JP2004513370A - ナビゲーション衛星信号品質監視用の装置 - Google Patents

ナビゲーション衛星信号品質監視用の装置 Download PDF

Info

Publication number
JP2004513370A
JP2004513370A JP2002541407A JP2002541407A JP2004513370A JP 2004513370 A JP2004513370 A JP 2004513370A JP 2002541407 A JP2002541407 A JP 2002541407A JP 2002541407 A JP2002541407 A JP 2002541407A JP 2004513370 A JP2004513370 A JP 2004513370A
Authority
JP
Japan
Prior art keywords
signal distortion
difference
deviation
correlation
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
JP2002541407A
Other languages
English (en)
Inventor
ブレナー,マッツ・エイ
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Honeywell International Inc
Original Assignee
Honeywell International Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Honeywell International Inc filed Critical Honeywell International Inc
Publication of JP2004513370A publication Critical patent/JP2004513370A/ja
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/22Multipath-related issues
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/24Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/30Acquisition or tracking or demodulation of signals transmitted by the system code related

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Navigation (AREA)

Abstract

測位システム衛星信号歪を検出するための装置が、相関曲線に沿った位置で複数の相関測定値を求める相関器を含む。相関測定値は、受信された衛星信号と基準との相関に基づいている。信号歪検出器は、相関曲線に沿った相関測定値間の差を求め、差から信号歪を検出する。

Description

