JP2000512380A

JP2000512380A - 未知の信号源の位置確定

Info

Publication number: JP2000512380A
Application number: JP09512454A
Authority: JP
Inventors: ディヴィッドパトリックハワース
Original assignee: UK Secretary of State for Defence
Current assignee: UK Secretary of State for Defence
Priority date: 1995-09-20
Filing date: 1995-09-20
Publication date: 2000-09-19
Anticipated expiration: 2015-09-20
Also published as: EP0852017A1; AU708274B2; KR100319993B1; DE69511118T2; CA2232396A1; GB9805264D0; CA2232396C; JP3556952B2; KR19990063678A; WO1997011383A1; EP0852017B1; DE69511118D1; US6018312A; GB2321356A; GB2321356B; GR3031185T3; ES2133806T3; BR9510640A; AU3612195A; MX9802116A

Abstract

(57)【要約】未知の信号を人工衛星中継器に伝達する未知の信号源の位置を確定する方法が、各受信器（１８）で中継器からの信号を受信することを含む。受信器（１８）は共通の信号源（２２）から各中継器を介して基準信号を受信する。各受信器（１８）によって受信された未知信号及び基準信号は、別の所で受信された信号とは独立し、且つそれらのタイミング及び位相情報を保存するように一緒に処理される。信号は周波数下げ変換され、デジタル化され、共通の処理コンピュータ（１５０）に送られる。コンピュータ（１５０）は基準信号の交差あいまい関数を達成して、差分時間オフセット（ＤＴＯ）及び差分周波数オフセット（ＤＦＯ）を決定する。同様のことを受信された未知信号に対して達成し、相対ＤＴＯ及びＤＦＯを決定する。相対時間及び周波数シフトが、基準信号ＤＴＯ及びＤＦＯを組み込むことによって、未知の信号間に導入され、基準信号ＤＦＯに対して、未知信号ＤＦＯが決定されることが可能とされる。基準信号は、未知信号の位相ノイズ及び周波数ドリフトの影響を相殺するためにも使用される。未知の信号ＤＴＯ及びＤＦＯから、未知信号源（１０）の位置が計算される。

Description

