(57)【特許請求の範囲】
【請求項１】
車両に搭載される車載装置であって、
センサーからの出力に基づき、前記車両の変動に関する観測量を観測する観測部と、
カルマンフィルタにより、前記車両の状態を示す状態量を推定する推定部と、を備え、
前記推定部は、
前記車両の前記状態量の予測値を算出し、
前記観測量と微積分の関係にある前記状態量の誤差として前記観測量の誤差が入力された前記カルマンフィルタにより、前記予測値の誤差を算出し、
算出した前記予測値と前記予測値の誤差とに基づいて、前記カルマンフィルタにより前記車両の前記状態量の推定値と前記推定値の誤差とを算出する、
ことを特徴とする車載装置。
【請求項２】
前記推定部は、
前回算出した前記状態量の前記推定値の誤差と、前記状態量の誤差として入力された前記観測量の誤差との共分散が入力された前記カルマンフィルタにより、前記予測値の誤差を算出する、
ことを特徴とする請求項１に記載の車載装置。
【請求項３】
前記観測部は、前記車両の速度を検出する車速センサーからの出力に基づき、前記車両の速度を観測し、
前記推定部は、前記車両の位置の誤差として前記車両の速度の誤差が入力された前記カルマンフィルタにより、前記予測値の誤差を算出する、
ことを特徴とする請求項１又は２に記載の車載装置。
【請求項４】
前記観測部は、前記車両の角速度を検出する角速度センサーからの出力に基づき、前記車両の角速度を観測し、
前記推定部は、前記車両の方位の誤差として前記車両の角速度の誤差が入力された前記カルマンフィルタにより、前記予測値の誤差を算出する、
ことを特徴とする請求項１又は２に記載の車載装置。
【請求項５】
前記観測部は、前記車両の加速度を検出する加速度センサーからの出力に基づき、前記車両の加速度を観測し、
前記推定部は、前記車両の速度の誤差として前記車両の加速度の誤差が入力された前記カルマンフィルタにより、前記予測値の誤差を算出する、
ことを特徴とする請求項１又は２に記載の車載装置。
【請求項６】
車両の状態を示す状態量をカルマンフィルタにより推定する推定方法であって、
前記車両に搭載される車載装置が、
前記車両の前記状態量の予測値を算出し、
センサーからの出力に基づき観測された前記車両の変動に関する観測量と微積分の関係にある前記状態量の誤差として、前記観測量の誤差が入力された前記カルマンフィルタにより、前記予測値の誤差を算出し、
算出した前記予測値と前記予測値の誤差とに基づいて、前記カルマンフィルタにより前記車両の前記状態量の推定値と前記推定値の誤差とを算出する、
ことを特徴とする推定方法。
【請求項７】
前記車載装置が、
前回算出した前記状態量の前記推定値の誤差と、前記状態量の誤差として入力された前記観測量の誤差との共分散が入力された前記カルマンフィルタにより、前記予測値の誤差を算出する、
ことを特徴とする請求項６に記載の推定方法。
【請求項８】
前記車載装置が、
前記車両の速度を検出する車速センサーからの出力に基づき、前記車両の速度を観測し、
前記車両の位置の誤差として前記車両の速度の誤差が入力された前記カルマンフィルタにより、前記予測値の誤差を算出する、
ことを特徴とする請求項６又は７に記載の推定方法。
【請求項９】
前記車載装置が、
前記車両の角速度を検出する角速度センサーからの出力に基づき、前記車両の角速度を観測し、
前記車両の方位の誤差として前記車両の角速度の誤差が入力された前記カルマンフィルタにより、前記予測値の誤差を算出する、
ことを特徴とする請求項６又は７に記載の推定方法。
【請求項１０】
前記車載装置が、
前記車両の加速度を検出する加速度センサーからの出力に基づき、前記車両の加速度を観測し、
前記車両の速度の誤差として前記車両の加速度の誤差が入力された前記カルマンフィルタにより、前記予測値の誤差を算出する、
ことを特徴とする請求項６又は７に記載の推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、車載装置、及び、推定方法に関する。
【背景技術】
【０００２】
本技術分野の背景技術として、特開２０００−５５６７８号公報（特許文献１）がある。特許文献１には、「車速センサ４、ジャイロ６での検出値から移動距離演算部１１、方位変化量演算部１２にて算出される移動距離、方位変化量に基づき、相対軌跡演算部１３、絶対位置演算部１４が算出する推定航法データ（車速、絶対方位、絶対位置）と、ＧＰＳ受信機８からのＧＰＳ測位データ（速度、方位、位置）との差を観測値とするカルマンフィルタ（誤差推定部１５）を備え、その状態量の一つとして、ジャイロ出力から角速度への変換ゲインの誤差（ゲイン誤差）を設定する。このカルマンフィルタにより求められるゲイン誤差の推定値により、ジャイロ出力及び変換ゲインを用いて算出された方位変化量を補正する。」と記載されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０００−５５６７８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ところで、車速パルスやジャイロ等のセンサーからの出力に基づき観測される速度や角速度等の観測量から求められる移動距離や方位変化量は、時間の経過により誤差が累積する場合がある。特許文献１では、累積する誤差を考慮したカルマンフィルタにより車両の状態を推定する技術について開示が無く、カルマンフィルタにより車両の状態を精度よく推定できない。
そこで、本発明は、カルマンフィルタにより車両の状態を精度よく推定できるようにすることを目的とする。
【課題を解決するための手段】
【０００５】
上記目的を達成するために、車載装置は、車両に搭載される車載装置であって、センサーからの出力に基づき、前記車両の変動に関する観測量を観測する観測部と、カルマンフィルタにより、前記車両の状態を示す状態量を推定する制御部と、を備え、前記制御部は、前記車両の前記状態量の予測値を算出し、前記観測量と微積分の関係にある前記状態量の誤差として前記観測量の誤差が入力された前記カルマンフィルタにより、前記予測値の誤差を算出し、算出した前記予測値と前記予測値の誤差とに基づいて、前記カルマンフィルタにより前記車両の前記状態量の推定値と前記推定値の誤差とを算出することを特徴とする。
【発明の効果】
【０００６】
本発明によれば、カルマンフィルタにより車両の状態を精度よく推定できる。
【図面の簡単な説明】
【０００７】
【図１】ナビゲーション装置の構成を示すブロック図である。
【図２】ナビゲーション装置の動作を示すフローチャートである。
【図３】車両状態推定処理における推定部の動作を示すフローチャートである。
【図４】推定した車両の状態に基づくマップマッチングの処理を説明するための図である。
【発明を実施するための形態】
【０００８】
図１は、ナビゲーション装置１（車載装置）の構成を示すブロック図である。
【０００９】
ナビゲーション装置１は、車両に搭載される車載装置であり、車両に搭乗しているユーザの操作に従って、地図の表示や、地図における車両の現在位置の表示、目的地までの経路探索、経路案内等を実行する。なお、ナビゲーション装置１は、車両のダッシュボード等に固定されてもよいし、車両に対し着脱可能なものであっても良い。
【００１０】
図１に示すように、ナビゲーション装置１は、制御部２と、記憶部３と、タッチパネル４と、ＧＰＳ受信部５と、車速センサー６（センサー）と、ジャイロセンサー７（角速度センサー、センサー）と、加速度センサー８（センサー）とを備える。
【００１１】
制御部２は、ＣＰＵや、ＲＯＭ、ＲＡＭ、その他の制御回路等を備え、ナビゲーション装置１の各部を制御する。制御部２は、ＲＯＭや記憶部３等に記憶された制御プログラムを実行することにより、後述する観測部２１、及び、推定部２２として機能する。
【００１２】
記憶部３は、ハードディスクや、ＥＥＰＲＯＭ等の不揮発性メモリを備え、データを書き換え可能に記憶する。記憶部３は、制御部２が実行する制御プログラムのほか、地図データ３ａを記憶する。地図データ３ａは、交差点やその他の道路網上の接続点を示すノードに関するノード情報や、ノードとノードとの道路区間を示すリンクに関するリンク情報、地図の表示に係る情報等を有する。リンク情報は、リンク毎に、少なくとも、リンクの位置に関する情報と、リンクの方位に関する情報とを含む。
【００１３】
タッチパネル４は、表示パネル４ａと、タッチセンサー４ｂとを備える。表示パネル４ａは、液晶ディスプレイやＥＬ（Electro Luminescent）ディスプレイ等により構成され、制御部２の制御の下、各種情報を表示パネル４ａに表示する。タッチセンサー４ｂは、表示パネル４ａに重ねて配置され、ユーザのタッチ操作を検出し、タッチ操作された位置を示す信号を制御部２に出力する。制御部２は、タッチセンサー４ｂからの入力に基づいて、タッチ操作に対応する処理を実行する。
【００１４】
ＧＰＳ受信部５は、ＧＰＳアンテナ５ａを介してＧＰＳ衛星からのＧＰＳ電波を受信し、ＧＰＳ電波に重畳されたＧＰＳ信号から、車両の位置と、車両の進行方向の方角（以下、「車両の方位」と表現する）とを演算により取得する。ＧＰＳ受信部５は、車両の位置を示す情報、及び、車両の方位を示す情報を制御部２に出力する。
【００１５】
車速センサー６は、車両の車速を検出して、検出した車速を示す信号を制御部２に出力する。
【００１６】
ジャイロセンサー７は、例えば振動ジャイロにより構成され、車両の回転による角速度を検出する。ジャイロセンサー７は、検出した角速度を示す信号を制御部２に出力する。
【００１７】
加速度センサー８は、車両に作用する加速度（例えば、進行方向に対する車両の傾き）を検出する。加速度センサー８は、検出した加速度を示す信号を制御部２に出力する。
【００１８】
図１に示すように、制御部２は、観測部２１と、推定部２２とを備える。
【００１９】
観測部２１は、車速センサー６、ジャイロセンサー７、及び、加速度センサー８から出力される信号に基づき、車両の変動に関する観測量を観測する。本実施形態では、観測部２１は、車速センサー６から出力される車速を示す信号に基づき、観測量として車両の速度（観測量）を観測する。観測部２１は、所定の演算により、車速を示す信号から、車両の速度と車両の速度の誤差とを算出する。この算出された車両の速度の誤差は、この速度から求めた移動距離に累積する誤差である。
【００２０】
また、観測部２１は、ジャイロセンサー７から出力される角速度を示す信号に基づき、観測量として車両の回転による角速度（観測量）を観測する。観測部２１は、所定の演算により、角速度を示す信号から、車両の角速度と車両の角速度の誤差とを算出する。この算出された車両の角速度の誤差は、この角速度から求めた方位変化量に累積する誤差である。
【００２１】
また、観測部２１は、加速度センサー８から出力される加速度（観測量）を示す信号に基づき、観測量として車両の加速度を観測する。観測部２１は、所定の演算により、加速度を示す信号から、車両の加速度と車両の加速度の誤差とを算出する。この算出された車両の加速度の誤差は、この加速度から求めた速度に累積する誤差である。
【００２２】
また、観測部２１は、ＧＰＳ受信部５から出力された車両の位置を示す情報と車両の方位を示す情報とに基づき、車両の位置と車両の方位とを観測する。
【００２３】
推定部２２は、観測部２１が観測する車両の速度、車両の加速度、車両の角速度、車両の位置、及び、車両の方位に基づき、カルマンフィルタにより、車両の状態を示す状態量を推定する。本実施形態では、推定部２２は、車両の位置、車両の方位、車両の速度、及び、車両の角速度を、車両の状態量として推定する。
【００２４】
後述するが、制御部２は、推定部２２が推定した車両の状態量に基づき、マップマッチングの対象となるリンクの評価を実行する。
【００２５】
なお、推定部２２が推定する車両の速度は、車両の状態量に相当する。また、観測部２１が観測する車両の速度は、観測量に相当する。同様に、推定部２２が推定する車両の角速度は、車両の状態量に相当する。また、観測部２１が観測する車両の角速度は、観測量に相当する。
【００２６】
ここで、カルマンフィルタによる基本的な車両の状態量の推定について説明する。
【００２７】
本実施形態では、カルマンフィルタにより推定する車両の状態量は、車両の位置、車両の方位、車両の速度、及び、車両の角速度である。以下に、推定する車両の各状態量を示す。
【００２８】
ｘ：車両の位置のｘ座標
ｙ：車両の位置のｙ座標
θ：車両の方位
ｖ：車両の速度
ω：車両の角速度
【００２９】
ここで、車両の状態量をベクトル表記した状態ベクトルを、（ｘ、ｙ、θ、ｖ、ω）とすると、車両の状態量についての状態方程式は、式（１）で表される。
【００３０】
【数１】
0006677533.tif 000002
【００３１】
添え字のｋは、時刻を示す。例えば、式（１）の左辺である（ｘｋ＋１、ｙｋ＋１、θｋ＋１、ｖｋ＋１、ωｋ＋１）は、時刻ｋ＋１における車両の状態量を示す。式（１）において、右辺の第２項であるｑｋは、システム雑音（平均０、誤差共分散行列であるＱｋを有する正規分布Ｎ（０、Ｑｋ））である。誤差共分散行列とは、分散と共分散との行列である。
【００３２】
本実施形態では、観測部２１は、観測対象として、車両の速度、車両の角速度、車両の加速度、車両の位置、及び、車両の方位を観測する。前述した通り、観測部２１は、車両の速度について、車速センサー６からの出力に基づき観測する。また、観測部２１は、車両の角速度について、ジャイロセンサー７からの出力に基づき観測する。また、観測部２１は、車両の加速度について、加速度センサー８からの出力に基づき観測する。また、観測部２１は、車両の位置、及び、車両の方位について、ＧＰＳ受信部５からの出力に基づき観測する。以下に、観測部２１が観測する観測対象を示す。なお、以下では、観測対象として車両の速度、車両の角速度、車両の位置、及び、車両の方位を例示する。
【００３３】
ｖＰＬＳ：車速センサー６からの出力に基づき観測される車両の速度
ωＧＹＲ：ジャイロセンサー７からの出力に基づき観測される車両の角速度
ｘＧＰＳ：ＧＰＳ受信部５からの出力に基づき観測される車両の位置のｘ座標
ｙＧＰＳ：ＧＰＳ受信部５からの出力に基づき観測される車両の位置のｙ座標
θＧＰＳ：ＧＰＳ受信部５からの出力に基づき観測される車両の方位
【００３４】
ここで、上記の観測対象をベクトル表記した観測ベクトルを（ｖＰＬＳ、ωＧＹＲ、ｘＧＰＳ、ｙＧＰＳ、θＧＰＳ）とすると、観測対象についての観測方程式は、式（２）で表される。
【００３５】
【数２】
0006677533.tif 000003
【００３６】
式（２）において、ｒｋは、観測雑音（平均０、誤差共分散行列であるＲｋを有する正規分布Ｎ（０、Ｒｋ））である。
【００３７】
以下、カルマンフィルタについて、車両の状態量を予測する予測処理と、車両の状態量を推定する推定処理とに分けて説明する。
【００３８】
なお、以下の説明において、添え字のｋ＋１｜ｋが付与された値は、時刻ｋまでの情報に基づき予測された時刻ｋ＋１における予測値を示す。また、添え字のｋ｜ｋ−１が付与された値は、時刻ｋ−１までの情報に基づき予測された時刻ｋにおける予測値を示す。また、添え字のｋ＋１｜ｋ＋１が付与された値は、時刻ｋ＋１までの情報に基づき推定された時刻ｋ＋１における推定値を示す。また、添え字のｋ｜ｋが付与された値は、時刻ｋまでの情報に基づき推定された時刻ｋにおける推定値を示す。また、添え字のｋ−１｜ｋ−１が付与された値は、時刻ｋ−１までの情報に基づき推定された時刻ｋ−１における推定値を示す。
【００３９】
＜予測処理＞
カルマンフィルタにおいて予測処理は、車両の状態量の予測値（以下、「車両状態予測値」と表現する）と、車両状態予測値の誤差共分散行列と、を算出する処理である。車両状態予測値は、式（３）に基づき算出される。なお、誤差共分散行列の算出とは、誤差共分散行列の各成分の値の算出を示す。
【００４０】
【数３】
0006677533.tif 000004
【００４１】
式（３）では、時刻ｋまでの情報に基づき予測された時刻ｋ＋１における車両状態予測値の算出を示す。この車両状態予測値は、式（３）の右辺に示すように、時刻ｋまでの情報に基づき推定された時刻ｋにおける車両の状態量の推定値により算出される。例えば、車両の位置のｘ座標の予測値を示すｘｋ＋１｜ｋは、時刻ｋまでの情報に基づき推定された時刻ｋにおける、車両の位置のｘ座標（ｘｋ｜ｋ）の推定値と、車両の方位（θｋ｜ｋ）の推定値と、車両の速度（ｖｋ｜ｋ）の推定値とに基づき算出される。なお、式（３）において、Ｔは、車速センサー６、ジャイロセンサー７、及び、加速度センサー８からの出力に基づき観測部２１が観測量を観測する間隔を示す。
【００４２】
なお、時刻ｋ−１までの情報に基づき予測された時刻ｋの車両状態予測値を算出する場合、式（３）において時刻を１ステップずつ下げた式により算出できる。つまり、時刻ｋ−１までの情報に基づき予測された時刻ｋにおける車両状態予測値は、時刻ｋ−１までの情報に基づき推定された時刻ｋ−１における車両の状態量の推定値に基づき算出される。
【００４３】
車両状態予測値の誤差共分散行列は、式（４）に基づき算出される。誤差共分散行列は、本実施形態において、車両の状態量に関する分散及び共分散の行列である。分散は、誤差を二乗したものである。つまり、車両の状態量の分散は、車両の状態量の誤差を二乗したものである。したがって、車両状態予測値の誤差共分散行列を算出することは、車両状態予測値の誤差を算出することに相当する。
【００４４】
【数４】
0006677533.tif 000005
【００４５】
式（４）において、Ｐは、誤差共分散行列を示す。式（４）の左辺は、時刻ｋまでの情報に基づき予測された時刻ｋ＋１における誤差共分散行列を示す。式（４）の左辺に示す誤差共分散行列は、時刻ｋまでの情報に基づき推定された時刻ｋにおける誤差共分散行列に基づき算出される。なお、Ｆは、式（１）の状態方程式から求められるヤコビ行列を示す。また、Ｆにおける添え字のＴは、転置行列を示す。また、Ｑは、システム雑音の誤差共分散行列を示す。
【００４６】
なお、時刻ｋ−１までの情報に基づき予測された時刻ｋにおける車両状態予測値の誤差共分散行列を算出する場合、式（４）においてＰの添え字の時刻を１ステップずつ下げた式により算出できる。つまり、時刻ｋ−１までの情報に基づき予測された時刻ｋにおける車両状態予測値の誤差共分散行列は、時刻ｋ−１までの情報に基づき推定された時刻ｋ−１における車両の状態量の推定値の誤差共分散行列により算出される。
【００４７】
このように、予測処理では、例えば、時刻ｋにおける車両の状態量を予測する場合、時刻ｋ−１までの情報に基づき算出された時刻ｋにおける車両状態予測値と、時刻ｋ−１までの情報に基づき算出された時刻ｋにおける車両状態予測値の誤差共分散行列とを算出する。つまり、予測処理は、車両の状態量の確立分布を予測している。
【００４８】
＜推定処理＞
次に、推定処理について説明する。
カルマンフィルタにおいて推定処理は、予測処理にて算出した車両状態予測値、及び、車両状態予測値の誤差共分散行列に基づいて、車両の状態量の推定値（以下、「車両状態推定値」と表現する）と、車両状態推定値の誤差共分散行列とを算出する処理である。
【００４９】
推定処理では、式（５）により観測残差が算出される。観測残差とは、観測対象の値と、車両状態予測値から算出される観測対象に対応する値との誤差である。
【００５０】
【数５】
0006677533.tif 000006
【００５１】
式（５）において、左辺は、観測残差をベクトル表記した観測残差ベクトルである。また、式（５）において、右辺の第１項は、観測部２１の観測対象とする観測対象ベクトルである。また、式（５）において、右辺の第２項は、予測処理で予測した車両状態予測値に、観測方程式から求められる観測行列であるＨを掛けたものである。
【００５２】
車両状態推定値は、式（５）に示す観測残差を用いて、式（６）により算出される。
【００５３】
【数６】
0006677533.tif 000007
【００５４】
式（６）では、時刻ｋ＋１までの情報に基づき予測された時刻ｋ＋１における車両状態推定値を示す。この車両状態推定値は、式（６）の右辺に示すように、時刻ｋまでの情報に基づき推定された時刻ｋ＋１における車両状態予測値を、観測残差より補正して算出される。なお、時刻ｋまでの情報に基づき予測された時刻ｋにおける車両状態推定値を算出する場合、時刻ｋ−１までの情報に基づき予測された時刻ｋにおける車両状態予測値を観測残差により補正することで算出される。
【００５５】
前述の通り、式（６）は、観測残差を用いて、車両状態予測値を補正し、車両状態推定値を算出する式を示す。式（６）に示すように、車両状態予測値を観測残差にて補正する際の補正係数としてＫｋが用いられる。Ｋｋは、カルマンゲインと呼ばれ、式（７）で表される。
【００５６】
【数７】
0006677533.tif 000008
【００５７】
式（７）において、Ｒｋは、観測雑音の誤差共分散行列である。また、添え字の「−１」は、逆行列を示す。式（７）に示すカルマンゲインＫｋは、時刻ｋまでの情報に基づく時刻ｋ＋１における車両状態予測値の誤差共分散行列に基づき算出される。
【００５８】
式（６）において、カルマンゲインＫｋは、車両状態推定値について、車両状態予測値を重視して算出するか、観測部２１が観測する観測対象の値を重視して算出するかを決定付ける係数である。
【００５９】
例えば、車速センサー６、ジャイロセンサー７、及び、ＧＰＳ受信部５からの出力に基づく観測対象の値の誤差が充分に小さい場合、車両状態推定値は、当該誤差が充分に小さいことから、観測対象の値となることが望ましい。これは、車両状態推定値が、誤差の充分に小さい値、つまり、精度の高い値になるためである。つまり、誤差共分散行列であるＲｋの値が充分に小さい場合、車両状態推定値は、観測対象の値となることが望ましい。ここで、Ｒｋの値が充分に小さい場合のカルマンゲインＫｋをＫｋ＝Ｈ−１として式（６）に与えると、式（６）の右辺の第１項は、無くなる。つまり、車両状態推定値は、誤差が充分に小さい観測対象の値となる。
【００６０】
また、例えば、車両状態予測値の誤差が、車速センサー６、ジャイロセンサー７、及び、ＧＰＳ受信部５からの出力に基づく観測対象の誤差より充分に小さい場合、車両状態推定値は、車両状態予測値となることが望ましい。これは、車両状態推定値が、観測対象の誤差より充分に誤差が小さい値、つまり、精度の高い値になるためである。つまり、車両状態予測値の誤差共分散行列が、Ｒｋの値より充分に小さい場合、車両状態推定値は、車両状態予測値となることが望ましい。ここで、車両状態予測値の誤差共分散行列がＲｋの値より充分に小さい場合のカルマンゲインＫｋとしてＫｋ＝０として式（６）に与えると、車両状態推定値は、車両状態予測値となる。
【００６１】
このように、カルマンゲインＫｋは、観測雑音の誤差共分散行列であるＲｋ、及び、車両状態予測値の誤差共分散行列に応じて、車両状態推定値が適切な値となるように設定する係数である。つまり、観測雑音の誤差共分散行列が観測対象の誤差共分散行列である場合、カルマンゲインＫｋは、車両状態予測値の誤差共分散行列が正確に予測されていれば、車速センサー６、ジャイロセンサー７、及び、ＧＰＳ受信部５からの出力に基づく観測対象の誤差共分散行列に応じて、車両状態推定値が適切な値となるように設定する。
【００６２】
車両状態推定値の誤差共分散行列は、式（８）により算出される。前述した通り、誤差共分散行列は、本実施形態では、車両の状態量に関する、分散及び共分散の行列である。車両の状態量の分散は、車両の状態量の誤差を二乗したものである。したがって、車両状態推定値の誤差共分散行列を算出することは、車両状態推定値の誤差を算出することに相当する。
【００６３】
【数８】
0006677533.tif 000009
【００６４】
式（８）において、Ｐは、式（４）と同様、誤差共分散行列を示す。また、式（８）において、Ｉは、単位行列を示す。式（８）の左辺は、時刻ｋ＋１までの情報に基づき推定された時刻ｋ＋１における誤差共分散行列を示す。式（８）の左辺に示す誤差共分散行列は、時刻ｋまでの情報に基づき予測された時刻ｋ＋１における誤差共分散行列に基づき算出される。
【００６５】
なお、時刻ｋまでの情報に基づき予測された時刻ｋにおける車両状態推定値の誤差共分散行列を算出する場合、式（８）においてＰの添え字の時刻を１ステップずつ下げた式により算出できる。つまり、時刻ｋまでの情報に基づき予測された時刻ｋにおける車両状態推定値の誤差共分散行列は、時刻ｋ−１までの情報に基づき推定された時刻ｋにおける車両状態予測値の誤差共分散行列により算出される。
【００６６】
式（８）は、車両状態推定値の誤差共分散行列は、車両状態予測値の誤差共分散行列に、（Ｉ−ＫｋＨ）を掛けた式である。式（８）に示すように、車両状態推定値の誤差共分散行列は、カルマンゲインＫｋの値に依存する。
【００６７】
例えば、観測雑音の誤差が充分に小さい場合のカルマンゲインＫｋとしてＫｋ＝Ｈ−１を式（８）に与えると、車両状態推定値の誤差共分散行列は零行列となる。これは、推定した車両状態推定値の誤差が充分に小さいことを示す。
【００６８】
また、例えば、観測雑音の誤差より車両状態予測値の誤差が充分に小さい場合のカルマンゲインの値としてＫｋ＝０を与えることにより、車両状態推定値の誤差共分散行列は、車両状態予測値の誤差共分散行列になる。これは、車両状態推定値の誤差が、誤差の充分に小さい車両状態予測値の誤差であることを示す。
【００６９】
このように、式（８）において、カルマンゲインＫｋは、観測雑音の誤差共分散行列と、車両状態予測値の誤差共分散行列に応じて、車両状態推定値の誤差共分散行列が適切になるように設定する。つまり、カルマンゲインＫｋは、観測雑音の精度と車両状態予測値の精度とに基づいて、車両状態推定値の誤差共分散行列が適切になるように設定する。
【００７０】
以上、推定処理では、例えば、時刻ｋにおける車両の状態量を推定する場合、時刻ｋまでの情報に基づき算出された時刻ｋにおける車両状態推定値と、時刻ｋまでの情報に基づき算出された時刻ｋにおける車両状態推定値の誤差共分散行列とを算出する。つまり、推定処理は、予測処理で予測した車両の状態量の確立分布に基づいて、車両の状態量の確立分布を推定している。
【００７１】
以上の算出により、推定部２２は、カルマンフィルタによって、車両状態推定値と、車両状態推定値の誤差共分散行列とを算出することにより、車両の状態を推定する。
【００７２】
前述した通り、カルマンゲインＫｋは、車両状態推定値と、車両状態推定値の誤差共分散行列とを適切に設定する係数である。つまり、車両の状態の推定の精度は、カルマンゲインＫｋに依存する。カルマンゲインＫｋは、前述した通り、観測雑音と、車両状態予測値の誤差共分散行列とに応じて、車両状態推定値と車両状態推定値の誤差共分散行列とを適切に設定する係数である。そのため、車両状態推定値と車両状態推定値の誤差共分散行列とを適切に算出するためには、車両状態予測値の誤差共分散行列が、正確に算出される必要がある。しかしながら、カルマンフィルタでは、システム雑音、及び、観測雑音が白色性の雑音であることが仮定される。つまり、カルマンフィルタにおいて誤差を雑音として扱う場合、誤差の平均が０であることが前提となる。そのため、平均が０とならない誤差、すなわち、累積する誤差を考慮して、車両状態予測値の誤差共分散行列を算出できない場合があり、この場合、カルマンゲインＫｋは、正確に算出されない。これは、車両状態推定値と車両状態推定値の誤差共分散行列とが精度よく算出できない、すなわち、車両の状態を精度よく推定できないことを示す。
【００７３】
そこで、本実施形態の推定部２２は、以下のカルマンフィルタにより、車両の状態を推定する。
【００７４】
以下、車両の状態を推定する際のナビゲーション装置１の動作を通して、推定部２２による車両の状態の推定について説明する。
【００７５】
図２は、ナビゲーション装置１の動作を示すフローチャートである。
【００７６】
ナビゲーション装置１の観測部２１は、車速センサー６、ジャイロセンサー７、及び、加速度センサー８から出力される信号に基づき、車両の速度、車両の角速度、及び、車両の加速度を観測する（ステップＳＡ１）。観測部２１は、車速センサー６、ジャイロセンサー７、及び、加速度センサー８が信号を出力する度に、車両の速度、車両の角速度、及び、車両の加速度を観測する。つまり、観測部２１がこれらを観測する間隔は、車速センサー６、ジャイロセンサー７、及び、加速度センサー８が検出する間隔と同じ間隔である。
【００７７】
次いで、観測部２１は、ＧＰＳ受信部５の出力に基づき、車両の位置、及び、車両の方位を観測する（ステップＳＡ２）。観測部２１は、車両の位置を示す情報と車両の方位を示す情報とがＧＰＳ受信部５から出力される度に、車両の位置、及び、車両の方位を観測する。つまり、観測部２１が車両の位置及び車両の方位を観測する間隔は、ＧＰＳ受信部５がＧＰＳ電波を受信する間隔となる。
【００７８】
次いで、ナビゲーション装置１の推定部２２は、観測部２１が観測した観測対象に基づき、車両の状態を推定する車両状態推定処理を実行する（ステップＳＡ３）。
【００７９】
図３は、車両状態推定処理における推定部２２の動作を示すフローチャートである。
【００８０】
推定部２２は、予測処理を実行する（ステップＳＢ１）。前述した通り、予測処理は、車両状態予測値、及び、車両状態予測値の誤差共分散行列を算出する処理である。車両状態予測値は、式（３）によって算出される。車両状態予測値の誤差共分散行列は、式（９）によって算出される。
【００８１】
【数９】
0006677533.tif 000010
【００８２】
式（９）において、Ｆｋは、式（１）の状態方程式から求められるヤコビ行列であり、式（１０）で表される。また、Ｑｋは、システム雑音の誤差共分散行列であり、式（１１）で表される。また、Ｃｋは、前回推定した車両状態推定値の誤差と、観測部２１が車速センサー６及びジャイロセンサー７からの出力に基づき観測する観測量の誤差との共分散行列であり、式（１２）で表される。Ｑｋ＋２・Ｃｋがシステム雑音である。
【００８３】
【数１０】
0006677533.tif 000011
【００８４】
【数１１】
0006677533.tif 000012
【００８５】
【数１２】
0006677533.tif 000013
【００８６】
式（１２）において、σｋｖＰＬＳは、観測部２１が算出する車両の速度の誤差を示す。また、σｋωＧＹＲは、観測部２１が算出する車両の角速度の誤差を示す。前述した通り、σｋｖＰＬＳは、時間の経過と共に、この速度から求めた移動距離に累積する誤差である。また、σｋωＧＹＲは、σｋｖＰＬＳと同様、時間と共に、この角速度から求めた方位変化量に累積する誤差である。
【００８７】
式（１２）に示すように、車両状態予測値の誤差共分散行列を算出する式（９）において、車速センサー６から出力に基づき観測部２１が算出する車両の速度の誤差は、車両の位置の誤差として入力される。すなわち、σｋｖＰＬＳは、ｃｏｓ（θｋ｜ｋ）σｋｖＰＬＳＴ、及び、ｓｉｎ（θｋ｜ｋ）σｋｖＰＬＳＴとして入力される。また、ジャイロセンサー７からの出力に基づき観測部２１が算出する車両の角速度の誤差は、車両の方位の誤差として入力される。すなわち、σｋωＧＹＲは、σｋωＧＹＲＴとして入力される。
【００８８】
このように、時間の経過と共に累積する誤差、すなわち、σｋｖＰＬＳとσｋωＧＹＲとについて、σｋｖＰＬＳを車両の位置の誤差として入力し、σｋωＧＹＲを車両の方位の誤差として入力したカルマンフィルタにより、推定部２２は、車両状態予測値の誤差共分散行列を算出する。
【００８９】
ここで、カルマンフィルタにσｋｖＰＬＳを車両の位置の誤差として入力し、σｋωＧＹＲを車両の方位の誤差として入力することについて説明する。以下では、カルマンフィルタにσｋｖＰＬＳを車両の位置の誤差として入力することについて詳述する。
【００９０】
ここで説明の便宜上、車両の状態量を、車両の位置のｘ座標と車両の速度として例示し、カルマンフィルタにσｋｖＰＬＳを車両の位置の誤差として入力することについて説明する。この例では、車両の位置のｘ座標は、車両の速度により決定する。この例では、観測対象は、車両の速度である。
【００９１】
したがって、車両の状態量についての状態方程式は、式（１３）で表される。また、観測対象についての観測方程式は、式（１４）で表される。
【００９２】
【数１３】
0006677533.tif 000014
【００９３】
【数１４】
0006677533.tif 000015
【００９４】
式（１３）及び式（１４）に基づき、時刻ｋまでの情報に基づき推定した時刻ｋにおける車両状態推定値の誤差共分散行列を、式（４）、式（６）、及び、式（８）に基づき算出する。ここで、車両状態推定値の誤差共分散行列を算出するに際し、Ｐ（ｋ｜ｋ）を式（１５）とする。Ｐ（ｋ｜ｋ）は、時刻ｋまでの情報に基づき推定した時刻ｋにおける車両状態推定値の誤差共分散行列である。
【００９５】
【数１５】
0006677533.tif 000016
【００９６】
式（１５）において、ｐ１１（ｋ｜ｋ）は、車両の位置の誤差分散を示す。また、ｐ２２（ｋ｜ｋ）は、車両の速度の誤差分散を示す。また、ｐ１２（ｋ｜ｋ）は、車両の位置と、車両の速度との誤差共分散を示す。
【００９７】
また、車両状態推定値の誤差共分散行列を算出するあたり、Ｆ、Ｑ、Ｈ、及び、Ｒのそれぞれを式（１６）〜式（１９）とする。ここで、Ｆは、状態方程式から求められるヤコビ行列である。また、Ｑは、システム雑音の分散行列である。また、Ｈは、観測方程式から求められるヤコビ行列である。また、Ｒは、観測雑音の分散行列である。観測量が車両の速度であることから、ｒｋは、車両の速度の誤差である。つまり、ｒｋ２は、観測部２１が観測する車両の速度の誤差分散を示す。
【００９８】
【数１６】
0006677533.tif 000017
【００９９】
【数１７】
0006677533.tif 000018
【０１００】
【数１８】
0006677533.tif 000019
【０１０１】
【数１９】
0006677533.tif 000020
【０１０２】
式（１５）〜式（１９）に基づき、車両状態推定値の誤差共分散行列、すなわち、Ｐ（ｋ｜Ｋ）を算出すると、車両状態推定値の誤差共分散行列は、式（２０）で表される。
【０１０３】
【数２０】
0006677533.tif 000021
【０１０４】
したがって、車両の位置の推定値の誤差分散は、式（２１）で表される。
【０１０５】
【数２１】
0006677533.tif 000022
【０１０６】
式（２１）は、共通因数であるｐ１１（ｋ｜ｋ−１）でくくると、式（２２）に変形される。ｐ１１（ｋ｜ｋ−１）は、時刻ｋ−１までの情報に基づく時刻ｋにおける車両の位置の予測値の誤差分散である。
【０１０７】
【数２２】
0006677533.tif 000023
【０１０８】
式（２２）に示すように、車両の位置の推定値の誤差分散であるｐ１１（ｋ｜ｋ）は、車両の位置の予測値の誤差分散であるｐ１１（ｋ｜ｋ−１）に、式（２３）を掛けることで求められる。
【０１０９】
【数２３】
0006677533.tif 000024
【０１１０】
式（２３）は、式（２２）に示すように、車両の位置の予測値の分散であるｐ１１（ｋ｜ｋ−１）に掛けられる係数である。この係数は、分母に車両の速度の誤差（ｒｋ）が含まれることを除くと、式（２４）と一致する。
【０１１１】
【数２４】
0006677533.tif 000025
【０１１２】
式（２４）に示す係数は、「１−（車両の位置の誤差と車両の速度の誤差との相関係数）２」である。これは、式（２４）に示す係数が、最小二乗法を用いた線形回帰モデルの非決定係数であることを示す。つまり、車両の位置の誤差と、車両の速度の誤差との相関性が小さければ小さいほど、当該係数は大きい値となり、車両の位置の予測値の分散は大きい値となる。一方で、当該相関性が大きければ大きいほど、当該係数は小さい値となり、車両の位置の予測値の分散は小さい値となる。また、式（２２）は、車両の位置の推定値の分散が、車両の速度の誤差に対して指数関数的に減少することを示す。
【０１１３】
したがって、車両の位置の推定値の誤差分散は、車両の位置の誤差と、車両の速度の誤差との相関性に基づき決定されていることを示している。また、車両の位置の推定値の誤差は、車両の速度の誤差によって累積しないことも示す。つまり、車両の速度の誤差の累積に伴い、車両の位置の推定値の誤差分散であるｐ１１（ｋ｜ｋ）を累積させるためには、車両の位置の誤差と、車両の速度の誤差との相関性を小さくする必要がある。つまり、ｐ１１（ｋ｜ｋ）を累積させるためには、カルマンフィルタにおいて、車両の位置の誤差と、車両の速度の誤差との相関性を小さくする必要がある。具体的には、移動距離に累積する速度の誤差に関しては、カルマンフィルタに車両の速度の誤差を入力し、車両の位置の誤差を推定するのではなく、カルマンフィルタに車両の位置の誤差として車両の速度の誤差を入力し、車両の位置の誤差を推定する必要がある。
【０１１４】
したがって、式（１２）に示すように、移動距離に累積する速度の誤差に関しては、車両の速度の誤差を、車両の速度と微積分の関係にある車両の位置の誤差としてカルマンフィルタに入力する。これによって、移動距離に累積する速度の誤差の分だけ、車両の速度の誤差と車両の位置の誤差との相関性を小さくでき、累積する車両の速度の誤差を累積する車両の位置の誤差として、車両の位置の予測値の誤差共分散行列を算出できる。つまり、累積する車両の速度の誤差を考慮した、車両の位置の予測値の誤差共分散行列を算出できる。
【０１１５】
しかしながら、車両の位置の誤差として車両の速度の誤差をカルマンフィルタに入力する際、式（４）をそのまま用いることができない。
【０１１６】
式（４）において、右辺の第１項は、前回に推定した車両状態推定値の誤差共分散行列を、状態方程式から求められるヤコビ行列により今回の車両の状態量の誤差共分散行列に変換している。一方、右辺の第２項は、システム雑音の誤差共分散行列である。つまり、式（４）は、前回の車両状態推定値の誤差共分散行列を今回の車両状態推定値の誤差共分散行列に変換したものに、システム雑音の誤差共分散行列を加算することで、車両状態予測値の誤差共分散行列を算出する式である。ここで、式（４）は、式（２５）に示す命題が用いられている。
【０１１７】
【数２５】
0006677533.tif 000026
【０１１８】
式（２５）は、確率変数Ｘと確率変数Ｙとが互いに独立である場合に、分散Ｖａｒ（）について成立する式を示す。つまり、カルマンフィルタにおいて車両状態予測値の誤差共分散行列を算出する際に用いられる式（４）は、車両の状態量の誤差共分散行列と、システム雑音の誤差共分散行列とが互いに独立であることが前提となっている。しかしながら、システム雑音の誤差共分散行列を累積する誤差の誤差共分散行列とする場合、式（４）では、車両状態予測値の誤差共分散行列が正確に算出できない。これは、車両の状態量の誤差共分散行列と、車両の状態量として入力した累積する誤差の誤差共分散行列とが独立でないためである。例えば、車両の位置の誤差が時間の経過で累積する速度の誤差に決定されることから、車両の位置の誤差と、車両の速度の誤差とは、一方の変化に伴い一方も変化する変数、すなわち、互いに独立しない変数である。そのため、確率変数Ｘと確率変数Ｙとが独立でない場合に成立する命題を用いる必要がある。当該命題は、式（２６）により表される。
【０１１９】
【数２６】
0006677533.tif 000027
【０１２０】
式（２６）において、Ｃｏｖ（）は共分散を示す。したがって、式（２４）と比較すると、式（２６）は、確率変数Ｘと確率変数Ｙとの共分散が入力されている。
【０１２１】
式（２６）に従った式（９）を用いることにより、車両の状態量の誤差共分散行列と、累積する誤差の誤差共分散行列との独立性をなくして、正確に車両状態予測値の誤差共分散行列を算出できる。
【０１２２】
以上、車両の速度の誤差を車両の位置の誤差としてカルマンフィルタに入力することを説明した。角速度の誤差を方位の誤差としてカルマンフィルタに入力することについての説明も、同様に説明できる。
【０１２３】
以上の説明のように、観測量と微積分の関係にある状態量の誤差として、観測量の誤差が入力され、また、観測量の誤差と、車両の状態量の誤差との独立性のない式により、累積する誤差を考慮した、車両状態予測値の誤差共分散行列を算出できる。
【０１２４】
図３に示すフローチャートの説明に戻り、推定部２２は、予測処理を実行すると、推定処理を実行する。推定部２２は、式（５）〜式（８）に基づいて、車両状態推定値と、車両状態推定値の誤差共分散行列とを算出する。
【０１２５】
前述の通り、カルマンゲインＫｋは、観測雑音の誤差共分散行列と、車両予測値の誤差共分散行列とに応じて、車両状態推定値と車両状態推定値の誤差共分散行列とを適切に設定する係数である。予測処理において、推定部２２は、累積する誤差を考慮した車両状態予測値の誤差共分散行列を算出できるため、カルマンゲインＫｋを正確に算出できる。したがって、推定部２２は、車両状態推定値と車両状態推定値の誤差とを精度よく算出できる。
【０１２６】
また、式（８）に示すように、車両状態推定値の誤差共分散行列は、車両状態予測値の誤差共分散行列に（Ｉ−ＫｋＨ）を掛けて算出される。したがって、車両状態予測値の誤差共分散行列は、累積する誤差を考慮した誤差共分散行列であるため、算出される車両状態推定値の誤差共分散行列も累積する誤差を考慮した誤差共分散行列である。
【０１２７】
このように、推定部２２は、累積する誤差を考慮した車両状態予測値の誤差共分散行列を算出できるため、カルマンゲインＫｋが正確に算出でき、カルマンフィルタにより精度よく車両の状態を推定できる。
また、カルマンフィルタで車両の状態を推定しているため、誤差を推定の対象としたカルマンフィルタに基づき車両の状態を推定する構成と比較し、誤差とは別に車両の状態を推定する必要がなく、容易に累積する誤差を考慮した車両の状態を推定できる。
【０１２８】
なお、状態量として車両の速度を推定する場合、推定部２２は、車両の加速度の誤差を車両の速度の誤差として入力したカルマンフィルタに基づき車両の速度の予測値の誤差共分散行列を算出し、車両の速度を推定する構成としてもよい。これにより、累積する加速度の誤差を考慮した車両の速度の予測値の誤差共分散行列が算出されるため、カルマンゲインＫｋを正確に算出でき、カルマンフィルタにより車両の速度を精度よく推定できる。
【０１２９】
上述した通り、推定部２２は、ＧＰＳ受信部５が出力する車両の位置を示す情報と車両の方位を示す情報とに基づき、観測部２１が観測した車両の位置と車両の方位とを加味して、車両状態推定値と車両状態推定値の誤差共分散行列とを算出し、車両の状態を推定する構成である。しかしながら、観測部２１が観測する車両の位置と車両の方位とを加味して、推定部２２が車両の状態を推定する構成に限定されない。すなわち、推定部２２は、ＧＰＳ受信部５からの出力に基づくことなく、車両の状態を推定してもよい。
【０１３０】
例えば、車両がＧＰＳ電波をＧＰＳ受信部５により受信できない環境下にいる場合、推定部２２は、ＧＰＳ受信部５からの出力に基づく車両の位置及び車両の方位の成分のない観測方程式に基づき車両の状態量を推定する。
【０１３１】
式（１２）に示すように、推定部２２は、車両状態予測値の誤差共分散行列を算出する際、車速センサー６及びジャイロセンサー７から出力される信号の間隔に基づき、車両状態予測値の誤差共分散行列を算出する。また、推定部２２は、式（３）に示すように、車両状態予測値を、前回推定した車両状態推定値と当該間隔に基づき算出する。したがって、推定部２２は、車両状態予測値、及び、車両状態予測値の誤差共分散行列の算出に際し、ＧＰＳ受信部５の受信間隔に依存しない。したがって、推定部２２は、式（５）において、ＧＰＳ受信部５からの出力に基づく観測対象の成分のない式に基づくことにより、車両状態推定値と車両状態推定値の誤差共分散行列とを算出できる。
【０１３２】
つまり、推定部２２は、ＧＰＳ受信部５のＧＰＳ電波の受信の間隔に依存せず、車両の状態を推定できる。これは、推定部２２が車速センサー６、ジャイロセンサー７、及び、加速度センサー８の検出間隔に基づき、車両の状態を推定することを示す。したがって、推定部２２が車速センサー６、ジャイロセンサー７、及び、加速度センサー８の検出間隔が、ＧＰＳ受信部５の受信間隔より短い場合において以下の効果を奏する。すなわち、当該受信間隔に基づいてカルマンフィルタにより車両の状態を推定する構成より、車両の状態を推定について高い頻度で実行でき、また、ＧＰＳ受信部５の受信環境に依存することなく、車両の状態を推定できる。
【０１３３】
前述した通り、ナビゲーション装置１の制御部２は、推定した車両の状態に基づき、マップマッチングの対象となるリンクの評価を実行する。
【０１３４】
図４は、推定した車両の状態に基づくマップマッチングの処理を説明するため、地図上の道路Ｒ１と、道路Ｒ２とを示す図である。図４において、道路Ｒ１は、方向Ｙ１に向かって延在する道路である。また、道路Ｒ２は、方向Ｙ１と並行でない方向Ｙ２に向かって延在する道路である。
【０１３５】
図４において、リンクＬ１は、道路Ｒ１に対応するリンクである。また、リンクＬ２は、道路Ｒ２に対応するリンクである。
【０１３６】
また、図４において、マークαは、推定部２２により推定された車両の位置を示すマークである。以下、図４を用いた説明では、推定部２２により、車両が、位置Ｍ１に位置し、進行方向Ｘ１に向かって走行しているものと推定されたとする。図４に示すように、位置Ｍ１は、道路Ｒ１上の位置であって、道路Ｒ１の幅方向の中心から、方向Ｙ１に向かって右側に離間した位置である。
【０１３７】
道路Ｒ１と、道路Ｒ２と、推定された車両の位置との関係が、図４に示す関係である場合、マップマッチングに際し、まず、制御部２は、地図データ３ａを参照し、マップマッチングの候補となるリンクであるマップマッチング候補リンクを取得する。制御部２は、マップマッチング候補リンクとして、推定された車両の位置から予め設定された所定の範囲内に位置し、かつ、車両とリンクとの方位の誤差が所定範囲内である１又は複数のリンクを取得する。図４の例では、リンクＬ１、及び、リンクＬ２は、推定された車両の位置である位置Ｍ１から所定範囲に位置し、また、車両とリンクとの方位の誤差が所定範囲内のリンクである。このため、制御部２は、リンクＬ１、及び、リンクＬ２をマップマッチング候補リンクとして取得する。
【０１３８】
次いで、制御部２は、取得したマップマッチング候補リンクのそれぞれについて、リンクを評価する評価量を算出する。
【０１３９】
評価量は、推定された車両の位置とマップマッチング候補リンクの位置との誤差、及び、推定された車両の方位とマップマッチング候補リンクの方位との誤差に基づいて、下記に示す式（２７）により算出される値である。
【０１４０】
τ＝δｘ２／Δｘ２＋δｙ２／Δｙ２＋δθ２／Δθ２・・・（２７）
【０１４１】
ここで、推定された車両の位置とマップマッチング候補リンクの位置との誤差とは、推定された車両の位置からマップマッチング候補リンクに対して垂線を延ばした場合における垂線及びマップマッチング候補リンクとの交点と、推定された車両の位置との、ｘ軸方向、及び、ｙ軸方向における差である。式（２７）において、δｘは、当該交点のｘ座標と、推定された車両の位置のｘ座標との差を示す。また、δｙは、当該交点のｙ座標と、推定された車両の位置のｙ座標との差を示す。
【０１４２】
図４の例において、推定された車両の位置とリンクＬ１との位置の誤差は、位置Ｍ１からリンクＬ１に対して延ばした垂線Ｓ１及びリンクＬ１の交点ＭＭ１と、位置Ｍ１とのｘ軸方向、及び、ｙ軸方向における差である。また、図４において、推定された車両の位置とリンクＬ２との位置の誤差は、位置Ｍ１からリンクＬ２に対して延ばした垂線Ｓ２及びリンクＬ２の交点ＭＭ２と、位置Ｍ１とのｘ軸方向、及び、ｙ軸方向における差である。
【０１４３】
また、推定された車両の方位とマップマッチング候補リンクの方位との誤差とは、推定された車両の方位に対応する角度と、マップマッチング候補リンクの方位に対応する角度との差である。
【０１４４】
前述した通り、車両の方位とは、車両の進行方向の方角を意味する。図４の例では、推定された車両の方位は、進行方向Ｘ１の方角である。また、車両の方位に対応する角度とは、東に向かう方向を基準として、東に向かう方向と、車両の方位との反時計回りの離間角度を意味する。図４の例では、推定された車両の方位に対応する角度は、仮想直線Ｋ１が東西に延びる仮想直線であるとした場合に、角度θ１である。
【０１４５】
また、リンクの方位とは、リンクが延在する方向の方角を意味する。リンクが延在する方向とは、リンクに沿った２つの方向のうち、リンクに対応する道路において車両が走行可能な方向に対応する方向を意味する。また、リンクの方位に対応する角度とは、東に向かう方向を基準として、東に向かう方向と、リンクの方位との反時計回りの離間角度を意味する。
【０１４６】
図４の例では、リンクＬ１の方位は、方向Ｚ１の方角である。また、リンクＬ１の方位に対応する角度は、仮想直線Ｋ２が東西に延びる仮想直線であるとした場合に、角度θ２である。また、図４の例では、リンクＬ２の方位は、方向Ｚ２の方角である。また、リンクＬ２の方位に対応する角度は、仮想直線Ｋ３が東西に延びる仮想直線であるとした場合に、角度θ３である。
【０１４７】
図４の例において、車両とリンクＬ１との方位誤差は、角度θ１と角度θ２との差である。また、車両とリンクＬ２との方位誤差は、角度θ１と角度θ３との差である。つまり、図４の例において、式（２７）では、δθは、角度θ１と角度θ２との差、または、角度θ１と角度θ３との差を示す。
【０１４８】
式（２７）に示すように、評価量τは、δｘを二乗した値にΔｘを二乗した値で割ったものと、δｙを二乗した値にΔｙを二乗した値で割ったものと、δθを二乗した値にΔθを二乗した値で割ったものと、の和により算出される。Δｘは、推定部２２が推定した車両のｘ軸方向の位置の誤差を示す。つまり、Δｘ２は、推定部２２が推定した車両のｘ軸方向の位置の分散を示す。また、Δｙは、推定部２２が推定した車両のｙ軸方向の位置の誤差を示す。つまり、Δｙ２は、推定部２２が推定した車両のｙ軸方向の位置の分散を示す。また、Δθは、推定部２２が推定した車両の方位の誤差を示す。つまり、Δθ２は、推定部２２が推定した車両の方位の分散を示す。Δｘ２、Δｙ２、及び、Δθ２は、推定部２２が推定した車両状態推定値の誤差共分散行列が含む分散である。このように、評価量は、車両の位置とリンクの位置との誤差を推定した車両の位置の誤差により無次元化した値と、車両の方位とリンクの方位との誤差を、推定した車両の位置の誤差により無次元化した値との和により算出される値である。
【０１４９】
図４の例では、制御部２は、マップマッチング候補リンクとして取得したリンクＬ１の評価量とリンクＬ２の評価量とを算出する。
【０１５０】
リンクＬ１の評価量をτ１とし、位置Ｍ１と交点ＭＭ１とのｘ座標の差をδｘ１とし、位置Ｍ１と交点ＭＭ１とのｙ座標の差をδｙ１とし、角度θ１と角度θ２との差をδθ１とした場合、評価量τ１は、下記に示す式（２８）で表される。
【０１５１】
τ１＝δｘ１２／Δｘ２＋δｙ１２／Δｙ２＋δθ１２／Δθ２・・・（２８）
【０１５２】
一方、リンクＬ２の評価量をτ２とし、位置Ｍ１と交点ＭＭ２とのｘ座標の差をδｘ２とし、位置Ｍ１と交点ＭＭ２とのｙ座標の差をδｙ２とし、角度θ１と角度θ３との差をδθ２とした場合、評価量τ２は、下記に示す式（２９）で表される。
【０１５３】
τ２＝δｘ２２／Δｘ２＋δｙ２２／Δｙ２＋δθ２２／Δθ２・・・（２９）
【０１５４】
図４の例では、距離ｌ２より距離ｌ１が小さい。すなわち、位置Ｍ１と交点ＭＭ２とにおけるｘ軸座標及びｙ座標の差より、位置Ｍ１と交点ＭＭ２とにおけるｘ座標及びｙ座標の差が小さい。また、図４の例では、角度θ１と角度θ３との差より角度θ１と角度θ２との差が小さい。このため、リンクＬ１の評価量であるτ１の値の方が、リンクＬ２の評価量であるτ２の値よりも小さくなる。したがって、制御部２は、評価量が小さいリンクＬ１を、車両の現在位置を対応付けるリンクとして決定する。
【０１５５】
このように、制御部２は、推定した車両の状態に基づいて、マップマッチングの対象となるリンクの評価を実行する。上述したように、制御部２は、リンクを評価する際、推定部２２が推定した車両状態推定値の誤差共分散行列が含む分散に基づき、リンクを評価する。当該誤差共分散行列は、累積する誤差を考慮して算出された車両状態予測値の誤差共分散行列に基づき算出された行列であり、精度よく算出された誤差を含む行列である。したがって、制御部２は、当該誤差共分散行列が含む分散をリンクの評価に用いることにより、正確にリンクを評価できる。
【０１５６】
以上、説明したように、ナビゲーション装置１（車載装置）は、車両の変動に関する観測量を観測する観測部２１と、カルマンフィルタにより、車両の状態を示す状態量を推定する推定部２２と、を備える。推定部２２は、車両の状態量の予測値を算出し、観測量と微積分の関係にある状態量の誤差として観測量の誤差が入力されたカルマンフィルタにより、予測値の誤差共分散行列（誤差）を算出し、算出した予測値と予測値の誤差共分散行列とに基づいて、カルマンフィルタにより車両の状態量の推定値と推定値の誤差共分散行列とを算出する。
【０１５７】
これにより、観測量の誤差を微積分の関係にある状態量の誤差として入力することにより、観測量の誤差と状態量の誤差との相関性を小さくし、累積する観測量の誤差を考慮して車両状態予測値の誤差共分散行列を算出できるため、カルマンゲインＫｋを正確に算出でき、カルマンフィルタにより車両の状態を精度よく推定できる。
【０１５８】
また、推定部２２は、前回算出した車両状態推定値の誤差と、車両の状態量の誤差として入力された観測量の誤差との共分散が入力されたカルマンフィルタにより、予測値の誤差共分散行列を算出する。
【０１５９】
前述した通り、式（４）は、状態量の誤差共分散行列と、システム雑音の誤差共分散行列とが互いに独立であることが前提となる命題が用いられている。しかしながら、システム雑音の誤差共分散行列を累積する誤差の誤差共分散行列とする場合、車両の状態量の誤差共分散行列と、累積する誤差の誤差共分散行列とが独立でない。そこで、式（２６）に従った式（９）、すなわち、前回推定した車両状態推定値の誤差と、状態量の誤差として入力された観測量の誤差との共分散が入力された式を用いることにより、車両の状態量の誤差共分散行列と、累積する誤差の誤差共分散行列との独立性をなくして、正確に車両状態予測値の誤差共分散行列を算出できる。
【０１６０】
また、観測部２１は、車両の速度を検出する車速センサー６からの出力に基づき、車両の速度を観測する。推定部２２は、車両の位置の誤差として車両の速度の誤差が入力されたカルマンフィルタにより、予測値の誤差共分散行列を算出する。
【０１６１】
これにより、車両の速度の誤差と車両の位置の誤差との相関性を小さくすることで、累積する車両の速度の誤差を考慮でき、カルマンフィルタにより車両の位置を精度よく推定できる。
【０１６２】
また、観測部２１は、車両の角速度を検出するジャイロセンサー７（角速度センサー）からの出力に基づき、車両の角速度を観測する。推定部２２は、車両の方位の誤差として車両の角速度の誤差が入力されたカルマンフィルタにより、予測値の誤差共分散行列を算出する。
【０１６３】
これにより、車両の角速度の誤差と車両の方位の誤差との相関性を小さくすることで、累積する車両の角速度の誤差を考慮でき、カルマンフィルタにより車両の方位を精度よく推定できる。
【０１６４】
観測部２１は、車両の加速度を検出する加速度センサー８からの出力に基づき、車両の加速度を観測する。推定部２２は、車両の速度の誤差として車両の加速度の誤差が入力されたカルマンフィルタにより、予測値の誤差共分散行列を算出する。
【０１６５】
これにより、車両の加速度の誤差と車両の速度の誤差との相関性を小さくすることで、累積する車両の加速度の誤差を考慮でき、カルマンフィルタにより車両の速度を精度よく推定できる。
【０１６６】
上述した実施形態は、あくまでも本発明の一態様を例示するものであって、本発明の趣旨を逸脱しない範囲で任意に変形、及び応用が可能である。
【０１６７】
例えば、図１は、本願発明を理解容易にするために、ナビゲーション装置１の機能構成を主な処理内容に応じて分類して示した概略図であり、ナビゲーション装置１の構成は、処理内容に応じて、さらに多くの構成要素に分類することもできる。また、１つの構成要素がさらに多くの処理を実行するように分類することもできる。
【０１６８】
また、例えば、図２、及び、図３のフローチャートの処理単位は、制御部２の処理を理解容易にするために、主な処理内容に応じて分割したものであり、処理単位の分割の仕方や名称によって、本発明が限定されることはない。制御部２の処理は、処理内容に応じて、さらに多くの処理単位に分割してもよい。また、１つの処理単位がさらに多くの処理を含むように分割してもよい。
【０１６９】
また、例えば、上述した実施形態では、ナビゲーション装置１として車載型を例示したが、ナビゲーション装置１の形態は任意であり、例えば歩行者が携帯するポータブル型の装置でも良い。
【符号の説明】
【０１７０】
１  ナビゲーション装置（車載装置）
６  車速センサー（センサー）
７  ジャイロセンサー（センサー）
８  加速度センサー（センサー）
２１  観測部
２２  推定部
【図１】
0006677533.tif 000028
【図２】
0006677533.tif 000029
【図３】
0006677533.tif 000030
【図４】
0006677533.tif 000031