JP2017514146A - 大気監視のための無線航法信号を処理するための方法および装置 - Google Patents

大気監視のための無線航法信号を処理するための方法および装置 Download PDF

Info

Publication number
JP2017514146A
JP2017514146A JP2017501498A JP2017501498A JP2017514146A JP 2017514146 A JP2017514146 A JP 2017514146A JP 2017501498 A JP2017501498 A JP 2017501498A JP 2017501498 A JP2017501498 A JP 2017501498A JP 2017514146 A JP2017514146 A JP 2017514146A
Authority
JP
Japan
Prior art keywords
measurement system
module
time
phase
signal
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.)
Granted
Application number
JP2017501498A
Other languages
English (en)
Other versions
JP6513179B2 (ja
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.)
European Union represented by European Commission
Original Assignee
European Union represented by European Commission
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 European Union represented by European Commission filed Critical European Union represented by European Commission
Publication of JP2017514146A publication Critical patent/JP2017514146A/ja
Application granted granted Critical
Publication of JP6513179B2 publication Critical patent/JP6513179B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/23—Testing, monitoring, correcting or calibrating of receiver elements
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining 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/40—Correcting position, velocity or attitude
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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
    • G01S1/00—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
    • G01S19/072—Ionosphere corrections
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/25—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/243—Demodulation of navigation message

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Atmospheric Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Electric Clocks (AREA)

Abstract

本発明は、GNSS（Global Navigation Satellite System：全地球航法衛星システム）の無線周波数信号の処理に基づいて、大気の監視及び測定を行うシステムを提供する。本発明の特徴は、受信した衛星信号から、振幅及び位相情報を抽出するためにオープンループ復調構成、および大気によって誘発される振幅および位相の変動に関連する統計情報を提供可能な信号処理技術である。

Description

（0001）
本発明は、地上ベースのGNSS（Global Navigation Satellite System：全地球航法衛星システム）受信機を用いる大気測定に関し、特に、電離層や電離層活動に関連する統計値の測定に関する。
（0002）
宇宙ベースの無線信号は、大気監視のために広く使用されている。これらの信号は、宇宙ベースの送信機から地球に伝播するため、大気により位相偏移、群遅延、および振幅の変化を受ける。
適切な方法でこれらの信号を処理する受信機は、これらの位相の推定量、遅延及び振幅の変化を抽出することができ、さらに大気に関するいくつかの情報を推測することができる。
GNSS信号は、その存在量の多さおよび地球規模のカバレッジ、ならびに複数の周波数で送信されるという事実により、これらの目的のために広く使用される。
これらの信号は電離層や対流圏の監視に使用される。その理由は、電離層や対流圏がＬバンドで送信された信号の伝播速度および伝播方向の変化を誘起するためである。
（0003）
これらの影響を測定するために、GNSS受信機が使用される。これらの受信機は、搬送波の周波数および位相と、これらの信号の変調された測距コードを追尾し、信号電力、搬送波位相および測距コード遅延の測定値を生成する。以後、未処理信号測定値(raw signal measurements)と呼ばれるこれらの値は、次に大気を通る信号伝播に関連する種々の特性（以下、大気測定値）の計算のために使用される。
従来、これらの未処理信号測定値を生成するための受信機は、対象のパラメータのクローズドループ追尾を行っている。典型的なシステムでは、測距コードのために遅延ロックループ（DLL）を使用し、搬送波のために位相ロックループ（PLL）を使用する。
他の多くのシステムも利用可能であるが、一般に、受信機は未処理信号の測定値を生成するために、再帰的なフィードバック/フィードフォワードメカニズムの並べ替えに依存する。
（0004）
大気測定値の計算は、利用可能性と未処理信号測定値の品質の両方に依存する。このため、受信機の追尾アルゴリズムが信号パラメータを正確に追尾することが困難になると、得られた大気測定値の品質が低下する。
受信機に実装された特定の追尾アルゴリズムが、結果として得られる大気測定値に影響を与える。例えば、追尾アルゴリズム内のフィルタリング作用または一時的なエラーは、大気測定値にアーチファクト(人工効果；artefacts)を生成することがある。
（0005）
大気中の異常（例えば、電離層シンチレーション）は、受信機の追尾アルゴリズムの正確な動作を困難にすることがある。したがって、受信機がこの異常を測定するために使用されるとき、追尾アルゴリズムが追尾に失敗し、未処理信号測定値の品質の低下またはそれらが利用できないとの理由で大気測定値の品質が大幅に低下する可能性がある。
受信機の追尾耐久力および測定値の利用可能性を高めるため、積分時間の拡張、追尾帯域の低減などの多くの技術が提案されている。しかし、これらは、未処理信号測定値の減少をもたらし、最終的には、大気測定値のアーチファクトが生成される結果となる。
（0006）
ある大気測定値の生成、例えばシグマ−ファイ（φ）として知られている電離層測定には、未処理信号測定値のフィルタリングが必要である。しばしばデトレンディング(de-trending)として知られるこのフィルタリングステージは、かなり長い収束時間を有する。未処理信号測定値が断続的に使用不能となった場合、その結果使用不能となる大気測定値の時間は、はるかに長くなり得る。
（0007）
現在のアプローチの弱点は推定段階にある。従来、受信機においては、大気測定値の計算に先立って、未処理信号パラメータの推定または追尾を行っている。理想的でない状態になると、この追尾段階に問題が発生する。従来のシステムの欠点について、理論的な要因を挙げながら、図１および図２（従来技術）に関連して説明する。
（0008）
典型的には地上受信機のアンテナ受信されるGNSS信号は次の式(１)のようにモデル化される。
Figure 2017514146
 ここで、Ssigは視野にある衛星からの信号のセット、sｉ(t)は可視の複数の衛星からのｉ番目の受信信号、n(t)は付加された熱雑音を意味する。式(１)における様々なパラメータは以下の信号特性を表す。即ち、Piは総受信電力（ワット）、ωiはラジアン／秒単位の公称のRF搬送波周波数、di(t)はバイポーダル(bi-podal)データ信号または二次コード、ci(t)は信号の拡散シーケンスおよび副搬送波、θi(t)は伝搬遅延、衛星対ユーザーダイナミックス(satellite-to-user dynamics)、大気の影響、および衛星時計の影響を含む総受信位相プロセス、プロセスτi(t)は伝搬遅延、衛星時計の影響、および大気の影響を含む受信機で観測された総遅延を表す。
（0009）
特に、搬送波の位相項、即ち式(１)におけるθi(t)は、異なる位相プロセスの数を表す。これは数学的には下記の式(２)に示す線形結合で表すことができる。
Figure 2017514146
 ここで、θ０はある任意の初期位相を表し、θLOS(t)は衛星と受信機との間の見通し線の長さ(line-of-sight geometry)/ダイナミックスによる誘起位相プロセスを表し、θSV Clk.(t)は衛星時計の誤差による位相プロセスを表し、θAtm.(t)は信号が伝播する大気による誘起位相プロセスを表す。
（0010）
図１（従来技術）は、搬送波位相（θi^）及び測距コード遅延（τi^）のローカルな推定値が相関値Yi[n]を生成するために使用される方法を示す、従来の受信機のデジタル整合フィルタ102-1のブロック図である。
（0011）
GNSS受信機は、一般に、受信RF信号のゼロまたはノンゼロの中間周波数へのダウンコンバージョンを実行し、次いで信号をサンプリングする。これらの信号サンプル(r)は、次に下記の式(３)に示される演算が実装されたデジタル整合フィルタ（DMF）102-1により処理される。
Figure 2017514146
 ここで、変数τi^およびθi^は、式(１)で定義された受信機における変数τiおよびθiの推定値であり、用語Yi[n]は相関値として知られている。
