JP7377866B2 - 利用不可能なgnss測定値を再作成するための方法およびシステム - Google Patents

利用不可能なgnss測定値を再作成するための方法およびシステム Download PDF

Info

Publication number
JP7377866B2
JP7377866B2 JP2021528385A JP2021528385A JP7377866B2 JP 7377866 B2 JP7377866 B2 JP 7377866B2 JP 2021528385 A JP2021528385 A JP 2021528385A JP 2021528385 A JP2021528385 A JP 2021528385A JP 7377866 B2 JP7377866 B2 JP 7377866B2
Authority
JP
Japan
Prior art keywords
frequency
measurements
gnss
carrier frequency
target
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.)
Active
Application number
JP2021528385A
Other languages
English (en)
Other versions
JP2021535401A (ja
Inventor
ドゥ・ビルドゥ,ビム
スリーバーゲン,ジャン-マリー
Original Assignee
セプテントリオ・エヌ・フェー
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 セプテントリオ・エヌ・フェー filed Critical セプテントリオ・エヌ・フェー
Publication of JP2021535401A publication Critical patent/JP2021535401A/ja
Application granted granted Critical
Publication of JP7377866B2 publication Critical patent/JP7377866B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/32Multimode operation in a single same satellite system, e.g. GPS L1/L2
    • 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/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/421Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system
    • G01S19/426Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system by combining or switching between position solutions or signals derived from different modes of operation in a single system
    • 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/21Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
    • 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/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/43Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry

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)

Description

