JP7329814B2 - 測位装置、測位方法、測位プログラム、測位プログラム記憶媒体、アプリケーション装置、および、測位システム - Google Patents

Description

本発明は、ＧＮＳＳの測位信号を用いて単独測位を行う測位技術に関する。

従来、単独測位を行う測位装置が各種考案されている。例えば、特許文献１に示すように、測位装置は、ＧＰＳ（Ｇｌｏｂａｌ Ｐｏｓｉｔｉｏｎｉｎｇ Ｓｙｓｔｅｍ）等のＧＮＳＳ（Ｇｒｏｂａｌ Ｎａｖｉｇａｔｉｏｎ Ｓａｔｅｌｌｉｔｅ Ｓｙｓｔｅｍｓ）の測位信号を受信し、擬似距離や搬送波位相積算値を用いて、位置の推定演算を行う。

この際、測位装置は、特許文献１に示すようなカルマンフィルタを用いて、位置の推定演算を行う。そして、カルマンフィルタを用いる従来の方法では、電離層遅延および対流圏遅延等の誤差要因によって、位置の推定精度が低下してしまう。このような電離層遅延および対流圏遅延等の誤差要因を、観測方程式から消去するため、従来では、測位装置は、衛星間一重位相差、受信機間一重位相差、または、衛星と受信機の二重位相差を用いていた。

特開２０１１－１３０３９９号公報

しかしながら、従来の構成および方法では、対となる複数の測位衛星から測位信号を継続的に受信する必要があったり、複数の受信機（アンテナ）を必要とする。

したがって、本発明の目的は、対となる複数の測位衛星からの測位信号を継続的に受信する必要が無く、少なくとも１つの受信機（アンテナ）によって、高精度な測位を実現する測位装置を提供することにある。

この発明の測位装置は、受信部、および、測位演算部を備える。受信部は、複数時刻において受信した測位信号の擬似距離および搬送波位相積算値を検出する。測位演算部は、擬似距離の時間差分値および搬送波位相積算値の時間差分値を算出する。測位演算部は、時間差分値を用いて、観測方程式および状態方程式を設定する。測位演算部は、観測方程式と状態方程式とを用いたフィルタ演算を実行することによって測位を行う。

この構成では、少なくとも１つの測位信号を複数の時刻で受信できていれば、観測方程式に含まれる未知数が少なくなり、測位精度が向上する。

また、この発明の測位装置は、次の構成であってもよい。測位演算部は、移動体の測位を行う場合に、前記測位する位置を、ＷＧＳ座標系から局地水平座標系へ座標変換する。測位演算部は、水平方向の動的モデルには、加速度の１次マルコフモデルを適用し、高さ方向の動的モデルには、速度の１次マルコフモデルを適用する。

この構成では、移動体に対する状態方程式が、より適切に設定される。

また、この発明の測位装置は、次の構成であってもよい。測位演算部は、水平方向の動的モデルに速度の項を含み、速度を、測位対象の速度として出力する。

この構成では、位置とともに速度が算出される。

また、この発明の測位装置は、次の構成であってもよい。測位演算部は、時間差分値算出部、二重差分値算出部、および、異常検出部を備える。時間差分値算出部は、擬似距離の時間差分値を複数時刻で算出する。二重差分値算出部は、複数時刻の擬似距離の時間差分値に対する二重時間差分値を算出する。異常検出部は、二重時間差分値が異常検出用閾値以上であれば、測位信号の受信の異常を検出する。

この構成では、観測異常、すなわち、測位信号の受信の異常が、容易に検出される。

また、この発明のアプリケーション装置は、上述のいずれかに記載の測位装置の各構成と、アプリケーション処理部と、を備える。アプリケーション処理部は、測位演算部で算出された位置を用いたアプリケーション情報を生成する。

この構成では、高精度な測位結果が得られることよって、より適切なアプリケーション情報が生成される。

