JPH03188318A - 方位検出装置 - Google Patents
方位検出装置Info
- Publication number
- JPH03188318A JPH03188318A JP32985389A JP32985389A JPH03188318A JP H03188318 A JPH03188318 A JP H03188318A JP 32985389 A JP32985389 A JP 32985389A JP 32985389 A JP32985389 A JP 32985389A JP H03188318 A JPH03188318 A JP H03188318A
- Authority
- JP
- Japan
- Prior art keywords
- angular velocity
- azimuth
- data
- sensor
- velocity sensor
- 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
- 238000004364 calculation method Methods 0.000 claims abstract description 13
- 238000012545 processing Methods 0.000 claims abstract description 10
- 238000001514 detection method Methods 0.000 claims description 18
- 238000000034 method Methods 0.000 description 13
- 230000006870 function Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 4
- 238000012886 linear function Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 238000009933 burial Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000005389 magnetism Effects 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 238000012887 quadratic function Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Landscapes
- Navigation (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〈産業上の利用分野〉
本発明は方位検出装置に関し、さらに詳細にいえば、地
磁気センサを使用し、この他に移動体の移動方向を検出
して旋回角速度を知る旋回角速度センサ(例えば、光フ
アイバジャイロ、機械式ジャイロ、振動ジャイロ、ガス
レートジャイロ)を使用して移動体のヨ一方向の方位を
検出する方位検出装置に関するものである。
磁気センサを使用し、この他に移動体の移動方向を検出
して旋回角速度を知る旋回角速度センサ(例えば、光フ
アイバジャイロ、機械式ジャイロ、振動ジャイロ、ガス
レートジャイロ)を使用して移動体のヨ一方向の方位を
検出する方位検出装置に関するものである。
〈従来の技術〉
従来から、道路交通網の任意の箇所を走行している車両
、あるいは航空路を航行する航空機、海路を航行する船
舶等(以下「車両」を想定して説明を進める。車両の「
走行」、航空機、船舶の「航行」をまとめていうときは
「走行」という)の位置を検出する方式として、距離セ
ンサと、方位センサ(地磁気センサ、または旋回角速度
センサのいずれか1つ)と、両センサからの出力信号に
必要な処理を施す処理装置とを具備し、車両の走行に伴
なって生ずる距離変化量Δノ、および方位θ(地磁気セ
ンサの場合)、もしくは方位変化量Δθ(旋回角速度セ
ンサの場合)を用いて移動体の現在位置データを得る推
測航法(DeadReckoning )が提案されて
いる。この方式は、Δノおよびθに基いて、例えばΔJ
の東西方向成分ΔX(−ΔJXcosθ)、および南北
方向成分Δy(−ΔJXsinθ)を算出し、従前の位
置出力データ(Px’、 Py’)に対して上記各成分
ΔX、Δyを加算することにより、現在の位置出力デー
タ(Px、Py)を求める方式であるが、方位センサが
必然的に有している誤差のために、得られる現在位置デ
ータに含まれる誤差が累積されてしまうという欠点があ
る。
、あるいは航空路を航行する航空機、海路を航行する船
舶等(以下「車両」を想定して説明を進める。車両の「
走行」、航空機、船舶の「航行」をまとめていうときは
「走行」という)の位置を検出する方式として、距離セ
ンサと、方位センサ(地磁気センサ、または旋回角速度
センサのいずれか1つ)と、両センサからの出力信号に
必要な処理を施す処理装置とを具備し、車両の走行に伴
なって生ずる距離変化量Δノ、および方位θ(地磁気セ
ンサの場合)、もしくは方位変化量Δθ(旋回角速度セ
ンサの場合)を用いて移動体の現在位置データを得る推
測航法(DeadReckoning )が提案されて
いる。この方式は、Δノおよびθに基いて、例えばΔJ
の東西方向成分ΔX(−ΔJXcosθ)、および南北
方向成分Δy(−ΔJXsinθ)を算出し、従前の位
置出力データ(Px’、 Py’)に対して上記各成分
ΔX、Δyを加算することにより、現在の位置出力デー
タ(Px、Py)を求める方式であるが、方位センサが
必然的に有している誤差のために、得られる現在位置デ
ータに含まれる誤差が累積されてしまうという欠点があ
る。
すなわち、方位センサが、地磁気を検出して移動体の絶
対方位を知る地磁気センサである場合には、地磁気方位
センサは微弱な地球磁界の強さを検出するものであり、
移動体本体が着磁してしまうとその出力データには誤差
が発生する。この誤差を打ち消すために地磁気方位セン
サの初期化処理が行われるが、車両であればその走行中
、特に踏切、電カケープル埋設場所、鉄橋、防音壁のあ
る高速道路や高層ビルの谷間を通過する時等にしばしば
外部からの強電磁界の影響を受けて車体の着磁量が変化
することにより、再度誤差が発生することがある。した
がって、このような磁界の乱れを含んだ地磁気センサ出
力データを的確に検出して排除しなければ、正しい方位
を求めることができない。
対方位を知る地磁気センサである場合には、地磁気方位
センサは微弱な地球磁界の強さを検出するものであり、
移動体本体が着磁してしまうとその出力データには誤差
が発生する。この誤差を打ち消すために地磁気方位セン
サの初期化処理が行われるが、車両であればその走行中
、特に踏切、電カケープル埋設場所、鉄橋、防音壁のあ
る高速道路や高層ビルの谷間を通過する時等にしばしば
外部からの強電磁界の影響を受けて車体の着磁量が変化
することにより、再度誤差が発生することがある。した
がって、このような磁界の乱れを含んだ地磁気センサ出
力データを的確に検出して排除しなければ、正しい方位
を求めることができない。
一方、旋回角速度センサを使用する場合には、例えば、
方位変化が所定値以上となった時、電源(イグニション
)オン時、極低速走行時、または、山道などの悪路走行
が検出された時等には、センサの出力データにノイズが
多く現れることが知られており、これらのデータをその
まま使用すると方位検出精度が低下する。
方位変化が所定値以上となった時、電源(イグニション
)オン時、極低速走行時、または、山道などの悪路走行
が検出された時等には、センサの出力データにノイズが
多く現れることが知られており、これらのデータをその
まま使用すると方位検出精度が低下する。
そこで、方位センサとして、旋回角速度センサと、地磁
気センサとの2つの方位センサを併用し、旋回角速度セ
ンサの出力データと、地磁気センサの出力データとのい
ずれかの信頼性が低下した場合に他のデータで補うよう
にした方位検出装置が考えられる。
気センサとの2つの方位センサを併用し、旋回角速度セ
ンサの出力データと、地磁気センサの出力データとのい
ずれかの信頼性が低下した場合に他のデータで補うよう
にした方位検出装置が考えられる。
この方位検出装置に適用しうる1つの方法は、旋回角速
度センサの出力データに含まれるノイズと、地磁気セン
サの出力データに含まれるノイズとを、ともに正規分布
に従うランダムノイズと考えて、カルマンフィルタの理
論を用いて、ノイズの予測値が少ないほうのデータすな
わち信頼性の多いと思われるデータに重みをおいて最終
的な方位を推定する方法である。
度センサの出力データに含まれるノイズと、地磁気セン
サの出力データに含まれるノイズとを、ともに正規分布
に従うランダムノイズと考えて、カルマンフィルタの理
論を用いて、ノイズの予測値が少ないほうのデータすな
わち信頼性の多いと思われるデータに重みをおいて最終
的な方位を推定する方法である。
ところが、旋回角速度センサでは、直線走行中でセンサ
出力が0であるべき時でも、温度や湿度の影響を受けて
幾らかの出力(オフセット)が発生するという傾向があ
る。このオフセット出力は、累積するという性質を有す
るので実際の走行方向からずれた方向を検知してしまう
ことになる。
出力が0であるべき時でも、温度や湿度の影響を受けて
幾らかの出力(オフセット)が発生するという傾向があ
る。このオフセット出力は、累積するという性質を有す
るので実際の走行方向からずれた方向を検知してしまう
ことになる。
しかも、オフセット値、およびその時間変化特性は、温
度や湿度、旋回角速度センサの種類によって異なる。ま
た、同種の旋回角速度センサであっても固体差を反映す
る。
度や湿度、旋回角速度センサの種類によって異なる。ま
た、同種の旋回角速度センサであっても固体差を反映す
る。
さらに、このオフセット出力は、旋回角速度センサ出力
データのランダムノイズとして検出することができない
。この理由は、ランダムノイズの場合、平均値が0にな
るのでデータの分散値をとることにより検出できるが、
オフセット出力では、平均値が0とならないからである
。
データのランダムノイズとして検出することができない
。この理由は、ランダムノイズの場合、平均値が0にな
るのでデータの分散値をとることにより検出できるが、
オフセット出力では、平均値が0とならないからである
。
したがって、カルマンフィルタの理論を用いて、ノイズ
の少ないほうのデータすなわち信頼性の多いほうのデー
タに重みをおいて最終的な方位を推定する場合にあって
は、オフセットに対する何らかの補正処理が望まれる。
の少ないほうのデータすなわち信頼性の多いほうのデー
タに重みをおいて最終的な方位を推定する場合にあって
は、オフセットに対する何らかの補正処理が望まれる。
本発明の目的は、旋回角速度センサと地磁気センサの出
力データを取り込み、それらの値と、過去の推定方位か
ら移動体の現在の方位を算出し、もって、移動体の現在
位置を求める方位検出装置において、上記旋回角速度セ
ンサの出力データの信頼性を判定する場合、実際に使わ
れている旋回角速度センサに生じるオフセットのずれを
含めた形で信頼性を判定し、移動体の現在の方位を正確
に推定することができる方位検出装置を提供することに
ある。
力データを取り込み、それらの値と、過去の推定方位か
ら移動体の現在の方位を算出し、もって、移動体の現在
位置を求める方位検出装置において、上記旋回角速度セ
ンサの出力データの信頼性を判定する場合、実際に使わ
れている旋回角速度センサに生じるオフセットのずれを
含めた形で信頼性を判定し、移動体の現在の方位を正確
に推定することができる方位検出装置を提供することに
ある。
く課題を解決するための手段〉
上記の目的を達成するための本発明の方位検出装置は、
第1図に示すように、旋回角速度センサのオフセット値
を移動体停止時の実測オフセット値と、停止時から次の
停止時までに計測した期間とに基づき、時間の関数o
(t)として予測するオフセット値予測手段15と、地
磁気センサの出力データθ8のノイズ成分σHを記憶す
るとともに、旋回角速度センサの出力データΔθGのノ
イズ成分σGを記憶し、上記予測されたオフセットo
(t)を用いて旋回角速度センサの出力データΔθGの
見掛けのノイズ成分を式 %式%()) により計算するノイズ算出手段12と、上記ノイズ成分
を用いてカルマンフィルタゲインKを算出するフィルタ
ゲイン算出手段13と、地磁気センサの出力から求めた
方位データθ”、および旋回角速度センサの出力から求
めた方位データθ6に、上記カルマンフィルタゲインK
に基づいた重み付け処理をすることにより、移動体の現
在の推定方位θを求める方位推定手段11とを有するも
のである。
第1図に示すように、旋回角速度センサのオフセット値
を移動体停止時の実測オフセット値と、停止時から次の
停止時までに計測した期間とに基づき、時間の関数o
(t)として予測するオフセット値予測手段15と、地
磁気センサの出力データθ8のノイズ成分σHを記憶す
るとともに、旋回角速度センサの出力データΔθGのノ
イズ成分σGを記憶し、上記予測されたオフセットo
(t)を用いて旋回角速度センサの出力データΔθGの
見掛けのノイズ成分を式 %式%()) により計算するノイズ算出手段12と、上記ノイズ成分
を用いてカルマンフィルタゲインKを算出するフィルタ
ゲイン算出手段13と、地磁気センサの出力から求めた
方位データθ”、および旋回角速度センサの出力から求
めた方位データθ6に、上記カルマンフィルタゲインK
に基づいた重み付け処理をすることにより、移動体の現
在の推定方位θを求める方位推定手段11とを有するも
のである。
なお、ノイズ算出手段が地磁気センサの出力データθH
に含まれるノイズ成分σ1、および旋回角速度センサの
出力データΔθGに含まれるノイズ成分σGを、走行デ
ータに基づいて算出するものであってもよい。
に含まれるノイズ成分σ1、および旋回角速度センサの
出力データΔθGに含まれるノイズ成分σGを、走行デ
ータに基づいて算出するものであってもよい。
く作用〉
地磁気センサの方位データθH1および旋回角速度セン
サの出力から求められる方位データθGは、それぞれラ
ンダムノイズ、オフセットを含めた形で次のように表わ
される。
サの出力から求められる方位データθGは、それぞれラ
ンダムノイズ、オフセットを含めた形で次のように表わ
される。
θH−θ+V
θG−〇+Δθ’+W+o(t)
ここに、v、Wは、地磁気センサ、旋回角速度センサの
出力に含まれるランダムノイズであり、それらの平均値
はOである。o (t)はオフセット値であり、時間と
ともに変化する。
出力に含まれるランダムノイズであり、それらの平均値
はOである。o (t)はオフセット値であり、時間と
ともに変化する。
本発明の方位検出装置によれば、ノイズ算出手段12に
より、上記ランダムノイズ成分v、Wとして、地磁気セ
ンサの出力データθH1旋回角速度センサの出力データ
ΔθGのノイズσH1σGを記憶しておき、オフセット
0(t)を、移動体停止時の実測オフセット値と、停止
から次の停止までの期間とに基づき具体化する。
より、上記ランダムノイズ成分v、Wとして、地磁気セ
ンサの出力データθH1旋回角速度センサの出力データ
ΔθGのノイズσH1σGを記憶しておき、オフセット
0(t)を、移動体停止時の実測オフセット値と、停止
から次の停止までの期間とに基づき具体化する。
すなわち、オフセット値は移動体が停止したときに実測
することができるので、上記関数□(1)は1−0で一
定値A1とおける。また、停止から次の停止までの期間
をTとし、今回の停止時の実測オフセット値をA2とす
るとo (T)−A2とおける。このことは、関数0(
t)に境界条件を与えることを意味する。関数形0(t
)を理論または実測により仮定できれば、これにより、
実際に用いられている旋回角速度センサのオフセット特
性を正確に反映するオフセット関数0(t)を求めるこ
とができる。
することができるので、上記関数□(1)は1−0で一
定値A1とおける。また、停止から次の停止までの期間
をTとし、今回の停止時の実測オフセット値をA2とす
るとo (T)−A2とおける。このことは、関数0(
t)に境界条件を与えることを意味する。関数形0(t
)を理論または実測により仮定できれば、これにより、
実際に用いられている旋回角速度センサのオフセット特
性を正確に反映するオフセット関数0(t)を求めるこ
とができる。
そして、式(σG+o(t))2により見掛けのノイズ
成分を計算する。これにより、各センサの出力データθ
Hの信頼性とともに、ΔθGの信頼性をオフセット分を
含めて評価することができる。沖!定されたノイズの値
は、フィルタゲイン算出手段13に供給され、カルマン
フィルタゲインKが算出される。
成分を計算する。これにより、各センサの出力データθ
Hの信頼性とともに、ΔθGの信頼性をオフセット分を
含めて評価することができる。沖!定されたノイズの値
は、フィルタゲイン算出手段13に供給され、カルマン
フィルタゲインKが算出される。
そして、カルマンフィルタ手段11において、所定の重
み付け処理が施され推定方位θが出力される。
み付け処理が施され推定方位θが出力される。
このようにして、ΔθGについて、オフセット値の変化
に応じた重み付けを行うことができるので、旋回角速度
センサのオフセット特性に応じた、正確な移動体の方位
を推定することができる。
に応じた重み付けを行うことができるので、旋回角速度
センサのオフセット特性に応じた、正確な移動体の方位
を推定することができる。
また、ノイズ算出手段がノイズ成分σHσGを走行デー
タに基づき収集するものであれば、実際の走行条件に応
じたσ8 σ6の評価をリアルタイムで行うことができ
、推定方位の信頼性をさらに向上させることができる。
タに基づき収集するものであれば、実際の走行条件に応
じたσ8 σ6の評価をリアルタイムで行うことができ
、推定方位の信頼性をさらに向上させることができる。
なお、上記走行中のノイズ成分を測定する理由は、走行
中の振動や、移動体の動特性等に起因するノイズを含ん
でいる走行中のデータによってノイズを評価することが
好ましいからである。
中の振動や、移動体の動特性等に起因するノイズを含ん
でいる走行中のデータによってノイズを評価することが
好ましいからである。
特に、旋回による角速度成分を含まない直線走行中のデ
ータを利用することがより好ましい。
ータを利用することがより好ましい。
〈実施例〉
以下実施例を示す添付図面によって詳細に説明する。
第3図は、車両に適用された本発明の方位検出装置の一
実施例を示すブロック図であり、・左右両輪の回転数を
検出する車速センサ41(このセンサは、距離センサと
して利用される。)・地磁気センサ42、 ・ジャイロ43(旋回角速度を干渉光の位相変化として
読み取る光フアイバジャイロ、ピエゾエレクトリック素
子の片持ちぼり振動技術を利用して旋回角速度を検出す
る振動ジャイロ、機械式ジャイロ等から選ばれたもの。
実施例を示すブロック図であり、・左右両輪の回転数を
検出する車速センサ41(このセンサは、距離センサと
して利用される。)・地磁気センサ42、 ・ジャイロ43(旋回角速度を干渉光の位相変化として
読み取る光フアイバジャイロ、ピエゾエレクトリック素
子の片持ちぼり振動技術を利用して旋回角速度を検出す
る振動ジャイロ、機械式ジャイロ等から選ばれたもの。
旋回角速度センサとして利用される。)、
・道路地図データを格納した道路地図メモリ2、・ジャ
イロ43、地磁気センサ42により検出された出力デー
タに基づいて車両の推定方位を算出するとともに、車速
センサ41のデータと合わせて車両の位置出力データを
演算しバッファメモリ3に格納するロケータ1、 ・当該読出した車両現在位置を地図に重ねてデイスプレ
ィ7に表示させるとともに、キーボード8とのインター
フェイスをとるナビゲーションコントローラ5 とから構成されている。
イロ43、地磁気センサ42により検出された出力デー
タに基づいて車両の推定方位を算出するとともに、車速
センサ41のデータと合わせて車両の位置出力データを
演算しバッファメモリ3に格納するロケータ1、 ・当該読出した車両現在位置を地図に重ねてデイスプレ
ィ7に表示させるとともに、キーボード8とのインター
フェイスをとるナビゲーションコントローラ5 とから構成されている。
上記ロケータ1は、例えば、車速センサ41からの出力
パルス信号の数をカウンタでカウントするこ六により車
輪の回転数を得、カウンタから出力されるカウント出力
データに対して、乗算器により1カウント当りの距離を
示す所定の定数を乗算することにより単位時間当りの走
行距離出力データを算出するとともに、ジャイロ43か
ら車両方位の相対変化を求め、これと地磁気センサ42
の絶対方位出力データとから車両の方位出力データを算
出するものである。
パルス信号の数をカウンタでカウントするこ六により車
輪の回転数を得、カウンタから出力されるカウント出力
データに対して、乗算器により1カウント当りの距離を
示す所定の定数を乗算することにより単位時間当りの走
行距離出力データを算出するとともに、ジャイロ43か
ら車両方位の相対変化を求め、これと地磁気センサ42
の絶対方位出力データとから車両の方位出力データを算
出するものである。
上記道路地図メモリ2は、所定範囲にわたる道路地図デ
ータが予め格納されているものであり、半導体メモリ、
カセットテープ、CD−ROM。
ータが予め格納されているものであり、半導体メモリ、
カセットテープ、CD−ROM。
ICメモリ、DAT等が使用可能である。
上記デイスプレィ7はCRT、液晶表示器等を使用して
、車両走行中の道路地図と車両位置とを表示するもので
ある。
、車両走行中の道路地図と車両位置とを表示するもので
ある。
上記ナビゲーション・コントローラ5は、図形処理プロ
セッサ、画像処理メモリ等から構成され、デイスプレィ
7上における地図の検索、縮尺切り替え、スクロール、
車両の現在位置の表示等を行わせる。
セッサ、画像処理メモリ等から構成され、デイスプレィ
7上における地図の検索、縮尺切り替え、スクロール、
車両の現在位置の表示等を行わせる。
上記の構成の装置による車両方位検出手順について説明
する。車両走行中は、上記ロケータ1に取り込んだ各セ
ンサ出力データに基づいて、車両の位置をデイスプレィ
7上に地図とともに表示しているが、その表示中も一定
時間ごとの割り込みにより各センサの出力データを取り
込み、車両方位を更新するようにしている。この割り込
み時の車両方位検出フローを第2図に示す。なお、この
割り込みは、車両の走行距離出力データを基にして求ま
る一定の走行距離ごとに行ってもよい。上記一定時間ま
たは一定走行距離の長さは、使用されるジャイロの種類
や地磁気センサの性能等により適宜設定される。
する。車両走行中は、上記ロケータ1に取り込んだ各セ
ンサ出力データに基づいて、車両の位置をデイスプレィ
7上に地図とともに表示しているが、その表示中も一定
時間ごとの割り込みにより各センサの出力データを取り
込み、車両方位を更新するようにしている。この割り込
み時の車両方位検出フローを第2図に示す。なお、この
割り込みは、車両の走行距離出力データを基にして求ま
る一定の走行距離ごとに行ってもよい。上記一定時間ま
たは一定走行距離の長さは、使用されるジャイロの種類
や地磁気センサの性能等により適宜設定される。
まず、ステップ■において、ジャイロ43の出カデータ
ΔθGと地磁気センサ42の出力データをθ8を取り込
む。
ΔθGと地磁気センサ42の出力データをθ8を取り込
む。
次に、車両が直線走行中かどうか判定する(ステップ■
)。直線道路を走行中かどうかは、比較的短期間のジャ
イロ43の出力データに基づいて簡単な比較計算をすれ
ば十分判断可能である。すなわちジャイロ43の出力デ
ータΔθGと前回の推定方位θとを用いて、 θ6−θ+ΔθG を計算する。この値と、例えば過去何回かの推定方位の
平均値くθ〉との差1θ6−くθ〉1が基準値eより少
なければ直線走行中と判断できる。
)。直線道路を走行中かどうかは、比較的短期間のジャ
イロ43の出力データに基づいて簡単な比較計算をすれ
ば十分判断可能である。すなわちジャイロ43の出力デ
ータΔθGと前回の推定方位θとを用いて、 θ6−θ+ΔθG を計算する。この値と、例えば過去何回かの推定方位の
平均値くθ〉との差1θ6−くθ〉1が基準値eより少
なければ直線走行中と判断できる。
このほかに、マツプマツチング手法により判断してもよ
い。
い。
直線走行中でなければ、方位検出に利用するため、上記
推定方位θGを登録して(ステップ■)、割り込み継続
回数を示すカウンタnを0とおき(ステップ■)、ステ
ップ[株]に進む。
推定方位θGを登録して(ステップ■)、割り込み継続
回数を示すカウンタnを0とおき(ステップ■)、ステ
ップ[株]に進む。
直線走行中であれば、カウンタnを1だけインクリメン
トしくステップ■)、センサ出力データΔθG θ”を
バッファメモリ3の第n番地に記憶しくステップ■;以
下、第n番地に記憶されたデータをΔθG n1θHn
で表わす)、nがバッファメモリ3の容量等から決まる
最大値mに達しているかどうか判定する(ステップ■)
。最大値mに達していれば、以下の計算をするのに十分
な数のデータが集められたと判断して、ステップ■以下
に進む。最大値mに達していなければ、割り込み前の状
態に戻り、上記の処理を繰り返す。
トしくステップ■)、センサ出力データΔθG θ”を
バッファメモリ3の第n番地に記憶しくステップ■;以
下、第n番地に記憶されたデータをΔθG n1θHn
で表わす)、nがバッファメモリ3の容量等から決まる
最大値mに達しているかどうか判定する(ステップ■)
。最大値mに達していれば、以下の計算をするのに十分
な数のデータが集められたと判断して、ステップ■以下
に進む。最大値mに達していなければ、割り込み前の状
態に戻り、上記の処理を繰り返す。
ステップ■では、過去m回で集められたデータΔθ’
n (n−ト1)を用いて平均値、−1 〈Δθ’ > = (1/m) Σ ΔθG n
。
n (n−ト1)を用いて平均値、−1 〈Δθ’ > = (1/m) Σ ΔθG n
。
n雪1
分散値、 。、。
σG 2 − (1/at) Σ (ΔθG n
−〈θG〉)2ri を求める。
−〈θG〉)2ri を求める。
ステップ■では、過去m回で集められたデータθ’ n
(n−ト1)を用いて平均値、n−m 〈θ’ > −(1/a+) Σ θHnr1 分散値、 −m σ” = (1/+) Σ (θHn <θH〉)
2n=1 を求める。
(n−ト1)を用いて平均値、n−m 〈θ’ > −(1/a+) Σ θHnr1 分散値、 −m σ” = (1/+) Σ (θHn <θH〉)
2n=1 を求める。
上記のようにして求めた分散値σG2 σH2はセンサ
の出力データの、正規分布に従うランダムノイズの大き
さに対応するものである。
の出力データの、正規分布に従うランダムノイズの大き
さに対応するものである。
次にジャイロ43の出力データΔθGのオフセットを求
める。第4図はジャイロ43の出力データの瞬時値を表
わすグラフであり、停車時t2から次の発進時t3まで
の間にオフセット値がAからOに補正され、発進後再び
オフセットo (t)が現れ、時間tとともに増加して
いく様子を示している。もし、オフセットの関数形o
(t)を時間の一次関数αtと仮定すれば、オフセット
補正値Aを前の発進時tlから停車時t2までに要した
時間T−t2−tlで割ることによりオフセットの時間
変化率αを求めることができる。
める。第4図はジャイロ43の出力データの瞬時値を表
わすグラフであり、停車時t2から次の発進時t3まで
の間にオフセット値がAからOに補正され、発進後再び
オフセットo (t)が現れ、時間tとともに増加して
いく様子を示している。もし、オフセットの関数形o
(t)を時間の一次関数αtと仮定すれば、オフセット
補正値Aを前の発進時tlから停車時t2までに要した
時間T−t2−tlで割ることによりオフセットの時間
変化率αを求めることができる。
上記時間変化率αは、停車、発進を繰り返すごとに求ま
るので、その都度αを更新していけばよい。これによっ
て、温度等によるオフセットのドリフト分を各時点に応
じて正確に評価することができる。
るので、その都度αを更新していけばよい。これによっ
て、温度等によるオフセットのドリフト分を各時点に応
じて正確に評価することができる。
上記のように関数形o (t)を時間の一次関数αtと
仮定したのは、オフセット値が第4図に示すように、時
間とともに増加することが多いので一次関数で十分近似
できると考えたからであるが、−次間数に限られるもの
でなく、旋回角速度センサの特性や実測データに応じて
最適な関数を決定すればよい。例えば、o (t)の形
として、二次関数αt2、指数関数α(e’−1)等を
仮定してもよい。
仮定したのは、オフセット値が第4図に示すように、時
間とともに増加することが多いので一次関数で十分近似
できると考えたからであるが、−次間数に限られるもの
でなく、旋回角速度センサの特性や実測データに応じて
最適な関数を決定すればよい。例えば、o (t)の形
として、二次関数αt2、指数関数α(e’−1)等を
仮定してもよい。
そこで、前回求めたオフセットの時間変化率αを用いて
、次式 %式%) により見掛けの分散値を計算する。
、次式 %式%) により見掛けの分散値を計算する。
なおこの実施例では、データΔθGのサンプリングは時
間Δtごととしており、ステップ[株]において、Δt
ごとに増加するオフセット増加率α′−α/Δtを用い
て (σ6+α′n)2 で計算している。以後この見掛けの分散値をジャイロ4
3の出力データΔθGの分散値と見なすことにし、改め
てσG2とおく。
間Δtごととしており、ステップ[株]において、Δt
ごとに増加するオフセット増加率α′−α/Δtを用い
て (σ6+α′n)2 で計算している。以後この見掛けの分散値をジャイロ4
3の出力データΔθGの分散値と見なすことにし、改め
てσG2とおく。
ステップ■では分散値σ62.σH2を考慮にいれた推
定方位を次式に基づいて計算する。
定方位を次式に基づいて計算する。
θi −(1−Ki ) (θi−1+Δθ’i)+
Ki θH1 ここに、θiは今回の割り込みで求める推定方位、θi
−1は前回に求めた方位である。ΔθGi。
Ki θH1 ここに、θiは今回の割り込みで求める推定方位、θi
−1は前回に求めた方位である。ΔθGi。
θHiは今回推定方位を計算するときに使用するセンサ
の出力データであり、例えば最新のデータΔθG m、
θHrt、を使用してもよいし上記平均値くΔθ’><
θ”〉を使用してもよい。K1は0<Ki<1なる変数
であって、カルマンゲインと呼ばれる。Klは前回に求
めたカルマンゲインK i−1を使って、 σG 2 +σG2+K1−1 σH2として求めら
れる。
の出力データであり、例えば最新のデータΔθG m、
θHrt、を使用してもよいし上記平均値くΔθ’><
θ”〉を使用してもよい。K1は0<Ki<1なる変数
であって、カルマンゲインと呼ばれる。Klは前回に求
めたカルマンゲインK i−1を使って、 σG 2 +σG2+K1−1 σH2として求めら
れる。
以上のようにして、特にジャイロ43のオフセットを考
慮に入れたノイズ量に応じてカルマンゲインKiをリア
ルタイムで求め、重み付けをした結果として、新しい車
両の方位θlを求め(ステップ@)、初期状態に戻る(
ステップ@)。
慮に入れたノイズ量に応じてカルマンゲインKiをリア
ルタイムで求め、重み付けをした結果として、新しい車
両の方位θlを求め(ステップ@)、初期状態に戻る(
ステップ@)。
この方位θlと車速センサ41の距離データとから車両
の推定位置を算出することができる。
の推定位置を算出することができる。
勿論この時に道路地図データと比較し、道路地図データ
との相関度を評価して車両の推定位置を補正し、車両の
現在位置を道路上に設定するマツプマツチング方式を採
用してもよい(特開昭63−148115号公報参照)
。
との相関度を評価して車両の推定位置を補正し、車両の
現在位置を道路上に設定するマツプマツチング方式を採
用してもよい(特開昭63−148115号公報参照)
。
以上、実施例に基づいて本発明の方位検出装置を説明し
てきたが、本発明は上記実施例に限るものではない。例
えば、上記実施例のステップ■■では、過去の走行で集
められたデータΔθG n。
てきたが、本発明は上記実施例に限るものではない。例
えば、上記実施例のステップ■■では、過去の走行で集
められたデータΔθG n。
θHnを用いて平均値、分散値を求めていたが、旋回角
速度センサ、地磁気センサの特性に応じた一定値を予め
記憶しておき、その値のみを用いて計算してもよい。こ
うすることによって、制御が簡単になるとともに、直線
走行中でない場合でも正確な方位を推定できるようにな
る。
速度センサ、地磁気センサの特性に応じた一定値を予め
記憶しておき、その値のみを用いて計算してもよい。こ
うすることによって、制御が簡単になるとともに、直線
走行中でない場合でも正確な方位を推定できるようにな
る。
また、車両の代わりに航空機、船舶等の移動体に適用す
ることができ、その他、本発明の要旨を変更しない範囲
内において、種々の設計変更を施すことが可能である。
ることができ、その他、本発明の要旨を変更しない範囲
内において、種々の設計変更を施すことが可能である。
〈発明の効果〉
以上のように、本発明の方位検出装置によれば、旋回角
速度センサのオフセット値を、移動体停止時の実測オフ
セット値と、停止時から停止時までに計測した期間とに
基づき、時間とともに変化する関数o (t)として予
測し、地磁気センサの出力データθHおよび旋回角速度
センサの出力データΔθGに含まれるノイズ成分と、上
記予測されたオフセットo (t)とを用いて、旋回角
速度センサの実際のオフセット特性を反映した移動体の
正確な方位を推定することができ、上記方位に基づいて
移動体の移動中の位置を正確に求めることができる。
速度センサのオフセット値を、移動体停止時の実測オフ
セット値と、停止時から停止時までに計測した期間とに
基づき、時間とともに変化する関数o (t)として予
測し、地磁気センサの出力データθHおよび旋回角速度
センサの出力データΔθGに含まれるノイズ成分と、上
記予測されたオフセットo (t)とを用いて、旋回角
速度センサの実際のオフセット特性を反映した移動体の
正確な方位を推定することができ、上記方位に基づいて
移動体の移動中の位置を正確に求めることができる。
また、地磁気センサの出力データθHおよび旋回角速度
センサの出力データΔθGに含まれるノイズ成分を実走
行データに基づきリアルタイムで求めることにすれば、
より正確な方位を推定することができる。
センサの出力データΔθGに含まれるノイズ成分を実走
行データに基づきリアルタイムで求めることにすれば、
より正確な方位を推定することができる。
第1図は本発明の方位検出装置の機能ブロック図、
第2図は方位検出手順を示すフローチャート、第3図は
方位検出装置のハードウェア構成を示すブロック図、 第4図は旋回角速度センサのオフセット値の時間変化を
示すグラフである。 1・・・ロケータ、3・・・バッファメモリ、11・・
・方位推定手段、12・・・ノイズ算出手段、13・・
・フィルタゲイン算出手段、 1 5・・・オフセラ ト予測手段、 2・・・地磁気センサ、 3・・・ジャイロ
方位検出装置のハードウェア構成を示すブロック図、 第4図は旋回角速度センサのオフセット値の時間変化を
示すグラフである。 1・・・ロケータ、3・・・バッファメモリ、11・・
・方位推定手段、12・・・ノイズ算出手段、13・・
・フィルタゲイン算出手段、 1 5・・・オフセラ ト予測手段、 2・・・地磁気センサ、 3・・・ジャイロ
Claims (1)
- 【特許請求の範囲】 1、地磁気センサの方位データ、および旋回角速度セン
サの出力から求められる方位データを取り込み、それら
の値と過去の推定方位とから移動体の現在の推定方位を
求める方位検出装置において、 旋回角速度センサのオフセット値を、移動 体停止時の実測オフセット値と、停止時から次の停止時
までに計測した期間とに基づき、時間の関数o(t)と
して予測するオフセット値予測手段と、 地磁気センサの出力データθ^Hのノイズ成分σ^Hを
記憶するとともに、旋回角速度センサの出力データΔθ
^Gのノイズ成分σ^Gを記憶し、上記予測されたオフ
セットo(t)を用いて旋回角速度センサの出力データ
Δθ^Gの見掛けのノイズ成分を次式により計算するノ
イズ算出手段と、 (σ^G+o(t))^2 上記ノイズ成分を用いてカルマンフィルタ ゲインKを算出するフィルタゲイン算出手段と、 地磁気センサの出力から求めた方位データ θ^H、および旋回角速度センサの出力から求めた方位
データθ^Gに、上記カルマンフィルタゲインKに基づ
いた重み付け処理をすることにより、移動体の現在の推
定方位θを求める方位推定手段とを有することを特徴と
する方位検出装置。 2、ノイズ算出手段が、地磁気センサの出力データθ^
Hに含まれるノイズ成分σ^H、および旋回角速度セン
サの出力データΔθ^Gに含まれるノイズ成分σ^Gを
実際の走行データに基づいて算出するものである請求項
1記載の方位検出装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP32985389A JPH03188318A (ja) | 1989-12-19 | 1989-12-19 | 方位検出装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP32985389A JPH03188318A (ja) | 1989-12-19 | 1989-12-19 | 方位検出装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH03188318A true JPH03188318A (ja) | 1991-08-16 |
Family
ID=18225967
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP32985389A Pending JPH03188318A (ja) | 1989-12-19 | 1989-12-19 | 方位検出装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH03188318A (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002333321A (ja) * | 2001-05-07 | 2002-11-22 | Pioneer Electronic Corp | 角速度検出方法及び装置、角度検出方法及び装置、ナビゲーションシステム及びコンピュータプログラム |
-
1989
- 1989-12-19 JP JP32985389A patent/JPH03188318A/ja active Pending
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002333321A (ja) * | 2001-05-07 | 2002-11-22 | Pioneer Electronic Corp | 角速度検出方法及び装置、角度検出方法及び装置、ナビゲーションシステム及びコンピュータプログラム |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5319561A (en) | Apparatus for estimating current heading using magnetic and angular velocity sensors | |
| US5235514A (en) | Apparatus for estimating current heading using magnetic and angular velocity sensors | |
| US5469158A (en) | Apparatus for correcting the detected heading of a vehicle | |
| EP0496538B1 (en) | Vehicle heading correction apparatus | |
| JP3018497B2 (ja) | 旋回角速度センサのオフセット補正装置 | |
| JPH04235310A (ja) | ジャイロのバイアス補正装置 | |
| JPH0731062B2 (ja) | ジャイロのオフセット補正方法及び装置 | |
| JPH02266221A (ja) | 自立航法装置に用いる角速度センサの較正装置 | |
| EP1315945A2 (en) | Calibration of multi-axis accelerometer in vehicle navigation system using gps data | |
| JPH04238216A (ja) | ジャイロのスケールファクタの算出方法 | |
| JPH02501855A (ja) | 電子コンパスを備えた車両用のナビゲーシヨン方法 | |
| US7032450B2 (en) | Method and apparatus for measuring speed of land vehicle using accelerometer | |
| JPH0833302B2 (ja) | 位置検出装置 | |
| JP2734150B2 (ja) | オフセット補正装置 | |
| JP2514766B2 (ja) | 車両方位補正装置 | |
| JPH03188317A (ja) | 方位検出装置 | |
| JPH03188318A (ja) | 方位検出装置 | |
| JPH03154820A (ja) | 車両位置検出装置 | |
| JPH03188316A (ja) | 方位検出装置 | |
| JP2723352B2 (ja) | 車載用ナビゲーションシステム | |
| JPH06288776A (ja) | 方位検出装置 | |
| JPH03221812A (ja) | 方位検出装置 | |
| JPH03154819A (ja) | 車両位置検出装置 | |
| JPH03206913A (ja) | 方位検出装置 | |
| JPH03221811A (ja) | オフセット補正機能を有する方位検出装置 |