JPH0868654A - Current position detector for vehicle - Google Patents

Current position detector for vehicle

Info

Publication number
JPH0868654A
JPH0868654A JP20630194A JP20630194A JPH0868654A JP H0868654 A JPH0868654 A JP H0868654A JP 20630194 A JP20630194 A JP 20630194A JP 20630194 A JP20630194 A JP 20630194A JP H0868654 A JPH0868654 A JP H0868654A
Authority
JP
Japan
Prior art keywords
azimuth
error
vehicle
gps
speed
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.)
Granted
Application number
JP20630194A
Other languages
Japanese (ja)
Other versions
JP3218876B2 (en
Inventor
Yoshitaka Ozaki
義隆 尾崎
Takahisa Ozaki
貴久 尾▲崎▼
Kaneyoshi Ueda
兼義 植田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Denso Corp
Original Assignee
NipponDenso Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by NipponDenso Co Ltd filed Critical NipponDenso Co Ltd
Priority to JP20630194A priority Critical patent/JP3218876B2/en
Publication of JPH0868654A publication Critical patent/JPH0868654A/en
Application granted granted Critical
Publication of JP3218876B2 publication Critical patent/JP3218876B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Abstract

PURPOSE: To reduce an error when the accuracy of azimuth/speed measured by GPS decreases when the estimation navigation is compensated by the GPS. CONSTITUTION: An offset error, a distance coefficient error, an absolute azimuth error, and an absolute position error in the estimation navigation are obtained for compensation by a Karman filter 6 according to the information of the position, azimuth, and vehicle speed obtained from the estimation navigation and the information of the position, azimuth, and vehicle speed outputted from GPS 3 for the estimation navigation by a distance sensor 1, a relative azimuth sensor 2, a relative trace operation part 4, and an absolute trace operation part 5. When the decrease in the accuracy of azimuth/speed measured by the GPS 3 is judged, an error is obtained only by the information of the position of the vehicle.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】本発明は、車両の方位および移動
距離により車両の現在位置を検出する車両用現在位置検
出装置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a vehicle current position detecting device for detecting the current position of a vehicle based on the direction and the moving distance of the vehicle.

【0002】[0002]

【従来の技術】従来、この種の車載用ナビゲーション装
置においては、車両の方位変化量を検出する相対方位セ
ンサ(ジャイロ、ステアリングセンサ、車輪センサ等)
と、車両の速度(距離)を検出する距離センサ(車速セ
ンサ、車輪センサ等)の出力から、車両の位置・方位・
車速等を検出する推測航法が用いられている。
2. Description of the Related Art Conventionally, in this type of vehicle-mounted navigation device, a relative direction sensor (gyro, steering sensor, wheel sensor, etc.) for detecting the amount of change in the direction of the vehicle.
And the output of a distance sensor (vehicle speed sensor, wheel sensor, etc.) that detects the speed (distance) of the vehicle,
Dead reckoning which detects vehicle speed etc. is used.

【0003】この推測航法の出力(位置・方位・車速
等)には、センサの誤差が含まれるため、誤差が生じ
る。特に、位置・方位は積分的に求められるため、誤差
が徐々に増大してしまう。これに対し、GPSは、絶対
的な位置・方位・車速を求めることができるため、GP
Sが測位した場合に推測航法の出力をGPSの出力に合
わせることにより補正が可能である。例えば、推測航法
で得られた位置をマップマッチングにより道路地図上の
道路位置に位置合わせした時の位置と、GPSで得られ
た位置との差が所定値より大きくなった時に、道路地図
上の位置をGPSで得られた位置に修正するようにする
ことができる。
The dead reckoning navigation output (position, direction, vehicle speed, etc.) includes an error of the sensor, and thus an error occurs. In particular, the position / azimuth is obtained in an integrated manner, so the error gradually increases. On the other hand, since GPS can determine the absolute position, direction, and vehicle speed, GP
The correction can be made by matching the dead reckoning output with the GPS output when S is positioned. For example, when the difference between the position obtained by dead reckoning and the position obtained by aligning the road position on the road map by map matching and the position obtained by GPS becomes larger than a predetermined value, The position may be corrected to the position obtained by GPS.

【0004】[0004]

【発明が解決しようとする課題】しかしながら、推測航
法で求められる位置については、GPSの出力により補
正することはできるが、センサの補正はできない。この
ため、GPS非受信時においてセンサ出力の誤差により
推測航法の出力精度が悪いという問題がある。本発明
は、センサ出力の誤差を修正して現在位置を求めるよう
にすることを第1の目的とする。
However, the position obtained by dead reckoning can be corrected by the output of GPS, but not by the sensor. Therefore, there is a problem that the dead-reckoning navigation output accuracy is poor due to an error in the sensor output when GPS is not received. A first object of the present invention is to correct an error in sensor output so as to obtain a current position.

【0005】本発明者等は、上記目的を達成するため、
後述するように、カルマンフィルタを用い、推測航法か
ら求められる車両の方位、位置、速度の情報とGPSか
ら出力される車両の方位、位置、速度の情報との差によ
り、オフセット誤差および距離係数誤差を求めて、オフ
セット補正量および距離係数を修正するようにしたもの
を考えた。
In order to achieve the above object, the present inventors have
As will be described later, using a Kalman filter, the offset error and the distance coefficient error are determined by the difference between the vehicle azimuth, position and speed information obtained from dead reckoning and the vehicle azimuth, position and speed information output from GPS. In consideration of the above, an offset correction amount and a distance coefficient were corrected.

【0006】このカルマンフィルタによる誤差修正にお
いて、GPSで計測する方位・距離係数は、例えば車速
が低下すると精度が低下する。この場合にそのまま誤差
修正を行っていくと推測航法の方位・距離係数に誤差が
生じる。そこで、本発明はGPSで計測する方位・距離
係数の精度が低下した時に、推測航法の方位・距離係数
の誤差を少なくすることを第2の目的とする。
In the error correction by the Kalman filter, the accuracy of the azimuth / distance coefficient measured by GPS decreases as the vehicle speed decreases. In this case, if the error is corrected as it is, an error occurs in the direction / distance coefficient of dead reckoning. Therefore, a second object of the present invention is to reduce the error in the azimuth / distance coefficient of dead reckoning when the accuracy of the azimuth / distance coefficient measured by GPS decreases.

【0007】[0007]

【課題を解決するための手段】本発明は上記目的を達成
するため、請求項1に記載の発明においては、車両の方
位変化量に応じた信号を出力する相対方位センサ(2)
と、車両の走行速度に応じた信号を出力する距離センサ
(1)と、前記相対方位センサからの信号に基づき、オ
フセット補正量によりオフセット補正して方位変化量を
求めるとともにこの方位変化量により車両の方位を特定
し、前記距離センサからの信号に距離係数を乗じて車両
の移動距離を求め、前記車両の方位と移動距離に基づい
て車両の位置を検出する位置検出手段(4,5およびそ
れに対する演算処理)とから構成される推測航法手段
(1,2,4,5およびそれに対する演算処理)と、G
PS衛星からの衛星電波を受信して車両の方位、位置、
速度に関する情報を出力するGPS(3)と、少なくと
もオフセット誤差、車両の方位誤差および距離係数誤差
を状態量X(t)とし、この誤差の時間的変化を与える
プロセス行列φおよび信号生成過程で発生する雑音ωに
より、状態量X(t+1)をφ・X(t)+ωにて関係
付ける信号生成過程と、状態量X(t)と観測値Y
(t)とを、観測行列Hおよび観測過程で発生する雑音
vにより、観測値Y(t)をH・X(t)+vにて関係
付ける観測過程を形成するモデルを基に、前記推測航法
手段における車両の方位、位置、速度の情報と前記GP
Sから出力される車両の方位、位置、速度の情報との差
により、上記観測値Y(t)を演算するとともに、その
演算値を基に上記モデルから状態量X(t)を演算する
状態量演算手段(6およびそれに対する401〜40
8,410,411,501〜503の処理)と、この
状態量演算手段にて求めた状態量X(t)によるオフセ
ット誤差および距離係数誤差により前記オフセット補正
量および前記距離係数を修正する修正手段(409)と
を備え、前記状態量演算手段は、前記GPSから出力さ
れる方位、速度の精度が低下している状態を判定する
と、前記推測航法手段からの車両の方位、速度に関する
情報と、前記GPSから出力される方位、速度に関する
情報を除外して、前記車両の位置の情報の差に基づき、
前記オフセット誤差および距離係数誤差の状態量を演算
する(501,503)ことを特徴としている。
In order to achieve the above object, the present invention provides a relative azimuth sensor (2) for outputting a signal according to an amount of azimuth change of a vehicle.
And a distance sensor (1) that outputs a signal according to the traveling speed of the vehicle, and based on a signal from the relative azimuth sensor, an offset correction amount is offset-corrected to obtain an azimuth change amount, and the vehicle is detected based on the azimuth change amount. Position detecting means (4, 5 and its) for detecting the position of the vehicle based on the azimuth and the moving distance of the vehicle by determining the azimuth of the vehicle, multiplying the signal from the distance sensor by the distance coefficient to obtain the moving distance of the vehicle. And a dead reckoning means (1, 2, 4, 5 and calculation processing for it), and G
By receiving satellite radio waves from the PS satellite,
A GPS (3) that outputs information related to speed and at least an offset error, a heading error of a vehicle, and a distance coefficient error as a state quantity X (t), and generate a process matrix φ that gives a temporal change of this error and a signal generation process. Signal ω, the state quantity X (t + 1) is related to φ · X (t) + ω, and the state quantity X (t) and the observed value Y
The dead reckoning navigation is based on a model that forms an observation process in which (t) and the observation matrix H and noise v generated in the observation process relate the observed value Y (t) with H · X (t) + v. Information of the direction, position and speed of the vehicle in the means and the GP
A state in which the observed value Y (t) is calculated from the difference between the vehicle azimuth, position, and speed information output from S, and the state quantity X (t) is calculated from the model based on the calculated value. Quantity calculating means (6 and 401 to 40 corresponding thereto)
8, 410, 411, 501-503) and the correction means for correcting the offset correction amount and the distance coefficient by the offset error and the distance coefficient error by the state quantity X (t) obtained by the state quantity calculating means. (409), the state quantity computing means, when determining a state in which the accuracy of the heading and speed output from the GPS is reduced, information on the heading and speed of the vehicle from the dead reckoning means, Excluding the information on the direction and speed output from the GPS, based on the difference in the position information of the vehicle,
The state quantities of the offset error and the distance coefficient error are calculated (501, 503).