また、この発明の測位システムは、基地局と測位衛星とを備える。基地局は、上述のいずれかに記載の測位装置の各構成、および、測位結果を用いて衛星軌道情報の補正情報を生成する補正情報生成部を備える。測位衛星は、補正情報生成部からの補正情報を測位信号とともに送信する。

この構成では、高精度な補正情報が得られることによって、より高精度な測位が可能になる。

この発明によれば、極簡素な構成および処理によって、高精度な測位を実現できる。

（Ａ）は、本発明の実施形態に係る測位装置の構成を示すブロック図であり、（Ｂ）は、測位演算部の構成を示すブロック図である。 本発明の実施形態に係る測位方法を示すフローチャートである。 異常検出を行う場合の測位演算部の構成を示すブロック図である。 観測値の異常検出のフローチャートである。 本発明の実施形態に係るナビゲーション装置の構成を示すブロック図である。 （Ａ）は、本発明の実施形態に係る補正システムの構成を示す図であり、（Ｂ）は、基地局の構成を示すブロック図である。

本発明の実施形態に係る測位装置および測位方法について、図を参照して説明する。図１（Ａ）は、本発明の実施形態に係る測位装置の構成を示すブロック図であり、図１（Ｂ）は測位演算部の構成を示すブロック図である。

図１（Ａ）に示すように、測位装置１０は、受信部３０、測位演算部４０、および、記憶部５０を備える。受信部３０および測位演算部４０は、それぞれに、１または複数の半導体素子やＩＣによって実現されている。なお、受信部３０および測位演算部４０は、情報処理装置と該情報処理装置によって実行されるプログラムによって実現できる。記憶部５０は、磁気デバイスや半導体デバイス等の既知の記憶媒体によって実現されている。なお、記憶部５０が半導体デバイスである場合には、記憶部５０は、測位演算部４０に含まれる構成であってもよい。

受信部３０は、アンテナ２０に接続されている。なお、アンテナ２０を測位装置１０に含めてもよい。

アンテナ２０は、複数の測位衛星ＳＡＴ（図１では１つの測位衛星のみを図示しており、他の測位衛星の図示は省略している。）からの測位信号ＳＳを受信し、受信部３０に出力する。測位信号ＳＳとは、例えば、既知のＧＰＳ（Ｇｒｏｂａｌ Ｐｏｓｉｔｉｏｎｉｎｇ Ｓｙｓｔｅｍ）信号等のＧＮＳＳ（Ｇｒｏｂａｌ Ｎａｖｉｇａｔｉｏｎ Ｓａｔｅｌｌｉｔｅ Ｓｙｓｔｅｍｓ）で用いられる測位信号である。

受信部３０は、測位信号ＳＳを捕捉、追尾して、測位信号ＳＳの擬似距離ρおよび搬送波位相積算値Φを検出する。受信部３０は、複数の時刻において、測位信号ＳＳの擬似距離ρおよび搬送波位相積算値Φを検出する。受信部３０は、各時刻での擬似距離ρおよび搬送波位相積算値Φを、測位演算部４０に出力する。

測位演算部４０は、擬似距離ρおよび搬送波位相積算値Φを記憶部５０に記憶する。図１（Ｂ）に示すように、測位演算部４０は、時間差分値算出部４１とフィルタ演算部４２とを備える。時間差分値算出部４１は、複数時刻の擬似距離ρの差分値（時間差分値）Δρを算出する。時間差分値算出部４１は、同じ複数時刻の搬送波位相積算値Φの差分値（時間差分値）ΔΦを算出する。

フィルタ演算部４２は、擬似距離ρの時間差分値Δρと、搬送波位相積算値Φの時間差分値ΔΦとを用いて、次に示す原理に則して状態方程式と観測方程式とを設定し、カルマンフィルタ処理を実行する。これにより、観測方程式に含まれる電離層遅延、対流圏遅延等の未知数は無くなり、測位精度が向上する。

なお、以下では、測位信号としてＧＰＳ信号のＬ１信号を用いる場合を、一例として説明する。

時刻ｔでの擬似距離ρ（ｔ）の観測モデルは、次の（式１）で表される。この観測モデルは、複数の測位衛星ＳＡＴのそれぞれに対して個別に設定される。

