JP5553909B2 - サイクルスリップ検出装置、サイクルスリップ検出のための集積回路、同集積回路を備えるgnss受信装置、同gnss受信装置を備える情報端末機器、サイクルスリップ検出方法及びプログラム - Google Patents

サイクルスリップ検出装置、サイクルスリップ検出のための集積回路、同集積回路を備えるgnss受信装置、同gnss受信装置を備える情報端末機器、サイクルスリップ検出方法及びプログラム Download PDF

Info

Publication number
JP5553909B2
JP5553909B2 JP2012548686A JP2012548686A JP5553909B2 JP 5553909 B2 JP5553909 B2 JP 5553909B2 JP 2012548686 A JP2012548686 A JP 2012548686A JP 2012548686 A JP2012548686 A JP 2012548686A JP 5553909 B2 JP5553909 B2 JP 5553909B2
Authority
JP
Japan
Prior art keywords
parity check
cycle slip
result
received data
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.)
Expired - Fee Related
Application number
JP2012548686A
Other languages
English (en)
Other versions
JPWO2012081275A1 (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.)
Furuno Electric Co Ltd
Original Assignee
Furuno Electric Co Ltd
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 Furuno Electric Co Ltd filed Critical Furuno Electric Co Ltd
Priority to JP2012548686A priority Critical patent/JP5553909B2/ja
Publication of JPWO2012081275A1 publication Critical patent/JPWO2012081275A1/ja
Application granted granted Critical
Publication of JP5553909B2 publication Critical patent/JP5553909B2/ja
Expired - Fee Related 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/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
    • G01S19/44Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method

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