【0008】請求項2に記載の発明では、請求項1に記
載の発明において、前記状態量演算手段は、車両の速度
が所定値以下の時に前記GPSから出力される方位、速
度の精度が低下している状態を判定することを特徴とし
ている。請求項3に記載の発明では、請求項1に記載の
発明において、前記状態量演算手段は、前記GPSから
の方位、速度に関する情報を基に、前記GPSから出力
される方位、速度の精度が低下している状態を判定する
ことを特徴としている。
According to a second aspect of the present invention, in the first aspect of the present invention, the state quantity computing means reduces the accuracy of the azimuth and speed output from the GPS when the speed of the vehicle is below a predetermined value. The feature is that the state in which it is operating is determined. According to a third aspect of the present invention, in the first aspect of the present invention, the state quantity computing means determines the accuracy of the azimuth and speed output from the GPS based on the information on the azimuth and speed from the GPS. It is characterized by determining the state of decrease.

【0009】なお、上記各手段のカッコ内の符号は、後
述する実施例記載の具体的手段との対応関係を示すもの
である。また、後述するフローチャート内の各ステップ
はそれぞれの機能を実現する機能実現手段を構成するも
のである。
The reference numerals in parentheses of the above-mentioned means indicate the correspondence with the concrete means described in the embodiments described later. Further, each step in the flowcharts described later constitutes a function realizing means for realizing each function.

【0010】[0010]

【発明の作用効果】請求項1乃至3に記載の発明によれ
ば、前記推測航法手段からの車両の方位、速度に関する
情報と、前記GPSから出力される方位、速度に関する
情報を除外して、前記車両の位置の情報の差により、オ
フセット誤差および距離係数誤差を求め、オフセット補
正量および距離係数を修正するようにしているから、相
対方位センサおよび距離センサに対するセンサ出力の誤
差修正を行うことができる。
According to the invention as set forth in claims 1 to 3, information relating to the direction and speed of the vehicle from the dead reckoning means and information relating to the direction and speed output from the GPS are excluded. Since the offset error and the distance coefficient error are obtained from the difference in the vehicle position information and the offset correction amount and the distance coefficient are corrected, the error in the sensor output with respect to the relative direction sensor and the distance sensor can be corrected. it can.

【0011】さらに、GPSから出力される方位、速度
の精度が低下している状態を判定すると、推測航法手段
およびGPSからの車両の方位、速度に関する情報を除
外して、車両の位置の情報の差に基づき、オフセット誤
差および距離係数誤差の状態量を演算するようにしてい
るから、GPSで計測する方位・距離係数の精度が低下
した状態でのそれらの誤差を少なくすることができる。
Further, when it is determined that the accuracy of the heading and speed output from the GPS is low, information regarding the heading and speed of the vehicle from the dead reckoning means and the GPS is excluded, and the position information of the vehicle is calculated. Since the state quantities of the offset error and the distance coefficient error are calculated based on the difference, those errors can be reduced in a state where the accuracy of the azimuth / distance coefficient measured by GPS is lowered.

【0012】[0012]

【実施例】以下、本発明を図に示す実施例について説明
する。この実施例は、推測航法とGPSとの複合化を図
るため、カルマンフィルタを用いている。このカルマン
フィルタの概要について説明する。このカルマンフィル
タにおいては、図7に示すように、信号生成過程と観測
過程に分けられる。図において、線形システム(φ)が
あり、そのシステムの状態X(t)に対して、観測行列
Hで関係付けられるX(t)の一部が観測できる場合
に、フィルタはX(t)の最適な推定値を与える。ここ
で、ωは信号生成過程にて発生する雑音であり、vは観
測過程にて発生する雑音である。このフィルタの入力は
Y(t)であり、出力はX(t)の最適推定値である。
Embodiments of the present invention will be described below with reference to the drawings. This embodiment uses a Kalman filter in order to combine dead reckoning and GPS. The outline of this Kalman filter will be described. This Kalman filter is divided into a signal generation process and an observation process, as shown in FIG. In the figure, if there is a linear system (φ) and a part of X (t) related by the observation matrix H can be observed for the state X (t) of the system, the filter of X (t) Gives an optimal estimate. Here, ω is noise generated in the signal generation process, and v is noise generated in the observation process. The input of this filter is Y (t) and the output is the optimal estimate of X (t).

【0013】時刻tまでの情報を用いた状態Xの最適推
定値、すなわち状態量X(t|t)は、数1により求め
られる。
The optimum estimated value of the state X using the information up to the time t, that is, the state quantity X (t | t) is obtained by the equation 1.

【0014】[0014]

【数1】X(t|t)=X(t|t−1)+K(t)
{Y(t)−HX(t|t−1)} ここで、X(t|t−1)は事前推定値、K(t)はカ
ルマンゲインであり、それぞれ数2、数3により表され
る。
## EQU1 ## X (t | t) = X (t | t-1) + K (t)
{Y (t) -HX (t | t-1)} where X (t | t-1) is a pre-estimated value and K (t) is a Kalman gain, which are represented by Equations 2 and 3, respectively. It