【発明の詳細な説明】未知の信号源の位置確定本発明は、複数の信号中継器によって受信される未知の信号源の位置を確定する方法と装置に関する。従来技術の確定システムとしては、確定されるべき信号源、即ち未知の送信器からの２つの信号の到達時間差（ＴＤＯ）の確定に基ずいた技術を採用したものが知られている。この技術では、送信器からの信号は、各経路において地球の軌道を回る各々の人工衛星によって２つの独立した信号経路に沿って地上基地の受信ステーションに中継される。各人工衛星は、地上の送信器からの信号（アップリンク）を一の周波数で受信して折り返し発振器の助けを受けて地上の受信器へ周波数シフトされた信号（ダウンリンク）を戻すトランスポンダーシステムを組み込んでいる。信号経路は長さが同じではなく、それは２つの人工衛星が異なる位置にあるためである。その結果、一方の経路に関連した受信ステーションにおける信号到達時間は、他方の経路に関連したものとは異なっている。２つの人工衛星の位置と受信ステーションの位置は既知であり、結果的に２つの異なる経路を経た信号の到達時間差（ＴＤＯＡ）は未知の送信器の位置情報を提供することになる。特定の人工衛星の位置に対して、一定のＴＤＯＡの点の軌跡は、位置ライン（ＬＯＰ）と称されるカーブを地上に描く。ＴＤＯＡの特定値を発生する未知の地上基地の送信器はその値に関連したＬＯＰ上のある点で確定される。ＴＤＯＡの正しい値を決定するために、人工衛星間の幾何学的関係は、人工衛星トランスポンダーや受信ステーション機器における信号伝播遅れと共に既知でなければならない。然し乍ら、ＬＯＰはただ曲線上の任意の点に対して未知の送信器を特定するだけである。米国特許第５，００８，６７９号は、前述のように２つの中継人工衛星を組み込んだ送信器位置確定システムを開示しており、このシステムは、未知の送信器からの信号のＴＤＯＡと、それらの到達周波数差（ＦＤＯＡ）と称されるものの両方を採用している。ＦＤＯＡは、中継人工衛星が地球に対して運動しているので生じ、それぞれの場合でその運動は、ドップラーシフトをダウンリンクの信号周波数に持ち込む。２つの人工衛星のドップラー効果は異なり、これは、各々の人工衛星からの受信ステーションに到達する信号間に周波数差やＦＤＯＡを発生させる。ＴＤＯＡのアプローチにおけるように、ＦＤＯＡ測定は、未知の送信器が存在している位置線（ＬＯＰ）を確定する。ＴＤＯＡとＦＤＯＡの位置線は、一般的に異なり、それらの交点は未知の送信器の位置を決める。然し乍ら、米国特許第５，００８，６７９号のシステムは、幾つかの拘束を受ける。それは、人工衛星の位置及び速度を高い精度で知る必要がある。もし人工衛星が地球の赤道面に対して相対的に０．１度より大きな軌道の傾きを有していれば、それは、５０ｋＨｚ以下の帯域幅の通信信号とは有効に使用され得ない。このシステムは、更に地上ステーションの局部発振器と人工衛星の信号折り返し発振器から非常に高い精度の位相の安定性を要求している。その後者の要求は、非常に低い、即ち０．０１度より小さい傾斜角度の人工衛星にはそのシステムを無効にするかも知れない；人工衛星の運動によって起こされるドップラーシフトは小さく且つそのような傾斜では測定が困難であるからである。更に、米国特許第５，００８，６７９号のシステムは、共通時間と周波数基準源を持つ共通配置された受信地上ステーションに使用できるだけである。それは、受信地上ステーションが幾何学的に離れている場合には変更無しでは使用できない。そのシステムは、更に位相の固定された周波数の下方変換鎖の使用を必要としている。送信器の位置確定システムの代わりの形態を提供することが、本発明の目的である。本発明は、複数の信号中継器によって受信された未知の信号源の位置確定方法を提供するもので、その方法は、次の工程を有している：−即ち、（ａ）各信号中継器から未知の信号を受信する複数の受信器を配置し；（ｂ）単一信号の複製か、又は周波数、時間及び位相において単一信号に対してロックされた信号の複製のいずれかであって、公知の場所の基準送信手段から信号中継器に送信される各基準信号を各信号中継器から受信する受信器を各々配置し；（ｃ）各受信器によって受信された各未知の信号と基準信号とを、これらの信号が他で受信された信号から独立して、相互にそれらのタイミングと位相の情報を保存するようにまとめて処理し；（ｄ）処理された基準信号と処理された未知信号の交差あいまい関数演算を実施し、次の（ｉ）未知信号の差分時間オフセット（ＤＴＯ）と差分周波数オフセット（ＤＦＯ）の値、（ii）異なった信号中継器位置に対応した未知信号のＤＴＯの値、（iii）異なった信号中継器位置に対応した未知信号のＤＦＯの値、（iv）信号中継器の異なった組合せに対応した未知信号のＤＴＯの値、（v）信号中継器の異なった組合せに対応した未知信号のＤＦＯの値の内の少なくとも一つを確定するために未知信号における位相ノイズと周波数ドリフトの作用を打ち消すように処理された基準信号を採用し；（ｅ）工程（ｄ）で確定されるＤＴＯ及び／若しくはＤＦＯの値から未知の信号源の位置を算定する上記工程を有している。本発明は、従来の技術で要求されていた精度で人工衛星の位置と速度を知るのを要求せず、また少なくとも３度迄の軌道の傾斜を持つ人工衛星に使用できると言う長所を与えてくれる。本発明は、従来の技術で要求された中継や受信器位相の安定度を必要としない。それは、地理学上共に位置決めされ、共通した時間と基準周波数を持ち、共通した精度の周波数標準を用いた位相の固定された周波数下方変換を採用した受信器を必要としない。本発明は、更に、他の欲しくない信号が存在していても未知の信号源の位置確定を実施できる。本発明の工程（ｃ）の処理は、好ましくは未知信号と基準信号とを所定の帯域幅を有した中間周波数（ＩＦ）信号に下げ変換することによって、またサンプル採取タイミングと下げ変換とを正確な周波数とタイミング基準に従って制御してそれらのディジタルサンプルを得ることによって各受信器で受信された信号について実施される。未知信号と基準信号は、好ましくは４Ｍｈｚより大きくないＩＦ帯域幅を有するように構成される。ＤＴＯの確定の為に、ＩＦ帯域幅は、好ましくは、基準信号の帯域幅も好ましくは合わせられる未知信号の帯域幅にできるだけ近く設定される。正確なＤＦＯの確定の為に、ＩＦ帯域幅は、好ましくは１００ＫＨｚより小さい、例えば１０ＫＨｚと小さい。これは、ＤＦＯ処理の精度を改善する。好適な実施例では、本発明の工程（ｄ）における交差あいまい関数演算は、ヒルベルト（Ｈｉｂｅｒｔ）変換手法から複素数データを発生する工程を有する。これは、そのような処理に対して未知信号と基準信号の同相及び直角成分とを発生させる必要を無くする。本発明の工程（ｄ）での交差あいまい関数演算は、好ましくは更に基準信号のＤＴＯを確定し、また基準信号の次のような処理工程を有している：−即ち、（ａ）関連した試行周波数オフセットが必要な仮の値であることを示す関数最大値が得られるまで試行周波数オフセットの範囲に対して交差あいまい関数の数値を求めることによって基準信号の仮の値を見い出し；（ｂ）各受信器に関連した第１と第２の基準信号をそれらの周波数変域等値に変換し；（ｃ）周波数シフトがＤＦＯの仮の値から構成されていて、周波数変域において第２基準信号に対する第１基準信号を周波数シフトし；（ｄ）周波数変域基準信号のいずれもの欲しくない周波数成分も削除し；（ｅ）周波数変域において、周波数成分結果を造り出すために第１基準信号の各周波数成分の複素数共役に第２基準信号の対応した周波数成分を掛け；（ｆ）周波数成分結果を時間変域に変換し、第１と第２の基準信号間の相対時間オフセットの値の各範囲に対して交差あいまい関数の各値を発生させ；（ｇ）交差あいまい関数の一組の最大強度値を選択し、それらの間の補間法によって交差あいまい関数の強度の最大値を得て、その最大値に対応するもので、第１と第２の基準信号間の相対時間オフセットとしての基準信号のＤＴＯを導き出す上記工程を有している。交差あいまい関数の演算による基準信号ＤＴＯの計測は、（計算によってそれを得るのとは反対に）未知源の確定精度を改善することになる。本発明の工程（ｄ）での処理は、次の工程を有した方法によって基準信号のＤＦＯを得る工程を含むことかできる：− （ａ）各受信器に関連した第１と第２の基準信号間のＤＴＯを確定し；（ｂ）基準信号のＤＴＯに等しい第１と第２の基準信号間の相対時間シフトを導入し；（ｃ）相対時間シフト後に基準信号のサンプリングを行い；（ｄ）時間成分結果を発生させるために、各第１基準信号サンプルの複素数共役に第２基準信号の対応したサンプルを掛け；（ｅ）時間成分結果を周波数変域にフーリエ変換し、第１と第２の基準信号間の相対周波数オフセットの値の各範囲に対して交差あいまい関数の各値を発生し；（ｆ）交差あいまい関数の一組の最大強度値を選択し、それらの間の補間法によって交差あいまい関数の強度の最大値を得て、その最大強度値に対応するもので、第１と第２の基準信号間の相対時間オフセットとしての基準信号ＤＦＯを導き出す上記工程を含むことができる。本発明の工程（ｄ）での処理は、次の工程によって未知信号のＤＴＯを確定する工程を含むことができる：− （ａ）時間と周波数のシフトは各々基準信号のＤＴＯとＤＦＯに等しいものであって、各受信器に関連した第１と第２の未知信号間に相対時間と周波数のシフトを導入し；（ｂ）相対時間と周波数のシフト後に未知信号をサンプリングし；（ｃ）第１と第２の未知信号をそれらの周波数変域等価物に変換し；（ｄ）試行周波数オフセットを変換された第２未知信号に適用し；（ｅ）変換された未知信号におけるいずれもの欲しくない周波数成分も削除し；（ｆ）周波数成分結果を発生させるために、第１の未知信号の各周波数成分の複素数共役に第２の未知信号の対応した周波数成分を掛け；（ｇ）第１と第２の未知信号間の相対時間オフセットの値の各範囲に対して交差あいまい関数の各値を発生させるために周波数成分結果を時間変域に変換し；（ｈ）交差あいまい関数の値を得るために或る範囲の試行周波数オフセットに対して工程（ｄ）から（ｇ）を繰り返し；（ｉ）交差あいまい関数の一組の最大強度値を選択し、それらの間の補間法によって交差あいまい関数の強度の最大値を得て、その最大強度値に対応するもので、第１と第２の未知信号間の相対時間オフセットとしての基準信号ＤＴＯを導き出す上記工程による。本発明の工程（ｄ）での処理は、未知信号のＤＦＯを確定することを含み、また基準信号処理と次の工程を組み込むことから構成される段階を含むことができる：− （ａ）シフトが基準信号のＤＴＯとＤＦＯに各々等しいもので、各受信器に関連した第１と第２の基準信号間に相対時間と周波数のシフトを導入し；（ｂ）相対時間と周波数シフト後に基準信号をサンプリングし；（ｃ）時間成分結果を発生させめために、各第１基準信号サンプルの複素数共役に第２信号の対応したサンプルを掛け；（ｄ）第１と第２の基準信号の間の相対周波数オフセットの値の各範囲に対して交差あいまい関数の各値を発生させるために時間成分結果を周波数変数にフーリエ変換し；（ｅ）交差あいまい関数からの欲しくない周波数変数成分を削除し；（ｆ）位相ノイズの乱れ情報を含んだ第１と第２の基準信号間に一組の瀘過された時間変域結果を発生させるために削除された交差あいまい関数を時間変域にフーリエ変換で戻す上記工程を組み込むことから構成できる。本発明の工程（ｄ）での処理は、未知信号のＤＦＯを確定することを含み、また基準信号処理と次の工程を組み込むことから構成される別の段階を含むことができる：− （ａ）基準信号のＤＴＯとＤＦＯに各々等しい第１と第２の未知信号間に相対時間と周波数のシフトを導入し；（ｂ）相対時間と周波数のシフト後に未知信号をサンプリングし；（ｃ）時間成分結果を発生させるために、各第１未知信号サンプルの複素数共役に第２信号の対応したサンプルを掛け；（ｄ）各時間成分結果に第１と第２の基準信号の瀘過された時間成分結果の複素数共役を掛け；（ｅ）時間成分結果を周波数変域にフーリエ変換し、第１と第２の未知信号間の相対周波数オフセットの或る範囲の値の各々に対して交差あいまい関数の各値を発生させ；（ｆ）交差あいまい関数の一組の最大強度値を選択し、それらの間の補間法によって交差あいまい関数の最大値を得て、その最大値に相当している第１と第２の未知信号間の相対周波数オフセットとして基準信号のＤＦＯに対する未知信号のＤＦＯを導き出す上記工程を組み込んだ別の段階を含むことができる。未知源の位置は、好ましくは次に基ずいた方法によって算定される：− （ａ）各ケースでの差分傾斜範囲（ＤＳＲ）が、各送信器から信号中継器までの経路の長さの差であって、未知信号ＤＴＯと基準信号ＤＴＯとから基準送信器と未知送信器のＤＳＲを含んだ式を得、（ｂ）各ケースでの差分傾斜範囲割合（ＤＳＲＲ）が各ＤＳＲの変化割合であって、未知信号のＤＦＯと基準信号のＤＦＯから基準送信器と未知送信器のＤＳＲＲを含んだ式を得、（ｃ）ＤＳＲ及びＤＳＲＲの式のテーラー展開分析を行って、未知及び基準のＤＴＯとＤＦＯから未知送信器の位置を導出し、中継器の位置及び速度と基準送信器及び受信器の位置を導出することに基づいている。工程（ｃ）において未知源に対して導出される位置は、その第１の値とされ、その一つ上の改善された値は、テーラー展開分析における第１の値の戻し代入によって、また収斂が得られるまで未知源の位置の導出を反復して導出できる。未知源の位置の片寄りは、次の工程によってより正確にできる：− （ａ）未知送信器の位置、中継器の位置及び速度並びに受信器の位置の導出された値から、未知信号のＤＴＯ及びＤＦＯの算定値を発生させ；（ｂ）未知信号のＤＴＯ及びＤＦＯの算定値と、交差あいまい処理の助けで計測されたものとの間の差から導出された未知送信器の位置に対して補正を発生させる工程による。第１に、未知信号のＤＳＲとＤＳＲＲの算定された値と計測された値との間の差は、２つの一次方程式によって未知源の位置における経度と緯度の誤差に関連付けられている。経度と緯度の誤差は、好ましくはこれらの式を解くことで確定される。基準信号は、信号中継器から得られる通信信号としたり、代わりに、互いに位相コーヒレンスで信号を発生する送信器から得られよう。送信器は、各信号中継器に向けられている。処理された基準信号及び処理された未知信号は、好ましくは交差あいまい関数処理に先立って記憶される。代わりの局面では、本発明は複数の信号中継器によって受信された未知信号源を確定する装置を提供し、その装置は：− （ａ）基準信号が単一信号の複製か又は単一信号に対して周波数、時間及び位相においてロックされた信号の複製のいずれかであり、基準信号が公知位置の基準送信手段から信号中継器に伝達されるものであって、各信号中継器から未知信号と各基準信号を受信する複数の受信器と；（ｂ）各受信器によって受信された各未知信号と基準信号を、これら信号が他で受信された信号から独立して相互にそれらのタイミングと位相の情報を保存するようにまとめて処理する手段と；（ｃ）（i）処理された基準信号と処理された未知信号について交差あいまい関数処理を実行し且つ、次の：（１）未知信号の差分時間オフセット（ＤＴＯ）と差分周波数オフセット（ＤＦＯ）の値と、（２）異なった信号中継器位置に対応した未知信号のＤＴＯ値と、（３）異なった信号中継器位置に対応した未知信号のＤＦＯ値と、（４）信号中継器の異なった組合わせに対応した未知信号のＤＴＯ値と、（５）信号中継器の異なった組合わせに対応した未知信号のＤＦＯ値の内の少なくとも一つを確定するために未知信号における位相ノイズや周波数ドリフト作用を相殺する基準信号を採用し、（ii）交差あいまい関数処理で確定されるＤＴＯ及び／若しくはＤＦＯの値から未知信号源の位置を算定する手段とを有している。本発明の装置は、ＤＴＯ及び／若しくはＤＦＯに対して得られた値に従って適切に時間及び／若しくは周波数において受信された信号を相殺したり、そのような相殺信号からそれらの別の値を得るように構成される。それは、ほぼ等しい強さの基準信号を信号中継器に伝達する構成の基準信号伝達手段を有し、これら信号は、変調されたまた変調されていない搬送波成分とすることができる。それは、更に交差あいまい関数演算に先立って、処理された基準信号と処理された未知信号を記憶する手段を有することができる。本発明がより完全に理解されるようにする為に、その実施例を今添付図面を参照にして例示のためだけに説明する。そこで：− 図１は、地上基地の送信器と、人工衛星中継器と、地上基地受信器との間の信号伝播を示しており；図２は、関連した地上基地送信器と人工衛星中継器を共の有した本発明の送信器位置確定システムの概略線図であり；図３は、人工衛星中継器の回路をより詳細に示しており；図４は、図２の送信器位置確定システムの一部である受信器と組み合わせた信号処理システムの回路の詳細をより詳細に示しており；図５は、図２の送信器位置確定システムに組り込まれたコンピュータ演算装置をを示しており；図６は、異なる人工衛星中継器から受信された信号の時間オフセットに対する相関のグラフであり：図７は、異なる時間オフセットの確定を改善するための補間法を図解しており；図８は、異なる人工衛星中継器から受信された信号の周波数オフセットに対する相関関係のグラフであり；図９及び図１０は、未知と基準の送信器に対する周波数変域相関を図解しており；図１１は、欲しく無い周波数成分の除去を図解しており；図１２は、相対ＤＦＯを確定するための交差あいまい関数の相関を図解したグラフであり；図１３は、全地球的位置確定手順の幾何学的構成を図解した概略図であり；図１４は、２つの中継人工衛星によって受信された基準信号のレベルを平衡させるために設計された基準送信器を図解した概略図であり；また図１５は、ＤＦＯとＤＴＯの測定を容易にする波形特性を持つ基準送信器の代替形の概略図である。図１を参照すると、アメリカ合衆国１１に存在した未知の送信器１０が地球１２の表面上に示されており、その北半球は、北極（図示されていない）を中央に配置して図示されている。未知送信器１０は、静止軌道の第１人工衛星１４に主要な放射強さのローブ（図示されていない）を向けている。それは、第１アップリンクｌ₁ ^u沿って人工衛星に伝播し且つ人工衛星を使っている未知の信号との干渉を起こす信号を送信している。未知信号の周波数は、人工衛星の未知のチャンネルを日常的に監視するスペクトル分析装置によって確定される。Ｋｕ帯域（１１−１４ＧＨｚ）で作動する一般的な通信衛星は、各々３６Ｍｈｚ幅で各１００通信信号を搬送できる１６チャンネルを有している。送信器１０は、更に放射型サイドローブ（図示されていない）を静止軌道の第２人工衛星１６に向けており、その信号は第２アップリンクｌ₂ ^uに沿ってその衛星に伝播している。経路の符号ｌ₁ ^uとｌ₂ ^uに対するスーパスクリプトの『ｕ』は、未知送信器１０に由来することを示している。第１人工衛星１４は、未知送信器１０からの信号を受信して、その衛星に向けられてイスラエルに配置された第１地上基地ステーションや受信器１８Ａに第１ダウンリンクｌ₁ ^mに沿ってそれを再送信している。第２人工衛星１６は、更に未知送信器の信号を受信して、それを南アメリカ２１に配置された第２地上基地受信器１８Ｂに第２ダウンリンクｌ₂ ^mに沿って再送信している。ここで、スーパスクリプト『ｍ』は、人工衛星を監視している地上基地受信器への経路を示している。地上基地の受信器１８Ａ、１８Ｂは、図示されたいずれかの又は両方に対して区別せずに、特定の場合に適するように１８Ａ又は１８Ｂのように符号１８によって参照される。送信器１０から第１受信器１８Ａへの全信号伝播経路長は、経路ｌ₁ ^uとｌ₁ ^m長さの合計に等しく、また送信器１０から第２受信器１８Ａへのそれは、経路ｌ₂ ^u とｌ₂ ^mの長さの合計に等しい。アフリカ２３の既知の地上位置の基準送信器２２は、基準信号を第１と第２の人工衛星１４、１６に各々第３と第４のアップリンクｌ₁ ^rとｌ₂ ^rに沿って送信するが、ここでスーパスブリクトの『ｒ』は、基準送信器２２からの送信を示している。基準送信器２２は、人工衛星１４と１６の一方に関連した通信チャンネルを使っているものから選択される。人工衛星１４と１６は、基準信号を受信器１８へダウンリンク経路ｌ₁ ^mとｌ₂ ^mに沿って各々再送信する。今図２を参照すると、本発明の送信器確定システムが、概略的に一般に３０で示されている。未知送信器１０と基準送信器２２と受信器１８とは、アンテナのシンボルで示されている。人工衛星１４と１６は、矩形によって示されている。受信器１８Ａと１８は、第１と第２の捕捉システム３２Ａと３２Ｂに各々接続されており、それらの各々は、後でより詳細に説明する別々のチャンネルで未知と基準の信号を処理する。捕捉システム３２は、離れた処理現場の中央制御・演算コンピュータ（図示されていない）に各最新のデータリンク３６Ａと３６Ｂによって接続されている。人工衛星１４と１６の回路は、図３に示されている。各々、受信（アップリンク）アンテナ５２と送信（ダウンリンク）アンテナ５４とが搭載されているコンテナ５０から構成されている。受信アンテナ５２は、低ノイズ増幅器５６に接続され、次いで、周波数折り返し発振器６０からの局部発振器入力を受信するミキサー５８に接続されている。局部発振器周波数は、両人工衛星１４と１６に対して１．５ＧＨｚである。ミキサー５８は、結果的に１．５ＧＨｚの周波数シフトダウンを発生する。ミキサー５８からの出力は、帯域幅フィルター６２に至り、その後に、送信アンテナ５４に信号供給を行う出力増幅器６４に達する。さて、更に図４を参照すると、各捕捉システム３２の回路がより詳細に示されている。各捕捉システムは、全地球的測位システム（ＧＰＳ）の受信器１００から構成されており、タイミング信号を送るためにそれを一つ以上のＧＰＳ人工衛星（図示されていない）に接続するアンテナ１０２を備えている。ＧＰＳは、宇宙に展開された幾つかの人工衛星から構成されており、それらから信号が入手される。ＧＰＳ受信器１００は、各々タイミング（ｔ）と周波数（ｆｒ）の信号の為の出力部１０６と１０８と共に制御入力部１０４を有している。出力部１０６は、実際は、後で説明する捕捉システム３２の各チャンネルに各々接続された２つの出力部を有している。各々２つのチャンネルを有した２つの捕捉システム３２が設けられており、結果的に、信号サンプリングが開始される異なった開始時間Ｔを各々有した４つのチャンネルが設けられることになる。２つの受信器１８Ａと１８Ｂに関連したタイミングと周波数の信号は、非常に類似しているが、然し、必ずしも同一である必要はなく、ｔ_Aとｆｒ_A、ｔ_Bとｆｒ_Bとなっている。これは、未知と基準の各送信器が、ＧＰＳの異なった部分にアクセスすることになるように地球表面上で遠く離れて配置される為である。結果的に、受信器１８Ａの信号は、受信器１８Ｂの信号と位相が揃っておらず、本発明はそのような一致を必要としないのが特長である。ＧＰＳ受信器１００の制御入力部１０４は、それに制御信号を送る局地ホストコンピュータに接続されている。周波数信号ｆｒは５ＭＨｚである。タイミング信号ｔは、後でより詳細に説明されるように、未知の送信器を確定する作業における信号サンプリングを制御する。周波数信号ｆｒのように、それは、ＧＰＳから受信する信号からＧＰＳ受信器１００によって発生される。未知送信器を確定する作業を開始するために、コンピュータ１０５は、開始時間を示す指令制御入力部１０４に送り；この時間が生じたことをＧＰＳが示すと、ＧＰＳ受信器１００は、隣接パルスが一定の時間差△ｔを有している一連のパルスのように、タイミング信号の発生を開始する。タイミング間隔△tは、両受信器１８Ａと１８Ｂにおいて同じである。コンピュータ１０５は、ｔ_o＋ｊ△ｔからタイミング信号に応答して採取されたいずれかの信号サンプルの時間を得る。但し、ｔ_oは開始時間で、ｊはサンプル番号である。受信器１８からの出力信号は、低ノイズ増幅器１１０に至り、次に、発振器１１４からの局部発振器入力信号を受信するミキサー１１２に達する。発振器１１４は、１１６でＧＰＳ受信器出力部１０８に接続され、１１．８０５ＧＨｚの周波数ｆｒに位相固定される。ミキサー１１２からの出力信号は、各々未知と基準の信号の２つのチャンネル１２０Ｕと１２０Ｒに至る。これらのチャンネルは、各ケースで関連したチャンネルを示すためにサフィクスＵ又はＲで同様に参照されている同じ要素を有している。チャンネル及び要素は、区別せずにいずれかを又は両方を示す為にＵ又はＲのサフィクス無しに、また特定が必要な時に関連サフィクスで参照されることになる。各チャンネル１２０では、ミキサー１１２からの信号は、調節可能な予選択フィルター１２２に至り、次いで、発振器１２６から６９８．６ＭＨｚの周波数で局地発振入力信号を受信するミキサー１２４に達する。発振器１２６は、１２８でＧＰＳ受信器の出力部１０８に接続されており、後者の周波数ｆｒに位相固定されている。予選択フィルター１２２の中央周波数と帯域幅及び発振器１２６の周波数は、局地ホストコンピュータ１０５の制御の下で整調される。ミキサー出力信号は、一定の中央周波数と、局地ホストコンピュータの制御の下で整調可能な帯域幅とを有した予選択フィルター１３０に至り、次いで可変ゲイン増幅器１３２に達し、そこから次に出力信号はアナログ−ディジタル変換器（ＡＤＣ）１３４に至る。ＡＤＣ１３４は、高速で高安定の８ビット装置である。それは、タイミング入力部１３６をＧＰＳ受信器出力部１０６に接続しており、そこからそれはタイミング信号ｔを受信する。タイミング信号の各パルスを受信すると、ＡＤＣ１３４は、可変ゲイン増幅器１３２から出力信号のディジタル化されたサンプルを発生する。信号のサンプリング割合は、出力信号の帯域幅の２倍と最小であり、局地ホストコンピュータの制御を受けている。ＡＤＣ１３４は、メモリ１３７に接続されている。未知と基準のチャンネル１２０Ｕと１２０Ｒのメモリ１３６Ｕと１３６Ｒは、各々両方共局地ホストコンピュータ１０５に接続されており、それは次に公記録保管部１４０に、上述のようにデータリンク３６に、接続部１４２によってＧＰＳ受信器制御入力部１０４に、接続部１４４によってＧＰＳ受信器出力部１０６に接続されている。さて更に図５を参照すると、遠隔地の処理現場３４のエレメントがより詳細に示されている。現場３４は、データリンク３６に、またＧＰＳシステムと通信するアンテナ１５４を備えた第３のＧＰＳ受信器１５２に接続された中央制御演算コンピュータ１５０を組込んでいる。コンピュータ１５０は、更にＤＳＰユニット１５６と公記録保管部１５８にも接続されている。送信器確定システム３０は、次のように作動する。未知送信器１０は、第１人工衛星１４の通信チャンネルの信号と干渉を起こす信号を送信する。未知信号の周波数は、人工衛星の通信チャンネルを監視しているスペクトル分析装置によって確定される。未知信号は人工衛星１４と１６に伝播し、そこでそれはミキサー５８によって１．５ＧＨｚだけ周波数がシフトダウンされ、第１と第２の受信器１８Ａと１８Ｂに各々再度送信される。基準信号が、そこで人の介入によって選択される。それは、第１人工衛星１４の通信チャンネルに存在して、また第２人工衛星１６にサイドローブを向けた送信器で発生して、好ましくは人工衛星１４のダウンリンクの監視から確定されるように未知信号のものと同じ帯域幅を有しているいずれかの未知信号である。それは、未知信号のものと十分に異なった周波数を有しており、これら信号がミキサー１１２における周波数の下方転換後に異なったチャンネルに分けられるようにしている。一例として、一般的な未知信号が１４．００５Ｇｈｚの中央周波数で送信され、１２８ｋｂ／ｓのデータから構成されている。この信号は、人工衛星の折り返し発振器６０によって１２．５０５Ｇｈｚに周波数がシフトダウンされる。隣接した信号は、未知信号よりも周波数が約１０Ｍｈｚ高いチャンネルで識別される２５６ｋｂ／ｓデータ信号のような人工衛星１４のダウンリンクのスペクトルを監視することによって、基準として選択される。そのような基準信号は、１４．０１５Ｇｈｚの送信器周波数に対応して１２．５１５Ｇｈｚの周波数を有することになろう。基準信号は、人工衛星１４と１６によって各受信器にリレーされる。未知送信器１０の主ローブ用の目標である第１人工衛星１４における信号対ノイズの割合は、単位元よりも大幅に大きくなりそうであり、５〜１５ｄＢの一般的な値を有している。しかし、第２人工衛星１６は、未知信号１０のサイドローブから低出力信号を受信するだけなので非常に低い信号対ノイズの割合を有した信号に関連付けられそうでである。そのような低い信号レベルは、従来の手段では検出できないし、また後で説明する信号の相関技術を使用する必要がある。受信器１８が受信した後は、未知と基準の信号は、１１０において増幅され、１１．８０５Ｇｈｚの局部発振器周波数と１１２においてミキシングされる。局部発振器周波数は、それらの各々の関連した未知または基準の周波数との差が７００Ｍｈｚの所定の中間周波数（ＩＦ）に接近しているように各局地ホストコンピュータによって整調される。ミキサー１１２におけるミキシングによって、次に未知と基準の信号を各予選択フィルター１２２Ｕと１２２Ｒに至るＩＦ信号に変換する。予選択フィルター１２２は、局地ホストコンピュータによって整調される帯域幅を有している。最初の組の信号データに対して、未知チャンネルの予選択フィルター１２２Ｕは、整調され、帯域幅を未知信号ＩＦ上に集中させている。未知チャンネルの後選択フィルター１３０は、下方変換信号の帯域幅を設定する。広い帯域幅は、他の誤差がより重要になる点まで計測時間における誤差を低減し、これで４ＭＨｚの限度を設定している。基準フィルター１２２Ｒは、基準信号ＩＦ上に帯域幅を集中させるように整調される。未知と基準のチャンネルフィルター１２２は、未知チャンネルのフィルター１２２Ｕが基準信号を拒絶し、基準チャンネルのフィルター１２２Ｒが未知信号を拒絶するのに適した帯域幅と周波数選択性とを有している。フィルター掛け後は、未知と基準のＩＦ信号は、ミキサー１２４においてより低い中間周波数（ＬＩＦ）に下方変換される。局部発振器１２６の周波数は、基準信号に対する未知信号の位相と周波数とが各捕捉システム３２に保存されるようにＧＰＳ信号に対して正確に位相固定される。これら周波数は、局地ホストコンピュータ１０５の制御を受け、それで一定のＬＩＦが後選択フィルター１３０に至り、そこで信号は引き続きもう一度フィルターに通されて最終的にフィルター掛けされた帯域幅を定義し、それは未知信号に対して４Ｍｈｚにも及ぶことができる。この後、ＬＩＦ信号は、ＡＤＣ１３４（８ビット）の全ダイナミックレインジを使うために増幅器１３２のゲインを適切に設定することで振幅が調節される。ＡＤＣ１３４による信号サンプリングは、次のように開始される。中央コンピュータ１５０は、局地ホストコンピュータ１０５の各々に開始時間を示し、コンピュータ１０５はそれを各ＧＰＳ受信器１００に中継器する。ＧＰＳが、開始時間が生じた旨示すると、各ＧＰＳの受信器１００は、タイミング信号ｔを開始する。基準信号の場所は既知なので、２つの人工衛星経路を介する伝播遅れは算定され、また２つの現場間の捕捉開始時間オフセットも異なった伝播遅れを考慮して算定される。捕捉は、０．００１秒のタイミング精度で実行される。上述のように、タイミング信号ｔは、１．９５３１２５μ秒の一連の一定サンプリング時間間隔△ｔの一連のタイミングパルスから構成されている。パルスは、ＧＰＳ周波数ｆｒに正確に位相固定されており、従って更に関連した捕捉システム３２における局部発振器１１４と１２６の周波数にも固定されている。各ＡＤＣ１３４は、各タイミング信号に応答して（チャンネルに応じて）未知又は基準の信号のディジタル信号サンプルを発生する。各メモリ１３６は、関連した開始時間と共に各ディジタル信号サンプルを一時的に記憶する。各局地コンピュータ１０５は、引き続いて、各サンプルとそれに関連したメモリ１３６Ｕと１３６Ｒからサンプリングする開始時間とから成るデータを読み取り、それらをその公記録保管部１４０に記憶する。未知の送信器位置を個々に確定する際、タイミング信号が切られる以前に、合計で１６．３８４×１０６のサンプルが４つのＡＤＣ１３４ＵＡ、１３４ＲＡ、１３４ＵＢ及び１３４ＲＢの各々によって採取される。いずれかのディジタル信号サンプルが採取される時間は、t０＋ｊ△ｔから得られるが、その場合ｔ0は開始時間であり、ｊはサンプル番号である。上述したように、ＤＡＣ１３４当り一つずつ、ｔ_OUA、ｔ_ORA、ｔ_OUB及びｔ_ORBによって与えられた４つまでの異なった開始時間が用意され、そこで時間は協定世界時（ＵＴＣ）に対して決められている。サンプリングが完了した後、（各々第１と第２の受信器１８Ａと１８Ｂに関連した）公記録保管部１４０Ａと１４０Ｂは、各々未知と基準の両送信器１０と２２に対するサンプルと開始時間を収容する。更に、各個々の受信器１８Ａ又は１８Ｂでは、ミキサー１１２と１２４及びＡＤＣ１３４が、局部発振器と、ＧＰＳの周波数と時間の信号ｆｒとｔに対して位相固定されたタイミング信号とを採用しているために、未知と基準の信号は、下方変換されて可干渉的にサンプル採用される。しかし、受信現場は、ＧＰＳの異なった部分へアクセスすることになるように遠く離れて地球表面上に位置しているので、ｆｒ、ｔ及びｔ_oは、それらが受信器１８Ｂの場合のように受信器１８Ａでは全く同じではない。公記録保管部１４０Ａと１４０Ｂにおけるディジタル信号のサンプルとそれらのタイミングの記憶は、都合が好いかもしれないずれかの不特定時間にサンプルされたデータから未知送信器の位置確定が得られようにするものである。ディジタル信号サンプルは、ディジタル処理のためにデータリンク３６Ａと３６Ｂに沿って２つの受信地から中央の制御演算コンピュータ１５０に転送される。その演算コンピュータによって実施される第１処理作業は、一定の（８ビット）正確な数値から浮動小数点数値にディジタルサンブルを変換するものである。これによって、引き続き処理を受けて信号品質の劣化を減じることになる。ディジタル信号サンプルに対して演算コンピュータ１５０によって実施される第２の処理作業は、それらを実数から複素数の形に変換するものである。この変換には、ヒルベルト（Ｈｉｌｂｅｒｔ）変換手法が採用されている。これを図解するために、ゼロ位相角を有した真正サイン波形を考慮することである。同相と直角成分のの２次元プロットで見ると、サイン波形は同相軸線に平行な直線として現れる。軸が直角成分と同相成分と時間となっている３次元プロットでは、サイン波形は、ゼロの矩象要素位置で時間／同相成分平面においてサイン波形として現れる。ヒルベルト変換後には、サイン波形は、回転する位相ベクトルに変換される。直角成分と同相成分と時間の軸を備えた３次元プロットでは、位相ベクトルは、軸方向において時間軸に平行な円形螺線として現れる。位相ベクトルは、直角成分と同相成分の両方を有しているので、ヒルベルト変換が信号を実数から複素数形に変換することが証明される。演算コンピュータ１５０によって実行される処理作業を先ず概略説明し、引き続いて数学的処理を提供する。演算コンピュータ１５０は、未知と基準の送信器からの信号を処理する交差あいまい関数（ＣＡＦ）を実行し、これら信号間の差分時間オフセット（ＤＴＯ）と差分周波数オフセット（ＤＦＯ）を確定する。ＤＴＯとＤＦＯは、次のように定義される：− ＤＴＯ：異なったルートを介して受信後に、元々同じ信号の２つの複製の受信間における差分時間オフセット、又は時間遅れ。ＤＦＯ：異なったルートを介して受信後に、元々同じ信号の２つの複製の受信間における差分周波数オフセット、又は相対周波数シフト。ＤＴＯと相対ＤＦＯから、基準送信器に対する未知送信器の位置は本発明に依って確定される。次の理論的説明において、『信号』の表現は、演算コンピュータ１５０におけるヒルベルト変換から浮動小数点への変換後にコンピュータ１０５Ａと１０Ｂからのディジタル信号サンプルを意味するものとして構成されている。演算コンピュータ１５０によって実行される次の処理作業は、基準信号の処理をする交差あいまい関数（ＣＡＦ）処理を実行して、各受信器１８Ａと１８Ｂと捕捉システム３２Ａと３２Ｂを介して得られた基準信号の２つの複製の間のＤＴＯとＤＦＯを確定するものである。ＣＡＦ処理は、１９８１年６月のＩＥＥＥ会報のＡＳＳＰ−２９Ｎｏ．３『あいまい関数処理のアルゴリズム』のタイトルでＳステイン（Ｓｔｅｉｎ）氏による文献に説明されている。交差あいまい関数又はＣＡＦＡ（τ、υ）は、次の式で定義される。その場合、s₁とs₂は、２つのアナログ信号であり、s₁ ^*はs₁の複素共役を表し；τとυは、各々ｓ₁に対するｓ₂に与えられる時間と周波数オフセットである。式（１）は、相互関係プラス周波数シフトの操作を表している。Ａ（τ、υ）は、一般に同相成分と直角成分の両方を有した複素数の形である。Ａ（τ、υ）の絶対値は、ｓ₁（ｔ）がs₂（τ＋υ）ｅｘｐ（−ｉ２πυｔ）と同一の時に最大である。ｓ₁とｓ₂が共通の源から伝播した後の同じ信号の複製である場合に、また伝播中に、これら信号の一方が他方に対して時間遅れを起こして周波数シフトを起こしており、Ａ（τ、υ）の絶対値は、適用されたオフセットτとυが伝播の途中で生じた時間遅れと周波数シフトを正確に打ち消すと最大となっている。中継人工衛星を介して受信された信号に対して、時間遅れと周波数シフトは、信号経路と中継人工衛星の運動の間の差及び人工衛星の折り返し発振器における差に依るものである。演算コンピュータ１５０は、後で説明するように式（１）のディジタル等値を実行する。第１例では、それは第１受信器１８Ａを介して受信されたものに対して第２受信器１８Ｂを介して受信されたその基準信号複製に一連の試行周波数シフトを加える。各試行周波数オフセットυに対して演算コンピュータ１５０は、或る範囲の時間オフセットτの値を適用し、｜Ａ（τ、υ）｜の数値をディジタルで求める。それは、｜Ａ（τ、υ）｜の識別可能な最大値に対するこの手続きの結果を探求し、またこの最大値が対応しているτとυの値は、受信器１８Ａと１８Ｂによって受信された基準信号の２つの複製間のＤＴＯとＤＦＯに各々等しい。この手続きで得られたＤＴＯの値は、本発明の目的に対して許容可能な精度を有しているが、しかしＤＦＯのそれは、仮の概算値であり、『粗いＤＦＯ』と称される。グラフで表されているように、離散値として表された時間オフセットτを変える関数としての信号s₁とｓ₂間の相互関係は、数個の個別の時間オフセットに渡った一連のスパイクとなっており、その場合各スパイクの寸法は、相関の度合を表している。しかし、相関度合を最大にする必要な時間オフセットは、一般に２つの適用された時間オフセットの間に入っている。必要な時間オフセットを確定するために、３つの最大相関が採用され、放物曲線は、τの関数としてプロットされたそれらの大きさの対数に適合されている。最大相関の点を示している放物線のピークに対応している補間された時間オフセットは、必要な時間オフセットとして、必要なＤＴＯとして採用される。処理の次の段階において、コンピュータ１５０は、第１と第２の受信器１８Ａと１８Ｂを介して受信された未知の信号の複製を使用する。それは、基準の粗いＤＦＯに等しい周波数オフセットを、第１受信器１８Ａを介して受信されたものに対して第２受信器１８Ｂを介して受信された未知信号の複製に適用する。これの理由は、周波数シフトの主要因が人工衛星の折り返し発振器における差から、また一方の人工衛星の他方に対する運動から生じるので、基準と未知のＤＦＯが同じく成り、またこれからの周波数シフトが未知と基準の両方の信号にとって同じであるためである。この関係で、隣の人工衛星上の折り返し発振器は一般に１Ｈｚより小さい。未知信号の複製を使って、演算コンピュータ１５０は、次に基準信号の複製に以前使用されたＣＡＦ処理手続きを反復する。この手続きの結果は、受信器１８Ａと１８Ｂによって受信された未知信号の２つの複製間にＤＴＯと粗いＤＦＯを与える。この場合、未知の粗いＤＦＯは、基準の粗いＤＦＯに対して確定される。未知の粗いＤＦＯの絶対値は、次いでその相対値を基準の粗いＤＦＯに加えることによって確定される。この点で、全ての基準と未知の信号が１００ＫＨｚの狭い帯域幅に１２２においてフィルター掛けされ、周波数オフセットが未知と基準の信号に対する粗いＤＦＯから確定されるように適用され、４つのＡＤＣ１３４ＵＡ、１３４ＲＡ、１３４ＵＢ、１３４ＲＢが異なった開始時間ｔ_OUA、ｔ_ORA、ｔ_OUB、ｔ_ORBを有するように開始時間が調節されることを除外して、新しい組の基準と未知の信号は前述のように得られる。演算コンピュータ１５０は、局地ホストコンピュータ１０５Ａと１０５Ｂに指令して、未知信号の粗いＤＦＯの絶対値に等しいオフセットを、第２受信器地１８Ｂにおける未知チャンネル１２０ＵＢの局部発振器１２６ＵＢの周波数に適用させる。このオフセットは、第１受信器現場１８Ａにおける未知チャンネル１２０ＵＡの局部発振器１２６ＵＡの周波数に相対している。同様に、基準の粗いＤＦＯに等しいオフセットは、第２受信器現場における基準チャンネル１２０ＲＢの局部発振器１２６ＲＢの周波数に適用される。ここで、再度そのオフセットは、第１受信器現場における基準チャンネル１２０ＲＡの局部発振器１２６ＲＡの周波数に相対している。この手続きの目的は、人工衛星１４との運動によって誘発される周波数オフセットを修正するものである。未知チャンネルの後選択フィルター１３０ＵＡと１３０ＵＢは、演算コンピュータ１５０からの指令に応答して、各々局地ホストコンピュータ１５０Ａと１５０Ｂによって１００ＫＨｚの帯域幅に戻される。開始時間の間の関係は、次のようになっている：ｔ_OUB−ｔ_OUA＝未知信号ＤＴＯ（アップリンク＋ダウンリンク）（2a）ｔ_ORB−ｔ_ORA＝基準信号ＤＴＯ（アップリンク＋ダウンリンク）（2b）式（２ａ）と（２ｂ）は、開始時間が各々がアップリンクとダウンリンクの信号経路の組合せから生じる未知信号ＤＴＯと基準信号ＤＴＯに各々等しいオフセットを採用することを示している。人工衛星１４と１６から受信された信号は、下方変換とディジタル信号サンプリングを受け、以前説明したように明白に規定された時間及び周波数のオフセットの影響を受ける新しい組の基準及び未知の信号データを提供する。そのように発生されたディジタル信号は、コンピュータ１５０においてＣＡＦ処理を受け、基準信号に対して未知信号の新しいより正確な値を発生する。この値と前に得られたＤＴＯの値は、引き続き（後で説明するように）コンピュータ１５０によって処理されて次の量：即ち異なった傾斜範囲（ＤＳＲ）と、異なった傾斜範囲割合（ＤＳＲＲ）と各々ｄ−ＤＳＲとｄ−ＤＳＲＲによって示されるＤＳＲとＤＳＲＲの各々における変化とを発生する。これらは、次のように定義される：− ＤＳＲ：地上の点から２つの人工衛星への経路の長さにおける差；ｄ−ＤＳＲ：或る期間に渡るＤＳＲにおける全体の変化；ＤＳＲＲ：時間に対する変化割合；ｄ−ＤＤＳＲ：或る期間に渡るＤＤＳＲにおける全体の変化；量ＤＳＲ、ｄ−ＤＳＲ、ＤＳＲＲ及びＤＤＳＲから、既知送信器２２対する未知送信器１０の位置が計算によって得られる。ＤＳＲとＤＳＲＲの組合せを使用した演算コンピュータ１５０によるＣＡＦ処理と未知送信器位置についての理論的分析についてより詳細に説明する量ＤＳＲ、ｄ−ＤＳＲ、ＤＳＲＲ及びＤＤＳＲから、既知送信器２２対する未知送信器１０の位置が計算によって得られる。ＤＳＲとＤＳＲＲの組合せを使用した演算コンピュータ１５０によるＣＡＦ処理と未知送信器位置についての理論的分析についてより詳細に説明するＤＴＯの確定の為に交差あいまい関数やＣＡＦＡ（τ、υ）の使用を参考にする場合、式（１）に従ってＣＡＦ処理を受けることになっている２つの信号s₁とｓ₂を考慮すること。もしｓ₁とｓ₂が時間変域アナログ信号ｓ₁（ｔ）とｓ₂（ｔ）であれば、次に一定間隔Δｔでのディジタルサンプリングの後にそれらはｓ₁ （ｊΔｔ）とｓ₂（ｊΔｔ）に成り、そこでｊはサンプリング番号である。Δｔは定数なので、ｊはただ変数であり、信号の表現はｓ₁（ｊ）とｓ₂（ｊ）に相当する。同様に、もしｓ１（ｊ）とｓ₂（ｊ）が周波数変域フーリエ変換Ｓ₁（ｆ）とＳ₂（ｆ）を有していれば、分離周波数間隔Δｆでの後者のディジタル表現は、各々Ｓ₁（ｋΔｆ）とＳ₂（ｋΔｆ）であり、Ｓ₁（ｋ）とＳ₂（ｋ）に相当しており、そこでｋは、周波数番号である。即ち、信号ｓ₁（ｔ）とそのフーリエ変換Ｓ₁（ｆ）のサンプリングされた表現は：− s₁(t)=s₁(jΔt)≡s₁(j)，S₁(kΔf)≡S₁(k)，ｊ=0，1，…,N-1，k=0,1，…,N-1，且つΔfΔt=1/N (3) 同様な式は、サブスクリプトのインデックストの１を２に変えることでｓ₂（ｔ）とＳ₂（ｆ）に対して得られる。信号と変換ｓ₁、ｓ₂、Ｓ₁とＳ₂は、複素数であり、即ち、それらは同相成分と直角成分を有している。分離した時間と周波数の変域表現（例えば、ｓ₁、Ｓ₁）は、次の分離フーリエ変換（ＤＦＴ）とこの変換の逆によって連結されている：− ＤＦＴは、例えば、１９９５年のＭａｔｈＣｏｍｐｕｔａｔｉｏｎ社発行のクールリー（Ｃｏｏｌｅｙ）ＪＷとトゥーキー（Ｔｕｋｅｙ）ＪＷの著『複素数フーリエシリーズの機械計算のアルゴリズム』の第１９巻の２９７−３０１ページに説明されているような早いフーリエ変換アルゴリズムを使って実行される。一般に、ｓ₂は、ｓ₁に対して周波数シフトを経験することになろう。良好な相関がs₁とｓ₂の間で達成されるようにこのシフトを補償する為には、周波数オフセットδfがｓ2(ｊ)に適用され、必然的にそのフーリエ変換Ｓ2(ｋ)の表現への変化を伴って、それをＳ₂ ^'(ｊ)に変換する。Ｓ₂ ^'(ｊ)は、次の式（６）から（８）によって定義される。式（８）において、δｋはδｆ／Δｆに等しい。式（６）から（８）は、ｓ₂がδｆによって周波数相殺され、Ｓ₂の対応シフトがδｋによってされていることを示している。サーチ処理中の演算速度を早める為に、粗いＤＦＯが、周波数変域において直接数単位の周波数段だけＳ₂を相殺することによって実行され、それで演算される必要があるフーリエ変換の数を減らしている。フーリエ変換Ｓ₁(ｋ)とＳ₂(ｋ)が発生されると、演算コンピュータは、後続の演算でノイズ汚染された、さもなければ欲しく無いその中の周波数成分を除去する。以前に示したように、一連の試行周波数オフセットを使用して、また各周波数オフセットに対して或る範囲の時間オフセットに渡ってＣＡＦの数値を求めることで大凡のＤＦＯを見出す必要がある；これは、関連試行周波数オフセットが必要なＤＦＯであることを示しているＣＡＦ大きさにおける有効最大値が得られるまで続行される。次の演算において、説明を簡単にする為に、このことが既に行われていると仮定しよう；即ち、ｓ₂は、ＣＡＦ大きさ｜Ａ（τ、υ）｜の最大値を発生するものに少なくとも等しいδｆの値だけ周波数相殺されていると仮定される；ＤＴＯ値を得る目的でＡ（τ、υ）の数値を求める為に、演算コンピュータは、２つのフーリエ変換Ｓ₁とＳ₂の結果を算定する。これは、式（９）で示されている：− 式（９）の逆変換は、式（１０）と（１１）に示されている：− Ｓ₁ ^*(k)とＳ₂(ｋ)を代入して:− 但し、１_tauは、演算コンピュータによって適用される時間間隔Δｔの離散数である。即ち、τ＝１_tauΔｔ式（１１）におけるｋに対する合計は、δ項を含まない等比級数として実行される。かくして、即ち：− Ｎを式（１１）に代入し、ｊとｊ’についての合計は、ただｊ’＝ｊ＋１_tau だけが（１３）から許容されるので、ｊについての合計に還元されることを喚起して：− 式（１４）は、正常化係数Ｎの項においてディジタルで再実行されている。前に説明してように、今ＤＴＯは、Ａ（１_tau）の大きさを最大にする時間オフセット数１_tauなので演算コンピュータ１５０によって確定される。演算コンピュータ１５０は有限時間分解を行う、即ち、それは、１_tauの有限段Δｔに時間オフセットπをを適用することである。Ａ（１_tau）の大きさを最大にするオフセットは、一般に２つの個々の時間オフセット間に存在している。コンピュータ１５０は、Ａ（１_tau）の３つの最大連続値に対応した３つの時間オフセット１_tau を確定し、それら値を対数スケールに変換し、次いで対数値の放物線補間法を実行して推定ピーク値を見出す。ピーク値が生じる時間オフセットは、必要なＤＴＯ値である。ＤＴＯを確定する手続きは、次のように要約される。（ａ）一連の試行周波数オフセットを使用し、それらの各々に対する或る範囲の１_tauの値に対してＡ（１_tau）の数値を求めることによって概略のＤＦＯを見出し；これは、関連試行周波数オフセットが必要なＤＦＯ値であることを示しているＡ（１_tau）の大きさにおける有効最大値が得られるまで続行される。（ｂ）概略のＤＦＯによってｓ₁に対してｓ₂を周波数シフトし；（ｃ）ｓ₁とｓ₂の双方をそれらの周波数変域等価Ｓ₁とＳ₂に変換し；（ｄ）どの欲しく無い周波数成分も除去し; （ｅ）Ｓ₂の各周波数成分にＳ₁の同じ周波数成分の複素数共役を掛け；（ｆ）周波数成分結果を時間変域に戻し変換し；（ｇ）（対数スケールにおいて）Ａ（１_tau）の３つの最大強度の値間で補間して、Ａ（１_tau）を最大にする時間オフセットを確定する。図６は、ＤＴＯを確定する為の時間オフセットの関数としての時間変域相関の又はＡ（１_tau）のグラフである。それは、約＋０．００３秒で最大２００を示している。図７は、適用された時間オフセットπに対してプロットされているＡ（１_tau）の最大値（２０ｌｏｇ｜Ａ｜として表現されている）についての補間法を示している。２０ｌｏｇ｜Ａ｜の内の３つの値は、それらが適合される放物線２１６と共に、２１０、２１２と２１４で示されている。２０ｌｏｇ｜Ａ｜のピーク値は、鎖線２１８によって示されており、そのτ軸上での交点はＤＴＯを与える。ＤＴＯをより正確に確定する為に、時間変域におけるｓ₁ ^*とｓ₂の結果が形成される。かくして、相対遅れを含めてｓ₂を有した式（５）と（６）から：− フーリエ変換式（１５）は、次のものを与える但し、１_nuは分離周波数オフセット数を示しいる、即ち、υ＝1_upsilon Δｆ（１３）に従って：− 式（１８）は、周波数変域交差相関を示している。ＤＦＯは、Ａ（1_upsilon）として表現された相関の大きさを最大にする周波数オフセットを見出すことによって確定される。従前通り、３つの最大相関の対数の放物線補間法が、ピーク相関とＤＦＯの推定値とを確定する為に採用される。ＤＦＯを確定する手順は、次のように要約される：− （ａ）粗いＤＦＯによってｓ₁に対してｓ₂を周波数シフトし；（ｂ）ＤＴＯによってｓ₁に対してｓ₂を時間シフトし；（ｃ）比較的狭い帯域におけるフィルター掛け後にｓ₁と時間シフトされたｓ2 とをサンプル採取し；（ｄ）各時間要素ｓ₂（ｊ）に同じ時間要素ｓ₁（ｊ）の複素数共役を掛け、各時間成分結果を発生し；（ｅ）時間成分結果を周波数変域にフーリエ変換し；（ｆ）Ａ（１_upsilon）の３つの最大の強度値の間で補間し、相関を最大にする周波数オフセットであるＤＦＯを確定する。周波数変域における相関応答は、サイドローブを有しているので、周波数スロット間における補間に誤差が存在する。理想的には、補間法は、ｓｉｎｘ／ｘの補間関数を使用して、リニアーな信号大きさについて補正結果を得るべきである。対数の信号大きさについて簡単な放物線補間法が使用されるので、劣化した補間精度に成っていまう。この補間法の問題を克服するために、時間変域信号は周波数変域への変換前にハミング（Ｈａｍｍｉｎｇ）窓で窓付けされる。この窓付け機能で、無視できるレベルに相関応答してサイドローブを減らし、かくして３つの最大成分の大きさの対数について放物線補間法を使用する正確な補間ができるようにする。図８は、ＤＦＯの確定の為の周波数オフセットνの関数として周波数変域相関の、又はＡ（１_upsilon）のグラフである。それは、約＋８０Ｈｚのυの値において最大２４０を示している。実際には、ＤＦＯの計測の精度には限度がある。これらの限度は、２つの人工衛星１４と１５に搭載されている折り返し発振器６０に存在している位相ノイズと、アンテナ送信器の幾何学的形状における変化と、地上ステーション受信器１８Ａと１８Ｂ上の位相ノイズによって設定される。図９と１０は、各々未知送信器１０と基準送信器２２のための周波数（即ち、周波数変域において）に対するＣＡＦのグラフである。これらのデータは、実際の信号と人工衛星から得られる。以前に説明したピーク補間技術が役立たなくなるようになんら単一ピークが存在していないことが分かる。図９と１０は、各々ピーク２５０（−９．２８１Ｈｚ）と２５２（−５９．２７６Ｈｚ）を有した同様な構成となっている。図９と１０における相関を洗練する手続きは、次のようになっている。基準信号のＣＡＦは、図１１に示されているように相関ピークから成分を除去する為に削り取られる。その目的は、逆変換された信号上のノイズを減らすことである。欲しく無い成分の削り取り後には、基準ＣＡＦは、時間変域に逆フーリエ変換される。基準送信器２２のものに対する未知送信器１０のＤＦＯを得るために、それらの時間変域波形は、共に掛けられる。これは、基準信号に対するフィルタ −掛けされた時間変域結果の対応成分の複素数共役が掛けられた未知信号に対する式（１５）によって特徴付けされた時間変域結果の各成分を含んでいる。この掛け算の結果は、次いで周波数変域に戻し変換される。結果として生じた相対ＣＡＦは、図１２に示されており、そこでは相関大きさが周波数に対してプロットされており、ＣＡＦは−０．０３８Ｈｚにおいて最大２７０を有している。図９と１０の相関スパイクにおける劣化は、図１２において実質的に除去されている。この技術は、受信器１８を介して得られた信号における位相ノイズ劣化を打ち消す。それは、周囲の人工衛星の折り返し発振器や受信器の作用によるかなりの位相ノイズ劣化が存在している状況下で本発明が使用され得るようにするものである。これは、未知の信号上の位相ノイズが、同じ受信器１８で受信された基準信号上のものと相関され、適当な処理による低減を受けやすいという理由による。更に、本発明は、地球の赤道面に対して０．５度より大きな軌道傾斜角度を有した静止軌道の人工衛星１４と１６にとって効果的である。そのような傾斜角度は、時間についてのＤＦＯに有効な変形を生じさせる。さて位相補償技術についてより詳細に説明する。式（１５）の左手側の結果の中の信号ｓ₁（ｊ）とｓ₂（ｊ）を考慮すること。これら信号は、ＤＦＯの粗い値によって周波数シフトされ、ＤＴＯの正確な値によって遅延される（未知又は基準のいずれかのＤＦＯは、どちらの信号が処理されているかに左右される）。ＤＦＯの周波数シフトは、ただＤＦＯの実際値に対する近似値でしかないので、信号は、残留未補償の周波数シフトを含んでいる。更に、その状況は、信号上に位相乱れを持ち込む。最後に、残留周波数シフトは、時間に左右される。これらの要因を考慮して、信号は記録される：− 但し：− ｓ₁ ^U（ｔ）は、経路ｌ₁ ^uとｌ₁ ^mに沿って人工衛星を介して受信された未知信号であり；ｓ₂ ^U（ｔ）は、経路ｌ₂ ^uとｌ₂ ^mに沿って人工衛星を介して受信された未知信号であり；ｓ₁ ^R（ｔ）は、経路ｌ₁ ^rとｌ₁ ^mに沿って人工衛星を介して受信された基準信号であり；ｓ₂ ^R（ｔ）は、経路ｌ₂ ^rとｌ₂ ^mに沿って人工衛星を介して受信された基準信号であり； ν₁ ^U（ｔ）は、経路ｌ₁ ^uとｌ₁ ^mを介した未知信号上の残留周波数シフトであり； ν₂ ^U（ｔ）は、経路ｌ₂ ^uとｌ₂ ^mを介した未知信号上の残留周波数シフトであり； ν₁ ^R（ｔ）は、経路ｌ₁ ^rとｌ₁ ^mを介した未知信号上の残留周波数シフトであり ν₂ ^R（ｔ）は、経路ｌ₂ ^rとｌ₂ ^mを介した未知信号上の残留周波数シフであり； φ₁ ^U,M（ｔ）は、経路ｌ₁ ^uとｌ₁ ^mを介した未知信号上の残留相乱れであり； φ₂ ^U,M（ｔ）は、経路ｌ₂ ^uとｌ₂ ^mを介した未知信号上の残留相乱れであり； φ₁ ^R,M（ｔ）は、経路ｌ₁ ^rとｌ₁ ^m介した未知信号上の残留位相乱れであり； φ₂ ^R,M（ｔ）は、経路ｌ₂ ^rとｌ₂ ^m介した未知信号上の残留位相乱れである。位相乱れは、次のように成分に分解する：− φ₂ ^U,M=φ₂ ^M,a+φ₂ ^T+φ₂ ^U,a+φ₂ ^U+φ^U φ₂ ^R,M=φ₂ ^M,a+φ₂ ^T+φ₂ ^R,a+φ₂ ^R+φ^R （20） φ₁ ^U,M=φ₁ ^M,a+φ₁ ^T+φ₁ ^u,a+φ₁ ^U+φ^U φ₁ ^R,M=φ1^M,a+φ₁ ^T+φ₁ ^R,a+φ₁ ^R+φ^R 但し； φ₁ ^M.aは、第１受信器１８Ａと第１衛星１４の間の大気に因る乱れであり； φ₂ ^M,aは、第２受信器１８Ｂと第２衛星１６の間の大気に因る乱れであり； φ₁ ^Tは、第１衛星１４上の折り返し発振器に因る乱れであり； φ₂ ^Tは、第２衛星１６上の折り返し発振器に因る乱れであり； φ₁ ^U,aは、未知送信器１０と第１衛星１４の間の大気に因る乱れであり； φ₂ ^U,bは、未知送信器１０と第２衛星１６の間の大気に因る乱れであり； φ₁ ^R,aは、基準送信器２２と第１衛星１４の間の大気に因る乱れであり； φ₁ ^Uは、未知送信器１０から第１衛星１４に経路ｌ₁ ^uを介して伝搬された未知信号の時間に関する残留周波数変化に因る位相乱れであり； φ₁ ^Rは、基準送信器２２から第１衛星１４に経路ｌ₁ ^rを介して伝搬された基準信号の時間に関する残留周波数変化に因る位相乱れであり； φ₂ ^Uは、未知送信器１０から第２衛星１６に経路ｌ₂ ^rを介して伝搬された未知信号についての時間に関する残留周波数変化に因る位相乱れであり； φ₂ ^Rは、基準送信器２２から第２衛星１６に経路ｌ₂ ^rを介して伝搬された基準信号についての時間に関する残留周波数変化に因る位相乱れであり； φ^Uは、未知信号上の残留位相であり； φ^Rは、基準信号上の残留位相である。式（19）から我々は次の結果を得る：− 但し φ＝｛２π〔（ν₂ ^U−ν₁ ^U）−（ν₂ ^R−ν₁ ^R）〕ｔ＋（φ₂ ^U,M −φ₁ ^U,M）−（φ₂ ^R,M−φ₁ ^R,M）｝（22）最後に式（２０）から：− φ＝｛２π〔（ν₂ ^U−ν₁ ^U）−（ν₂ ^R−ν₁ ^R）〕ｔ＋（φ₂ ^U−φ₁ ^U）−（φ₂ ^P−φ₁ ^P） −（φ₂ ^U,a−φ₁ ^U,a）−（φ₂ ^R,a−φ₁ ^R,a）｝（23）式（２１）に従った処理の後に残っている残留位相成分は、基準信号のものに対する未知信号の残留ＤＦＯと付加位相補正との組み合わせであることが式（２３）から理解されよう。これらの位相補正は、時間に対するＤＦＯの変化と、未知と基準の送信器１０と２２から２つの人工衛星１４と１６へのアップリンクに対する大気の異なった作用とによって生じる。干渉を被っている第１人工衛星１４から軌道で３度離れている第２人工衛星１６に対して、大気の影響は、式（２ 1）による処理が図１２に図解されているように基準信号のものに対して未知信号の残留ＤＦＯを発生するように実質的に無効になる。従って、本発明に従って基準信号を使用した処理は、位相ノイズと残留周波数ドリフトに対して補償する。前者は、０．０５度より小さい人工衛星軌道傾斜値において特に重要であり、後者は、０．５度より大きい人工衛星軌道傾斜値において特に重要である。この位相補償技術は、信号帯域幅が小さく（＜１００ｋＨｚ）且つ処理ゲインが大きい（＞６０ｄＢ）場合に未知信号に対するＤＴＯの確定にも適用される。大きい処理ゲインと狭い帯域幅が信号に対して、波形持続時間は、人工衛星の折り返し発振器の一般的な可干渉性時間である数秒を越えており、また結果的に相関強度は低下される。位相補償技術を適用することで、完全な相関強度が達成される。さて未知送信器の位置確定についてより詳細に説明する。第１段階は、所謂全地球的位置解答を得ることである。もし未知と基準のＤＴＯが各々ＤＴＯ（未知）とＤＴＯ（基準）であって、各々対応した送信器／人工衛星の分離の程度と方向とをそこで幾何学的に表しているベクトルとして図１で未知と基準のアップリンク経路インデックスｌ₁ ^U、ｌ₂ ^U、ｌ₁ ^rとｌ₂ ^rを扱うと：− 但し、ｌ₂₁（ｒ）は未知送信器からのＤＳＲであり、ｌ₂₁（ｒ₀）は基準送信器からのＤＳＲであり、ｃは光の速度である。更に、ｒとｒ₀は各場合における原点としての地球の中心からの未知送信器１０と基準送信器２２の位置ベクトルである。人工衛星を経た経路を移動する信号に対するドップラー周波数シフトは、アップリンク周波数におけるアップリンクＤＳＲＲとダウンリンク周波数におけるダウンリンクＤＳＲＲによるものである。アップリンクとダウンリンクの周波数は、人工衛星の折り返し発振器によって持ち込まれた周波数シフトによって異なる。かくして：− 但し、ｆ^Tは、人工衛星の折り返し発振器の周波数である。（２６）から（２５）を減算し、両側にｃ／ｆ^uを掛けることで得るものは：− 但し、ν₂₁（ｒ）は、未知送信器１０に対するＤＳＲＲであり； ν₂₁（ｒ₀）は、基準送信器２２に対するＤＳＲＲであり； ν₁（ｒ_m1）は、第１人工衛星１４から第１受信器１８Ａへの傾斜範囲率（ＳＲＲ）であり； ν（ｒ_m2）は、第２人工衛星１６から第２受信器１８ＢへのＳＲＲであり；ｆ^rは、最初に送信された基準信号の周波数であり；ｆ^uは、最初に送信された未知信号の周波数である。もし受信器１８Ａと１８Ｂが同じ側にあれば、式（２７）の右手側の少なくとも２つの項は、その側でＤＳＲＲとなる。式（２４）と（２７）は、基準送信器に対する未知送信器のＤＴＯとＤＦＯの計測値と、基準送信器に対する未知送信器のＤＳＲとＤＳＲＲの計測値とに関係している。ＤＴＯとＤＦＯは前述のように計測可能であり、ＤＳＲとＤＳＲＲは未知送信器を位置確定するのに使用される幾何学的値である。ＤＳＲと位置との関係は、未知送信器位置のＤＳＲを既知送信器位置のＤＳＲに関連させているテーラー式展開を介して近似的に把握される。かくして：− る。それは、最大変化率の方向を明快に示しているベクトルである。但し、ｅ_x、ｅ_y、ｅ_zは、ｘ、ｙ、ｚの直交軸に沿っ単位ベクトルである。式（２４）と（２８）が結合され、また式（２７）と（２９）が結合されて、次の２つの式を得る：− その場合：− Δｋ₂₁ ^l(r,r₀)＝ｃ〔ＤＴＯ（未知）−ＤＴＯ（基準）〕＝ｌ₂₁(r)−ｌ₂₁(r₀) （33）さて図１３を参照すると、未知送信器位置の確定の幾何学的描写が示されており、地球３００とその中心３０４とが図示されている。ＤＳＲとＤＳＲＲに対応した式（３２）と（３１）の左手側ベクトル勾配項は、各々ベクトル３１０と３１２として示されている。基準送信器現場は、３１４でベクトルｒ₀によって示されており、未知源は３１６でベクトルｒによって示されている。第１の半角度３１８は、地球表面上の第１円３２０によって地球中心３０４で限界が定められている。未知送信器は、円３２０と３２４の交点で特定される。ｒに対する式（３１）と（３２）とを解く為に、左手側と右手側とは全く等しいものとして扱われ、第１の概算値ｒ1は、左手側のrに置き換えられる。式（３１）と（３２）の右手側は、ＤＴＯ及びＤＦＯの計測値と、人工衛星１４と１６の位置及び速度と、基準送信器２２及び受信器１８Ａと１８Ｂの位置に左右される。式（３１）と（３２）の左手側のベクトル勾配の項は、人工衛星１４と１６の位置及び速度と基準送信器２２の場所に依存している。従って、全ての項は（未知送信器の現場の第１推定値を定義する）ｒに対して代入されたベクトルｒ₁ ３１６を除いて、式（３１）と（３２）では既知である。式（３１）と（３２）の右手側は、或る基準場所３１４と人工衛星の位置及び速度に対して一定であり、それらは、（各々）ＤＳＲとＤＳＲＲの勾配ベクトル３１０と３１２を伴った未知源の位置ベクトル３１６のベクトルドット結果は一定であると言う概念を構成する。もし地球が球形であるとしたら、未知源の位置ベクトル３１６は一定の大きさとなり、一定大きさの点軌跡は、地球の中心における一定の半角度３１８の範囲を定めている円３２０となろう。同様に、未知信号源は、地球の中心における一定の半角度３２２の範囲を定めている第２円３２２を（球形の）地球の表面上に描くように拘束されている。円３２０と３２４は、２つの位置３２６と３２８で交わる。好ましい人工衛星の幾何学形状としては、これら位置の一方３２６が両方の人工衛星１４と１６に見えているが、然し他方３２８が見えなく、度外視されいるものでる。両方の人工衛星に見える前者の位置は、未知源１０の必要位置である。これら２つの位置についての解答は、一対の連立方程式を簡単に解くことによって確定され、ＤＳＲとＤＳＲＲの計測に対応した２つの位置線の明確な確定を必要としない。更に精度を高めるに、第１の確定解答r₁が、新しい「基準」位置として使用され、また式（３１）と（３２）が計算し直されて次のものを提供する：− そこで：− Δｋ₂₁ ^l(r,r₁)＝Δｋ₂₁ ^I(r,r₀)−〔ｌ₂₁(r₁)−ｌ₂₁(r₀)〕（37） Δｋ₂₁ ^nu(r,r₁)＝Δｋ₂₁ ^nu(r,r₀)−〔ν₂₁(r₁)−ν₂₁(r₀)〕（38）式（３７）のΔｋ^l ₂₁（r、ｒ₁）は、計測されたアップリンクのＤＳＲと確定された位置ｒ₁に基ずいて算定されたものとの間の差であることか理解される。同様に、式（３８）のΔｋ^nu ₂₁（ｒ、ｒ₁）は、計測されたアップリンクのＤＳＲＲと確定された位置ｒ₁に基ずいて算定されたものとの間の差であることが理解される。ｒ_j-1からｒ_jを再算定するプロセスは、差|ｒ_j−ｒ_j-l｜が或る所定量（一般に５０ｋｍ）より小さくなるまで反復される。５０ｋｍの収斂を達成するのに必要とされる反復回数は、一般に３回であるが、人工衛星の軌道の都合に左右される。初期の（又は全地球的な）位置確定解答について収斂が達成されると、次の段階は、改良された（又は局地的な）位置確定を行うことである。この関係で、全地球的な位置確定は、とりわけ球形の地球と仮定した為に不正確である。全地球的な位置確定は、地球の本当の形状等、小さな不安定要素を考慮して、結果的により精度を高めることかできる。式（２４）と（２７）は、ＤＳＲとＤＳＲＲを観察されたＤＴＯとＤＦＯに関係付けている。初期の（又は全地球的な）位置確定解答は、式（２４）と（２７）の右手側に代入されて、これらの式の左手側の値の予測を行う。観察されたものと予測されたもののＤＴＯとＤＦＯ間の差は、位置の誤差に関する誤差項を引き出すのに使用される。式（３１）と（３２）において、ｒに対して僅かに増加変化ｄｒさせることで、次のものが示される：− 但し、式（３９）と（４０）の右手側は、基準値に対する観察されたものと予測されたもののＤＳＲとＤＳＲＲ間の差である。式（３７）と（３８）の左手側のｄｒ項は、未知送信器位置の残留誤差である。未知送信器は地球表面上になければならないと言う制約の為に、ｄｒ項は、ただ東と北の方向に沿った方位面に各々要素ｄＥとｄＮを有しているだけである。かくして：− ｒ＝ｒ_j＋ｄｒ＝ｒ_j＋ｄＥｅ_E＋ｄＮｅ_N （41）式（３９）と（４０）は、次のように式（４１）を使用している成分の項に記入される一対の連立方程式を構成する：− 但し、・ｅ_Nとｅ_E及びｅ_Nは、点ｒで位置的に各々束と北を指している単位ベクトルである。式（４２）から、ｄＥとｄＮの項が推定され、未知の送信器位置ｒの推定値を正確にするのに使用される。これは、新しい位置推定値ｒを提供し、局地位置確定の解答作業が反復され、更に収斂解答が得られるまで位置推定値の精度向上が行われる。局地位置確定の解答作業は正確な形の式を使うので、得られた解答は、計測や伝搬や天体暦の誤差によって設定された範囲までは正確なものとなろう。以上の説明は、図解の目的にＤＳＲとＤＳＲＲの計測の組合わせを使用した位置確定に集中して来た。この関連で、式（２４）と（２７）はＤＳＲをＤＴＯに、ＤＳＲＲをＤＦＯに関係付けていることに注目される。観察の次の組合わせを使って未知源を位置確定することも可能である。（ａ）人工衛星位置が効果的に変化している場合に異なった時間に一対の人工衛星から得られたＤＳＲ計測値；（ｂ）人工衛星位置が効果的に変化している場合に異なった時間に一対の人工衛星から得られたＤＳＲＲ計測値：（ｃ）人工衛星位置が効果的に変化している場合に同様な時間に、又は全体的に異なった時間のいずれかで、例えば３つの人工衛星の異なった一対の組合わせを使うことで第２対の人工衛星から得られたＤＳＲ計測値と共に、第１対の人工衛星から得られたＤＳＲ計測値；（ｄ）人工衛星位置が効果的に変化している場合に同様な時間に、又は全体的に異なった時間のいずれかで、例えば３つの人工衛星の異なった一対の組み合わせを使うことで第２対の人工衛星から得られたＤＳＲＲ計測値と共に、第１対の人工衛星から得られたＤＳＲＲ計測値；（ｅ）人工衛星位置が効果的に変化している場合に、同じ対の人工衛星から、然し異なった時間に得られたＤＳＲＲ計測値と共に、一対の人工衛星から得られたＤＳＲ計測値；（ｆ）人工衛星位置が効果的に変化している場合に同様な時間に、又は全体的に異なった時間のいずれかで、異なった一対の人工衛星から得られたＤＳＲＲ計測値と共に、一対の人工衛星から得られたＤＳＲ計測値；図１から図１３を参照して述べたＤＦＯとＤＳＲの計測値の組み合わせは、無変化のままか又は入れ替わるかもしれない一対の組合わせや人工衛星位置には無関係に通常適用可能なので、最も柔軟な選択である。それは、利用可能な異なった人工衛星の一対の組合わせに頼ることがないし、また２つの独立した計測値に成されるように十分に長い期間に渡って検出可能になっている未知信号にも頼らない。上記の代わりの選択（ａ）から（ｆ）は、一度以上同じタイプの計測値を使用し、付加的な位置情報は、中継人工衛星の一対の組合わせ又は位置を変えることで得られる。選択（ａ）〜（ｆ）からのデータ分析へのアプローチは、ここで得たものから容易に得られ；後者の分析はＤＴＯとＤＦＯの確定に基ずいた式を提供する。選択（ａ）から（ｆ）に必要な全ては、（そのケースがそうであるように）ＤＴＯ又はＤＦＯに関するそれら式を２回使うことである。例えば、式（２８）と（２９）では、ｌ₂₁（ｒ）とν₂₁（ｒ）に対する表現の代わりに、各時間に又は各一対の組合わせに対してｌ₂₁（ｒ）とν₂₁（ｒ）のいずれかに対して２つの表現が存在することになろう。式（４２）は、更に計測や他の誤差の影響によって位置確定誤差を推定するのにも使用される。この場合、ｄｋ項は誤差を表し、ｄｒ項は位置確定誤差である。計測や他の誤差は、しばしばランダムに述べられるので、統計的手法が使用され、ルートは引き出された平方位置誤差を意味している。表１と表２は、本発明に従って基準信号を使用しまた使用せずに各種の誤差項によって影響されるＤＴＯとＤＦＯの計測値における有り得る誤差レベルを示している。表１ＤＴＯ誤差における基準信号の効果表２ＤＦＯ誤差における基準信号の効果表１は、正常化後のＤＴＯにおける優勢的な誤差が人工衛星位置誤差であることを示しており、表２は、正常化後のＤＦＯにおける優勢的な誤差が人工衛星速度誤差であることを示している。これらの表は、本発明に従って基準信号を使用して達成可能な誤差低減において非常に重要な改良が成されていること、即ち人工衛星位置誤差に関して一桁以上に程度が優れ且つ人工衛星速度誤差に関して二桁以上に程度が優れていることを示している。表１と表２における結果は、±０．０５度の経度帯域と±０．０５度の緯度帯域以内に人工衛星を維持する静止軌道の日常的な自己位置保持を行っている人工衛星に対して位置と速度の誤差の内の一般的な値を使用して得られてたものである。２つの人工衛星は、東経７度と東経１０度に位置決めされ、「未知の」信号はフランスのパリ（約東経２．５度と北緯５０度）に位置決めされ、基準信号送信器と監視ステーションはＵＫのデフォード（約西経２．１４度と北緯５２．１度）に共に位置決めされていた。本発明に従って得られる長所を今説明する。本発明は、離れた受信器現場３２Ａと３２Ｂにおける未知の位置の信号に対する位相と周波数の可干渉性と同時に且つ同調して既知位置の基準送信器２２からの信号を捕捉し、サンプリングし、信号処理現場３４において信号サンプルを引き続き処理する。それは、局部発振器の位相ノイズや周波数オフセット、時間に対するドリフト等の人工衛星と地上ステーションからの影響を低減し、それで基準送信器のものに対するＤＦＯの計測値の精度を改善する。この相殺は、長時間に渡る（数十秒）サンプリングが未知の送信器から隣の人工衛星チャンネルに溢れ出した非常に弱い信号の検出や、ＤＴＯ及び相対ＤＦＯの引き続いた満足な計測を実現するのに必要とされる狭い帯域信号（１００Ｈｚ未満）にとって特に重要である。改善された相対ＤＦＯの計測精度は、人工衛星や地上ステーションの受信器発振器の位相ノイズの限界によって以前は起こり得た約１００ミリＨｚに比較して数ミリＨｚの精度でＤＦＯが計測され得るようにしている。同様に、相対ＤＦＯの計測は、傾斜した静止軌道での人工衛星の衛星運動に因る周波数ドリフト効果がほぼ相殺されるようにしている。この効果は、未知信号の所定の帯域幅と引き続いた集積時間Ｔに対して、ＤＦＯ処理の周波数スロットの最小分解能Δｆは、この集積時間によって定義され限定される。この集積時間に渡って、信号のＤＦＯは、周波数スロットによってドリフトしてはならない。次いで、これは赤道面に対して静止軌道の最大可能な傾斜角度を限定している。これらを考慮して、生のＤＦＯ計測値の最大傾斜が定義される。同様に、基準送信器に対するＤＦＯ計測値の最大傾斜が定義される。この後者の傾斜は、未知送信器に対する基準送信器の位置に依存しており、表３の目的に対して、緯度５２度の未知送信器の真南の緯度２度の基準送信器が想定されている。表３生及び相対ＤＦＯ測定値に対する最大許容傾斜角度表３は、一般的な従来技術のシステムが、地球の赤道に対して０．１度より大きな傾斜角度の軌道を有した人工衛星とは効果的に働くことができないことを示している。図１から図１３を参照して説明した本発明の例は、共通したタイミングと、周波数及び位相の基準とを採用している。このことは、自から共通したタイミングと周波数及び位相の基準とに整合されていない地球上の離れた場所で未知信号が捕捉されるようにするものであり、引き続きのＤＴＯと相対ＤＦＯの計測の質低下は回避される。このことは、従来技術とは違って共通した地域をカバーする必要が無く、従って単一の地上ステーションの現場で監視できる必要がない人工衛星１４と１６のダウンリンク担当域内に受信現場が配置されるようにしている。未知信号は、基準信号のように両人工衛星１４と１６によって受信可能でなければならないが、このことがケースとなっている多くの実際の周囲条件が存在しているが故に、このことは重大な限界にはなっていない。しかし、単一の地上ステーション現場から両人工衛星を監視すると言う従来技術の特徴は、人工衛星のダウンリンクが必ずしも共通した担当域を持つ必要がないので重大な限界となっている点である。図１から図１３を参照して説明した本発明の例は、更に、人工衛星の天体暦誤差に依る効果を削減可能としている。人工衛星の天体暦誤差のＤＦＯとＤＴＯに関するインパクトは、未知と基準の送信器の間では広範囲に渡って共通しており、従って実質的な相殺が可能となっている。例えば、表２と３に示されているように、東経７度と東経１０度に配置された一対の静止人工衛星は、ＤＦＯの計測値に２Ｈｚの誤差を、ＤＴＯの計測値に２μｓの誤差を持ち込む一般的な自己位置保持誤差を有している。このことは、ＤＴＯに関して約１０ｋｍの、ＤＦＯに関しては２３４０ｋｍの位置誤差を持ち込むことになろう。既知位置に送信器を使用することで、天体暦誤差の実質的な修正を可能にする。例えば、パリ（東経２．５度、北緯５０度）に配置された「未知」の送信器に対して、未知現場の真西に１０度のところに配置された基準信号に対するＤＴＯとＤＦＯの減算は、天体暦誤差による相対ＤＴＯ誤差を０．２９μｓまで、天体暦誤差による相対ＤＦＯ誤差を３５ｍＨｚまで減じる。これは、１．５ｋｍのＤＴＯ位置誤差と４２ｋｍのＤＦＯ位置誤差に対応している。未知現場の真南１０度に配置された基準信号のＤＴＯとＤＦＯの減算は、天体暦誤差による相対ＤＴＯを０．０６５μｓま、天体暦誤差による相対ＤＦＯ誤差を５０ｍＨｚまで減じる。これらの誤差は、ＤＴＯとＤＦＯに対する０．３ｋｍと６０ｋｍの位置誤差に各々対応している。もし、基準送信器が未知信号により接近して利用できれば、天体暦誤差のインパクトは比例して減ぜられる。かくして、未知位置の真西１度の基準送信器１に対しては、天体暦誤差によるＤＴＯ誤差は、０．１５ｋｍの位置誤差に対応して０．０３μｓとなり、天体暦誤差によるＤＦＯ誤差は、４．３ｋｍの位置誤差に対応した３．６ｍＨｚとなり、それらの誤差の０．１は、経度において１０度離れた基準信号で達成される。本発明の副次的な長所は、次のようなものである：− （ａ）テーラー級数展開を使って、位置線の演算を必要とせずに未知源に急速に収斂できる。（ｂ）時間／周波数の変域を含んだディジタルな信号処理技術を使用したＣＡＦの実行は、欲しくない信号成分が周波数又は時間の変域で容易に除去され得るようにする。もし除去されなければ、これらの欲しくない成分は、信号が検出されないように及び／若しくはＤＴＯ計測精度が低減されるように大きな質低下を結果のＣＡＦに生ぜしめる。（ｃ）信号のサンプリングとディジタルな信号処理のアプローチは、位置確定を行うのに使用されるデータからの信号を明確に再構築できるようにし、これによって特定の信号が特定の位置に関連されていることを証明可能としている。これは、データが処理で壊れたり、特定の確定位置に信号を関連させるのに他に手段を必要としている従来技術に対しての改良となる。さて、図１４を参照すると、一般に４００で示された基準送信器が設けられている。基準送信器４００は、２つのミキサー４０４Ｘと４０４Ｙの各々に変調信号を与える変調された波形の発生器４０２を組み込んでいる。これらのミキサーは、基準発振器４０８からの信号に対して両方共位相固定されている各発振器４０６Ｘと４０６Ｙからの１４．０１５ＧＨｚの呼称上等しい周波数を有した局部発振器信号を受信する。ミキサー４０４Ｘと４０４Ｙからの出力信号は、各増幅器４１０Ｘと４１０Ｙに至り、次に第１と第２の送信アンテナ４１２Ｘと４１２Ｙに各々に達する。基準送信器４００は、つぎのように作動する。ミキサー４０４Ｘと４０４Ｙは基準発振器信号に対して両方共位相固定された１４．０１５ＧＨｚに集中された変調信号から成る出力信号を発生する。結果的に生じる信号は、４１０Ｘと４１０Ｙで増幅され、次いで各アンテナ４１２Ｘと４１２Ｙから送信される。第１アンテナ４１２Ｘは、第１人工衛星１４に主送信ローブを向けており、第２アンテナ４１２Ｙは、第２人工衛星１６に主送信ローブを向けている。アンテナ送信信号は、人工衛星１４と１６を使用した通信交換との衝突を回避するために一般的な通信信号のものよりも実質的により低い振幅に設定されている。しかし、その振幅は、受信器１８Ａと１８Ｂでの検出にとって十分に高く、また通信送信器のサイドローブの振幅よりもより高くなっている。従って、それらは前に述べた同等のサイドローブに比較してノイズ量に対して改良された信号を提供する。この構成の長所は、２つの部分から成っている。第１に、基準信号の周波数は、人工衛星１４又は１６のいずれかに、例えばチャンネル性能が正常な通信信号に対して受け入れ難く、従って通常は採用されないようなトランスポンダーの縁部にアクセスすることが無い周波数を占有するように選択される。第２に、基準信号のレベルは、信号チャンネル負荷のインパクトが無視できるように両方の人工衛星１４と１６において低いレベルに低減される。図１５を参照すると、一般に５００で示された基準送信器が示されている。基準送信器５００は、２つのミキサー５０４Ｘと５０４Ｙの各々に変調信号を与える疑似ランダムバイナリシーケンス（ＰＲＢＳ）発生器５０２を組み込んでいる。これらのミキサーは、基準送信器５０８からの信号に対して両方共位相固定されている各発振器５０６Ｘと５０６Ｙからの１４．０１５ＧＨｚの名目上等しい周波数を有した局部発振器信号を受信する。それらは、更にミキサー５０４Ｘと５０４からの信号も送られる電力結合器５１２Ｘと５１２Ｙに可変減衰器５１０Ｘと５１０Ｙを介して搬送波信号も与える。電力結合器５１２Ｘと５１２Ｙからの出力信号は、各増幅器５１４Ｘと５１４Ｙに至り、次に第１と第２の送信アンテナ５１６Ｘと５１６Ｙに各々達する。基準送信器５００は、次のように作動する。ミキサー５０４Ｘと５０４Ｙは、基準発振器信号に対して両方共位相固定された疑似ランダムバイナリシーケンス（ＰＲＢＳ）の変調された１４．０１５ＧＨｚ搬送波から成る出力信号を発生する。ミキサーの出力信号は、電力結合器５１２Ｘと５１２Ｙにおいて付加搬送波成分を受け取るが、それら成分は、可変減衰器５１０Ｘと５１０Ｙによって都合の良い振幅に調節されている。結果的に生じる結合された信号は、５１４Ｘと５１４Ｙで振幅され、次に各アンテナ５１６Ｘと５１６Ｙから送信される。第１アンテナ５１６Ｘは、第１人工衛星１４に主送信ローブを向けており、また第２アンテナ５１６Ｙは、第２人工衛星１６に主送信ローブを向けている。アンテナ送信信号は、人工衛星１４と１６を使用した通信交換との衝突を回避するために、一般的な通信信号のものよりも実質的により低い振幅に設定されている。しかし、それら振幅は、受信器１８Ａと１８Ｂでの検出のために十分に高く、通信送信器のサイドローブの振幅よりも高い。従って、それらは、前に説明したサイドローブの等値に比較して信号対ノイズの比率を改善している。各アンテナ送信信号における付加搬送波成分は、受信器１８Ａと１８で受信されたＣＷ信号成分の観察から直接決められるので、基準信号のＤＦＯの確定を容易にする。更に、それは、位相の質低下は、各チャンネルで受信されたＣＷ信号で直接観察されるので、位相ノイズの相殺処理を容易にする。引き続きのＣＡＦ処理の作用は、受信されたダウンリンクの基準信号が基準信号のＤＦＯシフトを包含しているので、基準信号に対して未知信号のＤＦＯを引き出すことである。信号のＰＲＢＳ成分は、基準ＤＴＯにおいて単一のピークを有してうまく限定された時間変域を提供してくれる。これは、隠された反復性を含んで、従ってＤＴＯの本当の値として混乱を起こしながら時間変域ＣＡＦに一つ以上の相関ピークを発生させる。上記例では、基準ＤＴＯとＤＦＯの計測を未知ＤＴＯとＤＦＯを確定する目的で説明した。これは精度を最大限高めることになるが、基準ＤＴＯを計測するのは実は必須ではない。ＣＡＦ処理によって直接未知ＤＴＯを計測することは或る目的で可能であり適切なことである。代わりに、基準ＤＴＯは、基準送信器２２と中継人工衛星１４と１６の位置から幾何学的に確定されよう。しかし、全ての実際の目的のために、基準ＤＦＯを計測して未知ＤＦＯを確定することは必要なことであり、これは、ＤＦＯ計測の必要な精度がＨｚ×１０^-3のオーダであるからであり、また計測された基準ＤＦＯを使うことで、中継人工衛星１４と１６によって、導入された数Ｈｚのオーダでのこの計測値の誤差を補償することになる。