【0001】
発明の技術分野
本発明は、一般に、衛星航法システム(GPS)など、衛星ベースの測位システムに関し、より詳細には、衛星ベースの測位システムにおいて衛星によって伝送される信号の品質の監視に関する。
発明の背景および従来技術
衛星ベースの測位システムは、受信器の位置を求めるために使用されており、典型的には、衛星制御設備と、複数の衛星と、受信器と、1つまたは複数の局所または地域地上局とを含む。各衛星が、受信器の位置を求める際に受信器に有用なコードおよび何らかの所定の情報を含む信号を伝送する。受信器は、少なくとも4つの衛星のコードにそれ自体同期して、これらの衛星からの信号の情報を使用して三角測量のような処置を行って、地球の中心やGPS標準時間などの基準に対する座標および時間オフセットを求める。
【0002】
受信器は、特定の位置に束縛されず、したがって可変位置を表す。実際、衛星ベースの測位システムの目的は、受信器の位置にかかわらず、受信器がその位置を求めることができるようにすることである。他方、局所または地域地上局は固定位置にあり、これを用いて、衛星によって伝送される信号を監視する。衛星によって伝送される信号は、例えば大気状態によって悪影響を受けることがあり、これにより、受信器が不適切に位置を決定する可能性がある。したがって、地上局は、受信器がより正確な位置計算を行えるように、必要な信号補正を受信器に知らせる。この構成は、ディファレンシャル測位(differential positioning)と呼ばれる。
【0003】
本発明の地上局はまた、衛星内部の欠陥を検出するために、衛星によって伝送される信号を監視する。GPSでは、これらの欠陥が測位システム信号欠陥から使用者を保護するために厳しい要件を課しているFAAによって指定される。より大きな欠陥の少なくともいくつかを表すために、一組の試験波形がFAAによって選択されている。これらの波形は、地上局機器の認証試験のために使用される。
【0004】
従来技術は、様々な相関器スペーシングで従来のコード追跡ディスクリミネータ(discriminator)を比較することによって欠陥を判定する。図1に示されるように、衛星から受信されるコードを、衛星によって伝送されるコードの時間シフトされた複製であるコード基準の組と相関させることによって、相関曲線が確立される。例えば、図1に示されるように、7つの相関測定値を計算することができる。同相測定値IPが、受信されたコードと、受信されたコードに対してゼロ時間シフトを有する基準コードとの相関量を表す(この測定値は時間に正確な測定値と呼ばれる)。同相測定値IEは、受信されたコードと、受信されたコードに対して先になるように第1の所定の時間シフトを有する基準コードとの相関量を表す。同相測定値ILは、受信されたコードと、受信されたコードに対して後になるように第2の所定の時間シフトを有する基準コードとの相関量を表す。同様に、同相測定値IEが、第3の所定の時間シフトを用いて導出され、同相測定値ILが、第4の所定の時間シフトを用いて導出され、同相測定値IEが、第5の所定の時間シフトを用いて導出され、同相測定値ILが、第6の所定の時間シフトを用いて導出される。第1の所定の時間シフトの大きさは、第2の所定の時間シフトの大きさに等しい場合があり、第3の所定の時間シフトの大きさは、第4の所定の時間シフトの大きさに等しい場合があり、第5の所定の時間シフトの大きさは、第6の所定の時間シフトの大きさに等しくあり得る。測定される相関が単に時間シフトの関数であって、受信された衛星信号の絶対パワーの関数とはならないように、全ての測定値が正規化されると仮定する。
【0005】
このとき、第1、第2、および第3のディスクリミネータが、以下の式に従って形成される。
【0006】
【数23】
Figure 2004513370
【0007】
【数24】
Figure 2004513370
【0008】
【数25】
Figure 2004513370
【0009】
その後、以下の式に従って量d1,2、d1,3、およびd2,3を形成することによって、これらのディスクリミネータを互いに比較する。
【0010】
【数26】
Figure 2004513370
【0011】
【数27】
Figure 2004513370
【0012】
【数28】
Figure 2004513370
【0013】
量d1,2、d1,3、およびd2,3は、対応するしきい値D1,2、D1,3、およびD2,3と比較され、第1のディスクリミネータd1,2がD1,2を超える場合、第2のディスクリミネータd1,3がD1,3を超える場合、または第3のディスクリミネータd2,3がD2,3を超える場合に、欠陥が存在すると考えられる。衛星航法システムの通常の動作中、各衛星から受信される信号に対してこの試験が行われる。認証中、欠陥検出が生じたことを実証するために、FAAによって選択された各試験波形を使用して試験が行われる。
【0014】
この方法での問題の少なくとも1つは、3つの量d1,2、d1,3、およびd2,3を求めるために6つの相関器を必要とすることであり、これは、提供される有用なデータの量に関して非常に大きなハードウェアである。
【0015】
また、ピークが何らかの定められた標準から所定量だけずれているかどうか判定するために、相関ピーク全体(すなわち、時間に正確な同相測定値IPの周りの相関曲線の部分)をスキャンすることによって、地上局が欠陥を判定することが知られている。しかし、この欠陥検出構成は、かなりの量の計算力を必要とし、精度に欠ける。
【0016】
従来技術での第3の方法は、測定値IE、IE、IE、IL、IL、およびILの間の以下の割合を使用する。
【0017】
【数29】
Figure 2004513370
【0018】
【数30】
Figure 2004513370
【0019】
【数31】
Figure 2004513370
【0020】
これらの割合がそれぞれ、対応する所定の値と比較される。
本発明は、従来技術の1つまたは複数の問題を克服する構成を対象とする。
発明の概要
本発明の一態様によれば、測位システム衛星信号欠陥を検出するための装置が、相関器および欠陥検出器を備える。相関器は、相関曲線に沿った点で複数の相関測定値を求め、各相関測定値が、受信された衛星信号と基準との相関に基づいている。欠陥検出器は、相関曲線に沿った相関測定値間の差を求め、差から欠陥を検出する。
【0021】
本発明の別の態様によれば、測位システム衛星によって伝送される信号に影響を及ぼす欠陥を検出する方法が、伝送された信号を第1の基準と相関させて、相関曲線に沿った第1の点で第1の相関測定値を求めるステップと、伝送された信号を第2の基準と相関させて、相関曲線に沿った第2の点で第2の相関測定値を求めるステップと、伝送された信号を第3の基準と相関させて、相関曲線に沿った第3の点で第3の相関測定値を求めるステップと、第1の相関測定値と第2の相関測定値から第1の差を求めるステップと、第2の相関測定値と第3の相関測定値から第2の差を求めるステップと、第1の差を第1のしきい値と直接比較するステップと、第2の差を第2のしきい値と直接比較するステップと、第1および第2の差と第1および第2のしきい値との比較に基づいて、衛星での欠陥を検出するステップとを含む。
【0022】
本発明のさらに別の態様によれば、測位システム衛星によって伝送される信号に影響を及ぼす欠陥を検出する方法が、伝送された信号を基準と相関させて、複数の相関測定値を、相関曲線に沿った対応する点で求めるステップと、n>2として、n対の相関測定値から単一値を求めるステップと、単一値をしきい値と比較するステップと、比較に基づいて衛星での欠陥を検出するステップとを含む。
【0023】
これらおよびその他の特徴および利点は、本発明の詳細な考察を、図面と関連付けて見ればより明らかになろう。
詳細な説明
本発明に関連する地上局10の一部を図2に示す。地上局は、相関器12−E,・・・,12−E,12−E,12−E,12−P,12−L,12−L,12−L,・・・,12−Lを有し、ここでn+mは2より大きく、nは、後の相関測定値の数であり、mは、先の相関測定値の数であり、欠陥を判定する際に使用される。相関器12−Pは、受信された信号における通常のコードを基準14−Pと相関させて、時間に正確な相関出力IPを生成し、相関器12−Lは、受信された信号におけるコードを基準14−Lと相関させて、後の相関出力ILを生成し、相関器12−Lは、受信された信号におけるコードを基準14−Lと相関させて、後の相関出力ILを生成し、相関器12−Lは、受信された信号におけるコードを基準14−Lと相関させて、後の相関出力ILを生成し、・・・、相関器12−Lは、受信された信号におけるコードを基準14−Lと相関させて、後の相関出力ILを生成する。
【0024】
さらに、相関器12−Eは、受信された信号におけるコードを基準14−Eと相関させて、先の相関出力IEを生成し、相関器12−Eは、受信された信号におけるコードを基準14−Eと相関させて、先の相関出力IEを生成し、相関器12−Eは、受信された信号におけるコードを基準14−Eと相関させて、先の相関出力IEを生成し、・・・、相関器12−Eは、受信された信号におけるコードを基準14−Eと相関させて、先の相関出力IEを生成する。
【0025】
地上局10はプロセッサ16を有し、プロセッサ16は、本明細書で以後開示する時間に正確な相関出力IPおよび後の相関出力IL、IL、ILを使用して、欠陥が存在するかどうか判定する。別法として、または追加として、プロセッサ16は、本明細書で以後開示する先の相関出力IE、IE、IE、・・・、IEを使用して、欠陥が存在するかどうか判定する。
【0026】
時間に正確な相関出力IPを生成するために、プロセッサ16は、最適な相関が得られるまで基準14−Pをシフトさせる。この基準14−Pは、受信された信号に含まれるコードの複製であってよい。次いで、プロセッサ16は、基準14−Lが基準14−Pの複製となり、かつ第1の所定の時間量だけ基準14−Pに関して時間シフトされるように、基準14−Lを制御する。したがって、相関器12−Lが後の相関出力ILを生成する。プロセッサ16はまた、基準14−Lが基準14−Pの複製となり、かつ第2の所定の時間量だけ基準14−Pに関して時間シフトされるように、基準14−Lを制御し、ここで第2の所定の時間量は第1の所定の時間量よりも大きい。したがって、相関器12−Lが後の相関出力ILを生成する。同様に、プロセッサ16は、基準14−Lが基準14−Pの複製となり、かつ第3の所定の時間量だけ基準14−Pに関して時間シフトされるように、基準14−Lを制御し、ここで第3の所定の時間量は第1および第2の所定の時間量よりも大きい。したがって、相関器12−Lは後の相関出力ILを生成する。ILまでの残りの後の相関出力も同様に生成される。第1、第2、第3、およびそれ以降の所定の時間量は全て、後の相関出力IL〜ILが全て、図3に示される相関曲線の下向きの、または後の傾斜の上にあるように選択される。
【0027】
追加として、または別法として、相関器12−E、12−E、12−E、・・・、12−Eを、先の相関出力IE、IE、IE、・・・、IEを生成するように配置することができる。また、受信された信号におけるコードを、基準14−Pの時間シフトされた直角位相形態に相関させることによって、直角位相相関出力QE、・・・、QE、QP、QL、・・・、QLを生成することができる。この後者の代替形態によれば、欠陥表示を生成するために使用される各測定値を、対応する同相測定値と直角位相測定値のRMS(根二乗平均)として形成することができる。
【0028】
集合IE,・・・,IE,IE,IE,IP,IL,IL,IL,・・・,ILを、I,・・・,I−3,I−2,I−1,I,I,I,I,・・・,Iと表すことができ、それに対応するRMS値の集合
【0029】
【数32】
Figure 2004513370
【0030】
を、R,・・・,R−3,R−2,R−1,R,R,R,R,・・・,Rと表すことができる。
先の相関出力および後の相関出力が誤差検出のために使用される場合、プロセッサ16は、先の相関出力IE〜IE、時間に正確な相関出力IP、および/または後の相関出力IL〜ILを処理し、それにより1つまたは複数の測定差di,jを導出する。これらの測定差di,jは、以下の式に従って生成される。
【0031】
【数33】
Figure 2004513370
【0032】
または
【0033】
【数34】
Figure 2004513370
【0034】
ここで、i=−m,・・・,nであり、j=−m,・・・,nであり、負の符号は、先の傾斜での測定値を示し、正の符号は、相関曲線の後の傾斜での測定値を示す。
【0035】
この時点で、これらの測定差di,jの全てまたは部分集合から、期待される差を引き、結果として得られる差の偏差を、対応するしきい値と比較して、欠陥の存在を判定することができる。例えば、これらの差の偏差全てが使用されると仮定すると、これらの差の偏差を、以下の式に従って、対応するしきい値と比較することができる。
【0036】
【数35】
Figure 2004513370
【0037】
ここで、Edi,jは、欠陥がないときに、対応する測定差di,jそれぞれに関して期待される差である。
いくつかの場合に、測定差di,jは、熱およびマルチパスノイズによって影響を及ぼされる場合があり、これは、欠陥検出装置の感度、すなわちしきい値Di,jの大きさによっては、誤った欠陥検出を導く可能性がある。したがって、これらの場合、欠陥が実際に存在しないときに欠陥が検出される可能性があり、または存在する欠陥が全く検出されない可能性がある。
【0038】
i,jの熱ノイズの内容は、基準コード14−E、・・・、14−E、14−E、14−E、14−P、14−L、14−L、14−L、・・・、14−L間の遅延hi,jの関数として求めることができる。遅延hi,jは、lおよびlを生成するために受信信号と相関される2つの基準の間の遅延である。典型的には、hi,j=0.025〜0.05チップであり、しかしこの範囲とは異なる場合もある。di,jでの熱ノイズth1は、信号対雑音比およびth1の標準偏差(1−sigma)に依存し、以下の式によって与えられる。
【0039】
【数36】
Figure 2004513370
【0040】
ここで、Bは、ノイズの両側帯域幅である。さらに、時間に正確な基準(すなわち、基準14−P)の変動による熱ノイズに対する別の寄与th2が存在する。したがって、総熱ノイズはth=th1+th2である。マルチパスノイズmpは、アンテナ利得パターンに依存し、そのオーバーバウンドする1−σ
【0041】
【数37】
Figure 2004513370
【0042】
は、衛星高度
【0043】
【数38】
Figure 2004513370
【0044】
の関数として表される。thおよびmpの統計学的な性質は通常、地上局の設置時に識別され、統計学的な情報がパラメータ化され、その後メモリに格納される。
【0045】
熱およびマルチパスノイズの悪影響を最小限に抑える一方法は、欠陥の検出に使用される、測定差di,jそれぞれに関する複数の測定値を生み出すことである。次いで、測定差di,jそれぞれに関する複数の測定値を平均化する、またはフィルタリングすることができる。熱ノイズと、いくつかのマルチパスノイズとは、特に、1つの測定値から次の測定値に相関されないので、平均化は、熱およびマルチパスノイズの影響を低減する傾向がある。
【0046】
一例として、時間に正確な相関出力IPおよび後の相関出力ILおよびILを使用して欠陥を検出すると仮定する。したがって、以下の測定差が求められる。d0,1=IP−IL;d0,2=IP−IL;およびd1,2=IL−IL。しかし、熱およびマルチパスノイズの影響を低減するために、測定差d0,1の複数の計算値が複数の相関測定値に基づいて作成され、その結果、複数の時間に正確な相関出力IPおよび複数の後の相関出力ILが生じる。次いで、測定差d0,1の全てのそのような計算値が平均化される。同様に、測定差d0,2の複数の計算値が複数の相関測定値に基づいて生み出され、その結果、複数の時間に正確な相関出力IPおよび複数の後の相関出力ILが生じる。前と同様に、測定された差d0,2の全てのそのような計算値が平均化される。同様に、測定差d1,2の複数の計算値が複数の相関測定値に基づいて生み出され、その結果、複数の後の相関出力ILおよび複数の後の相関出力ILが生じる。やはり、測定差d1,2の全てのそのような計算値が平均化される。次いで、これらの平均値を、それらの対応するしきい値D0,1、D0,2、およびD1,2と比較して、誤差の存在を求めることができる。
【0047】
熱およびマルチパスノイズの影響を低減するための別の方法は、測定差di,j、または時間に正確な相関出力IP、後の相関出力IL〜IL、および先の相関出力IE〜IEを、例えば低域フィルタを用いて適切にフィルタリングすることである。
【0048】
熱およびマルチパスノイズの影響を低減するためのさらに別の方法は、以下の手順を実施することによるものである。この手順を説明する際、以下の式に従って共分散行列Pを定義することが有用である。
【0049】
【数39】
Figure 2004513370
【0050】
ここで、下線はベクトルを示し、E[A]は、Aの統計学的な期待値であり、ベクトルは、ベクトルの平均値であり、ベクトルは、以下の式に従って求められる。
【0051】
【数40】
Figure 2004513370
【0052】
ここで、l、・・・、l−3、l−2、l−1、l、l、l、l、・・・、lなどN+1個の相関測定値を仮定すると、k=−m,・・・,n−1に関してd=l−lk−1−Edであり、またはk=−m,・・・,n−1に関してd=R−Rk−1−Edである。上三角行列Uおよび対角行列Dは、以下の式に従って求められる。
【0053】
【数41】
Figure 2004513370
【0054】
ここで、Pは、式(6)によって与えられる共分散行列である。式(6)からわかる共分散行列Pに関して、上三角行列Uおよび対角行列Dは、例えば、コレスキー分解を用いることによって求めることができる。したがって、以下の関係を、以下の式に従って定義することができる。
【0055】
【数42】
Figure 2004513370
【0056】
ここで、
【0057】
【数43】
Figure 2004513370
【0058】
は、ベクトルを生成する無相関化偏差を表すベクトルである。式(9)は、以下の式に従って書き直すことができる。
【0059】
【数44】
Figure 2004513370
【0060】
次いで、式(6)と(10)を組み合わせて、以下の式を生成する。
【0061】
【数45】
Figure 2004513370
【0062】
式(6)と(11)を比較することによって、Dが以下の式
【0063】
【数46】
Figure 2004513370
【0064】
によって与えられること、および上で定義したDが以下のフォーマット
【0065】
【数47】
Figure 2004513370
【0066】
を有する対角行列であることがわかる。このとき、分散
【0067】
【数48】
Figure 2004513370
【0068】
が、対角行列Dから求められる。上の式からわかるように、
【0069】
【数49】
Figure 2004513370
【0070】
(ここでiは1からNまで変化する)の偏差は相関されず、分散
【0071】
【数50】
Figure 2004513370
【0072】
を有する。
欠陥を判定するための最終的な
【0073】
【数51】
Figure 2004513370
【0074】
値は、以下の式
【0075】
【数52】
Figure 2004513370
【0076】
によって得られる。
【0077】
【数53】
Figure 2004513370
【0078】
の定義で必要な
【0079】
【数54】
Figure 2004513370
【0080】
への正規化は、式(14)で行われる。値
【0081】
【数55】
Figure 2004513370
【0082】
は、小さな熱およびマルチパスノイズを有する単一値であり、複数の相関測定値に関する情報を表し、欠陥の存在を判定するためにしきい値Dと比較することができる。
【0083】
本発明の特定の修正形態を上述した。他の修正形態も、当業者には想定されよう。例えば、上述したように、
【0084】
【数56】
Figure 2004513370
【0085】
分布は、全ての関連する分布がガウス分布であるという仮定に基づいている。dの分布は、この仮定から逸脱する場合があり、ここで与えた公式に対する適切な修正が必要な場合がある。
【0086】
さらに、本発明を、FAAによって指定されるものなど衛星信号欠陥の検出に関連して上述した。これらの欠陥は、本発明の使用によって検出可能な信号歪をもたらす。また、頭記の特許請求の範囲によって具体化される本発明を使用して、マルチパスと衛星コードの相互相関の影響から生じるものなど他の信号歪を検出することもできる。
【0087】
したがって、本発明の説明は、例示としてのみ解釈すべきであり、本発明を実施するための最良形態を当業者に教示するためのものである。本発明の精神を逸脱することなく詳細を大幅に変えることができ、頭記の特許請求の範囲内にある全ての修正形態の独占的な使用権が確保される。
【図面の簡単な説明】
【図1】
衛星航法システムにおいて地上局で実施される従来技術欠陥検出を説明するのに有用な相関図を示す波形である。
【図2】
本発明に関係する地上局受信器の一部の概略図である。
【図3】
本発明による衛星航法システムにおいて地上局によって実施される欠陥検出を説明するのに有用な相関図を示す波形である。