本発明は、GPS衛星からの受信信号において発生するサイクルスリップを検出する技術に関する。
従来、GPSによる測位において、例えば、数センチ程度の高い測位精度が要求される場合に、搬送波位相を用いた相対測位が行われる。固定点に設置された基準局と移動体に搭載された移動局で同一のGPS信号を受信することによって前記移動体の位置を算出する移動体の相対測位においては、
同一のGPS衛星から送信されるGPS信号の搬送波位相を基準局と移動局とで同時に積算することが必要である。GPS信号の搬送波位相の積算中にGPS信号が途切れた場合、GPS信号を追尾している受信機の発振器の発振周波数をカウントすることにより前記積算が継続されるが、この発振器がいわゆるフリーラン状態になるため、中断の長さやそのときの状況によっては、前記搬送波位相の積算値がずれを起こす。これをサイクルスリップという。観測量にサイクルスリップが生じると測位結果に誤差が生じる。なお、このサイクルスリップは1つの受信機で測位される単独測位においても発生する。
このサイクルスリップを適切に検出する方法として、様々な方法が提案されている。
一つの従来例によれば、受信機と慣性装置とから算出された速度差を利用してサイクルスリップを検出する方法が開示されている(例えば、特許文献1参照)。ここでは、受信機で計測された移動局の速度V1及び位置P1と、慣性装置(INS)で計測された移動局の速度V2及び位置P2とをフィルタリングして、より精度のよい補正された移動局の速度V3及び位置P3を計算し、外部に出力する。そして、速度V3(またはV2)を受信機に与え、また速度V2に含まれる誤差ΔV2=V2−V3を演算して、慣性装置に与える。GPS受信機は、速度V1を補正された速度V3(またはV2)と比較して、サイクルスリップの検出とその補正を行い、慣性装置は、速度誤差ΔV2を用いて速度V2及び位置P2を補正する。
特開2000−65593号公報
特許文献1に開示されたサイクルスリップ検出方法では、受信機及び慣性装置により計測される2つの観測量の差分を基に検出するため、閾値の設定によっては誤判断の可能性がある。
そこで、本発明は、簡単な構成で、精度の高いサイクルスリップの検出を実現することを目的とする。
上記目的を達成するため、本発明の一つの観点によれば、データ取得部と、パリティチェック演算部と、サイクルスリップ判定部と、を備えるサイクルスリップ検出装置が提供される。データ取得部は、測位用衛星からの受信データを取得する。パリティチェック演算部は、受信データに対しパリティチェックを行なう。サイクルスリップ判定部は、パリティチェック演算部によるパリティチェック結果に基づき受信データにおけるサイクルスリップの有無を判定する。
ここでは、パリティチェック演算部によりパリティチェックを行い、そのパリティチェック結果に基づき受信データにおけるサイクルスリップの有無を判定するため、従来に比して簡単且つ正確にサイクルスリップの検出が可能である。
上記サイクルスリップ検出装置は、次のような構成であってもよい。パリティチェック演算部は、第1パリティチェック演算部と第2パリティチェック演算部とを有する。第1パリティチェック演算部は、パリティチェックにより第1のパリティチェック結果である値を取得する。第2パリティチェック演算部は、第1のパリティチェック結果の値と、受信データの各ビット位置においてサイクルスリップが発生した場合のパリティチェック結果である第2のパリティチェック結果の値とを比較する。そして、サイクルスリップ判定部は、第2パリティチェック演算部による上記比較の結果、第1のパリティチェック結果の値が、上記各ビット位置で取得した第2パリティチェック結果の値のいずれかと一致すると判断した場合、受信データにおいてサイクルスリップが有ると判定する。
また、上記サイクルスリップ検出装置は、データ取得部により取得された受信データにおけるサイクルスリップの検出タイミングを推定するサイクルスリップ推定部を更に備えていてもよい。この場合、パリティチェック演算部は、パリティチェックにより第1のパリティチェック結果である値を取得する第1パリティチェック演算部と、サイクルスリップの検出タイミングに応じた第2のパリティチェック結果である値を取得する第2パリティチェック演算部と、を有し、サイクルスリップ判定部は、第1のパリティチェック結果である値と第2のパリティチェック結果である値とを比較し、一致する場合は受信データにおいてサイクルスリップが有ると判定する。
更に、上記サイクルスリップ検出装置においては、第1パリティチェック演算部は、所定の式を用いて第1のパリティチェック結果である値を取得し、第2パリティチェック演算部は、受信データの各ビット位置においてサイクルスリップが発生した場合のパリティチェック結果である値からなるテーブルを用いて第2のパリティチェック結果である値を取得するようにしてもよい。
上記サイクルスリップ検出装置のパリティチェック演算部は、受信データのビット値d−d31を数式1に適用して第1のパリティチェック結果である値P〜Pを取得するようにしてもよい。
Figure 0005553909
また、上記サイクルスリップ検出装置は、サイクルスリップ判定部による判定結果に応じて、位置の測位を行なう測位演算部を更に備えていてもよい。
本発明の別の観点によれば、データ取得部と、パリティチェック演算部と、サイクルスリップ判定部と、を備える集積回路が提供される。データ取得部は、測位用衛星からの受信データを取得する。パリティチェック演算部は、受信データに対しパリティチェックを行なう。サイクルスリップ判定部は、パリティチェック演算部によるパリティチェック結果に基づき受信データにおけるサイクルスリップの有無を判定する。
上記集積回路は、サイクルスリップ判定部による判定結果に応じて、位置の測位を行なう測位演算部を更に備えていてもよい。
上記集積回路は更に、アンテナからの受信信号をデジタル信号に変換し、データ取得部に入力するRF部、を備えていてもよい。
本発明の別の観点によれば、上記集積回路と、測位用衛星からの信号を受信するアンテナと、を備えるGNSS受信装置が提供される。
本発明の更に別の観点によれば、上記GNSS受信装置を備える情報端末機器が提供される。
本発明の更に別の観点によれば、測位用衛星からの受信データを取得するデータ取得ステップと、受信データに対しパリティチェックを行なうパリティチェック演算ステップと、パリティチェック演算部によるパリティチェック結果に基づき受信データにおけるサイクルスリップの有無を判定するサイクルスリップ判定ステップと、備えるサイクルスリップ検出方法が提供される。
本発明の更に別の観点によれば、測位用衛星からの受信データを取得するデータ取得ステップと、受信データに対しパリティチェックを行なうパリティチェック演算ステップと、パリティチェック演算部によるパリティチェック結果に基づき受信データにおけるサイクルスリップの有無を判定するサイクルスリップ判定ステップと、をコンピュータに実行させるプログラムが提供される。
本発明によれば、簡単な構成で、精度の高いサイクルスリップの検出を実現することができる。
第1実施形態に係るGPS受信機を含むシステムの概略構成を示すブロック図である。 第1実施形態に係るGPS受信機のサイクルスリップ検出部の機能を説明するブロック図である。 航法メッセージデータの構成を示す図である。 航法メッセージデータにおけるワードの構成を示す図である。 サイクルスリップが生じた場合のデータ状態を説明する図である。 第1実施形態に係るサイクルスリップ検出部の動作を示すフローチャートである。 第2実施形態に係るGPS受信機のサイクルスリップ検出部の機能を説明するブロック図である。 第2実施形態に係るサイクルスリップ検出部の動作を示すフローチャートである。
<1 第1実施形態>
以下、本発明のGNSS受信装置の一実施形態を図面を参照して説明する。本実施形態に係るGNSS受信装置は、例えばGPS(Global Positioning System)に用いられるGPS受信機1である。
<1.1 GPS受信機1>
図1は、GPS受信機1の概略構成を示すブロック図である。
図1に示すGPS受信機1は、アンテナ10、RF部20、及びベースバンド部30を備える。GPS受信機1は情報端末機器40に接続されている。情報端末機器40とは、所定の通信や情報処理を行なう情報通信装置を備えた、例えば携帯端末や移動体端末を含むコンピュータ機器等(情報端末機器)である。ここで接続とは有線、無線のいずれによって接続されてもよく、またGPS受信機1が情報端末機器40に内蔵されている場合も含む。
アンテナ10は、GPS衛星(図示省略)から送出される高周波信号(衛星信号)を受信する。
RF部20は、アンテナ10からの高周波信号を中間周波数にダウン・コンバートしたり、アナログ信号をデジタル信号に変換して、次のベースバンド部30で処理できるような信号フォーマットに変換したりする。RF部20は、例えば、アンテナ10からの入力信号を増幅する高周波増幅部、入力信号の周波数を低下させるミキサ、バンドパス・フィルタ、A/D変換器、基準周波数源を含むPLL回路等を備える。なお、基準周波数源としては、例えば、TCXO(Temperature Compensated Crystal Oscillator)等が用いられる。この基準周波数源は、RF部20の信号生成のために用いられる他、ベースバンド部30のサンプリングクロックを生成してベースバンド部30へ出力する。
ベースバンド部30は、RF部20からデータ信号(受信データ)とクロック信号を入力し、データ信号をデコードし、GPS受信機1の位置の測位を行う。以下、ベースバンド部30の詳細な構成について説明する。
なお、本実施形態においては、RF部20及びベースバンド部30は1チップの集積回路で構成されていてもよいし、それぞれ別々のチップで構成されていてもよい。
<1.2 ベースバンド部30>
図1に示すベースバンド部30は、一例であり、相関器31と、PLL32と、CPU33と、UART34と、メモリインターフェース35と、ROM36と、RAM37と、RTC38と、を備える。
相関器31は、CPUの設定にしたがって動作する複数組の相関器から構成され、相関器31の出力にピークが出るように位相制御を行う。この位相制御によりGPS信号の捕捉や追尾(トラッキング)が行われる。
PLL32は、RF部20から入力されたクロック信号をCPU33に出力する。
CPU33は、PLL32から入力されたクロック信号に基づき動作し、GPS信号の捕捉/トラッキング、後述する航法メッセージデータのデコードや管理、測位計算、衛星の位置予測計算等を行う。
UART34は、GPS受信機1に接続された情報端末機器40との間のデータの入出力等を行う。メモリインターフェース35は、外部メモリとのインターフェースとして機能する。ROM36は、マスクROMやフラッシュROMとして構成され、所定のプログラム等を格納する。RAM37は、バッテリバックアップ領域等を有する。RTC38(Realtime Clock)は、電源投入時に参照され、この時刻を利用してGPS衛星の位置予測計算を行う。
なお、上記ベースバンド部30の構成は一例に過ぎず、その他の構成であってもよい。
<1.3 サイクルスリップ検出部300>
図2は、ベースバンド部30により実現されるサイクルスリップ検出部(サイクルスリップ検出装置)300の機能ブロック図である。サイクルスリップ検出部300は、上述したベースバンド部30のCPUがROM36等に記憶された所定のプログラムを実行することによってサイクルスリップの検出を実現するものである。
このサイクルスリップの検出は、データ信号補足/追尾部310より取得される、位相制御されているデータ信号(受信データ)に対して行われる。このデータ信号は、GPS衛星から送信される航法メッセージデータを含む。航法メッセージデータは、次のような構成を有する。
<1.3.1 航法メッセージデータの構成>
航法メッセージデータの1サイクルは、フレームという単位の集合であり、図3に示す構造を有する。1フレームは1500ビットであり、これを送信するには30秒の時間がかかる。フレームは5つのサブフレーム(各300ビット)から構成されており、サブフレーム1から順番に送信を始め、サブフレーム5を送信し終わると再びサブフレーム1の送信に戻る。
サブフレームの内部は、図3に示すようにワードという単位に分割されている。1ワードは30ビットである。1サブフレームは10ワードに対応し、各ワードは24ビットのデータ部とパリティチェック用の6ビットから構成される。
図4は、1ワード(30ビット)の構成を示す。dからd31はそれぞれ1ビットの情報に相当する。d及びdは、前ワードの終端2ビットであり当該ワードのパリティチェック計算に用いられる。dからd25は、当該ワードのデータビットであり、d26からd31は当該ワードのパリティビットである。
なお、航法メッセージデータにサイクルスリップが生じた場合、データは図5に示すような状態になる。図5(a)は追尾衛星が送信したデータ、つまりサイクルスリップが生じていない状態のデータを示し、図5(b)受信機1が収集したデータでサイクルスリップが生じているデータを示す。同図に示すように、サイクルスリップが生じた時点でビット値が反転し、その後も反転し続けている。このビット値が最初に反転したタイミングがサイクルスリップが検出されたタイミングである。
<1.3.2 サイクルスリップ検出部300の構成>
図2に戻ると、サイクルスリップ検出部300は、サイクルスリップ推定部301、メッセージ収集部(第1パリティチェック演算部)302、第2パリティチェック演算部303、及びパリティチェック比較部(サイクルスリップ判定部)304を備える。
サイクルスリップ推定部301は、データ信号補足/追尾部310より位相制御されているデータ信号を取得し、その受信データに基づき、公知の手法によってサイクルスリップの有無を予備的に判定する。サイクルスリップ推定部301は、サイクルスリップが起こったと判定した場合はその検出したタイミングを取得する。
ここでのサイクルスリップの判定手法の一例としては、例えば、次のようなものがある。GPS信号の搬送波の位相値を一定時間積算し、この位相積算値に基づいて移動体の前記一定時間の平均速度を求め、GPS衛星−移動体の相対速度による送信周波数のずれを表す搬送波追尾周波数に基づいて前記移動体の瞬時速度を求める。そして、上記平均速度と瞬時速度とを比較し、これらがほぼ一致していればサイクルスリップは生じていないと判断することができ、一方、これらの差が大きければ、位相積算値にサイクルスリップが生じていると判断することができる(例えば、特開平09−145816号公報)。或いは、GPS受信装置と慣性装置から算出された速度差を基にサイクルスリップを検出する(例えば、特許文献1)こともできる。その他にもサイクルスリップの検出手法は様々な公知の手法があり、本実施形態においてはいずれの手法も適用可能である。
サイクルスリップ推定部301は、サイクルスリップがあると判定した場合は、その旨をメッセージ収集部302に通知するとともに、そのサイクルスリップを検出したタイミングを第2パリティチェック演算部303に出力する。
メッセージ収集部302は、サイクルスリップ推定部301によるサイクルスリップ有りとの判定を受けて、データ信号補足/追尾部310より入力された、位相制御されているデータ信号に対し、所定の式を用いてパリティチェックを行い、その結果(第1のパリティチェック結果)をパリティチェック比較部304に出力する。所定の式とは、例えば次のような式である。
Figure 0005553909
上記数式1においては、当該ワード及び前ワードに含まれるビット値d〜d31を用いて、PからPのビット値を算出し、パリティチェック比較部304に出力する。ここでは、出力されたパリティチェック結果の値は、例えば、P=0、P=1、P=1、P=1、P=1、P=1であったとする。
一方、第2パリティチェック演算部303は、サイクルスリップ推定部301より入力されたサイクルスリップ検出タイミングに基づき、データ信号においてサイクルスリップが生じたビット位置を特定する。第2パリティチェック演算部303は、例えば、データ信号の各ビット位置においてサイクルスリップが発生した場合のパリティチェック結果の値を予めメモリ等に記憶しておき、その値を参照する。このサイクルスリップが発生した場合のパリティチェック結果の値は、例えば上記数式1を用いて算出され、例えば次の表1に示すようなテーブルとして予め作成されメモリに記憶されている。
なお、表1は、1番目のビット位置として、前ワードの終端ビット(図4のdに相当)を含む。このビット位置でサイクルスリップが検出されることにより、それ以降に続く当該ワード(図4のd〜d30)においてサイクルスリップが発生しているかどうかを特定でできる。一方、表1は、当該ワードの終端ビット((図4のd31に相当)を含まないが、これは、このビット位置では、当該ワードに関する限り単独エラーなのかサイクルスリップなのか判別できないためである。
Figure 0005553909
第2パリティチェック演算部303は、サイクルスリップが検出されたタイミングに相当するビット位置に対応する値P〜Pをパリティチェック結果(第2のパリティチェック結果)として取得する。ここでは、例えば、サイクルスリップが検出されたビット位置がd13であるとする。この場合、d13に対応する値(ここでは、P=0、P=1、P=1、P=1、P=1、P=1)をパリティチェック結果としてパリティチェック比較部304に出力する。
パリティチェック比較部304は、メッセージ収集部302から入力されたパリティチェック結果の値と、第2パリティチェック演算部303から入力されたパリティチェック結果の値とが一致するかどうか比較し、その結果を測位演算部305に出力する。ここで、例えば、メッセージ収集部302及び第2パリティチェック演算部303共、出力された値はP=0、P=1、P=1、P=1、P=1、P=1となるため、一致すると判断されその結果が測位演算部320に出力される。
測位演算部320は、パリティチェック比較部304より入力された比較結果に応じて、GPS受信機1の位置を測位するための演算を行う。一致するとの比較結果であった場合は、サイクルスリップが正確に検出されたことになる。したがって、測位演算部320は該当する航法メッセージデータを破棄又は修正等し、同航法メッセージデータに基づく観測量を測位に使用しないようにする。一方、一致しないとの比較結果があった場合は、航法メッセージデータを破棄又は修正する。
<1.3.3 サイクルスリップ検出部300の動作>
以下、サイクルスリップ検出部300の動作について、図6のフローチャートを参照しながら説明する。
ステップS101:サイクルスリップ推定部301は、データ信号補足/追尾部310より位相制御されている航法メッセージデータを含むデータ信号を取得する。
ステップS102:サイクルスリップ推定部301は、データ信号に基づき、公知の手法によってサイクルスリップの有無を予備的に判定する。サイクルスリップ有りと判定した場合はステップS103に進み、サイクルスリップ無しと判定した場合はステップS107に進む。
ステップS103:サイクルスリップ推定部301は、サイクルスリップの検出タイミングを第2パリティチェック演算部303に出力する。サイクルスリップ推定部301はまた、メッセージ収集部302にサイクルスリップ有りとの判定結果を出力する。
ステップS104:メッセージ収集部302の第1パリティチェック演算部は、データ信号に対し、上述のように数式1によりパリティチェック結果P〜Pを算出する。一方第2パリティチェック演算部303は、上述のように表1を参照して、サイクルスリップの検出タイミングであるビット位置に対応するP〜Pの値を取得する。第1及び第2パリティチェック演算部によるパリティチェック結果はパリティチェック比較部304に入力される。
ステップS105:パリティチェック比較部304は、第1及び第2パリティチェック演算部が取得した第1及び第2パリティチェック結果P〜Pが一致するかどうかを判定する。一致する場合は、ステップS106に進み、一致しない場合はステップS107に進む。
ステップS106:パリティチェック比較部304は、データ信号にサイクルスリップがあると最終的に判定し、その判定結果を測位演算部305に出力する。サイクルスリップが有ると判定された場合、測位演算部305は、そのデータ信号に基づく観測量を測位に使用しないようにするため、そのデータ信号を破棄又は修正等する。
ステップS107:パリティチェック比較部304は、データ信号にサイクルスリップが無いと判定し、その判定結果を測位演算部305に出力する。サイクルスリップが無いと判定された場合、測位演算部305は、そのデータ信号に基づく観測量について少なくともサイクルスリップを考慮することなく測位に使用することができる。
なお、上記サイクルスリップ検出部の動作の流れは一例であり、各処理の順序やタイミング等は上記に限定されるものではない。
<1.4 第1実施形態の特徴>
本実施形態に係るサイクルスリップ検出部300によれば、サイクルスリップ推定部301で予備的にサイクルスリップの判定を行なった後、パリティチェックを用いて最終的にサイクルスリップの判定を行なっているため、簡単且つ精度よくサイクルスリップの検出が可能である。
<2.第2実施形態>
第1実施形態のGPS受信機1のサイクルスリップ検出部300では、サイクルスリップ推定部301によりサイクルスリップの有無を予備的に判定した後、パリティチェックを用いてその判定の当否を最終的に判断しているが、本発明はこれに限定されない。
第2実施形態に係るサイクルスリップ検出部2300は、パリティチェックを用いて直接サイクルスリップの有無を検出する一例である。なお、本実施形態においては、受信データに1ビットの単独エラーはないことを前提に、サイクルスリップの判定を行うものとする。
なお、図1に示す第1実施形態の受信機1の構成は本実施形態においても同様であり、以下の説明において、第1実施形態と同様の構成については説明を省略し、同様の構成要素ついては同一の符号を付して説明する。
<2.1 サイクルスリップ検出部2300の構成>
サイクルスリップ検出部(サイクルスリップ検出装置)2300は、メッセージ収集部(パリティチェック演算部)2302、及びサイクルスリップ判定部2305を備える。
メッセージ収集部(パリティチェック演算部)2302は、第1パリティチェック演算部2302aと第2パリティチェック演算部2302bとを有する。第1パリティチェック演算部2302aは、データ信号補足/追尾部310より位相制御されている航法メッセージデータを含むデータ信号(受信データ)を取得し、所定の式を用いてパリティチェック結果(第1のパリティチェック結果)を算出する。所定の式とは第1実施形態と同じである次の数式1である。
Figure 0005553909
上記数式1においては、当該ワード及び前ワードに含まれるビットd〜d31を用いて、PからPの値を算出する。ここでは、出力された値は、例えば、P=0、P=1、P=1、P=1、P=1、P=1であったとする。
次に、メッセージ収集部2302の第2パリティチェック演算部2302bは、第1パリティチェック演算部2302aにより算出された上記パリティチェック結果の値P〜Pを取得し、同値が予め記憶された値と一致するかどうかを判定する。
ここで、予め記憶された値とは、例えば、メッセージ収集部2302がデータ信号の各ビット位置においてサイクルスリップが発生した場合のパリティチェック結果の値を予めメモリ等に記憶したものである。この各ビット位置においてサイクルスリップが発生した場合のパリティチェック結果の値は、例えば上記数式1を用いて算出され、例えば第1実施形態と同じ表1に示すようなテーブルとして予め作成されている。
ここでは、第1パリティチェック演算部2302aにより算出されたパリティチェックの値P〜P(ここでは、P=0、P=1、P=1、P=1、P=1、P=1)は、次の表1を参照すると、ビット位置d13とd27におけるパリティチェック結果の値P〜Pと一致していることが分かる。
Figure 0005553909
なお、表1におけるd13とd27の値は同じであるため、ここではサイクルスリップを検出したタイミング(ビット位置)を一意には取得できないが、当該航法メッセージデータにおいてサイクルスリップが発生していることは明らかになる。
メッセージ収集部2302は、このチェック結果をサイクルスリップ判定部2305に出力する。
サイクルスリップ判定部2305は、メッセージ収集部2302より入力されたパリティチェック結果に基づき、サイクルスリップの有無を判定する。具体的には、メッセージ収集部2302の第1パリティチェック演算部2302aによりによって算出されたパリティチェック結果の値P〜Pが、第2パリティチェック演算部2302bの上記表1のいずれかの値P〜Pと一致した場合、サイクルスリップ有りと判定し、一致しない場合はサイクルスリップ無しと判定し、その判定結果を測位演算部320に出力する。
測位演算部320は、パリティチェック判定部2305より入力された判定結果に応じて、GPS受信機1の位置を測位するための演算を行う。例えば、サイクルスリップが有りと判定された場合は該当する航法メッセージデータを破棄又は修正等し、サイクルスリップ無しと判定された場合は、サイクルスリップが検出されなかったものとして、当該航法メッセージデータを扱う。
<2.2 サイクルスリップ検出部2300の動作>
以下、第2実施形態に係るサイクルスリップ検出部2300の動作について、図8のフローチャートを参照しながら説明する。
ステップS201:メッセージ収集部2302は、データ信号補足/追尾部310より位相制御されている航法メッセージデータを含むデータ信号を取得する。
ステップS202:メッセージ収集部2302の第1パリティチェック演算部2302aは、取得したデータ信号に対し、上述の通り数式1によりパリティチェック結果P〜Pを算出する。
ステップS203:メッセージ収集部2302の第2パリティチェック演算部2302bは、上述のように表1を参照して、第1パリティチェック演算部2302aにより算出したP〜Pと一致する値が表1にあるかどうかを判定する。一致する値がある場合は、ステップS204に進み、一致する値がない場合はステップS205に進む。
ステップS204:メッセージ収集部2302より、一致する値があるとのチェック結果を受けたサイクル判定部2305は、サイクルスリップ有りと判定する。この判定を受けた測位演算部305は、当該データ信号を破棄又は修正等して、そのデータ信号に基づく観測量を測位に使用しないようにする。
ステップS205:メッセージ収集部2302より、一致する値が無いとのチェック結果を受けたサイクル判定部2305は、サイクルスリップ無しと判定する。この判定を受けた測位演算部305は、当該データ信号に基づく観測量について少なくともサイクルスリップを考慮することなく測位に使用することができる。
なお、上記サイクルスリップ検出部の動作の流れは一例であり、各処理の順序やタイミング等は上記に限定されるべきでない。
<2.3 第2実施形態の特徴>
本実施形態に係るサイクルスリップ検出部2300によれば、パリティチェックを利用してサイクルスリップを検出することができるため、従来に比して簡単で、正確にサイクルスリップの検出が可能である。
また、追尾する一チャンネルの観測量から、閾値等を用いることなく、サイクルスリップの検出を行なうことが可能となる点においても、従来に比べて簡単にサイクルスリップの検出が可能である。
<2.4 第2実施形態の変形例>
上記第2実施形態では、メッセージ収集部2302の第1パリティチェック演算部2302aがパリティチェック結果P〜Pを算出した後、第2パリティチェック演算部2302bにより表1を参照しているがこれに限定されない。例えば、第1パリティチェック演算部2302aによるパリティチェック結果P〜Pが全て0の値を示していれば、サイクルスリップは生じていないと判断し、表1等の参照はせずに、その結果をサイクルスリップ判定部2305に入力し、サイクルスリップは無しと判定するようにしてもよい。
<3 その他実施形態>
<3.1>
上記第1及び第2実施形態においては、サイクルスリップが生じたことを示す所定の値P〜Pは表1として格納され、参照されているが、これに限定されない。かかるサイクルスリップが生じたことを示す所定の値はリアルタイムに算出するようにしてもよい。
<3.2>
上記第1及び第2実施形態に係るGPS受信機1は例示に過ぎず、これに基づいて以下に示すような変更を行うことが可能である。
上記各実施形態では、本発明のGNSS受信機の一例としてGPS受信機1を適用した場合について説明したが、受信機の他の例として、測位システムからの衛星信号を利用可能な受信機であればよく、例えば、ガリレオ(Galileo Positioning System)に適用されるGNSS受信機などでもよい。
<3.3>
上記第1及び第2実施形態では、アンテナ10、RF部20およびベースバンド部30を備えるGPS受信機1について示しているが、本発明は、例えばベースバンド部30のみを備える形態(例えば、集積回路、半導体、モジュール、チップ)、または、RF部20およびベースバンド部30のみを備える形態(例えば、集積回路、半導体、チップ、モジュール)と適用することも可能である。
<3.4>
上記第1及び第2実施形態においては、サイクルスリップ検出部300は装置や集積回路として説明したが、本発明は、サイクルスリップ検出方法やコンピュータプログラムとしても実現し得るものである。
1 GPS受信機(GNSS受信装置)
10 アンテナ
20 RF部
30 ベースバンド部
31 相関器
32 PLL
33 CPU
34 UART
35 メモリインターフェース
36 ROM
37 RAM
38 RTC
40 情報端末機器
300 サイクルスリップ検出部(サイクルスリップ検出装置)
301 サイクルスリップ推定部
302 メッセージ収集部(第1パリティチェック演算部)
303 第2パリティチェック演算部
304 パリティチェック比較部(サイクルスリップ判定部)
320 測位演算部
2300 サイクルスリップ検出部(サイクルスリップ検出装置)
2302 メッセージ収集部(パリティチェック演算部)
2305 サイクルスリップ判定部

Claims (12)

  1. 測位用衛星からの受信データを取得するデータ取得部と、
    前記受信データに対しパリティチェックを行なうパリティチェック演算部と、
    前記パリティテックの結果に基づき前記受信データにおけるサイクルスリップの有無を判定するサイクルスリップ判定部と、を備えるサイクルスリップ検出装置において、
    前記パリティチェック演算部は、
    前記パリティチェックにより第1のパリティチェック結果である値を取得する第1パリティチェック演算部と、
    前記第1のパリティチェック結果の値と、前記受信データの各ビット位置においてサイクルスリップが発生した場合のパリティチェック結果である第2のパリティチェック結果の値とを比較する第2パリティチェック演算部と、を有し、
    前記サイクルスリップ判定部は、
    前記第2パリティチェック演算部による前記比較の結果、前記第1のパリティチェック結果の値が前記各ビット位置で取得した前記第2パリティチェック結果の値のいずれかと一致すると判断した場合、前記受信データにおいてサイクルスリップが有ると判定する、ことを特徴とするサイクルスリップ検出装置。
  2. 測位用衛星からの受信データを取得するデータ取得部と、
    前記受信データに対しパリティチェックを行なうパリティチェック演算部と、
    前記パリティテックの結果に基づき前記受信データにおけるサイクルスリップの有無を判定するサイクルスリップ判定部と、
    前記データ取得部により取得された前記受信データにおけるサイクルスリップの検出タイミングを推定するサイクルスリップ推定部と、
    を備えるサイクルスリップ検出装置において、
    前記パリティチェック演算部は、
    前記パリティチェックにより第1のパリティチェック結果である値を取得する第1パリティチェック演算部と、
    前記サイクルスリップの検出タイミングに応じた第2のパリティチェック結果である値を取得する第2パリティチェック演算部と、を有し、
    前記サイクルスリップ判定部は、
    前記第1のパリティチェック結果である値と前記第2のパリティチェック結果である値とを比較し、一致する場合は前記受信データにおいてサイクルスリップが有ると判定する、ことを特徴とするサイクルスリップ検出装置。
  3. 前記第1パリティチェック演算部は、所定の式を用いて前記第1のパリティチェック結果である値を取得し、
    前記第2パリティチェック演算部は、前記受信データの各ビット位置においてサイクルスリップが発生した場合のパリティチェック結果である値からなるテーブルを用いて前記第2のパリティチェック結果である値を取得する、請求項1または2に記載のサイクルスリップ検出装置。
  4. 前記パリティチェック演算部は、前記受信データのビット値d−d31を数式1に適用して前記第1のパリティチェック結果である値P〜Pを取得する、請求項1から3のいずれかに記載のサイクルスリップ検出装置。
    Figure 0005553909
  5. 測位用衛星からの受信データを取得するデータ取得部と、
    前記受信データに対しパリティチェックを行なうパリティチェック演算部と、
    前記パリティテックの結果に基づき前記受信データにおけるサイクルスリップの有無を判定するサイクルスリップ判定部と、を備える集積回路において、
    前記パリティチェック演算部は、
    前記パリティチェックにより第1のパリティチェック結果である値を取得する第1パリティチェック演算部と、
    前記第1のパリティチェック結果の値と、前記受信データの各ビット位置においてサイクルスリップが発生した場合のパリティチェック結果である第2のパリティチェック結果の値とを比較する第2パリティチェック演算部と、を有し、
    前記サイクルスリップ判定部は、
    前記第2パリティチェック演算部による前記比較の結果、前記第1のパリティチェック結果の値が前記各ビット位置で取得した前記第2パリティチェック結果の値のいずれかと一致すると判断した場合、前記受信データにおいてサイクルスリップが有ると判定する、ことを特徴とする集積回路。
  6. 測位用衛星からの受信データを取得するデータ取得部と、
    前記受信データに対しパリティチェックを行なうパリティチェック演算部と、
    前記パリティテックの結果に基づき前記受信データにおけるサイクルスリップの有無を判定するサイクルスリップ判定部と、
    前記データ取得部により取得された前記受信データにおけるサイクルスリップの検出タイミングを推定するサイクルスリップ推定部と、
    を備える集積回路において、
    前記パリティチェック演算部は、
    前記パリティチェックにより第1のパリティチェック結果である値を取得する第1パリティチェック演算部と、
    前記サイクルスリップの検出タイミングに応じた第2のパリティチェック結果である値を取得する第2パリティチェック演算部と、を有し、
    前記サイクルスリップ判定部は、
    前記第1のパリティチェック結果である値と前記第2のパリティチェック結果である値とを比較し、一致する場合は前記受信データにおいてサイクルスリップが有ると判定する、ことを特徴とする集積回路。
  7. 請求項5または6のいずれかに記載の集積回路と、
    前記測位用衛星からの信号を受信するアンテナと、を備える、GNSS受信装置。
  8. 請求項7に記載のGNSS受信装置、を備える、情報端末機器。
  9. 測位用衛星からの受信データを取得するデータ取得ステップと、
    前記受信データに対しパリティチェックを行なうパリティチェック演算ステップと、
    前記パリティテックの結果に基づき前記受信データにおけるサイクルスリップの有無を判定するサイクルスリップ判定ステップと、からなるサイクルスリップ検出方法において、
    前記パリティチェック演算ステップは、
    前記パリティチェックにより第1のパリティチェック結果である値を取得する第1パリティチェック演算ステップと、
    前記第1のパリティチェック結果の値と、前記受信データの各ビット位置においてサイクルスリップが発生した場合のパリティチェック結果である第2のパリティチェック結果の値とを比較する第2パリティチェック演算ステップと、を有し、
    前記サイクルスリップ判定ステップは、
    前記第2パリティチェック演算ステップによる前記比較の結果、前記第1のパリティチェック結果の値が前記各ビット位置で取得した前記第2パリティチェック結果の値のいずれかと一致すると判断した場合、前記受信データにおいてサイクルスリップが有ると判定する、ことを特徴とするサイクルスリップ検出方法。
  10. 測位用衛星からの受信データを取得するデータ取得ステップと、
    前記受信データに対しパリティチェックを行なうパリティチェック演算ステップと、
    前記パリティテックの結果に基づき受信データにおけるサイクルスリップの有無を判定するサイクルスリップ判定部と、
    前記データ取得部により取得された前記受信データにおけるサイクルスリップの検出タイミングを推定するサイクルスリップ推定ステップと、からなるサイクルスリップ検出方法において、
    前記パリティチェック演算ステップは、
    前記パリティチェックにより第1のパリティチェック結果である値を取得する第1パリティチェック演算ステップと、
    前記サイクルスリップの検出タイミングに応じた第2のパリティチェック結果である値を取得する第2パリティチェック演算ステップと、を有し、
    前記サイクルスリップ判定ステップは、
    前記第1のパリティチェック結果である値と前記第2のパリティチェック結果である値とを比較し、一致する場合は前記受信データにおいてサイクルスリップが有ると判定する、ことを特徴とするサイクルスリップ検出方法。
  11. 測位用衛星からの受信データを取得するデータ取得ステップと、
    前記受信データに対しパリティチェックを行なうパリティチェック演算ステップと、
    前記パリティテックの結果に基づき前記受信データにおけるサイクルスリップの有無を判定するサイクルスリップ判定ステップと、をコンピュータに実行させるプログラムにおいて、
    前記パリティチェック演算ステップは、
    前記パリティチェックにより第1のパリティチェック結果である値を取得する第1パリティチェック演算ステップと、
    前記第1のパリティチェック結果の値と、前記受信データの各ビット位置においてサイクルスリップが発生した場合のパリティチェック結果である第2のパリティチェック結果の値とを比較する第2パリティチェック演算ステップと、を有し、
    前記サイクルスリップ判定ステップは、
    前記第2パリティチェック演算ステップによる前記比較の結果、前記第1のパリティチェック結果の値が前記各ビット位置で取得した前記第2パリティチェック結果の値のいずれかと一致すると判断した場合、前記受信データにおいてサイクルスリップが有ると判定する、ことを特徴とするプログラム。
  12. 測位用衛星からの受信データを取得するデータ取得ステップと、
    前記受信データに対しパリティチェックを行なうパリティチェック演算ステップと、
    前記パリティテックの結果に基づき受信データにおけるサイクルスリップの有無を判定するサイクルスリップ判定部と、
    前記データ取得部により取得された前記受信データにおけるサイクルスリップの検出タイミングを推定するサイクルスリップ推定ステップと、をコンピュータに実行させるプログラムにおいて、
    前記パリティチェック演算ステップは、
    前記パリティチェックにより第1のパリティチェック結果である値を取得する第1パリティチェック演算ステップと、
    前記サイクルスリップの検出タイミングに応じた第2のパリティチェック結果である値を取得する第2パリティチェック演算ステップと、を有し、
    前記サイクルスリップ判定ステップは、
    前記第1のパリティチェック結果である値と前記第2のパリティチェック結果である値とを比較し、一致する場合は前記受信データにおいてサイクルスリップが有ると判定する、ことを特徴とするプログラム。
JP2012548686A 2010-12-13 2011-07-14 サイクルスリップ検出装置、サイクルスリップ検出のための集積回路、同集積回路を備えるgnss受信装置、同gnss受信装置を備える情報端末機器、サイクルスリップ検出方法及びプログラム Expired - Fee Related JP5553909B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2012548686A JP5553909B2 (ja) 2010-12-13 2011-07-14 サイクルスリップ検出装置、サイクルスリップ検出のための集積回路、同集積回路を備えるgnss受信装置、同gnss受信装置を備える情報端末機器、サイクルスリップ検出方法及びプログラム

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2010276976 2010-12-13
JP2010276976 2010-12-13
JP2012548686A JP5553909B2 (ja) 2010-12-13 2011-07-14 サイクルスリップ検出装置、サイクルスリップ検出のための集積回路、同集積回路を備えるgnss受信装置、同gnss受信装置を備える情報端末機器、サイクルスリップ検出方法及びプログラム
PCT/JP2011/066132 WO2012081275A1 (ja) 2010-12-13 2011-07-14 サイクルスリップ検出装置、サイクルスリップ検出のための集積回路、同集積回路を備えるgnss受信装置、同gnss受信装置を備える情報端末機器、サイクルスリップ検出方法及びプログラム

Publications (2)

Publication Number Publication Date
JPWO2012081275A1 JPWO2012081275A1 (ja) 2014-05-22
JP5553909B2 true JP5553909B2 (ja) 2014-07-23

Family

ID=46244391

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012548686A Expired - Fee Related JP5553909B2 (ja) 2010-12-13 2011-07-14 サイクルスリップ検出装置、サイクルスリップ検出のための集積回路、同集積回路を備えるgnss受信装置、同gnss受信装置を備える情報端末機器、サイクルスリップ検出方法及びプログラム

Country Status (2)

Country Link
JP (1) JP5553909B2 (ja)
WO (1) WO2012081275A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111190200A (zh) * 2019-12-09 2020-05-22 北京时代民芯科技有限公司 一种动态环境下单频周跳检测与修复方法

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140062778A1 (en) * 2012-08-31 2014-03-06 Javad Gnss, Inc. Cycle slip detection
US9231690B2 (en) * 2013-09-17 2016-01-05 Qualcomm Incorporated Signal tracking and decoding in GNSS
CN106199659B (zh) * 2016-07-04 2018-08-21 武汉大学 基于模糊数学的gnss单站双频观测数据周跳探测及处理方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11237463A (ja) * 1998-02-20 1999-08-31 Sokkia Co Ltd Gps受信機の同期判定装置
JP2002530627A (ja) * 1998-05-06 2002-09-17 スナップトラック・インコーポレーテッド 衛星ポジショニングシステムにおける信号処理方法及び装置
JP2011052987A (ja) * 2009-08-31 2011-03-17 Casio Computer Co Ltd 測位データ受信機、エラー訂正装置およびプログラム
JP2011069629A (ja) * 2009-09-24 2011-04-07 Seiko Epson Corp 擬似衛星信号受信方法及び受信装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11237463A (ja) * 1998-02-20 1999-08-31 Sokkia Co Ltd Gps受信機の同期判定装置
JP2002530627A (ja) * 1998-05-06 2002-09-17 スナップトラック・インコーポレーテッド 衛星ポジショニングシステムにおける信号処理方法及び装置
JP2011052987A (ja) * 2009-08-31 2011-03-17 Casio Computer Co Ltd 測位データ受信機、エラー訂正装置およびプログラム
JP2011069629A (ja) * 2009-09-24 2011-04-07 Seiko Epson Corp 擬似衛星信号受信方法及び受信装置

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111190200A (zh) * 2019-12-09 2020-05-22 北京时代民芯科技有限公司 一种动态环境下单频周跳检测与修复方法

Also Published As

Publication number Publication date
WO2012081275A1 (ja) 2012-06-21
JPWO2012081275A1 (ja) 2014-05-22

Similar Documents

Publication Publication Date Title
US9739888B2 (en) Method, system and device for position determination with predicted ephemeris
WO2013016800A1 (en) System, method, and computer program for a low power and low cost gnss receiver
US11709278B2 (en) Method and system for signal detection including positioning signals
JP5553909B2 (ja) サイクルスリップ検出装置、サイクルスリップ検出のための集積回路、同集積回路を備えるgnss受信装置、同gnss受信装置を備える情報端末機器、サイクルスリップ検出方法及びプログラム
JP5968998B2 (ja) 測位信号検出方法、測位信号検出プログラム、測位信号受信装置、測位装置および情報機器端末
US9287921B2 (en) Method and apparatus for frame synchronization in a positioning system
US20120022779A1 (en) Navigation receivers
JP2016532119A (ja) クラウドオフロード型全地球衛星測位
JP2008170338A (ja) 衛星航法装置、外部基地局、及び衛星測位システム
CN112327335A (zh) Gnss接收机和卫星捕获跟踪方法
US20120027133A1 (en) Demodulation of data collected prior to bit edge detection
JP4888110B2 (ja) 相関演算制御回路及び相関演算制御方法
JPWO2020100455A1 (ja) Gnss受信装置
US9184786B2 (en) Systems and methods for clock calibration for satellite navigation
US20120154217A1 (en) Method and program of acquiring navigation message, gnss receiving apparatus, and mobile terminal
JP2015001526A (ja) 測位ビット同期を実施し検証するための方法
CN108732598B (zh) 一种gnss接收机及其时间确定方法
JPH11223669A (ja) Gps受信装置及びそのzカウント抽出方法
US9014304B2 (en) Demodulation method, demodulation device, and electronic apparatus
US11841444B2 (en) Resilient ephemeris decoding of GNSS satellite information
US10578745B2 (en) Positioning signal receiving method and positioning signal receiving device
JP2009236517A (ja) Gpsレシーバ
FI111482B (fi) Menetelmä sijainninmäärityksen suorittamiseksi ja elektroniikkalaite
US9091756B2 (en) Receiving method and receiving apparatus
US20230168384A1 (en) Methods and apparatus for monitoring the quality of correcton data associated with a navigation satellite or the propagation of signals transmitted thereby

Legal Events

Date Code Title Description
A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140314

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140527

R150 Certificate of patent or registration of utility model

Ref document number: 5553909

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

LAPS Cancellation because of no payment of annual fees