JP2995852B2 - 車両用gps航法システム - Google Patents
車両用gps航法システムInfo
- Publication number
- JP2995852B2 JP2995852B2 JP29424290A JP29424290A JP2995852B2 JP 2995852 B2 JP2995852 B2 JP 2995852B2 JP 29424290 A JP29424290 A JP 29424290A JP 29424290 A JP29424290 A JP 29424290A JP 2995852 B2 JP2995852 B2 JP 2995852B2
- Authority
- JP
- Japan
- Prior art keywords
- positioning
- error
- accuracy
- value
- vehicle
- 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 - Lifetime
Links
Landscapes
- Navigation (AREA)
- Fittings On The Vehicle Exterior For Carrying Loads, And Devices For Holding Or Mounting Articles (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Description
に適したGPS航法システムに関する。
3−198887号公報に示されているように、各人工衛星か
らの衛星軌道及び時刻に関する各情報の受信に基いて各
人工衛星の位置と自己位置との距離を測定して自己位置
を測位し、さらに、各人工衛星の位置から誤差の生じる
方向を求め、この誤差方向との関連により道路地図情報
を利用して測位値を補正するようにしたものがある。
を利用して測位する場合には、高度値を与えた上で同一
高度面内で二次元測位を行うこととなるため、測位結果
には、上述の距離の測定誤差に加え、高度値の設定誤差
も含まれることとなる。しかしながら、上述の構成で
は、高度設定誤差が測位結果に与える影響までは考慮さ
れておらず、その結果、二次元測位における測位誤差予
測の精度が低いという不具合があった。
両用GPS航法システムにおいて、二次元測位をするにあ
たり、高度設定誤差をも考慮して、測位誤差の予測精度
を高めようとするものである。
て示すごとく、 車両に装備されて複数の人工衛星S1,S2,…,Siからの
送信電波を受信する受信手段1と、 この受信手段1の受信電波に基いて車両の位置を測位
するとともにこの測位値の誤差の予測精度を演算する測
位演算手段2とを備えたGPS航法システムにおいて、 測位演算手段2は、高度設定誤差を予測する予測手段
2aと、この予測手段2aによって予測された高度設定誤差
により生じる測距精度を演算する測距精度演算手段2bと
を具備して、二次元測位時に、少なくとも前記予測精度
に前記測距精度を合成して新たな測位誤差の予測精度を
演算することを特徴とする車両用GPS航法システムが提
供される。
定誤差が測位結果に与える影響を考慮した構成を採るの
で、測位誤差の予測精度を高めることができる。その結
果、当該予測精度の活用により、当該車両の位置に対す
るより一層正確な把握が可能となる。
2図は、本発明が車両用ナビゲータに適用された例を示
している。このナビゲータは、方位センサ10を有してお
り、この方位センサ10は、当該車両の方位を検出し方位
検出信号を発生する。車速センサ20は、当該車両の車速
を検出しこの検出結果に比例する車速パルスを順次発生
する。GPS受信装置30は、四つ又は三つの人工衛星から
の各送信電波を受信部30aにて受信し、この受信電波に
基き第3図〜第5図に示す各フローチャートに従い測位
プログラムをマイクロコンピュータ30bにより実行し
て、三次元測位又は二次元測位及びその予測精度等のた
めの演算処理をする。但し、測位プログラムはマイクロ
コンピュータ30bのROMに予め記憶されている。
素子からなるもので、この地図メモリ40には、所定範囲
の地図データ(道路形状データを含む)が予め記憶され
ている。コントロールスイッチ50は、地図メモリ40の記
憶地図データ中から表示地図を選択するとき、或いは初
期値をマイクロコンピュータ60に入力するときに操作さ
れる。マイクロコンピュータ60は、表示プログラムを、
第8図に示すフローチャートに従い、方位センサ10、車
速センサ20、GPS受信装置30、地図メモリ40及びコント
ロールスイッチ50との協働により実行し、この実行中に
おいて、CRT80に接続したCRTコントローラ70の制御に必
要な演算処理をする。但し、表示プログラムはマイクロ
コンピュータ60のROMに予め記憶されている。
走行状態にて本発明装置を作動させれば、マイクロコン
ピュータ30bが第3図のフローチャートに従い測位プロ
グラムの実行をステップ100にて開始するとともに、マ
イクロコンピュータ60が第8図のフローチャートに従い
ステップ200にて表示プログラムの実行を開始する。ス
テップ100での実行開始後、マイクロコンピュータ30b
が、第3図及び第4図にて示すごとく、測位プログラム
を第1測位誤差演算ルーティン110に進める。
いては、そのステップ111にて、GPS受信装置30の受信部
30aにより受信された四つの人工衛星からの各軌道デー
タを表わす各送信電波に基き各人工衛星の各三次元位置
が演算され、また衛星電波の到達時間より各衛星とGPS
受信装置30との間の距離が演算され、かつこれら演算結
果に応じてGPS受信装置30の三次元位置(即ち、経度、
緯度及び高度)が測位置として演算される。但し、送信
電波を受信できている人工衛星が三つしかない場合は、
高度値を任意に設定することにより、三つの人工衛星を
用いてGPS受信装置30の二次元位置(即ち、経度、緯
度)が測位値としてステップ111にて演算される。従っ
て、以下の各ステップ112〜115の演算は、三次元測位を
行った場合と、二次元測位を行った場合で異なる。
2において、GPS受信装置30の測位値を中心として、東、
北及び真上を三軸とするx,y,z直交座標系を定め、GPS受
信装置30の測位値から各人工衛星の位置に向かう単位ベ
クトルをそれぞれ求め、これら各単位ベクトルを上述の
直交座標系の各座標軸成分に分割し、各人工衛星の時計
の時刻とGPS受信装置30の時刻との間の時計ズレ量tに
対する係数「1」を加えてA行列を次のように完成す
る。
=∂ri/∂Lg,βi=∂ri/∂Lt,γi=∂ri/∂hが成立
する。ここに、γiはGPS受信装置30と衛星との間の距
離を表わす。また、Lg,Lt,hは、GPS受信装置30の三次元
位置を特定する経度、緯度及び高度を表わす。また、L
g,Lt,hはx,y,zの各座標にそれぞれ対応する。
A行列に基き、以下のように、分散・共分散行列(AT・
A)-1を演算する。
行列(AT・A)-1の固有値λを演算する。かかる場合、
この固有値λは、以下のように、誤差成分の最大となる
方向の分散λと、誤差成分の最小となる方向の分散λ
sとの二つの値として得られる。
としたときの測位誤差の分散を与えるものであることを
考慮し、λ及びλsに関する各式(3)、(4)にお
いては、User Equivalent Range Error(以下、UEREと
いう)で、測位に使用している衛星のうちの最大値UERE
maxを乗じることにより、λ及びλsを実際の測位誤
差の分散(1σ)として求める。なお、「1σ」は68.3
%の確率を表わす。
れかに基き、URA及びGPS受信装置30の内部で発生する測
距離誤差(以下、受信装置内測距誤差という)に応じて
演算される。
をしていないことを表わす。また、URAは、各人工衛星
からの受信電波で特定される人工衛星毎の測距誤差の予
測値を表わすパラメータであって、各人工衛星の時計及
び詳細軌道データを原因として生じる測距誤差成分を含
む。
すると、次のステップ115において、測位誤差の分散が
最大となる方向、即ち、三次元測位時の測位値が確率1
σにて存在する楕円の長軸長さσLの方向θL、及び測
位誤差の分散が最小となる方向、即ち、前記楕円の短軸
長さσSの方向θSが次の各式(7),(8)によって
演算される。
ュータ30bのROMに予め記憶されている。
の演算について示す。ステップ112では、二次元側位の
場合は高度の測定を行わないのでA行列の要素からγの
項がなくなる。更に測位に使用する衛星が三つであるの
で四番目の衛星の項がなくなり以下のようにA行列を完
成する。
下の様になる。
全く同様に求める事ができる。三次元測位の場合は、以
上で測位演算は、終了する。二次元測位の場合は、測位
誤差の要因として更に、高度設定値の誤差が考えられる
ので、これが測位値に与える誤差を求めるため、測位プ
ログラムを第2測位誤差演算ルーティン120(第3図及
び第5図参照)に進める。このことは、ルーティン120
は二次元測位時のみ実行されることを意味する。
にて、三次元測位時のステップ111にて求めた高度を高
度設定値として採用し、その精度(VDOP×UERE)を高度
設定誤差Δhとする。但し、VDOPは、高さ方向の測位精
度の劣化を表わす指数である。なお、二次元測位が連続
する場合には、上述の高度設定誤差Δhに加え、高度値
の設定に採用した三次元測位値と現在の二次元測位値と
から求められるGPS受信装置30の移動量に比例した高度
設定誤差の増加をも考慮する。この増加は、土地の勾配
に基く値であって、測位精度の予測値よりも大きな誤差
をもつ測位値を出力しないようにするために、山岳地等
を想定した上で、大きめの値とする必要がある。
て、高度設定誤差Δhに起因して生じる各人工衛星毎の
測距誤差Δriを次式(9)に基きGPS受信装置30から各
人工衛星へ向かう各単位ベクトルの高さ方向への方向余
弦(ステップ112におけるγ0,γ1,γ2に相当)及び高
度設定誤差Δhに応じて演算する。
おいて、測距誤差Δriに起因して生じる測位誤差、即
ち、経度方向のズレΔLg、緯度方向のズレΔLt,及びGPS
受信装置30の時計と各人工衛星の時計とのズレΔBを次
式(10)に基きステップ112におけるA行列の二次元測
位相当成分の行列(式(1a))に応じて演算する。
の測位誤差はズレの方向が同一であっても、その大きさ
は高度設定誤差Δhに比例して変化することが、式(1
0)から理解される。しかして、高度設定誤差の予測値
としてその1σの大きさを予測したとすれば、測位値
は、真の位置を中心とし(ΔLt,ΔLg)と(−ΔLt,−Δ
Lg)とを結ぶ線上に確率1σにて得られると予測され
る。なお、上述した各式(9),(10)はマイクロコン
ピュータ30bのROMに予め記憶されている。
行が終了すると、マイクロコンピュータ30bが、ステッ
プ130(第3図参照)にて、第1及び第2の測位誤差演
算ルーティン110,120で得た各測位誤差を第6図にて示
すように合成する。この合成によれば、測距誤差による
測位誤差の分布を示す楕円の中心を、GPS受信装置30の
真の位置を中心とする直交座標系において、(ΔLt,ΔL
g)から(−ΔLt,−ΔLg)だけ移動したときにできる楕
円図形の範囲内に、測位点が1σで含まれると予測され
る。三次元測位の場合は、ステップ110で求めたとお
り、第7図に示した楕円図形の範囲内に、測位点が1σ
で含まれると予測される。このことは、GPS受信装置30
の二次元測位及び三次元測位における測位精度が高精度
にて予測できたことを意味する。但し、現実には真の位
置が不明なため、測位値を中心に予測精度の範囲を描け
ば、この範囲内に真の位置が確率1σで含まれる。な
お、マイクロコンピュータ30bは、ステップ130にて、上
述の測位値及びその予測精度をマイクロコンピュータ30
に付与する。
イクロコンピュータ60が、表示プログラムを推測航法演
算処理ルーティン210に進める。しかして、この推測航
法演算処理ルーティン210においては、マイクロコンピ
ュータ60内に記憶された当該車両の最終位置或いはコン
トロールスイッチ50により入力された当該車両の位置を
初期値とし、方位センサ10からの方位検出信号の値及び
車速センサ20からの各車速パルスに基く車速に応じて当
該車両の相対的移動量を演算して同車両の現在位置を推
測演算する。
後、マイクロコンピュータ60がマップマッチング処理ル
ーティン220の実行に移行する。このマップマッチング
処理ルーティン220においては、方位センサ10からの方
位検出信号及び車速センサ20からの各車速パルスに誤差
が混入している場合には、推測航法演算処理ルーティン
210における当該車両の推測演算位置にも誤差が生じる
とともにこの誤差が累積的に増大してゆくことを考慮し
て、かかる累積位置誤差を、地図メモリ40の道路形状デ
ータに基いて補正する。かかる場合、この補正による現
在位置の精度は、道路形状を参照しているため、非常に
高く維持される。
おいて、マップマッチング処理ルーティン220にて補正
した推測演算処理とマイクロコンピュータ30bからの入
力測位値との間のズレがマイクロコンピュータ30bによ
る予測精度内にあるか否かにつき判別する。しかして、
予測精度内になければ、マイクロコンピュータ60がステ
ップ230にて「NO」と判別し、ステップ240にて、マイク
ロコンピュータ30bによる測位値を当該車両の現在位置
とセットし、ステップ250にて、同現在位置を地図メモ
リ40からの地図データにデータとして合成して表示デー
タを作成するとともにこの表示データを表示出力信号と
して発生する。一方、ステップ230での判別が「YES」に
なる場合には、マイクロコンピュータ60が、ステップ25
0にて、上述のマップマッチングで求めた現在位置を上
述と同様に地図データに合成して表示データを作成し表
示出力信号を発生する。なお、GPS受信装置30による測
位が、人工衛星からの送信電波の遮断等で不可能な場合
には、上述の表示データの作成は、マップマッチングで
求めた現在位置を利用してなされる。
を発生すると、CRTコントローラ70が前記表示データの
内容をCRT80に表示させる。かかる場合、上述のように
ステップ230での判別が「NO」となった場合には、GPS受
信装置30による測位値で特定される当該車両の現在位置
が地図上に表示される。一方、上述のようにステップ23
0での判別が「YES」となった場合には、マップマッチン
グにより補正した現在位置が地図上に表示される。
0の実行のもとに、二次元測位時の高度設定誤差による
測位誤差の影響をも考慮して当該車両の現在位置を決定
するようにしたので、従来のように、二次元測位時にHD
OPが小さいときに測位値の精度を良いものとして誤判断
してしまうこともなく、現在位置の予測精度を高めるこ
とができる。従って、CRT80の表示内容によって当該車
両の位置を精度よく認識できる。かかる場合、意図的な
精度劣化(Selective availability)が行なわれたり、
或いは、人工衛星の故障や経年変化等による精度劣化が
生じた場合でも、URAの値を上述のように予測精度のパ
ラメータとしてとりこむことにより正確な予測が可能と
なる。
定にあたり、三次元測位時の高度値を採用した例につい
て説明したが、これに代えて、予め地図データに含めた
高度値、気圧の変化等を利用した高度計の出力、或いは
サインポスト等の位置情報をもつ標識で特定される高度
値を利用して実施するようにしてもよい。
める方法について説明したが、これに限らず、第2測位
誤差演算ルーティン120の活用により、正確な経度・緯
度が得られていることを前提に、高度値をも求めるよう
にしてもよい。かかる場合、例えば、サインポストから
の情報或いはマップマッチングで交差点の右左折等によ
り確度の高い位置決定ができたとき等において、正確な
経度・緯度が分かったときにGPS受信装置30の二次元測
位が終了しておりその測位値が一定方向にずれている場
合、高度設定値に誤差が存在すると考えられる。これに
対して、実際の高度値を求めるにあたっては、ステップ
121におけるΔhを「1」としたときの測位誤差 を演算する。そして、GPS受信装置30の測位値と、その
他の情報から得られた真の位置との距離を求め、この距
離を前記演算測位誤差で除算すれば、高度の設定誤差が
得られる。なお、測位値の真値に対するズレの方向がΔ
Lg,ΔLtの示す方向であれば、高度の設定値が大きい。
一方、逆方向のズレの場合には、高度の設定値が小さ
い。
定した現在位置が、GPS受信装置30の予測精度1σの外
に存在する場合には、GPS受信装置30の測位置を当該車
両の現在位置とするようにしたが、地図が詳細に作成さ
れているために、正確なマップマッチングが可能な場合
には、GPS受信装置30の測位値の使用範囲を、予測精度
1σから2σ或いは3σに拡大するようにして実施して
もよい。
は本発明を適用した車両用ナビゲータのブロック図、第
3図〜第5図は第2図のGPS受信装置のマイクロコンピ
ュータの作用を示すフローチャート、第6図は二次元測
位時の測位誤差の合成説明図、第7図は三次元測位時の
測位誤差の説明図、及び第8図は第2図のマイクロコン
ピュータ60の作用を示すフローチャートである。 符号の説明 30……GPS受信装置、30a……受信部、30b……マイクロ
コンピュータ。
Claims (1)
- 【請求項1】車両に装備されて複数の人工衛星からの送
信電波を受信する受信手段と、 この受信手段の受信電波に基いて車両の位置を測位する
とともにこの測位値の誤差の予測精度を演算する測位演
算手段とを備えたGPS航法システムにおいて、 前記測位演算手段は、高度設定誤差を予測する予測手段
と、この予測手段によって予測された高度設定誤差によ
り生じる測距精度を演算する測距精度演算手段とを具備
して、二次元測位時に、少なくとも前記予測精度に前記
測距精度を合成して新たな測位誤差の予測精度を演算す
ることを特徴とする車両用GPS航法システム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP29424290A JP2995852B2 (ja) | 1990-10-30 | 1990-10-30 | 車両用gps航法システム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP29424290A JP2995852B2 (ja) | 1990-10-30 | 1990-10-30 | 車両用gps航法システム |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH04166781A JPH04166781A (ja) | 1992-06-12 |
JP2995852B2 true JP2995852B2 (ja) | 1999-12-27 |
Family
ID=17805191
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP29424290A Expired - Lifetime JP2995852B2 (ja) | 1990-10-30 | 1990-10-30 | 車両用gps航法システム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2995852B2 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4738944B2 (ja) * | 2005-09-02 | 2011-08-03 | クラリオン株式会社 | Gpsレシーバ |
JP6546730B2 (ja) * | 2014-10-09 | 2019-07-17 | 日本無線株式会社 | 衛星信号受信装置 |
-
1990
- 1990-10-30 JP JP29424290A patent/JP2995852B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JPH04166781A (ja) | 1992-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US4924699A (en) | Position measuring method using a satellite | |
US5400254A (en) | Trace display apparatus for a navigation system | |
JP3062301B2 (ja) | Gpsナビゲーション装置 | |
EP0527558B1 (en) | GPS navigation system with local speed direction sensing and PDOP accuracy evaluation | |
JP3483962B2 (ja) | 航法装置 | |
US4837700A (en) | Method and apparatus for processing data in a GPS receiving device in a road vehicle | |
US5483456A (en) | Navigation system and a method of calculating GPS measuring deviation | |
US5442559A (en) | Navigation apparatus | |
JP3164659B2 (ja) | 航法装置 | |
JPH0531924B2 (ja) | ||
JP2995852B2 (ja) | 車両用gps航法システム | |
JP2783924B2 (ja) | 車両位置検出装置 | |
JPH065169B2 (ja) | 相対位置ナビゲーションシステム | |
JPH0816698B2 (ja) | 移動体の衛星航法装置 | |
JPH07248230A (ja) | 航法装置 | |
JP3019719B2 (ja) | 位置測定方法 | |
JPH0782085B2 (ja) | 衛星航法装置 | |
JP2921355B2 (ja) | 車載用ナビゲーション装置 | |
JPS63109381A (ja) | Gps受信装置のデ−タ処理方法 | |
JPH06281721A (ja) | ナビゲーション装置 | |
JP2786309B2 (ja) | 車両位置検出装置 | |
JPH11142167A (ja) | 車両位置計測装置 | |
JP3263437B2 (ja) | 車両用ナビゲーション装置 | |
JP2593224B2 (ja) | 車両搭載型ナビゲータ装置 | |
JPH0820263B2 (ja) | 車両方位修正装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 9 Free format text: PAYMENT UNTIL: 20081029 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091029 Year of fee payment: 10 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101029 Year of fee payment: 11 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 11 Free format text: PAYMENT UNTIL: 20101029 |
|
EXPY | Cancellation because of completion of term | ||
FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 12 Free format text: PAYMENT UNTIL: 20111029 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111029 Year of fee payment: 12 |