Claims (43)

  1. 測位システム衛星信号歪を検出するための装置であって、
    相関曲線に沿った点で複数の相関測定値を求める相関器であって、各相関測定値が、受信された衛星信号と基準との相関に基づいている相関器と、
    相関曲線に沿った相関測定値間の差を求め、差から信号歪を検出する信号歪検出器と
    を備える装置。
  2. 相関測定値がそれぞれ、基準と衛星信号との異なる時間シフトを表す請求項1に記載の装置。
  3. 全ての異なる時間シフトが、後の時間シフトである請求項2に記載の装置。
  4. 全ての異なる時間シフトが、先の時間シフトである請求項2に記載の装置。
  5. 異なる時間シフトが、後の時間シフトと先の時間シフトとを含む請求項2に記載の装置。
  6. 異なる時間シフトが、ゼロ時間シフトも含む請求項2に記載の装置。
  7. 信号歪検出器が、差と、前記差の期待値との偏差を生成し、偏差をしきい値と比較して、信号歪の存在を検出する請求項1に記載の装置。
  8. 信号歪検出器が、各差と、それに対応する前記差の期待値との偏差を生成し、各偏差をそれに対応するしきい値と比較して、信号歪の存在を検出する請求項1に記載の装置。
  9. 信号歪検出器が、各差と、それに対応する前記差の期待値との偏差を生成し、信号歪検出器が、前記偏差に基づいて単一偏差値を求め、単一偏差値をしきい値と比較して、信号歪の存在を検出する請求項1に記載の装置。
  10. 信号歪検出器が、各差と、それに対応する前記差の期待値との偏差を生成し、信号歪検出器が、偏差の統計学的な性質に基づいて共分散行列を求め、共分散行列を使用して偏差に対して
    Figure 2004513370
    処置を行って、信号歪を示す信号偏差を生み出す請求項1に記載の装置。
  11. 信号歪検出器が、式
    Figure 2004513370
    を形成することによって信号歪を検出し、ここで各di,jが、一対の相関測定値iとjの差であり、Edi,jが、信号歪が存在しないときの差di,jの期待値である請求項1に記載の装置。
  12. 信号歪検出器が、式
    Figure 2004513370
    に従って信号歪を検出し、ここでdi,jが、一対の相関測定値iとjの差であり、Edi,jが、信号歪が存在しないときの差di,jの期待値であり、Di,jがしきい値である請求項1に記載の装置。
  13. 信号歪検出器が、熱およびマルチパスノイズの影響を低減するために平均化を行う請求項12に記載の装置。
  14. 信号歪検出器が、熱およびマルチパスノイズの影響を低減するためにフィルタリングを行う請求項12に記載の装置。
  15. 信号歪検出器が、熱およびマルチパスノイズの影響を低減するために平均化を行う請求項1に記載の装置。
  16. 信号歪検出器が、熱およびマルチパスノイズの影響を低減するためにフィルタリングを行う請求項1に記載の装置。
  17. 測位システム衛星によって伝送される信号に影響を及ぼす信号歪を検出する方法であって、
    伝送された信号を第1の基準と相関させて、相関曲線に沿った第1の点で第1の相関測定値を求めるステップと、
    伝送された信号を第2の基準と相関させて、相関曲線に沿った第2の点で第2の相関測定値を求めるステップと、
    伝送された信号を第3の基準と相関させて、相関曲線に沿った第3の点で第3の相関測定値を求めるステップと、
    第1の相関測定値と第2の相関測定値から第1の差を求めるステップと、
    第2の相関測定値と第3の相関測定値から第2の差を求めるステップと、
    第1の差を第1のしきい値と直接比較するステップと、
    第2の差を第2のしきい値と直接比較するステップと、
    第1および第2の差と第1および第2のしきい値との比較に基づいて、衛星での信号歪を検出するステップと
    を含む方法。
  18. さらに、第1および第3の相関測定値から第3の差を求めるステップと、第3の差を第3のしきい値と直接比較するステップとを含み、信号歪の検出が、第1、第2、および第3の差と第1、第2、および第3のしきい値との比較に基づいて、衛星での信号歪を検出するステップを含む請求項17に記載の方法。
  19. 第1、第2、および第3の相関測定値が、基準と伝送される信号との異なる時間シフトを表す請求項17に記載の方法。
  20. 全ての異なる時間シフトが、後の時間シフトである請求項19に記載の方法。
  21. 全ての異なる時間シフトが、先の時間シフトである請求項19に記載の方法。
  22. 異なる時間シフトが、後の時間シフトと先の時間シフトとを含む請求項19に記載の方法。
  23. 異なる時間シフトが、ゼロ時間シフトも含む請求項19に記載の方法。
  24. 信号歪の検出が、
    第1および第2の差それぞれと、それに対応する前記差に関して期待値との偏差を生成するステップと、
    各偏差をそれに対応するしきい値と比較して、信号歪の存在を検出するステップと
    を含む請求項17、18、19、20、21、22、または23に記載の方法。
  25. 信号歪の検出が、
    各差と、それに対応する前記差の期待値との偏差を生成するステップと、
    前記偏差に基づいて単一偏差値を求めるステップと、
    単一偏差値をしきい値と比較して、信号歪の存在を検出するステップと
    を含む請求項17に記載の方法。
  26. 信号歪の検出が、
    各差と、それに対応する前記差の期待値との偏差を生成するステップと、
    偏差の統計学的な性質に基づいて共分散行列および平均値を求めるステップと、
    共分散行列および平均値を使用して偏差に対して
    Figure 2004513370
    処置を行って、信号歪を示す単一偏差値を生み出す請求項17に記載の方法。
  27. 信号歪の検出が、単一偏差値をしきい値と比較して、信号歪の存在を検出するステップを含む請求項26に記載の方法。
  28. 信号歪の検出が、式
    Figure 2004513370
    に従って信号歪を検出するステップを含み、ここでdi,jが、相関測定値iとjの差であり、Edi,jが、信号歪が存在しないときの差di,jの期待値である請求項17、18、19、20、21、22、または23に記載の方法。
  29. 信号歪の検出が、式
    Figure 2004513370
    に従って信号歪を検出するステップを含み、ここでdi,jが、相関測定値iとjの差であり、Edi,jが、信号歪が存在しないときの差di,jの期待値であり、Di,jがしきい値である請求項17、18、19、20、21、22、または23に記載の方法。
  30. 信号歪の検出が、熱およびマルチパスノイズの影響を低減するために平均化を行うステップを含む請求項17、18、19、20、21、22、または23に記載の方法。
  31. 信号歪の検出が、熱およびマルチパスノイズの影響を低減するためにフィルタリングを行うステップを含む請求項17、18、19、20、21、22、または23に記載の方法。
  32. 測位システム衛星によって伝送される信号に影響を及ぼす信号歪を検出する方法が、
    伝送される信号を基準と相関させて、複数の相関測定値を、相関曲線に沿った対応する点で求めるステップと、
    N>2として、それぞれが一対の相関測定値に基づいて形成されているN個の値から単一値を求めるステップと、
    単一値をしきい値と比較するステップと、
    比較に基づいて衛星での信号歪を検出するステップと
    を含む方法。
  33. 各相関測定値が、基準と、伝送される信号との異なる時間シフトを表す請求項32に記載の方法。
  34. 全ての時間シフトが、後の時間シフトである請求項33に記載の方法。
  35. 全ての時間シフトが、先の時間シフトである請求項33に記載の方法。
  36. 時間シフトの少なくとも1つが後の時間シフトであり、時間シフトの少なくとも1つが前の時間シフトである請求項33に記載の方法。
  37. 異なる時間シフトが、ゼロ時間シフトも含む請求項33に記載の方法。
  38. 単一値を求めるステップが、
    相関測定値の対の間でN個の差を形成するステップと、
    N個の差と、それに対応するN個の差の期待値との偏差から単一値を求めるステップと
    を含む請求項32に記載の方法。
  39. 単一値を求めるステップが、
    相関測定値の対の間でN個の差を形成するステップと、
    N個の差と、N個の差の期待値との偏差を生成するステップと、
    偏差の統計学的な性質に基づいて共分散行列および平均値を求めるステップと、
    共分散行列および平均値を使用して偏差を無相関化して、相関していない新たな偏差を生成するステップと、
    無相関化した偏差に対して
    Figure 2004513370
    処置を行って、単一値を求めるステップと
    を含む請求項32に記載の方法。
  40. 共分散行列および平均値を使用して偏差を無相関化するステップが、共分散を使用して、無相関化されて正規化された偏差を生成するステップを含み、無相関化されて正規化された偏差に対して
    Figure 2004513370
    処置を行って、単一値を求める請求項39に記載の方法。
  41. さらに、熱およびマルチパスノイズの影響を低減するために平均化を行うステップを含む請求項32に記載の方法。
  42. さらに、熱およびマルチパスノイズの影響を低減するためにフィルタリングを行うステップを含む請求項32に記載の方法。
  43. N対の相関測定値から単一値を求めるステップが、

    Figure 2004513370
    に従って共分散行列Pを定義するステップを含み、ここで下線はベクトルを示し、E[A]は、Aの統計学的な期待値であり、ベクトルは、ベクトルの平均値であり、ベクトルは、式
    Figure 2004513370
    に従って求められ、ここでNが偏差の数であり、偏差dが、式
    Figure 2004513370
    に従って相関測定値lとlの対から形成され、Edがdの期待値であり、
    さらに、

    Figure 2004513370
    に従って、上三角行列Uおよび対角行列Dを求めるステップと、

    Figure 2004513370
    に従って
    Figure 2004513370
    を定義するステップとを含み、ここで
    Figure 2004513370
    は、ベクトルを生成する無相関化偏差を表すベクトルであり、さらに、
    上述した式から式
    Figure 2004513370
    を生成するステップと、
    上述した式から式
    Figure 2004513370
    を求めるステップとを含み、ここでDが、以下のフォーマット
    Figure 2004513370
    を有し、さらに、
    上のDから分散
    Figure 2004513370
    を求めるステップと、

    Figure 2004513370
    に従って最終的な
    Figure 2004513370
    値を求めるステップと、
    Figure 2004513370
    をしきい値Dと比較して、信号歪の存在を求めるステップと
    を含む方法。
