JP4988028B2 - Gps受信機 - Google Patents
Gps受信機 Download PDFInfo
- Publication number
- JP4988028B2 JP4988028B2 JP2010258962A JP2010258962A JP4988028B2 JP 4988028 B2 JP4988028 B2 JP 4988028B2 JP 2010258962 A JP2010258962 A JP 2010258962A JP 2010258962 A JP2010258962 A JP 2010258962A JP 4988028 B2 JP4988028 B2 JP 4988028B2
- Authority
- JP
- Japan
- Prior art keywords
- pseudorange
- pseudo distance
- pseudo
- error
- gps
- 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.)
- Active
Links
- 238000009499 grossing Methods 0.000 claims description 108
- 238000012545 processing Methods 0.000 claims description 79
- 230000005855 radiation Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 description 78
- 238000012937 correction Methods 0.000 description 18
- 230000000694 effects Effects 0.000 description 18
- 238000000034 method Methods 0.000 description 15
- 239000005433 ionosphere Substances 0.000 description 13
- 238000011156 evaluation Methods 0.000 description 12
- 239000005436 troposphere Substances 0.000 description 9
- 230000005856 abnormality Effects 0.000 description 5
- 230000007423 decrease Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 102100036848 C-C motif chemokine 20 Human genes 0.000 description 2
- 101000713099 Homo sapiens C-C motif chemokine 20 Proteins 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 102100035353 Cyclin-dependent kinase 2-associated protein 1 Human genes 0.000 description 1
- 101000661816 Homo sapiens Suppression of tumorigenicity 18 protein Proteins 0.000 description 1
- 206010027339 Menstruation irregular Diseases 0.000 description 1
- 102100029860 Suppressor of tumorigenicity 20 protein Human genes 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000000593 degrading effect Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Description
三角測量を行うには、3つ以上のGPS衛星の位置と、個々のGPS衛星とGPS受信機間の距離(以下、擬似距離と呼ぶ)とが必要になる。
GPS衛星の位置については、GPS受信機がGPS衛星から放送される航法メッセージを受信すれば、その航法メッセージに含まれている軌道情報から時々刻々の位置を把握することができる。
一方、擬似距離については、GPS衛星から放射された電波が地上に届く伝播時間τから算出することができる。
ただし、GPS受信機が伝播時間τを算出するには、GPS受信機の内蔵時計をGPSの時系である“GPS Time”に同期させる必要がある。
なお、測位計算に使用するGPS衛星の数が3つであれば、GPS受信機の3次元位置を計測することができる。また、GPS衛星の数が4つ以上になれば、さらにGPS受信機の内蔵時計を“GPS Time”に同期させることができる。
これらの誤差が未補正な状態で擬似距離ρCτを測位計算に使用すると、測位結果である位置に誤差を生じるほか、GPS受信機の内蔵時計を“GPS Time”に上手く同期させることができなくなるので、電波伝搬遅延誤差とGPS衛星搭載時計誤差を除去する必要がある。
ρ=ρCτ+dTSATC−δρion−δρtro−δρmp [m] (1)
ただし、Cは光速である。
マルチパスなどの擬似距離誤差を低減して、測位精度を高めているGPS受信機が以下の特許文献1〜3に開示されている。
即ち、特許文献1では、GPS衛星と自動車の相対運動によるドップラー偏位周波数fdop-estを算出する。式(2)において、V0はGPS衛星の進行方向に対する速度[m/s]、θ1はGPS衛星と自動車とを結ぶ線と、GPS衛星の進行方向とがなす角度[rad]である。
fdop-est=fL1V0cosθ1/(C−V0cosθ1) (2)
そして、GPS受信機により計測されたドップラー偏位周波数fdopと、GPS衛星と自動車の相対運動によるドップラー偏位周波数fdop-estとの差異が所定値以上になると、そのGPS衛星の擬似距離ρCτに大きなマルチパス誤差が含まれていると判断し、この擬似距離ρCτを測位計算に使用しないようにしている。
即ち、特許文献2では、下記の式(3)を用いて、GPS衛星と自動車の相対運動によるドップラー偏位周波数fdop-estを算出する。
fdop-est={LOSx(VSxi−VCxi)+LOSy(VSyi−VCyi)
+LOSz(VSzi−VCzi)}fL1/C [Hz] (3)
LOSx=(PSxi−PCxi)/R
LOSy=(PSyi−PCyi)/R
LOSz=(PSzi−PCzi)/R
R={(PSxi−PCxi)2+(PSyi−PCyi)2+(PSzi−PCzi)2}1/2
ただし、(PSxi,PSyi,PSzi)はGPS衛星の位置のxyz成分、(VSxi,VSyi,VSzi)はGPS衛星の速度のxyz成分、(PCxi,PCyi,PCzi)は自動車の位置のxyz成分、(VCxi,VCyi,VCzi)は自動車の速度のxyz成分である。
LOSは地上の自動車から上空のGPS衛星を見る視線方向ベクトル(line−of−site vector;LOS)である。
そして、下記の式(6)のように、擬似距離ρCτiから擬似距離誤差δρiを除去し、擬似距離誤差除去後の擬似距離ρiを測位計算に使用するようにする。
ρi=ρCτi−δρi (6)
具体的には、下記の通りである。
まず、キャリアスムージングは、下記の式(7)に示すように、搬送波位相(ドップラー偏位周波数fdopi)の積算値φiの差分(φi−φi-1)を用いて、電波伝播時間から計測された擬似距離ρCτiを平滑化(重み付け平均)する方法であり、このキャリアスムージングは、擬似距離の誤差が周期的に変動する場合に有効である。
ρSi=ρCτi/N+(N−1)/N{ρSi-1+(φi−φi-1)λΔT} [m]
(7)
ただし、Nはスムージングのサンプル数であり、通常100に設定される。λは搬送波の波長[m]、ΔTは測位インターバル時間[s]である。
図1はこの発明の実施の形態1によるGPS受信機を示す構成図であり、図において、GPSアンテナ1は例えば自動車に取り付けられるGPS受信機のアンテナであり、例えば、自動車の上空に位置する複数個のGPS衛星から放射された電波を受信する。なお、GPSアンテナ1は電波受信手段を構成している。
信号処理部2はGPSアンテナ1により受信された電波に含まれているGPS衛星のコードを参照して電波の放射元であるGPS衛星を識別し、GPS衛星毎にGPS衛星の擬似距離と擬似距離変化率を計測する処理を実施する。また、GPSアンテナ1により受信された電波からGPS衛星の軌道情報などが記述されている航法メッセージを抽出し、その航法メッセージや電波の受信時刻を出力する処理を実施する。なお、信号処理部2は信号処理手段を構成している。
擬似距離誤差算出部4は信号処理部2から出力された電波の受信時刻と航法メッセージ中の補正パラメータ(電離層通過時の擬似距離誤差の補正パラメータと、GPS衛星に搭載されている時計の誤差の補正パラメータ)などを用いて、各GPS衛星から放射された電波が大気中(電離層、対流圏)を通過する際に生じる擬似距離誤差と衛星搭載時計誤差を算出する処理を実施する。なお、擬似距離誤差算出部4は擬似距離誤差算出手段を構成している。
擬似距離誤差評価部6は信号処理部2により計測された擬似距離と擬似距離平滑化部5により平滑化された擬似距離の誤差分散を算出する処理を実施する。なお、擬似距離誤差評価部6は誤差分散算出手段を構成している。
図2はこの発明の実施の形態1によるGPS受信機の処理内容を示すフローチャートある。
信号処理部2は、GPSアンテナ1が自動車の上空に位置するGPS衛星から放射された電波を受信すると、例えば所定の測位周期(1Hz)毎に、GPSアンテナ1により受信された電波に含まれているGPS衛星のコードを参照して電波の放射元であるGPS衛星を識別する(図3、図4を参照)。
信号処理部2は、電波の放射元であるGPS衛星を識別すると、そのGPS衛星の擬似距離ρCτiと擬似距離変化率Δρrateiを計測し、その擬似距離ρCτiと擬似距離変化率Δρrateiを擬似距離平滑化部5及び擬似距離誤差評価部6に出力する。
また、信号処理部2は、GPSアンテナ1により受信された電波からGPS衛星の軌道情報などが記述されている航法メッセージを抽出し、その航法メッセージや電波の受信時刻を衛星位置算出部3及び擬似距離誤差算出部4に出力する。
擬似距離誤差算出部4は、信号処理部2から電波の受信時刻と航法メッセージを受けると、電波の受信時刻と航法メッセージ中の補正パラメータ(電離層通過時の擬似距離誤差の補正パラメータ)などを用いて、各GPS衛星から放射された電波が大気中(電離層、対流圏)を通過する際に生じる擬似距離誤差を算出する(ステップST2)。
また、擬似距離誤差算出部4は、電波の受信時刻と航法メッセージ中の補正パラメータ(GPS衛星に搭載されている時計誤差の補正パラメータ)などを用いて、衛星搭載時計誤差を算出する(ステップST2)。
Δρratei=fdopiC/fL1 [m/s] (8)
ΔρCτi=(ρCτi−ρCτi-1)/Δt (9)
k=1/|ΔρCτi−Δρratei| (10)
ただし、fL1は搬送波周波数、fdopはドップラー偏位周波数、Cは光速、Δtは測位インターバル時間(例えば、1秒)である。
ρSi=ρCτik+{ρSi-1+ΔρrateiΔt}(1−k) (11)
なお、擬似距離ρCτiと擬似距離変化率Δρrateiの双方に誤差が含まれていない場合、擬似距離の差分ΔρCτiと擬似距離変化率Δρrateiは同じ値になる。
ΔδρCτi=|ρCτi−ρSi| (12)
補正後のρSi=ρSi+dTSATC−δρion−δρtro [m] (13)
ただし、(PCxi,PCyi,PCzi)は自動車の位置のxyz成分、(PSxi,PSyi,PSzi)はGPS衛星の位置のxyz成分、Aは航法行列、Wは擬似距離誤差の重み付き行列である。
一方、収束していなければ、ステップST1〜ST6の処理を繰り返し実施させる制御信号を衛星位置算出部3や擬似距離平滑化部5などに出力する。
なお、位置算出部7は、電波を受信できているGPS衛星の個数が4衛星未満、または、収束計算の回数が所定回数以上であれば、自動車の位置を測位することができないと判断して、測位周期毎の動作を終了する。
また、電波伝搬時間τから計測した擬似距離の誤差が小さい場合、この擬似距離を多く使うような重み係数kとなるため、擬似距離変化率Δρrateiを不必要に多用しなくなり、安価なGPS受信機を使用しても、精度の高い測位結果を得ることができる。
上記実施の形態1では、信号処理部2により計測された擬似距離ρCτiと擬似距離変化率Δρrateiから重み係数kを求め、その重み係数kに応じた擬似距離ρCτiと擬似距離変化率Δρrateiの重み加算を実施して、その擬似距離ρCτiを平滑化するものについて示したが、信号処理部2により前回計測された擬似距離ρCτi-1(第二の擬似距離)と今回計測された擬似距離ρCτi(第一の擬似距離)の差分と、信号処理部2により今回計測された擬似距離変化率Δρratei(第一の擬似距離変化率)との差異が所定値より大きい場合、今回計測された擬似距離ρCτiを使用しないで、擬似距離の平滑化を実施するようにしてもよい。
図5はこの発明の実施の形態2によるGPS受信機の処理内容を示すフローチャートある。
なお、この実施の形態2では、擬似距離平滑化部5の処理内容のみが、上記実施の形態1と異なるので、擬似距離平滑化部5の処理内容のみを説明する。
即ち、擬似距離平滑化部5は、GPSアンテナ1が所定のGPS衛星から放射された電波を受信しているか否かを判定する(ステップST11)。
擬似距離平滑化部5は、所定のGPS衛星から放射された電波を受信していれば、さらに、連続受信であるか否かを判定する(ステップST12)。即ち、電波を途切れ途切れに受信している状況であるか否かを判定する。
即ち、擬似距離平滑化部5は、信号処理部2により前回計測された擬似距離ρCτi-1と今回計測された擬似距離ρCτiの差分と、信号処理部2により今回計測された擬似距離変化率Δρrateiとの差異が所定値より大きい場合、擬似距離の誤差が過大であると判断する。
擬似距離平滑化部5は、擬似距離の誤差が過大であると判断すると、下記の式(15)に示すように、信号処理部2により今回計測された擬似距離ρCτiを使用しないで、擬似距離の平滑化を実施する(ステップST14)。
ρSi=ρSi-1+ΔρrateiΔt (15)
擬似距離平滑化部5は、擬似距離の誤差が過大でないと判断すると、上記実施の形態1と同様に、式(11)を用いて、擬似距離の平滑化を実施する(ステップST15)。
擬似距離平滑化部5は、非受信時間tが所定時間に満たない場合、下記の式(16)に示すように、電波が途絶える前の擬似距離変化率Δρratei-tを用いて擬似距離の平滑化を再開する(ステップST17)。
ρSi=ρSi-t+Δρratei-tt (16)
擬似距離平滑化部5は、非受信時間tが所定時間を超えた後、同じGPS衛星の電波を再受信した場合、下記の式(17)に示すように、擬似距離の平滑値を初期化する(ステップST18)。
ρSi=ρCτi (17)
擬似距離平滑化部5は、所定のGPS衛星から放射された電波を受信していなければ、所定のGPS衛星の電波が連続して非受信となっている時間tを計測し、測位周期毎の擬似距離の平滑化処理を終了する(ステップST20)。
上記実施の形態1では、信号処理部2により計測された擬似距離ρCτiと擬似距離変化率Δρrateiから重み係数kを求め、その重み係数kに応じた擬似距離ρCτiと擬似距離変化率Δρrateiの重み加算を実施して、その擬似距離ρCτiを平滑化するものについて示したが、GPSアンテナ1により受信されているGPS衛星の電波が途絶えても、その電波が途絶えてから所定時間を経過するまでは、電波が途絶える前の擬似距離変化率を用いて擬似距離の平滑化を継続するようにしてもよい。
図6はこの発明の参考例1によるGPS受信機の処理内容を示すフローチャートある。
なお、この参考例1では、擬似距離平滑化部5の処理内容のみが、上記実施の形態1と異なるので、擬似距離平滑化部5の処理内容のみを説明する。
即ち、擬似距離平滑化部5は、GPSアンテナ1が所定のGPS衛星から放射された電波を受信しているか否かを判定する(ステップST21)。
擬似距離平滑化部5は、所定のGPS衛星から放射された電波を受信していれば、所定のGPS衛星の電波が連続して非受信となっている時間tが所定時間に到達しているか否かを判定する(ステップST22)。
即ち、擬似距離平滑化部5は、信号処理部2により前回計測された擬似距離ρCτi-1と今回計測された擬似距離ρCτiの差分と、信号処理部2により今回計測された擬似距離変化率Δρrateiとの差異が所定値より大きい場合、擬似距離の誤差が過大であると判断する。
擬似距離平滑化部5は、擬似距離の誤差が過大であると判断すると、上記実施の形態2と同様に、式(15)を用いて、信号処理部2により今回計測された擬似距離ρCτiを使用しないで、擬似距離の平滑化を実施する(ステップST24)。
擬似距離平滑化部5は、擬似距離の誤差が過大でないと判断すると、上記実施の形態1と同様に、式(11)を用いて、擬似距離の平滑化を実施する(ステップST25)。
擬似距離平滑化部5は、所定のGPS衛星の電波が連続して非受信となっている時間tが所定時間を超えていれば、上記実施の形態2と同様に、式(17)を用いて、擬似距離の平滑値を初期化する(ステップST26)。
擬似距離平滑化部5は、上記のようにして、擬似距離を平滑化すると、所定のGPS衛星の電波が連続して非受信となっている時間tをゼロにして、測位周期毎の擬似距離の平滑化処理を終了する(ステップST27)。
そして、擬似距離平滑化部5は、所定のGPS衛星の電波が連続して非受信となっている時間tが所定時間に到達しているか否かを判定する(ステップST29)。
擬似距離平滑化部5は、所定のGPS衛星の電波が連続して非受信となっている時間tが所定時間に到達していなければ、下記の式(18)に示すように、電波が途絶える前の擬似距離変化率Δρratei-tを用いて擬似距離の平滑化を継続する(ステップST30)。
ρSi=ρSi-1+Δρratei-tΔt (18)
なお、自動車の上空のGPS衛星数が4衛星未満の場合に、電波が途絶える前の擬似距離変化率を用いて擬似距離の平滑化を継続するようにすれば、電波受信が途切れがちとなり、測位が間欠的になっていた場所(例えば、市街地、橋の下)でも、自動車の位置を連続的に測位して、アベイラビリティを高めることができる効果を奏する。
図7はこの発明の参考例2によるGPS受信機を示す構成図であり、図において、図1と同一符号は同一または相当部分を示すので説明を省略する。
衛星位置・速度算出部11は信号処理部2から出力された電波の受信時刻と航法メッセージ中の軌道情報(エフェメリスなど)を用いて、各GPS衛星の時々刻々の位置と速度を算出する処理を実施する。なお、衛星位置・速度算出部11は衛星位置特定手段を構成している。
擬似距離変化推定部12はGPS衛星に対する相対運動のドップラー効果を考慮して、衛星位置・速度算出部11により算出されたGPS衛星の位置及び速度と、速度算出部13により算出された自動車の速度と、位置算出部7により算出された自動車の位置とから擬似距離変化率を推定する処理を実施する。
なお、衛星位置・速度算出部11、擬似距離変化推定部12及び速度算出部13から擬似距離変化率推定手段が構成されている。
擬似距離平滑化部14は信号処理部2により計測された擬似距離と擬似距離変化推定部12により推定された擬似距離変化率から重み係数を求め、その重み係数に応じた擬似距離と擬似距離変化率の重み加算を実施して、各GPS衛星の擬似距離を平滑化する処理を実施する。擬似距離平滑化部14は平滑化手段を構成している。
図8はこの発明の参考例2によるGPS受信機の処理内容を示すフローチャートある。
信号処理部2は、GPSアンテナ1が自動車の上空に位置するGPS衛星から放射された電波を受信すると、例えば所定の測位周期(1Hz)毎に、GPSアンテナ1により受信された電波に含まれているGPS衛星のコードを参照して電波の放射元であるGPS衛星を識別する(図3、図4を参照)。
信号処理部2は、電波の放射元であるGPS衛星を識別すると、そのGPS衛星の擬似距離ρCτiと擬似距離変化率Δρrateiを計測し、その擬似距離ρCτiを擬似距離平滑化部14及び擬似距離誤差評価部6に出力する。また、その擬似距離変化率Δρrateiを速度算出部13に出力する。
また、信号処理部2は、上記実施の形態1と同様に、GPSアンテナ1により受信された電波からGPS衛星の軌道情報などが記述されている航法メッセージを抽出し、その航法メッセージや電波の受信時刻を衛星位置算出部3及び擬似距離誤差算出部4に出力する。
擬似距離誤差算出部4は、信号処理部2から電波の受信時刻と航法メッセージを受けると、上記実施の形態1と同様に、電波の受信時刻と航法メッセージ中の補正パラメータ(電離層通過時の擬似距離誤差の補正パラメータ)などを用いて、各GPS衛星から放射された電波が大気中(電離層、対流圏)を通過する際に生じる擬似距離誤差を算出する(ステップST32)。
また、擬似距離誤差算出部4は、電波の受信時刻と航法メッセージ中の補正パラメータ(GPS衛星に搭載されている時計の誤差の補正パラメータ)などを用いて、衛星搭載時計誤差を算出する(ステップST32)。
Δρratei-mi={LOSxi(VSxi−VCxi)+LOSyi(VSyi−VCyi)
+LOSzi(VSzi−VCzi)} [m/s] (21)
LOSxi=(PSxi−PCxi)/R
LOSyi=(PSyi−PCyi)/R
LOSzi=(PSzi−PCzi)/R
R={(PSxi−PCxi)2+(PSyi−PCyi)2+(PSzi−PCzi)2}1/2
Δρratei-mi=fdopiC/fL1 [m/s] (22)
ΔρCτi=(ρCτi−ρCτi-1)/Δt (23)
k=1/|ΔρCτi−Δρratei-mi| (24)
ただし、fL1は搬送波周波数、fdopはドップラー偏位周波数、Cは光速、Δtが測位インターバル時間(例えば、1秒)である。
ρSi=ρCτik+{ρSi-1+Δρratei-miΔt}(1−k) (25)
なお、擬似距離ρCτiと擬似距離変化率Δρratei-miの双方に誤差が含まれていない場合、擬似距離の差分ΔρCτiと擬似距離変化率Δρratei-miは同じ値になる。
また、位置算出部7は、擬似距離誤差算出部4により算出された衛星搭載時計誤差を用いて、GPS受信機が内蔵している時計誤差を算出し、その時計誤差を信号処理部2に出力する。
これにより、信号処理部2は、電波の受信時刻を衛星位置・速度算出部11や擬似距離誤差算出部4に出力する際、位置算出部7により算出された時計誤差を考慮して、電波の受信時刻を補正することになる。
一方、収束していなければ、ステップST31〜ST38の処理を繰り返し実施させる制御信号を衛星位置・速度算出部11や擬似距離平滑化部5などに出力する。
なお、位置算出部7は、電波を受信できているGPS衛星の個数が4衛星未満、または、収束計算の回数が所定回数以上であれば、自動車の位置を測位することができないと判断して、測位周期毎の動作を終了する。
即ち、この参考例2では、重み付き最小二乗法で自動車位置を収束計算する過程で、擬似距離変化推定部12により推定された擬似距離変化率と自動車速度が徐々に正確になり、自動車位置が収束した時点では、GPS衛星毎の擬似距離誤差変化率の誤差が均等になる。また、擬似距離変化推定部12により推定された擬似距離変化率を用いて平滑化された擬似距離も同様に、GPS衛星毎の擬似距離の誤差が均等になる。
このような擬似距離を用いて擬似距離誤差を求めるので、GPS受信機により計測された擬似距離の誤差をより正確に検出することができる。
このような擬似距離を用いて自動車の位置計算を行えば、擬似距離誤差の大きさが大幅に小さくなるだけでなく、三角測量の原理で算出した自動車位置の3次元空間の重心が、より正確な自動車の位置を示すようになる(測位精度が向上する)。
また、この参考例2では、式(19)又は式(21)で擬似距離変化率を推定する際に、GPS測位で算出した自動車位置と自動車速度を使用しているが、例えば、GPS受信機をカーナビなどに組み込んで使用する場合には、自律測位やマップマッチングで算出した位置や速度を適宜使うように切り替えることで、マルチパスの影響をさらに軽減したり、GPS衛星の電波受信状態が悪い場所でも、自動車の位置や速度を連続して計測したりすることができる。また、高い精度を安定して維持することができる効果があるなど、測位精度をさらに向上するための拡張性を有する。
安価なGPS受信機を自動車に搭載して市街地を走行する場合、GPSアンテナ1により受信された複数個のGPS衛星の電波の中には、異常を含む擬似距離ρCτがある(図9を参照)。同じGPS衛星の擬似距離変化率Δρrateについては、擬似距離ρCτが異常となった同じ時間帯に、程度は小さいが異常を生じる(図10を参照)。
図10には、擬似距離の差分ΔρCτを擬似距離変化率Δρrateと合わせてグラフ化しているが、擬似距離変化率Δρrateと比較すれば、擬似距離ρCτの変動が大きいことがわかる。このような擬似距離ρCτと擬似距離変化率ΔρrateがGPS受信機から出力された場合には、従来の方法では、次のような問題がある。
擬似距離変化率Δρrateには、擬似距離ρCτと同じ誤差が生じていないので、上記の特許文献2の方法では、擬似距離ρCτの大きな誤差を検出できない。これは測位精度の低下につながる。
上記の特許文献3であれば、平滑化した擬似距離にマルチパス誤差が残留して擬似距離誤差が有色化する。これは測位精度の低下につながる。
式(19)より、図10で擬似距離変化率が異常である時間帯を含めて、自動車が停車していた場合の擬似距離変化率Δρrate-siを高精度に推定することができる(図11を参照)。
また、擬似距離変化推定部12により推定された擬似距離変化率Δρrate-siは、複数個のGPS衛星と自動車の相対運動より算出したものなので、他のGPS衛星の擬似距離変化率についても、同様に推定することができる。
擬似距離変化推定部12により推定された擬似距離変化率Δρrate-siを用いれば、図示していないが、式(20)により、自動車速度も高精度に算出することができる。
また、擬似距離変化推定部12により推定された擬似距離変化率Δρrate-miは、複数個のGPS衛星と自動車の相対運動より算出したものなので、他のGPS衛星の擬似距離変化率についても、同様に推定することができる。
擬似距離変化推定部12により推定された擬似距離変化率Δρrate-miを基準にすれば、下記の式(26)に示すように、信号処理部2により計測された擬似距離変化率Δρrateの誤差も検出することができる。
δΔρrate=|Δρrate−Δρrate-mi| [m/s] (26)
擬似距離平滑化部14により平滑化された擬似距離ρsiの差分ΔρSiと、信号処理部2により計測された擬似距離ρCτiの差分ΔρCτiを比較すれば、図14より、平滑化された擬似距離ρSiの誤差が大幅に低減したことがわかる。
式(21)による擬似距離変化率Δρrate-miと比較すれば、平滑化された擬似距離ρSiに約10m以下で変動する誤差が残っているが、この誤差は、下記の式(27)により検出することができる。
δΔρSi=|ΔρSi−Δρrate-mi| [m/s] (27)
自動車上空の複数個のGPS衛星(PRN:5,6,9,14,18,22,25,30)について、式(21)で推定した擬似距離変化率Δρrate-miを基準にして、擬似距離ρCτiの差分ΔρCτiとの差異、平滑化した擬似距離ρSiの差分ΔρSiとの差異のそれぞれの標準偏差を計算すれば、擬似距離平滑化部14により平滑化された擬似距離ρSiには、各GPS衛星の不均等な擬似距離誤差が均等になる改善効果があることがわかる。なお、図中において、平滑化なしの欄は、下記の式(28)、平滑化ありの欄は上記の式(27)の結果である。
δΔρCτi=|ΔρCτi−Δρrate-mi| [m/s] (28)
Claims (1)
- 複数のGPS衛星から放射された電波を受信する電波受信手段と、
上記電波受信手段により受信された電波の放射元であるGPS衛星を識別して、そのGPS衛星から放射された電波が地上に届く伝播時間に基づいてGPS衛星の擬似距離を算出し、その電波の搬送波周波数がドップラー偏位した周波数から、擬似距離変化率を算出するとともに、その電波に含まれている航法メッセージを抽出する信号処理手段と、
上記信号処理手段により抽出された航法メッセージからGPS衛星の位置を特定する衛星位置特定手段と、
上記信号処理手段が計測した擬似距離変化率により、該擬似距離変化率とともに計測した擬似距離を平滑化する平滑化手段と、
該平滑化手段により平滑化された擬似距離と上記衛星位置特定手段により特定されたGPS衛星の位置から現在位置を測位する位置測位手段とを備え、
上記平滑化手段は、上記信号処理手段が計測した第一の擬似距離と、該第一の擬似距離を計測する前に該信号処理手段が計測した第二の擬似距離との差分を求め、上記第一の擬似距離とともに上記信号処理手段が計測した第一の擬似距離変化率と上記差分との差異が、所定値より小さいかまたは該所定値と等しい場合、上記第二の擬似距離、上記第一の擬似距離および上記第一の擬似距離変化率から重み係数を求め、その重み係数に応じた擬似距離と擬似距離変化率の重み加算を用いて、上記GPS衛星の擬似距離を平滑化し、一方、上記第一の擬似距離とともに上記信号処理手段が計測した第一の擬似距離変化率と上記差分との差異が、所定値より大きい場合、上記第一の擬似距離を使用しないで上記第一の擬似距離変化率を用いて、上記GPS衛星の擬似距離を平滑化するGPS受信機。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010258962A JP4988028B2 (ja) | 2010-11-19 | 2010-11-19 | Gps受信機 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010258962A JP4988028B2 (ja) | 2010-11-19 | 2010-11-19 | Gps受信機 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2005147083A Division JP2006322846A (ja) | 2005-05-19 | 2005-05-19 | Gps受信機 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2011053223A JP2011053223A (ja) | 2011-03-17 |
JP4988028B2 true JP4988028B2 (ja) | 2012-08-01 |
Family
ID=43942356
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010258962A Active JP4988028B2 (ja) | 2010-11-19 | 2010-11-19 | Gps受信機 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4988028B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10761215B2 (en) | 2016-03-09 | 2020-09-01 | Mitsubishi Electric Corporation | Positioning device and positioning method |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112230254B (zh) * | 2020-10-29 | 2024-05-14 | 新乡医学院三全学院 | 一种gps载波相位多径误差的校正方法及装置 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07198821A (ja) * | 1994-01-06 | 1995-08-01 | Japan Radio Co Ltd | Gps受信機及びその測位方法 |
JPH11118903A (ja) * | 1997-10-16 | 1999-04-30 | Matsushita Electric Ind Co Ltd | 位置検出装置 |
JP2000193733A (ja) * | 1998-12-25 | 2000-07-14 | Japan Radio Co Ltd | 測位装置及び測位方法 |
JP2000193732A (ja) * | 1998-12-25 | 2000-07-14 | Furuno Electric Co Ltd | 測位装置 |
JP2001004734A (ja) * | 1999-06-16 | 2001-01-12 | Matsushita Electric Ind Co Ltd | Gps受信機 |
JP2002196060A (ja) * | 2000-12-25 | 2002-07-10 | Furuno Electric Co Ltd | キャリアスムージングディファレンシャル測位装置 |
JP2003021672A (ja) * | 2001-07-06 | 2003-01-24 | Matsushita Electric Ind Co Ltd | 航法衛星システム受信機 |
JP2004239841A (ja) * | 2003-02-07 | 2004-08-26 | Furuno Electric Co Ltd | 波浪計測システムおよび波浪計測方法 |
JP2006010573A (ja) * | 2004-06-28 | 2006-01-12 | Furuno Electric Co Ltd | キャリアスムージング測位装置 |
JP2006322846A (ja) * | 2005-05-19 | 2006-11-30 | Mitsubishi Electric Corp | Gps受信機 |
-
2010
- 2010-11-19 JP JP2010258962A patent/JP4988028B2/ja active Active
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10761215B2 (en) | 2016-03-09 | 2020-09-01 | Mitsubishi Electric Corporation | Positioning device and positioning method |
Also Published As
Publication number | Publication date |
---|---|
JP2011053223A (ja) | 2011-03-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2006322846A (ja) | Gps受信機 | |
US7498979B2 (en) | Fast decimeter-level GNSS positioning | |
JP7153427B2 (ja) | 衛星測位システムを用いた測位方法および測位装置 | |
US7982667B2 (en) | Post-processed accuracy prediction for GNSS positioning | |
JP3408593B2 (ja) | 衛星をベースとするナビゲーションシステムにおいて衛星の位置を予測する方法及び装置 | |
US8174437B2 (en) | System and method for augmenting DGNSS with internally-generated differential correction | |
US8378887B2 (en) | Systems and methods for synthesizing GPS measurements to improve GPS location availability | |
US9244174B2 (en) | Mitigation of scintillations in signals of global navigation satellite systems caused by ionospheric irregularities | |
KR101437346B1 (ko) | Gps 수신기의 위치 계산 방법, 그 방법을 구현하기 위한 프로그램이 기록된 기록매체 및 gps 수신기 | |
JPH075240A (ja) | Gps位置推定の精度を改善する方法及び装置 | |
KR101949274B1 (ko) | 정지 궤도 위성의 궤도 결정 장치 및 방법 | |
JP2010151725A (ja) | Gnss受信装置及び測位方法 | |
EP3030922A1 (en) | Detection of scintillations in signals of global navigation satellite systems caused by lonospheric irregularities | |
JPH0777570A (ja) | 衛星基地局航法システムにおけるサイクルスリップ検出方法及び装置 | |
EP2864810A1 (en) | Selection of a subset of global navigation satellite system measurements based on relation between shifts in target parameters and sum of residuals | |
US20130002485A1 (en) | Satellite Subset Selection | |
CN104316943A (zh) | 一种伪距离和多普勒组合差分定位系统及方法 | |
JP5564275B2 (ja) | Gps受信機の時刻決定方法 | |
US6704650B1 (en) | Technique for accurate distance and velocity calculations using the global positioning system (GPS) | |
JP5413118B2 (ja) | 測位システム | |
JP4988028B2 (ja) | Gps受信機 | |
JP2008051572A (ja) | ナビゲーション装置及び、その方法、並びにそのプログラム | |
US20230266480A1 (en) | Method for Determining an Integrity Datum regarding a GNSS-Based Localization of a Vehicle | |
JP7148039B2 (ja) | 移動体情報推定装置及びプログラム | |
US20220026583A1 (en) | Enhancing Sensitivity to Reflected GNSS Signals |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110628 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110826 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20120117 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20120306 |
|
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: 20120417 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20120425 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 Ref document number: 4988028 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150511 Year of fee payment: 3 |
|
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 |
|
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 |
|
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 |