【0015】[0015]

【数2】X(t|t−1)=φX(t−1|t−1)## EQU00002 ## X (t | t-1) =. Phi.X (t-1 | t-1)

【0016】[0016]

【数3】K(t)=P(t|t−1)HT (HP(t|
t−1)HT +V)-1 ここで、Pは状態量Xの誤差共分散であり、P(t|t
−1)は誤差共分散の予測値、P(t−1|t−1)は
誤差共分散であり、それぞれ数4、数5により表され
る。
[Number 3] K (t) = P (t | t-1) H T (HP (t |
t−1) H T + V) −1 where P is the error covariance of the state quantity X, and P (t | t
-1) is the predicted value of the error covariance, and P (t-1 | t-1) is the error covariance, which are represented by equations 4 and 5, respectively.

【0017】[0017]

【数4】 P(t|t−1)=φP(t−1|t−1)φT +W[Number 4] P (t | t-1) = φP (t-1 | t-1) φ T + W

【0018】[0018]

【数5】P(t−1|t−1)=(I−K(t−1)
H)P(t−1|t−2) なお、Vは観測過程で発生する雑音vの分散、Wは信号
過程で発生する雑音ωの分散である。また、A(i|
j)は時刻jまでの情報に基づく時刻iでのAの推定値
を表す。なお、添字のT は転置行列を意味し、-1は逆行
列を意味する。Iは単位行列である。
## EQU5 ## P (t-1 | t-1) = (I-K (t-1)
H) P (t-1 | t-2) V is the variance of noise v generated in the observation process, and W is the variance of noise ω generated in the signal process. In addition, A (i |
j) represents the estimated value of A at time i based on the information up to time j. The subscript T means a transposed matrix, and -1 means an inverse matrix. I is a unit matrix.

【0019】さらに、VとWは平均0の白色ガウス雑音
であり、互いに無相関である。上記のようなカルマンフ
ィルタにおいて、状態量Xと誤差共分散Pの初期値に適
当な誤差を与えてやり、新しい観測が行われる度に以上
の計算を繰り返し行うことにより、状態量Xの精度が向
上する。このようなカルマンフィルタを推測航法へ適用
したのが本実施例である。
Further, V and W are white Gaussian noises having an average of 0 and are uncorrelated with each other. In the Kalman filter as described above, the accuracy of the state quantity X is improved by giving an appropriate error to the initial values of the state quantity X and the error covariance P and repeating the above calculation each time a new observation is made. To do. This embodiment applies such a Kalman filter to dead reckoning.

【0020】まず、上記の信号生成過程の定義について
説明する。推測航法でのカルマンフィルタは、推測航法
の誤差の補正を目的とするので、状態量Xは以下の5つ
の誤差値を定義する。この誤差値の時間的な変化を与え
るものがプロセス行列φである。 オフセット誤差(εG)
First, the definition of the above signal generation process will be described. The Kalman filter in dead-reckoning aims to correct the error in dead-reckoning, so the state quantity X defines the following five error values. It is the process matrix φ that gives a temporal change of this error value. Offset error (εG)

【0021】[0021]

【数6】εGt =εGt-1 +ω0 確定的な変化はなく、前回の誤差にノイズが付加され
る。 絶対方位誤差(εA)
[Equation 6] εG t = εG t-1 + ω 0 There is no definite change, and noise is added to the previous error. Absolute heading error (εA)

【0022】[0022]

【数7】εAt =T×εGt-1 +εAt-1 +ω1 前回の誤差に、オフセット誤差に前回からの経過時間を
かけて求める方位誤差とノイズが付加される。 距離係数誤差(εK)
In Equation 7] εA t = T × εG t- 1 + εA t-1 + ω 1 last error, the azimuth error and noise to obtain over the elapsed time from the previous to the offset error is added. Distance coefficient error (εK)

【0023】[0023]

【数8】εKt =εKt-1 +ω2 確定的な変化はなく、前回の誤差にノイズが付加され
る。 絶対位置北方向誤差(εY)
[Equation 8] εK t = εK t-1 + ω 2 There is no definite change, and noise is added to the previous error. Absolute position north direction error (εY)

【0024】[0024]

【数9】εYt =sin(AT +εAt-1 +εGt-1 ×
T/2)×L×(1+εKt-1 )−sin(AT )×L
+εYt-1 前回の誤差に方位誤差・距離誤差によって生じる誤差が
付加される。 絶対位置東方向誤差(εX)
## EQU9 ## εY t = sin (A T + εA t-1 + εG t-1 ×
T / 2) × L × (1 + εK t-1 ) −sin (A T ) × L
+ ΕY t-1 The error generated by the azimuth error / distance error is added to the previous error. Absolute position east error (εX)

【0025】[0025]

【数10】εXt =cos(AT +εAt-1 +εGt-1
×T/2)×L×(1+εKt-1 )−cos(AT )×
L+εXt-1 前回の誤差に方位誤差・距離誤差によって生じる誤差が
付加される。上記の定義において、AT は真の絶対方
位、Lは前回からの移動距離、Tは前回からの経過時間
である。
ΕX t = cos (A T + εA t-1 + εG t-1
× T / 2) × L × (1 + εK t-1 ) -cos (A T ) ×
L + εX t-1 The error caused by the azimuth error / distance error is added to the previous error. In the above definition, AT is the true absolute azimuth, L is the moving distance from the previous time, and T is the elapsed time from the previous time.

【0026】上記の各式を状態量で偏微分し線形化する
と信号生成過程は以下のように定義される。
When the above equations are partially differentiated by the state quantity and linearized, the signal generation process is defined as follows.

【0027】[0027]

【数11】 [Equation 11]

【0028】上記Aは、絶対方位AT +εAt-1 +εG
t-1 ×T/2を意味する。この値は真の絶対方位AT
センサ誤差が加わったものであり、後述するように、方
位変化量から求められる絶対方位Aとする。また、ω0
は、オフセット雑音(温度ドリフト等によるオフセット
の変動分)、ω1 は絶対方位雑音(ジャイロのゲイン的
な誤差)、ω2 は距離係数雑音(経年変化)を意味す
る。
The above A is the absolute azimuth A T + εA t-1 + εG
It means t-1 × T / 2. This value is obtained by adding the sensor error to the true absolute azimuth A T , and is the absolute azimuth A obtained from the amount of azimuth change, as will be described later. Also, ω 0
Indicates offset noise (change in offset due to temperature drift etc.), ω 1 indicates absolute direction noise (gyro gain error), and ω 2 indicates distance coefficient noise (aging).

【0029】次に、上記観測過程の定義について説明す
る。観測値は推測航法の出力と、GPSの出力の差より
求める。それぞれの出力には誤差が含まれるため、観測
値において、推測航法の誤差とGPSの誤差の和が得ら
れる。この観測値Yと状態量Xを関係付け、数12のよ
うに定義される。
Next, the definition of the above observation process will be described. The observation value is obtained from the difference between the dead reckoning output and the GPS output. Since each output includes an error, the sum of the dead reckoning error and the GPS error is obtained in the observed value. The observed value Y and the state quantity X are related to each other and are defined as in Expression 12.

【0030】[0030]

【数12】 [Equation 12]

【0031】但し、観測過程で発生する雑音vはGPS
の雑音であり、数13のように定義される。
However, the noise v generated in the observation process is GPS
Noise, which is defined as

【0032】[0032]

【数13】 [Equation 13]