（0012）
式(３)により記述される演算は、受信機内の追尾アルゴリズム内で搬送波位相および測距コードの追尾ループの一部として実施される。
（0013）
図２（従来技術）は、搬送波追尾ループ104および測距コード追尾ループの双方のためのループを描いた典型的なクローズドループ追尾構成のブロック図である。当業者に理解されるように、DMFバンク102は複数のインスタンス102-1，102-2，102-3（即ちチャンネルあたり１つ；ここには３つのみ図示）を実装している。
（0014）
典型的な実装は図２に示されるブロック図に従う。ここで、相関値Yi[n]は、信号パラメータτi^（ｍＴＳ）およびθi^（ｍＴＳ）の推定値を生成するために、搬送波追尾ブロック104および測距コード追尾ブロック106により処理される。次いで、これらの推定値は後一連の相関値Yi[n]）を生成するために使用される。
（0015）
GNSS信号が伝播している大気の特性や属性を推定するために使用される特定の複数のアルゴリズムは、‘大気監視アルゴリズム’108との名称のブロックに実装されている。これらのアルゴリズムは、次の量を処理する。即ち、DMF102によって生成されるYi[n]、および追尾アルゴリズムにより推定されるτi^とθi^である。監視アルゴリズムの性能は未処理信号測定値の品質に直接的に影響されるので、搬送波および測距コードの双方の追尾アルゴリズムの適正な動作は大気監視受信機にとって非常に重要である。
問題は、伝搬路における高い大気活性状態の変動により、これらの追尾アルゴリズムが不十分に実行されたり、失敗したりするおそれがあることである。
（0015a）
YORK J ET AL：“Development of a Prototype Texas Ionospheric Ground Receiver（TIGR）”,ITM2012-PROCEEDINGS OF THE 2012 INTERNATIONAL TECHNICAL MEETING OF THE INSTITUTE OF NAVIGATION,THE INSTITUTE OF NAVIGATION,8551 RIXLEW LANE SUITE 360 MANASSAS,VA 20109,USA,1 February 2012（2012-02-01）, 1526-1556頁,XP056000936は、衛星信号から電離層測定を行うように設計されたソフトウェア受信機を開示している。
RFデータは、直接毎秒２ギガでADCによりサンプリングされ、FPGAに渡され、ここでデジタル的にフィルタリングされ、各々が２０MHzの帯域を有する３つの可変帯域にダウンサンプリングされる。
データ量の低減したデジタルデータストリームは第２のFPGAに渡され、ここで個々のチャンネルは、複数の狭い信号帯域にフィルタリングされる。この信号帯域の中心は、予測されたドップラーシフトを補償するために調整された衛星信号の周波数である。このデータ中の信号の位相および強度の予測は、汎用CPU上で動作するオンボードソフトウェアの使用により達成される。
（0015b）
LULICH TD ET AL：“Open Loop Tracking of Radio Occultation Signals from an Airbome Platform”,GNSS 2010-PROCEEDINGS OF THE 23RD INTERNATIONAL TECHNICAL MEETING OF THE SATELLITE DIVISION OF THE INSTITUTE OF NAVIGATION(ION GNSS 2010),THE INSTITUTE OF NAVIGATION,8551 RIXLEW LANE SUITE 360 MANASSAS,VA 20109,USA、24 September 2010（2010-09-24）,1049-1060頁, XP056000217は、リモートセンシング技術をベースとする無線掩蔽（RO）を開示している。
この無線掩蔽（RO）の技術は、電離層内の電子密度を決定するために全地球測位システム（GPS）からの信号を使用し、モデルベース（model-based）のドップラー周波数の推定値およびGPSデータビットの記録(レコード)を使用したオープンループ（OL）追尾を行う。
（0015c）
G. BEYERLE ET AL：“Observation and Simulation of receiver-induced refractivity biases in GPS radio occultation”,JOURNAL OF GEOPHYSICAL RESERCH,Vol.111,no.D12,1 January 2006（2006-01 -01）,XP055158431,ISSN：0148-0227、DOI：10.1029 / 2005J D006673は、GPS電波掩蔽内の受信機に誘発される屈折率の偏り(バイアス)の観測とシミュレーションを開示している。
（0015d）
NIU F ET AL：“GPS Carrier Phase Detrending Methods and Performances for Ionospheric Scintillation Studies”,ITM2012-PROCEEDINGS OF THE 2012 INTERNATIONAL TECHNICAL MEETING OF THE INSTITUTE OF NAVIGATION,THE INSTITUTE OF NAVIGATION,8551, RD LEW,LANE SUITE 360,MANASSAS,VA 20109,USA,1 February 2012（2012-02-01）、1462-1467頁,XP056000934は、GPS搬送波の位相のデトレンド方法および電離層シンチレーションの研究のための講演を開示している。
（0016）
本発明は、さらに無線航法（例えばGNSS）信号から、純粋にオープンループの方法で生成された未処理信号測定値から取得される大気測定値を生成することを目的とする。
（0017）
本発明は、低品質の未処理信号測定値および/または高い大気活動の状況で大気測定値を生成することを目的とする。
（0018）
本発明の一態様によれば、無線航法システムの衛星搭載送信機からの少なくとも１つの無線航法信号に基づいて、大気測定値を生成するための測定システムであって、
オープンループ構成で配置された、データ取得モジュール、復調器モジュールおよび大気監視アルゴリズムモジュールを含み；
前記データ取得モジュールは基準クロックを含み、前記データ取得モジュールは、前記無線航法信号を受信して、そこから複数の中間周波数(IF)サンプル(r)を生成するように適合されており、各サンプルは前記基準クロックから生成される関連時刻タグ(TOW)を有しており；
前記復調器モジュールは、前記IFサンプル(r)と、関連時刻タグ(TOW)と、前記衛星システムに関連する補助データとを受信するよう適合されており、かつそこから相関値(Yi)を生成するように適合されており；
前記大気監視アルゴリズムモジュールは、前記相関値(Yi)を受信して、そこから前記大気測定値を生成するよう適合されている、
測定システムが提供される。
（0018a）
データ取得モジュールは、タグ付きIFサンプル(r)としての各IFサンプル(r)を出力するように適合されてもよい。各タグ付きIFサンプル（r）は、それぞれの時刻タグ(TOW)が付加されたIFサンプル(r)を含む。
（0018b）
復調器モジュールは、上記タグ付きIFサンプル(r)を受信するように適合されてもよい。ここで、復調器モジュール(32)により生成される各相関値(Yi)は、それぞれの時刻タグ(TOW)に関連している。
（0018c）
データ取得モジュールは、IFサンプル(r)を生成するためのアナログ-デジタルコンバータ（変換器；ADC）を含んでもよい。このADCは、基準クロック(418)に結合されており、そのタイミングに基づいてIFサンプル(r)を生成する。
（0018d）
データ取得モジュールは、それぞれのIFサンプル(r)に同期して前記時刻タグ(TOW)を出力するように適合された時刻付加モジュール(420)を含んでもよい。
（0018e）
データ取得モジュールは、基準クロックに結合された時刻付加モジュールを含んでもよい。この時刻付加モジュールは、ローカル時刻を表すローカルカウンタとして作動するように適合されている。このローカルカウンタのカウント値は、各IFサンプル(r)が生成されたときインクリメントされる。
（0018f）
データ取得モジュールは、無線航法信号からアナログIF信号を生成するように適合されたダウンコンバータを含んでもよい。このダウンコンバータは、基準クロックの出力から得られる変換信号に基づいて動作する。
（0018g）
データ取得モジュールは、基準クロックの出力を受信するように結合されたPLLを含んでもよい。ここで、PLLは、ダウンコンバータに変換信号を供給するためのVCOを駆動する。
（0018h）
基準クロックは、無線航法信号の時間フレームに対して、無線航法信号の測距コードのチップの期間の１０分の１未満の整合度を有してもよい。
（0018i）
基準クロックは、約１ナノ秒の精度で時刻タグの生成のために、時刻の推定値をプロパゲート(propagate)することに使用可能な程度の、無線航法信号の時間フレームに対する整合度を有してもよい。
（0018j）
基準クロックは、初期の時刻タグに対応する初期の同期点から、進んだ時刻にプロパゲートするように適合されてもよい。
（0019）
１つの実施形態では、基準クロックは、完全にモデル化されたクロックを含む。ここで、上記時刻タグは、過去の信号同期点から所定のサンプル周期を用いてプロパゲートする時刻信号から得られる。
（0019a）
別の実施形態では、基準クロックは、統制(disciplined)発振器を含む。この統制発振器は、外部の周波数標準から統制クロック信号を受信するように適合された内部発振器を含む。
このシステムは、開始段階とデータ取得段階で動作可能である。ここで、基準クロックは、統制発振器による統制が開始段階では有効となり、データ取得段階で無効となるように動作可能である。外部の周波数基準は、GNSS信号とGPS統制発振器(GPSDO)から供給されてもよい。
（0019b）
別の実施形態では、基準クロックは、未知のパラメータのライブモデリングがクロック推定アルゴリズムを使用して実行される自走クロックを含む。
基準クロックは、推定されたクロックのパラメータを測定し、推定されたクロックパラメータに基づいて初期同期点から時刻タグを正確にプロパゲートするように動作可能であってもよい。
この推定されたクロックパラメータは、第１の衛星搭載送信機のセットからの衛星航法信号から測定されてもよい。また、上記データ取得モジュールにより取得される衛星航法信号は、第２セットの１つ以上の衛星搭載送信機のセットからのものである。第１のセットと第２のセットは、共通の衛星搭載送信機を持っていない。
（0020）
補助データは、受信機に関連するパラメータ（Rec.）、即ち、地球中心、地球固定フレームにおける受信機アンテナの区分的連続軌道を表すパラメータを含んでもよい。
（0021）
補助データは、送信機が搭載された衛星の軌道パラメータ（SV）を含んでもよい。軌道パラメータ（SV）は、ブロードキャストエフェメリス(放送暦；broadcast ephemerides)を含むことができる。ブロードキャストエフェメリスは、GNSSエフェメリスや精密(precise)エフェメリスを含んでもよい。
（0022）
補助データは、エフェメリス情報（ATM）を含んでもよい。
（0023）
復調器モジュールは、時刻タグ(TOW)と受信機に関連するパラメータ(Rec.)を受信し、かつ受信機に関連する時間遅延(δtRX)を出力するためのユーザー受信機モデルを含んでもよい。
（0024）
復調器モジュールは、宇宙飛行体モデル(506)を含んでもよい。宇宙飛行体モデルは、時刻タグ(TOW)および軌道パラメータ（S.V.）を受信し、そして宇宙飛行体に関連する時間遅延(δtSV)を生成するように適合されている。
（0025）
復調器モジュールは、大気モデルを含んでもよい。大気モデルは、時刻タグ（TOW）、受信機に関連する時間遅延(δtRX)、宇宙飛行体に関連する時間遅延(δtSV)、およびエフェメリス情報(Atm.)を受信し、かつ大気に関連する時間差(δtA)を出力するように適合されている。
（0026）
復調器モジュールは、信号に関連する時間遅延(tSIG)を生成するために、受信機に関連する時間差(δtRX)と宇宙飛行体に関連する時間差(δtSV)の和を含む第１の和、および、第１の和と大気に関連する時間差(δtA)との和を含む第２の和を生成するよう適合されている。
（0027）
復調器モジュールは、さらに、コードおよび搬送波MCOを含んでもよい。このコードおよび搬送波MCOは、信号遅延(tSIG)を受信し、DMFへの入力の推定値(Si^)を生成するように適合されている。
（0028）
大気監視アルゴリズムモジュールは、位相プロセス再構成アルゴリズムを含んでもよい。この位相プロセス再構成アルゴリズムは、
以前の位相推定値により現在の相関値(Yi)をデローテート(回転補正：de-rotate)すること、
弁別器(ディスクリミネータ)を用いて、残留位相(Φ)を推定すること、および
以前の位相と現在の残留位相(Φ)の和としての現在の位相(θ)を計算することを含む。
（0029）
大気監視アルゴリズムモジュールは位相差プロセス（処理）アルゴリズムを含んでもよい。この位相差プロセスアルゴリズムは、以下の式を用いて位相差を再構成するように動作可能である。
Figure 2017514146
（0030）
１つの実施形態では、σφはφの値を生成するために、下記の式で与えられるフィルタを用いて、Δ(n)の値から計算されてもよい。
Figure 2017514146
（0031）
１つの実施形態では、σφは下記の関数を用いて、σφ ２から取得されてもよい。
Figure 2017514146
（0032）
本発明の別の態様によれば、無線航法システムの衛星搭載送信機からの少なくとも１つの無線航法信号に基づいて、大気測定値を生成するための測定方法であって、
オープンループ構成で配置された、データ取得モジュールであって基準クロックを含む前記データ取得モジュール、復調器モジュールおよび大気監視アルゴリズムモジュールを提供し；
前記データ取得モジュールを使用して、前記無線航法信号を受信し、そこから各々前記基準クロックから生成される関連時刻タグ(TOW)を有する複数の中間周波数(IF)サンプル(r)を生成し；
前記復調器モジュールを使用して、前記IFサンプル(r)と関連時刻タグ(TOW)と前記衛星システムに関連する補助データを受信し、そこから相関値(Yi)を生成し；
前記大気監視アルゴリズムモジュールを使用して、前記相関値(Yi)を受信し、そこから前記大気測定値を生成することを含む、
測定方法が提供される。
（0033）
本発明の別の態様によれば、処理回路による実行のため、および少なくとも請求項３２におけるステップに対応する、命令を定義するかまたは変換可能なデータが記録または貯蔵された、記録可能、書き換え可能または貯蔵可能な媒体が提供される。
（0034）
本発明の別の態様によれば、通信装置と記憶装置を含み、オンデマンドまたは他の方法で、処理回路による実行のため、および少なくとも請求項３２におけるステップに対応する、命令を定義するかまたは変換可能なデータを送信するように適合されたコンピューターサーバーが提供される。
（0035）
利点は、受信機の時間と位置に関する支援情報の使用、および大気測定値を算出する方法を修正することにより、追尾段階の問題が完全に回避されることで、より強固(ロバスト)な監視システムを構築できることである。
（0036）
さらなる利点は、オープンループモードで動作させ、かつ連続かつ完全な位相推定のための必要性を回避することによって、本発明の実施形態は、弱信号環境での回復力(resilience)、極めて厳しい電離層活動下での強固な追尾、および伝統的な受信機構成よりも高い測定値の利用可能性を提供する。
（0037）
本発明の実施形態について添付図面を参照して説明する。
（0038）
図１（従来技術）は、従来の受信機のデジタル整合フィルタ102-1のブロック図であり、相関値Yi[n]の生成に搬送波位相(θi^)と測距コード(τi^)がどのように使用されるかを示す。
（0039）
図２（従来技術）は、搬送波追尾ループ104および測距コード追尾106の両方のためのループを描いた典型的なクローズドループ追尾構成のブロック図である。
（0040）
図３は、本発明の実施形態に係る測定システムのシステム構成のブロック図である。
（0041）
図４は、図３のデータ取得モジュール31のブロック図である。
（0042）
図５は、図３のオープンループ復調器モジュール32のブロック図である。
（0043）
図６は、図３のAMAモジュール33で行われる一連の相関値Yiが与えられるS4の値の計算処理のブロック図である。
（0044）
図７は、図３のAMAモジュール33で行われる一連の位相推定値θi^が与えられるσφの計算処理のブロック図である。
（0045）
図８は、図３のAMAモジュール33で行われる現在の位相を計算する相関値Yiからの位相の再構成処理のブロック図である。
（0046）
図９は、図３のAMAモジュール33で行われる一連の位相差推定値Δ^が与えられるσφの計算処理のブロック図である。
（0047）
図１０は、図６に示されるアルゴリズムを使用して、方法１と方法２の両方で計算され、かつ市販の大気監視受信機によって提供された基準値と比較したS4を時間毎にプロットして示したものである。
（0048）
図１１は、図７と図９に示されるアルゴリズムを用いて、図３のAMAモジュール33で実行された方法１および方法２の両方で計算され、かつ市販の大気監視受信機によって提供された基準値と比較したσφを時間毎にプロットして示したものである。
（0049）
図３は、本発明の実施形態に係る測定システムのシステム構成のブロック図である。
（0050）
このシステム構成は、３つの主要なブロックである、データ取得（DAQ）31、オープンループ復調器（OLD）32と大気監視アルゴリズム（AMA）33で構成することができる。
（0051）
GNSS受信機におけるフィードフォワードクローズドループ構成に関連する問題および欠点を克服するために、本発明の実施形態は、オープンループ構成を利用している。
図３に示されているように、本発明の実施形態は以下の構成要素を含む。即ち、データ取得モジュール（DAQ）31、オープンループ復調器（OLD）32、および大気監視アルゴリズムのセットを組み込んだ大気監視アルゴリズムモジュール（AMA）33を含む。
以下でより詳細に説明するように、DAQモジュール31は、未処理（例えばGNSS）信号を受信し、それぞれが関連する時刻タグ（TOW）を有するIFサンプル(r)を出力する。復調器32は、rとTOWを受け取り、これら、および補助データ(一般に35で示される)に基づいて、相関値Yjを生成する。相関値YjはAMAモジュール33に出力される。
（0052）
図４は、図３のデータ取得モジュール31のブロック図である。
（0053）
図４に示されているDAQモジュール33は、アンテナ402からの１つ以上のGNSSバンドのための中間周波数サンプル(r)を取得するタスクを実行する。
アンテナ402で受信された信号は、低雑音増幅器(LNA)404で増幅され、バンドパスRFフィルタ406およびさらなる増幅器408を通り、IFに変換するダウンコンバータ410に送られる。
ダウンコンバータ410からのIF出力信号は、ADC416でデジタル(IF)サンプルに変換される前に、ゲインブロック412でさらに増幅され、IFバンドパスフィルタ414を通る。
UTCまたは特定のGNSSシステムの時刻に調整された、ここではTOWと呼ばれる時刻タグを、各IFサンプルrに提供するため、正確なクロック418がサンプルをダウンコンバートし、デジタル化するために使用される。
基準クロック418からのクロック信号422は、PLL424に供給される。PLL424は、VCO426とともにループを形成し、VCO426の出力は、IF信号を生成するために、ダウンコンバータ410に使用される(当業者には理解されるように、IFデータは、リアルタイム処理のためにOLDモジュール32(図３)に直接流されるか、または後処理のためにストレージディスクに流される。)。
DAQの重要な特徴は基準クロック418である。好ましくは、基準クロック418は、UTC、または特定のGNSSシステムの時刻に正確に同期しているか、またはそのオフセット（基準値からの偏差）がよく知られている。
しかしながら、好ましい実施形態では、データ取得中は基準クロック418がGNSSシステムから直接統制されないことが必要である。それが、実施形態のシステムが動作するオープンループの原理に反するからである。DAQ31の基準クロック418の特定の実施形態については、以下に説明する。
（0054）
図３に戻ると、OLDモジュール32は、DAQ31からのIFサンプルおよび関連する時刻タグ(TOWおよびr)、ならびに補助情報35を入力として受け取る；本実施形態において、後者は、それぞれが図３にS.V.、Rec.およびAtm.と示されている衛星エフェメリス、受信機の位置情報、および大気情報を含むが、それらの１つ、いくつか、または全てが使用される。
好ましくは、各サンプルrについて、τi^およびθi^と呼ばれる受信信号パラメータの値を予測するため、関連するTOWは、S.V.、Rec.およびAtm.情報35とともに使用される。
本質的に、この情報は、幾何学的距離、相対論的な影響、任意の公知のおよび/または決定論的な大気の影響、衛星クロックとハードウェアバイアスおよび公知の受信機バイアスからなる伝搬路を予測するために使用されてもよい。
次にこれらの信号パラメータは、相関値Y[ｎ]を生成するためにDMFに渡される。DMFは図１に示されている標準的なDMF102でよい。
このDMF102は、相関値Y[n]を生成するために、搬送波および測距コードのワイプオフ（除去）およびそれに続く積分およびダンプ(dump)操作を実行する。
式（２）で表される位相プロセスを参照すると、そのプロセスθLOS(t)およびθSV Clk.(t)は、大気と衛星エフェメリスの誤差に対応する残留誤差のみを残して除去される。
測距コード遅延又は搬送波位相中の何れか又は両方の特定のオフセットに対応する余分な相関値が生成される可能性があるが、IFデータ内で捕捉された各信号、および各可視衛星のために、少なくとも１つの相関値Y[n]が生成される。
（0055）
AMAモジュール33は、OLDモジュール32により生成される相関値Y[n]を入力として受け取り、それらを様々な大気の測定のために使用する。このモジュールは、大気の状態または大気の活動のレベルに関する情報を提供する様々なアルゴリズムを実装してもよい。
特に、AMAモジュール33は、シンチレーションを記述または定量化するものを含む電離層活動に関連する測定値の生成を行ってもよい。
多くの受信機は、アルゴリズムの特性を利用することにより、相関値Y [n]および完全な位相推定値の両方を使用しているが、AMAモジュール33は、Yiのみの標準電離層測定アルゴリズムを実装してもよい。
（0056）
図４を参照すると、DAQモジュール31は、受信アンテナ402を内蔵し、受信したRF GNSS信号を処理する。このブロックの出力は、正確に時刻タグが付与されたデジタルGNSSデータ（r,TOW）のストリームである。
後者は、それがアンテナによって受信されるとき、各サンプルrは正確な時刻タグTOWと対になっており、本質的にGNSS信号のデジタル表現である。
時刻は、ローカル時間フレーム（受信機）またはリモート時間フレーム（宇宙飛行体、送信機）で表すことができる。いずれの場合も、現在の時刻を示すローカルカウンタは、IF信号の各サンプルrがアナログ-デジタル変換器４１６により記録(record)された時にインクリメントされてよい。このカウンタの値は、このサンプルrと対になって、時刻タグTOWを表す。
（0057）
好ましくは、IFデータを取得し、時刻タグ（TOW）を生成するために十分に正確なクロック418が使用される。より好ましくは、シンチレーションの監視ため、クロック418は、非常に低い位相ノイズを示すことが必要である。
また、本発明の実施形態は、クロックのドリフトおよびドリフトのレートに対して、データ取得動作の全期間にわたって既知のモデルに適合するのに十分に安定なことであるとの第２の要件を課している。
与えられたGNSS時間フレームに対して必要な時刻の整合(アライメント)の度合は、監視された信号の特性により決定され、公称の測距コードチップの１０分の１未満に設定される。複数の実施形態によれば、以下に示すように、異なる実装を用いることができる。
（0058）
例えば、本発明の１つのライブ実装は、GNSSフロントエンドとデジタイザを使用している。
典型的には、このようなシステムは、１つの信号事前調整(pre-conditioning)ブロック（プリアンプとフィルタ）、１つ以上の合成器(シンセサイザー)と混合器(ミキサー)、アンチエイリアスフィルタ、およびアナログ-デジタル変換器（ADC）を含む。
局部発信器を調整するクロックはADCを駆動する必要がある。後のミッション(post-mission)の実装は、GNSSデジタイザで以前に捕捉されたファイルからのストリーミングデータを含む。
（0059）
図４に示すように、GNSSサンプルストリームにタグを付ける処理は、単一の基準クロック418によって駆動される。このクロックは、GNSS時間フレームに整合されること、および非常に高い精度（≒１.−ナノ秒）で、時刻タグを生成するために、時刻の推定値をプロパゲートするために使用できることが重要である。
（0060）
一実施形態では、完全な（すなわち完全にモデル化された）クロック418が使用される。この場合、DAQモジュール31は、予測可能なサンプル周期を用いて、過去の単一同期点からの時刻を単純にプロパゲートする。
（0061）
別の実施形態は統制発振器（DO）に依存することを含む。この場合、DAQモジュール31は、データが取得される期間中に基準発信器がアライメント（整合）を維持することを必要とする。
統制発振器は、一般に内部発振器を含み、内部発振器に対し、修正または操縦（steering）を提供するために、外部の周波数標準を利用する。
外部の周波数標準の一つの共通のソースはGNSS信号である。GNSS(即ちGPS統制発振器(GPSDO；GPS Disciplined Oscillator）の使用）から統制される場合には、データ取得処理を開始する前に統制が行われることを保証しなければならない。
その上、装置は、データ取得処理の全体を通して、アライメントを維持するために十分に良好な持続(ホールドオーバー)性能を有していることが保証されなければならない。
もちろん、このような実装は、標準、既製のGPSDOを基準クロックとして採用し、データ取得期間中の統制動作を無効にすることを選択してもよい。
（0062）
自走クロックが使用される別の実施形態では、その未知のパラメータ（すなわちドリフトとドリフトレート）のライブモデリングは、オーダーメイド(bespoke)のクロック推定アルゴリズムを用いて実行される。
本質的に、これは、GNSSシステムが基準クロックの未知のパラメータの推定に対する基準として使用されるという意味で、前項の実施形態に相当する。
しかし、推定されたクロックパラメータは、発振器を操縦または統制しようとするというよりも、むしろ、単に測定され、より正確な初期同期点からの時刻タグをプロパゲートするために使用される。
この特殊なケースは、以前に取得したデータの後処理である。
このデータは、後処理され、テスト中の大気の異常によって、影響を受けていないか、または軽度の影響を受けただけに見える1つ（またはそれ以上）の衛星信号は、未知のクロックパラメータを推定するために利用することができる。
未知のクロックパラメータを推定するために使用する衛星のセットと、大気を監視するために使用する衛星のセットとは、相互に排他的でなければならないことに留意されたい。
このアプローチは、正確な時刻タグ付け機能を使用することなく、時間的に先に記録されたシンチレーションファイルに対する結果を検証するために使用されている。
（0063）
上述した基準クロック418の３つの実施形態は、初期時刻タグから進んだ時刻に初期同期ポイントをプロパゲートするための手段を表す。全ての場合に、この初期の時刻同期が達成されなければならない。実施形態においては、これはDAQ31の出力rを処理することによって行われる。
これらのサンプルは、従来の典型的なGPSSで行われるように処理することができ、それによって、伝統的な、クローズドループアルゴリズムがサンプルrに適用され、位置と固定時刻（time-fix）が計算される。この固定時刻は時刻タグ処理のための初期同期ポイントを提供する。
（0064）
図５は、図３のオープンループ復調器モジュール32のブロック図である。
（0065）
オープンループ復調器（OLD）モジュール32は、相関値Yiを生成するために、DAQモジュール31によって生成されたIFサンプルrを処理する。
OLDモジュール32に必要な入力は、好ましくは、正確に時刻タグが付加されたデジタルGNSSデータ（r、TOW）のストリーム；ユーザーダイナミックス、宇宙飛行体ダイナミックス（オンボードクロックを含む）、および大気遅延を含む。
受信機側で誘起される遅延は、IFサンプルの時刻タグTOWへの吸収によってモデル化されていることに注意されたい。このモジュールは、それぞれの可視衛星について相関Yi値の出力として生成する。典型的なコヒーレント積分期間は１ミリ秒であるが、より一般的な実装では短いか長い期間を使用することができる。
（0066）
OLD32は、図５に示すように、２つの主要なサブブロックから構成されている：GNSS信号シンセサイザ502とDMF102である。
このシンセサイザは、ユーザーの特定の位置および時刻において受信アンテナ402(図４)で観測される衛星信号を再生する。
それを行うためには、ユーザーのアンテナ402の軌道がゼロまたは正確に知られている必要がある。これは、図５に“Rec.”として示されている、地球中心（earth-centered）、地球固定フレームにおける受信アンテナ402の区分的連続軌道を表す入力を提供することにより達成される。OLD32は、適切に時刻タグ（TOW）と入力Rec.を受信するためのユーザー受信機モデルを含み；ユーザー受信機モジュール504は時間遅延Δ(t)RXを出力する。
（0067）
衛星信号ダイナミックスは、軌道パラメータおよび正確な時刻基準と結合されたクロック補正を使用して生成される。この軌道モデルのパラメータは、図５に変数“S.V.”として示されている。
これらの軌道パラメータは、しばしばエフェメリス（ephemerides）と呼ばれ、GNSS、いわゆるサードパーティ(第三者)精密エフェメリス、または任意の他の適切な軌道モデルのセットによって提供されるように、一連のブロードキャストエフェメリスの形態をとることができる。
典型的な例は、GPS、ガリレオ(Galileo)と北斗(BeiDou)に使用されるケプラーパラメータ、またはGLONASSのために使用されるデカルト位置および派生(デリバティブ)モデルのブロードキャストセットである。他のオプション、国際GNSSサービス（IGS）などのサードパーティが提供する正確なエフェメリスモデルを含む。
OLD32は、適切には、宇宙飛行体モデル506も含み；かつ飛行体モデル506は、時刻タグ(TOW)および軌道パラメータ(S.V.)を受信し、宇宙飛行体に関連する時間遅延δTSVを生成するように適合されている。
最終的には、大気遅延の推定値は、それぞれ別のモデルを使用して生成され、既知の値、または補助エフェメリス情報を用いて適切にパラメータ化される。この補助エフェメリス情報は、図５に変数“Atm.”として示されている。
典型的なモデルは、ほんの数例挙げるだけでも、電離層のためのKlobucharまたはNeQuickモデル、および対流圏のためのSaastamoinenモデルを含む。
大気モデル508は、時刻タグ（TOW）、受信機に関連する時間遅延(δtRX)、宇宙飛行体に関する時間遅延(δtSV)、およびエフェメリス情報(Atm.)を受信するように適合されている。大気モデル508は、これらから大気に関連する時間遅延(δtA)を生成する。
（0068）
OLD32は、信号に関連する遅延(δtRX)と宇宙飛行体に関連する遅延(δtSV)の和からなる第１の和512を生成するための加算要素510を含む。さらに、第２の加算要素514が、時間遅延(δtSIG)を生成するために、第１の和512と大気に関連する遅延(δtA)の和からなる第２の和を生成するように適合されている。この信号に関連する遅延(δtRX)は、DMF102への入力の推定値(Si^)を生成するために、コードおよび搬送波MCO516に入力される。
（0069）
DMF102は、図１に示されているように、シミュレートされたレプリカ信号に対する入力信号ストリームの逆拡散と累積を行う。このステップは、フィードバック動作がないとの違いはあるが、図２に示されている一般的なクローズドループGNSS受信機内で相関器エンジンが行う動作と極めて類似している。
有利なことに、本発明は、これらの異常の影響を受ける相関値Yiを観測するいかなるフィードバックシステムにも依存せず、受信信号において観測可能な振幅、遅延および位相変化に正確に依存する。
DMF102の一実施形態においては、衛星信号毎に唯一の相関器が使用されるが、コンセプトはより一般的であり、任意の数を使用することができる。別の実施形態では、受信信号パラメータの一層の耐久性及び正確な予測を提供するために、公称の周波数および遅延に調整された相関器に対して、周波数および遅延の両方が周辺で一定の間隔を置いた多くの相関器が実装される。各実施形態では、実装は、古典的なリアルタイムまたは後処理を含む。
（0070）
OLD32の内部において、２つのサブブロック(信号シンセサイザ502およびDMF102)の間のインタフェースは、正確に計時されたコード遅延および搬送波位相のシーケンス、それぞれτi^およびθi^である。どちらも、ユーザーと航法衛星のアンテナとの間の見通し線(LOS:Line Of Sight)の幾何学的な距離に関連する。しかし、両者の間には、数ある要素の中で、特に分散性媒体である大気に起因する実質的な差異が存在する。
正確なストリーム時間を使用することで、アキュムレータは何時でもダンプされることが可能であるが、通常は拡散されたコードのエッジに同期してダンプされる。
そのうえ、ブロードキャストされた航法メッセージにアクセス権が与えられると、受信信号に対するビットワイプオフ操作を許容し、航法メッセージの重要な部分を予測することができる。可能な場合、次のセクションで説明するように、それは位相再構成を著しく強化する。
（0071）
大気監視アルゴリズム（AMA）モジュール33（図3）は、複素相関値Y（n）をOLDモジュール32から入力として受け取り、かつ大気測定値を出力として生成する。以下、以後で定義するS4およびσφを含む一般的なシンチレーション指数を生成する電離層シンチレーション監視モジュールの実施形態について論じる。
（0072）
シンチレーション特性のための広く受け入れられているパラメータは、振幅のためのS4と搬送波位相のためのσφである。一般的には、S4の計算には、Iと呼ばれる中間パラメータを相関値Y[n]から計算することが必要である。多くの場合、これは、狭帯域パワー(NBP)と広帯域パワー(WBP)と呼ばれるさらなる中間変数を介して実装される。
あるいは、図６に示されているように、他の実装では、相関値Yiから直接この値を計算する。
Iの値は、一般に0.1Hzの帯域幅を持つ６次のバターワース型のローパスフィルタを用いてデトレンドされる。このローパスフィルタは、数値安定性の理由から、２段のカスケードフィルタとして実装され、図６にLPF62として示されている。
（0073）
図７は、図３のAMAモジュール33で行われる一連の位相推定値θi^が与えられるσφの計算処理のブロック図である。
（0074）
σφを計算する従来の方法は以下の３つのステップを含む。最初に、受信されたGNSS信号の位相は、等価プロセスθ（n]＝θ(nTI)を生成するために、再構成され、固定のサンプルレートでサンプリングされる。
二番目に、それは0.1Hzの帯域幅を持つ６次のバターワース型のハイパスフィルタを用いてデトレンドされる。
最後に、デトレンドされた位相プロセスの変動が、有限の重複しないブロック期間τ秒で計算される(ブロック74)。
したがって、搬送波位相測定が位相異常を測定するために使用されるのに対し、プロンプトレプリカ（またはむしろそれらの副産物）の複素相関器の出力は、振幅異常を評価するために使用される。
（0075）
標準的なGNSS受信機は、PLLのようなクローズドループ位相追尾アルゴリズムを介して、搬送波位相観測値を導出することに留意されたい。対照的に、本発明の実施形態においては、σφを計算するための２つの新規なアプローチが採用されている。
受信信号の位相の追尾をクローズドループで行わないことの利点は；シンチレーション指数の計算を追尾ループのフィルタのアーチファクトから分離すること、及び非常に悪い信号条件下（例えば、サイクルスリップ）での位相追尾障害の問題を防ぐことである。
（0076）
σφの算出方法の両方の実施形態は、再構成された位相が直ちにデトレンドされるという事実を利用する。
注意喚起したように、OLDモジュール32は、信号合成を経由して、見通し線（LOS）ダイナミック、衛星クロックおよび既知の大気の寄与などを含む、搬送波位相に対する既知の決定的寄与のすべてを除去する。残るのは、エフェメリスと大気モデルと局部発振器からのさらなる位相の寄与から生じる残留誤差である。
デトレンドプロセスの目的は、θAtm.(t)以外の全ての要因の寄与を除去することであり、したがって、位相再構成処理に必要なことは、この用語（term）を表すことだけである。２つの実施形態は以下のとおりである。
［第１の方法：位相プロセスの再構成］
（0077）
図８は、図7のAMAモジュール33で行われる、相関値Yiからの位相再構成、即ち、以前の位相と、以前の位相推定値により現在の相関器をデローテート（de-rotating；回転解除）したときに観察される残留位相の和としての現在の位相を算出する処理のブロック図である。
符号ｅ-jθは、実軸に対しての角度−θをもつ複素平面上の単位ベクトルを表し、符号Ｚ-1は１サンプルの遅延を表す。
この場合、位相プロセスの有効部分(useful part)は、連続する対をなす相関値間で計算される位相差の累積和によって再構成される（すなわち、Y[N-1]とY[N]）。
図８に示すように、相関値Yiは、以前の位相推定値によってデローテートされ(ブロック82を介して第１の乗算要素84へ入力する)、かつ残留位相の誤差が測定される。弁別ブロック84は、式（４）で示した位相弁別器の機能を表す。
総位相推定値θは、これらの残留位相測定値を積算することによって与えられる。結果として得られた位相プロセスは、典型的なGNSS受信機によって生成されるのと同様の方法で処理することができる。
搬送波周波数測定および位相測定に最も重要な寄与は、ユーザーと衛星アンテナとの間のLOS速度だけではなく、ユーザーと衛星クロックのドリフトの差に起因する。
従来の受信機からの位相測定値を使用する場合、一般に、この寄与は、無視できない収束時間を犠牲にしてデトレンディングフィルタによって除去される。実際、このようなフィルタのこの収束時間は、数分のオーダーになる可能性がある。この技術は、当然、この位相の寄与を除去する。その結果は、公称の雑音条件では、同じデトレンディングフィルタは、従来の位相プロセスと第１の方法の一つに正確に同じ結果をもたらすが、後者の場合がはるかに速く収束することを示す。
（0078）
詳細には、位相プロセスは以下のように再構成される。まず、現在の相関値が以前の位相推定値によってデローテートされる（最初のサンプルのためにゼロに初期化される）。
二番目に、残留位相φは、現在のサンプルに対応するデータビットの符号が既知であるか否かに応じて、コヒーレント（四象限アークタンジェント）であるかまたは非コヒーレント（アークタンジェント）であるかが弁別器86により推定される。
現在の位相の値φは、次に、以前の位相(ブロック88のアプリケーションにより取得される)と現在の残留位相φの和(加算要素87を介して)として計算される。
数学的には、このプロセスは以下の式（４）によって表される。
Figure 2017514146
 ここで、d[n]は、現在のデータビットの値であり、かつ位相プロセスは、位相差の再帰的な合計に分解されると仮定する。即ち、
θ[n]＝θ[n−1]＋Δ[n]、符号R(x)およびI(x)は、それぞれxの実部および虚部、そして関数atan(x)およびatan2(y,x)は、それぞれアークタンジェントおよび４象限アークタンジェント関数を表す。このσφを生成するための位相推定のプロセスは、次の図７に示されている標準的なアルゴリズムに続く。
［第２の方法：変更されたデトレンディングフィルタ］
（0079）
図９は、図３のAMAモジュール33で行われる、一連の位相差推定値Δ^が与えられσφの値が計算される処理のブロック図である。
この第２のアプローチは、本質的に0.1Hzの帯域幅を持つ６次のハイパスフィルタが搬送波位相の全てのゆっくりと変化するダイナミックスの除去のてこ入れをする。このことは、位相プロセスの再構成を完全に回避し、絶対的な位相に代えて位相差が用いられ得ることを示唆している。
デトレンディングフィルタの次数を１低減することができ、従来のアプローチに従って、等価でより簡単なフィルタが、全く同じ結果を生成する。
この新しいフィルタは以下のように見出される。
下記の式（５）で示されるｚ領域の伝達関数を有する伝統的な６次のバターワース型ハイパスフィルタを想定する。
Figure 2017514146
 ここで、式（５）の出力はデトレンドされた位相プロセスであり、ここにθＤ[n]で示されるものは、再構成された位相θ[n]にH(z)を適用することで見出される。このフィルタは、フィルタB(z)をθ(t)に対して適用し、次いでその結果に対して1/A(z)を適用するカスケードプロセスと完全に等価である。したがって、この位相プロセスθ[n]は、B(z)によりのみ行われる。
（0080）
エポックnでは、ここにθB[n]として示されるθ[n]に対して適用されたフィルタB（z）の出力は、
Figure 2017514146
 で与えられる。
（0081）
また、その位相プロセスを位相差の総和に分解することを想定すると、式（６）は、
Figure 2017514146
 により与えられる。
（0082）
ここでH(z)はハイパスフィルタであり、直流ゲインが０である。そうするとH(1)＝０であるから、
Figure 2017514146
 であり、式（７）は、
Figure 2017514146
 に削減される。
ここで、
Figure 2017514146
 式（８）を思い起こすと、c６＝０であるから、
Figure 2017514146
（0083）
したがって、フィルタH(z)を用いて、元の位相プロセスθを処理することにより得られたデトレンディング位相プロセスθDは、
Figure 2017514146
 で与えられる新しいフィルタを用いて、位相差プロセスΔ[n]＝θ[n]−θ[n−1]を処理することと完全に等価であることがわかる。
（0084）
式（１２）は、位相プロセスの再構成が全く必要ないこと、および２つの隣接する相関値の間の位相差の推定のみが必要であることを示唆する。このことは、σφの生成を大幅に簡素化し、結果として、著しくより強固な監視アルゴリズムになる。具体的には、一実施形態では、位相差プロセスの再構成は、
Figure 2017514146
 を含む。
ここで、相関値の両方のビットが既知か否かに依存して、２つの位相差推定値の一方が使用される。この位相差プロセスの生成の後、変更されたフィルタが適用可能であり、最後にσφの推定値を計算することができる。
Δ(n)からσφを計算するこのプロセスは、図９に示されている。この図において、“M-FLT”と名付けられたブロックが、ここで述べた変更されたフィルタを表している。
（0085）
式（１３）で与えられるような位相差の単純な直接的推定は、以下に想定されるが、代替的な実施形態では、Δ[n]の推定値の精度を高めるために、相関値のさらなる処理を行ってもよい。
例えば、この推定を、位相プロセスのためのアプリオリ（演繹的な）統計モデルを条件としてもよい。即ち、外れ値または失敗と思われる推定値；即ち、３つまたはそれより多い一連の推定値が観測されるような推定値、を拒否するかまたは修正する仕組みを実装してもよい。
（0086）
S4およびσφについて、第１の方法、第２の方法より計算された値を比較したもの、および市販の大気監視(基準)受信機をそれぞれ図１０および図１１に示す。これらより、オープンループアルゴリズムは、伝統的なクローズドループ受信機のそれを完全に再生可能であることが判る。
（0087）
それぞれの実施において種々の構成要素を有する様々な実施形態を参照して、実施形態について説明してきたが、他の実施形態では、これらの構成要素および他の構成要素の組合せおよび置換を利用することが理解されるであろう。
（0088）
また、実施形態のいくつかは、コンピュータシステムのプロセッサによってまたは機能を実行する他の手段によって実装可能な方法または方法の要素の組合せとして本明細書中に記載されている。したがって、このような方法または方法の要素を実行するために必要な命令を有するプロセッサは、方法または方法の要素を実行するための手段を構成する。さらに、装置の実施形態の本明細書に記載の要素は、本発明を実施する目的で要素によって実行される機能を実行するための手段の一例である。
（0089）
本明細書で提供される説明において、多数の特定の細部が記載されている。しかしながら、本発明の実施形態はこれらの特定の細部なしに実施できることが理解される。さらに、周知の方法、構造、および技術は、この説明に対する理解を不明瞭にしないために詳細には示されていない。
（0090）
このように、本発明の好ましい実施形態であると考えられるものを説明してきたが、当業者は、他の変形または更なる変形が本発明の精神および範囲から逸脱することなくなされ得ることを理解するであろう。また、本発明の範囲内に入るような全ての変更及び変形は、請求することを意図したものである。例えば、上記の任意の式は、使用可能な手順の単なる代表例である。
ブロック図において機能を追加または削除することができるし、機能ブロック間で動作を交換することもできる。また、記載した方法において、本発明の範囲内でステップを追加または削除することができる。

Claims (35)

  1. 無線航法システムの衛星搭載送信機からの少なくとも１つの無線航法信号に基づいて、大気測定値を生成するための測定システムであって、
    オープンループ構成で配置された、データ取得モジュール(31)、復調器モジュール（32）および大気監視アルゴリズムモジュール（33）を含み；
    前記データ取得モジュール（31）は基準クロック(418)を含み、前記データ取得モジュール(31)は、前記無線航法信号を受信して、そこから複数の中間周波数(IF)サンプル(r)を生成するように適合されており、各IFサンプルは前記基準クロック(418)から生成される関連時刻タグ(TOW)を有しており；
    前記復調器モジュール（32）は、前記IFサンプル(r)と、関連時刻タグ(TOW)と、前記衛星システムに関連する補助データ(35)とを受信するように適合されており、かつそこから相関値(Yi)を生成するように適合されており；
    前記大気監視アルゴリズムモジュールは、前記相関値(Yi)を受信して、そこから前記大気測定値を生成するよう適合されている、
    測定システム。
  2. 請求項１に記載された測定システムにおいて、
    前記データ取得モジュール（31）は、各IFサンプル(r)をタグ付きIFサンプル(r)として出力するように適合されており、各タグ付きIFサンプル(r)は、それぞれの時刻タグ(TOW)が付加されたIFサンプル(r)を含む、測定システム。
  3. 請求項１または２に記載された測定システムにおいて、
    前記復調器モジュール（32）は、前記タグ付きIFサンプル(r)を受信するように適合されており、それによって前記復調器モジュール(32)により生成された各相関値(Yi)はそれぞれの時刻タグ(TOW)に関連付けられる、測定システム。
  4. 先行する請求項のいずれかに記載された測定システムにおいて、
    前記データ取得モジュール（31）は、IFサンプル(r)を生成するアナログ-デジタルコンバータ(ADC)(416）を含み、前記ADCは、前記基準クロック(418)に結合され、そのタイミングで前記IFサンプル(r)を生成する、測定システム。
  5. 先行する請求項のいずれかに記載された測定システムにおいて、
    前記データ取得モジュール（31）は、それぞれのIFサンプル(r)に同期して前記時刻タグ(TOW)を出力するように適合された時刻付加モジュール(420)を含む、測定システム。
  6. 先行する請求項のいずれかに記載された測定システムにおいて、
    前記データ取得モジュール(31)は、前記基準クロック(418)に結合された時刻付加モジュール(420)を含み、時刻付加モジュール(420)は、ローカル時刻を表すローカルカウンタとして動作するように適合されており、前記ローカルカウンタのカウント値は各IFサンプル(r)が生成されたときインクリメントされる、測定システム。
  7. 先行する請求項のいずれかに記載された測定システムにおいて、
    前記データ取得モジュール(31)は、前記無線航法信号からアナログIF信号を生成するように適合されたダウンコンバータ(410)を含み、前記ダウンコンバータ(410)は、基準クロック(418)の出力から得られた変換信号に基づいて動作する、測定システム。
  8. 請求項７に記載された測定システムにおいて、
    前記データ取得モジュール(31)は、前記基準クロック(418)の出力を受信するように結合されたPLL(424)を含み、前記PLL(424)は、前記ダウンコンバータ(410)に変換信号を供給するVCO(426)を駆動する、測定システム。
  9. 先行する請求項のいずれかに記載された測定システムにおいて、
    前記基準クロック(418)は、無線航法信号の時間フレームに対して、前記無線航法信号の測距コードのチップの期間の１０分の１未満の整合度を有する、測定システム。
  10. 先行する請求項のいずれかに記載された測定システムにおいて、
    前記基準クロック(418)は、約１ナノ秒の精度の時刻タグの生成のために、時刻の推定値をプロパゲートすることに使用可能な程度の無線航法信号の時間フレームに対する整合度を有する、測定システム。
  11. 先行する請求項のいずれかに記載された測定システムにおいて、
    前記基準クロック(418)は、初期の時刻タグに対応する初期の同期点から進んだ時刻にプロパゲートするように適合されている、測定システム。
  12. 先行する請求項のいずれかに記載された測定システムにおいて、
    前記基準クロック(418)は、完全にモデル化されたクロックを含み、前記時刻タグは、所定のサンプル周期を用いて過去の信号同期点からプロパゲートする時刻信号から生成される、測定システム。
  13. 請求項１乃至１１のいずれかに記載された測定システムにおいて、
    前記基準クロック(418)は、統制発振器を含み、前記統制発振器は内部発振器を含み、統制クロック信号を外部の周波数標準から受信するように適合されている、測定システム。
  14. 請求項１３に記載された測定システムにおいて、
    開始段階およびデータ取得段階で動作可能であり、前記基準クロック(418)は、前記統制発振器による統制が、前記開始段階中は有効であり、前記データ取得段階中は無効となるように動作可能である、測定システム。
  15. 請求項１４に記載された測定システムにおいて、
    前記外部の周波数標準は、GNSS信号とGPS統制発振器(GPSDO)信号のいずれか一方によって提供される、測定システム。
  16. 請求項１乃至１１のいずれかに記載された測定システムにおいて、
    前記基準クロック(418)は、未知のパラメータのライブモデリングがクロック推定アルゴリズムを使用して実行される自走クロックを含む、測定システム。
  17. 請求項１６に記載された測定システムにおいて、
    前記基準クロック(418)は、
    推定されたクロックパラメータを測定し、
    前記推定されたクロックパラメータに基づく初期同期点から正確に時刻タグをプロパゲートするように動作可能である、測定システム。
  18. 請求項１７に記載された測定システムにおいて、
    前記推定されたクロックパラメータは、第１の衛星搭載送信機のセットからの衛星航法信号から測定され、かつ前記データ取得モジュール(31)により取得される前記衛星航法信号は、第２の衛星搭載送信機のセットの１つ以上の衛星搭載送信機からのものであり、前記第１のセットおよび第２のセットは共通の衛星搭載送信機を備えていない、測定システム。
  19. 先行する請求項のいずれかに記載された測定システムにおいて、
    前記補助データ(35)は、受信機に関連するパラメータ(Rec.)を含み、前記受信機に関連するパラメータは、地球中心、地球固定フレームにおける受信機アンテナ（402）の区分的連続軌道を表す、測定システム。
  20. 先行する請求項のいずれかに記載された測定システムにおいて、
    前記補助データ(35)は、前記送信機が搭載された衛星の軌道パラメータ(S.V.)を含む、測定システム。
  21. 請求項２０に記載された測定システムにおいて、
    前記軌道パラメータ(S.V.)は、ブロードキャストエフェメリスを含む、測定システム。
  22. 請求項２１に記載された測定システムにおいて、
    前記ブロードキャストエフェメリスは、GNSSエフェメリスまたは精密エフェメリスを含む、測定システム。
  23. 先行する請求項のいずれかに記載された測定システムにおいて、
    前記補助データ(35)は、エフェメリス情報(ATM)を含む、測定システム。
  24. 先行する請求項のいずれかに記載された測定システムにおいて、
    前記復調器モジュール(32)は、時刻タグ(TOW)および前記受信機に関連するパラメータ(Rec.)を受信するため、ならびに受信機に関連する時間遅延(δtRX)を出力するためのユーザー受信機モデル(504)を含む、測定システム。
  25. 先行する請求項のいずれかに記載された測定システムにおいて、
    前記復調器モジュールは、宇宙飛行体モデル(506)を含み、前記宇宙飛行体モデルは、時刻タグ(TOW)および前記軌道パラメータ(S.V.)を受信し、宇宙飛行体に関連する時間遅延(δtSV)を生成するように適合されている、測定システム。
  26. 請求項２４に従属する請求項２５に記載された測定システムにおいて、
    前記復調器モジュール(32)は、大気モデル(508)を含み、前記大気モデルは、前記時刻タグ(TOW)、前記受信機に関連する時間遅延(δtRX)、前記飛行体に関連する時間遅延(δtSV)、および前記エフェメリス情報(Atm.)を受信し、大気に関連する時間差(δtA)を出力するように適合されている、測定システム。
  27. 請求項２４に従属する請求項２５に記載された測定システムにおいて、
    前記復調器モジュール(32)は、前記受信機に関連する時間差(δtRX)と前記宇宙飛行体に関連する時間差(δtSV)の和を含む第１の和を生成し、かつ信号に関連する時間遅延(tSIG)を生成するために、前記第１の和(512)と前記大気に関連する時間差(δtA)との和を含む第２の和を生成するように適合されている、測定システム。
  28. 請求項２７に記載された測定システムにおいて、
    前記復調器モジュールは、前記信号遅延(tSIG)を受信し、DMF(102)の入力のための推定値(Si^)を生成するように適合されているコードおよび搬送波MCOをさらに含む、測定システム。
  29. 先行する請求項のいずれかに記載された測定システムにおいて、
    大気監視アルゴリズムモジュール(33)は、位相プロセス再構成アルゴリズムを含み、前記位相プロセス再構成アルゴリズムは、
    以前の位相推定値により現在の相関値(Yi)をデローテートすること；
    弁別器を使用して残留位相(Φ)を推定すること；および
    以前の位相と現在の残留位相(Φ)との和としての現在の位相(θ)を計算すること；
    を含む、測定システム。
  30. 請求項１乃至２８のいずれかに記載された測定システムにおいて、
    前記大気監視アルゴリズムモジュール(33)は、位相差プロセスアルゴリズムを含み、前記位相差プロセスアルゴリズムは、
    Figure 2017514146
     を用いて、位相差を再構成するように動作可能である、測定システム。
  31. 請求項３０に記載された測定システムにおいて、
    σφは、
    Figure 2017514146
     で与えられるφの値を生成するためのフィルタを用いることでΔ[n]の値から計算される、測定システム。
  32. 請求項３１に記載された測定システムにおいて、
    σφは、関数
    Figure 2017514146
     を用いて、σφ ２から得られる、測定システム。
  33. 無線航法システムの衛星搭載送信機からの少なくとも１つの無線航法信号に基づいて、大気測定値を生成するための測定方法であって、
    オープンループ構成で配置された、データ取得モジュール(31)であって基準クロック(418)を含むデータ取得モジュール(31)、復調器モジュール（32）および大気監視アルゴリズムモジュール（33）を提供し；
    前記データ取得モジュール(31)を使用して前記無線航法信号を受信し、そこから、各IFサンプルが前記基準クロック(418)から生成される関連時刻タグ(TOW)を有する、複数の中間周波数(IF)サンプル(r)を生成し；
    前記復調器モジュール(32)を使用して、前記IFサンプル(r)と、関連時刻タグ(TOW)と、前記衛星システムに関連する補助データ(35)とを受信し、そこから相関値(Yi)を生成し；かつ
    前記大気監視アルゴリズムモジュール(33)を使用して、前記相関値(Yi)を受信し、そこから前記大気測定値を生成する、
    測定方法。
  34. 処理回路による実行のため、および少なくとも請求項３２におけるステップに対応する命令を定義するかまたは変換可能なデータが記録または貯蔵された、記録可能、書き換え可能または貯蔵可能な媒体。
  35. 通信装置と記憶装置を含み、オンデマンドまたは他の方法で、処理回路による実行のため、および少なくとも請求項３２におけるステップに対応する命令を定義するかまたは変換可能なデータを送信するように適合されたコンピューターサーバー。
JP2017501498A 2014-03-28 2015-03-27 大気監視のための無線航法信号を処理するための方法および装置 Active JP6513179B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP14162381.9 2014-03-28
EP14162381.9A EP2924468A1 (en) 2014-03-28 2014-03-28 Method and apparatus for processing radionavigation signals for atmospheric monitoring
PCT/EP2015/056795 WO2015144914A1 (en) 2014-03-28 2015-03-27 Method and apparatus for processing radionavigation signals for atmospheric monitoring

Publications (2)

Publication Number Publication Date
JP2017514146A true JP2017514146A (ja) 2017-06-01
JP6513179B2 JP6513179B2 (ja) 2019-05-15

Family

ID=50391051

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017501498A Active JP6513179B2 (ja) 2014-03-28 2015-03-27 大気監視のための無線航法信号を処理するための方法および装置

Country Status (11)

Country Link
US (1) US10663601B2 (ja)
EP (2) EP2924468A1 (ja)
JP (1) JP6513179B2 (ja)
KR (1) KR102299105B1 (ja)
CN (1) CN106415314B (ja)
AU (1) AU2015238231B2 (ja)
BR (1) BR112016022524B1 (ja)
CA (1) CA2940330C (ja)
ES (1) ES2713148T3 (ja)
RU (1) RU2680711C2 (ja)
WO (1) WO2015144914A1 (ja)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SG11201701730WA (en) * 2014-09-05 2017-04-27 Fugro Intersite B V Ionospheric scintillation prediction
CN106707306B (zh) * 2016-12-09 2019-04-02 北京无线电计量测试研究所 一种gnss接收机载波相位测量值错误检测方法
RU2664019C9 (ru) * 2017-04-10 2019-08-15 Ашот Эрнстович Кочарян Устройство и способ адаптивной компенсации искажений и восстановления несущей сигнала для когерентных приёмников
JP6906683B2 (ja) * 2017-08-01 2021-07-21 ヴィアサット, インコーポレイテッドViaSat, Inc. 地理的に分布した衛星アクセスノードにおけるシンチレーション緩和
KR102450580B1 (ko) 2017-12-22 2022-10-07 삼성전자주식회사 금속 배선 하부의 절연층 구조를 갖는 반도체 장치
US11506796B2 (en) * 2018-11-21 2022-11-22 Unicore Communications, Inc. Method, apparatus and mobile device for extending real-time kinematic positioning during reference data outage
US11038606B1 (en) * 2019-04-19 2021-06-15 The Government of the Unites States of Ameroica as represented by the Secretary of the Air Force Source discrimination of fluctuating radio signals
US20200363536A1 (en) * 2019-05-16 2020-11-19 Orolia Usa Inc. Methods for enhancing non-global navigation satellite system location and timing pseudorange positioning calculations and systems thereof
RU2713211C1 (ru) * 2019-07-23 2020-02-04 Ашот Эрнстович Кочарян Устройство и способ приема оптического сигнала
CN111339676A (zh) * 2020-03-10 2020-06-26 中国科学院国家空间科学中心 一种近地空间环境综合数据分析系统
JP2022104765A (ja) * 2020-12-29 2022-07-11 財團法人工業技術研究院 信号センシングのためのシステムと方法
RU2769113C1 (ru) * 2021-01-11 2022-03-28 Федеральное государственное автономное образовательное учреждение высшего образования "Сибирский федеральный университет" Способ коррекции погрешности дискретизации дальномерного кода
KR102580909B1 (ko) 2021-08-24 2023-09-20 세종대학교산학협력단 위성항법 반송파 측정치의 사이클 슬립 검출장치 및 그 방법
WO2023224508A1 (en) * 2022-05-16 2023-11-23 Limited Liability Company "Topcon Positioning Systems" Method and apparatuses of measuring phase ionosphere scintillations
CN115987309B (zh) * 2023-03-15 2023-06-20 湖南卫导信息科技有限公司 一种非相干多路干扰信号模拟器及方法
CN116338735A (zh) * 2023-05-11 2023-06-27 天津云遥宇航科技有限公司 一种基于巴特沃斯滤波的电离层掩星闪烁指数s4计算方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11304903A (ja) * 1998-04-17 1999-11-05 Toshiba Corp シインチレーション監視装置
JP2006349672A (ja) * 2005-05-20 2006-12-28 Nec Toshiba Space Systems Ltd Ｇｐｓｒ多周計測型電離層遅延補正装置、ｇｐｓｒ多周計測型電離層遅延補正方法及びプログラム
JP2012137448A (ja) * 2010-12-28 2012-07-19 Japan Radio Co Ltd 電離層遅延評価装置および航法装置

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4894662A (en) * 1982-03-01 1990-01-16 Western Atlas International, Inc. Method and system for determining position on a moving platform, such as a ship, using signals from GPS satellites
US5594453A (en) * 1994-11-01 1997-01-14 Trimble Navigation, Ltd GPS receiver having a rapid acquisition of GPS satellite signals
US6633255B2 (en) * 1995-10-09 2003-10-14 Qualcomm Inc. Method for open loop tracking GPS signals
US6691031B2 (en) * 2002-05-31 2004-02-10 Magellan Corporation Method and apparatus for substituting sensor data for satellite signal in GPS receiver
US7912422B2 (en) * 2006-07-31 2011-03-22 University Corporation For Atmospheric Research Method and system for demodulation of open-loop GPS radio occultation signals
CN101424734B (zh) 2007-10-31 2012-01-04 中国科学院微电子研究所 降低全球定位系统接收机的晶振成本的方法
CN103064091B (zh) * 2011-10-20 2015-02-11 神讯电脑(昆山)有限公司 定位装置与其信号处理方法
US9910158B2 (en) * 2012-12-28 2018-03-06 Trimble Inc. Position determination of a cellular device using carrier phase smoothing
RU2523912C1 (ru) * 2013-02-20 2014-07-27 Межрегиональное общественное учреждение "Институт инженерной физики" Устройство пеленгации исскуственных ионосферных образований
RU2560525C1 (ru) * 2014-06-25 2015-08-20 Александр Васильевич Тертышников Способ определения положения эпицентральной зоны источника и скорости распространения перемещающихся ионосферных возмущений

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11304903A (ja) * 1998-04-17 1999-11-05 Toshiba Corp シインチレーション監視装置
JP2006349672A (ja) * 2005-05-20 2006-12-28 Nec Toshiba Space Systems Ltd Ｇｐｓｒ多周計測型電離層遅延補正装置、ｇｐｓｒ多周計測型電離層遅延補正方法及びプログラム
JP2012137448A (ja) * 2010-12-28 2012-07-19 Japan Radio Co Ltd 電離層遅延評価装置および航法装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
LULICH T D; GARRISON J L; HAASE J S; Y-M YANG; VOO J; XIE F; MURADYAN P: "OPEN LOOP TRACKING OF RADIO OCCULTATION SIGNALS FROM AN AIRBORNE PLATFORM", GNSS 2010 - PROCEEDINGS OF THE 23RD INTERNATIONAL TECHNICAL MEETING OF THE SATELLITE 以下備考, JPN5017002211, 24 September 2010 (2010-09-24), US, pages 1049 - 1060, XP056000217, ISSN: 0003946487 *
NIU F; MORTON Y; WANG J; PELGRUM W: "GPS CARRIER PHASE DETRENDING METHODS AND PERFORMANCES FOR IONOSPHERE SCINTILLATION STUDIES", ITM 2012 - PROCEEDINGS OF THE 2012 INTERNATIONAL TECHNICAL MEETING OF THE INSTITUTE OF NAVIGATION, JPN5017002213, 1 February 2012 (2012-02-01), US, pages 1462 - 1467, XP056000934, ISSN: 0003946489 *
WANG K-N; MURADYAN P; GARRISON J L; HAASE J S; MURPHY B; ACIKOZ U; LULICH T: "OPEN-LOOP TRACKING OF RISING AND SETTING GNSS RADIO-OCCULTATION SIGNALS FROM AN 以下備考", 2013 IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM - IGARSS, JPN5017002212, 21 July 2013 (2013-07-21), pages 3367 - 3370, XP032558256, ISSN: 0003946488, DOI: 10.1109/IGARSS.2013.6723550 *
YORK J; MUNTON D; MAYO M; SUMP T: "DEVELOPMENT OF A PROTOTYPE TEXAS IONOSPHERIC GROUND RECEIVER (TIGR)", ITM 2012 - PROCEEDINGS OF THE 2012 INTERNATIONAL TECHNICAL MEETING OF THE INSTITUTE OF NAVIGATION, JPN5017002210, 1 February 2012 (2012-02-01), US, pages 1526 - 1556, XP056000936, ISSN: 0003946486 *

Also Published As

Publication number Publication date
CN106415314A (zh) 2017-02-15
ES2713148T3 (es) 2019-05-20
JP6513179B2 (ja) 2019-05-15
AU2015238231A1 (en) 2016-10-20
EP3123201B1 (en) 2019-01-30
RU2680711C2 (ru) 2019-02-26
BR112016022524A2 (ja) 2017-08-22
CA2940330C (en) 2022-01-04
BR112016022524B1 (pt) 2022-10-04
RU2016142308A3 (ja) 2018-08-28
RU2016142308A (ru) 2018-05-03
EP2924468A1 (en) 2015-09-30
AU2015238231B2 (en) 2019-08-01
CN106415314B (zh) 2019-06-18
US10663601B2 (en) 2020-05-26
KR20160147775A (ko) 2016-12-23
WO2015144914A1 (en) 2015-10-01
CA2940330A1 (en) 2015-10-01
EP3123201A1 (en) 2017-02-01
US20170139050A1 (en) 2017-05-18
KR102299105B1 (ko) 2021-09-07

Similar Documents

Publication Publication Date Title
JP6513179B2 (ja) 大気監視のための無線航法信号を処理するための方法および装置
JP7066187B2 (ja) 低地球軌道衛星通信システムにおける位置推定
CN108646543B (zh) 一种具有高稳定度性能的gnss驯服时钟方法
O'Hanlon et al. CASES: A smart, compact GPS software receiver for space weather monitoring
JP3657004B2 (ja) 衛星クロックの無線リンク制御
BR112013017960B1 (pt) Método e sistema para determinar correções de relógio de satélite
JP6349362B2 (ja) 無線通信ユニット、集積回路、衛星通信システム及び電離層群遅延の補正方法
US20150042511A1 (en) System and method for continuous carrier wave reconstruction
KR20170099343A (ko) Gnss 수신기에 의해 널 영역 내의 신호를 프로세싱 및 검출하는 방법 및 시스템
JP4848146B2 (ja) 測位信号を送信するための装置、その装置を備える測位システムおよび測位信号を送信するシステム
Dong et al. Implementation and verification of a software-based IF GPS signal simulator
CN106338748A (zh) 基于卡尔曼滤波的gps接收机跟踪环路
JP2008249402A (ja) 周波数間バイアス推定装置及び周波数間バイアス推定方法
Curran et al. An open-loop vector receiver architecture for GNSS-based scintillation monitoring
CN116224746A (zh) 一种星地原子钟融合的高稳定时间基准建立方法
CN115903444A (zh) 分布式星座时间基准建立的本地时间生成方法及装置
Stevanovic et al. Coasting through wideband interference events using robust carrier phase tracking
EP3748401A1 (en) System and method for position determination of a stationary gnss receiver using a distributed time signal
Lashley et al. Vector processing
Gogoi et al. NavIC receiver clock offsets estimation with common view master clock method
Kanouj et al. A new approach to designing a positioning system for operation under the conditions of dynamic stress and high noise level
Vani et al. Evaluation of GPS data for navigational solution and error reduction using kalman filter
Zhao et al. GPS Signal Tracking and Positioning through Long Duration Wide Band Interference with Multiple Model Kalman Filter
Zhao et al. Inertial Aiding for Moving GPS Receiver Carrier Phase Estimation Through Wideband Interference
Abbasiannik Multichannel dual frequency GLONASS software receiver in combination with GPS L1 CA:(URL: httpwww. geomatics. ucalgary. caresearchpublications).

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180202

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20181031

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20181106

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20181212

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20181225

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190325

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: 20190402

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190409

R150 Certificate of patent or registration of utility model

Ref document number: 6513179

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250