ρ（ｔ）＝ｒ（ｔ）＋δＩ（ｔ）＋δＴ（ｔ）＋ｃδｔ（ｔ）
＋δｂ_ＣＡ＋ｅ_ＣＡ（ｔ） （式１）
（式１）において、ｒ（ｔ）は測位装置１０と測位衛星ＳＡＴとの幾何学距離であり、δＩ（ｔ）は電離層遅延であり、δＴ（ｔ）は対流圏遅延であり、δｔ（ｔ）は測位装置１０および測位衛星ＳＡＴの時計誤差であり、δｂ_ＣＡは測位装置１０および測位衛星ＳＡＴのコードに対するバイアス誤差であり、ｅ_ＣＡ（ｔ）は観測誤差である。また、ｃは光速である。

また、時刻ｔでの搬送波位相積算値Φ（ｔ）の観測モデルは、次の（式２）で表される。この観測モデルは、複数の測位衛星ＳＡＴのそれぞれに対して個別に設定される。

Φ（ｔ）＝ｒ（ｔ）－δＩ（ｔ）＋δＴ（ｔ）＋ｃδｔ（ｔ）
＋δｂ_Ｌ１＋λＮ＋λε_Ｌ１（ｔ） （式２）
（式２）において、ｒ（ｔ）は測位装置１０と測位衛星ＳＡＴとの幾何学距離であり、δＩ（ｔ）は電離層遅延であり、δＴ（ｔ）は対流圏遅延であり、δｔ（ｔ）は測位装置１０および測位衛星ＳＡＴの時計誤差であり、δｂ_Ｌ１は測位装置１０および測位衛星ＳＡＴの搬送波位相に対するバイアス誤差であり、ＮはＬ１信号に対する整数値バイアスであり、ε_Ｌ１（ｔ）は観測誤差である。また、ｃは光速である。

ここで、時刻ｔ１と時刻ｔ２を設定し、時刻ｔ１と時刻ｔ２との時間間隔を短く設定する。例えば、時間間隔を約１０秒以下とする。なお、時間間隔は、測位衛星の速度および測位装置の概略の移動速度に応じて適宜設定すればよく、一例として３秒から５秒程度であるとよい。時刻ｔ１と時刻ｔ２との時間間隔が短いと、測位衛星および測位装置の移動距離は小さい。したがって、時刻ｔ１での電離層遅延δＩ（ｔ１）と、時刻ｔ２での電離層遅延δＩ（ｔ２）とは、一致しているものとして考えることが可能である。同様に、時刻ｔ１での対流圏遅延δＴ（ｔ１）と、時刻ｔ２での対流圏遅延δＴ（ｔ２）とは、一致しているものとみなせる。

したがって、電離層遅延δＩ（ｔ１）と電離層遅延δＩ（ｔ２）との差分値は略０である。同様に、対流圏遅延δＴ（ｔ１）と対流圏遅延δＴ（ｔ２）との差分値も略０である。

また、時刻ｔ１と時刻ｔ２との時間間隔が短く、測位衛星および測位装置の移動距離は小さいので、時刻ｔ１での整数値バイアスと時刻ｔ２での整数値バイアスとは同じとして考えることができる。

また、測位装置および測位衛星の時計誤差やバイアス誤差は、時刻ｔ１と時刻ｔ２で変化しないと考えることができる。

したがって、時刻ｔ１での擬似距離ρ（ｔ１）と時刻ｔ２での擬似距離ρ（ｔ２）との時間差分値Δρ（ｔ１，ｔ２）は、次の（式３）で近似できる。

Δρ（ｔ１－ｔ２）＝ｒ（ｔ１）－ｒ（ｔ２）＋Ｅ_ＣＡ（ｔ１，ｔ２） （式３）
なお、Ｅ_ＣＡ（ｔ１，ｔ２）は、擬似距離の時間差分値Δρに対するコード位相の観測誤差である。