【0033】以上の定義を基に、カルマンフィルタを用
いた推測航法について説明する。図1に本実施例におけ
る概略構成を示す。この図に示すように、車速センサ
1、ジャイロ2からの信号を基に、相対軌跡演算部4、
絶対位置演算部5での演算が行われ、それらの演算(推
測航法演算)により、車速、相対軌跡、絶対位置、絶対
方位が出力される。また、GPS3からは位置・方位・
車速の出力が得られる。カルマンフィルタ6は、推測航
法により得られた車速、絶対位置・絶対方位の情報およ
びGPS3からの車速、位置・方位の情報を基に、車速
センサ1の距離係数補正、ジャイロ2のオフセット補
正、絶対位置補正、絶対方位補正を行う。
The dead reckoning using the Kalman filter will be described based on the above definitions. FIG. 1 shows a schematic configuration of this embodiment. As shown in this figure, based on the signals from the vehicle speed sensor 1 and the gyro 2, the relative trajectory calculation unit 4,
The absolute position calculation unit 5 performs the calculation, and the vehicle speed, the relative trajectory, the absolute position, and the absolute azimuth are output by the calculation (dead-reckoning navigation calculation). In addition, from GPS3
Output of vehicle speed can be obtained. The Kalman filter 6 corrects the distance coefficient of the vehicle speed sensor 1, the offset correction of the gyro 2, and the absolute position based on the vehicle speed, the absolute position / absolute direction information and the vehicle speed / position / direction information from the GPS 3 obtained by dead reckoning. Correction, absolute azimuth correction.

【0034】このような車載用ナビゲーション装置へカ
ルマンフィルタを適用すると、車速センサ1の距離係数
補正、ジャイロ2のオフセット補正、および絶対方位補
正、絶対位置補正により、数2に示す、事前推定X(t
|t−1)は0となる。従って、数1は数14に示すよ
うになる。
When the Kalman filter is applied to such a vehicle-mounted navigation device, the pre-estimation X (t shown in Formula 2 is obtained by the distance coefficient correction of the vehicle speed sensor 1, the offset correction of the gyro 2, the absolute azimuth correction, and the absolute position correction.
| T−1) becomes 0. Therefore, equation 1 becomes as shown in equation 14.

【0035】[0035]

【数14】X(t|t)=K(t)Y(t) 従って、上記信号生成過程にて定義された5つの誤差値
による状態量Xは、数3〜数5によって求められるカル
マンゲインK(t)および観測値Y(t)により求めら
れる。
[Mathematical formula-see original document] X (t | t) = K (t) Y (t) Therefore, the state quantity X by the five error values defined in the signal generation process is the Kalman gain obtained by the equations 3 to 5. It is determined by K (t) and the observed value Y (t).

【0036】ここで、数3における誤差共分散Pは、数
15により定義される。
Here, the error covariance P in Equation 3 is defined by Equation 15.

【0037】[0037]

【数15】 (Equation 15)

【0038】この誤差共分散PにおけるσGG 2 はオフセ
ット誤差の大きさの見積もりを表し、σAA 2は絶対方位
誤差の大きさの見積もりを表し、σKK 2 は距離係数誤差
の大きさの見積もりを表し、σYY 2 は絶対位置北方向誤
差の大きさの見積もりを表し、σXX 2 は絶対位置東方向
誤差の大きさの見積もりを表す。それら以外のσij 2
i行とj列の相互相関値を表す。例えばσAG 2はオフセ
ット誤差と絶対方位誤差の相互相関値を表す。
In this error covariance P, σ GG 2 represents an estimate of the magnitude of the offset error, σ AA 2 represents an estimate of the magnitude of the absolute azimuth error, and σ KK 2 is an estimate of the magnitude of the distance coefficient error. , Σ YY 2 represents an estimate of the absolute position northward error magnitude, and σ XX 2 represents an absolute position eastward error magnitude estimate. Other than that, σ ij 2 represents the cross-correlation value of the i-th row and the j-th column. For example, σ AG 2 represents the cross-correlation value between the offset error and the absolute azimuth error.

【0039】この誤差共分散Pの値は、数4の計算によ
って更新される。なお、初期値においては、σGG 2 、σ
AA 2、σKK 2 、σYY 2 、σXX 2 の各値を誤差が最大とな
る値に設定しておき、また相互相関値については全て0
に設定しておく。また、数3におけるHは数12で示さ
れる行列を用い、Vについては数13に示されるものを
用いる。また、数4におけるWは数11に示されるωの
分散を用いる。
The value of this error covariance P is updated by the calculation of equation (4). In addition, in the initial value, σ GG 2 , σ
The values of AA 2 , σ KK 2 , σ YY 2 , and σ XX 2 are set to values that maximize the error, and the cross-correlation values are all 0.
Set to. In addition, the matrix shown in Expression 12 is used for H in Expression 3, and the matrix shown in Expression 13 is used for V. Further, W in Expression 4 uses the variance of ω shown in Expression 11.

【0040】観測過程における観測値Yとしては、数1
2に示すように、εADRt −εAGP St、εKDRt −εK
GPSt、εYDRt −εYGPSt、εXDRt −εXGPStを用い
ている。ここで、添字のDRt は時刻tにおいて車速セン
サ1、ジャイロ2からの信号に基づく推測航法にて求め
られた値を意味し、GPStは時刻tにおいてGPS3から
出力される値を意味する。
As the observed value Y in the observation process,
As shown in FIG. 2, εA DRt −εA GP St , εK DRt −εK
GPSt , εY DRt −εY GPSt , εX DRt −εX GPSt are used. Here, the subscript DRt means a value obtained by dead-reckoning navigation based on signals from the vehicle speed sensor 1 and the gyro 2 at time t, and GPSt means a value output from GPS 3 at time t.

【0041】εADRt −εAGPStは、推測航法により求
められた絶対方位とGPS3から出力される方位の差、
すなわち推測航法により求められた絶対方位には真の絶
対方位とその誤差εADRt が含まれており、またGPS
3から出力される方位には真の絶対方位とその誤差εA
GPStが含まれているため、それらの差を取ることにより
εADRt −εAGPStが得られる。
ΕA DRt −εA GPSt is the difference between the absolute azimuth determined by dead reckoning and the azimuth output from GPS3,
That is, the absolute azimuth obtained by dead-reckoning contains the true absolute azimuth and its error εA DRt.
The true absolute azimuth and its error εA
Since GPSt is included, εA DRt −εA GPSt can be obtained by taking the difference between them.

【0042】同様に、εKDRt −εKGPStは、推測航法
により求められる速度とGPS3から出力される速度の
差から求まる距離係数誤差であり、具体的には、(推測
航法による速度−GPSによる速度)/(推測航法によ
る速度)により求められる。また、εYDRt −εYGPSt
は、推測航法により求められる絶対位置のY成分とGP
S3から出力される位置のY成分の誤差の差であり、ε
DRt −εXGPStは、推測航法により求められる絶対位
置のX成分とGPS3から出力される位置のX成分の誤
差の差である。
Similarly, εK DRt −εK GPSt is a distance coefficient error obtained from the difference between the speed obtained by dead reckoning and the speed output from GPS3. Specifically, (speed due to dead reckoning−speed due to GPS ) / (Speed by dead reckoning). Also, εY DRt −εY GPSt
Is the Y component of absolute position and GP obtained by dead reckoning
It is the difference in error of the Y component of the position output from S3, and ε
X DRt −εX GPSt is the difference between the error between the X component of the absolute position obtained by dead reckoning and the X component of the position output from the GPS 3.