───────────────────────────────────────────────────── フロントページの続き (81)指定国ＥＰ(ＡＴ，ＢＥ，ＣＨ，ＤＥ，ＤＫ，ＥＳ，ＦＲ，ＧＢ，ＧＲ，ＩＥ，ＩＴ，ＬＵ，ＭＣ，ＮＬ，ＰＴ，ＳＥ)，ＡＭ，ＡＴ，ＡＵ，ＢＢ，ＢＧ，ＢＲ，ＢＹ，ＣＡ，ＣＨ，ＣＮ，ＣＺ，ＤＥ，ＤＫ，ＥＥ，ＥＳ，ＦＩ，ＧＢ，ＧＥ，ＨＵ，ＩＳ，ＪＰ，ＫＥ，ＫＧ，ＫＰ，ＫＲ，ＫＺ，ＬＫ，ＬＲ，ＬＴ，ＬＵ，ＬＶ，ＭＤ，ＭＧ，ＭＮ，ＭＷ，ＭＸ，ＮＯ，ＮＺ，ＰＬ，ＰＴ，ＲＯ，ＲＵ，ＳＤ，ＳＥ，ＳＧ，ＳＩ，ＳＫ，ＴＪ，ＴＭ，ＴＴ，ＵＡ，ＵＧ，ＵＳ，ＵＺ，ＶＮ【要約の続き】ら、未知信号源（１０）の位置が計算される。