同様に、時刻ｔ１での搬送波位相積算値Φ（ｔ１）と時刻ｔ２での搬送波位相積算値Φ（ｔ２）との時間差分値ΔΦ（ｔ１，ｔ２）は、次の（式４）で近似できる。

ΔΦ（ｔ１－ｔ２）＝ｒ（ｔ１）－ｒ（ｔ２）＋Ｅ_Ｌ１（ｔ１，ｔ２） （式４）
なお、Ｅ_Ｌ１（ｔ１，ｔ２）は、搬送波位相積算値の時間差分値ΔΦに対する搬送波位相の観測誤差である。

ここで、ｒ（ｔ）は、測位装置１０と測位衛星ＳＡＴとの距離であり、測位装置１０の位置をｕ（ｔ）、測位衛星ＳＡＴの位置をＳｐ（ｔ）とし、距離演算を示す記号をＤＩＳとすると、次の（式５）となる。

ｒ（ｔ）＝ＤＩＳ（ｕ（ｔ）－Ｓｐ（ｔ）） （式５）
測位衛星ＳＡＴの位置Ｓｐ（ｔ）は、測位信号に重畳された航法メッセージ、他の精密軌道情報が提供されるシステム（例えば、ＭＡＤＯＣＡ）からの情報等によって取得が可能である。

したがって、上述の（式３）および（式４）に、（式５）を代入することによって、測位装置１０の位置ｕ（ｔ）に対する、擬似距離の時間差分値Δρおよび搬送波位相差積算値の時間差分値ΔΦを用いた観測モデルを設定できる。そして、この観測モデルを用いることによって、既知の式変形から、測位装置１０の位置に対する観測方程式を設定できる。

また、状態方程式は、次のように設定できる。

（Ａ）静止状態（Ｓｔａｔｉｃ）の場合
静止状態の場合、測位装置１０は移動していない。したがって、ｕ（ｔ_ｋ）＝ｕ（ｔ_ｋ－１）となり、状態方程式を設定できる。

したがって、静止状態に対して、カルマンフィルタを適用でき、測位装置１０の位置ｕ（ｔ）を推定演算できる。なお、観測誤差が白色雑音で仮定できる場合は、一般的なカルマンフィルタの理論を用いればよく、観測誤差が有色雑音の場合は、状態・観測雑音に相関のあるカルマンフィルタの理論を用いればよい。

（Ｂ）移動状態（ｋｉｎｅｍａｔｉｃ）の場合
移動状態の場合、測位装置１０の移動モデルを、加速度または速度を用いた１次近似モデルで設定する。具体的には、例えば、次のように設定する。まず、位置ｕ（ｔ）をＷＧＳ座標系からＥＮＵ座標系に変換する。この座標変換は、既知の方法によって実現可能である。

そして、ＥＮＵ座標系におけるＥａｓｔ方向とＮｏｒｔｈ方向の移動モデルは、加速度１次マルコフモデルに設定する。また、ＥＮＵ座標系におけるＵｐ方向の移動モデルは、速度１次マルコフモデルに設定する。

これにより、移動状態に対する状態方程式を設定できる。

したがって、移動状態に対して、カルマンフィルタを適用でき、測位装置１０の位置ｕ（ｔ）を推定演算できる。

以上のような構成および処理を用いることによって、観測方程式における、電離層遅延、対流圏遅延、時計誤差、およびバイアス誤差の誤差要因を無くすことができ、推定精度を向上できる。また、１重位相差、２重位相差を用いないことによって、演算が簡素化される。さらに、少なくとも１つのアンテナで測位演算が可能である。

なお、移動状態の場合においては、カルマンフィルタに用いる二時刻ｔ_ｋ－１、ｔ_ｋの時間差が極小さいと、距離差ｒ（ｔ_ｋ）－ｒ（ｔ_ｋ－１）に相当する値を算出するための位置ｕ（ｔ_ｋ）、ｕ（ｔ_ｋ－１）に乗算する勾配ベクトルＧ（ｔ_ｋ－１）、Ｇ（ｔ_ｋ）が殆ど変化しない。この場合、Ｇ（ｔ_ｋ）－Ｇ（ｔ_ｋ－１）によって決まるカルマンゲインが小さくなり、推定誤差が大きくなってしまうことがある。