【0043】また、数13に示す、観測過程で発生する
雑音vはGPS3の雑音であり、以下のようにして求め
られる。GPS3における擬似距離の計測誤差(UER
E)とHDOP(Horizontal Dilution of Presision)
の関係により測位精度が、UERE×HDOPで求めら
れ、この測位精度を2乗することにより、v2t、v3t
求められる。また、ドップラー周波数の計測誤差とHD
OPの関係より速度精度が、ドップラー周波数の計測誤
差×HDOPで求められ、この速度精度/車速にて距離
係数計測誤差が求められ、これを2乗することによりv
1tが求められる。さらに、車両の速度Vc と速度精度か
ら方位精度がtan-1(速度精度/Vc)で求められ、
この方位精度を2乗するこによりv0tが求められる。
Further, the noise v generated in the observation process shown in the equation 13 is the noise of the GPS 3 and is obtained as follows. Pseudo-range measurement error in GPS3 (UER
E) and HDOP (Horizontal Dilution of Presision)
Therefore, the positioning accuracy is calculated by UARE × HDOP, and by squaring the positioning accuracy, v 2t and v 3t are calculated. In addition, Doppler frequency measurement error and HD
The speed accuracy is calculated from the relationship of OP by the measurement error of the Doppler frequency x HDOP, and the distance coefficient measurement error is calculated by this speed accuracy / vehicle speed.
1t is required. Further, the azimuth accuracy is calculated as tan −1 (speed accuracy / Vc) from the vehicle speed Vc and the speed accuracy,
The square of this azimuth accuracy gives v 0t .

【0044】従って、観測過程におけるεADRt −εA
GPSt、εKDRt −εKGPSt、εYDR t −εYGPSt、εX
DRt −εXGPStおよび上記雑音Vを入力とし、数3〜数
5および数1を実行することにより、信号生成過程にて
定義された5つの誤差値による状態量Xが求められ、こ
れらにより車速センサ1の距離係数補正、ジャイロ2の
オフセット補正、絶対位置補正、絶対方位補正が行われ
る。
Therefore, εA DRt −εA in the observation process
GPSt, εK DRt -εK GPSt, εY DR t -εY GPSt, εX
DRt- εX GPSt and the noise V are input, and by executing the equations 3 to 5 and the equation 1, the state quantity X based on the five error values defined in the signal generation process is obtained. The distance coefficient correction of 1, the offset correction of the gyro 2, the absolute position correction, and the absolute azimuth correction are performed.

【0045】上記の相対軌跡演算、絶対位置演算、カル
マンフィルタはマイクロコンピュータによる演算処理に
て行われるため、以下これについて説明する。図2に推
測航法のメインルーチンの演算処理を示す。ステップ1
00にて方位変化量・移動距離の演算を行う。この処理
の詳細を図3に示す。まず、ステップ101にてジャイ
ロ2の出力角速度にメインルーチンの起動周期TM を掛
けて方位変化量を算出する。次のステップ102にて、
その方位変化量から、オフセット補正量(この補正量に
ついては後述する)にメインルーチンの起動周期TM
掛けたものを引き、方位変化量のオフセット補正を行
う。次のステップ103では、車速センサ1からの車速
パルス数に距離係数(この距離係数についても後述す
る)を掛けて移動距離を算出する。
The above-mentioned relative locus calculation, absolute position calculation and Kalman filter are carried out by the calculation processing by the microcomputer, so that they will be described below. FIG. 2 shows the arithmetic processing of the main routine of dead reckoning. Step 1
At 00, the azimuth change amount / movement distance is calculated. Details of this processing are shown in FIG. First, in step 101, the output angular velocity of the gyro 2 is multiplied by the start cycle T M of the main routine to calculate the bearing change amount. In the next step 102,
The offset correction amount (the correction amount will be described later) multiplied by the start cycle T M of the main routine is subtracted from the direction change amount to perform the offset correction of the direction change amount. In the next step 103, the moving distance is calculated by multiplying the number of vehicle speed pulses from the vehicle speed sensor 1 by a distance coefficient (this distance coefficient will also be described later).

【0046】このステップ100の次に、ステップ20
0の相対軌跡演算処理を行う。この処理の詳細を図4に
示す。まず、ステップ201にて、方位変化量(ステッ
プ102にて求めたもの)を基に相対方位を更新する。
この更新した相対方位およびステップ103にて求めた
移動距離によりステップ202にて相対位置座標の更新
を行う。この更新は、移動距離に対する相対方位のX,
Y成分をそれまでの相対位置座標に加算することにより
行う。この相対位置座標は相対軌跡を求めるたに行うも
ので、その相対軌跡と道路形状との関係により、いわゆ
るマップマッチングが行われる。
After this step 100, step 20
A relative locus calculation process of 0 is performed. Details of this processing are shown in FIG. First, in step 201, the relative azimuth is updated based on the azimuth change amount (obtained in step 102).
The relative position coordinates are updated in step 202 based on the updated relative azimuth and the moving distance obtained in step 103. This update is based on the relative azimuth X,
This is performed by adding the Y component to the relative position coordinates up to that point. This relative position coordinate is used to determine the relative trajectory, and so-called map matching is performed based on the relationship between the relative trajectory and the road shape.

【0047】このステップ200の次に、ステップ30
0の絶対方位・絶対位置の演算処理を行う。この処理の
詳細を図5に示す。まず、ステップ301にて、方位変
化量(ステップ102にて求めたもの)を基に絶対方位
を更新する。この更新した絶対方位およびステップ10
3にて求めた移動距離によりステップ202にて絶対位
置座標の更新を行う。このステップ200の処理にて更
新された絶対方位Aと絶対位置は後述するGPSとの複
合化処理にて利用される。
After step 200, step 30
The absolute azimuth and absolute position of 0 are calculated. Details of this processing are shown in FIG. First, in step 301, the absolute azimuth is updated based on the azimuth change amount (obtained in step 102). This updated absolute orientation and step 10
In step 202, the absolute position coordinates are updated based on the movement distance obtained in 3. The absolute azimuth A and the absolute position updated in the processing of this step 200 are used in the composite processing with GPS described later.

【0048】このGPSとの複合化処理を行うステップ
400の詳細を図6に示す。まず、ステップ401にて
前回の測位又は予測計算からT1 秒経過したか否かを判
定する。これはGPS3の測位が行われる毎にステップ
403〜409にて推測航法の誤差を補正する処理を行
うが、GPS3の測位ができない場合には誤差が大きく
なるため、それに対応した誤差の予測計算をステップ4
10、411にて定期的に行うために設けられている。
FIG. 6 shows details of step 400 for carrying out the compounding process with the GPS. First, in step 401, it is determined whether or not T1 seconds have elapsed since the last positioning or prediction calculation. This performs the process of correcting the dead reckoning error in steps 403 to 409 every time the GPS3 positioning is performed. However, when the GPS3 positioning cannot be performed, the error becomes large. Therefore, the prediction calculation of the error corresponding thereto is performed. Step 4
It is provided to perform regularly at 10 and 411.

【0049】ステップ401の判定がNOになると、ス
テップ402にてGPS3からの測位データがあるか否
かを行う。GPS3からの測位データがあると、ステッ
プ403以降のカルマンフィルタの演算処理に進む。ま
ず、ステップ403にて観測値Yの計算を行う。これ
は、GPS3から出力される速度、位置、方位データお
よび推測航法におけるステップ300の処理にて求めた
絶対方位、絶対位置および図示しない速度演算処理によ
り車速センサ1からの車速パルスに基づく車両の速度と
から、数12に示した、εADRt −εA GPSt、εKDRt
−εKGPSt、εYDRt −εYGPSt、εXDRt −εXGPSt
を計算するとともに、数13に示す、観測過程で発生す
る雑音vをGPS3の測位データ等を基に計算する。
When the determination in step 401 is NO, the scan
Whether there is positioning data from GPS3 at step 402
Do that. If there is positioning data from GPS3,
The processing proceeds to the Kalman filter calculation processing after step 403. Well
First, in step 403, the observation value Y is calculated. this
Is the speed, position, and direction data output from GPS3.
And obtained in the process of step 300 in dead reckoning
Absolute azimuth, absolute position and speed calculation processing not shown
Based on the vehicle speed pulse from the vehicle speed sensor 1,
Therefore, εA shown in Equation 12DRt-ΕA GPSt, ΕKDRt
-ΕKGPSt, ΕYDRt-ΕYGPSt, ΕXDRt-ΕXGPSt
As well as calculating
Noise v is calculated based on GPS3 positioning data and the like.