Claims

【特許請求の範囲】 1.（ａ）各信号中継器から未知の信号を受信する複数の受信器を配置し；（ｂ）単一信号の複製か、又は周波数、時間及び位相において単一信号に対してロックされた信号の複製のいずれかであって、公知の場所の基準送信手段から信号中継器に送信される各基準信号を各信号中継器から受信する受信器を各々配置し；（ｃ）各受信器によって受信された各未知の信号と基準信号とを、これらの信号が他で受信された信号から独立して、相互にそれらのタイミングと位相の情報を保存するようにまとめて処理し；（ｄ）処理された基準信号と処理された未知信号の交差あいまい関数演算を実施し、次の（i）未知信号の差分時間オフセット（ＤＴＯ）と差分周波数オフセット（ＤＦＯ）の値、（ii）異なった信号中継器位置に対応した未知信号のＤＴＯの値、（iii）異なった信号中継器位置に対応した未知信号のＤＦＯの値、（iv）信号中継器の異なった組合せに対応した未知信号のＤＴＯの値、（v）信号中継器の異なった組合せに対応した未知信号のＤＦＯの値の内の少なくとも一つを確定するために未知信号における位相ノイズと周波数ドリフトの作用を打ち消すように処理された基準信号を採用し；（ｅ）工程（ｄ）で確定されるＤＴＯ及び／若しくはＤＦＯの値から未知の信号源の位置を算定する上記工程を有していることを特徴とする複数の信号中継器によって受信された未知信号源を確定する方法。 2．工程（ｃ）の処理は、未知信号と基準信号とを所定の帯域幅を有した中間周波数（ＩＦ）信号に下げ変換することによって、またサンプル採取タイミングと下げ変換とを正確な周波数とタイミング基準に従って制御してそれらのディジタルサンプルを得ることによって各受信器で受信された信号について実施される請求の範囲第１項記載の方法。 3．未知信号と基準信号は、４Ｍｈｚより大きくないＩＦ帯域幅を有している請求の範囲第２項記載の方法。 4．未知信号のＤＴＯが確定され、未知信号が下げ変換に先立って少なくともその帯域幅に接近しているＩＦ帯域幅を有している請求の範囲第３項記載の方法。 5．未知信号のＤＦ０が確定され、未知信号の帯域幅が、基準信号の帯域幅に合わせられる請求の範囲第２項又は第３項記載の方法。 6．未知信号のＤＦＯが確定され、未知信号のＩＦが、１００ＫＨｚより小さい帯域幅を有している請求の範囲第２項又は第３項記載の方法。 7．未知信号のＩＦは、１０ＫＨｚにほぼ等しい帯域幅を有している請求の範囲第６項記載の方法。 8．請求の範囲第１項の工程（ｄ）における交差あいまい関数演算は、ヒルベルト（Ｈｉｂｅｒｔ）変換手法から複素数データを発生する工程を有している上記請求の範囲のいずれかの項に記載の方法。 9．請求の範囲第１項の工程（ｄ）における交差あいまい関数は、基準信号のＤＴＯを確定するために実行される上記請求の範囲のいずれかの項に記載の方法。１０．（ａ）関連した試行周波数オフセットが必要な仮の値であることを示す関数最大値が得られるまで試行周波数オフセットの範囲に対して交差あいまい関数の数値を求めることによって基準信号の仮の値を見い出し；（ｂ）各受信器に関連した第１と第２の基準信号をそれらの周波数変域等値に変換し；（ｃ）周波数シフトがＤＦＯの仮の値から構成されていて、周波数変域において第２基準信号に対する第１基準信号を周波数シフトし；（ｄ）周波数変域基準信号のいずれもの欲しくない周波数成分も削除し；（ｅ）周波数変域において、周波数成分結果を造り出すために第１基準信号の各周波数成分の複素数共役に第２基準信号の対応した周波数成分を掛け；（ｆ）周波数成分結果を時間変域に変換し、第１と第２の基準信号間の相対時間オフセットの値の各範囲に対して交差あいまい関数の各値を発生させ；（ｇ）交差あいまい関数の一組の最大強度値を選択し、それらの間の補間法によって交差あいまい関数の強度の最大値を得て、その最大値に対応するもので、第１と第２の基準信号間の相対時間オフセットとしての基準信号のＤＴＯを導き出す上記工程を有している請求の範囲第９項記載の方法。１１．請求の範囲第１項の工程（ｄ）の処理は、基準信号のＤＦＯを得るために実行され、（ａ）各受信器に関連した第１と第２の基準信号間のＤＴＯを確定し；（ｂ）基準信号のＤＴＯに等しい第１と第２の基準信号間の相対時間シフトを導入し；（ｃ）相対時間シフト後に基準信号のサンプリングを行い；（ｄ）時間成分結果を発生させるために、各第１基準信号サンプルの複素数共役に第２基準信号の対応したサンプルを掛け；（ｅ）時間成分結果を周波数変域にフーリエ変換し、第１と第２の基準信号間の相対周波数オフセットの値の各範囲に対して交差あいまい関数の各値を発生し；（ｆ）交差あいまい関数の一組の最大強度値を選択し、それらの間の補間法によって交差あいまい関数の強度の最大値を得て、その最大強度値に対応するもので、第１と第２の基準信号間の相対時間オフセットとしての基準信号ＤＦＯを導き出す上記工程を有している上記請求の範囲のいずれかの項に記載の方法。１２．請求の範囲第１項の工程（ｄ）の処理は、基準信号のＤＴＯを得るために実行され、（ａ）時間と周波数のシフトは各々基準信号のＤＴＯとＤＦＯに等しいものであって、各受信器に関連した第１と第２の未知信号間に相対時間と周波数のシフトを導入し；（ｂ）相対時間と周波数のシフト後に未知信号をサンプリングし；（ｃ）第１と第２の未知信号をそれらの周波数変域等価物に変換し；（ｄ）試行周波数オフセットをこの変換された第２未知信号に適用し；（ｅ）変換された未知信号におけるいずれもの欲しくない周波数成分も削除し；（ｆ）周波数成分結果を発生させるために、第１の未知信号の各周波数成分の複素数共役に第２の未知信号の対応した周波数成分を掛け；（ｇ）第１と第２の未知信号間の相対時間オフセットの値の各範囲に対して交差あいまい関数の各値を発生させるために周波数成分結果を時間変域に変換し；（ｈ）交差あいまい関数の値を得るために或る範囲の試行周波数オフセットに対して工程（ｄ）から（ｇ）を繰り返し；（ｉ）交差あいまい関数の一組の最大強度値を選択し、それらの間の補間法によって交差あいまい関数の強度の最大値を得て、その最大強度値に対応するもので、第１と第２の未知信号間の相対時間オフセットとしての基準信号ＤＴＯを導き出す工程とを有している上記請求の範囲のいずれかの項に記載の方法。１３．請求の範囲第１項の工程（ｄ）の処理は、未知信号のＤＦＯを得るために実行され、基準信号処理から構成され、且つ（ａ）シフトが基準信号のＤＴＯとＤＦＯに各々等しいもので、各受信器に関連した第１と第２の基準信号間に相対時間と周波数のシフトを導入し；（ｂ）相対時間と周波数シフト後に基準信号をサンプリングし；（ｃ）時間成分結果を発生させめために、各第１基準信号サンプルの複素数共役に第２信号の対応したサンプルを掛け；（ｄ）第１と第２の基準信号の間の相対周波数オフセットの値の各範囲に対して交差あいまい関数の各値を発生させるために時間成分結果を周波数変数にフーリエ変換し；（ｅ）交差あいまい関数からの欲しくない周波数変数成分を削除し；（ｆ）位相ノイズの乱れ情報を含んだ第１と第２の基準信号間に一組の瀘過された時間変域結果を発生させるために削除された交差あいまい関数を時間変域にフーリエ変換で戻す上記行程を組み込んだ段階を有している上記請求の範囲のいずれかの項に記載の方法。１４．請求の範囲第１項の工程（ｄ）の処理は、未知信号のＤＦＯを得るために実行され、未知信号処理から構成され、且つ（ａ）基準信号のＤＴＯとＤＦＯに各々等しい第１と第２の未知信号間に相対時間と周波数のシフトを導入し；（ｂ）相対時間と周波数のシフト後に未知信号をサンプリングし；（ｃ）時間成分結果を発生させめために、各第１未知信号サンプルの複素数共役に第２信号の対応したサンプルを掛け；（ｄ）各時間成分結果に第１と第２の基準信号の瀘過された時間成分結果の複素数共役を掛け；（ｅ）時間成分結果を周波数変域にフーリエ変換し、第１と第２の未知信号間の相対周波数オフセットの或る範囲の値の各々に対して交差あいまい関数の各値を発生させ；（ｆ）交差あいまい関数の一組の最大強度値を選択し、それらの間の補間法によって交差あいまい関数の最大値を得て、その最大値に相当している第１と第２の未知信号間の相対周波数オフセットとして基準信号のＤＦＯに対する未知信号のＤＦＯを導き出す上記工程を組み込んだ別の段階を有している上記請求の範囲のいずれかの項に記載の方法。１５．上記請求の範囲第１項の工程（ｃ）における、未知送信器位置の算定は、（ａ）各ケースでの差分傾斜範囲（ＤＳＲ）が、各送信器から信号中継器までの経路の長さの差であって、未知信号ＤＴＯと基準信号ＤＴＯとから基準送信器と未知送信器のＤＳＲを含んだ式を得、（ｂ）各ケースでの差分傾斜範囲割合（ＤＳＲＲ）が各ＤＳＲの変化割合であって、未知信号のＤＦＯと基準信号のＤＦＯから基準送信器と未知送信器のＤＳＲＲを含んだ式を得、（ｃ）ＤＳＲ及びＤＳＲＲの式のテーラー展開分析を行って、未知及び基準のＤＴＯとＤＦＯから未知送信器の位置を導出し、中継器の位置及び速度と基準送信器及び受信器の位置を導出することに基づいている上記請求の範囲のいずれかの項に記載の方法。１６．その請求の範囲の工程（ｃ）において導出される位置は、その第１値であり、その一つ以上の改善された値は、テーラー展開分析における第１値の戻し代入によって、また未知送信器の位置の導出の反復によって導出される請求の範囲第１５項記載の方法。１７．（ａ）未知送信器の位置、中継器の位置及び速度並びに受信器の位置の導出された値から、未知信号のＤＴＯ及びＤＦＯの算定値を発生させ；（ｂ）未知信号のＤＴＯ及びＤＦＯの算定値と、交差あいまい処理の助けで計測されたものとの間の差から導出された未知送信器の位置に対して補正を発生させる工程によって未知送信器の位置の導出を改良することを含んでいる請求の範囲第１５項又は第１６項記載の方法。１８．処理された基準信号及び処理された未知信号は、交差あいまい関数処理に先立って記憶される上記請求の範囲第のいずれかの項に記載の方法。１９．（ａ）基準信号が単一信号の複製か又は単一信号に対して周波数、時間及び位相においてロックされた信号の複製のいずれかであり、基準信号が公知位置の基準送信手段から信号中継器に伝達されるものであって、各信号中継器から未知信号と各基準信号を受信する複数の受信器と；（ｂ）各受信器によって受信された各未知信号と基準信号を、これら信号が他で受信された信号から独立して相互にそれらのタイミングと位相の情報を保存するようにまとめて処理する手段と；（ｃ）（i）処理された基準信号と処理された未知信号について交差あいまい関数処理を実行し且つ、次の：（１）未知信号の差分時間オフセット（ＤＴＯ）と差分周波数オフセット（ＤＦＯ）の値と、（２）異なった信号中継器位置に対応した未知信号のＤＴＯ値と、（３）異なった信号中継器位置に対応した未知信号のＤＦＯ値と、（４）信号中継器の異なった組合わせに対応した未知信号のＤＴＯ値と、（５）信号中継器の異なった組合わせに対応した未知信号のＤＦＯ値の内の少なくとも一つを確定するために未知信号における位相ノイズや周波数ドリフト作用を相殺する基準信号を採用し、（ii）交差あいまい関数処理で確定されるＤＴＯ及び／若しくはＤＦＯの値から未知信号源の位置を算定する手段とを有していることを特徴とする複数の信号中継器によって受信された未知信号源を確定する装置。２０．ＤＴＯ及び／若しくはＤＦＯに対して得られる値に従って受信された信号を時間及び周波数に関してオフセットし、このオフセット信号から別の値を得るように構成されている請求の範囲第１９項記載の装置。２１．ほぼ等しい強さの基準信号を信号中継器に伝達する構成の基準信号伝達手段を有している請求の範囲第１９項又は第２０項記載の装置。２２．変調されたまた変調されていない搬送波成分を組み込んだ基準信号を信号中継器に伝達する構成の基準信号伝達手段を有している請求の範囲第１９項又は第２０項記載の装置。２３．交差あいまい関数演算に先立って、処理された基準信号と処理された未知信号を記憶する手段を有している請求の範囲第１９項又は第２０項記載の装置。