JP2015021900A - 衛星時計補正パラメータ生成装置、その方法及び測位装置 - Google Patents
衛星時計補正パラメータ生成装置、その方法及び測位装置 Download PDFInfo
- Publication number
- JP2015021900A JP2015021900A JP2013151945A JP2013151945A JP2015021900A JP 2015021900 A JP2015021900 A JP 2015021900A JP 2013151945 A JP2013151945 A JP 2013151945A JP 2013151945 A JP2013151945 A JP 2013151945A JP 2015021900 A JP2015021900 A JP 2015021900A
- Authority
- JP
- Japan
- Prior art keywords
- satellite
- satellite clock
- correction parameter
- clock correction
- gnss
- 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.)
- Pending
Links
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
【課題】エフェメリス・データが正常に取得できず、前回受信したエフェメリス・データの有効期限が切れた場合、従来より精度の高い測位が行える補正パラメータを得ることができる衛星時計補正パラメータ生成装置を提供する。【解決手段】本発明の衛星時計補正パラメータ生成装置は、GNSSにおける測位装置に設けられ、GNSS衛星から衛星時計補正パラメータが正常に受信できなかった場合に、測位に用いるGNSS衛星と測位装置との観測距離を求める際に使用するGNSS衛星の衛星時計の誤差を求める衛星時計補正パラメータを推定する装置であり、過去の衛星時計補正パラメータであるバイアス値及び誤差値の変化率が書き込まれて記憶されている衛星時計補正パラメータ記憶部と、過去の複数のバイアス値から求めた傾きを変化率の推定値である推定変化率として出力する衛星時計変化率算出部とを備える。【選択図】図1
Description
本発明は、衛星時計補正パラメータ生成装置、その方法及び測位装置に関する。
GNSS(Global Navigation Satellite System:グローバル・ナビゲーション・サテライト・システム)は、地球上空を高速で飛行する人工衛星を利用して、このGNSS衛星から送出される電波を、GNSS受信機にて受信することにより、自己の位置(GNSS受信機の位置)を高精度で測位することができる位置検出システムである。このGNSS測位に用いられるシステムとしては、GPS(Global Positioning System)、GLONASS(Global Navigation Satellite System)及びGALILEO(ガリレオ、EU(European Union)の全地球航法衛星システム)などがある。
例えば、GPSシステムの場合、合計32個のGPS衛星が、地球上空約2万kmの6つの円軌道上に、5機及び6機ずつ配置され、地球上の如何なる場所においても、常に数個のGPS衛星からの電波を受信することができ、これら数個のGPS衛星から同時にGPS電波を受信することにより、位置検出装置側では自己の正確な位置を確認できるようになっている。
GNSSは、上述したGPS等を含んでおり、全地球的な位置測定、測量などに用いられている。地球上空の軌道を周回するGNSS衛星の各々は、アルマナック(Almanac)データとエフェメリス(Ephemeris)・データの二種類の軌道情報を含む航法メッセージを送信している。ここで、GNSS衛星とは、GPS衛星、GLONASS衛星(GLONASSで用いられる衛星)及びGALILEO衛星(GALILEOで用いられる衛星)などを含む、GNSS測位において用いられる衛星である。
カーナビゲーション装置などのGNSS受信機(測位装置)は、この航法メッセージを用いてGNSS衛星からGNSS受信機までの信号到達時間を測定する。そしてGNSS受信機は、信号到達時間に光速を掛けることでGNSS衛星とGNSS受信機間の距離(擬似距離)を求めている。
この作業を4機以上の衛星に対して行うことで、受信機の位置を三次元上で算出している。アルマナック・データは、軌道上における全ての衛星に関する軌道情報であり、少なくとも6日に1回、データが更新される。エフェメリス・データは、各衛星の正確な位置情報であり、2時間毎にデータが更新される。信号を発射した時刻情報は、6秒毎にデータ更新される。この時刻情報は、衛星に搭載されている高精度の原子時計(以下、衛星時計とする)により得られた時刻を示している。また、エフェメリス・データには、上記衛星時計の時間を補正する補正パラメータも含まれている。
この作業を4機以上の衛星に対して行うことで、受信機の位置を三次元上で算出している。アルマナック・データは、軌道上における全ての衛星に関する軌道情報であり、少なくとも6日に1回、データが更新される。エフェメリス・データは、各衛星の正確な位置情報であり、2時間毎にデータが更新される。信号を発射した時刻情報は、6秒毎にデータ更新される。この時刻情報は、衛星に搭載されている高精度の原子時計(以下、衛星時計とする)により得られた時刻を示している。また、エフェメリス・データには、上記衛星時計の時間を補正する補正パラメータも含まれている。
GNSS受信機は、受信した各情報を自機の記憶部に、例えば電源を切る前に受信したアルマナック・データとエフェメリス・データを保存している。これらの情報には、各々有効期限がある。エフェメリス・データは、有効期限が受信してから2時間から4時間程度とされている。しかし、エフェメリス・データは、30秒の長さのデータで構成されており、周囲の建物などの環境によって、電波の受信が阻害されて受信に時間を要する。GNSS受信機の電源がオフ状態で、GNSS衛星からの送信情報をエフェメリス・データの有効時間(例えば4時間)以上で受信できていない場合など、GNSS受信機は、再度、送信情報を受信する必要がある。エフェメリス・データの有効期限が切れたことにより再受信している期間、GNSS受信機は、位置測定を行えない。このため、エフェメリス・データの有効期限の大幅な延長が重要な課題となっている。
エフェメリス・データが受信できず、エフェメリス・データの有効期限切れにより再受信している期間、GNSS受信機は、位置測定を行えない。このため、エフェメリス・データの有効期限が切れた場合、GNSS受信機は、衛星位置については、以前受信したエフェメリス・データを用いて、軌道予測を行うことで算出し、衛星時計については、以前受信したエフェメリス・データを用いて、現在のGNSS衛星の衛星時計の時間を補正し、GNSS衛星とGNSS受信機との距離を求め、GNSS受信機の位置測定を行う方法が提案されている(例えば、特許文献1参照)。
しかしながら、特許文献1などの従来の手法において、衛星時計については、最後に取得した古い補正パラメータ(過去に取得した衛星時計補正パラメータ)を使用しており、現時点の補正パラメータと異なる数値のため、上述したGNSS衛星とGNSS受信機との距離の計算精度が低下してしまう。
このため、従来の手法においては、有効期限が切れた後に新たなエフェメリス・データが得られない場合、古い補正パラメータを使用することにより、GNSS衛星とGNSS受信機との距離が正確に求められず、GNSS受信機の測位の精度が大幅に低下してしまう。
このため、従来の手法においては、有効期限が切れた後に新たなエフェメリス・データが得られない場合、古い補正パラメータを使用することにより、GNSS衛星とGNSS受信機との距離が正確に求められず、GNSS受信機の測位の精度が大幅に低下してしまう。
本発明は、このような事情に鑑みてなされたもので、エフェメリス・データが正常に取得できず、前回受信したエフェメリス・データの有効期限が切れた場合においても、従来に比較してより精度の高い位置測定を行うために用いる補正パラメータを得ることができる衛星時計補正パラメータ生成装置、その方法及び測位装置を提供することを目的とする。
本発明の衛星時計補正パラメータ生成装置は、GNSSにおける測位装置に設けられ、GNSS衛星から衛星時計補正パラメータが正常に受信できなかった場合に、測位に用いる前記GNSS衛星と前記測位装置との観測距離を求める際に用いる前記GNSS衛星の衛星時計の誤差を求めるための衛星時計補正パラメータを推定する衛星時計補正パラメータ生成装置であり、過去の衛星時計補正パラメータであるバイアス値及び誤差値の変化率が書き込まれて記憶されている衛星時計補正パラメータ記憶部と、過去の複数の前記バイアス値から求めた傾きを前記変化率の推定値である推定変化率として出力する衛星時計変化率算出部とを備えることを特徴とする。
本発明の衛星時計補正パラメータ生成装置は、前記バイアス値を2つ用いる場合、前記衛星時計変化率算出部が、前記GNSS衛星から衛星時計補正パラメータを正常に受信できなかった場合に、直前とその前のバイアス値を結ぶ線分の傾き前記推定変化率とすることを特徴とする。
本発明の測位装置は、上記に記載の衛星時計補正パラメータ生成装置を備える測位装置であり、GNSS衛星から衛星時計補正パラメータが正常に受信できなかった場合に前記観測距離を求める際、前記衛星時計の時刻誤差を、直前に正常に受信した前記バイアス値と、前記衛星時計補正パラメータ生成装置が推定した前記推定変化率とを用いて算出することを特徴とする。
本発明の衛星時計補正パラメータ生成方法は、GNSSにおける測位装置において、GNSS衛星から衛星時計補正パラメータが正常に受信できなかった場合に、測位に用いる前記GNSS衛星と前記測位装置との観測距離を求める際に用いる前記GNSS衛星の衛星時計の誤差を求めるための衛星時計補正パラメータを推定する衛星時計補正パラメータ生成方法であり、衛星時計変化率算出部が、過去の衛星時計補正パラメータであるバイアス値及び誤差値の変化率が書き込まれて記憶されている衛星時計補正パラメータ記憶部から、過去の複数の前記バイアス値から求めた傾きを前記変化率の推定値である推定変化率として出力する衛星時計変化率算出過程を有することを特徴とする。
この発明によれば、過去に正常に受信された複数の衛星時計補正パラメータのバイアス値の変化より、誤差値の変化率を推定するため、従来のように衛星時計補正パラメータを取得できない場合に、直前の変化率を代用する従来に比較し、より衛星時計の誤差の変化を反映することができるため、GNSS衛星の観測距離をより高い精度で求めることが可能となる。
また、本実施形態によれば、GNSS衛星の観測距離を、従来に比較してより高い精度で求めることができるため、従来に比較して、測位装置の測位の結果も高い精度で得ることが可能となる。
また、本実施形態によれば、GNSS衛星の観測距離を、従来に比較してより高い精度で求めることができるため、従来に比較して、測位装置の測位の結果も高い精度で得ることが可能となる。
以下、図面を用いて、本発明の実施形態について説明する。
まず、一般的なGNSS衛星の軌道予測の概略原理について、GPS(Global Positioning System)衛星の例にとり説明する。
地球上空の6つの軌道を24機以上のGPS衛星が周回し、各衛星は、航法メッセージを送信している。航法メッセージは、フレームを単位に送信している。1フレームの送信には、30秒を要する。1フレームは、5組のサブフレームで構成され、サブフレーム1から5が順番に送信される。また、サブフレーム1から3は、各衛星に固有の情報を含み、毎回同じ内容で構成される。サブフレーム4と5は、全衛星が同じ内容の情報で構成され、サブフレーム毎にページ1から25で構成されている。このため、サブフレーム4と5の全ての情報を送信するには、25フレームが必要である。このため、GPS受信機(測位装置)が航法メッセージの全情報を得るには、12分30秒の時間がかかる。サブフレーム1は、各衛星の状態(正常に動作しているか否か)、衛星が送信する衛星のクロック誤差を補正するための係数であるクロック補正係数などにより構成されている。サブフレーム2は、各衛星の軌道情報(エフェメリス(Ephemeris)・データ)1/2で構成されている。
サブフレーム3は、各衛星の軌道情報(エフェメリス・データ)2/2で構成されている。サブフレーム4は、GPS受信機が受信する信号が、電離層により遅延する量を補正するための係数である電離層遅延補正係数、大気により遅延する量を補正するための係数である大気遅延補正係数、GPS時刻(GNSS時刻)とUTC(協定世界時間;Universal Time, Coordinated))との関係を表す情報であるUTC関係情報、全衛星の軌道情報(アルマナック(Almanac)・データ)1/2などで構成されている。サブフレーム5は、全衛星の軌道情報(アルマナック・データ)2/2で構成されている。また、各サブフレームの先頭には、GPS時刻を示す情報が含まれている。なお、GPS時刻とは、1週間を単位としてGPS衛星で管理されている時間(衛星時計のカウントしている時間)であり、毎週、日曜日の0時からの経過時間で表される情報である。
まず、一般的なGNSS衛星の軌道予測の概略原理について、GPS(Global Positioning System)衛星の例にとり説明する。
地球上空の6つの軌道を24機以上のGPS衛星が周回し、各衛星は、航法メッセージを送信している。航法メッセージは、フレームを単位に送信している。1フレームの送信には、30秒を要する。1フレームは、5組のサブフレームで構成され、サブフレーム1から5が順番に送信される。また、サブフレーム1から3は、各衛星に固有の情報を含み、毎回同じ内容で構成される。サブフレーム4と5は、全衛星が同じ内容の情報で構成され、サブフレーム毎にページ1から25で構成されている。このため、サブフレーム4と5の全ての情報を送信するには、25フレームが必要である。このため、GPS受信機(測位装置)が航法メッセージの全情報を得るには、12分30秒の時間がかかる。サブフレーム1は、各衛星の状態(正常に動作しているか否か)、衛星が送信する衛星のクロック誤差を補正するための係数であるクロック補正係数などにより構成されている。サブフレーム2は、各衛星の軌道情報(エフェメリス(Ephemeris)・データ)1/2で構成されている。
サブフレーム3は、各衛星の軌道情報(エフェメリス・データ)2/2で構成されている。サブフレーム4は、GPS受信機が受信する信号が、電離層により遅延する量を補正するための係数である電離層遅延補正係数、大気により遅延する量を補正するための係数である大気遅延補正係数、GPS時刻(GNSS時刻)とUTC(協定世界時間;Universal Time, Coordinated))との関係を表す情報であるUTC関係情報、全衛星の軌道情報(アルマナック(Almanac)・データ)1/2などで構成されている。サブフレーム5は、全衛星の軌道情報(アルマナック・データ)2/2で構成されている。また、各サブフレームの先頭には、GPS時刻を示す情報が含まれている。なお、GPS時刻とは、1週間を単位としてGPS衛星で管理されている時間(衛星時計のカウントしている時間)であり、毎週、日曜日の0時からの経過時間で表される情報である。
また、エフェメリス・データは、衛星の位置を計算するために必要な軌道の6要素(昇交点赤経、軌道傾斜角、近地点引数、軌道長半径、離心率、真近点角)のデータ、各補正値、及び軌道のエポック時刻toe(軌道の元期(ephemeris reference time))などにより構成されている。エフェメリス・データは、2時間毎に更新される。また、エフェメリス・データの有効期間は、受信時より2時間から4時間である。
アルマナック・データは、各衛星の軌道情報である。また、アルマナック・データは、少なくとも6日に1度の頻度で更新されている。また、アルマナック・データの有効期間は、6日である。
アルマナック・データは、各衛星の軌道情報である。また、アルマナック・データは、少なくとも6日に1度の頻度で更新されている。また、アルマナック・データの有効期間は、6日である。
GPS受信機は、受信したこれらのデータを記憶部に記憶する。記憶されているアルマナック・データ及びエフェメリス・データ共に有効な状態から、GPS受信機が測位を開始することをホット・スタートという。また、エフェメリス・データが無効でありアルマナック・データのみ有効な状態から、GPS受信機が測位を開始することをウォーム・スタートという。アルマナック・データ及びエフェメリス・データ共に無効な状態から、GPS受信機が測位を開始することをコールド・スタートという。
そして、ホット・スタートの場合、GPS受信機(測位装置)は、電源投入後に最初の位置情報が出力されるまでの時間は、記憶されているエフェメリス・データを用いて測位できるので数秒程度である。一方、ウォーム・スタートの場合、GPS受信機は、エフェメリス・データを再取得する必要があるので、取得に30秒以上の時間がかかり、さらに取得したエフェメリス・データの検査にも時間がかかる。このため、ウォーム・スタートの場合、電源投入後に最初の位置情報が出力されるまでの時間は、30秒以上である。さらに、コールド・スタートの場合、アルマナック・データが無効なため,現時刻における可視衛星を知ることができず、受信可能である衛星をランダムに探す必要が有る。この衛星を探す時間に加え、エフェメリス・データを受信するための時間を要するため、電源投入後に最初の位置情報が出力されるまでの時間は、ウォーム・スタートよりも長くなる。
そして、ホット・スタートの場合、GPS受信機(測位装置)は、電源投入後に最初の位置情報が出力されるまでの時間は、記憶されているエフェメリス・データを用いて測位できるので数秒程度である。一方、ウォーム・スタートの場合、GPS受信機は、エフェメリス・データを再取得する必要があるので、取得に30秒以上の時間がかかり、さらに取得したエフェメリス・データの検査にも時間がかかる。このため、ウォーム・スタートの場合、電源投入後に最初の位置情報が出力されるまでの時間は、30秒以上である。さらに、コールド・スタートの場合、アルマナック・データが無効なため,現時刻における可視衛星を知ることができず、受信可能である衛星をランダムに探す必要が有る。この衛星を探す時間に加え、エフェメリス・データを受信するための時間を要するため、電源投入後に最初の位置情報が出力されるまでの時間は、ウォーム・スタートよりも長くなる。
軌道予測とは、以前収集されたエフェメリス・データを利用し、未来の衛星位置を予測する技術である。軌道予測装置は、この予測した衛星位置を測位に用いることで、ウォーム・スタート時にエフェメリス・データの受信完了を待たずに測位開始することができるため、GPS初期位置算出時間(TTFF)を短縮することができる。GPS受信機が衛星からの情報取得後、電源オフ状態の期間が短ければ、以前取得したデータを用いて、衛星の位置を把握できるため、データの再取得までの時間が短くなる。
そして、軌道予測装置は、エフェメリス・データから計算される衛星位置・速度をバックアップ・メモリに保存する。次に、軌道予測装置は、軌道予測計算に必要なデータが揃い次第、軌道予測パラメータを推定し、推定した軌道予測パラメータを使用して軌道を予測する。次に、軌道予測装置は、予測した軌道をバックアップ・メモリへ保存する。そして、軌道予測装置は、エフェメリス・データが無効の時、バックアップ・メモリに保存された予測した軌道を使用して衛星位置を計算する。
そして、軌道予測装置は、エフェメリス・データから計算される衛星位置・速度をバックアップ・メモリに保存する。次に、軌道予測装置は、軌道予測計算に必要なデータが揃い次第、軌道予測パラメータを推定し、推定した軌道予測パラメータを使用して軌道を予測する。次に、軌道予測装置は、予測した軌道をバックアップ・メモリへ保存する。そして、軌道予測装置は、エフェメリス・データが無効の時、バックアップ・メモリに保存された予測した軌道を使用して衛星位置を計算する。
衛星の位置は、衛星の初期位置と初期速度が与えられ、衛星にかかる加速度が計算され、積分が行われることで、計算される。衛星にかかる加速度ベクトルasatは、地球重力定数をGMe、地球中心から衛星までの位置ベクトルをrとし、地球重力の地球非球面による加速度ベクトルをa1、太陽重力による加速度ベクトルをa2、月重力による加速度ベクトルをa3、太陽輻射圧による加速度ベクトルをa4とすると、次式(1)のように表される。なお、太陽輻射圧とは、太陽から飛んでくる光子(Photon)がGPS衛星の表面に衝突して発生する力である。
一般的な軌道予測プログラムでは、式(1)の加速度asatを数値積分により時間積分することで、GPS衛星の位置と速度を算出する。用いる数値積分法は、例えば、RKF45(ルンゲクッタフェールベルグ法;Runge-Kutta-Fehlberg method)及びCowellである。
次に、本実施形態に係る衛星時計補正パラメータ生成装置を用いた測位装置によるGNSSの構成について、図1を用いて説明する。図1は、本実施形態に係る衛星時計補正パラメータ生成装置を用いた測位装置によるGNSSの概略ブロック図である。
図1に示すように、GNSSは、測位装置1、GNSS衛星2から構成されている。測位装置(GNSS受信機)1は、軌道予測装置10、測位算出部20、衛星時計補正パラメータ生成装置30を備えている。なお、GNSS衛星2は、1機ではなく、GPS衛星だけでも24機以上の衛星から構成される。
軌道予測装置10は、情報取得部110、記憶部115、アルマナック取得部120、衛星型抽出部125、エフェメリス取得部130、衛星位置速度算出部135、衛星位置・速度記憶部140、演算開始判定部145、パラメータ推定部150、軌道予測演算部155、予測結果記憶部160を備えている。
図1に示すように、GNSSは、測位装置1、GNSS衛星2から構成されている。測位装置(GNSS受信機)1は、軌道予測装置10、測位算出部20、衛星時計補正パラメータ生成装置30を備えている。なお、GNSS衛星2は、1機ではなく、GPS衛星だけでも24機以上の衛星から構成される。
軌道予測装置10は、情報取得部110、記憶部115、アルマナック取得部120、衛星型抽出部125、エフェメリス取得部130、衛星位置速度算出部135、衛星位置・速度記憶部140、演算開始判定部145、パラメータ推定部150、軌道予測演算部155、予測結果記憶部160を備えている。
情報取得部110は、GNSS衛星2が送信する航法メッセージから、UTC関係情報、GNSS時刻、電離層遅延補正係数を取得し、取得したGNSS時刻、UTC関係情報と電離層遅延補正係数を記憶部115に記憶させる。
記憶部115には、GNSS時刻、UTC関係情報と電離層遅延補正係数、大気遅延補正係数が記憶されている。
記憶部115には、GNSS時刻、UTC関係情報と電離層遅延補正係数、大気遅延補正係数が記憶されている。
アルマナック取得部120(情報取得部)は、GNSS衛星2が送信したアルマナック・データを取得し、取得したアルマナック・データを衛星型抽出部125に出力する。
衛星型抽出部125は、アルマナック取得部120が出力するアルマナック・データから衛星型番号を抽出し、抽出した衛星型番号を衛星位置・速度記憶部140に記憶させる。なお、衛星型番号とは、Block IIA,Block IIR,Block IIR−M,Block IIFなどのGNSS衛星2の型を示す番号である。
衛星型抽出部125は、アルマナック取得部120が出力するアルマナック・データから衛星型番号を抽出し、抽出した衛星型番号を衛星位置・速度記憶部140に記憶させる。なお、衛星型番号とは、Block IIA,Block IIR,Block IIR−M,Block IIFなどのGNSS衛星2の型を示す番号である。
エフェメリス取得部130(情報取得部)は、GNSS衛星2が送信するエフェメリス・データを取得し、取得したエフェメリス・データを衛星位置速度算出部135に出力する。
衛星位置速度算出部135は、エフェメリス取得部130が出力したエフェメリス・データから衛星の状態(Health)を抽出する。衛星位置速度算出部135は、抽出した衛星の状態が正常であるか否かを判別する。衛星の状態とは、衛星の状態を表すコードで、0以外は何らかの異常が衛星にあることを示している。
このコードが正常を示す0であり、衛星の状態が正常であると判別した場合、衛星位置速度算出部135は、取得したエフェメリス・データをエポック時刻toeに近い偶数正時(例えば、0時、2時、4時、6時などの2で除算して割り切れる時間)TOE2に変換する。衛星位置速度算出部135は、衛星位置・速度記憶部140に記憶されているGNSS衛星2の位置と速度情報(以下、位置速度情報という)の内、算出した時刻TOE2から3日を減算した、時刻TOE2から3日より以前(時刻TOE2の示す日時より3日前)の情報を保存領域から削除する。時刻TOE2より3日以前の情報を保存領域から削除する理由は、エフェメリス・データを均等に3日分記憶させるためである。後述するように、このように均等に3日分、記憶させたエフェメリス・データを用いて軌道予測装置10が軌道予測を行った場合、予測誤差が少なくなるためである。
衛星位置速度算出部135は、エフェメリス取得部130が出力したエフェメリス・データから衛星の状態(Health)を抽出する。衛星位置速度算出部135は、抽出した衛星の状態が正常であるか否かを判別する。衛星の状態とは、衛星の状態を表すコードで、0以外は何らかの異常が衛星にあることを示している。
このコードが正常を示す0であり、衛星の状態が正常であると判別した場合、衛星位置速度算出部135は、取得したエフェメリス・データをエポック時刻toeに近い偶数正時(例えば、0時、2時、4時、6時などの2で除算して割り切れる時間)TOE2に変換する。衛星位置速度算出部135は、衛星位置・速度記憶部140に記憶されているGNSS衛星2の位置と速度情報(以下、位置速度情報という)の内、算出した時刻TOE2から3日を減算した、時刻TOE2から3日より以前(時刻TOE2の示す日時より3日前)の情報を保存領域から削除する。時刻TOE2より3日以前の情報を保存領域から削除する理由は、エフェメリス・データを均等に3日分記憶させるためである。後述するように、このように均等に3日分、記憶させたエフェメリス・データを用いて軌道予測装置10が軌道予測を行った場合、予測誤差が少なくなるためである。
衛星位置速度算出部135は、衛星位置・速度記憶部140に時刻TOE2の位置速度情報が記憶されているか否かを判別する。衛星位置・速度記憶部140に時刻TOE2の位置速度情報が記憶されていると判別した場合、後述するように、衛星位置速度算出部135は、TOE2の位置速度情報を削除後、衛星位置・速度記憶部140の保存領域の空きが、2個分の位置速度情報を記憶できる領域になるまで情報を削除する。衛星位置速度算出部135は、時刻TOE2と時刻TOE2+2時間(時刻TOE2に対して2時間加算した結果)の各時刻の各GPS衛星2の各位置と各速度を算出し、算出した2つの時刻の位置と速度を位置速度情報として衛星位置・速度記憶部140に記憶させる。
衛星位置・速度記憶部140に時刻TOE2の位置速度情報が記憶されていないと判別した場合、後述するように、衛星位置速度算出部135は、衛星位置・速度記憶部140の保存領域の空きが、3個分の位置速度情報を記憶できる領域になるまで情報を削除する。衛星位置速度算出部135は、時刻TOE2、時刻TOE2−2時間(時刻TOE2から2時間減算した結果)、及び時刻TOE2+2時間(時刻TOE2に対して2時間加算した結果)の各時刻のGPS衛星2の各位置と各速度を算出し、算出した2つの時刻の位置と速度を位置速度情報として衛星位置・速度記憶部140に記憶させる。
衛星位置・速度記憶部140に時刻TOE2の位置速度情報が記憶されていないと判別した場合、後述するように、衛星位置速度算出部135は、衛星位置・速度記憶部140の保存領域の空きが、3個分の位置速度情報を記憶できる領域になるまで情報を削除する。衛星位置速度算出部135は、時刻TOE2、時刻TOE2−2時間(時刻TOE2から2時間減算した結果)、及び時刻TOE2+2時間(時刻TOE2に対して2時間加算した結果)の各時刻のGPS衛星2の各位置と各速度を算出し、算出した2つの時刻の位置と速度を位置速度情報として衛星位置・速度記憶部140に記憶させる。
衛星位置・速度記憶部140には、衛星型番号、各GPS衛星2における各時刻の位置速度情報が記憶されている。
演算開始判定部145は、衛星位置・速度記憶部140から過去の位置速度情報、衛星型番号を読み出し、記憶部115からGNSS時刻、UTC関係情報及び電離層遅延補正係数を読み出す。演算開始判定部145は、読み出した情報に基づき、予め設定されている演算開始条件を満たしているか否かを判定する。演算開始条件を満たしていると判別した場合、演算開始判定部145は、パラメータ推定部150に読み出した情報と判定結果を出力する。
演算開始判定部145は、衛星位置・速度記憶部140から過去の位置速度情報、衛星型番号を読み出し、記憶部115からGNSS時刻、UTC関係情報及び電離層遅延補正係数を読み出す。演算開始判定部145は、読み出した情報に基づき、予め設定されている演算開始条件を満たしているか否かを判定する。演算開始条件を満たしていると判別した場合、演算開始判定部145は、パラメータ推定部150に読み出した情報と判定結果を出力する。
パラメータ推定部150は、演算開始判定部145が出力する読み出された情報と判定結果を用いて、推定パラメータsを推定する。パラメータ推定部150は、推定した推定パラメータを軌道予測演算部155に出力する。なお、推定パラメータsは、例えば11個の成分を有するベクトルである。
軌道予測演算部155は、パラメータ推定部150が出力する推定パラメータsを使用し、後述するように、例えば将来の3日間までの2時間ごとの時刻te(k)(k=0,・・・,36)におけるGNSS衛星2_i(i=0,・・・,110)の衛星位置pe(i)(i=0,・・・,110)を軌道計算する。軌道予測演算部155は、算出した予測軌道値を予測結果記憶部160に記憶させる。
予測結果記憶部160には、GNSS衛星2_iの各々の予測軌道値である衛星位置pe(i)が記憶されている。
軌道予測演算部155は、パラメータ推定部150が出力する推定パラメータsを使用し、後述するように、例えば将来の3日間までの2時間ごとの時刻te(k)(k=0,・・・,36)におけるGNSS衛星2_i(i=0,・・・,110)の衛星位置pe(i)(i=0,・・・,110)を軌道計算する。軌道予測演算部155は、算出した予測軌道値を予測結果記憶部160に記憶させる。
予測結果記憶部160には、GNSS衛星2_iの各々の予測軌道値である衛星位置pe(i)が記憶されている。
測位算出部20は、予測結果記憶部160に記憶されている予測軌道値である衛星位置をpe(i)読み出し、読み出した予測軌道値を外挿して衛星位置を算出して測位を行う。すなわち、測位算出部20は、算出された位置速度情報、及び実測により取得された位置速度情報、UTC関係情報、GNSS時刻、衛星時計補正パラメータ、電離層遅延補正係数を用いて、公知の手法で位置測定を行う。ここで、測位算出部20は、航法メッセージを用いてGNSS衛星から測位装置1までの信号到達時間を測定する。測位算出部20は、信号到達時間に光速を掛けることでGNSS衛星と測位装置1間の距離(擬似距離)を求め、複数のGNSS衛星から距離に基づいて測位装置1の位置を三次元上で算出している。ここで、測位算出部20は、以下の(2)式を用いて、GNSS衛星と測位装置1との間の距離を求めている。
ρ(tu)=r(tu,ts)+c{δu(tu)−δs(ts)}
+I(tu)+T(tu)+ep(tu) …(2)
ρ(tu):時刻tuで求めた観測距離
r(tu,ts):tuでの測位装置1の位置とtsでの衛星位置との幾何距離
c:光の速度
δu(tu):測位装置1の時計の誤差
δs(ts):衛星時計の誤差
I(tu):電離層遅延による誤差距離
T(tu):大気遅延による誤差距離
ep(tu):マルチパス誤差
tu:測位装置1が航法メッセージを受信した時間
ts:GNSS衛星2が航法メッセージを送信した時間(GNSS時刻)
ここで、測位算出部20は、記憶部115から電離層遅延補正係数及び大気遅延補正係数を読み出す。測位算出部20は、電離層遅延による誤差距離I(tu)を電離層遅延補正係数により、また大気遅延による誤差距離T(tu)を大気遅延補正係数により、それぞれ予め設定された関係式により算出し、(2)式に用いる。
ρ(tu)=r(tu,ts)+c{δu(tu)−δs(ts)}
+I(tu)+T(tu)+ep(tu) …(2)
ρ(tu):時刻tuで求めた観測距離
r(tu,ts):tuでの測位装置1の位置とtsでの衛星位置との幾何距離
c:光の速度
δu(tu):測位装置1の時計の誤差
δs(ts):衛星時計の誤差
I(tu):電離層遅延による誤差距離
T(tu):大気遅延による誤差距離
ep(tu):マルチパス誤差
tu:測位装置1が航法メッセージを受信した時間
ts:GNSS衛星2が航法メッセージを送信した時間(GNSS時刻)
ここで、測位算出部20は、記憶部115から電離層遅延補正係数及び大気遅延補正係数を読み出す。測位算出部20は、電離層遅延による誤差距離I(tu)を電離層遅延補正係数により、また大気遅延による誤差距離T(tu)を大気遅延補正係数により、それぞれ予め設定された関係式により算出し、(2)式に用いる。
なお、軌道予測装置10は、航法メッセージを受信し、受信した航法メッセージを復調して、復調した航法メッセージから各データを抽出する。情報取得部110とアルマナック取得部120及びエフェメリス取得部130は、この抽出された各データを取得している。例えば、情報取得部110が航法メッセージを復調し、復調した航法メッセージから各データを抽出し、抽出した各情報を衛星型抽出部125及び衛星位置速度算出部135に出力し、記憶部115に書き込むようにしてもよい。また、記憶部115、衛星位置・速度記憶部140、及び予測結果記憶部160に分ける例を説明したが、これらの記憶部をまとめて1つの記憶部に記憶させるようにしてもよい。
軌道予測装置10が行う処理について説明する。まず、軌道予測装置10は、航法データの取得を開始する。次に、アルマナック取得部120は、GNSS衛星2が送信するアルマナック・データを取得し、取得したアルマナック・データを衛星型抽出部125に出力する。次に、衛星型抽出部125は、アルマナック取得部120が出力するアルマナック・データから衛星型番号を抽出し、抽出した衛星型番号を衛星位置・速度記憶部140に記憶させる。次に、情報取得部110は、GNSS衛星2が送信する航法メッセージのサブフレーム4のデータから、電離層遅延補正係数、大気遅延補正係数、大気遅延補正係数、GNSS時刻(ts)、及びUTC関係情報を取得し、取得した電離層遅延補正係数、GNSS時刻、及びUTC関係情報を記憶部115に記憶させる。
衛星時計補正パラメータ生成装置30は、衛星時計補正パラメータ記憶部310と衛星時計変化率算出部320から構成されている。
衛星時計補正パラメータ記憶部310には、衛星時計補正パラメータ(バイアス値Af0及び変化率Af1)が時系列に、所定の期間の複数のデータが衛星時計補正パラメータのエポック時刻とともに書き込まれて記憶されている。例えば、本実施形態の場合、現在から3つ前までの衛星時計補正パラメータが記憶されている。ここで、衛星時計補正パラメータにおけるバイアス値Af0は、エポック時刻における衛星時計の誤差である。変化率Af1は、単位時間あたりの衛星時計の誤差の変化率である。
衛星時計補正パラメータ記憶部310には、衛星時計補正パラメータ(バイアス値Af0及び変化率Af1)が時系列に、所定の期間の複数のデータが衛星時計補正パラメータのエポック時刻とともに書き込まれて記憶されている。例えば、本実施形態の場合、現在から3つ前までの衛星時計補正パラメータが記憶されている。ここで、衛星時計補正パラメータにおけるバイアス値Af0は、エポック時刻における衛星時計の誤差である。変化率Af1は、単位時間あたりの衛星時計の誤差の変化率である。
また、上述したエフェメリス取得部130は、エフェメリス・データを受信した際、衛星時計補正パラメータ記憶部310に記憶されている衛星時計補正パラメータの最も古い時刻のもの消去し、新たに受信したエフェメリス・データにおける衛星時計補正パラメータを、衛星時計補正パラメータ記憶部310に書き込んで記憶させる。すなわち、エフェメリス取得部130は、過去のエフェメリス・データにおける衛星時計補正パラメータ(バイアス値Af0及び変化率Af1)を時系列に所定の範囲で衛星時計補正パラメータ記憶部310に対して書き込んで記憶させる。
衛星時計変化率算出部320は、以下の(3)式により、推定変化率Af1(ts)’の推定を行う。
Af1(ts)’=
{Af0(ts−1)−Af0(ts−2)}/{(ts−1)−(ts−2)} …(3)
上記(3)式において、(ts−1)は、時刻tsの直前(一つ前)に受信したパラメータのエポック時刻である。Af0(ts−1)は、エポック時刻ts−1におけるバイアス値である。時刻(ts−2)は、時刻tsの2つ前に受信したパラメータのエポック時刻である。Af0(ts−2)は、エポック時刻ts−2におけるバイアス値である。すなわち、2つの過去のバイアス値Af0を用いて求めた推定変化率Af1(ts)’は、この2個のバイアス値Af0を結ぶ線分の傾きである。
衛星時計変化率算出部320は、以下の(3)式により、推定変化率Af1(ts)’の推定を行う。
Af1(ts)’=
{Af0(ts−1)−Af0(ts−2)}/{(ts−1)−(ts−2)} …(3)
上記(3)式において、(ts−1)は、時刻tsの直前(一つ前)に受信したパラメータのエポック時刻である。Af0(ts−1)は、エポック時刻ts−1におけるバイアス値である。時刻(ts−2)は、時刻tsの2つ前に受信したパラメータのエポック時刻である。Af0(ts−2)は、エポック時刻ts−2におけるバイアス値である。すなわち、2つの過去のバイアス値Af0を用いて求めた推定変化率Af1(ts)’は、この2個のバイアス値Af0を結ぶ線分の傾きである。
また、測位算出部20は、エフェメリス・データが受信された場合、(2)式で用いる衛星時計の誤差δs(ts)を以下の(4)式から求める。
δs(ts)=Af0(ts)+Af1(ts)・(ts−t0) …(4)
この(4)式において、t0は、衛星時計補正パラメータの元期であり、航法メッセージに含まれている。
一方、測位算出部20は、エフェメリス・データが受信されなかった場合、衛星時計補正パラメータ記憶部310から、前回のバイアス値Af0(ts−1)を読み出し、かつ衛星時計補正パラメータ生成装置30から供給される推定変化率Af1(ts)’とにより、衛星時計の推定誤差時間δ’s(ts)を、以下の(5)式により求める。
δ’s(ts)=Af0(ts−1)+Af1(ts)’・(ts−ts−1) …(5)
測位算出部20は、エフェメリス・データが受信されなかった際において、(5)式から求めた推定誤差時間δ’s(ts)を、エフェメリス・データが受信できた際の誤差時間δs(ts)に変えて、(2)式に代入し、観測距離ρ(tu)を算出して求める。
δs(ts)=Af0(ts)+Af1(ts)・(ts−t0) …(4)
この(4)式において、t0は、衛星時計補正パラメータの元期であり、航法メッセージに含まれている。
一方、測位算出部20は、エフェメリス・データが受信されなかった場合、衛星時計補正パラメータ記憶部310から、前回のバイアス値Af0(ts−1)を読み出し、かつ衛星時計補正パラメータ生成装置30から供給される推定変化率Af1(ts)’とにより、衛星時計の推定誤差時間δ’s(ts)を、以下の(5)式により求める。
δ’s(ts)=Af0(ts−1)+Af1(ts)’・(ts−ts−1) …(5)
測位算出部20は、エフェメリス・データが受信されなかった際において、(5)式から求めた推定誤差時間δ’s(ts)を、エフェメリス・データが受信できた際の誤差時間δs(ts)に変えて、(2)式に代入し、観測距離ρ(tu)を算出して求める。
図2は、指定誤差時間δ’s(ts)の算出を説明する図である。縦軸がδsを示し、横軸がGNSS時刻を示している。ここで、衛星時計補正パラメータ記憶部310には、衛星時計補正パラメータ(Af0(ts−3),Af1(ts−3))、(Af0(ts−2),Af1(ts−2))、(Af0(ts),Af1(ts))が記憶されている。
エフェメリス取得部130は、GNSS時刻tsにGNSS衛星2が送信したエフェメリス・データが正常に受信されなかった場合、衛星時計補正パラメータ生成装置30に対して、エフェメリス・データが得られなかったことを示す異常信号を出力する。
これにより、衛星時計変化率算出部320は、エポック時刻がGNSS時刻ts−1及びts−2である衛星時計補正パラメータ(Af0(ts−3),Af1(ts−3))、(Af0(ts−2),Af1(ts−2))を衛星時計補正パラメータ記憶部310から読み出し、(3)式に代入して、推定変化率Af1(ts)’を算出し、測位算出部20に対して出力する。また、衛星時計変化率算出部320は、上記推定変化率Af1(ts)’を求める際、2つの時刻における衛星時計補正パラメータではなく、3つ以上の複数の時刻の衛星時計補正パラメータのバイアス値Af0を用いて、例えば各時刻のバイアス値との二乗誤差が最小となる線分の傾きを変化率として求めるように構成しても良い。
エフェメリス取得部130は、GNSS時刻tsにGNSS衛星2が送信したエフェメリス・データが正常に受信されなかった場合、衛星時計補正パラメータ生成装置30に対して、エフェメリス・データが得られなかったことを示す異常信号を出力する。
これにより、衛星時計変化率算出部320は、エポック時刻がGNSS時刻ts−1及びts−2である衛星時計補正パラメータ(Af0(ts−3),Af1(ts−3))、(Af0(ts−2),Af1(ts−2))を衛星時計補正パラメータ記憶部310から読み出し、(3)式に代入して、推定変化率Af1(ts)’を算出し、測位算出部20に対して出力する。また、衛星時計変化率算出部320は、上記推定変化率Af1(ts)’を求める際、2つの時刻における衛星時計補正パラメータではなく、3つ以上の複数の時刻の衛星時計補正パラメータのバイアス値Af0を用いて、例えば各時刻のバイアス値との二乗誤差が最小となる線分の傾きを変化率として求めるように構成しても良い。
測位算出部20は、推定変化率Af1(ts)’が供給されると、衛星時計補正パラメータ記憶部310から直前の衛星時計補正パラメータのバイアス値Af0(ts−1)を読み出す。そして、測位算出部20は、衛星時計補正パラメータ生成装置30から供給された推定変化率Af1(ts)’と、直前に受信された衛星時計補正パラメータのバイアス値Af0(ts−1)とを(5)式に代入して推定誤差時間δ’s(ts)の算出を行う。
上述したように本実施形態によれば、過去に正常に受信された複数の衛星時計補正パラメータのバイアス値Af0の変化より、推定変化率Af1’を推定するため、従来のように直前の推定値を代用する従来に比較し、より衛星時計の誤差の変化を反映することができるため、GNSS衛星2の観測距離ρ(tu)をより高い精度で求めることが可能となる。
また、本実施形態によれば、GNSS衛星2の観測距離ρ(tu)を、従来に比較してより高い精度で求められるため、従来に比較して、測位装置1の測位の結果も高い精度で得ることが可能となる。
また、本実施形態によれば、GNSS衛星2の観測距離ρ(tu)を、従来に比較してより高い精度で求められるため、従来に比較して、測位装置1の測位の結果も高い精度で得ることが可能となる。
また、図1における測位装置1の機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより、エフェメリス・データが得られない際における衛星時計補正パラメータの変化率Af1を求め、測位装置1の位置を測位する測位処理を行ってもよい。なお、ここでいう「コンピュータシステム」とは、OSや周辺機器等のハードウェアを含むものとする。
また、「コンピュータシステム」は、WWWシステムを利用している場合であれば、ホームページ提供環境(あるいは表示環境)も含むものとする。
また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ROM、CD−ROM等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含むものとする。また上記プログラムは、前述した機能の一部を実現するためのものであっても良く、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであっても良い。
また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ROM、CD−ROM等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含むものとする。また上記プログラムは、前述した機能の一部を実現するためのものであっても良く、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであっても良い。
以上、この発明の実施形態を図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。
1…測位装置 2…GNSS衛星 10…軌道予測装置 20…測位算出部 30…衛星時計補正パラメータ生成装置 110…情報取得部 115…記憶部 120…アルマナック取得部 125…衛星型抽出部 130…エフェメリス取得部 135…衛星位置速度算出部 140…衛星位置・速度記憶部 145…演算開始判定部 150…パラメータ推定部 155…軌道予測演算部 160…予測結果記憶部 310…衛星時計補正パラメータ記憶部 320…衛星時計変化率算出部
Claims (4)
- GNSSにおける測位装置に設けられ、GNSS衛星から衛星時計補正パラメータが正常に受信できなかった場合に、測位に用いる前記GNSS衛星と前記測位装置との観測距離を求める際に用いる前記GNSS衛星の衛星時計の誤差を求めるための衛星時計補正パラメータを推定する衛星時計補正パラメータ生成装置であり、
過去の衛星時計補正パラメータであるバイアス値及び誤差値の変化率が書き込まれて記憶されている衛星時計補正パラメータ記憶部と、
過去の複数の前記バイアス値から求めた傾きを前記変化率の推定値である推定変化率として出力する衛星時計変化率算出部と
を備えることを特徴とする衛星時計補正パラメータ生成装置。 - 前記バイアス値を2つ用いる場合、
前記衛星時計変化率算出部が、前記GNSS衛星から衛星時計補正パラメータを正常に受信できなかった場合に、直前とその前のバイアス値を結ぶ線分の傾き前記推定変化率とすることを特徴とする請求項1に記載の衛星時計補正パラメータ生成装置。 - 請求項1または請求項2に記載の衛星時計補正パラメータ生成装置を備える測位装置であり、
GNSS衛星から衛星時計補正パラメータが正常に受信できなかった場合に前記観測距離を求める際、前記衛星時計の時刻誤差を、直前に正常に受信した前記バイアス値と、前記衛星時計補正パラメータ生成装置が推定した前記推定変化率とを用いて算出する
ことを特徴とする測位装置。 - GNSSにおける測位装置において、GNSS衛星から衛星時計補正パラメータが正常に受信できなかった場合に、測位に用いる前記GNSS衛星と前記測位装置との観測距離を求める際に用いる前記GNSS衛星の衛星時計の誤差を求めるための衛星時計補正パラメータを推定する衛星時計補正パラメータ生成方法であり、
衛星時計変化率算出部が、過去の衛星時計補正パラメータであるバイアス値及び誤差値の変化率が書き込まれて記憶されている衛星時計補正パラメータ記憶部から、過去の複数の前記バイアス値から求めた傾きを前記変化率の推定値である推定変化率として出力する衛星時計変化率算出過程
を有することを特徴とする衛星時計補正パラメータ生成方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013151945A JP2015021900A (ja) | 2013-07-22 | 2013-07-22 | 衛星時計補正パラメータ生成装置、その方法及び測位装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013151945A JP2015021900A (ja) | 2013-07-22 | 2013-07-22 | 衛星時計補正パラメータ生成装置、その方法及び測位装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2015021900A true JP2015021900A (ja) | 2015-02-02 |
Family
ID=52486462
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2013151945A Pending JP2015021900A (ja) | 2013-07-22 | 2013-07-22 | 衛星時計補正パラメータ生成装置、その方法及び測位装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2015021900A (ja) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105158780A (zh) * | 2015-07-24 | 2015-12-16 | 北京跟踪与通信技术研究所 | 一种基于多种导航卫星可互换的导航定位方法 |
US20160187488A1 (en) * | 2014-12-30 | 2016-06-30 | Mediatek Inc. | Satellite-based positioning method and associated apparatus |
CN108229746A (zh) * | 2018-01-08 | 2018-06-29 | 唐山学院 | 卫星钟差数据的多项式预报模型的最优化改正方法 |
CN112051597A (zh) * | 2020-09-27 | 2020-12-08 | 西安领航遥感技术有限公司 | 一种精密单点定位方法与装置 |
CN112782734A (zh) * | 2021-02-05 | 2021-05-11 | 中国空间技术研究院 | 一种导航星座自主导航电文系统差检测与处理方法 |
CN113810146A (zh) * | 2021-08-05 | 2021-12-17 | 西安理工大学 | 一种无外部参考时钟的数据同步传输方法 |
JP2023521983A (ja) * | 2020-04-09 | 2023-05-26 | テレフオンアクチーボラゲット エルエム エリクソン(パブル) | 非地上ネットワークシステム情報のgnssデータ |
CN116931026A (zh) * | 2023-09-15 | 2023-10-24 | 北京神导科技股份有限公司 | 一种卫星导航信号的异常判定方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0627216A (ja) * | 1992-07-08 | 1994-02-04 | Japan Radio Co Ltd | Gps受信機 |
US20080238765A1 (en) * | 2007-03-30 | 2008-10-02 | Gengsheng Zhang | Satellite Clock Prediction |
JP2010156631A (ja) * | 2008-12-29 | 2010-07-15 | Seiko Epson Corp | 衛星軌道データの圧縮方法、衛星軌道データの提供方法、衛星軌道データの展開方法及びデータ圧縮装置 |
JP2012013505A (ja) * | 2010-06-30 | 2012-01-19 | Japan Radio Co Ltd | 測位支援装置 |
JP2012502298A (ja) * | 2008-09-11 | 2012-01-26 | カリフォルニア インスティチュート オブ テクノロジー | Gnss暦の自律的な受信機内予測方法及び装置 |
-
2013
- 2013-07-22 JP JP2013151945A patent/JP2015021900A/ja active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0627216A (ja) * | 1992-07-08 | 1994-02-04 | Japan Radio Co Ltd | Gps受信機 |
US20080238765A1 (en) * | 2007-03-30 | 2008-10-02 | Gengsheng Zhang | Satellite Clock Prediction |
JP2012502298A (ja) * | 2008-09-11 | 2012-01-26 | カリフォルニア インスティチュート オブ テクノロジー | Gnss暦の自律的な受信機内予測方法及び装置 |
JP2010156631A (ja) * | 2008-12-29 | 2010-07-15 | Seiko Epson Corp | 衛星軌道データの圧縮方法、衛星軌道データの提供方法、衛星軌道データの展開方法及びデータ圧縮装置 |
JP2012013505A (ja) * | 2010-06-30 | 2012-01-19 | Japan Radio Co Ltd | 測位支援装置 |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160187488A1 (en) * | 2014-12-30 | 2016-06-30 | Mediatek Inc. | Satellite-based positioning method and associated apparatus |
US10185035B2 (en) * | 2014-12-30 | 2019-01-22 | Mediatek Inc. | Satellite-based positioning method and associated apparatus |
CN105158780A (zh) * | 2015-07-24 | 2015-12-16 | 北京跟踪与通信技术研究所 | 一种基于多种导航卫星可互换的导航定位方法 |
CN108229746A (zh) * | 2018-01-08 | 2018-06-29 | 唐山学院 | 卫星钟差数据的多项式预报模型的最优化改正方法 |
JP2023521983A (ja) * | 2020-04-09 | 2023-05-26 | テレフオンアクチーボラゲット エルエム エリクソン(パブル) | 非地上ネットワークシステム情報のgnssデータ |
JP7549675B2 (ja) | 2020-04-09 | 2024-09-11 | テレフオンアクチーボラゲット エルエム エリクソン(パブル) | 無線デバイスによって実行される方法、ネットワークノードによって実行される方法、無線デバイス、及びネットワークノード |
CN112051597A (zh) * | 2020-09-27 | 2020-12-08 | 西安领航遥感技术有限公司 | 一种精密单点定位方法与装置 |
CN112782734B (zh) * | 2021-02-05 | 2024-05-31 | 中国空间技术研究院 | 一种导航星座自主导航电文系统差检测与处理方法 |
CN112782734A (zh) * | 2021-02-05 | 2021-05-11 | 中国空间技术研究院 | 一种导航星座自主导航电文系统差检测与处理方法 |
CN113810146A (zh) * | 2021-08-05 | 2021-12-17 | 西安理工大学 | 一种无外部参考时钟的数据同步传输方法 |
CN113810146B (zh) * | 2021-08-05 | 2023-09-22 | 西安理工大学 | 一种无外部参考时钟的数据同步传输方法 |
CN116931026B (zh) * | 2023-09-15 | 2023-12-12 | 北京神导科技股份有限公司 | 一种卫星导航信号的异常判定方法 |
CN116931026A (zh) * | 2023-09-15 | 2023-10-24 | 北京神导科技股份有限公司 | 一种卫星导航信号的异常判定方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6191074B2 (ja) | 衛星時計精度判定装置、その方法及び測位装置 | |
JP6426222B2 (ja) | エフェメリス拡張システムとgnssでの使用方法 | |
JP2015021900A (ja) | 衛星時計補正パラメータ生成装置、その方法及び測位装置 | |
KR100958090B1 (ko) | 향상된 자율 위치 확인 시스템을 위한 방법 및 장치 | |
CN100397094C (zh) | 卫星定位系统接收机中的时间确定及其方法 | |
EP2309290B1 (en) | Satellite positioning system receivers and methods | |
KR101000984B1 (ko) | 위성 위치 확인 시스템에서 거리 변화율 측정값을 이용하여위치를 검증하는 방법 및 장치 | |
EP2056122B1 (en) | Method and apparatus for backround decoding of a satellite navigation message to maintain integrity of long term orbit information in a remote receiver | |
US20070299609A1 (en) | Method and system for ephemeris extension for GNSS applications | |
US7355551B2 (en) | Method for a space-efficient GPS receiver | |
US20080111738A1 (en) | Method and apparatus in standalone positioning without broadcast ephemeris | |
US20100198512A1 (en) | Method and apparatus for providing reliable extended ephemeris quality indicators | |
WO2001014903A1 (en) | Positioning receiver and related method | |
JP4723932B2 (ja) | 測位システム | |
JP2011080928A (ja) | 測位システム | |
JP4322829B2 (ja) | サイクルスリップ検出装置及びサイクルスリップ検出方法 | |
JP2013044576A (ja) | 人工衛星ベースの受信機に対するsbas支援情報を持ったagpsサーバ | |
KR20100093346A (ko) | 위치 결정 시스템 수신기 및 그것의 위성 시간 보정 방법 | |
US8487811B1 (en) | Sync feedback for time to first fix | |
JP5801111B2 (ja) | 軌道予測装置、及び軌道予測方法 | |
US20150070213A1 (en) | Location determination in multi-system gnns environment using conversion of data into a unified format | |
Zhang et al. | A Novel Method to Provide URE in Broadcast or Extended Ephemeris | |
JP2013003092A (ja) | 軌道予測装置、及び軌道予測方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20160707 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20170425 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20170509 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20171107 |