【0050】ステップ404では、プロセス行列φの計
算を行う。これは、前回のプロセス行列の計算時点から
の移動距離L、経過時間T(これらは図示しない計測処
理により別途求められている)およびステップ301に
て求めた絶対方位Aにより、数11に示すプロセス行列
φを求める。このようにして計算した観測値Yおよびプ
ロセス行列φを基に、上述した数3〜数5の計算を行っ
て数14に示す状態量Xを求める。すなわち、ステップ
405では、数3により誤差共分散Pの予測計算を行
う。ステップ406では、数4によりカルマンゲインK
の計算を行う。ステップ407では、数5により誤差共
分散Pの計算を行う。この後、カルマンゲインKおよび
観測値Yに基づき、ステップ408にて、数14の計算
により状態量Xを求める。この状態量Xは、数11の左
辺に示すように、オフセット誤差(εG)、絶対方位誤
差(εA)、距離係数誤差(εK)、絶対位置北方向誤
差(εY)、絶対位置東方向誤差(εX)を表してい
る。
In step 404, the process matrix φ is calculated. This is the process shown in Formula 11 according to the moving distance L from the time when the previous process matrix was calculated, the elapsed time T (these are separately obtained by the measurement process not shown), and the absolute azimuth A obtained in step 301. Find the matrix φ. Based on the observation value Y and the process matrix φ calculated in this way, the above-described calculations of the expressions 3 to 5 are performed to obtain the state quantity X shown in the expression 14. That is, in step 405, the prediction calculation of the error covariance P is performed by the equation 3. In step 406, the Kalman gain K is calculated by the equation 4.
Calculate. In step 407, the error covariance P is calculated by the equation 5. After that, based on the Kalman gain K and the observed value Y, in step 408, the state quantity X is obtained by the calculation of Expression 14. This state quantity X is, as shown on the left side of the equation 11, offset error (εG), absolute heading error (εA), distance coefficient error (εK), absolute position north direction error (εY), absolute position east direction error ( εX) is represented.

【0051】これらの誤差により、ステップ409に
て、図に示す計算にて推測航法誤差の修正、すなわちジ
ャイロ2のオフセット補正、車速センサ1の距離係数補
正、絶対方位補正、絶対位置補正が行われる。ジャイロ
2のオフセット補正により、ステップ102にて用いら
れるオフセット補正量が修正され、車速センサ1の距離
係数補正により、ステップ103にて用いられる距離係
数が修正され、絶対方位補正により、ステップ301に
て用いられる絶対方位Aが修正され、絶対位置補正によ
りステップ302にて用いられる絶対位置が修正され
る。
Due to these errors, in step 409, the dead reckoning error is corrected by the calculation shown in the figure, that is, the gyro 2 offset correction, the vehicle speed sensor 1 distance coefficient correction, absolute azimuth correction, and absolute position correction are performed. . The offset correction amount used in step 102 is corrected by the offset correction of the gyro 2, the distance coefficient used in step 103 is corrected by the distance coefficient correction of the vehicle speed sensor 1, and the absolute direction correction is performed in step 301 by the absolute direction correction. The absolute orientation A used is modified, and the absolute position correction modifies the absolute position used in step 302.

【0052】上記の処理を、GPS3からの測位データ
が有る毎に繰り返し行い、上記誤差修正を行って、より
正確なる推測航法データを得ることができる。また、G
PS3の測位が長時間できない場合で、ステップ401
の判定がYESとなると、ステップ410、411に進
み、プロセス行列φの計算および誤差共分散Pの予測計
算を行う。これによって、GPS3の測位ができない場
合の誤差に対応した誤差共分散の予測計算を行い、その
後にGPS3が測位できた時に行われるカルマンフィル
タの処理を正確に行えるようにする。
The above process is repeated every time there is positioning data from the GPS 3 to correct the above error, and more accurate dead reckoning data can be obtained. Also, G
If positioning of PS3 is not possible for a long time, step 401
If the determination is YES, the process proceeds to steps 410 and 411 to calculate the process matrix φ and predict the error covariance P. This makes it possible to perform an error covariance prediction calculation corresponding to an error when the GPS 3 cannot be positioned, and to accurately perform the Kalman filter process performed when the GPS 3 is subsequently positioned.

【0053】上記したカルマンフィルタによる誤差修正
において、GPS3で計測する方位・速度は、車速が低
下すると精度が低下する。カルマンフィルタの理論によ
れば、このような時には観測雑音を実際に見合った量ま
で大きくすれば正常に動作可能である。しかしながら、
カルマンフィルタの理論は観測雑音が白色ガウス雑音で
あるという前提に立っており、実際のGPSの誤差はS
A(故意の精度劣化)等の影響がある時は白色ガウス雑
音になっておらず、偏った分布となっている。この場合
に上記のように観測雑音だけで処理を続けると、推測航
法の方位・距離係数に誤差が生じる。
In the error correction by the Kalman filter described above, the accuracy of the azimuth / speed measured by the GPS 3 decreases as the vehicle speed decreases. According to the Kalman filter theory, in such a case, it is possible to operate normally if the observation noise is increased to an actually commensurate amount. However,
The Kalman filter theory is based on the assumption that the observation noise is white Gaussian noise, and the actual GPS error is S
When there is an influence of A (intentional deterioration of accuracy), white Gaussian noise does not occur, and the distribution is biased. In this case, if the processing is continued only with the observation noise as described above, an error occurs in the direction / distance coefficient of dead reckoning.

【0054】また、車速が低い場合にGPSによる補正
を行わない構成とすると、GPSの測位精度が低い市街
地等での性能が低下してしまう。そこで、車速が低下し
GPSで計測する方位・速度の精度が低下した場合は、
これを補正のリファレンスとして利用せず、また車速が
高くGPSで計測する方位・速度の精度が高い場合に
は、これを利用するようにする。
If the GPS correction is not performed when the vehicle speed is low, the performance in an urban area or the like where GPS positioning accuracy is low will deteriorate. Therefore, if the vehicle speed decreases and the accuracy of the direction and speed measured by GPS decreases,
This is not used as a reference for correction, and is used when the vehicle speed is high and the azimuth / speed accuracy measured by GPS is high.

【0055】具体的には、図8に示す処理を行う。図8
において、ステップ501にて車速は規定値以上か否か
を判定する。車速は車速センサ1からの出力により求め
られる。この判定がYESになると、ステップ502に
進み、数12により観測値Yを計算してカルマンフィル
タによる補正を行う。この場合は、上述したものと同様
である。
Specifically, the processing shown in FIG. 8 is performed. FIG.
In step 501, it is determined whether the vehicle speed is equal to or higher than a specified value. The vehicle speed is obtained from the output from the vehicle speed sensor 1. If this determination is YES, the process proceeds to step 502, and the observed value Y is calculated by Expression 12, and the Kalman filter is used for correction. In this case, it is the same as that described above.

【0056】一方、車速が規定値より低く、ステップ5
01の判定がNOになると、ステップ503に進み、数
16により観測値Yを計算してカルマンフィルタによる
補正を行う。
On the other hand, if the vehicle speed is lower than the specified value, step 5
When the determination in 01 is NO, the process proceeds to step 503, the observed value Y is calculated by the equation 16, and the Kalman filter is used for correction.

