【手続補正書】
【提出日】平成15年8月19日(2003.8.19)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【発明の名称】ＧＰＳ衛星ドリフト監視装置
【特許請求の範囲】
【請求項１】  ＧＰＳ受信機の視野に入り、そして離れる複数の衛星から情報を受信するＧＰＳ受信機と共に使用される装置であって、
視野に入ってくる複数の衛星から前記情報を受信し、時々最小値へ変化する識別子を決定し、複数の前記最小値を用いて所定の信頼性レベルでドリフトエラーの上限を決定する衛星ドリフト監視装置。
【請求項２】  前記ドリフトエラーの上限が、時間的に可変の閾値と比較されて、視野に入ってくる新しい衛星のドリフトを示す選択取消信号を発生する、請求項１の衛星ドリフト監視装置。
【請求項３】  衛星エラー限界ＳＥＬが、ドリフトエラー範囲及び現在の衛星エラー範囲に基き１−ｐmd信頼性レベルで再帰的に計算され、水平保護限界ＨＰＬが解行列Ｔに基きＳＥＬから計算される、請求項１の衛星ドリフト監視装置。
【請求項４】  前記衛星エラー範囲が、式
【式１】
2002520623.tif 000001
に従って計算され、ここで｜ｄn｜は識別子ｄnの絶対値であり、Ｋmdは検出ミス確率ｐmdに相当する統計的シグマ数であり、σSAは選択的有効性ノイズシグマであり、
【式２】
2002520623.tif 000002
は衛星の幾何的配置に依存する係数パラメータの絶対値であり、その上側と下側の添字が同一であり、衛星指標ｎに等しく、前記ドリフトエラー範囲が以前に計算されたドリフトエラー範囲の最小値として計算され、現在観測されたドリフトエラー範囲が、式
ｒ（ｔ）＝ｍｉｎ（ｒ（ｔ−Δｔ），ε／ｔ）
に従って再帰的に決定され、ここでｔは新しい衛星が最初に受信された後の時間であり、Δｔは時間ステップである、請求項２の衛星ドリフト監視装置。
【請求項５】  前記衛星エラー範囲が、式
【式３】
2002520623.tif 000003
に従って計算され、ここで｜ｄn｜は識別子ｄnの絶対値であり、Ｋmdは検出ミス確率ｐmdに相当する統計的シグマ数であり、σSAは選択的有効性ノイズシグマであり、
【式４】
2002520623.tif 000004
は衛星の幾何的配置に依存する係数パラメータの絶対値であり、その上側と下側の添字が同一であり、衛星指標ｎに等しく、前記ドリフトエラー範囲は以前に計算されたドリフトエラー範囲の最小値として計算され、現在観測されたドリフトエラー範囲が、式
ｒ（ｔ）＝ｍｉｎ（ｒ（ｔ−Δｔ），ε／ｔ）
に従って再帰的に決定され、ここでｔは新しい衛星が最初に受信された後の時間であり、Δｔは時間ステップであり、ｒは１−ｐmd信頼性レベルで予測されたドリフト速度であり、ｔは衛星が最初に受信された後の時間である、請求項３の衛星ドリフト監視装置。
【請求項６】  衛星エラー限界（ＳＥＬ）が、式
ＳＥＬ（ｔ）＝ｍｉｎ（ε，ＳＥＬ（ｔ＋Δｔ））
に従って再帰的に決定され、ここでΔｔは時間ステップであり、衛星ｎに関連する水平保護限界（ＨＰＬｎ）が、式
ＨＰＬn＝ｔhnＳＥＬ＋ＫmdｒｆσSA
に従って計算され、ここでＫmdは検出ミス確率ｐmdに相当する統計的シグマ数であり、ｒｆは１より小さな縮小係数であり、
【式５】
2002520623.tif 000005
であって、ｔn1及びｔn2がこの最小二乗解行列Ｔの要素である、請求項５の衛星ドリフト監視装置。
【発明の詳細な説明】
【０００１】
（技術分野）
この発明は航空機ナビゲーション、特に地球測位システム（ＧＰＳ）及び慣性基準システム（ＩＲＳ）を採用し、衛星信号が信頼できなくなるとき、並びに新しく得られた衛星信号が信頼できないときを早期に決定可能なシステムに関する。
【０００２】
（背景技術）
従来ＧＰＳは複数の衛星からの信号を受信することにより航空機の位置を決定するために使用されている。各信号は衛星の位置及び伝送時間に関する情報を有しており、航空機上のＧＰＳ受信機はそれ自体の位置を演算可能である。４個の変数（３軸の位置及び時間）が存在するので、受信機の位置を決定するためには少なくとも４個の衛星からの信号が必要になる。良好な幾何的配置を有する少なくとも５個の衛星が存在する場合、４個の信号からなる各サブセットが位置決めするために利用可能であり、これらの信号の内の１つがエラーであるか否かを決定するために、これらの信号を互いに比較可能である（フェールセーフ）。衛星が少なくとも６個視野内に存在しエラー信号が存在した場合、どの信号がエラーであるかを決定するために、この群の配置を使用可能である（フェールオペレーショナル）。この手順は従来受信機自律完全性監視装置（ＲＡＩＭ）と呼ばれるシステムにより達成された。このシステムについては、バージニア州アレキサンドリアにあるナビゲーション協会（Ｉｎｓｔｉｔｕｔｅ  ｏｆ  Ｎａｖｉｇａｔｉｏｎ）により出版された「グローバルポジショニングシステム」第５巻と題する出版物内のハネウエル・インコーポレーテツドの「ＧＰＳ受信機内のＲＡＩＭ監視装置及び統合ＧＰＳ／ＩＲＳの実現」と題した記事に開示されている。
【０００３】
ＲＡＩＭは４個の衛星ではエラーの有る衛星信号を検出できず、また５個の衛星のみではどの信号がエラーであるかが識別できないので、パイロットはこれらの状況で受信する位置情報を頼ることはできず、ＧＰＳ入力は無視する必要がある。従って４個の衛星のみが視野内に有るとき信頼性のある出力を発生可能なシステムを提供することが望まれている。
【０００４】
１９９６年９月１１日に出願され本発明の譲受人に譲受された第０８／７１２，２３２号「精度不足を検出する解分離装置を有するナビゲーションシステム」と題した特許出願に、この問題の１解決法が開示される。更に詳述するにこの出願には複数の遠隔送信機から信号を受信し、精度不足を検出し、保護限界（即ち最遠の統計的距離であり、これにより装置によって決定された位置が１つの衛星がエラーである場合である可能性が極めて高い）を決定する完全性監視装置及び位置決めアルゴリズムを開示している。この位置決めアルゴリズムは好ましくはカルマンフィルタリング技術を利用し、検出能力を高めるため内部に慣性基準データを内蔵する。このシステムにより、限られた時間の間４個の衛星からの信号を信頼性高く使用できるが、カルマンフィルタバンクにより与えられる保護限界は時々飛行の最も厳しい時期（非精密進入及び境界領域ナビゲーション）に対して必要な保護の限界を超え、従って十分有効でない。更にカルマンフィルタのバンクは慣性基準装置以外の装置内で実施するためには計算上高価であり複雑である。従ってより簡潔で低廉であり、メーカの装置内に容易に内蔵可能で且つ動作の早期の時点で衛星信号エラーを検出可能なシステムの必要性が生じている。
【０００５】
本発明者により本出願と同日付で出願され本発明の譲受人に譲渡された「ＧＰＳ信号欠陥分離監視装置」と題する第０９／１１８，０４６号出願は、ＧＰＳ信号から所定の軸に沿った第１の加速出力を、また慣性基準装置（ＩＲＵ）から同じ所定軸に沿った第２の加速出力を発生する新たな加速監視装置を提供することにより、これらの必要性に対応している。
【０００６】
ＩＲＵ出力は位置を決定するにはそれほど正確ではないが、所定軸に沿った加速の測定は極めて精度が高い。従ってＩＲＵ加速信号は所定軸に沿ったＧＰＳ加速信号と比較される。所定量より大きなエラーが存在する場合、不完全な衛星信号が識別され、計算から排除可能である。従って連続動作には４個の衛星のみが必要である。ランプ（傾斜）出力として考えられる信号ドリフトの場合、ランプの初期変化により、過渡的加速が発生され、加速監視装置により検出されてこの状態の早期の識別が得られる。この特徴は上述の出願に開示されている。しかしながら、ＧＰＳ受信機により観察されるより前にドリフトし始める衛星の信号ドリフトの検出は、過渡的加速がそれ以前に生じているので、新たな加速監視装置によっては検出されない。
【０００７】
（発明の開示）
本発明はこの問題を、衛星エラー範囲（所定の信頼度レベルに関連する範囲又は限界）を決定するために用いる識別子（ｄｉｓｃｒｉｍｉｎａｔｏｒ）を発生し、この衛星エラー範囲をドリフトエラー範囲を決定するために使用し、このドリフトエラー範囲がドリフトが生じたことを示すことにより、新たに捕捉した衛星のドリフトを検出する新規な衛星ドリフト監視装置（時としてゼロクロッシングＲＡＩＭあるいはＺ−ＲＡＩＭとも呼ぶ）を提供することにより解決する。
【０００８】
加速監視装置及び本発明の新規な衛星ドリフト監視装置の両方について以下に詳述するが、加速監視装置は上述した出願で請求され、新規な衛星ドリフト監視装置が本明細書で請求される。
【０００９】
（発明を実施するための最良の形態）
図１を参照すると、ＧＰＳ受信機１２と慣性基準装置１４とを備え、航空機ナビゲーションに使用されるシステム１０が示される。ＧＰＳ受信機１２は当業者には周知なものであり、複数の衛星、例えば図１の衛星Ｓ１〜Ｓ６からそれらの位置及び伝達時間を示す信号を受信するよう動作する。６個の衛星を示したが、ＧＰＳ受信機１２の視野内に６個以上あるいはそれ以下であっても衛星を使用可能であり、その数は時間的期間にわたり変更可能であることは理解されよう。ＧＰＳ受信機１２は衛星Ｓ１〜Ｓ６からの信号で動作して、受信機の位置を決定し、各衛星に対しＧＰＳコードに基く擬似レンジ測定値ｐｒを発生する。また各衛星からの搬送波のサイクル（サイクルの端数を含む）を追跡しカウントすることにより、ＧＰＳ受信機１２は以下の式（１）に従って、時間Δｔにわたる擬似レンジの変化を得るために使用される累積搬送波カウント値ｐｃを与える。
Δｐｒ＝Δｐｃ＝ｐｃ（ｔ＋Δｔ）−ｐｃ（ｔ）      （１）
搬送波カウント値の差分をとることにより得られる擬似レンジの変化は正確であり、センチメートル（０．０１メートル）内にあって且つ擬似レンジ測定値の差分をとることにより形成される擬似レンジ変化は正確で２〜５メートル内である。ＧＰＳコードに基く擬似レンジの変化Δｐｒ、搬送波に基く擬似レンジの変化Δｐｃ、あるいはその組合せ（以下平滑化コード測定値）が本発明で使用される。以下にこの点について説明する。
【００１０】
ｐｒ信号及びｐｃ信号はＧＰＳ受信機１２からの経路１６により示される出力として与えられ、経路１７、１８を介しボックス２０で示す標準ＲＡＩＭへ供給される。標準ＲＡＩＭ２０はまた経路２４を介して高度検出装置２２から気圧高度信号を受け取る。信頼できない衛星信号に遭遇した事が標準ＲＡＩＭ２０により検出された場合、選択取消信号が経路２６により示される出力として標準ＲＡＩＭ２０から与えられる。これは良好な幾何的配置を有する少なくとも６個の衛星が視野に有る場合に生じ、このとき５個の衛星信号のグループはどの衛星信号が不完全かを決定するように比較できる。その衛星信号の内の１つがこのように識別されると、その情報は経路１６から（経路２６及び経路２８を介し）ＧＰＳｐｒ及びｐｃを受け取るボックス３０により示される衛星選択機能へ送られる。衛星選択機能３０は不完全な衛星信号を排除し、経路３２により示される出力として残りの有効な信号をボックス３４により示される位置決定機能へ送る。位置決定機能３４はまた経路２４を介し高度検出装置２２から気圧高度信号も受け取り、指示器若しくは飛行管理システム（図示せず）のような下流の装置への航空機の位置を示す出力を（経路３６により示されるように）発生する。良好な幾何的配置を有する衛星が５個だけ視野に有る場合、４個の衛星信号のグループが依然として少なくとも１つの他の衛星信号と比較でき、何れかの衛星信号が不完全である事を決定可能であるが、不完全な衛星信号の識別ができないので選択取消はできない。このような場合、下流の装置及びパイロットはＧＰＳ信号を使用すべきでないことを通知される（この信号は図１には図示していない）。
【００１１】
上述した出願によれば、経路１６により示される出力信号は経路１７及び経路１８を介し、ボックス４４により示される加速監視装置へ与えられる。加速監視装置４４は慣性基準装置１４（これは当業者には周知であり、複数のジャイロと加速度計からなる）からも入力を受信し、慣性基準装置１４は高度ｈ、高度行列Ｃ、緯度／経度／外れ角行列Ｄ、速度ｖ、及び速度変化Δｖを示すライン４６の出力を発生する。最小３個のジャイロと３個の加速度計が使用されるが、フェールセーフやフェールオペレーショナル動作及び高い信頼性を保証するため、２あるいは３個の冗長システムが使用されることが好ましい。経路４６からのΔｖ信号及び経路１６、１７、１８からのｐｒ信号、ｐｃ信号が本発明の加速監視装置４４により使用され、この動作は図２に沿ってより良好に理解されよう。
【００１２】
図２において、図１に相当する矢印４６により示される経路を経て慣性信号ｈ、Ｃ、Ｄ、ｖ、及びΔｖを受け取り、また図１に相当する矢印４０により示される経路を経てＧＰＳ入力、ｐｒ及びｐｃを受け取る加速監視装置４４が示される。
【００１３】
慣性信号Δｖは速度１０〜６０Ｈｚでフィルタリング処理された速度の増分からなり、「Ｌフレームへ変換」と書かれた機能ボックス５４に与えられる。機能ボックス５４はまた矢印４６で示す姿勢入力Ｃを受け取り、フィルタリング処理された速度増分を局部垂直フレームＬに変換し、この変換された増分を「第２の１Ｈｚ時間差の形成」と書かれた機能ボックス５８に出力する。機能ボックス５８では、各軸に沿った「ダブル位置差信号」（即ち、現在のΔｔ期間での位置変化と以前のΔｔ期間の位置変化との差を示す信号）が慣性加速度（即ち変換された高速速度増分）を積分することにより形成される。機能ボックス５８の出力はＧＰＳ信号（以下に説明する）から引き出される加速度内に存在しない地球の自転による加速度成分を含む、慣性により測定された加速度を反映する。従って機能ボックス５８の出力は望ましくない地球回転成分を除去する「地球回転による加速度の除去」と書かれた機能ボックス６０への入力として与えられる。最後にＧＰＳ信号から引き出される加速度は衛星の方向において（視線に沿って）だけ有効なので、機能ボックス６０からの慣性による基準信号は「視線に沿って投影」と書かれた機能ボックス６２に与えられ、機能ボックス６２はまた矢印６４により示される信号ＬＯＳ（衛星に対する視線に沿っての単位ベクトル）を受け取り、機能ボックス６２の出力は衛星に対する視線に沿って投影された慣性ダブル位置差である。この信号は加速度監視の識別子を形成する「識別子及び時間除去変換」と書かれた機能ボックス６６に与えられる。
【００１４】
ＧＰＳ信号に関し、すべての追跡された衛星に対する擬似レンジ測定値ｐｒ及び累積された搬送波サイクルカウント値ｐｃは経路１８から「範囲差を形成」と書かれた機能ボックス７０に与えられる。機能ボックス７０は初期化時のＧＰＳ位置である矢印７２により示される入力として図１の位置決め機能３４からの初期ＧＰＳ位置ｒinitを受け取る。上述したように、信号ｐｃは累積された搬送波サイクル（各サイクルは通常の衛星信号内の約０．１９ｍ位置変化に相当する）からなり、信号ｐｒはＧＰＳコードに基く擬似レンジ測定値である。信号ｐｒ、ｐｃあるいはその組合せ（平滑化測定値）は衛星の運動を含み、機能ボックス７０は衛星運動成分を除去し、その結果を「第２の１Ｈｚ時間差の形成」と書かれた機能ボックス７４に与えるように動作する。ダブル差信号（即ち現在のΔｔ期間のサイクルカウント値（あるいは平滑化擬似レンジ）の変化と以前のΔｔ期間のサイクルカウント値（あるいは平滑化擬似レンジ）の変化との差）が機能ボックス７４により形成され、（現在位置での視線ベクトル方向での変化及び初期の基準位置における視線ベクトル方向に関連する成分も含むが）ＧＰＳ受信機の加速度を示す機能ボックス７４の出力は「視線補償」と書かれた機能ボックス７６に与えられる。この時点で衛星加速度が現在の位置から見たときと初期の基準位置から見たときが異なるよう投影される。従って機能ボックス７６により形成される出力は視線方向に関連する成分が除去されたものである。機能ボックス７６の出力は「衛星加速度補償」と書かれた機能ボックス７８に与えられる。サイクルカウント値は衛星の運動を含み、機能ボックス７８は衛星運動加速度成分を除去するように動作する。最終の結果は視線に沿ったＧＰＳ信号に基くダブル位置差であり、これは機能ボックス６６に与えられる。「識別子及び時間除去変換」の機能ボックス６６は機能ボックス７８、６２からのＧＰＳ信号及び慣性基準加速度信号で動作し、この２信号の差分を求めることにより識別子を形成する。この機能はまた各衛星特定の識別子からすべての識別子の（すべての衛星に対する）平均を減算して、受信機のクロックオフセットを除去する。この信号は、時間に対し識別子の出力を平均化し固定の閾値と比較して、ライン８２に選択取消信号を発生する機能ボックス８０に与えられ、この選択取消信号は加速度が閾値を超える衛星の選択を取り消すために使用される。上述した機能ボックスにより実行される機能のすべてはコンピュータプログラムにより実行でき、各機能はプログラミング分野における当業者により容易にプログラミングできることは理解されよう。
【００１５】
図１において、図２のライン８２上の選択取消信号がライン２６、２８を介して衛星選択機能３０へ与えられる事が示され、選択取消しされた衛星信号は位置決め機能３４へ送られない。
【００１６】
図１には、経路１６、１７、１８を経てＧＰＳ信号を受け取り、また経路２４を経て高度検出装置２２から気圧高度信号を受け取る受信監視装置（Ａｃｃｅｐｔａｎｃｅ  Ｍｏｎｉｔｏｒ）８６が示される。受信監視装置８６は、例えば衛星擬似レンジが適正な範囲から遠く離間するため明らかに正しくないＧＰＳ信号を検出すべくこの分野で使用される。このようなエラーの有る衛星信号が受信されると、受信監視装置８６は経路２６に選択取消信号を発生し、経路２８を経て衛星選択機能３０に与え、衛星選択機能３０はこのエラーの有る信号の位置決め機能３４による使用を防止するよう動作する。
【００１７】
本発明によれば、新規な「衛星ドリフト監視装置」即ち「Ｚ−ＲＡＩＭ」が経路１６、１７、１８を経て、ＧＰＳ信号を、及び経路２４を経て気圧高度信号を受け取るボックス９２として示される。Ｚ−ＲＡＩＭ９２は水平保護限界ＨＰＬを示す出力９５を経路９５を経てパイロットあるいは飛行管理システム（図示せず）のような下流の航空機装置に与える。
【００１８】
本発明の新規な衛星ドリフト監視装置即ちＺ−ＲＡＩＭ９２は、新たな衛星からの信号が最初に得られるときドリフトが既に開始されているか否かを決定するよう動作する。これについては図３Ａ及び図３Ｂに沿って説明する。
【００１９】
良好な幾何的配置を有する衛星がＮ個存在すると仮定すると、ｎ番目の衛星に対し当業者に周知のＲＡＩＭ識別子が以下の式により形成される。
【式６】
2002520623.tif 000006
ここで、
【式７】
2002520623.tif 000007
は上述した出版物に見られるように、周知な衛星配置に依存する係数であり、Ｎは衛星の数、ｋはｋ番目の衛星を示し、Δｐｒkはｋ番目の衛星に対する測定された擬似レンジ（あるいは平滑化擬似レンジ）と予想擬似レンジとの差である。選択的有効性（ＳＡ、ＤＯＤによりＧＰＳの出力に重畳された意図的なノイズ信号）のため、識別子は図３Ａに見られるように変化する。
【００２０】
１−ｐmdの所定の信頼性レベル（通常９９．９％）を有するよう定義されたエラー範囲εは以下の式により定義される。
【式８】
2002520623.tif 000008
ここで、｜ｄn｜は識別子ｄnの絶対値であり、Ｋmdは検出ミス確率ｐmdに相当する統計的シグマ数、σSAは選択的有効性ノイズシグマであり、
【式９】
2002520623.tif 000009
は式（２）内の衛星配置依存係数の絶対値であり上側と下側の添字は同一であり衛星指標ｎに等しい。
【００２１】
ＳＡノイズの特性に基き、識別子はゼロにクロスするか、時間を置いて即ち少なくとも５分毎に最小値となる。識別子がゼロにクロスするか最小の絶対値｜ｄn｜minに達すると、１−ｐmdエラー範囲の推定のエラーは最小となる。
【式１０】
2002520623.tif 000010
また、１−ｐmd信頼性レベルでの推定されたドリフト速度は次のように表される。
ｒ＝εmin／ｔ        （５）
ここで、ｔは衛星が最初に受信された後の時間である。ｔが増加するに伴い、各最小値｜ｄn｜minに対しより正確なドリフトエラー範囲が得られる。ドリフトエラー範囲が以下の式に従い再帰的に最小値として計算されると、
ｒ（ｔ）＝ｍｉｎ（ｒ（ｔ−Δｔ），ε／ｔ）        （６）
ここで、Δｔは時間ステップであり、その結果は同じであることに留意する必要がある。
【００２２】
衛星エラー限界ＳＥＬは現在の衛星エラー範囲ε及び１−ｐmdドリフトエラー範囲が加えられた、以前に推定された衛星エラー範囲の最小値である。これは以下の式により再帰的に決定される。
ＳＥＬ（ｔ）＝ｍｉｎ（ε，ＳＥＬ（ｔ＋Δｔ））        （７）
ここで、Δｔは時間ステップである。
【００２３】
衛星ｎと関連する水平保護限界ＨＰＬ（航空電子工学装置に於ける所要出力）は以下のように計算できる。
ＨＰＬn＝ｔhnＳＥＬ＋ＫmdｒｆσSA        （８）
ここで、Ｋmdは検出ミス確率ｐmdに相当する統計的シグマ数であり、ｒｆは１より小さな縮小係数で、
【式１１】
2002520623.tif 000011
は当業者には周知な最小二乗解行列Ｔであり、ｔn1及びｔn2はこの行列Ｔの要素である。初期のドリフトを有する衛星を排除するアクションが講じられない場合でも、ドリフトが一定である、即ちフェールセーフ動作を保証する限りにおいて、この衛星と関連する水平保護限界ＨＰＬnがこのエラーを常に制限する。ドリフトが一定でない場合、上述した出願の加速監視装置４４は加速度を検出する必要がある。速度ｒが所定の時間依存しきい値を越えるときは、フェールセーフ動作能力を達成するために経路２２、２３を経て新たな衛星の選択を解除し、かつ図１の機能ブロック２４を実行する必要がある。既に位置決めに使用されている衛星は加速監視装置４４により監視されているので、新しく捕捉された衛星だけこのように監視する必要がある。
【００２４】
ここで示すボックス及び経路は好ましい実施形態の動作を説明する目的のために使用されたが、これらのボックス及び経路はソウフトウエア若しくはハードウエア、あるいはその組合せにより実現可能であることは理解されよう。
【００２５】
従って新しく捕捉された衛星が捕捉前にドリフトを開始するときを検出可能なドリフト監視装置が提供されたことが理解されよう。多くの別の実施形態が当業者において推測でき、ここで開示される本発明は本発明の精神から離れることなく他の特定形態で実施可能である。従ってここで説明した実施形態はすべての点で例示のためと考えるべきであり、これに制限されるものではなく、本発明の範囲は上述の説明によるものに限られるものではなく添付の請求項により表わされ、請求項の等価物の意味及び範囲内に入るすべての設計変更が本発明に含まれる。
【図面の簡単な説明】
【図１】
図１は本発明を利用するシステムの概略ブロック図である。
【図２】
図２は上述した出願の加速監視装置の動作を説明するためのブロック図である。
【図３Ａ】
図３Ａは本発明の新規な衛星ドリフト監視装置の動作を示すグラフである。
【図３Ｂ】
図３Ｂは本発明の新規な衛星ドリフト監視装置の動作を示すグラフである。