この場合、位置差ｕ（ｔ_ｋ）－ｕ（ｔ_ｋ－１）を未知数として、勾配ベクトルＧ（ｔ_ｋ－１）、Ｇ（ｔ_ｋ）が殆ど同じであることを利用して、勾配ベクトルの平均値（Ｇ（ｔ_ｋ－１）＋Ｇ（ｔ_ｋ））／２を用いて観測方程式を設定すればよい。これにより、安定して、高精度な測位が可能になる。

また、移動状態の場合、状態方程式に含まれる速度の項を用いることによって、測位装置１０は、速度を算出することも可能である。

また、上述の構成では、１つのアンテナを用いる態様を示したが、複数のアンテナを用いてもよい。この場合、測位装置１０は、複数のアンテナのそれぞれに対して、上述の測位演算を実行する。そして、測位装置１０は、これらの複数のアンテナの位置に対する平均値等の統計処理を行ってもよい。これにより、位置の推定精度をさらに向上することが可能になる。

また、上述の説明では、各処理をそれぞれの機能部で分担して実行する態様を示した。しかしながら、上述の各処理をプログラム化して記憶部５０に記憶しておき、ＣＰＵ等の演算装置で実行してもよい。この場合、図２に示す処理を実行すればよい。図２は、本発明の実施形態に係る測位方法を示すフローチャートである。

演算装置は、第１時刻ｔ１での擬似距離ρ１と搬送波位相積算値Φ１を検出する（Ｓ１０１）。演算装置は、第２時刻ｔ２での擬似距離ρ２と搬送波位相積算値Φ２を検出する（Ｓ１０２）。

演算装置は、第１時刻ｔ１の擬似距離ρ１と第２時刻ｔ２の擬似距離ρ２との時間差分値Δρ１２＝ρ１－ρ２を算出する。また、演算装置は、第１時刻ｔ１の搬送波位相積算値Φ１と第２時刻ｔ２の搬送波位相積算値Φ２との時間差分値ΔΦ１２＝Φ１－Φ２を算出する（Ｓ１０３）。

演算装置は、擬似距離の時間差分値Δρ１２、および、搬送波位相積算値の時間差分値ΔΦ１２を観測値として、測位装置の位置を推定するための状態方程式および観測方程式を設定する（Ｓ１０４）。

演算装置は、ステップＳ１０４で設定した状態方程式および観測方程式に対してカルマンフィルタ処理を実行し、測位装置の位置を算出する（Ｓ１０５）。

なお、上述の測位装置および測位方法では、観測値に異常検出も可能である。図３は、異常検出を行う場合の測位演算部の構成を示すブロック図である。

測位演算部４０Ａは、時間差分値算出部４１、二重差分値算出部４３、および、異常検出部４４を備える。時間差分値算出部４１は、複数時刻の擬似距離ρの差分値（時間差分値）Δρを算出する。二重差分値算出部４３は、異なる時刻で得られた時間差分値Δρａと時間差分値Δρｂとの二重差分値ΔΔρａｂを算出する。異常検出部４４は、二重差分値Δρａｂを異常検出用の閾値Ｔｈと比較する。異常検出部４４は、二重差分値Δρａｂが異常検出用の閾値Ｔｈを超えていれば、観測値の異常、すなわち、測位信号の受信の異常を検出する。

なお、この異常検出処理は、上述のように演算処理装置とプログラムによって実現することもできる。図４は、観測値の異常検出のフローチャートである。