【0057】[0057]

【数16】 [Equation 16]

【0058】すなわち、車速が低い場合は方位・速度の
精度が低下しているため、方位・距離係数を利用しない
ようにして、観測値Yを計算する。そして、このように
して得られた観測値Yに基づき状態量Xを計算してステ
ップ409による補正を行う。なお、ステップ501の
判定は、GPS3からの受信データにより求められる車
両の方位・速度の精度が低下している状態を判定すれば
よいため、車速が低下していることを判定する以外に、
方位精度、速度精度を示すv0t、v1tを所定の比較値と
比較して、車両の方位・距離係数の精度低下を判定する
ようにしてもよい。
That is, when the vehicle speed is low, the azimuth / speed accuracy is low, so the observed value Y is calculated without using the azimuth / distance coefficient. Then, the state quantity X is calculated based on the observation value Y obtained in this way, and the correction in step 409 is performed. Note that the determination in step 501 may be performed by determining the state in which the accuracy of the azimuth / speed of the vehicle obtained from the reception data from the GPS 3 is reduced, and therefore, in addition to determining that the vehicle speed is reduced,
It is also possible to compare v 0t and v 1t indicating the azimuth accuracy and the speed accuracy with a predetermined comparison value to determine the accuracy decrease of the azimuth / distance coefficient of the vehicle.

【図面の簡単な説明】[Brief description of drawings]

【図1】本発明の一実施例を示す概略構成図である。FIG. 1 is a schematic configuration diagram showing an embodiment of the present invention.

【図2】推測航法のメインルーチンの演算処理を示すフ
ローチャートである。
FIG. 2 is a flowchart showing a calculation process of a main routine of dead reckoning.

【図3】方位変化量・移動距離の演算処理を示すフロー
チャートである。
FIG. 3 is a flowchart showing a calculation process of a direction change amount / movement distance.

【図4】相対軌跡の演算処理を示すフローチャートであ
る。
FIG. 4 is a flowchart showing a relative locus calculation process.

【図5】絶対方位・絶対位置の演算処理を示すフローチ
ャートである。
FIG. 5 is a flowchart showing a calculation process of absolute azimuth and absolute position.

【図6】GPSとの複合化処理を示すフローチャートで
ある。
FIG. 6 is a flowchart showing a combination process with GPS.

【図7】カルマンフィルタのモデルを示す構成図であ
る。
FIG. 7 is a configuration diagram showing a model of a Kalman filter.

【図8】方位・距離係数の精度判定に基づいて行うカル
マンフィルタの処理を示すフローチャートである。
FIG. 8 is a flowchart showing a process of a Kalman filter performed based on accuracy determination of azimuth / distance coefficient.

【符号の説明】[Explanation of symbols]

1 距離センサ 2 相対方位センサ 3 GPS 4 相対軌跡演算部 5 絶対位置演算部 6 カルマンフィルタ 1 Distance Sensor 2 Relative Direction Sensor 3 GPS 4 Relative Trajectory Calculation Section 5 Absolute Position Calculation Section 6 Kalman Filter

Claims (3)

【特許請求の範囲】[Claims] 【請求項1】 車両の方位変化量に応じた信号を出力す
る相対方位センサと、 車両の走行速度に応じた信号を出力する距離センサと、 前記相対方位センサからの信号に基づき、オフセット補
正量によりオフセット補正して方位変化量を求めるとと
もにこの方位変化量により車両の方位を特定し、前記距
離センサからの信号に距離係数を乗じて車両の移動距離
を求め、前記車両の方位と移動距離に基づいて車両の位
置を検出する位置検出手段とから構成される推測航法手
段と、 GPS衛星からの衛星電波を受信して車両の方位、位
置、速度に関する情報を出力するGPSと、 少なくともオフセット誤差、車両の方位誤差および距離
係数誤差を状態量X(t)とし、この誤差の時間的変化
を与えるプロセス行列φおよび信号生成過程で発生する
雑音ωにより、状態量X(t+1)をφ・X(t)+ω
にて関係付ける信号生成過程と、状態量X(t)と観測
値Y(t)とを、観測行列Hおよび観測過程で発生する
雑音vにより、観測値Y(t)をH・X(t)+vにて
関係付ける観測過程を形成するモデルを基に、前記推測
航法手段における車両の方位、位置、速度の情報と前記
GPSから出力される車両の方位、位置、速度の情報と
の差により、上記観測値Y(t)を演算するとともに、
その演算値を基に上記モデルから状態量X(t)を演算
する状態量演算手段と、 この状態量演算手段にて求めた状態量X(t)によるオ
フセット誤差および距離係数誤差により前記オフセット
補正量および前記距離係数を修正する修正手段とを備
え、 前記状態量演算手段は、前記GPSから出力される方
位、速度の精度が低下している状態を判定すると、前記
推測航法手段からの車両の方位、速度に関する情報と、
前記GPSから出力される方位、速度に関する情報を除
外して、前記車両の位置の情報の差に基づき、前記オフ
セット誤差および距離係数誤差の状態量を演算すること
を特徴とする車両用現在位置検出装置。
1. A relative azimuth sensor that outputs a signal in accordance with the amount of change in the azimuth of the vehicle, a distance sensor that outputs a signal in accordance with the traveling speed of the vehicle, and an offset correction amount based on the signal from the relative azimuth sensor. Offset correction is performed to obtain the azimuth change amount, and the azimuth change amount of the vehicle is specified, and the signal from the distance sensor is multiplied by the distance coefficient to obtain the moving distance of the vehicle. Dead reckoning means composed of position detecting means for detecting the position of the vehicle based on the GPS, GPS for receiving satellite radio waves from GPS satellites and outputting information on the azimuth, position and speed of the vehicle, at least an offset error, The heading error and the distance coefficient error of the vehicle are taken as the state quantity X (t), and the process matrix φ that gives a temporal change of this error and the coarseness generated in the signal generation process. By the sound ω, the state quantity X (t + 1) is calculated as φ · X (t) + ω
, The state quantity X (t) and the observed value Y (t) are related to each other by the observation matrix H and the noise v generated in the observed process. ) + V based on a model forming an observation process, which is related to the direction, position, and speed information of the vehicle in the dead reckoning means and the difference between the direction, position, and speed information of the vehicle output from the GPS. , The observation value Y (t) is calculated, and
A state quantity calculating means for calculating the state quantity X (t) from the model based on the calculated value, and the offset correction by the offset error and the distance coefficient error due to the state quantity X (t) obtained by the state quantity calculating means. An amount and a correction means for correcting the distance coefficient, and when the state quantity calculation means determines a state in which the accuracy of the azimuth and speed output from the GPS is reduced, Information about bearing and speed,
Current position detection for vehicle characterized in that the state quantities of the offset error and the distance coefficient error are calculated based on the difference in the information on the position of the vehicle, excluding the information on the azimuth and speed output from the GPS. apparatus.
【請求項2】 前記状態量演算手段は、車両の速度が所
定値以下の時に前記GPSから出力される方位、速度の
精度が低下している状態を判定することを特徴とする請
求項1に記載の車両用現在位置検出装置。
2. The state quantity calculating means determines a state in which the accuracy of the azimuth and speed output from the GPS is reduced when the speed of the vehicle is equal to or lower than a predetermined value. The vehicle current position detection device described.
【請求項3】 前記状態量演算手段は、前記GPSから
の方位、速度に関する情報を基に、前記GPSから出力
される方位、速度の精度が低下している状態を判定する
ことを特徴とする請求項1に記載の車両用現在位置検出
装置。
3. The state quantity calculating means determines a state in which the accuracy of the azimuth and speed output from the GPS is reduced based on the information on the azimuth and speed from the GPS. The vehicle current position detection device according to claim 1.
JP20630194A 1994-08-31 1994-08-31 Current position detection device for vehicles Expired - Lifetime JP3218876B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP20630194A JP3218876B2 (en) 1994-08-31 1994-08-31 Current position detection device for vehicles

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP20630194A JP3218876B2 (en) 1994-08-31 1994-08-31 Current position detection device for vehicles