発明の分野
本発明は、全地球航法衛星システム(Global Navigation Satellite Systems:GNS
S)に関し、特に、衛星信号と基地局信号との組合わせに基づいて位置パラメータを判定するシステムおよび方法に関する。
背景
全地球航法衛星システム(GNSS)は、地球の表面上またはその付近の任意のユーザに位置、速度および時間を提供するために広く用いられている。GNSSは、特に、「測位信号」とも称される専用のナビゲーション信号を送信する複数のナビゲーション衛星(GNSS衛星)の配列を備える。米国によって開発された全地球測位システム(Global Positioning System:GPS)が最もよく知られている。他のシステムには、欧州のガリ
レオシステム、ロシアのGLONASS、および中国のBeiDouシステムがある。
大多数のGNSS衛星は、さまざまな周波数帯域で複数の測位信号を送信する。GNSS受信機では、まず、アナログフロントエンド回路においてさまざまな周波数帯域が増幅されてベースバンドにダウンコンバートされるとともに、さまざまな衛星から信号が取得されて、専用の追跡チャネルで追跡される。受信機は、各衛星の各追跡信号ごとに、いわゆる擬似距離測定値および搬送波位相測定値を計算する。これらの測定値は、受信機クロックと衛星クロックとのオフセットによって、かつ、複数の誤差原因(大気の影響、衛星の偏り、マルチパスなど)によって偏りはあるものの、衛星までの距離を表している。搬送波位相は、擬似距離よりも正確であるが、搬送波波長の整数の不定性によって影響を受ける。
高精度GNSS測位は、典型的には、2つ以上の周波数帯域からの擬似距離測定値および搬送波位相測定値に依拠するとともに、差動モードで動作することを必要とする。たとえば、リアルタイムキネマティックモード(Real-Time Kinematic mode:RTK)では、ローバーGNSS受信機は、既知の座標を用いて、基地局に対するその位置を計算する。この目的のために、基地局は、その位置ならびにその擬似距離測定値および搬送波位相測定値をローバに送信する。ローバーは、基地局データを利用して、それ自体の測定値に影響を及ぼす誤差の大部分を補償することができる。多周波RTKモードでの相対測位精度は、センチメートルレベルまたはさらにはミリメートルレベルである。
精密単独測位(Precise Point Positioning)では、衛星クロックならびに軌道の誤差
および偏りがGNSS受信機のネットワークによって計算されて、補正値がローバに送信される。ローバは、これらを用いてそれ自体の測定値を補正する。
現在の最先端技術では、特定セットの周波数帯域を追跡するローバーは、当該セットの周波数帯域に適用可能な測定値または補正値を受信しなければならない。基地局がデータを送信する周波数帯域のセットが、ローバによって追跡される周波数帯域のセットと一致しない場合、共通の周波数帯域しか補正されず、精度の低下につながることとなる。
基地局データがGNSS受信機によって追跡されるすべての周波数帯域のために利用できない場合、より優れた精度を得るGNSS受信機が必要となる。
米国特許出願公開番号US2005/203702A1は、二重周波数ナビゲーションが依拠する2つの周波数のうちの1つが利用できない短期間の間にバックアップ二重周波数ナビゲーションを実行するための方法を開示している。当該方法は、保持された周波数についての搬送波位相測定値と、電離層屈折作用のモデルとを用いて、利用不可能な周波数上でのコードおよび搬送波位相測定値を合成することを含む。電離層屈折作用のモデルは、両方の周波数についての測定値が利用可能である場合に更新されるものである。当該出願の方法は、通常動作中に両方の周波数で測定値にアクセスできるシステムにおいてのみ用いることができる。なぜなら、当該方法は、これら周波数のうちの1つが利用できない短い期間中に推定することに依拠しているからである。したがって、上記公報の方法は、受信機が一般に上記周波数のうちの1つしか追跡しない場合には採用することができない。
概要
本発明の一局面に従うと、目標搬送波周波数(f)で少なくとも1つのGNSSパラメータ推定値
Figure 0007377866000001
を生成することによって、GNSSシステムにおいて利用不可能な測定値を再作成するための方法が提供される。当該方法は、それぞれの他の搬送波周波数(f)での1つ以上の利用可能な擬似距離測定値(P)から、当該目標搬送波周波数(f)での擬似距離推定値
Figure 0007377866000002
を導出するステップと、それぞれの他の搬送波周波数(f)での当該1つ以上の利用可能な擬似距離測定値(P)および1つ以上の利用可能な搬送波位相測定値(φ)から、当該目標搬送波周波数(f)での搬送波位相推定値
Figure 0007377866000003
を導出するステップとを含む。
本発明は、特に、衛星での信号生成に関していくつかの条件が満たされた場合、1つ以上の他の周波数帯域上で得られた測定値から所与の周波数帯域に関する擬似距離および搬送波位相を導出することができるという発明者らの見識に基づいている。当該1つ以上の「他の周波数帯域」は、所与の目標周波数帯域とは異なる。本発明の方法では、目標周波数帯域内の以前の測定値にアクセスする必要はない。
本発明に従った方法の利点は、基地受信機およびローバー受信機が同じセットの信号を追跡しない場合に差動動作を有効にする点である。本発明に従った方法のさらなる利点は、何らかの理由で(一時的に)ジャミングされているかまたは利用できない周波数帯域に関する測定値を入手し続けることを可能にする点である。
本発明に従った方法の一実施形態においては、当該擬似距離推定値
Figure 0007377866000004
を導出するステップ(1030)は、当該目標周波数(f)に対する対応する他の搬送波周波数(f)の比の関数である第1の補正項を、当該1つ以上の擬似距離測定値(P)のうちの1つに加えるステップを含む。
本実施形態は、目標擬似距離(P)と利用可能な擬似距離(P)との差を、許容可能な近似度にまで、比
Figure 0007377866000005
の衛星非依存関数として表現することができるという見識に基づいている。
特定の実施形態では、第1の補正項は、対応する電離層遅延(I)に比例するとともに、比の二乗から1を差引いたものとして規定される係数である。
発明者らは、目標擬似距離(P)と利用可能な擬似距離(P)との差を、許容可能な近似度にまで、積
Figure 0007377866000006
として表現することができることを見出した。この場合、Iは適用可能な電離層遅延である。
本発明に従った方法の一実施形態では、当該搬送波位相推定値
Figure 0007377866000007
を導出するステップは、当該目標搬送波周波数(f)と対応する他の搬送波周波数(f)との差の関数である第2の補正項を、当該1つ以上の搬送波位相測定値(φ)のうちの1つに加えるステップを含む。
本実施形態は、目標搬送波位相(φ)と利用可能な搬送波位相(φ)との差を、許容可能な近似度にまで、差f-fの衛星非依存関数として表現することができるという見識に基づいている。
特定の実施形態では、当該第2の補正項は、対応する擬似距離測定値(P)と、当該目標搬送波周波数(f)と当該対応する他の搬送波周波数(f)との差と、光の速度(c)の逆数と、の積である項を含む。
特定の実施形態では、当該第2の補正項は、対応する電離層遅延(I)と、当該目標搬送波周波数(f)と当該対応する他の搬送波周波数(f)との差の二乗と、目標搬送波周波数(f)と光の速度(c)との積の逆数と、の積である項を含む。
発明者らは、目標搬送波位相(φ)と利用可能な搬送波位相(φ)との差を、許容可能な近似度にまで、
Figure 0007377866000008
として表現することができることを見出した。
本発明に従った方法の一実施形態では、導出するステップは、いずれの先験的な衛星依存項にも依拠しない。
本実施形態の一利点は、導出が、利用可能な帯域で(すなわち、周波数fで)衛星によって送信される(複数の)信号からは得られない任意の衛星依存項に依拠していないので、関与する衛星についての先験的知識が不要である点である。
本発明に従った方法の一実施形態では、当該目標搬送波周波数(f)は、他の搬送波周波数(f)から100MHz以内にある。
発明者らは、本発明の変換式に関して、たとえば、信号がAltBOC変調の一部である場合と同様に、P項およびI項の係数が特に小さいので、送信元周波数fと目標周波数fとが互いに近い場合に特に興味深いものであることを見出した。
本発明に従った方法の一実施形態では、目標搬送波周波数(f)での信号とそれぞれの他の搬送波周波数(f)での信号との間の周波数間バイアスは衛星に依存しない。
本実施形態の利点は、別々の衛星の送信元追跡信号/目標追跡信号の間で適用できることである。
本発明の一局面に従うと、GNSSを用いて位置を判定するための方法が提供される。当該方法は、第1のセットの周波数帯域でGNSS衛星からそれぞれの追跡信号を受信するステップと、第2のセットの周波数帯域で基地局からそれぞれの補正信号を受信するステップと、当該第1のセットに存在するが当該第2のセットには存在しない少なくとも1つの目標周波数帯域について、上述の導出するステップを適用することによって、異なる周波数帯域で受信された少なくとも1つの補正信号からシミュレート済み補正信号を計算するステップと、当該目標周波数帯域で受信された当該追跡信号および当該対応するシミュレート済み補正信号に基づいて位置パラメータを計算するステップとを含む。
本発明の一局面に従うと、GNSSを用いて位置を判定するための方法が提供される。当該方法は、第1のセットの周波数帯域でGNSS衛星からそれぞれの追跡信号を受信するステップと、第2のセットの周波数帯域で基地局からそれぞれの補正信号を受信するステップと、当該第2のセットに存在するが当該第1のセットには存在しない少なくとも1つの目標周波数帯域について、前述の請求項のいずれか1項に記載の導出するステップを適用することによって、異なる周波数帯域で受信された少なくとも1つの追跡信号からシミュレート済み追跡測定値を計算するステップと、当該シミュレート済み追跡測定値および当該目標周波数帯域で受信された当該対応する補正信号に基づいて、位置パラメータを計算するステップとを含む。
本発明の一局面に従うと、GNSSシステムのジャミングに対抗するための方法が提供される。当該方法は、第1のセットのジャミングされていない周波数帯域でGNSS衛星からそれぞれの追跡信号を受信するステップと、少なくとも1つのジャミングされている周波数帯域について、上述の導出するステップを適用することによって、当該セットのジャミングされていない周波数帯域における1つ以上の受信済み追跡信号からシミュレート済み追跡測定値を計算するステップとを含む。
本発明の一局面に従うと、コード手段を備えるコンピュータプログラムプロダクトが提供される。当該コード手段は、実行されると、プロセッサに上述の方法のステップを実行させるように構成されている。
本発明の一局面に従うと、GNSSを用いて位置を判定するためのシステムが提供される。当該システムは、第1のセットの周波数帯域でGNSS衛星から追跡信号を受信するように構成された第1の受信手段と、第2のセットの周波数帯域で基地局からそれぞれの補正信号を受信するように構成された第2の受信手段と、当該追跡信号および当該補正信号に対して上述の方法を実行するように構成された処理手段とを備える。
本発明に従ったコンピュータプログラムプロダクトおよびGNSS受信機についての実施形態の技術的効果および利点は、本発明に従った方法についての対応する実施形態の技術的効果および利点に必要な変更を加えたものに対応している。
図面の簡単な説明
次に、本発明の実施形態についてのこれらおよび他の技術的特徴および利点を、添付の図面を参照してより詳細に説明する。
基地局(ガリレオ(Galileo)E1およびガリレオE5a)によって追跡される周波数とローバー(ガリレオE1およびガリレオE5b)によって追跡される周波数との間の不一致を示す図である。 ローバーによって計算された測定値に対する基地から受信した測定値のアライメントを概略的に示す図である。 基地から受信した測定値に対するローバからの測定値のアライメントを概略的に示す図である。 本発明に従った方法の一実施形態のフローチャートを示す図である。 本発明に従ったシステムの一実施形態を概略的に示す図である。 真のGPS L5搬送波位相と、GPS L2搬送波位相から導出されるGPS L5搬送波位相との間のバイアスの値を示す図である。 真のガリレオE5a搬送波位相とガリレオE5b搬送波位相から導出されるガリレオE5a搬送波位相との間のバイアスの値を示す図である。 電離層寄与を無視することの影響を説明する値を示す図である。 真のE5a擬似距離と、E5b擬似距離から推定されるE5a擬似距離との差を説明する値を示す図である。
発明の詳細な説明
上述したように、大多数のGNSS衛星は、さまざまな周波数帯域で複数の測位信号を送信する。たとえば、最新世代のGPS衛星は、L1搬送波、L2搬送波およびL5搬送波について、それぞれ1575.42MHz、1227.6MHzおよび1176.45MHzで送信する。ガリレオ衛星は、E1搬送波、E5a搬送波、E5b搬送波、およびE6搬送波について、それぞれ1575.42MHz、1176.45MHz、1207.14MHz、および1278.75MHzで送信する。この場合、E5aおよびE5bは単一の広帯域AltBOC変調として送信される。
費用、サイズまたは電力消費を最適化するために、多くのGNSS受信機は、所与の衛星によって送信される周波数帯域をすべて追跡するわけではない。
たとえば、ガリレオ受信機は、ガリレオE1信号およびガリレオE5b信号のみを追跡するように設計されてもよい。このような受信機は、E1およびE5bについての搬送波位相および擬似距離のみを生成することとなるだろう。
周波数帯域iおよびjを追跡するローバは、周波数帯域iおよびjに適用可能な測定値または補正値を受信しなければならない。たとえば、RTKローバは、ガリレオE1信号およびガリレオE5b信号を追跡している場合にはE1信号およびE5b信号に関する基地局データを受信しなければならない。基地局が代わりにE1およびE5aに関するデータを送信している場合、E1測定値のみが使用可能となり、ローバーは、より低精度の単一周波数でのソリューションに戻らなければならなくなるだろう。これを図1に示す。
本発明の実施形態の一目的は、この問題を少なくとも部分的に克服することである。発明者らは、衛星における信号生成時にいくつかの条件が満たされた場合、1つ以上の他の周波数帯域に関して得られた測定値から所与の周波数帯域に関する擬似距離および搬送波位相を導出できることを見出した。
たとえば、本発明に従った受信機の一実施形態は、図2に示されるように、E5a帯域からの基地局測定値をE5b帯域からの測定値に変換することができるだろう。図示した実施形態では、本発明の核心は、基地局が受信したE5a測定値からE5b測定値をエミュレートすることで二重周波数PVTを計算することを可能にする「E5aからE5bを導出する」ブロックによって表わされる。
図3に示される別の例示的な実施形態では、ローバからのE5b測定値はE5aに変換される。
物理的には、1つ以上の他の周波数帯域に関して得られた測定値から所与の周波数帯域に関する擬似距離および搬送波位相を導出できるということは、所与の衛星の別々の周波数帯域に関して得られた測定値間に強い相関関係が存在するという事実に関係している。別々の周波数帯域での測定値間の差をもたらす主な原因は、分散媒体である電離層、ならびに受信機および衛星の周波数間バイアスであることが知られている。
電離層は、しばしば、予測可能もしくは測定可能であるか、または場合によっては無視されてもよい。受信機バイアスは、受信機ハードウェアにおける遅延および位相シフトによってもたらされる。これらはすべての衛星に共通であり、このため、測位誤差に寄与しない。したがって、異なる周波数帯域における測定値間の相関関係を利用するという観点から、衛星の周波数間バイアスは、観察された偏差に対して最も厄介な寄与となる。
当技術分野において主流となっている見知によれば、衛星の周波数間バイアス(すなわち、衛星ハードウェアにおける遅延および位相シフト)は、衛星に依存するものであり、場合によっては時間とともに変動可能であり、各周波数帯域に特有のものであると考えらられている。衛星の周波数間バイアスが存在しているということは、ローバが追跡しているものと同じ周波数帯域から測定値または補正値を受信する必要があることの主な理由となる。
しかしながら、発明者らは、一般的な見解とは逆に、衛星の周波数間バイアスが常に衛星に依存するわけではないことを見出した。たとえば、ガリレオE5a信号とガリレオE5b信号とガリレオE5-AltBOC信号との間の周波数間バイアスは、驚くべきことに、本質的に衛星に依存していない。このことは、これらの信号が単一の広帯域変調として送信されるとともに衛星において決定論的にデジタル方式で生成されるという事実に起因する。
本発明の実施形態はこの特性を利用している。したがって、信号iおよび信号kの周波数間バイアスが衛星に依存しないようにこれらの信号iおよび信号kが衛星によって送信される場合、信号iからの測定値のみが利用可能であっても、信号kからの測定値がどのようなものになるかを予測することができる。本発明は、ローバーが基地局またはネットワークから利用可能な信号と同じ信号を追跡していない場合に、ローバーからの測定値を基地局またはネットワークから受信された測定値に合わせるのに特に有用である。
図4は、本発明に従った方法の一実施形態についてのフローチャートを提供する。当該方法は、それぞれの他の搬送波周波数(f)での1つ以上の利用可能な擬似距離測定値(P)および搬送波位相測定値(φ)から、目標搬送波周波数(f)での少なくとも1つのGNSSパラメータ推定値
Figure 0007377866000009
を導出する方法を提供する。
当該方法は、特に、目標周波数(f)に対する対応する他の搬送波周波数(f)の比の関数である第1の補正項を1つ以上の擬似距離測定値(P)のうちの1つに加えることによって、目標搬送波周波数(f)での擬似距離推定値
Figure 0007377866000010
を導出するステップ1030を含み得る。
当該方法は、特に、目標搬送波周波数(f)と対応する他の搬送波周波数(f)との差の関数である第2の補正項を、1つ以上の搬送波位相測定値(φ)のうちの1つに加えることによって、目標搬送波周波数(f)での搬送波位相推定値
Figure 0007377866000011
を導出するステップ1040を含み得る。
上述の少なくとも1つのGNSSパラメータ推定値を導出するステップは、GNSSを用いて位置を判定するための方法において用いることができる。
図2に概略的に示されるとともに図4においてさらに詳述されるシナリオでは、そのような方法は、第1のセットの周波数帯域でGNSS衛星からそれぞれの追跡信号を受信するステップ1010と、第2のセットの周波数帯域で基地局からそれぞれの補正信号を受信するステップ1020と、当該第1のセットに存在するが当該第2のセットには存在しない少なくとも1つの目標周波数帯域について、上述の導出するステップ1030;1040を適用することによって、異なる周波数帯域で受信された少なくとも1つの補正信号からシミュレート済み補正信号を計算するステップと、当該目標周波数帯域で受信された当該追跡信号および当該対応するシミュレート済み補正信号に基づいて、位置パラメータを計算するステップ1050とを含み得る。
図3に概略的に示される(図4には示されない)シナリオでは、そのような方法は、第1のセットの周波数帯域でGNSS衛星からそれぞれの追跡信号を受信するステップ1010と、第2のセットの周波数帯域で基地局からそれぞれの補正信号を受信するステップ1020と、当該第2のセットに存在するが当該第1のセットには存在しない少なくとも1つの目標周波数帯域について、上述の導出するステップ1030;1040を適用することによって、異なる周波数帯域で受信された少なくとも1つの補正信号からシミュレート済み追跡測定値を計算するステップと、当該シミュレート済み追跡測定値および当該目標周波数帯域で受信された当該対応する補正信号に基づいて位置パラメータを計算するステップ1050とを含み得る。
上述の少なくとも1つのGNSSパラメータ推定値を導出するステップは、GNSSシステムにおけるジャミングに対抗するための方法において用いることができる。当該用途では、一実施形態は、1セットのジャミングされていない周波数帯域でGNSS衛星からそれぞれの追跡信号を受信するステップと、少なくとも1つのジャミングされている周波数帯域について、上述の導出するステップ1030;1040を適用することによって、当該1セットのジャミングされていない周波数帯域における1つ以上の受信済み追跡信号からシミュレート済み追跡測定値を計算するステップとを含み得る。
図5は、本発明に従ったシステムの一実施形態を概略的に示す。当該システムは特にGNSS受信機であり得る。
GNSSを用いて位置を判定するための例示のシステムは、第1のセットの周波数帯域でGNSS衛星から追跡信号を受信するように構成された第1の受信手段110と、第2のセットの周波数帯域で基地局からそれぞれの補正信号を受信するように構成された第2の受信手段120と、当該追跡信号および当該補正信号に対して上述の方法を実行するように構成された処理手段130とを備える。
代替的なGNSSシステムは、上述のように、ジャミングに対抗するための方法を実行するように構成され得る。この場合、処理手段130は、ジャミングされていない周波数帯域に関して受信された上記追跡信号に対して上述の方法を実行するように構成される。
本発明に従ったシステムの処理手段に起因する機能は、適切なソフトウェア、構成可能なハードウェア構成要素(たとえばFPGA)、専用のハードウェア構成要素(たとえばASIC)、またはそれらの任意の組合せを含むプログラム可能なハードウェア構成要素(たとえば、汎用プロセッサ、DSP等)によって実行され得る。同じ構成要素(複数の構成要素)が他の機能を提供してもよい。
本発明はまた、コード手段を備えるコンピュータプログラムプロダクトに関する。当該コード手段は、実行されると、本発明に従った方法のステップをプロセッサに実行させるように構成されている。当該コンピュータプログラムプロダクトは、コンピュータ可読媒体(たとえば、光学メモリ、磁気メモリもしくは半導体ベースのメモリ)に格納されたソフトウェアまたはダウンロード可能なソフトウェアアプリケーションを含み得るがこれらに限定されない。
ここで、或る搬送波周波数でのGNSS測定値(擬似距離、搬送波位相、ドップラー、およびCN0)が、1つ以上の他の搬送波周波数で得られた測定値からどのように推測され得るかについて、より詳細に説明する。
所与のGNSS衛星に関する2つの搬送波周波数fおよびfの擬似距離測定値および搬送波位相測定値は、以下のように表すことができる。
Figure 0007377866000012
これらの式においては、Pはメートル単位の周波数iでの擬似距離であり、φは搬送波サイクル単位の周波数iでの搬送波位相であり、ρは衛星-受信機間の距離(および本発明に関連しない追加のクロック誤差項)であり、Nは整数位相不定性であり、Iはメートル単位の周波数fでの電離層における遅延であり、δP は擬似距離バイアスであり、δφは搬送波位相バイアスであり、cは光の速度である。
電離層遅延Iは、典型的には数メートルのオーダであり、極端な電離層作用を受けた状態では最大で200mにまで達し得る。
擬似距離バイアスδPは、衛星依存バイアスと受信機依存バイアスとの和である。両方のバイアスは典型的にはメートルレベルである。衛星依存バイアスは、通常、衛星によって送信されるビゲーションメッセージに含まれる情報を用いて補償することができる。補償後、残りの衛星バイアスは数デシメートルにまで低減される。受信機バイアスは、すべての衛星に共通であり、典型的には位置計算の際に相殺される。このため、本記載ではそれ以上の懸念はない。
搬送波位相バイアスδφは、衛星依存バイアスと受信機依存バイアスとの和である。整数部分が整数不定性Nに吸収されるので、δφの分数部分だけが関連している。受信機依存バイアスは、典型的には、位置計算の際に相殺されるので、本発明には関連していないが、衛星依存位相バイアスは、補償されない場合、位置アルゴリズムに直接影響を及ぼす。
式(1)を配列し直すことにより、φおよびPを以下のように表すことができる。
Figure 0007377866000013
この場合、Nki=N-N、δφki=δφ-δφ、および、δPki=δP-δPである。
式(2)および式(3)は、周波数fでの搬送波位相および擬似距離ならびに電離層遅延から、周波数fでの搬送波位相および擬似距離を導出できることを示唆している。
ここで、fでの搬送波位相は、以下の式を用いて、周波数fの擬似距離および搬送波位相から計算され得る。
Figure 0007377866000014
ここで、
Figure 0007377866000015
は、導出された搬送波位相測定値である。電離層遅延Iを得るためのいくつかの方法を以下に説明する。
(4)と(2)とを組合わせることにより、周波数fでの真の搬送波位相測定値と周波数fでの導出された測定値との差が以下の通りであることが分かるだろう。
Figure 0007377866000016
この場合、
Figure 0007377866000017
である。
真の搬送波位相測定値と導出された搬送波位相測定値との差は、整数不定性項Nkiとバイアス項Bk,iとの和である。搬送波位相測定値が常に未知の整数不定性で定義され
るので、不定性項は無関係であるが、バイアス項Bk,iはそのとおりではない。Bk,iは、任意の整数部分が整数不定性Nkiに吸収されるので、0から1までの値である。Bk,iは、多くの場合、衛星に依存するかまたは時間に応じて変動可能である。図6は、
一例として、(4)を用いてGPS L5測定値がGPS L2測定値(BL5,L2
から導出される場合に対応するBk,i項を示す。
図6は、1日間にわたって観察される、中緯度位置からのすべてのGPS衛星についてのBL5,L2値を示す。各トレース線は衛星の通過に対応する。BL5,L2がすべての衛星に対して異なっており、衛星の通過に関して一定ではないことが分かるだろう。
この衛星依存性により、衛星依存の影響が補償されない限り、(4)を用いて搬送波位相測定値を導出することが非実用的となる。
しかしながら、いくつかのGNSS信号の場合、Bk,iは衛星に依存していないよう
に見える。これは、たとえば、周波数fの信号および周波数fの信号が衛星によって広帯域変調として送信される場合に該当する。その場合、導出された搬送波位相は、整数サイクル分および衛星非依存バイアス分だけ真の搬送波位相から異なっているに過ぎない。衛星非依存バイアスは、受信器バイアスとは区別できないので、測位アルゴリズムにおいては不適切である。
衛星によって広帯域変調で送信される信号の例として、1176.45MHzではE5aガリレオ信号および1207.14MHzではE5bガリレオ信号がある。一例として、中緯度位置から丸一日の間にわたって観察されるすべてのガリレオ衛星についてのBE5a,E5b項が図7にグラフ化されている。
図6と比較すると、図7は、E5b搬送波位相からのE5a搬送波位相の予測可能性が、GPS L2からのGPS L5の予測可能性よりもはるかに優れていることを示している。グローバルオフセット(図7では約0.63サイクル)が存在しているが、そのようなオフセットは搬送波位相測定に固有のものであり、測位アルゴリズムには関係がない。
周波数fおよび周波数fが互いに近い(たとえば、信号がAltBOC変調の一部である場合と同様に、分離<100MHzである)場合、式(4)は特に興味深いものとなる。これは、この場合、P項およびI項の係数が小さいからである。
たとえば、fがE5b周波数であり、fがE5a周波数である場合、(4)は以下のとおりである。
Figure 0007377866000018
これは、E5b測定値および電離層遅延からE5a搬送波位相を推定する際に、E5b擬似距離における任意の誤差に対する結果の感度が小さくなる(係数0.1024分だけ減衰される)とともに、電離層遅延に対する感度がさらに小さくなる(係数0.00267)だろうことを示している。これは、電離層遅延を高い精度で認識しておく必要がないことを意味する。一例として、(7)では、IE5bに関する10メートルの誤差は、ほとんどの用途では無視できる
Figure 0007377866000019
の0.0267サイクルの誤差にしか寄与しないだろう。
電離層遅延Iは、いくつかの方法で推定することができる。これは、たとえば、以下の周知の式を用いて、2つの周波数fおよびf(fは、典型的にはfとは異なる)での2つの擬似距離の差から得ることができる。
Figure 0007377866000020
(8)では、PおよびPは周波数fおよび周波数fの擬似距離であり、BGD()は、衛星によって送信される任意の補正項である。たとえば、E5bガリレオ周波数およびE1ガリレオ周波数からの測定値を有する受信機は、E5b(Pi)およびE1(Pj)に関する擬似距離を用いて電離層遅延IE5bを推定し得るとともに、さらに、(7)を用いてE5a搬送波位相を推定し得る。
(8)を用いてIを判定することができる精度は数メートルのレベルであり、これは、低い倍率を考慮すると十分過ぎるほどである。Pi擬似距離とPj擬似距離との間の任意の受信機バイアス作用は、導出された位相のグローバルオフセットに寄与し得るものであって、関連性がない。
本発明に適した精度レベルでI遅延を得るための他の方法は、電離層マップ、たとえばSBAS衛星によって送信されるマップを用いること、または、GPS衛星によって送信されるKlobucharモデルなどのモデルを適用すること、を含む。
本発明の一変形例では、周波数fでの搬送波位相は、周波数fでの搬送波位相から、かつ、周波数fおよび周波数fの擬似距離から、(8)および(4)を用いて得られる。
擬似距離は、任意には、ノイズを減らすために周知の方法で搬送波平滑化され得る。
電離層に対する感度が低いので、(4)において完全にI項を無視することさえも可能である。たとえば、図8は、E5b搬送波位相からE5a搬送波位相を導出する場合の、中緯度局に関する丸一日にわたる導出済み搬送波位相に対するI項の寄与を示す。当該寄与は、電離層作用の増大に応じて午後に1ビット増加させるが、常に数百サイクルのレベルのままで保たれているので、差し支えなく無視することができる。
本発明の別の変形例では、導出された搬送波位相測定値は、単一の他の周波数のコードおよび搬送波位相の測定値からのみ、好ましくは近接周波数分離で、計算される。この場合、Iは(4)においてゼロに設定される。
なお、上記導出が受信側の位相ワインドアップ(phase wind-up)の影響に明確に対処
していないことに留意されたい。位相ワインドアップがすべての衛星に同じように影響を及ぼす共通のバイアスであるので、測位アルゴリズムには不適切である。
(3)から開始すると、周波数fの擬似距離からの、周波数fの擬似距離についての好適な推定量は以下のとおりとなる。
Figure 0007377866000021
ここで、fでの真の擬似距離測定値と(9)を用いて導出された測定値との差は以下のとおりである。
Figure 0007377866000022
バイアス項δPk,iは、衛星依存成分および衛星非依存成分を有する。ここでは衛星
依存成分だけが関連している。これは、衛星が送信するグループ遅延情報を用いて正確に補償することができる。代替的には、ガリレオE5a、ガリレオE5bおよびガリレオE5-AltBOCなどの広帯域変調として送信される信号の場合、δPk,iについては
衛星からの非依存性が大きいので、このδPk,iは無視しても差し支えない。δPk,iが衛星に依存していないという事実を図9に示す。グレーの色調ごとに異なる衛星に対応している。顕著な衛星依存オフセットは見られない。
式(9)は、周波数fおよび周波数fが互いに近い(たとえば、ガリレオAltBOC信号の場合と同様に、分離<100MHzである)場合に特に興味深いものとなる。というのも、この場合、I項の係数が小さいからである。たとえば、fがE5b周波数に対応しており、fがE5a周波数である場合、(9)は以下のとおりである。
Figure 0007377866000023
搬送波位相の場合と同様に、式(11)は、電離層依存項が結果に対して限られた作用しか及ぼさないことを示している。
本発明に従った受信機では、周波数fの擬似距離は(9)を用いて計算されてもよく、I値は(8)によって計算されてもよい。
いくつかのGNSS受信機はまた、さまざまな周波数でドップラーを計算する。周波数fでのドップラーから周波数fでのドップラーを推量することは当技術分野では周知である。これは以下の式を用いて行うことができる。
Figure 0007377866000024
C/N0値に関しては、通常、CN0が周波数に依存していないと想定しても差し支えないだろう。このため、以下のとおりとなる。
C/N0=C/N0
特定の実施形態を参照して本発明を上述してきたが、これは、本発明を限定するためではなく明確にするためになされたものであり、その範囲は添付の特許請求の範囲によって特定される。

Claims (14)

  1. 目標搬送波周波数(f)で少なくとも1つのGNSSパラメータ推定値
    Figure 0007377866000025

    を生成することによって、GNSSシステムにおいて利用不可能な測定値を再作成するための方法であって、
    1つ以上の利用可能な擬似距離測定値(P)から、前記目標搬送波周波数(f)での擬似距離推定値
    Figure 0007377866000026

    を導出するステップ(1030)と、
    前記1つ以上の利用可能な擬似距離測定値(P)および1つ以上の利用可能な搬送波位相測定値(φ)から、前記目標搬送波周波数(f)での搬送波位相推定値
    Figure 0007377866000027

    を導出するステップ(1040)とを含み、
    記擬似距離推定値
    Figure 0007377866000028

    および前記搬送波位相推定値
    Figure 0007377866000029

    は前記1つ以上の利用可能な擬似距離測定値(P)および前記1つ以上の利用可能な搬送波位相測定値(φ)から導出されるものであり、前記1つ以上の利用可能な擬似距離測定値(P)および前記1つ以上の利用可能な搬送波位相測定値(φ)は、前記目標搬送波周波数(f)とは異なるそれぞれの周波数(f)での測定値である、方法。
  2. 記擬似距離推定値
    Figure 0007377866000030

    を導出するステップ(1030)は、前記目標搬送波周波数(f)に対する対応する他の搬送波周波数(f)の比の関数である第1の補正項を、前記1つ以上の擬似距離測定値(P)のうちの1つに加えるステップを含む、請求項1に記載の方法。
  3. 前記第1の補正項は、対応する電離層遅延(I)に比例するとともに、前記比の二乗から1を差引いたものとして規定される係数である、請求項2に記載の方法。
  4. 前記搬送波位相推定値
    Figure 0007377866000031

    を導出するステップ(1040)は、前記目標搬送波周波数(f)と、対応する他の搬送波周波数(f)との差の関数である第2の補正項を、前記1つ以上の搬送波位相測定値(φ)のうちの1つに加えるステップを含む、請求項1~請求項3のいずれか1項に記載の方法。
  5. 前記第2の補正項は、
    対応する擬似距離測定値(P)と、
    前記目標搬送波周波数(f)と前記対応する他の搬送波周波数(f)との差と、
    光の速度(c)の逆数と、
    の積である項を含む、請求項4に記載の方法。
  6. 前記第2の補正項は、
    対応する電離層遅延(I)と、
    前記目標搬送波周波数(f)と前記対応する他の搬送波周波数(f)との差の二乗と、
    目標搬送波周波数(f)と光の速度(c)との積の逆数と、
    の積である項を含む、請求項4または5に記載の方法。
  7. 前記導出するステップ(1030;1040)は、いずれの先験的な衛星依存項にも依拠しない、請求項1~請求項6のいずれか1項に記載の方法。
  8. 前記目標搬送波周波数(f)は、他の搬送波周波数(f)から100MHz以内にある、請求項1~請求項7のいずれか1項に記載の方法。
  9. 前記目標搬送波周波数(f)での信号とそれぞれの他の搬送波周波数(f)での信号との間の周波数間バイアスが衛星に依存しない、請求項1~請求項8のいずれか1項に記載の方法。
  10. GNSSを用いて位置を判定するための方法であって、
    第1のセットの周波数帯域でGNSS衛星からそれぞれの追跡信号を受信するステップ(1010)と、
    第2のセットの周波数帯域で基地局からそれぞれの補正信号を受信するステップ(1020)と、
    前記第1のセットに存在するが前記第2のセットには存在しない少なくとも1つの目標周波数帯域について、請求項1~請求項9のいずれか1項に記載の導出するステップ(1030;1040)を適用することによって、異なる周波数帯域で受信された少なくとも1つの補正信号からシミュレート済み補正信号を計算するステップと、
    前記目標周波数帯域で受信された前記追跡信号および対応する前記シミュレート済み補正信号に基づいて、位置パラメータを計算するステップ(1050)とを含む、方法。
  11. GNSSを用いて位置を判定するための方法であって、
    第1のセットの周波数帯域でGNSS衛星からそれぞれの追跡信号を受信するステップ(1010)と、
    第2のセットの周波数帯域で基地局からそれぞれの補正信号を受信するステップ(1020)と、
    前記第2のセットに存在するが前記第1のセットには存在しない少なくとも1つの目標周波数帯域について、請求項1~請求項10のいずれか1項に記載の導出するステップ(1030;1040)を適用することによって、異なる周波数帯域で受信された少なくとも1つの追跡信号からシミュレート済み追跡測定値を計算するステップと、
    前記シミュレート済み追跡測定値および前記目標周波数帯域で受信された対応する前記補正信号に基づいて、位置パラメータを計算するステップ(1050)とを含む、方法。
  12. GNSSシステムのジャミングに対抗するための方法であって、
    第1のセットのジャミングされていない周波数帯域でGNSS衛星からそれぞれの追跡信号を受信するステップと、
    少なくとも1つのジャミングされている周波数帯域について、請求項1から9のいずれか1項に記載の導出するステップ(1030;1040)を適用することによって、前記セットのジャミングされていない周波数帯域における1つ以上の受信済み追跡信号からシミュレート済み追跡測定値を計算するステップとを含む、方法。
  13. コード手段を備えるコンピュータプログラムプロダクトであって、前記コード手段は、実行されると、請求項1~請求項12のいずれか1項に記載の方法のステップをプロセッサに実行させるように構成されている、コンピュータプログラムプロダクト。
  14. GNSSを用いて位置を判定するためのシステムであって、
    第1のセットの周波数帯域でGNSS衛星から追跡信号を受信するように構成された第1の受信手段(110)と、
    第2のセットの周波数帯域で基地局からそれぞれの補正信号を受信するように構成された第2の受信手段(120)と、
    前記追跡信号および前記補正信号に対して請求項10に記載の方法を実行するように構成された処理手段(130)とを備える、システム。
JP2021528385A 2018-11-21 2019-11-21 利用不可能なgnss測定値を再作成するための方法およびシステム Active JP7377866B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP18207645.5 2018-11-21
EP18207645.5A EP3657218A1 (en) 2018-11-21 2018-11-21 Method and system for recreating unavailable gnss measurements
PCT/EP2019/082099 WO2020104594A1 (en) 2018-11-21 2019-11-21 Method and system for recreating unavailable gnss measurements

Publications (2)

Publication Number Publication Date
JP2021535401A JP2021535401A (ja) 2021-12-16
JP7377866B2 true JP7377866B2 (ja) 2023-11-10

Family

ID=64426813

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021528385A Active JP7377866B2 (ja) 2018-11-21 2019-11-21 利用不可能なgnss測定値を再作成するための方法およびシステム

Country Status (4)

Country Link
US (1) US11906635B2 (ja)
EP (2) EP3657218A1 (ja)
JP (1) JP7377866B2 (ja)
WO (1) WO2020104594A1 (ja)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2020114350A (ru) * 2020-04-22 2021-10-22 Тримбл Инк. Приемник gnss, приспособленный для создания, использования и передачи программно-сформированных данных спутникового сигнала
EP3978955A1 (en) 2020-10-04 2022-04-06 Septentrio N.V. Method and system for determining a phase of a gnss carrier
WO2023009463A1 (en) 2021-07-24 2023-02-02 Swift Navigation, Inc. System and method for computing positioning protection levels
US12013468B2 (en) 2022-09-01 2024-06-18 Swift Navigation, Inc. System and method for determining GNSS corrections
US12019163B2 (en) 2022-09-12 2024-06-25 Swift Navigation, Inc. System and method for GNSS correction transmission
US20240094403A1 (en) * 2022-09-12 2024-03-21 Qualcomm Incorporated Cross-frequency correction for precise positioning

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007529010A (ja) 2004-03-12 2007-10-18 ナヴコム テクノロジー インコーポレイテッド 2周波数の一方で測定データが利用できない場合に短期間バックアップ2周波数ナビゲーションを行なう方法
JP2010504523A (ja) 2006-09-22 2010-02-12 ナヴコム テクノロジー インコーポレイテッド 3つのGPS周波数を用いて整数値サイクル(whole−cycle)搬送波位相アンビギュイティを解消する方法
US20100141510A1 (en) 2008-12-09 2010-06-10 Dai Liwen L Methods and systems to increase accuracy in the navigation of single frequency receivers
US20110279314A1 (en) 2009-02-22 2011-11-17 Trimble Navigation Limited GNSS Signal Processing Methods and Apparatus with Ionospheric Filters

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103163533B (zh) * 2013-03-27 2014-10-22 武汉大学 一种gnss全球与区域电离层延迟无缝融合表达和改正方法
US10524226B2 (en) * 2017-10-06 2019-12-31 Skycatch, Inc. Determining the location of a UAV in flight utilizing real time kinematic satellite navigation and precise point positioning

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007529010A (ja) 2004-03-12 2007-10-18 ナヴコム テクノロジー インコーポレイテッド 2周波数の一方で測定データが利用できない場合に短期間バックアップ2周波数ナビゲーションを行なう方法
JP2010504523A (ja) 2006-09-22 2010-02-12 ナヴコム テクノロジー インコーポレイテッド 3つのGPS周波数を用いて整数値サイクル(whole−cycle)搬送波位相アンビギュイティを解消する方法
US20100141510A1 (en) 2008-12-09 2010-06-10 Dai Liwen L Methods and systems to increase accuracy in the navigation of single frequency receivers
US20110279314A1 (en) 2009-02-22 2011-11-17 Trimble Navigation Limited GNSS Signal Processing Methods and Apparatus with Ionospheric Filters

Also Published As

Publication number Publication date
EP3884303A1 (en) 2021-09-29
JP2021535401A (ja) 2021-12-16
US11906635B2 (en) 2024-02-20
US20220011443A1 (en) 2022-01-13
WO2020104594A1 (en) 2020-05-28
EP3657218A1 (en) 2020-05-27

Similar Documents

Publication Publication Date Title
JP7377866B2 (ja) 利用不可能なgnss測定値を再作成するための方法およびシステム
RU2354991C2 (ru) Способ использования трех частот gps для разрешения целочисленных неоднозначностей фазы несущей
JP7153427B2 (ja) 衛星測位システムを用いた測位方法および測位装置
EP1762824B1 (en) Position determination using carrier phase measurements of satellite signals
TWI424184B (zh) 信號處理方法和系統
US7576690B2 (en) Position determination with reference data outage
KR101000984B1 (ko) 위성 위치 확인 시스템에서 거리 변화율 측정값을 이용하여위치를 검증하는 방법 및 장치
US10816670B2 (en) Navigation satellite system positioning with enhanced satellite-specific correction information
US6999027B1 (en) Accommodation of anti-jamming delays in GNSS receivers
KR20130014447A (ko) 의사거리 수정을 이용하는 위치 신호의 수신 방법 및 장치
US11125890B2 (en) Advanced navigation satellite system positioning method and system using seeding information
US11846715B2 (en) DGNSS using reference station carrier phase measurements
CN114935770B (zh) 一种多历元加快精密单点定位收敛速度的方法及装置
CN105044741A (zh) 一种伪距相位综合广域差分改正值的求解方法
EP3223038A1 (fr) Procédé de geopositionnement par satellites et terminal associé
Melgard et al. Interchangeable integration of GPS and Galileo by using a common system clock in PPP
JPWO2006121023A1 (ja) 測位装置および測位システム
WO2020145839A1 (en) Method and apparatus for improving the quality of position determination
Duan et al. Performance of Galileo satellite products determined from multi-frequency measurements
CN110568464A (zh) 基于bds/gnss多模芯片的精密定位方法及装置
JP2005069866A (ja) 相対測位装置
WO2022201606A1 (ja) 測位装置、測位プログラムおよび測位方法
EP2929369A1 (fr) Procede et dispositif de determination du temps de propagation d'un filtre a ondes acoustiques de surface
JP2007127579A (ja) キャリア位相相対測位装置
KR20110127274A (ko) 보조 데이터를 이용한 위치결정 방법

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20221017

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230725

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230912

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20231003

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20231030

R150 Certificate of patent or registration of utility model

Ref document number: 7377866

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150