(57)【要約】  【課題】ＧＰＳ信号の搬送波位相のサイクルスリップを正確かつ経済的に検知し、迅速に補正量を算出して搬送波位相データを補正することにより、高精度の測位・航法を継続的に可能とするＧＰＳサイクルスリップ自動補正手段、方法及びプログラムを提供する。  【解決手段】ＧＰＳ受信部１０で受信したＧＰＳ信号の搬送波位相データに対して外部接続したＧＰＳサイクルスリップ自動補正装置２０の、ＭＰＵ２１で実行する自動補正プログラムにて１次スプラインフィルタ処理を施し、更新前予測位相データと観測位相データとの比較からサイクルスリップ現象の発生を検知すると共に補正量を算出し、該補正量により搬送波位相データを自動補正する。なお、ＧＰＳ受信機にＧＰＳサイクルスリップ自動補正装置２０そのものを内蔵するように構成しても良い。【選択図】  図１
2006275860.tif 000002
【特許請求の範囲】
【請求項１】
受信したＧＰＳ信号の搬波位相のサイクルスリップを検出して、搬送波位相を補正するＧＰＳサイクルスリップ自動補正装置において、
前記受信したＧＰＳ信号に対して1次スプラインフィルタリング処理を施すことによって前記サイクルスリップを検出するとともに搬送波位相を補正することを特徴とするＧＰＳサイクルスリップ自動補正装置。
【請求項２】
前記1次スプラインフィルタリング処理として、観測した搬送波位相データと、現在までの受信状況に基づいて予測した予測位相データとの比較結果が予め設定した閾値に収まっているか否かに基づいて、サイクルスリップの発生の有無を検出すると共に、施すべく補正量を算出して、算出した補正量により、観測した搬送波位相データを補正することを特徴とするＧＰＳサイクルスリップ自動補正装置。
【請求項３】
ＧＰＳ信号を受信するＧＰＳ受信機において、請求項１または２に記載のＧＰＳサイクルスリップ自動補正装置を外部接続する接続手段を備えて、該ＧＰＳサイクルスリップ自動補正装置により、受信したＧＰＳ信号の搬送波位相のサイクルスリップの検出と補正を行うことを特徴とするＧＰＳ受信機。
【請求項４】
受信したＧＰＳ信号の搬送波位相のサイクルスリップを検出して、搬送波位相を補正するために、前記受信したＧＰＳ信号に対して1次スプラインフィルタリング処理を施す自動補正手段を内蔵していることを特徴とするＧＰＳ受信機。
【請求項５】
前記1次スプラインフィルタリング処理として、観測した搬送波位相データと、現在までの受信状況に基づいて予測した予測位相データとの比較結果が予め設定した閾値に収まっているか否かに基づいて、サイクルスリップの発生の有無を検出すると共に、施すべく補正量を算出して、算出した補正量により、観測した搬送波位相データを補正することを特徴とするＧＰＳ受信機。
【請求項６】
受信したＧＰＳ信号の搬送波位相のサイクルスリップを検出して、搬送波位相を補正するＧＰＳサイクルスリップ自動補正方法において、
前記受信したＧＰＳ信号に対して1次スプラインフィルタリング処理を施すことによって前記サイクルスリップを検出するとともに搬送波位相を補正することを特徴とするＧＰＳサイクルスリップ自動補正方法。
【請求項７】
前記1次スプラインフィルタリング処理として、観測した搬送波位相データと、現在までの受信状況に基づいて予測した予測位相データとの比較結果が予め設定した閾値に収まっているか否かに基づいて、サイクルスリップの発生の有無を検出すると共に、施すべく補正量を算出して、算出した補正量により、観測した搬送波位相データを補正することを特徴とするＧＰＳサイクルスリップ自動補正方法。
【請求項８】
請求項６または7に記載のＧＰＳサイクルスリップ自動補正方法をコンピュータにより実行することが可能なプログラムとして実現することを特徴とするＧＰＳサイクルスリップ自動補正プログラム。
【請求項９】
請求項８に記載のＧＰＳサイクルスリップ自動補正プログラムをコンピュータにより読み取り可能な記録媒体に記録することを特徴とするプログラム記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ＧＰＳサイクルスリップ自動補正装置、ＧＰＳ受信機、ＧＰＳサイクルスリップ自動補正方法、ＧＰＳサイクルスリップ自動補正プログラム及びプログラム記録媒体に関し、特に、ＧＰＳＲ(ＧＰＳ受信機: Global Positioning System Receiver)による搬送波位相データ計測におけるサイクルスリップ現象の発生の検出とこれに対する補正を可能とするＧＰＳサイクルスリップ自動補正装置に関する。本発明は、ＧＰＳＲによる測位・航法の分野において好適に適用可能であり、搬送波位相データ計測によりサイクルスリップの補正を経済的かつ正確に行うことができる。
【背景技術】
【０００２】
従来のＧＰＳ受信機（ＧＰＳＲ）では、搬送波位相の計測時に、障害物などで電波が遮断されて生じる、ＧＰＳ受信機の瞬時的な位相ロック外れに起因して、搬送波位相データが波長の整数倍(位相表現では、３６０degの整数倍）飛躍する現象が発生する。これを、サイクルスリップ現象と呼んでいる。
【０００３】
一般に、ＧＰＳＲが計測する他の測距データである擬似距離データのノイズレベルが数センチから数十センチであるのに比べて、搬送波位相データのノイズは３ミリから数ミリ程度であることから、ＧＰＳＲが出力する搬送波位相データの測位・航法への利用が、測位・航法の高精度化を行う上で重要な鍵となっている。すでに、ＲＴＫ(Real Time Kinematic)航法としても実際に利用されているが、この場合に、前述のようなサイクルスリップ現象が発生すると、搬送波位相による測位・航法を一旦リセットする必要が生じる。
【０００４】
この結果、サイクルスリップの発生によるリセットに伴い、データ処理に不連続を生じてしまうので、ＧＰＳＲは、搬送波位相のアンビギュイティ（Ambiguity：曖昧さ、つまり整数値バイアス）を最初から解き直さねばならないが、これには、相応の時間を掛けた擬似距離データの蓄積が必要となる。
【０００５】
このようなサイクルスリップの検出と補正を高速に行う従来技術として、例えば、特許文献１に示す特開２０００−６５５９３号公報「ハイブリッド形航法システムに記載されたものがある。該特許文献1に記載の技術は、ＧＰＳ受信機と慣性装置とをハイブリッドにして、それぞれで計測された移動速度と位置とを用いて、より正確な移動速度を算出して、サイクルスリップの検出とその補正を正確に行おうとするものである。
【特許文献１】特開２０００−６５５９３号公報（第３−４頁、図１）
【発明の開示】
【発明が解決しようとする課題】
【０００６】
しかしながら、従来の技術では、ＧＰＳ受信機能のみを活用して、サイクルスリップの発生を正確に検知することが困難であり、例えば、前記特許文献１に記載の技術でも、ＧＰＳ受信機のみならず慣性装置を併用する構成とする必要があり、移動機器の小型化が困難であり、かつ、コスト低減も困難であり、更には、取得した観測データの補正量の算出に要する計算量が多く、補正を行うまでにかなりの時間を要する。
【０００７】
つまり、従来技術における第１の課題は、観測データの補正の有無を認知するために必須である、サイクルスリップの発生を、経済的に、誤り無く正確に検知することが不可能な点にある。
【０００８】
第２の課題は、サイクルスリップを正確に検知したとしても、補正量を迅速に正しく計算して、観測データの補正を行い、データの連続性を確保することにより、補正された搬送波位相データを用いて、高精度の測位・航法を継続的に可能にするというＧＰＳ測位システムの本質的な課題を解決することができないという点にある。
【０００９】
（本発明の目的）
従って、本発明の目的は、サイクルスリップの発生を経済的かつ精度良く検知し、同時に、そのサイクルスリップ量を迅速に計算して、観測したＧＰＳ搬送波の位相データを補正することにより、ＧＰＳ受信機の搬送波位相計測に関する信頼性の向上と性能向上とを図ることが可能なＧＰＳサイクルスリップ自動補正装置、ＧＰＳ受信機、ＧＰＳサイクルスリップ自動補正方法、ＧＰＳサイクルスリップ自動補正プログラム及びプログラム記録媒体を提供することにある。
【課題を解決するための手段】
【００１０】
前述の課題を解決するため、本発明によるＧＰＳサイクルスリップ自動補正装置、ＧＰＳ受信機、ＧＰＳサイクルスリップ自動補正方法、ＧＰＳサイクルスリップ自動補正プログラム及びプログラム記録媒体は、次のような特徴的な構成を採用している。
【００１１】
（１）受信したＧＰＳ信号の搬送波位相のサイクルスリップを検出して、搬送波位相を補正するＧＰＳサイクルスリップ自動補正装置において、
前記受信したＧＰＳ信号に対して1次スプラインフィルタリング処理を施すことによって前記サイクルスリップを検出するとともに搬送波位相を補正するＧＰＳサイクルスリップ自動補正装置。
（２）前記1次スプラインフィルタリング処理として、観測した搬送波位相データと、現在までの受信状況に基づいて予測した予測位相データとの比較結果が予め設定した閾値に収まっているか否かに基づいて、サイクルスリップの発生の有無を検出すると共に、施すべく補正量を算出して、算出した補正量により、観測した搬送波位相データを補正するＧＰＳサイクルスリップ自動補正装置。
（３）ＧＰＳ信号を受信するＧＰＳ受信機において、請求項１または２に記載のＧＰＳサイクルスリップ自動補正装置を外部接続する接続手段を備えて、該ＧＰＳサイクルスリップ自動補正装置により、受信したＧＰＳ信号の搬送波位相のサイクルスリップの検出と補正を行うＧＰＳ受信機。
（４）受信したＧＰＳ信号の搬送波位相のサイクルスリップを検出して、搬送波位相を補正するために、前記受信したＧＰＳ信号に対して1次スプラインフィルタリング処理を施す自動補正手段を内蔵しているＧＰＳ受信機。
（５）前記1次スプラインフィルタリング処理として、観測した搬送波位相データと、現在までの受信状況に基づいて予測した予測位相データとの比較結果が予め設定した閾値に収まっているか否かに基づいて、サイクルスリップの発生の有無を検出すると共に、施すべく補正量を算出して、算出した補正量により、観測した搬送波位相データを補正するＧＰＳ受信機。
（６）受信したＧＰＳ信号の搬送波位相のサイクルスリップを検出して、搬送波位相を補正するＧＰＳサイクルスリップ自動補正方法において、
前記受信したＧＰＳ信号に対して1次スプラインフィルタリング処理を施すことによって前記サイクルスリップを検出するとともに搬送波位相を補正するＧＰＳサイクルスリップ自動補正方法。
（７）前記1次スプラインフィルタリング処理として、観測した搬送波位相データと、現在までの受信状況に基づいて予測した予測位相データとの比較結果が予め設定した閾値に収まっているか否かに基づいて、サイクルスリップの発生の有無を検出すると共に、施すべく補正量を算出して、算出した補正量により、観測した搬送波位相データを補正するＧＰＳサイクルスリップ自動補正方法。
（８）上記（６）または（7）のＧＰＳサイクルスリップ自動補正方法をコンピュータにより実行することが可能なプログラムとして実現するＧＰＳサイクルスリップ自動補正プログラム。
（９）上記（８）のＧＰＳサイクルスリップ自動補正プログラムをコンピュータにより読み取り可能な記録媒体に記録するプログラム記録媒体。
【発明の効果】
【００１２】
本発明によるＧＰＳサイクルスリップ自動補正装置、ＧＰＳ受信機、ＧＰＳサイクルスリップ自動補正方法、ＧＰＳサイクルスリップ自動補正プログラム及びプログラム記録媒体によれば、以下のような効果を得ることができる。 
【００１３】
すなわち、本発明によれば、サイクルスリップの検知と観測データに対する補正に１次スプラインフィルタを用いることにより、小型化が可能で、かつ、経済的な装置構成が可能になると共に、取得データを瞬時に正確な搬送波位相データに補正することが可能になり、ＧＰＳ受信機による測位・航法の精度および信頼性を大幅に向上させることができる。その理由は、正確な補正を施した結果として、取得した搬送波位相データを継続的に使用することが可能になるからである。
【００１４】
また、サイクルスリップ自動補正用の自動補正プログラムを提供することにより、該ＧＰＳサイクルスリップ自動補正プログラムをＭＰＵ（Ｍicro Processing Unit）に組み込んで実行することにより、経済的かつ融通性に富むＧＰＳサイクルスリップ自動補正装置あるいはＧＰＳ受信機を提供することも可能であり、カーナビや携帯電話機などの小型の移動機器にも好適に適用することができる。
【発明を実施するための最良の形態】
【００１５】
以下、本発明によるＧＰＳサイクルスリップ自動補正装置、ＧＰＳ受信機、ＧＰＳサイクルスリップ自動補正方法、ＧＰＳサイクルスリップ自動補正プログラム及びプログラム記録媒体の好適な実施例について添付図を参照して説明する。なお、以下の実施例の説明においては、本発明によるＧＰＳ受信機の具体的な構成例やＧＰＳサイクルスリップ自動補正方法及び該ＧＰＳサイクルスリップ自動補正方法をコンピュータにより実行可能なＧＰＳサイクルスリップ自動補正プログラムについて、図面を参照しながら詳細に説明する。斯かる説明により、本発明によるＧＰＳサイクルスリップ自動補正装置の具体的な実施例や、さらには、該ＧＰＳサイクルスリップ自動補正プログラムをコンピュータにより読み取り可能な記録媒体に記録したプログラム記録媒体の実施例についても、容易に想到することができるので、これらの発明に関する具体的な実施例についての説明は割愛する。
【００１６】
（本発明の最良の形態）
搬送波位相データのサイクルスリップ現象の発生を瞬時に検出して、短時間で、測位システムに用いるＧＰＳ衛星／ＧＰＳ受信機の距離の変化を直線的な変化として十分な精度で近似できる技術として、ＧＰＳ受信機の搬送波位相計測回路が出力する観測データ（搬送波位相データ）に対して１次スプラインフィルタを適用することにより、更新前の予測データと取得した観測データとの比較から、サイクルスリップ現象の発生を瞬時に検知すると共に、サイクルスリップによるデータの飛躍量が量子化されていることを利用して、補正量を瞬時に算出して、該補正量を観測した実際の搬送波位相データに加える補間処理を施すことにより、正確に補正した搬送波位相データを出力するように構成する。
【００１７】
ここで、本実施例におけるＧＰＳ受信機（ＧＰＳＲ）とは、このようなアルゴリズムを実現する自動補正機能例えばＧＰＳサイクルスリップ自動補正プログラムを搭載したＭＰＵ(Micro Processing Unit)部分を装備したＧＰＳ受信機である。
【００１８】
もしくは、本実施例におけるＧＰＳ受信機（ＧＰＳＲ）とは、このようなアルゴリズムを実現するＧＰＳサイクルスリップ自動補正プログラムを搭載したＭＰＵ(Micro Processing Unit)部分を有するＧＰＳサイクルスリップ自動補正装置を既存のＧＰＳＲに外部接続した形態からなるＧＰＳ受信機である。あるいは、該ＧＰＳサイクルスリップ自動補正装置を内蔵したＧＰＳ受信機として構成しても良い。
【００１９】
また、本実施例のアルゴリズムとして適用する1次スプラインフィルタ（linear_spline
filter）とは、有限要素法における定義域全体を最も単純に１つの区分多項式により近似することにより、取得した搬送波位相データを線形補間して出力するフィルタであり、計算処理負荷が軽く高速化を図ることができる。一般に、スプラインフィルタは、ガウシアンフィルタよりも鋭い遮断特性を示し、データ始点あるいはデータ終点領域における変形の発生を抑えることができる他、測定データに含まれる周期の長いうねり成分に対する追随性も良いことが知られており、より正確にサイクルスリップの発生を検出して補正を行うことができる。
【００２０】
本実施例におけるＧＰＳＲでは、受信した搬送波位相データを読み込んで、これに対して、前述の１次スプラインフィルタリングアルゴリズムを適用したＧＰＳサイクルスリップ自動補正プログラムを動作させることにより、搬送波位相データの補正計算を実施し、適切に補正した搬送波位相データを継続的に出力することを可能にし、もって、ＧＰＳ受信機による測位・航法の精度および信頼性が向上させようとするものである。
【００２１】
すなわち、本実施例におけるＧＰＳＲでは、受信した搬送波位相データを読み込み、該搬送波位相データの取得前の挙動に基づいて１次スプラインフィルタが推定した予測値と、実際の搬送波位相データとの比較を行い、予め設定した閾値を超えているか否かにより、両者の差異が観測ノイズレベル及びフィルタの推定精度に比べて有意であるか否かを判定し、差異が有意であると判定されたものについては、サイクルスリップが発生したものと判定し、搬送波位相データの飛躍量に関する量子化条件から、サイクルスリップによる搬送波位相データの飛躍量を決定し、該飛躍量を用いて、搬送波位相データのアンビギュイティの値を補正して更新する。
【００２２】
（構成の説明）
図１は、本発明によるＧＰＳ受信機の構成の一実施例を示す構成概念図である。図１に例示するように、本実施例によるＧＰＳ受信機１は、既存のＧＰＳ受信機と全く同様の構成のＧＰＳ受信部１０に対して、本発明によるＧＰＳサイクルスリップ自動補正装置２０を外部接続して構成する例を示している。従って、ＧＰＳ受信部１０は、方式のバリエーションは兎も角、従来のＧＰＳ受信機と同様の回路ブロックを備えている。
【００２３】
すなわち、ＧＰＳ受信機１０では、受信部１１にて、ＧＰＳ信号を受信し、周波数変換(down convert)して復調部１２へ出力し、復調部１２にて、原振・シンセサイザ１４からの信号によって搬送波をロックして相関器・位相カウンタ１３とＧＰＳサイクルスリップ自動補正装置２０とに出力する。相関器・位相カウンタ１３では、まず、相関器にて、拡散符号のレプリカとの相関を計算して、拡散符号をロックして、航法信号を復調して、位相カウンタにて、拡散符号のロックを保持するための時間差、搬送波のロックを保持するための位相差をそれぞれ計測して、バイアス擬似距離と搬送波位相とを出力する。
【００２４】
一方、ＧＰＳ受信機１０に外部接続するＧＰＳサイクルスリップ自動補正装置２０は、１次スプラインフィルタ機能を有するＧＰＳサイクルスリップ自動補正プログラムを格納して、該ＧＰＳサイクルスリップ自動補正プログラムを実行するＭＰＵ２１を少なくとも内蔵している。すなわち、ＧＰＳ受信機１０の復調部１２からの搬送波と相関器・位相カウンタ１３からの搬送波位相データ及び今までの挙動を示すデータとを入力として、ＭＰＵ２１にてＧＰＳサイクルスリップ自動補正プログラムを実行することにより、搬送波位相データに対してスプラインフィルタリング処理を施してサイクルスリップの有無を自動検出すると共に、サイクルスリップが検出された場合には、搬送波位相データを自動的に補正して、補正結果を外部に出力する。また、ＧＰＳサイクルスリップ自動補正プログラムによる演算結果は、原振・シンセサイザ１４にも出力され、次の受信搬送波のロックのために用いられる。
【００２５】
つまり、図１に例示するＧＰＳサイクルスリップ自動補正装置２０では、ＧＰＳ受信部１０の動作を制御する各種のデータの生成に係わるソフトウェアを動作させるが、このソフトウェアの一部にＧＰＳサイクルスリップ自動補正プログラムとして１次スプラインフィルタ処理を組み込んで、搬送波位相データに対する自動サイクルスリップ補正機能を持たせることができる。すなわち、ＭＰＵ２１に１次スプラインフィルタリングアルゴリズムを実行するＧＰＳサイクルスリップ自動補正プログラムを格納し、ＧＰＳ受信機１０から出力される搬送波位相データを１次スプラインフィルタによって補間処理して出力する。
【００２６】
（１次スプラインフィルタの動作）
次に、本発明によるＧＰＳサイクルスリップ自動補正プログラムに適用する１次スプラインフィルタのアルゴリズムについて、図２−図７に示すフローチャートを参照して説明する。ここに、図２は、本発明によるＧＰＳサイクルスリップ自動補正プログラムに適用する１次スプラインフィルタの動作の一例を説明するためのフローチャートであり、図３は、本発明によるＧＰＳサイクルスリップ自動補正プログラムの初期設定処理の一例を説明するためのフローチャートである。また、図４は、本発明によるＧＰＳサイクルスリップ自動補正プログラムの１次スプラインフィルタ動作としてデータ観測時刻まで伝搬させる観測時刻伝搬処理の一例を説明するためのフローチャートである。また、図５は、本発明によるＧＰＳサイクルスリップ自動補正プログラムにおいてサイクルスリップ発生の有無を判定するサイクルスリップ判定処理の一例を説明するためのフローチャートであり、図６は、本発明によるＧＰＳサイクルスリップ自動補正プログラムにおいてサイクルスリップによる搬送波位相データの補正を行うサイクルスリップ補正処理の一例を説明するためのフローチャートである。また、図７は、本発明によるＧＰＳサイクルスリップ自動補正プログラムにおける１次スプラインフィルタリング状態の更新を行うスプラインフィルタ状態更新処理の一例を説明するためのフローチャートである。
【００２７】
まず、図２のフローチャートを用いて、本発明によるＧＰＳサイクルスリップ自動補正プログラムの全体の流れの一例について説明する。最初に、１次スプラインフィルタ処理に使用するパラメータを初期設定する（ステップＳ１）。その後、１次スプラインフィルタの微分方程式を数値積分することにより解いて、データ観測時刻まで伝搬させる処理を行う（ステップＳ２）。
【００２８】
次に、１次スプラインフィルタの状態更新前の搬送波位相データの予測値と、実際に得られた搬送波位相データとから、実測搬送波位相データが基準誤差の範囲内に収まっているか否かを判定することにより、サイクルスリップの発生の有無を判定する（ステップＳ３）。サイクルスリップが発生していると判定した場合は（ステップＳ４のＹＥＳ）、サイクルスリップによる搬送波位相データのアンビグイティを補正し（ステップＳ５）、一方、サイクルスリップが発生していないと判定した場合は（ステップＳ４のＮＯ）、搬送波位相データの補正処理をスキップする。しかる後、次の観測データの１次スプラインフィルタ処理に備えるために、１次スプラインフィルタリングの状態を更新する（ステップＳ６）。しかる後、ステップＳ２の観測時刻への伝搬処理に復帰し、次の観測データに対するサイクルスリップの判定・補正処理へ移行する。
【００２９】
次に、図２で説明した各ステップの処理について、図３−図７のフローチャートを用いて、その一例を更に詳細に説明する。まず、図３のフローチャートに示す初期設定処理について説明する。
【００３０】
図３の初期設定処理において、まず、搬送波位相データのアンビギュイティ（ambiguity）を示す変数「Ａ」を０に初期設定する（ステップＳ１１）。次に、搬送波位相（位相変化率）のアプリオリ推定における期待値（予測値）を、変数「ｘ」に設定する（ステップＳ１２）。次に、１次スプラインフィルタの共分散行列のアプリオリ推定における共分散値を、変数「Ｓ」に設定する（ステップＳ１３）。次に、実際に観測される搬送波位相データに含まれる観測ノイズの共分散値を、変数「Ｒ」に設定する（ステップＳ１４）。
【００３１】
更に、搬送波位相データの飛躍量（補正量）を量子化するために用いる量子化単位（「２π」が１波長に相当）を、変数「Ｔ」に設定する（ステップＳ１５）。最後に、サイクルスリップの発生の有無を判定する際の安全係数を、ユーザが任意に定義して、変数「ｋ」に設定する（ステップＳ１６）。
【００３２】
次に、図４のフローチャートに示す１次スプラインフィルタの観測時刻伝搬処理について説明する。まず、次の式１に示す２つの微分方程式を、数値積分を用いて解く（ステップＳ２１）。
【００３３】
【数１】
2006275860.tif 000003
【００３４】
この結果を、搬送波位相データの観測時刻ｔまで伝搬させることにより、次の式２に示す搬送波位相値ｘ(ｔ)と共分散値Ｓ（ｔ）との双方の観測時刻ｔにおける数値を算出する（ステップＳ２２）。これにより、観測した搬送波位相データに関するサイクルスリップの発生の有無と補正量の算出とを迅速に行うことを可能としている。
【００３５】
【数２】
2006275860.tif 000004
【００３６】
次に、図５のフローチャートに示す１次スプラインフィルタにおけるサイクルスリップ判定処理について説明する。まず、サイクルスリップ判定用の基準誤差σを、次の式３を用いて算出する（ステップＳ３１）。
【００３７】
【数３】
2006275860.tif 000005
【００３８】
次に、図４の観測時刻伝搬処理にて算出した結果を次の式４に適用して、１次スプラインフィルタの状態更新前における観測残差Ｒｅｓを算出する（ステップＳ３２）。
【００３９】
【数４】
2006275860.tif 000006
【００４０】
ここで、式３で求めた基準誤差σとユーザが定義した安全係数ｋとの積と、観測残差Ｒｅｓの絶対値とを比較して、前者の値が後者以上の値であれば（ステップＳ３３のＹＥＳ）、サイクルスリップは発生していないものと判定して、その旨を出力し（ステップＳ３４）、逆に、前者の値が後者よりも小さい値であれば（ステップＳ３３のＮＯ）、サイクルスリップが発生しているものと判定して、その旨を出力する（ステップＳ３５）。
【００４１】
次に、図６のフローチャートに示す１次スプラインフィルタにおけるサイクルスリップ補正処理について説明する。まず、図５のサイクルスリップ判定処理にて算出した観測残差Ｒｅｓを、量子化単位Ｔで除算した結果を整数化（量子化）して、変数ｍに設定する（ステップＳ５１）。しかる後、変数ｍの量子化数値に量子化単位Ｔを乗算して、観測データを補正すべき補正量（ｍ・Ｔ）を算出して、観測した搬送波位相データのアンビグイティＡから減算することにより、補正後の搬送波位相データのアンビギュイティＡを求める（ステップＳ５２）。この処理により、サイクルスリップ発生時でも、１次スプラインフィルタによる補間処理により瞬時にほぼ正確なアンビギュイティＡを出力することができ、観測データの連続性を確保することができる。
【００４２】
次に、図７のフローチャートに示す１次スプラインフィルタにおけるスプラインフィルタ状態更新処理について説明する。まず、搬送波位相データの期待値ｘ（ｉ）を、次の式５に示すように、図４の観測時刻伝搬処理において算出された共分散相関に従って更新する（ステップＳ６１）。
【００４３】
【数５】
2006275860.tif 000007
【００４４】
次に、共分散行列Ｓを、次の式６を用いて更新し、次の観測データに関するサイクルスリップ判定用の準備を完了する（ステップＳ６２）。
【００４５】
【数６】
2006275860.tif 000008
【００４６】
以上のように、ＧＰＳサイクルスリップ自動補正装置２０に搭載のＭＰＵ２１は、１次スプラインフィルタ処理を含むＧＰＳサイクルスリップ自動補正プログラムを実行させて、ＧＰＳ受信部１０から読み込んだ搬送波位相データについて、該搬送波位相データ取得前の挙動に基づいたアプリオリ推定の予測値と比較を行い、両者の差異が、観測ノイズレベル及び１次スプラインフィルタの推定精度に比べて、有意であるか否かを判定し、差異が有意であると判定されたものについては、取得した該搬送波位相データの飛躍量に関する量子化条件から飛躍量を決定し、該搬送波位相データのアンビギュイティＡの値を補正することにより、データの連続性を確保し、高精度の測位・航法を継続的に実施することが可能である。
【００４７】
また、本実施例に示すように、サイクルスリップの検知及び補正に用いる手段としては、ＧＰＳ受信部で受信した搬送波データに関する１次スプラインフィルタ処理を行うＧＰＳサイクルスリップ自動補正装置２０のみを装備すればよく、従来技術のように、他の手段（例えば、前記特許文献１に記載のような慣性装置やジャイロなど）を別途装備する必要はなく、経済的に、誤り無く正確にサイクルスリップの発生を検知することができる。
【００４８】
（本発明の他の実施例）
更には、前述したように、１次スプラインフィルタリングアルゴリズムを含むＧＰＳサイクルスリップ自動補正プログラムを、ＧＰＳ受信機そのもののソフトウェアとして内蔵するか、図１に例示のように該ＧＰＳサイクルスリップ自動補正プログラムを装備したＧＰＳサイクルスリップ自動補正装置２０を従来のＧＰＳ受信機に外部接続して構成するか、あるいは、該ＧＰＳサイクルスリップ自動補正装置２０そのものをＧＰＳ受信機に内蔵させて構成するか、あるいは、１次スプラインフィルタ処理アルゴリズムをプログラム論理ではなく、ゲートアレイ論理などのワイヤードロジックで構成するか、など、補正後の搬送波位相データを出力するまでのＧＰＳ受信機の回路構成に関しては、種々のバリエーションを採用することが可能である。
本発明の本質は、観測した搬送波位相データに対して、１次スプラインフィルタを利用したサイクルスリップの検知及び補正機構を適用する手段を備えていることである。
【００４９】
なお、本発明は、例えば、ＧＰＳ受信機能を有する、カーナビや携帯電話機などを含む各種の小型の移動用機器に好適に適用することが可能であり、一例として斯かる移動用機器に内蔵されたソフトウェアに本発明によるＧＰＳサイクルスリップ自動補正プログラムを組み込むことにより、ＧＰＳ信号のサイクルスリップに対する迅速かつ正確な自動補正を経済的に実現することができる。
【００５０】
以上、本発明の好適実施例の構成を説明した。しかし、斯かる実施例は、本発明の単なる例示に過ぎず、何ら本発明を限定するものではないことに留意されたい。本発明の要旨を逸脱することなく、特定用途に応じて種々の変形変更が可能であることが、当業者には容易に理解できよう。
【図面の簡単な説明】
【００５１】
【図１】本発明によるＧＰＳ受信機の構成の一実施例を示す構成概念図である。
【図２】本発明によるＧＰＳサイクルスリップ自動補正プログラムに適用する１次スプラインフィルタの動作の一例を説明するためのフローチャートである。
【図３】本発明によるＧＰＳサイクルスリップ自動補正プログラムの初期設定処理の一例を説明するためのフローチャートである。
【図４】本発明によるＧＰＳサイクルスリップ自動補正プログラムの１次スプラインフィルタ動作としてデータ観測時刻まで伝搬させる観測時刻伝搬処理の一例を説明するためのフローチャートである。
【図５】本発明によるＧＰＳサイクルスリップ自動補正プログラムにおいてサイクルスリップ発生の有無を判定するサイクルスリップ判定処理の一例を説明するためのフローチャートである。
【図６】本発明によるＧＰＳサイクルスリップ自動補正プログラムにおいてサイクルスリップによる搬送波位相データの補正を行うサイクルスリップ補正処理の一例を説明するためのフローチャートである。
【図７】本発明によるＧＰＳサイクルスリップ自動補正プログラムにおける１次スプラインフィルタリング状態の更新を行うスプラインフィルタ状態更新処理の一例を説明するためのフローチャートである。
【符号の説明】
【００５２】
１      ＧＰＳ受信機
１０    ＧＰＳ受信部
１１    受信部
１２    復調部
１３    相関器・位相カウンタ
１４    原振・シンセサイザ
２０    ＧＰＳサイクルスリップ自動補正装置
２１    ＭＰＵ(Micro Processing Unit)
【図１】
2006275860.tif 000009
【図２】
2006275860.tif 000010
【図３】
2006275860.tif 000011
【図４】
2006275860.tif 000012
【図５】
2006275860.tif 000013
【図６】
2006275860.tif 000014
【図７】
2006275860.tif 000015