Publications (2)

Publication Number Publication Date
JPH0868654A true JPH0868654A (en) 1996-03-12
JP3218876B2 JP3218876B2 (en) 2001-10-15

Family

ID=16521043

Family Applications (1)

Application Number Title Priority Date Filing Date
JP20630194A Expired - Lifetime JP3218876B2 (en) 1994-08-31 1994-08-31 Current position detection device for vehicles

Country Status (1)

Country Link
JP (1) JP3218876B2 (en)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6292752B1 (en) 1997-11-06 2001-09-18 Daimlerchrysler Ag Device for acquiring lane path indicative data
JP2001337150A (en) * 2000-03-24 2001-12-07 Clarion Co Ltd Gps receiver outputting 2drms added with error estimation value of kalman filter
JP2007248271A (en) * 2006-03-16 2007-09-27 Seiko Epson Corp Positioning device and positioning method
JP2008175786A (en) * 2007-01-22 2008-07-31 Zhencheng Hu Moving object position detecting method and detector
JP2008275530A (en) * 2007-05-02 2008-11-13 Alpine Electronics Inc Position detector and position detecting method
US7502688B2 (en) 2004-11-12 2009-03-10 Mitsubishi Electric Corporation System for autonomous vehicle navigation with carrier phase DGPS and laser-scanner augmentation
JP2011128138A (en) * 2010-09-30 2011-06-30 Seiko Epson Corp Method and system for calculating position
JP2011128071A (en) * 2009-12-18 2011-06-30 Seiko Epson Corp Method and system for calculating position
JP2013205348A (en) * 2012-03-29 2013-10-07 Toyota Mapmaster:Kk Apparatus and method for estimating position of mobile terminal device, computer program for estimating position of mobile terminal device, and recording medium with computer program recorded thereon
JP5667731B1 (en) * 2014-03-06 2015-02-12 ジオサーフ株式会社 Field guidance system, field guidance method, software, and storage medium storing software
JP2016211943A (en) * 2015-05-08 2016-12-15 株式会社デンソー Vehicle position detection system and vehicle position detection method
US9704391B2 (en) 2010-10-28 2017-07-11 Panasonic Intellectual Property Management Co., Ltd. Traffic accident detection device and method of detecting traffic accident
JP2021167727A (en) * 2020-04-08 2021-10-21 株式会社デンソー Positioning device for vehicle and program for positioning device for vehicle
JP2021167728A (en) * 2020-04-08 2021-10-21 株式会社デンソー Positioning device for vehicle and program for positioning device for vehicle

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6015459B2 (en) 2013-01-18 2016-10-26 株式会社デンソー Vehicle navigation device
JP2014153215A (en) 2013-02-08 2014-08-25 Denso Corp Navigation device for vehicle

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6292752B1 (en) 1997-11-06 2001-09-18 Daimlerchrysler Ag Device for acquiring lane path indicative data
JP2001337150A (en) * 2000-03-24 2001-12-07 Clarion Co Ltd Gps receiver outputting 2drms added with error estimation value of kalman filter
US7502688B2 (en) 2004-11-12 2009-03-10 Mitsubishi Electric Corporation System for autonomous vehicle navigation with carrier phase DGPS and laser-scanner augmentation
JP2007248271A (en) * 2006-03-16 2007-09-27 Seiko Epson Corp Positioning device and positioning method
JP2008175786A (en) * 2007-01-22 2008-07-31 Zhencheng Hu Moving object position detecting method and detector
JP2008275530A (en) * 2007-05-02 2008-11-13 Alpine Electronics Inc Position detector and position detecting method
JP2011128071A (en) * 2009-12-18 2011-06-30 Seiko Epson Corp Method and system for calculating position
JP2011128138A (en) * 2010-09-30 2011-06-30 Seiko Epson Corp Method and system for calculating position
US9704391B2 (en) 2010-10-28 2017-07-11 Panasonic Intellectual Property Management Co., Ltd. Traffic accident detection device and method of detecting traffic accident
JP2013205348A (en) * 2012-03-29 2013-10-07 Toyota Mapmaster:Kk Apparatus and method for estimating position of mobile terminal device, computer program for estimating position of mobile terminal device, and recording medium with computer program recorded thereon
JP5667731B1 (en) * 2014-03-06 2015-02-12 ジオサーフ株式会社 Field guidance system, field guidance method, software, and storage medium storing software
JP2016211943A (en) * 2015-05-08 2016-12-15 株式会社デンソー Vehicle position detection system and vehicle position detection method
JP2021167727A (en) * 2020-04-08 2021-10-21 株式会社デンソー Positioning device for vehicle and program for positioning device for vehicle
JP2021167728A (en) * 2020-04-08 2021-10-21 株式会社デンソー Positioning device for vehicle and program for positioning device for vehicle

Also Published As

Publication number Publication date
JP3218876B2 (en) 2001-10-15

Similar Documents

Publication Publication Date Title
US8510044B2 (en) Position sensing device and method
US8234091B2 (en) Angular velocity sensor correcting apparatus for deriving value for correcting output signal from angular velocity sensor, angular velocity calculating apparatus, angular velocity sensor correcting method, and angular velocity calculating method
JP4952518B2 (en) Vehicle current position detection device and program
US9395191B2 (en) Navigation apparatus
US7286933B2 (en) GPS/dead-reckoning combination system and operating method thereof
JP3218876B2 (en) Current position detection device for vehicles
CN109343095B (en) Vehicle-mounted navigation vehicle combined positioning device and combined positioning method thereof
US8401788B2 (en) Angular velocity sensor correcting apparatus and method
JP3473117B2 (en) Current position detection device for vehicles
JP2000055678A (en) Current position detecting device for vehicle
WO2017150106A1 (en) In-vehicle device and estimation method
JP2011122921A (en) Position location apparatus, position location method, position location program, velocity vector calculation apparatus, velocity vector calculation method, and velocity vector calculation program
CN110133694A (en) The vehicle positioning method and system assisted based on the course double antenna GNSS and wheel speed
JP5273127B2 (en) Angular velocity sensor correction apparatus and angular velocity sensor correction method
KR100948089B1 (en) Method for deciding car position by pseudo dead reckoning and car navigation system using the same
JP6248559B2 (en) Vehicle trajectory calculation device
JP3628046B2 (en) Current position detection device for vehicle
JPH07301541A (en) Navigation device
JP5365606B2 (en) Angular velocity sensor correction apparatus and angular velocity sensor correction method
WO2016098703A1 (en) Angular velocity sensor correction device, angular velocity sensor correction method, azimuth estimation device and azimuth estimation method
JP3321096B2 (en) Vehicle position measuring device and speed measuring device
JP5222814B2 (en) Positioning method and apparatus
JP3503204B2 (en) Current position detection device for vehicles
JP3367461B2 (en) Moving body attitude angle detection device
JP5569276B2 (en) Current position detection device for vehicle

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: 20100810

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100810

Year of fee payment: 9

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110810

Year of fee payment: 10

FPAY Renewal fee payment (prs date is renewal date of database)

Year of fee payment: 11

Free format text: PAYMENT UNTIL: 20120810

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130810

Year of fee payment: 12

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

EXPY Cancellation because of completion of term