(57)【要約】      （修正有）【課題】位相シンチレーションに対してリアル・タイム選別を実施する方法を提供する。【解決手段】位相シンチレーション・モニタにおいてサンプル時間期間中に位相シンチレーション・イベントを検出するステップと、位相シンチレーション・モニタにおける位相シンチレーション・イベントの検出に基づいて、関連する衛星測定データを今後の使用から除外するステップと、位相シンチレーション・モニタにおいて位相シンチレーション・イベントの終了を検出するステップと、位相シンチレーション・モニタによって検出された位相シンチレーション・イベントの終了後に収集された関連衛星測定データを再度受け入れるステップとを含む。【選択図】図３
2015180868.tif 000002
【特許請求の範囲】
【請求項１】
位相シンチレーションにリアル・タイム選別を実施する方法であって、
位相シンチレーション・モニタ（１３）においてサンプル時間期間中に位相シンチレーション・イベント（１６−１）を検出するステップと、
前記位相シンチレーション・モニタ（１３）における前記位相シンチレーション・イベント（１６−１）の検出に基づいて、関連する衛星測定データを今後の使用から除外するステップと、
前記位相シンチレーション・モニタ（１３）において前記位相シンチレーション・イベントの終了を検出するステップと、
前記位相シンチレーション・モニタ（１３）によって判定された前記位相シンチレーション・イベント（１６−１）の終了後に収集された関連衛星測定データを再度受け入れるステップと、
を含む、方法。
【請求項２】
請求項１記載の方法であって、更に、
蓄積デルタ距離データが入手可能である基準受信機衛星対（ＲＲｉ／ＳＶｊ）に対して、衛星（ＳＶ）運動およびＳＶクロック訂正キャリア・レート（ｃｒ＿ｃｏｒｒｅｃｔｅｄｉｊ（ｋ））を計算するステップと、
前記ＳＶ運動およびＳＶクロック訂正キャリア・レートから、他の全てのＳＶ運動およびＳＶクロック訂正キャリア・レート（ｃｒ＿ｃｏｒｒｅｃｔｅｄｉｊ（ｋ））を減算することによって、基準受信機クロックに対するＳＶ運動およびＳＶクロック訂正キャリア・レートを補償するステップと、
ｋ番目のサンプル時間期間における基準受信機／衛星対ＲＲｉ／ＳＶｊに対する基準受信機トレンド除去ＳＶ運動およびＳＶクロック訂正キャリア・レート（ｒｒＣＲｄｅｔｒｅｎｄ）ｉｊ（ｋ）を計算するステップと、
数値積分を使用して、キャリア位相推定値（Ｐ２ｉｎ）ｉｊ（ｋ）を計算するステップと、
基準受信機／衛星対（ＲＲｉ／ＳＶｊ）に対するキャリア位相推定値のサンプル平均（Ｐ２ａｖｇ）ｉｊ（ｋ）を計算するステップと、
現在（ｋ番目）のサンプル時間期間における基準受信機／衛星対（ＲＲｉ／ＳＶｊ）に対するキャリア位相推定値Ｐ２ｉｊ（ｋ）のばらつきを計算するステップと、
を含む、方法。
【請求項３】
請求項１記載の方法において、前記位相シンチレーション・モニタ（１３）において前記サンプル時間期間中に前記位相シンチレーション・イベント（１６−１）を検出するステップが、
少なくとも１つの基準受信機／衛星対（ＲＲｉ／ＳＶｊ）に対する現在の（ｋ番目）サンプル時間期間における位相シンチレーション・モニタ（１３）の識別子が、除外判断基準を満たすステップを含み、
前記位相シンチレーション・イベント（１６−１）の終了を検出するステップが、
除外基準受信機／衛星対（ＲＲｉ／ＳＶｊ＿ｅｘｃｌｕｄｅｄ）に対する現在の（ｋ番目）サンプル時間期間における位相シンチレーション・モニタ（１３）の識別子が、再受け入れ判断基準を満たすと判定するステップを含む、方法。
【発明の詳細な説明】
【背景技術】
【０００１】
[0001]  全地球ナビゲーション衛星システム（ＧＮＳＳ）用の地上型補強システム（ＧＢＡＳ）は、１つ以上のパラメータをブロードキャストし、ＧＮＳＳ受信機が、受信した衛星信号における種々の誤差を訂正するために、これらのパラメータを受信することができる。受信した衛星信号（全地球測位システム（ＧＰＳ）信号）において電離層シンチレーション攪乱が誤差の原因となるとき、シンチレーションの脅威が発生する。電離層（イオン）シンチレーションは、ＧＰＳ信号の振幅および／または位相の変動であり、信号が電離層における電離圏全電子数（ＴＥＣ）に局在する異常を通過するときに発生する。電離層（イオン）位相シンチレーションは、電離層を通過する信号におけるＧＰＳ信号振幅の急激な変動の原因となる。位相におけるこれらの擾乱は、地上局によってサービスされる空港において着陸しつつある飛行機について許容度外（out of tolerance）高精度接近情報を検出する、衛星システムの地上型補強システム（ＧＢＡＳ）の能力に影響を及ぼす虞がある。
【０００２】
[0002]  シンチレーションの発生は希ではないので、電離層シンチレーションは、地上局のサブシステムにおける衛星障害モニタの一部に対して「幻惑する」（blinding）運行シナリオまたは偽りの運行シナリオを生成することによって、衛星システムを脅威に晒す。加えて、電離層シンチレーションは、完全性モニタの基礎となる検査統計およびその関連する限界シグマ（bounding sigma）（σ）が、予期されるものとは著しく異なり、このためにその要求性能を満たすことができなくなる原因にもなり得る。
【発明の概要】
【０００３】
[0003]  本願は、位相シンチレーションのリアル・タイム選別プロセスを実施する方法に関する。この方法は、位相シンチレーション・モニタにおいてサンプル時間期間中において位相シンチレーション・イベントを検出するステップと、位相シンチレーション・モニタにける位相シンチレーション・イベントの検出に基づいて、関連する衛星測定を今後の使用（further use）から除外するステップと、位相シンチレーション・モニタにおける位相シンチレーション・イベントの終了を検出するステップと、位相シンチレーション・モニタによって判定された位相シンチレーション・イベントの終了後に収集された関連衛星測定データを再度受け入れるステップとを含む。
【図面の簡単な説明】
【０００４】
[0004]  図面は実施形態例を図示するに過ぎず、したがって範囲を限定すると解釈されないことを理解した上で、添付図面の使用により、更に具体性および詳細を加えて、実施形態例について説明する。
【図１】図１は、衛星システムに通信可能に結合された航空機、および本願にしたがって位相シンチレーションに対するリアル・タイム選別を行う地上型補強システムを示す。
【図２】図２は、本願にしたがって位相シンチレーション監視識別子（discriminator）を計算する方法の流れ図を示す。
【図３】図３は、本願にしたがって、位相シンチレーション監視識別子測定に関連する衛星測定データを除外すべきか否か判定するために使用される差分補正プロセッサ（ＤＣＰ）の処理図を示す。
【図４】図４は、本願にしたがって、一旦位相シンチレーション監視識別子（Ｐ２）測定値が除外レベルに達したなら、今後の使用から除外する方法の流れ図を示す。
【図５】図５は、本発明にしたがって、一旦位相シンチレーション監視識別子（Ｐ２）測定値がその再受け入れ判断基準を満たしたなら、衛星測定データの使用を再度受け入れる方法の流れ図を示す。
【０００５】
[0010]  共通の慣例にしたがって、説明する種々の特徴は同一の拡縮率で描かれるのではなく、実施形態例に関連する具体的な特徴を強調するように描かれるものとする。
【発明を実施するための形態】
【０００６】
[0011]  電離層位相シンチレーションは、ＧＰＳ衛星信号受信に影響を及ぼす可能性がある環境条件である。先に引用した問題に取り組むため、そして衛星システム・モニタを保護するために、本明細書において説明する選別プロセスおよびシステムは、潜在的にモニタの動作を変転させる可能性がある衛星測定データを排除する。本明細書において説明する選別プロセスは、シンチレーション・イベントの間に得られた衛星測定データを除外し、シンチレーション・イベントが終了したと判定された後に得られた衛星測定データを再度受け入れる。具体的には、本明細書では、位相シンチレーション・モニタについて記載し、更に位相シンチレーション・イベントによって影響を受けた衛星測定データを排除するための衛星測定データ除外および再受け入れプロセスについて記載する。この除外および再受け入れプロセスは、結果的に得られるＧＰＳ信号がシステムの完全性を保護し続けること、即ち、既存のモニタがそれらの完全性要件を満たすことに変わりはないことを確保する。この選別プロセスは、ＧＢＡＳの一部として使用される。ＧＢＡＳは、正確な着陸サービスを航空共同体（aviation community）に提供するために使用される。選別プロセスは、好ましくない衛星測定データを検出し、モニタにおけるその使用を禁止しつつ有効なデータが通過することを可能にする。有効な衛星測定データは、この選別プロセスを通過することを許され、システムの完全性を確保するために、残りの機能／モニタが適正に動作することを可能にする。
【０００７】
[0012]  本明細書において説明する位相シンチレーション・モニタは、位相シンチレーションについての衛星測定データを選別し（screen）、その使用がシステムの完全性を保護する能力を低下させる場合、その測定データを除去するように設計される。位相シンチレーションは、ＧＢＡＳにおいて使用される衛星測定毎に監視される。選別は、衛星測定におけるキャリア位相推定値の分散（variance）に基づく。
【０００８】
[0013]  図１は、衛星システム２０に通信可能に結合された航空機４０を示す。衛星システム２０は、複数の衛星２０（１〜Ｎ）と、通信可能に結合されたＧＢＡＳ１０とを含む。ＧＢＡＳ１０は、本願にしたがって、位相シンチレーションに基づくリアル・タイム選別を行う。航空機４０は、ここでは「空中運航体４０」（airborne vehicle）とも呼ばれる。ＧＢＡＳ１０は、受信した衛星信号を利用して、ある種の用途における位置精度を補強する。例えば、地上型補強システムは、高精度航空機接近用途（precision aircraft approach application）における位置精度を高めることができる。地上型補強システム１０は、飛行場における滑走路４６と関連がある地上局１１の一部である。地球の水平線を全体的に４５で表す。電離層１５は、衛星２０（１〜Ｎ）と地球４５との間にある大気の１つの層である。図１に示すように、全体的に１６−１および１６−２によって表される電離圏全電子数（ＴＥＣ）における局在異常が、電離層１５の部分にある。ＴＥＣ１６（１〜２）における局在異常は、全てのシンチレーション・イベントを示すが、ここでは「位相シンチレーション・イベント１６（１〜２）」のみをそのように呼ぶことにする。位相シンチレーション・イベント１６（１〜２）は、電離層１５の信号摂動の様相（signal-perturbing feature）であり、信号位相に影響を及ぼし、時間的に変化し、異なる時点において電離層１５の異なる部分で発見される。
【０００９】
[0014]  地上型補強システム１０は、少なくとも１つのモニタ７０、複数の基準受信機（ＲＲ）６０（１〜４）、および位相シンチレーション・モニタ１３を含む。位相シンチレーション・モニタ１３は、プロセッサ５０、記憶媒体８０、ソフトウェア８５、および任意のメモリ５５を含む。プロセッサ５０は、複数の基準受信機（ＲＲ）６０（１〜４）を通じて複数の衛星２０（１〜Ｎ）からの入力を受信するように、通信可能に結合される。当業者には周知であろうが、完全性モニタ７０およびＧＢＡＳ１０もプロセッサとインターフェースし、ソフトウェアおよびメモリを有する。図面を簡素化するために、これらは図１には示されない。
【００１０】
[0015]  第１基準受信機６０−１は、第２基準受信機６０−２から距離Ｌ１−２のところにある。第２基準受信機６０−２は、第３基準受信機６０−３から距離Ｌ２−３のところにある。第３基準受信機６０−３は、第４基準受信機６０−４から距離Ｌ３−４のところにある。第４基準受信機６０−４は、第１基準受信機６０−１から距離Ｌ４−１のところにある。基準受信機６０（１〜４）間の距離をここでは基準線と呼ぶ。
【００１１】
[0016]  また、少なくとも１つのモニタ７０は、ここでは「完全性モニタ７０」とも呼ぶ。前述のように、位相シンチレーションは、ある完全性モニタ７０の検査統計を低下させる可能性がある。本願は、これらのモニタ７０の位相シンチレーション脅威の影響の軽減を中心に据える。
【００１２】
[0017]  シンチレーション脅威を軽減するには、衛星測定値から望ましくないレベルの位相シンチレーションを選別する（screen）ことによって、完全性モニタ７０の変転を防止するリアル・タイム位相シンチレーション選別プロセスを実施する。位相シンチレーション監視識別子をここではＰ２と呼ぶ。また、リアル・タイム選別プロセスは、位相シンチレーション・イベントがいつ終了したかについても判定し、以前に排除された測定ソースから収集された衛星測定データを再度受け入れる。再受け入れが行われるのは、位相シンチレーション監視識別子が、以下で説明するように、再受け入れ判断基準を満たすときである。
【００１３】
[0018]  位相シンチレーション監視識別子Ｐ２を作成するために使用されるプロセスは、基準受信機（ＲＲ）と蓄積デルタ距離測定値が入手可能な被追跡衛星との全ての対について計算および評価される。蓄積デルタ距離という用語は、基準受信機／衛星対から受信された測定データを指し、最初に衛星信号を追跡して以来受信されてきたキャリア位相サイクルの蓄積から成る。キャリア位相サイクルはメートル単位に換算される。衛星から信号を受信するあらゆる基準受信機は、基準受信機／衛星対ＲＲｉ／ＳＶｊであると定義する。「宇宙船」（ＳＶ）および衛星という用語は、本明細書では相互交換可能に使用されるものとする。現在のサンプル時間期間における基準受信機／衛星対ＲＲｉ／ＳＶｊをＳＶｉｊとも呼ぶ。例えば、図１に示すように、衛星信号は、ワイヤレス通信リンク２２−１を通じて、第２基準受信機６０−２において第１衛星２０−１から受信されるので、第２基準受信機６０−２および第１衛星２０−１は基準受信機／衛星対ＲＲ２／ＳＶ１となる。勿論、第１基準受信機６０−１が第１衛星２０−１の受信エリア内にある場合、同じ衛星信号が第１基準受信機６０−１においても、ワイヤレス通信リンクを通じて、受信されるので、第１基準受信機６０−１および第１衛星２０−１は基準受信機／衛星対ＲＲ１／ＳＶ１を形成する。衛星の受信エリアとは、衛星から信号を受信することができる（露出される）地球のエリアである。受信エリアは、衛星が地球の軌道を進むにつれて、そして地球が自転するにつれて変化する。ある高度よりも下においてこれらの衛星信号の受け入れを制限するために、受信マスクを適用することもできる。ブロードキャスト信号はワイヤレス通信リンク２２−５を通じて航空機４０に送られ、位相シンチレーション監視識別子１３に基づいて除外されない衛星からのデータだけを含む。
【００１４】
[0019]  図２は、本願にしたがって位相シンチレーション監視識別子（Ｐ２）を計算する方法２００の流れ図を示す。方法２００は、位相シンチレーション・モニタ１３（図１）においてソフトウェア８５を実行するプロセッサ５０によって、各サンプル時間期間において、基準受信機／衛星対ＲＲｉ／ＳＶｊの各々について実施される。位相シンチレーション監視識別子（Ｐ２）は、ここでは、「位相シンチレーション監視識別子（Ｐ２）測定値」、および「キャリア位相推定値のばらつき」、ＲＲｉおよびＳＶｊに対するサンプル時刻ｋにおけるＰ２ｉｊ（ｋ）、ならびに「キャリアの位相のばらつき」とも呼ぶことにする。「キャリア・レート」（carrier rate）という用語は、蓄積デルタ距離の変化率であり、ここでは、「デルタ蓄積デルタ距離」（delta accumulated delta range）、または「デルタ・キャリア・レート」、または「Δキャリア」として相互交換可能に使用することにする。位相シンチレーション監視識別子（Ｐ２）ｉｊ（ｋ）は、蓄積デルタ距離（ＡＤＲ）データが入手可能な基準受信機／衛星対ＲＲｉ／ＳＶｊ毎に計算および評価される。
【００１５】
[0020]  ブロック２０２において、基準受信機クロック訂正キャリア・レートを、蓄積デルタ距離データが入手可能なそれぞれの基準受信機／衛星対について計算する。プロセッサ５０（図１）は、基準受信機／衛星対についての蓄積デルタ距離データが入手可能か否か判定する。一般に、以下で示すように、ＳＶｊならびにＲＲｉのＳＶ運動およびＳＶクロック訂正キャリア・レートを時間期間ｋにおいて計算する。ＳＶ運動およびＳＶクロック訂正キャリア・レート（SV motion and SV clock corrected carrier rate）（ｃｒ＿ｃｏｒｒｅｃｔｅｄｉｊ（ｋ））を受信エリア内にある基準受信機／衛星対ＲＲｉ／ＳＶｊ毎に計算する。
【００１６】
[0021]  以下のように、ＳＶ運動およびＳＶクロック訂正キャリア・レート値を計算する。
【００１７】
【数１】
2015180868.tif 000003
【００１８】
ここで、
ｃｒ＿ｃｏｒｒｅｃｔｅｄｉｊ（ｋ）は、時刻ｋにおけるＲＲｉ／ＳＶｊに対するＳＶ運動およびＳＶクロック訂正キャリア・レートである。
【００１９】
ａｃｃｉｊ（ｋ）は、時刻ｋにおけるＲＲｉ／ＳＶｊ上の蓄積デルタ距離である。
ａｃｃｉｊ（ｋ−１）は、時刻ｋ−１におけるＲＲｉ／ＳＶｊ上の蓄積デルタ距離である。
【００２０】
ｒａｎｇｅｉｊ（ｋ）は、時刻ｋにおけるＳＶｊからＲＲｉまでの距離である。
ｒａｎｇｅｉｊ（ｋ−１）は、時刻ｋ−１におけるＳＶｊからＲＲｉまでの距離である。
【００２１】
ｓｖｃｌｋｉｊ（ｋ）は、時刻ｋにおいてＲＲｉが見るエフェメリス・クロック値である。
ｓｖｃｌｋｉｊ（ｋ−１）は、時刻ｋ−１においてＲＲｉが見るエフェメリス・クロック値である。
【００２２】
ｔｒｏｐｏｉｊ（ｋ）は、時刻ｋにおけるＲＲｉ／ＳＶｊに対する対流圏遅延推定値である。
ｔｒｏｐｏｉｊ（ｋ−１）は、時刻ｋ−１におけるＲＲｉ／ＳＶｊに対する対流圏遅延推定値である。
【００２３】
[0022]  このように、時刻ｋにおけるＳＶｊおよＲＲｉのＳＶ運動およびＳＶクロック訂正キャリア測定値を計算する。
[0023]  このキャリア・レートから、このＲＲ上における全ての他のＳＶ運動およびＳＶクロック訂正キャリア・レートの平均を減算して、基準受信機クロック訂正キャリア・レート（ｒｒＣＲ）ｉｊ（ｋ）を求めることによって、ＳＶ運動およびＳＶクロック訂正キャリア・レートを基準受信機クロックに対して補償する。ｋ番目のサンプル時間期間における各基準受信機／衛星対ＲＲｉ／ＳＶｊに対する基準受信機クロック訂正キャリア・レートを、次のように計算する。
【００２４】
【数２】
2015180868.tif 000004
【００２５】
ここで、
ｒｒＣＲｉｊ（ｋ）は、時刻ｋにおけるＲＲｉ／ＳＶｊに対する基準受信機クロック訂正キャリア・レートである。
【００２６】
ｃｒ＿ｃｏｒｒｅｃｔｅｄｉｊ（ｋ）は、時刻ｋにおけるＲＲｉ／ＳＶｊに対するＳＶ運動およびＳＶクロック訂正キャリア・レートである。
Ｎｓｖは、時刻ｋにおいてｃｒ＿ｃｏｒｒｅｃｔｅｄｉｊ値を有するＲＲｉにおける衛星の数である。
【００２７】
[0024]  ｋ番目のサンプル時間期間における各基準受信機／衛星対ＲＲｉ／ＳＶｊに対する基準受信機クロック訂正キャリア・レートの移動平均は、次のように計算する。この実施形態の一実施態様では、各サンプル時間期間の長さを１／２秒として、４０サンプルにわたって移動平均を計算するが、他のサンプル数および長さも可能である。
【００２８】
【数３】
2015180868.tif 000005
【００２９】
ここで、
ＡＶＧ＿ＣＲｉｊ（ｋ）は、時刻ｋにおけるＲＲｉ／ＳＶｊに対する基準受信機クロック訂正キャリア・レートの移動平均である。
【００３０】
ｒｒＣＲｉｊ（ｋ−ｎ）は、時刻ｋ−ｎにおけるＲＲｉ／ＳＶｊに対する基準受信機クロック訂正キャリア・レートである。
[0025]  時刻ｋにおける基準受信機トレンド除去ＳＶ運動およびＳＶクロック訂正キャリア・レート（ｒｒＣＲｄｅｔｒｅｎｄ）ｉｊ（ｋ）は、基準受信機クロック訂正キャリア・レートｒｒＣＲｉｊ（ｋ）、および平均基準受信機クロック訂正キャリア・レート（ＡＶＧ＿ＣＲ）ｉｊ（ｋ）の差である。基準受信機トレンド除去ＳＶ運動およびＳＶクロック訂正キャリア・レート（ｒｒＣＲｄｅｔｒｅｎｄ）ｉｊ（ｋ）は、以下のように、ｋ番目のサンプル時間期間における基準受信機／衛星対ＲＲｉ／ＳＶｊ毎に計算する。
【００３１】
【数４】
2015180868.tif 000006
【００３２】
ここで、
（ｒｒＣＲｄｅｔｒｅｎｄ）ｉｊ（ｋ）は、時刻ｋにおけるＲＲｉ／ＳＶｊに対する基準受信機トレンド除去ＳＶ運動およびＳＶクロック訂正キャリア・レートである。
【００３３】
ｒｒＣＲｉｊ（ｋ−ｎ）は、時刻ｋ−ｎにおけるＲＲｉ／ＳＶｊに対する基準受信機クロック訂正キャリア・レートである。
ＡＶＧ＿ＣＲｉｊ（ｋ）は、時刻ｋにおけるＲＲｉ／ＳＶｊに対する基準受信機クロック訂正キャリア・レートの移動平均である。
【００３４】
[0026]  このように、基準受信機トレンド除去ＳＶ運動およびＳＶクロック訂正キャリア・レート（ｒｒＣＲｄｅｔｒｅｎｄ）ｉｊ（ｋ）は、時刻ｋにおいて有効な基準受信機／衛星対ＲＲｉ／ＳＶｊ毎に計算する。
【００３５】
[0027]  ブロック２０４において、台形積分を使用してキャリア位相推定値（Ｐ２ｉｎ）ｉｊ（ｋ）を計算する。ブロック２０２において計算した基準受信機クロック訂正キャリア・レートを、基準受信機／衛星対ＲＲｉ／ＳＶｊの各々について積分する。このように、現在のサンプル時間期間ｋにおける基準受信機ｉ上における衛星測定値ｊのキャリア位相推定値が、（Ｐ２ｉｎ）ｉｊ（ｋ）で表され、現在のサンプル時間期間における基準受信機／衛星対ＲＲｉ／ＳＶｊの各々について計算する。「キャリア位相推定値」および「キャリア位相」は、本明細書では相互交換可能に使用されるものとし、双方共（Ｐ２ｉｎ）ｉｊとして表される。キャリア位相推定値（Ｐ２ｉｎ）ｉｊ（ｋ）は、基準受信機／衛星対ＲＲｉ／ＳＶｊ毎に、基準受信機クロック訂正キャリア・レートを積分することによって生成される。
【００３６】
[0028]  キャリア位相推定値（Ｐ２ｉｎ）ｉｊ（ｋ）は、以下のように数値積分（台形積分を使用する）によって生成する。
【００３７】
【数５】
2015180868.tif 000007
【００３８】
ここで、
（Ｐ２ｉｎ）ｉｊ（ｋ）は、時刻ｋにおけるＲＲｉ／ＳＶｊに対するキャリア位相推定値である。
【００３９】
（ｒｒＣＲｄｅｔｒｅｎｄ）ｉｊ（ｋ）は、時刻ｋにおけるＲＲｉ／ＳＶｊに対する基準受信機トレンド除去ＳＶ運動およびＳＶクロック訂正キャリア・レートである。
Δｔは、時間期間の長さ（この実施形態では、１／２秒）である。
【００４０】
ＣＮＶは、１／２秒毎のメートルを１秒毎のメートルに変換する（Δｔが０．５秒のとき２．０）。
[0029]  ブロック２０６において、キャリア位相推定値の平均（Ｐ２ａｖｇ）ｉｊ（ｋ）を、以下のように、基準受信機／衛星対ＲＲｉ／ＳＶｊ毎に計算する。
【００４１】
【数６】
2015180868.tif 000008
【００４２】
ここで、
（Ｐ２ａｖｇ）ｉｊ（ｋ）は、時刻ｋにおけるＲＲｉ／ＳＶｊに対する平均キャリア位相推定値である。
【００４３】
（Ｐ２ｉｎ）ｉｊ（ｋ−ｎ）は、時刻ｋ−ｎにおけるＲＲｉ／ＳＶｊに対するキャリア位相推定値である。
Ｋｐは、平均キャリア位相推定値におけるサンプル数である（この実施形態では３０であるが、他の値も可能である）。
【００４４】
[0030]  ブロック２０８において、現在のサンプル時間期間において、基準受信機／衛星対ＲＲｉ／ＳＶｊ毎に、キャリア位相推定値のばらつきを計算する。「キャリア位相推定値のばらつき」および「位相シンチレーション監視識別子測定値」という用語は、本明細書では相互交換可能に使用されるものとし、双方ともｋ番目のサンプルに対して、（Ｐ２）ｉｊ（ｋとして表される。
【００４５】
[0031]  現在のサンプル時間期間におけるＲＲｉ上におけるＳＶｊに対するキャリア位相推定値の各衛星のばらつきは、以下のように計算する。
【００４６】
【数７】
2015180868.tif 000009
【００４７】
ここで、
（Ｐ２）ｉｊ（ｋ）は、時刻ｋにおけるＲＲｉ／ＳＶｊに対するキャリア位相推定値のばらつきである。
【００４８】
（Ｐ２ａｖｇ）ｉｊ（ｋ）は、時刻ｋにおけるＲＲｉ／ＳＶｊに対する平均キャリア位相推定値である。
（Ｐ２ｉｎ）ｉｊ（ｋ−ｎ）は、時刻ｋ−ｎにおけるＲＲｉ／ＳＶｊに対するキャリア位相推定値である。
【００４９】
Ｋｐは、平均キャリア位相推定値におけるサンプル数である（この実施形態では３０であるが、他の値も可能である）。
[0032]  図３は、本願にしたがって、関連する位相シンチレーション監視識別子測定値を除外すべきか否か判定するために使用される差分訂正プロセッサ（ＤＣＰ）の処理図１５０を示す。サンプル毎に、基準受信機６０（１〜４）（図１）は、生の蓄積デルタ距離キャリアに基づく測定値および測定時刻をプロセッサ５０（図１）に衛星（宇宙船（ＳＶ））２０（１〜Ｎ）から供給する。また、プロセッサ５０は、入力として、Ｐ２除外閾値、Ｐ２再受け入れ閾値、保持時間、およびキュー・サイズも受信する。この実施形態の一実施態様では、この入力データは、プロセッサ５０に入力される前に、メモリ５５に格納される。Ｐ２除外閾値をここでは「測定除外閾値」、「除外閾値」、および「ＳＶ測定除外閾値」とも呼ぶ。これらは全て、「Ｐ２＿ＴＨＲＥＳＨＯＬＤ」として表される。Ｐ２再受け入れ閾値をここでは「測定値再受け入れ閾値」、「再受け入れ閾値」、および「ＳＶ測定値再受け入れ閾値」とも呼ぶ。これらは全て、Ｐ２＿ＲＥＡＤＭＩＴ  ＴＨＲＥＳＨＯＬＤとして表される。
【００５０】
[0033]  各サンプル時間期間の間、プロセッサ５０（図１）は以上の入力を受信し、位相シンチレーション監視ソフトウェア１６０を実行して、関連する衛星測定データを除外すべきか否か判定する。関連する衛星測定データに対する位相シンチレーション監視識別子測定値を除外すべきである基準受信機／衛星対ＲＲｉ／ＳＶｊを、ここでは、除外基準受信機／衛星対ＲＲｉ／ＳＶｊ（またはＲＲｉ／ＳＶｊ＿ｅｘｃｌｕｄｅｄ）と呼ぶ。図３に示すように、位相監視ソフトウェア１６０は、除外される衛星測定値のリストを出力する。これらの衛星測定値は、入手可能な測定値を判定するために使用される。入手可能な測定値は、現行の１組の測定値から、位相シンチレーション監視ソフトウェア１６０から出力される除外測定値のリストを差し引いたものとなる。
【００５１】
[0034]  位相シンチレーション監視ソフトウェア１６０の実行のフローを、図４および図５の方法４００および５００として示す。図４は、本願にしたがって、位相シンチレーション監視識別子（Ｐ２）測定値と関連する衛星測定データを除外する方法４００の流れ図を示す。「位相シンチレーション監視識別子（Ｐ２）測定値」および「キャリア位相推定値のばらつき」、Ｐ２ｉｊ（ｋ）測定値は、ここでは相互交換可能に使用されるものとする。方法４００は、サンプル時間期間ｋ毎に実行される。この実施形態の一実施態様では、サンプル時間期間は１／２秒の長さを有するが、サンプル時間期間には他の長さも可能である。
【００５２】
[0035]  ブロック４０２において、本プロセスが開始する。ブロック４０４において、現在のサンプル時間期間における基準受信機／衛星対ＲＲｉ／ＳＶｊ（即ち、ＳＶｉｊ）が、現在除外されているか否か判定する。ＳＶｉｊが現在除外されている場合、現在のサンプル時間期間ｋでは基準受信機ｉにおいて衛星測定値ｊの衛星測定データＳＶｉｊ（ｋ）が現在除外されていることになる。
【００５３】
[0036]  ＳＶｉｊ（ｋ）が現在除外されている場合、以前のサンプル時間期間（例えば、（ｋ−１）番目のサンプル時間期間から（ｋ−ｎ）番目のサンプル期間までの間（ここで「ｎ」は正の整数である）においてＳＶｉｊが除外されたことになり、フローはブロック４０６に進み、再受け入れチェックを実行する。この場合、フローは、ブロック４０６から図５における方法５００のブロック５０２に進む。図５については以下で説明する。
【００５４】
[0037]  ＳＶｉｊ（ｋ）が現在除外されていない場合、フローはブロック４０４からブロック４０８に進む。ブロック４０８において、位相シンチレーション監視識別子測定値Ｐ２ｉｊ（ｋ）がその除外判断基準を満たしているか否か判定する。このプロセスは、基準受信機／衛星対ＲＲｉ／ＳＶｊの各々に対して繰り返される。
【００５５】
[0038]  この実施形態の一実施態様では、除外判断基準を満たすのは、Ｐ２ｉｊ（ｋ）が除外閾値（例えば、Ｐ２閾値）よりも大きい場合である。この実施形態の他の実施態様では、除外判断基準を満たすのは、Ｐ２ｉｊ（ｋ）がＰ２閾値以上である場合である。Ｐ２閾値（除外閾値）はメモリ５５（図１）および／またはプロセッサ５０（図１）に格納されている。いずれの場合でも、それぞれの基準受信機／衛星対ＲＲｉ／ＳＶｊに対する位相シンチレーション監視識別子測定値Ｐ２ｉｊ（ｋ）は、各々、除外閾値と比較される。他の除外判断基準も可能である。
【００５６】
[0039]  Ｐ２ｉｊが基準受信機／衛星対ＲＲｉ／ＳＶｊの１つ（以上）についてその除外判断基準を満たしている場合、この１つ（以上）の基準受信機／衛星対ＲＲｉ／ＳＶｊは、除外される基準受信機／衛星対ＲＲｉ／ＳＶｊ＿ｅｘｃｌｕｄｅｄとなる。この実施形態の一実施態様では、基準受信機／衛星対ＲＲｉ／ＳＶｊの２つ以上が、Ｐ２閾値以上のＰ２ｉｊ（ｋ）を有する。この場合、これらの基準受信機／衛星対ＲＲｉ／ＳＶｊは、除外される基準受信機／衛星対ＲＲｉ／ＳＶｊ＿ｅｘｃｌｕｄｅｄであると定義され、この除外される基準受信機／衛星対ＲＲｉ／ＳＶｊ＿ｅｘｃｌｕｄｅｄにおける衛星から送られるデータは、下流側のモニタ７０において使用されず、ワイヤレス通信リンク２２−５（図１）を通じて送られるブロードキャストの一部にはならない。次いで、フローは１つ以上の除外される基準受信機／衛星対ＲＲｉ／ＳＶｊ＿ｅｘｃｌｕｄｅｄのためにブロック４１０に進む。ブロック４１０において、この１つ（以上）の基準受信機／衛星対ＲＲｉ／ＳＶｊが、除外される基準受信機／衛星対ＲＲｉ／ＳＶｊ＿ｅｘｃｌｕｄｅｄとして定義される。このように、除外される基準受信機／衛星対ＲＲｉ／ＳＶｊ＿ｅｘｃｌｕｄｅｄについての位相シンチレーション監視識別子測定値が除外閾値（即ち、Ｐ２＿ＴＨＲＥＳＨＯＬＤ）以上であると判定されたとき（除外閾値との比較のときに）、基準受信機／衛星対における衛星からの、除外される基準受信機／衛星対ＲＲｉ／ＳＶｊ＿ｅｘｃｌｕｄｅｄについてのデータは、もはや下流側のモニタ７０において使用されず、ワイヤレス通信リンク２２−５（図１）を通じて送られるブロードキャストの一部とはならない。次いで、フローはブロック４１２に進み、フローは現在のサンプル時間期間に対するプロセスを終了する。
【００５７】
[0040]  Ｐ２ｉｊがその除外判断基準を満たしていない場合、フローはブロック４０８からブロック４１２に進み、この１つ（以上）の非除外基準受信機／衛星対ＲＲｉ／ＳＶｊに対するフローは、現在のサンプル時間期間に対するプロセスを終了する。この実施形態の一実施態様では、Ｐ２ｉｊ（ｋ）が基準受信機／衛星対ＲＲｉ／ＳＶｊの１つ（以上）に対してＰ２除外閾値未満である場合、基準受信機／衛星対ＲＲｉ／ＳＶｊのこの１つ（以上）は除外されず、衛星測定データは、下流側のモニタ７０において使用され続け、ワイヤレス通信リンク２２−５（図１）を通じて送られるブロードキャストの一部のままになっている。
【００５８】
[0041]  図５は、本願にしたがって位相シンチレーション監視識別子（Ｐ２）測定値に関連する衛星測定データを再度受け入れる方法５００の流れ図を示す。衛星測定データが再度受け入れられるのは、位相シンチレーション・イベントが終了したと判定されたときである。方法５００は、位相シンチレーション・モニタ１３（図１）においてソフトウェア８５を実行するプロセッサ５０によって、基準受信機／衛星対ＲＲｉ／ＳＶｊの各々に対して実施される。
【００５９】
[0042]  ブロック５０２において、本プロセスが開始する。ブロック５０４において、ＳＶｉｊ（ｋ）が現在除外されているか否か判定する。ＳＶｉｊ（ｋ）が現在除外されていない場合（例えば、基準受信機／衛星対ＲＲｉ／ＳＶｊに対するＳＶｉｊがｋ番目のサンプル時間期間において除外されていない場合）、フローはブロック５０６に進み、除外チェックを実行する。現在のサンプル時間期間は、ｋ番目のサンプル時間期間である。この場合、フローは、以上で説明した処理のために、ブロック５０６から図４における方法４００のブロック４０２に進む。
【００６０】
[0043]  ＳＶｉｊ（ｋ）が現在除外されている場合、フローはブロック５０８に進む。ブロック５０８において、除外基準受信機／衛星対ＲＲｉ／ＳＶｊ＿ｅｘｃｌｕｄｅｄに対する現在のサンプル時間期間におけるＰ２ｉｊ（ｋ）が、再受け入れ判断基準を満たすか否か判定する。
【００６１】
[0044]  この実施形態の一実施態様では、再受け入れ判断基準が満たされるのは、現在のサンプル時間期間における位相シンチレーション監視識別子Ｐ２ｉｊ（ｋ）がＰ２再受け入れ閾値未満である場合である。Ｐ２再受け入れ閾値は、予め選択され、メモリ５５（図１）またはプロセッサ５０（図１）に格納される。プロセッサ５０は、除外基準受信機／衛星対ＲＲｉ／ＳＶｊ＿ｅｘｃｌｕｄｅｄに対する位相シンチレーション監視識別子Ｐ２ｉｊ（ｋ）をＰ２再受け入れ閾値と比較する。他の再受け入れ判断基準も可能である。
【００６２】
[0045]  この実施形態の他の実施態様では、再受け入れ判断基準が満たされるのは、Ｐ２ｉｊ（ｋ）が、現在のサンプル時間期間の前で基準受信機／衛星対ＲＲｉ／ＳＶｊが除外された後に発生した、予め選択された数Ｍ（Ｍは正の整数）個のサンプルにわたって、現在のサンプル時間期間におけるＰ２再受け入れ閾値未満である場合である。この実施形態では、プロセッサ５０は、除外基準受信機／衛星対ＲＲｉ／ＳＶｊ＿ｅｘｃｌｕｄｅｄに対する位相シンチレーション監視識別子Ｐ２ｉｊ（ｋ）をＰ２再受け入れ閾値と比較し、基準受信機／衛星対ＲＲｉ／ＳＶｊが除外されて以降にこの閾値よりも小さいサンプルの数を計数し、基準受信機／衛星対ＲＲｉ／ＳＶｊが除外されて以降に計数した閾値未満のサンプル数を、格納されている予め選択された数Ｍと比較する。予め選択された数Ｍは、メモリ５５（図１）またはプロセッサ５０（図１）に格納されている。
【００６３】
[0046]  この実施形態の更に他の実施態様では、再受け入れ判断基準が満たされるのは、現在のサンプル時間期間におけるＰ２ｉｊ（ｋ）が、現在のサンプル時間期間の直前に連続的に発生した、予め選択された数Ｍ個の連続サンプルについて、Ｐ２再受け入れ閾値未満である場合である。この実施形態の更に他の実施態様では、再受け入れ判断基準が満たされるのは、予め選択された数Ｍ個のサンプルが、これらの連続サンプルに対するＰ２ｉｊ（ｋ）の値には関係なく、計数されたときである。これらの任意の再受け入れ判断基準は、位相シンチレーション・イベントが終了したことの指示として使用される。位相シンチレーション・イベントが終了したことを示す他の再受け入れ判断基準も可能である。
【００６４】
[0047]  Ｐ２ｉｊ（ｋ）が除外基準受信機／衛星対ＲＲｉ／ＳＶｊ＿ｅｘｃｌｕｄｅｄに対する再受け入れ判断基準を満たす場合、フローはブロック５１０に進む。ブロック５１０において、除外基準受信機／衛星対の衛星測定値はもはや除外されない。フローはブロック５１０からブロック５１２に進み、フローは現在のサンプル時間期間に対するプロセス（ブロック５１２）を終了する。
【００６５】
[0048]  ブロック５０８において、Ｐ２ｉｊ（ｋ）が、除外基準受信機／衛星対ＲＲｉ／ＳＶｊ＿ｅｘｃｌｕｄｅｄに対して再受け入れ判断基準を満たしていないと判定した場合、除外基準受信機／衛星対ＲＲｉ／ＳＶｊ＿ｅｘｃｌｕｄｅｄの衛星測定値は相変わらず除外される（除外されたままである）。フローはブロック５１２に進む。このように、シンチレーション・イベントの間に基準受信機／衛星対ＲＲｉ／ＳＶｊから得られた衛星測定データは、ナビゲーション・システムにおける使用から除外され、この除外基準受信機／衛星対ＲＲｉ／ＳＶｊ＿ｅｘｃｌｕｄｅｄからの測定データは、シンチレーション・イベントがこの除外基準受信機／衛星対ＲＲｉ／ＳＶｊ＿ｅｘｃｌｕｄｅｄに対して終了した後に、再度受け入れられる。
【００６６】
[0049]  この実施形態の一実施態様では、キャリア位相推定値を使用して、位相シンチレーションの存在／不在が、４機の基準受信機から判定される。この実施形態の他の実施態様では、位相シンチレーションを検出するために使用されるキャリア位相推定値は、１０個の連続サンプルに対する０．０５（単位なし）のＰ２除外閾値および０．０１３メートルの再受け入れ閾値に基づく。
【００６７】
[0050]  本明細書において説明した方法および技法は、ディジタル電子回路において、あるいは少なくとも１つのプロセッサ（例えば、プログラマブル・プロセッサ、特殊目的プロセッサ、コンピュータのような汎用プロセッサ、または図１のプロセッサ５０）、ファームウェア、ソフトウェア、あるいはこれらの任意の組み合わせによって実現することができる。これらの技法を具体化する装置は、しかるべき入力および出力デバイス、プロセッサ、ならびにプロセッサによる実行のためのプログラム命令を実体的に（tangibly）具体化する記憶媒体を含めばよい。これらの技法を具体化するプロセスは、入力データに対して動作してしかるべき出力を生成することによって所望の機能を実行するために命令のプログラムを実行する少なくとも１つのプロセッサによって実行すればよい。有利なのは、これらの技法が、プログラマブル・システム上で実行可能な１つ以上のプログラムにおいて実現できることである。プログラマブル・システムは、データ記憶システムからデータおよび命令を受け、データおよび命令をデータ記憶システムに送信するように結合された少なくとも１つのプロセッサ、少なくとも１つの入力デバイス、ならびに少なくとも1つの出力デバイスを含む。一般に、プロセッサは、読み取り専用メモリおよび／またはランダム・アクセス・メモリから命令およびデータを受ける。
【００６８】
[0051]  コンピュータ・プログラムおよびデータを実体的に具体化するのに適した記憶デバイスには、一例としてＥＰＲＯＭ、ＥＥＰＲＯＭ、およびフラッシュ・メモリ・デバイスのような半導体メモリ・デバイスを含むあらゆる形態の不揮発性メモリ、内部ハード・ディスクおよびリムーバブル・ディスクのような磁気ディスク、およびＤＶＤディスクが含まれる。以上の内いずれもが、特別に設計された特定用途集積回路（ＡＳＩＣ）によって補足されても、またはこれらに組み込まれてもよい。
実施形態例
[0052]  例１は、位相シンチレーションに対してリアル・タイム選別プロセスを実施する方法を含む。この方法は、位相シンチレーション・モニタにおいてサンプル時間期間中に位相シンチレーション・イベントを検出するステップと、位相シンチレーション・モニタにおける位相シンチレーション・イベントの検出に基づいて、関連する衛星測定データを今後の使用から除外するステップと、位相シンチレーション・モニタにおいて位相シンチレーション・イベントの終了を検出するステップと、位相シンチレーション・モニタによって判定される位相シンチレーション・イベントの終了後に収集された関連衛星測定データを再度受け入れるステップとを含む。
【００６９】
[0053]  例２は、例１の方法を含み、更に、蓄積デルタ距離データが入手可能な基準受信機／衛星対に対して衛星（ＳＶ）運動およびＳＶクロック訂正キャリア・レートを計算するステップと、他の全てのＳＶ運動およびＳＶクロック訂正キャリア・レートの平均をＳＶ運動およびＳＶクロック訂正キャリア・レートから減算することによって、ＳＶ運動およびＳＶクロック訂正キャリア・レートを基準受信機クロックに対して補償するステップと、サンプル時間期間における基準受信機／衛星対に対する基準受信機トレンド除去ＳＶ運動およびＳＶクロック訂正キャリア・レートを計算するステップとを含む。
【００７０】
[0054]  例３は、例２の方法を含み、更に、数値積分を使用してキャリア位相推定値を計算するステップと、基準受信機／衛星対に対してキャリア位相推定値のサンプル平均を計算するステップと、現在のサンプル時間期間における基準受信機／衛星対に対するキャリア位相推定値のばらつきを計算するステップとを含む。
【００７１】
[0055]  例４は、例１〜３のいずれかの方法を含み、位相シンチレーション・モニタにおいてサンプル時間期間中に位相シンチレーション・イベントを検出するステップが、少なくとも１つの基準受信機／衛星対に対する現サンプル時間期間における位相シンチレーション監視識別子が除外判断基準を満たすと判定するステップを含む。
【００７２】
[0056]  例５は、例４の方法を含み、少なくとも１つの基準受信機／衛星対に対する現在のサンプル時間期間における位相シンチレーション監視識別子が、除外判断基準を満たすことを判定するステップが、少なくとも１つの基準受信機／衛星対に対する現在のサンプル時間期間における位相シンチレーション監視識別子が、除外閾値よりも大きい、または除外閾値に等しい、の内一方であることを判定するステップを含む。
【００７３】
[0057]  例６は、例１〜例５のいずれかの方法を含み、位相シンチレーション・イベントの終了を検出するステップが、除外基準受信機／衛星対に対する現在のサンプル時間期間にける位相シンチレーション監視識別子が、再受け入れ判断基準を満たすことを判定するステップを含む。
【００７４】
[0058]  例７は、例６の方法を含み、除外基準受信機／衛星対に対する現在のサンプル時間期間における位相シンチレーション監視識別子が、再受け入れ判断基準を満たすことを判定するステップが、除外基準受信機／衛星対に対する現在のサンプル時間期間における位相シンチレーション監視識別子が、再受け入れ閾値未満であることを判定するステップを含む。
【００７５】
[0059]  例８は、例６および例７のいずれかの方法を含み、除外基準受信機／衛星対に対する現在のサンプル時間期間にける位相シンチレーション監視識別子が、再受け入れ判断基準を満たすことを判定するステップが、除外基準受信機／衛星対に対する現在のサンプル時間期間における位相シンチレーション監視識別子が、予め選択された数のサンプルにわたって再受け入れ閾値未満であることを判定するステップを含む。
【００７６】
[0060]  例９は、例６〜例８のいずれかの方法を含み、除外基準受信機／衛星対に対する現在のサンプル時間期間にける位相シンチレーション監視識別子が、再受け入れ判断基準を満たすことを判定するステップが、除外基準受信機／衛星対に対する現在のサンプル時間期間における位相シンチレーション監視識別子が、予め選択された数の連続サンプルにわたって再受け入れ閾値未満であることを判定するステップを含む。
【００７７】
[0061]  例１０は、位相シンチレーションに対してリアル・タイム選別を行う位相シンチレーション・モニタを含む。この位相シンチレーション・モニタは、複数の基準受信機から入力を受信するように通信可能に結合された少なくとも１つのプロセッサと、少なくとも１つのプロセッサによる実行のためにプログラム命令を有形に具体化する記憶媒体とを含み、このプログラム命令は、少なくとも１つのプロセッサによって実行されたときに、サンプル時間期間中に位相シンチレーション・イベントを検出し、位相シンチレーション・イベントの検出に基づいて、関連する衛星測定データを今後の使用から除外し、位相シンチレーション・イベントの終了を検出し、位相シンチレーション・イベントの終了後に収集された関連衛星測定データを再度受け入れるように動作可能である。
【００７８】
[0062]  例１１は、例１０の位相シンチレーション・モニタを含み、サンプル測定の間位相シンチレーション・イベントを検出するプログラム命令が、更に、少なくとも１つのプロセッサによって実行されたときに、少なくとも１つの基準受信機／衛星対に対して現在のサンプル時間期間における位相シンチレーション監視識別子が、除外判断基準を満たすと判定するように動作可能である。
【００７９】
[0063]  例１２は、例１１の位相シンチレーション・モニタを含み、少なくとも１つの基準受信機／衛星対に対して現在のサンプル時間期間における位相シンチレーション監視識別子が、除外判断基準を満たすと判定するように動作可能なプログラム命令が、少なくとも１つのプロセッサによって実行されたときに、少なくとも１つの基準受信機／衛星対に対して現在のサンプル時間期間における位相シンチレーション監視識別子が、除外閾値よりも大きい、または除外閾値に等しい、の内一方であると判定するように動作可能なプログラム命令を含む。
【００８０】
[0064]  例１３は、例１０〜例１２のいずれかの位相シンチレーション・モニタを含み、位相シンチレーション・イベントの終了を検出するように      動作可能なプログラム命令が、少なくとも１つのプロセッサによって実行されたときに、除外基準受信機／衛星対に対する現在のサンプル時間期間における位相シンチレーション監視識別子が、再受け入れ判断基準を満たすと判定するように動作可能なプログラム命令を含む。
【００８１】
[0065]  例１４は、例１３の位相シンチレーション・モニタを含み、除外基準受信機／衛星対に対する現在のサンプル時間期間における位相シンチレーション監視識別子が、再受け入れ判断基準を満たすと判定するように動作可能なプログラム命令が、少なくとも１つのプロセッサによって実行されたときに、除外基準受信機／衛星対に対する現在のサンプル時間期間における位相シンチレーション監視識別子が、再受け入れ閾値未満であると判定するように動作可能なプログラム命令を含む。
【００８２】
[0066]  例１５は、例１３および例１４のいずれかの位相シンチレーション・モニタを含み、除外基準受信機／衛星対に対する現在のサンプル時間期間における位相シンチレーション監視識別子が、再受け入れ判断基準を満たすと判定するように動作可能なプログラム命令が、少なくとも１つのプロセッサによって実行されたときに、除外基準受信機／衛星対に対する現在のサンプル時間期間における位相シンチレーション監視識別子が、予め選択された数のサンプルにわたって、再受け入れ閾値未満であると判定するように動作可能なプログラム命令を含む。
【００８３】
[0067]  例１６は、例１０〜例１５のいずれかの位相シンチレーション・モニタを含み、前述のプログラム命令が、更に、少なくとも１つのプロセッサによって実行されたときに、位相シンチレーション・モニタに、少なくとも１つの基準受信機／衛星対に対して現在のサンプル時間における位相シンチレーション監視識別子を計算させるように動作可能である。
【００８４】
[0068]  例１７は、例１６の位相シンチレーション・モニタを含み、少なくとも１つの基準受信機／衛星対に対して現在のサンプル時間における位相シンチレーション監視識別子を計算するプログラム命令が、少なくとも１つのプロセッサによって実行されたときに、蓄積デルタ距離データが入手可能な基準受信機／衛星対に対して衛星（ＳＶ）運動およびＳＶクロック訂正キャリア・レートを計算し、他の全てのＳＶ運動およびＳＶクロック訂正キャリア・レートの平均をＳＶ運動およびＳＶクロック訂正キャリア・レートから減算することによって、ＳＶ運動およびＳＶクロック訂正キャリア・レートを基準受信機クロックに対して補償し、サンプル時間期間における基準受信機／衛星対に対する基準受信機トレンド除去ＳＶ運動およびＳＶクロック訂正キャリア・レートを計算するように動作可能なプログラム命令を含む。
【００８５】
[0069]  例１８は、例１７の位相シンチレーション・モニタを含み、少なくとも１つの基準受信機／衛星対に対して現在のサンプル時間における位相シンチレーション監視識別子を計算するプログラム命令が、更に、少なくとも１つのプロセッサによって実行されたときに、数値積分を使用してキャリア位相推定値を計算し、基準受信機／衛星対に対してキャリア位相推定値のサンプル平均を計算し、現在のサンプル時間期間における基準受信機／衛星対に対するキャリア位相推定値のばらつきを計算するように動作可能なプログラム命令を含み、計算されたキャリア位相推定値のばらつきが、位相シンチレーション監視識別子となる。
【００８６】
[0070]  例１９は、位相シンチレーション監視識別子を計算する方法を含む。この方法は、蓄積デルタ距離データが入手可能な基準受信機／衛星対に対して衛星（ＳＶ）運動およびＳＶクロック訂正キャリア・レートを計算するステップと、他の全てのＳＶ運動およびＳＶクロック訂正キャリア・レートの平均をＳＶ運動およびＳＶクロック訂正キャリア・レートから減算することによって、ＳＶ運動およびＳＶクロック訂正キャリア・レートを基準受信機クロックに対して補償するステップとを含む。
【００８７】
[0071]  例２０は、例１９の方法を含み、更に、サンプル時間期間における基準受信機／衛星対に対する基準受信機トレンド除去ＳＶ運動およびＳＶクロック訂正キャリア・レートを計算するステップと、数値積分を使用してキャリア位相推定値を計算するステップと、基準受信機／衛星対に対してキャリア位相推定値のサンプル平均を計算し、現在のサンプル時間期間における基準受信機／衛星対に対するキャリア位相推定値のばらつきを計算するステップとを含み、計算されたキャリア位相推定値のばらつきが、位相シンチレーション監視識別子となる。
【００８８】
[0072]  以上、本明細書では具体的な実施形態について例示し説明したが、同じ目的を達成するために計算される任意の構成を、以上で示した具体的な実施形態と置換してもよいことは、当業者には認められよう。本願は、本発明のあらゆる改変および変形をも権利範囲に含むことを意図している。したがって、本発明は特許請求の範囲およびその均等物によってのみ限定されることを明白に意図している。
【図１】
2015180868.tif 000010
【図２】
2015180868.tif 000011
【図３】
2015180868.tif 000012
【図４】
2015180868.tif 000013
【図５】
2015180868.tif 000014
【外国語明細書】
2015180868000001.pdf