測位装置１０の測位演算部４０は、擬似距離の時間差分値Δρの更なる時間差分値である二重差分値ΔΔρを算出する（Ｓ２０２）。具体的には、測位演算部４０は、上述のように、第１時刻ｔ１の擬似距離ρ１と第２時刻ｔ２の擬似距離ρ２との時間差分値Δρ１２を算出する。また、測位演算部４０は、同様に、第２時刻ｔ２の擬似距離ρ２と第３時刻ｔ３の擬似距離ρ３との時間差分値Δρ２３を算出する。測位演算部４０は、時間差分値Δρ２３と時間差分値Δρ１２との二重差分値ΔΔρを算出する。なお、この説明では、２つの時間差分値に同じ時刻の擬似距離が含まれているが、これに限るものではない。

測位演算部４０は、二重差分値ΔΔρを、異常検出用の閾値Ｔｈと比較する（Ｓ２０２）。異常検出用の閾値Ｔｈは、既に取得している測位結果等を参照して、必要とする測位精度よりも劣化したときの時間差分値の変化量に基づいて設定されている。

測位演算部４０は、二重差分値ΔΔρが閾値Ｔｈを超えれば（Ｓ２０２：ＹＥＳ）、観測値である擬似距離が異常であると判定する（Ｓ２０３）。

測位演算部４０は、二重差分値ΔΔρが閾値Ｔｈを超えなければ（Ｓ２０２：ＮＯ）、観測値である擬似距離が正常であると判定する（Ｓ２０４）。

このように、本実施形態の構成および処理を用いることによって、観測値の異常を検出できる。なお、この説明では、擬似距離の場合を示したが、搬送波位相積算値を用いても異常の検出は可能である。

上述の測位装置は、自装置の位置を用いた各種のアプリケーション装置に適用できる。例えば、図５は、本発明の実施形態に係るナビゲーション装置の構成を示すブロック図である。

図５に示すように、ナビゲーション装置１は、測位装置１０、アンテナ２０、および、ナビゲーション処理部６１を備える。ナビゲーション装置１が、本発明の「アプリケーション装置」に対応し、ナビゲーション処理部６１が、本発明の「アプリケーション処理部」に対応する。

測位装置１０およびアンテナ２０は、上述の構成を備える。測位演算部４０は、測位結果を、ナビゲーション処理部６１に出力する。

ナビゲーション処理部６１は、測位結果を適用した既知の方法から、自装置が搭載される移動体の目的位置までのルート照会、ルート案内、または、現在位置の周辺情報等のナビゲーション情報を生成する。ナビゲーション処理部６１は、ナビゲーション情報を、図示しない表示部等に表示する。ここで、本実施形態の測位装置１０を用いることによって、自装置位置を高精度に得られるので、ナビゲーション装置１は、高精度なルート照会、ルート案内や、より適切な周辺情報を提供できる。

上述の説明では、測位装置の位置を算出する態様を示したが、測位装置を基地局に備えて、上記の推定演算を用いて、測位衛星の位置を高精度に算出することも可能である。そして、この算出結果を用いることによって、衛星軌道情報に対する補正情報を算出することもできる。図６（Ａ）は、本発明の実施形態に係る補正システムの構成を示す図であり、図６（Ｂ）は、基地局の構成を示すブロック図である。

図６（Ａ）、図６（Ｂ）に示すように、基地局２は、測位装置１０、アンテナ２０、送信用アンテナ２０Ｔ、補正情報生成部７１、および、送信部７２を備える。アンテナ２０および測位装置１０は、上述の構成を備える。測位演算部４０は、測位結果を、補正情報生成部７１に出力する。

補正情報生成部７１は、基地局２の位置（補正用の基準位置）を予め記憶している。補正情報生成部７１は、補正用の基準位置と測位位置とを比較する。補正情報生成部７１は、補正用の基準位置と測位位置との比較結果から、既知の方法を用いて、衛星位置の誤差を算出する。補正情報生成部７１は、この衛星位置の誤差から衛星軌道情報の補正情報を生成する。補正情報生成部７１は、補正情報を送信部７２に出力する。

送信部７２は、送信用アンテナ２０Ｔを介して、補正情報を測位衛星ＳＡＴに送信する。測位衛星ＳＡＴは、受信した補正情報によって、補正情報の更新を行う。

このようなシステムの構成を用いることによって、補正情報を高精度に更新でき、測位精度が向上する。