JP2002541407A 2000-10-18 2001-10-16 ナビゲーション衛星信号品質監視用の装置 Withdrawn JP2004513370A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US09/691,690 US6987820B1 (en) 2000-10-18 2000-10-18 Apparatus for navigation satellite signal quality monitoring
PCT/US2001/032163 WO2002039136A2 (en) 2000-10-18 2001-10-16 Apparatus for navigation satellite signal quality monitoring

Publications (1)

Publication Number Publication Date
JP2004513370A true JP2004513370A (ja) 2004-04-30

Family

ID=24777562

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2002541407A Withdrawn JP2004513370A (ja) 2000-10-18 2001-10-16 ナビゲーション衛星信号品質監視用の装置

Country Status (4)

Country Link
US (1) US6987820B1 (ja)
EP (1) EP1332378B1 (ja)
JP (1) JP2004513370A (ja)
WO (1) WO2002039136A2 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012093134A (ja) * 2010-10-25 2012-05-17 Seiko Epson Corp 受信信号信頼度判定方法、コード位相誤差算出方法及び受信信号信頼度判定装置
WO2013088528A1 (ja) * 2011-12-14 2013-06-20 古野電気株式会社 Gnss信号処理方法、測位方法、gnss信号処理プログラム、測位プログラム、gnss信号処理装置、測位装置、および、移動端末
WO2013088529A1 (ja) * 2011-12-14 2013-06-20 古野電気株式会社 Gnss信号処理方法、測位方法、gnss信号処理プログラム、測位プログラム、gnss信号処理装置、測位装置、および、移動端末

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1288672A1 (en) * 2001-08-08 2003-03-05 Septentrio N.V. Method and apparatus for processing signals for ranging applications
US7657230B2 (en) * 2002-08-29 2010-02-02 Qualcomm Incorporated Procedure for detecting interfering multi-path condition
WO2004031797A1 (en) 2002-09-24 2004-04-15 Honeywell International Inc. Radio frequency interference monitor
US7555262B2 (en) 2002-09-24 2009-06-30 Honeywell International Inc. Radio frequency interference monitor
US6809684B2 (en) 2002-09-24 2004-10-26 Honeywell International Inc. Signal deformation monitor
US6873910B2 (en) 2002-10-22 2005-03-29 Qualcomm Incorporated Procedure for searching for position determination signals using a plurality of search modes
EP2463685B1 (en) * 2010-12-07 2017-08-02 Airbus DS GmbH Method of synchronizing a receiver with a received ambiguous signal having a known number of at least two peaks
FR2984524B1 (fr) * 2011-12-19 2014-11-21 Sagem Defense Securite Systeme de positionnement par satellites avec attenuation post-correlation des interferences bande etroite
FR3012620A1 (fr) * 2013-10-31 2015-05-01 Thales Sa Procede de detection d'interferences dans un signal de radio-navigation par satellite par detection d'une deformation de la fonction de correlation
CN103869334B (zh) * 2013-12-27 2016-02-24 中国科学院国家授时中心 Gnss空间信号畸变的自动识别与处理方法
US10215862B2 (en) * 2014-04-07 2019-02-26 Honeywell International Inc. Systems and methods for a code carrier divergence high-pass filter monitor
HUE053985T2 (hu) 2018-08-08 2021-08-30 Thales Man & Services Deutschland Gmbh Eljárás GNSS vevõk mûködtetésére mûhold jelek deformációjának meghatározásához
IT202000011794A1 (it) * 2020-05-20 2021-11-20 St Microelectronics Srl Procedimento per rilevare spoofing in un ricevitore gnss (global navigation satellite system), corrispondenti apparecchiatura ricevente e prodotto informatico
CN112526559B (zh) * 2020-12-03 2024-05-10 北京航空航天大学 一种多工况条件下的系统关联性状态监测方法
US11874382B2 (en) * 2021-02-11 2024-01-16 Mitre Corporation GNSS spoofing detection using peak suppression monitor
CN114879229B (zh) * 2022-05-25 2023-07-18 中国民用航空飞行学院 一种适用于gbas系统的卫星导航信号畸变检测方法

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5402450A (en) 1992-01-22 1995-03-28 Trimble Navigation Signal timing synchronizer
US5414729A (en) 1992-01-24 1995-05-09 Novatel Communications Ltd. Pseudorandom noise ranging receiver which compensates for multipath distortion by making use of multiple correlator time delay spacing
KR970005598B1 (ko) * 1994-06-29 1997-04-18 삼성전자 주식회사 대역확산통신방식의 비동기 디지탈 수신장치
KR960015842B1 (ko) * 1994-07-14 1996-11-22 삼성전자 주식회사 적응 임계치를 이용한 의사잡음 코드 동기장치
US5630208A (en) 1994-07-19 1997-05-13 Trimble Navigation Limited Adaptive multipath equalization
US5600329A (en) 1995-06-30 1997-02-04 Honeywell Inc. Differential satellite positioning system ground station with integrity monitoring
FR2742612B1 (fr) * 1995-12-15 1998-02-06 Sextant Avionique Procede et circuit de reception de signaux de positionnement par satellites avec elimination des erreurs de multitrajets
US6198765B1 (en) * 1996-04-25 2001-03-06 Sirf Technologies, Inc. Spread spectrum receiver with multi-path correction
US6047017A (en) 1996-04-25 2000-04-04 Cahn; Charles R. Spread spectrum receiver with multi-path cancellation
JP2002506606A (ja) * 1998-04-30 2002-02-26 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ コード分割多重アクセス送信機及び受信機
US6345068B1 (en) * 1998-09-16 2002-02-05 Infineon Technologies Ag Hierarchical delay lock loop code tracking system with multipath correction
US6603803B1 (en) * 1999-03-12 2003-08-05 Navcom Technology, Inc. Global positioning system receiver for monitoring the satellite transmissions and for reducing the effects of multipath error on coded signals and carrier phase measurements
US6868110B2 (en) 1999-12-01 2005-03-15 The Board Of Trustees Of The Leland Stanford Junior University Multipath and tracking error reduction method for spread-spectrum receivers
US6658048B1 (en) * 2000-04-07 2003-12-02 Nokia Mobile Phones, Ltd. Global positioning system code phase detector with multipath compensation and method for reducing multipath components associated with a received signal

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012093134A (ja) * 2010-10-25 2012-05-17 Seiko Epson Corp 受信信号信頼度判定方法、コード位相誤差算出方法及び受信信号信頼度判定装置
WO2013088528A1 (ja) * 2011-12-14 2013-06-20 古野電気株式会社 Gnss信号処理方法、測位方法、gnss信号処理プログラム、測位プログラム、gnss信号処理装置、測位装置、および、移動端末
WO2013088529A1 (ja) * 2011-12-14 2013-06-20 古野電気株式会社 Gnss信号処理方法、測位方法、gnss信号処理プログラム、測位プログラム、gnss信号処理装置、測位装置、および、移動端末
JPWO2013088529A1 (ja) * 2011-12-14 2015-04-27 古野電気株式会社 Gnss信号処理方法、測位方法、gnss信号処理プログラム、測位プログラム、gnss信号処理装置、測位装置、および、移動端末
JPWO2013088528A1 (ja) * 2011-12-14 2015-04-27 古野電気株式会社 Gnss信号処理方法、測位方法、gnss信号処理プログラム、測位プログラム、gnss信号処理装置、測位装置、および、移動端末