１：ナビゲーション装置
２：基地局
１０：測位装置
２０：アンテナ
３０：受信部
４０：測位演算部
４１：時間差分値算出部
４２：フィルタ演算部
４３：二重差分値算出部
４４：異常検出部
５０：記憶部
６１：ナビゲーション処理部
７１：補正情報生成部
７２：送信部
２０Ｔ：送信用アンテナ
ＳＡＴ：測位衛星
ＳＳ：測位信号

Claims (9)

1. 整数値バイアスの値が同じ複数時刻において受信した測位信号の擬似距離および搬送波位相積算値を検出する受信部と、
   前記擬似距離の時間差分値および前記搬送波位相積算値の時間差分値を用いて、観測方程式および状態方程式を設定し、前記観測方程式と前記状態方程式とを用いたフィルタ演算を実行することによって測位を行う、測位演算部と、
   を備える、測位装置。
2. 前記測位演算部は、
   移動体の測位を行う場合に、前記測位する位置を、ＷＧＳ座標系から局地水平座標系へ座標変換し、
   水平方向の動的モデルには、加速度の１次マルコフモデルを適用し、高さ方向の動的モデルには、速度の１次マルコフモデルを適用し、
   前記観測方程式において、前記複数時刻の位置差を未知数として、前記複数時刻の位置の勾配ベクトルの平均値を用いる、
   請求項１に記載の測位装置。
3. 前記測位演算部は、
   前記水平方向の動的モデルに速度の項を含み、
   該速度を、測位対象の速度として出力する、
   請求項２に記載の測位装置。
4. 前記測位演算部は、
   前記擬似距離の時間差分値を複数時刻で算出する時間差分値算出部と、
   該複数時刻の擬似距離の時間差分値に対する二重時間差分値を算出する二重差分値算出部と、
   該二重時間差分値が、異常検出用閾値以上であれば、前記測位信号の受信の異常を検出する異常検出部と、
   を備える、
   請求項１乃至請求項３のいずれかに記載の測位装置。
5. 請求項１乃至請求項４のいずれかに記載の測位装置の各構成と、
   前記測位演算部で算出された位置を用いたアプリケーション情報を生成するアプリケーション処理部と、
   を備える、アプリケーション装置。
6. 請求項１乃至請求項４のいずれかに記載の測位装置の各構成、および、前記測位の結果を用いて衛星軌道情報の補正情報を生成する補正情報生成部を備える基地局と、
   前記補正情報生成部からの補正情報を前記測位信号とともに送信する測位衛星と、
   を備える、測位システム。
7. 整数値バイアスの値が同じ複数時刻において受信した測位信号の擬似距離および搬送波位相積算値を検出し、
   前記擬似距離の時間差分値および前記搬送波位相積算値の時間差分値を算出し、
   前記時間差分値を用いた観測方程式および状態方程式を設定し、
   前記観測方程式と前記状態方程式とを用いたフィルタ演算を実行することによって測位を行う、
   測位方法。
8. 整数値バイアスの値が同じ複数時刻において受信した測位信号の擬似距離および搬送波位相積算値を検出する処理と、
   前記擬似距離の時間差分値および前記搬送波位相積算値の時間差分値を算出する処理と、
   前記時間差分値を用いた観測方程式および状態方程式を設定する処理と、
   前記観測方程式と前記状態方程式とを用いたフィルタ演算を実行することによって測位を行う処理と、
   を、演算処理装置に実行させる、測位プログラム。
9. 整数値バイアスの値が同じ複数時刻において受信した測位信号の擬似距離および搬送波位相積算値を検出する処理と、
   前記擬似距離の時間差分値および前記搬送波位相積算値の時間差分値を算出する処理と、
   前記時間差分値を用いた観測方程式および状態方程式を設定する処理と、
   前記観測方程式と前記状態方程式とを用いたフィルタ演算を実行することによって測位を行う処理と、
   を、演算処理装置に実行させるプログラムが記憶された、測位プログラム記憶媒体。

Images