Also Published As

Publication number Publication date
US6987820B1 (en) 2006-01-17
EP1332378A2 (en) 2003-08-06
WO2002039136A2 (en) 2002-05-16
WO2002039136A3 (en) 2003-04-17
EP1332378B1 (en) 2012-09-26

Similar Documents

Publication Publication Date Title
JP2004513370A (ja) ナビゲーション衛星信号品質監視用の装置
US7663532B2 (en) Method for measuring distance and position using spread spectrum signal, and an equipment using the method
US11409002B2 (en) Method for operating a plurality of GNSS receivers for detecting satellite signal deformation
US5828694A (en) Determination of multipath tracking error
US7702003B2 (en) Multi-path detection method for CDMA receivers
US20020101912A1 (en) Multipath and tracking error reduction method for spread-spectrum receivers
US20100104048A1 (en) Time delay measurement
JP2001356162A (ja) 統合型ワイヤレス全地球測位システムでパイロット位相オフセット時間遅延パラメータを決定する方法
CN104155662B (zh) 基于gnss相关峰值探测器的自适应互干扰抑制方法
JP2006506830A (ja) 到着誤差の無線信号時間の検出と補償のためのシステムと方法
WO1998029755A1 (en) Gps multipath detection method and system
US20100135365A1 (en) Method and apparatus for multipath mitigation
KR100687243B1 (ko) 다중 경로 오차 제거를 위한 코드추적루프 및 다중 경로오차 제거 방법
US20060133461A1 (en) Method and apparatus for processing signals for ranging applications
US20050147191A1 (en) Extended frequency error correction in a wireless communication receiver
US7876738B2 (en) Preventing an incorrect synchronization between a received code-modulated signal and a replica code
JP2005077318A (ja) 全地球測位システム
Akos et al. Signal quality monitoring: test results
US9236903B2 (en) Multi-path detection
JP2003329757A (ja) 無線位置測定方法および装置
CN108627858B (zh) 一种基于码相位调整量的多径修复方法及装置
Jakab et al. Analysis of the multipath meter performance in the presence of anomalous satellite signals
EP4361662A1 (fr) Procédé de radiodiffusion de temps précis intègre
RU2248669C2 (ru) Способ определения местоположения абонентской мобильной станции
JP2002006027A (ja) 無線網グローバルポジショニング一体化(wgp)システムの無線端末に対する往復遅延時間(rtd)パラメータを得るための方法。

Legal Events

Date Code Title Description
A300 Application deemed to be withdrawn because no request for examination was validly filed

Free format text: JAPANESE INTERMEDIATE CODE: A300

Effective date: 20050104