(57)【要約】【課題】車重の推定に利用可能であり、かつ推定の精度および速度を高めることが可能な、車両の状態量を推定する技術を実現する。【解決手段】状態量推定装置は、データ格納部（１０１）、予測量演算部（１０２）、取得部（１０７）、カルマンゲイン演算部（１０３）、推定状態量および推定共分散を算出する推定量演算部（１０４）およびプロセスノイズ共分散を修正するプロセスノイズ共分散修正部（１０６）を備える。算出または修正された推定状態量、推定共分散およびプロセスノイズ共分散は、それぞれ、状態量、状態共分散およびプロセスノイズ共分散としてデータ格納部（１０１）に書き込まれ、次回の状態量推定のための演算に用いられる。【選択図】図１
2020064015.tif 000002
【特許請求の範囲】
【請求項１】
車両に作用する状態量、前記状態量の共分散である状態共分散、および、演算ノイズの共分散であるプロセスノイズ共分散、を格納するデータ格納部と、
前記状態量から予測状態量を演算し、前記状態共分散および前記プロセスノイズ共分散から前記予測状態量の共分散である予測共分散を演算する予測量演算部と、
車両のセンサ値を取得する取得部と、
前記センサ値および前記予測共分散を用いてカルマンゲインを演算するカルマンゲイン演算部と、
前記カルマンゲイン、前記予測状態量および前記センサ値、からカルマンフィルタを用いて、推定状態量を演算し、前記カルマンゲイン、前記予測共分散および前記センサ値、からカルマンフィルタを用いて、推定共分散を演算する推定量演算部と、
前記推定状態量を用いて前記プロセスノイズ共分散を修正するプロセスノイズ共分散修正部と、を備え、
前記推定状態量が前記状態量として、前記推定共分散が前記状態共分散として、修正された前記プロセスノイズ共分散が前記プロセスノイズ共分散として、前記データ格納部に書き込まれる、状態量推定装置。
【請求項２】
前記推定状態量から前記車両の重量成分を抽出する、重量成分抽出部をさらに備える、請求項１に記載の状態量推定装置。
【請求項３】
前記重量成分抽出部は、前記重量成分の平均と、前記重量成分の標準偏差と、を演算する、請求項２に記載の状態量推定装置。
【請求項４】
前記重量成分抽出部は、
初期状態をオフとする安定フラグと、
前記推定状態量の演算回数、前記標準偏差の値、および、前記安定フラグのオンオフ、からなる群から選ばれる一以上を表す情報に応じて前記重量成分を決定する重量決定部と、
を含み、
前記重量決定部は、
前記推定状態量の演算回数が所定回数以下の場合には、最新の抽出した前記重量成分を抽出した前記重量成分と決定し、
前記推定状態量の演算回数が前記所定回数を超え、かつ、前記標準偏差が所定の閾値未満の場合には、前記平均を抽出した前記重量成分と決定して前記安定フラグをオンにし、
前記安定フラグがオンであるが、前記推定状態量の演算回数が前記所定回数を超え、かつ前記標準偏差が所定の閾値以上である場合には、前回以前の前記安定フラグをオンにしたときに抽出した前記重量成分を抽出した前記重量成分と決定し、
前記推定状態量の演算回数が前記所定回数を超え、前記標準偏差が所定の閾値以上であり、かつ、前記安定フラグがオフの場合には、最新の抽出した前記重量成分を抽出した前記重量成分と決定する、
請求項３に記載の状態量推定装置。
【請求項５】
前記プロセスノイズ共分散修正部は、前記重量成分抽出部が抽出した前記重量成分を用いて前記プロセスノイズ共分散を修正する、請求項２に記載の状態量推定装置。
【請求項６】
前記プロセスノイズ共分散修正部は、前記重量成分抽出部が演算した前記標準偏差を用いて前記プロセスノイズ共分散を修正する、請求項３または４に記載の状態量推定装置。
【請求項７】
前記重量成分抽出部は、抽出した前記重量成分を前記車両の現在の重量として出力する、請求項２〜４のいずれか一項に記載の状態量推定装置。
【請求項８】
懸架装置を有する車両に作用する状態量を推定して、前記懸架装置の減衰力を前記状態量に応じて制御する制御装置であって、
車両に作用する状態量、前記状態量の共分散である状態共分散、および、演算ノイズの共分散であるプロセスノイズ共分散、を格納するデータ格納部と、
前記状態量から予測状態量を演算し、前記状態共分散および前記プロセスノイズ共分散から前記予測状態量の共分散である予測共分散を演算する予測量演算部と、
車両のセンサ値を取得する取得部と、
前記センサ値および前記予測共分散を用いてカルマンゲインを演算するカルマンゲイン演算部と、
前記カルマンゲイン、前記予測状態量および前記センサ値、からカルマンフィルタを用いて、推定状態量を演算し、前記カルマンゲイン、前記予測共分散および前記センサ値、からカルマンフィルタを用いて、推定共分散を演算する推定量演算部と、
前記推定状態量を用いて前記プロセスノイズ共分散を修正し、するプロセスノイズ共分散修正部と、を備え、
前記推定状態量が前記状態量として、前記推定共分散が前記状態共分散として、修正された前記プロセスノイズ共分散が前記プロセスノイズ共分散として、前記データ格納部に書き込まれる、制御装置。
【請求項９】
車両に作用する状態量、前記状態量の共分散である状態共分散、および、演算ノイズの共分散であるプロセスノイズ共分散、をデータ格納部に格納するステップと、
前記状態量から予測状態量を演算し、前記状態共分散および前記プロセスノイズ共分散から前記予測状態量の共分散である予測共分散を演算するステップと、
車両のセンサ値を取得するステップと、
前記センサ値および前記予測共分散を用いてカルマンゲインを演算するステップと、
前記カルマンゲイン、前記予測状態量および前記センサ値、からカルマンフィルタを用いて、推定状態量を演算するステップと、
前記カルマンゲイン、前記予測共分散および前記センサ値、からカルマンフィルタを用いて、推定共分散を演算するステップと、
前記推定状態量を用いて前記プロセスノイズ共分散を修正するステップと、
前記推定状態量を前記状態量とし、前記推定共分散を前記状態共分散として、修正された前記プロセスノイズ共分散を前記プロセスノイズ共分散として、前記データ格納部に書き込むステップと、
を含む、状態量推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、状態量推定装置、制御装置、および状態量推定方法に関する。
【背景技術】
【０００２】
近年、自動車などの車両では、例えばその運転の安全性および快適性を高めるための種々の制御が行われている。このような車両の運転に係る制御を適切に実施するために、走行中の車両の質量（車重）を精度よく把握することが求められている。車重を推定する技術には、例えば、各種センサによりその車両に依存する状態量を取得し、運動方程式に当て嵌め、さらに２つの運動方程式を減算することによりノイズを相殺して車重を推定する技術が知られている。各種センサで取得した状態量のうち、状態量の推定に適さない数値は、通常、閾値を用いて破棄される（例えば、特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１３−１５２１７０号公報（２０１３年８月８日公開）
【発明の概要】
【発明が解決しようとする課題】
【０００４】
上記の技術では、センサで取得したセンサ値のうち、状態量の推定に適さない一つでも数値が含まれると、その他の適する数値も破棄されることがある。そのため、車両の状態量を推定する機会が少なくなり、少ない機会で状態量を推定しなければならない場合がある。よって、状態量の推定の精度が不十分になることがあり、また状態量の推定の安定化に時間がかかることがある。このように、上記の技術では、車両の状態量の推定における精度および速度を高める観点から検討の余地が残されている。
【０００５】
本発明の一態様は、より高い精度で迅速に車両の状態量を推定する技術を実現することを目的とする。
【課題を解決するための手段】
【０００６】
上記の課題を解決するために、本発明の一態様に係る状態量推定装置は、車両に作用する状態量、前記状態量の共分散である状態共分散、および、演算ノイズの共分散であるプロセスノイズ共分散、を格納するデータ格納部と、前記状態量から予測状態量を演算し、前記状態共分散および前記プロセスノイズ共分散から前記予測状態量の共分散である予測共分散を演算する予測量演算部と、車両のセンサ値を取得する取得部と、前記センサ値および前記予測共分散を用いてカルマンゲインを演算するカルマンゲイン演算部と、前記カルマンゲイン、前記予測状態量および前記センサ値、からカルマンフィルタを用いて、推定状態量を演算し、前記カルマンゲイン、前記予測共分散および前記センサ値、からカルマンフィルタを用いて、推定共分散を演算する推定量演算部と、前記推定状態量を用いて前記プロセスノイズ共分散を修正するプロセスノイズ共分散修正部と、を備え、前記推定状態量が前記状態量として、前記推定共分散が前記状態共分散として、修正された前記プロセスノイズ共分散が前記プロセスノイズ共分散として、前記データ格納部に書き込まれる。
【発明の効果】
【０００７】
本発明の一態様によれば、より高い精度で迅速に車両の状態量を推定する技術を実現することができる。
【図面の簡単な説明】
【０００８】
【図１】本発明の実施形態１に係る状態量推定装置の機能的構成の一例を示すブロック図である。
【図２】実施形態１に係る状態量を推定する処理の流れの一例を示すフローチャートである。
【図３】実施形態１においてプロセスノイズ共分散を修正する処理における第一の態様の流れを示すフローチャートである。
【図４】実施形態１においてプロセスノイズ共分散を修正する処理における第二の態様の流れを示すフローチャートである。
【図５】実施形態１においてプロセスノイズ共分散を修正する処理における第三の態様の流れを示すフローチャートである。
【図６】実施形態１においてプロセスノイズ共分散を修正する処理における第四の態様の流れを示すフローチャートである。
【図７】実施形態１に係る状態量推定装置が適用される車両の構成の一例を模式的に示す図である。
【発明を実施するための形態】
【０００９】
〔実施形態１〕
以下、本発明の一実施形態について、詳細に説明する。
【００１０】
［状態量推定のロジックの説明］
車両に発生する運動量は、運動方程式Ｆ＝ｍａで表すことができる。ここで、Ｆはエンジン出力から各種走行抵抗値を減算したものであり、ｍは車両の重量（車重）であり、ａは車両の加速度である。ａは、車速Ｖの時間微分値で求めることができるため、下記式１で表すことができる。式１中、Ｆｘはエンジン由来出力を表し、Ｆａｉｒは空気抵抗成分を表し、Ｆｓｕｒｆａｃｅは路面抵抗成分を表す。
【００１１】
【数１】
2020064015.tif 000003
【００１２】
ここで、エンジン由来出力であるＦｘは、車輪半径で車輪トルクを除することで求められ、路面の傾き分の重力加速度を受けることから、下記式２で表すことができる。下記式２中、τは車輪トルクを表し、Ｒは車輪半径を表し、θｒは路面の傾斜角度を表す。
【００１３】
【数２】
2020064015.tif 000004
【００１４】
本明細書でいうところの「車輪トルク」とは、走行中の車両の駆動源が発生させる、車輪に働く、車両を加減速させる方向に働くトルクである。たとえば、内燃機関を駆動源とする車両では、車輪トルクとは、該車両の車輪に掛けられているトルクである。内燃機関の車輪トルクは、内燃機関が発生させたトルクを、空燃比、外気温およびスロットルバルブ解放量などから推定し、得られた推定値に、車両毎に設定されている伝達損失係数と、各減速機構による所定の減速比とを乗ずることにより求められる。また、車輪内に駆動源として電動モータを独立して持つ車両では、車輪トルクは、該車両の各車輪に掛けられているトルクの合計である。当該車両の車輪トルクは、各モータの仕事効率および当該モータに印加される電圧から推定されるトルクに、車両毎に設定されている伝達損失係数と、各減速機構による所定の減速比とを乗ずることにより求められる。また、車両が電気式デフなどの差動制御装置（ＬＳＤ）をさらに有する場合では、これらの装置の作動状況をさらに参照して上記車輪トルクを算出してもよい。
【００１５】
空気抵抗成分であるＦａｉｒは、その車両に特有の空気抵抗係数に車速の二乗を乗じることで求められるから、下記式３で表すことができる。下記式３中、Ｃｄはその車両の空気抵抗係数を表し、Ｖは車速を表す。
【００１６】
【数３】
2020064015.tif 000005
【００１７】
路面抵抗成分であるＦｓｕｒｆａｃｅは、路面の摩擦係数に車重を乗じたものの路面角度の重力加速度分であるから、下記式４で表すことができる。下記式４中、μはその車両の路面抵抗係数を表し、ｍは車重を表し、θｒは路面の傾斜角度を表す。
【００１８】
【数４】
2020064015.tif 000006
【００１９】
したがって、車両に発生する運動量の運動方程式Ｆ＝ｍａは、式１〜式４により、下記式５、式６で表される。
【００２０】
【数５】
2020064015.tif 000007
【００２１】
車両の走行状態を制御する場合では、車両の前後Ｇを検知する前後Ｇセンサが一般に用いられる。この前後Ｇセンサの値Ｇｓｅｎｓは、車速を時間微分した値の車体ピッチ角分と、路面角度と車体ピッチ角の合計の重力加速度成分であるから、下記式７、式８で表すことができる。下記式７、式８中、θｐは車両のピッチ角度を表し、θｒは路面の傾斜角度を表す。
【００２２】
【数６】
2020064015.tif 000008
【００２３】
式６および式８より、式９が導き出され、式９から式１０が導き出される。
【００２４】
【数７】
2020064015.tif 000009
【００２５】
ここで、θｐは十分小さいとして、θｐ≒０として近似を行うと、式１０から下記式１１が導き出される。
【００２６】
【数８】
2020064015.tif 000010
【００２７】
式１１において、左辺をＹに、ＧｓｅｎｓをＧに、ｍｇｃｏｓθｒ（μ−ｓｉｎθｐ）をＣにそれぞれ置き換えることにより、下記式１２が導き出される。
【００２８】
【数９】
2020064015.tif 000011
【００２９】
式１２においてＹおよびＧは、前述の通りセンサの値を参照可能である。また、ｍは車重であり、車載物の重量が実質的に変わらない適当な期間であれば実質的に一定であるから、定数とする。また、式１２のＣについて、式１１より、ｍは前述の通り定数であり、ｇは重力加速度であるから、これも定数である。車体ピッチ角ｓｉｎθｐは、前述の通りθｐ≒０であるから、定数（０）とする。よって、本実施形態では、路面角度ｃｏｓθｒおよび路面抵抗μは、変化分をすべて後述するカルマンフィルタのノイズとして扱うことができることから、ここでは定数として扱う。したがって、式１２においてＣは定数と考える。
【００３０】
一方で、下記式１３は、カルマンフィルタの予測関数である。式１３において、ｔは演算回数を表し、ｘａt−１は前回の演算により推定した、車両に作用する状態量を表し、ｘｂtは今回の演算による、車両に作用する状態の予測量（予測状態量）を表す。車両に作用する状態量は、車両に係る物理量であり、その例には、車両の重量（車重）および道路の勾配成分が含まれる。
【００３１】
【数１０】
2020064015.tif 000012
【００３２】
式１４は、式１３からｘｂを行列で表示したものである。前述のようにｍおよびＣは定数であるから、式１３は下記式１４に基づき変化しないと考えられる。したがって、前回推定の状態量ｘａt−１と今回の予測状態量ｘｂtとは同一であり、下記式１５が成り立つ。
【００３３】
【数１１】
2020064015.tif 000013
【００３４】
また、予測状態量ｘｂtの共分散（予測共分散）Ｐ−tは、下記式１６で表される。式１６において、Ｐt−１は前回の演算により推定された共分散を（状態共分散）表し、Ｑt−１は前回の演算により推定されたプロセスノイズ共分散を表す。状態共分散は、例えば、車両に作用する状態量とプロセスノイズ共分散とを二変数とする共分散である。プロセスノイズ共分散は、例えば、車重ｍの推定に係るノイズと路面傾斜等Ｃの推定に係るノイズとを二変数とする共分散である。
【００３５】
【数１２】
2020064015.tif 000014
【００３６】
ここで、カルマンフィルタの観測方程式として、式１２を下記式１７に対応させる。式１７中、ｈtはＹに対応し、ｍｂはｍに対応し、ＣｂはＣに対応する。
【００３７】
【数１３】
2020064015.tif 000015
【００３８】
式１１の左辺は、前述したようにセンサによって観測可能であり、また式１２におけるＹである。これをセンサ値（「観測量」とも言う）ｚとすると、下記式１８が導き出される。
【００３９】
【数１４】
2020064015.tif 000016
【００４０】
カルマンゲインは、下記式１９、２０から求められる。当該カルマンゲインは、カルマンフィルタのゲインであり、観測誤差が修正された（「最適化された」とも言える）ものとなっている。
【００４１】
【数１５】
2020064015.tif 000017
【００４２】
下記式２１によってカルマンフィルタ演算を実施し、推定状態量ｘａｔを得る。また、下記式２２によってカルマンフィルタ演算を実施し、推定共分散Ｐｔを得る。
【００４３】
【数１６】
2020064015.tif 000018
【００４４】
［機能的構成例］
図１は、本実施形態に係る状態量推定装置の機能的構成の一例を示すブロック図である。図１に示されるように、状態量推定装置１００は、データ格納部１０１、予測量演算部１０２、カルマンゲイン演算部１０３、推定量演算部１０４、重量成分抽出部１０５、プロセスノイズ共分散修正部１０６、取得部１０７およびセンサ類１０８を含む。
【００４５】
データ格納部１０１は、状態量（ｘａ）、状態共分散（Ｐ）およびプロセスノイズ共分散（Ｑ）の各データを格納している。
【００４６】
予測量演算部１０２は、状態量ｘａから予測状態量（ｘｂ）を演算する。また、予測量演算部１０２は、状態共分散Ｐおよびプロセスノイズ共分散Ｑから予測共分散（Ｐ−）を演算する。
【００４７】
取得部１０７は、車両のセンサ値（ｚ）を取得する。取得部１０７は、車両に配置されている種々のセンサ類１０８からセンサ値ｚを取得する。センサ類１０８は、例えば、車両の状態を表す物理量を検出するための種々のセンサ、および、センサが検出した値から当該物理量を算出する装置を含む。
【００４８】
カルマンゲイン演算部１０３は、センサ値ｚおよび予測共分散Ｐ−を用いてカルマンゲイン（Ｋ）を演算する。
【００４９】
推定量演算部１０４は、カルマンゲインＫ、予測状態量Ｐ−およびセンサ値ｚからカルマンフィルタを用いて、推定状態量ｘａを演算する。推定状態量ｘａは、予測状態量ｘｂから推定される状態量の推定値である。また、推定量演算部１０４は、カルマンゲインＫ、予測共分散Ｐ−およびセンサ値ｚからカルマンフィルタを用いて、推定共分散Ｐを演算する。
【００５０】
重量成分抽出部１０５は、推定状態量ｘａから車両の重量成分（ｍ）を抽出する。
【００５１】
プロセスノイズ共分散修正部１０６は、重量成分抽出部１０５が抽出した重量成分を用いてプロセスノイズ共分散Ｑを修正する。
【００５２】
また、推定量演算部１０４は、ある状態量ｘａに基づく推定状態量ｘａを次の状態量ｘａとしてデータ格納部１０１に格納する。また、推定量演算部１０４は、ある状態量ｘａに基づく推定共分散Ｐを次の状態共分散Ｐとしてデータ格納部１０１に格納する。さらに、プロセスノイズ共分散修正部１０６は、ある推定状態量から抽出された重量成分ｍを用いてあるプロセスノイズ共分散Ｑを修正し、次のプロセスノイズ共分散Ｑとしてデータ格納部１０１に格納する。このように、状態量推定装置１００は、今回演算した推定状態量ｘａが次回の演算における状態量ｘａとしてデータ格納部１０１に書き込まれ、今回演算した推定共分散Ｐが次回の演算における状態共分散Ｐとして、データ格納部１０１に書き込まれる。
【００５３】
なお、状態量推定装置１００は、重量成分抽出部１０５を含まず、プロセスノイズ共分散修正部１０６は、重量成分ｍではなく、その抽出前のデータである推定状態量ｘａを推定量演算部１０４から取得し、それを用いてプロセスノイズ共分散Ｑを修正してもよい。本発明においては、プロセスノイズ共分散を修正する態様には様々な態様が含まれる。以下、状態量の推定に係る処理について説明するが、プロセスノイズ共分散を修正する種々の態様については、特に後述する。
【００５４】
［処理例］
図２は、本実施形態に係る状態量を推定する処理の流れの一例を示すフローチャートである。
【００５５】
データ格納部１０１は、状態量ｘａｔ、状態共分散Ｐｔおよびプロセスノイズ共分散Ｑｔを格納している（ステップＳ２０１）。ｔは演算回数であり、初期値であれば０である。
【００５６】
ステップＳ２０２において、状態量推定装置１００は、演算回数ｔを一つ増やすように更新する。たとえば、今回の演算回数がｔであるとすると、前回の演算回数はｔ−１で表され、前々回の演算回数はｔ−２で表され、次回の演算回数はｔ＋１で表される。
【００５７】
次いで、予測量演算部１０２は、前回演算時または初期値の状態量から予測状態量を演算し、前回演算時または初期値の状態共分散およびプロセスノイズ共分散から予測状態量の共分散である予測共分散を演算する。たとえば、予測量演算部１０２は、予測状態量と予測共分散とを所定の間隔で連続して演算する。
【００５８】
すなわち、ステップＳ２０３において、予測量演算部１０２は、予測状態量ｘｂを前述の式１３に基づき演算し、前回の演算回数における状態量ｘａｔ−１を求めるべき予測状態量ｘｂｔとする。
【００５９】
また、ステップＳ２０４において、予測量演算部１０２は、予測共分散Ｐ−を前述の式１６に基づき演算し、前回の演算回数における状態共分散Ｐｔ−１およびプロセスノイズ共分散Ｑｔ−１との和により予測共分散Ｐ−ｔを求める。
【００６０】
一方、取得部１０７は、センサ類１０８から車両の状態を検出する各種センサ値（「観測量ｚ」とも言う）を取得する。たとえば、取得部１０７は、センサ値を所定の間隔で連続して取得する。観測量ｚは、個々のセンサのセンサ値そのものであってもよく、当該センサ値から算出された、走行状態にある車両の、その車両に依存する物理量の値であってもよい。たとえば、ステップＳ２０５において、取得部１０７は、前後Ｇセンサのセンサ値であるＧｓｅｎｓｔ、車速Ｖｔおよび車両におけるエンジン由来の出力値Ｆｘｔを取得する。ここでＧｓｅｎｓｔは、式２０におけるＨｔである。
【００６１】
ステップＳ２０６において、カルマンゲイン演算部１０３は、観測量ｚｔおよび予測共分散Ｐ−ｔを用いて、前述の式１９に基づいてカルマンゲインＫｔを演算する。たとえば、カルマンゲイン演算部１０３は、カルマンゲインＫを所定の間隔で連続して演算する。
【００６２】
ステップＳ２０７において、推定量演算部１０４は、カルマンゲインＫｔ、予測状態量ｘｂｔおよび観測量ｚｔ（Ｈｔ）からカルマンフィルタを用いて、前述の式２１に基づいて推定状態量ｘａｔを演算する。
【００６３】
また、ステップＳ２０８において、推定量演算部１０４は、カルマンゲインＫｔ、予測共分散Ｐ−ｔおよび観測量ｚｔ（Ｈｔ）からカルマンフィルタを用いて、前述の式２２に基づいて推定共分散Ｐｔを演算する。たとえば、推定量演算部１０４は、推定状態量と推定共分散とを所定の間隔で連続して演算する。
【００６４】
ステップＳ２０９において、プロセスノイズ共分散修正部１０６は、推定状態量ｘａｔを用いてプロセスノイズ共分散Ｑｔ−１を修正する。たとえば、プロセスノイズ共分散修正部１０６は、所定の間隔で連続してプロセスノイズ共分散Ｑｔを修正する。修正されたプロセスノイズ共分散を「Ｑｔ」とする。以下、プロセスノイズ共分散を修正する種々の態様について説明する。
【００６５】
図３は、本実施形態においてプロセスノイズ共分散を修正する処理における第一の態様の流れを示すフローチャートである。ステップＳ３０１において、プロセスノイズ共分散修正部１０６は、推定量演算部１０４が求めた推定状態量ｘａｔとプロセスノイズ共分散Ｑｔ−１とに基づいて、プロセスノイズ共分散Ｑｔ−１を修正し、修正されたプロセスノイズ共分散Ｑｔを演算する。
【００６６】
通常、車両に作用する状態量を推定するための演算の周期は、当該状態量が継時的に変化する期間に比べて十分に短いことから、状態量ｘａは、通常、一定と見なすことができることがある。この場合、今回の推定状態量ｘａｔが前回の状態量ｘａｔ−１に対して変化しているとすると、それは、ノイズに起因すると考えられる。プロセスノイズ共分散Ｑｔ−１は、状態量ｘａとプロセスノイズ共分散Ｑとを対応させたマップを車両に応じて予め作成しておき、当該マップを参照することによって適宜に修正することが可能である。
【００６７】
なお、当該第一の態様では、状態量推定装置１００は、重量成分抽出部１０５を含まなくてよい。車両に作用する状態量には車重が含まれるが、状態量推定装置１００によれば、カルマンフィルタの使用とプロセスノイズ共分散の修正とによって、重量成分抽出部１０５を有していなくても、車重を含む状態量を適切かつ迅速に推定することが可能である。
【００６８】
図４は、本実施形態においてプロセスノイズ共分散を修正する処理における第二の態様の流れを示すフローチャートである。当該第二の態様では、状態量推定装置１００は、重量成分抽出部１０５を含む。当該第二の態様では、プロセスノイズ共分散修正部１０６は、重量成分抽出部１０５が抽出した重量成分ｍを用いてプロセスノイズ共分散Ｑを修正する。
【００６９】
ステップＳ４０１において、重量成分抽出部１０５は、推定状態量ｘａｔから車両の重量成分ｍｔを、例えば所定の間隔で連続して抽出する。推定状態量ｘａは、前述した式１２における車両の重量成分ｍと定数Ｃの行列で表され、重量成分ｍは、当該行列中で表される値として決めることができる。
【００７０】
ステップＳ４０２において、プロセスノイズ共分散修正部１０６は、重量成分抽出部１０５が抽出した重量成分ｍｔとプロセスノイズ共分散Ｑｔ−１とに基づいてプロセスノイズ共分散Ｑｔ−１を修正し、修正されたプロセスノイズ共分散Ｑｔを演算する。修正されたプロセスノイズ共分散Ｑｔは、例えば、重量成分ｍｔとプロセスノイズ共分散Ｑｔ−１とを、重量成分ｍとプロセスノイズ共分散Ｑとについて上記のように車両ごとに予め作成されたマップを参照することによって求められる。上記の修正処理によれば、車両の重量成分の推定に対するプロセスノイズ共分散の不要な変動の影響をより抑制することが可能である。
【００７１】
図５は、本実施形態においてプロセスノイズ共分散を修正する処理における第三の態様の流れを示すフローチャートである。当該第三の態様では、重量成分抽出部１０５は、重量成分ｍの平均ｍａｖｅと、重量成分ｍの標準偏差ｍσとを演算する。そして、プロセスノイズ共分散修正部１０６は、重量成分抽出部１０５が演算した標準偏差ｍσを用いてプロセスノイズ共分散Ｑを修正する。
【００７２】
ステップＳ４０１において重量成分ｍｔを抽出すると、ステップＳ５０１において、重量成分抽出部１０５は、演算回数の所定の範囲において（例えば二つ前の演算回数までの重量成分ｍｔ−２を用いて）、抽出した重量成分ｍｔの平均値ｍａｖｅｔと標準偏差ｍσｔとを演算する。
【００７３】
ステップＳ５０２において、プロセスノイズ共分散修正部１０６は、当該標準偏差ｍσｔとプロセスノイズ共分散Ｑｔ−１とに基づいてプロセスノイズ共分散Ｑｔ−１を修正し、修正されたプロセスノイズ共分散Ｑｔを演算する。この場合、修正されたプロセスノイズ共分散Ｑｔは、例えば以下のようにして求められる。
【００７４】
プロセスノイズ共分散Ｑは、下記式２３の２×２行列で表される。下記式中、Ｑ１は車重ｍの推定に係るノイズの分散を表し、Ｑ２は路面勾配等の他の状態量Ｃの推定に係るノイズの分散を表す。プロセスノイズ共分散Ｑは、標準偏差ｍσに応じて当該Ｑ１およびＱ２の一方または両方を適宜に変更することにより修正することができる。
【００７５】
【数１７】
2020064015.tif 000019
【００７６】
より詳しくは、推定状態量ｘａｔの重量成分ｍｔの標準偏差ｍσｔが小さい場合では、ｍ推定に有るノイズの分散Ｑ１も小さいと言える。よって、Ｑ１を小さくする修正を行うか、またはＱ２を大きくする修正を行ってＱ１を相対的に小さくすることにより、プロセスノイズ共分散Ｑが適宜に修正され得る。
【００７７】
推定状態量ｘａｔの重量成分ｍｔの標準偏差ｍσｔが大きい場合では、ｍ推定に有るノイズの分散Ｑ１も大きいと言える。よって、Ｑ１を大きくする修正を行うか、またはＱ２を小さくする修正を行ってＱ１を相対的に大きくすることにより、プロセスノイズ共分散Ｑが適宜に修正され得る。
【００７８】
標準偏差ｍσの大小は、複数の閾値によって判断することができ、当該閾値に応じて上記のようにＱ１またはＱ２を段階的に修正することができる。たとえば、ｍσに対応するＱ１またはＱ２をマップから読み出すことにより、ｍσに応じた適当なＱ１またはＱ２の値を決めることができ、それに基づきＱ１またはＱ２を上記のように修正することが可能である。このような修正方法によれば、重量成分ｍのそれまでの挙動に応じてプロセスノイズ共分散Ｑが修正されるため、より安定した状態量の推定が可能となる。
【００７９】
図６は、本実施形態においてプロセスノイズ共分散を修正する処理における第四の態様の流れを示すフローチャートである。当該第四の態様では、重量成分抽出部１０５は、安定フラグと重量決定部とを含む。状態量の推定処理における安定フラグの初期状態はオフである。安定フラグには、オンの場合に１が入力され、オフの場合には０が入力される。重量決定部は、推定状態量ｘａの演算回数ｔ、標準偏差ｍσの値、および、安定フラグのオンオフ、からなる群から選ばれる一以上を表す情報に応じて重量成分ｍを決定する。
【００８０】
ステップＳ４０１において、重量成分抽出部１０５は、推定状態量ｘａｔから重量成分ｍｔを抽出する。
【００８１】
ステップＳ６０１において、重量決定部は、演算回数ｔが所定の回数(n)を超えているか否かを判定する。ステップＳ６０１において演算回数ｔが所定の回数(n)を超えている場合には、ステップＳ５０２において、重量成分抽出部１０５は、重量成分ｍｔの平均ｍａｖｅｔおよび標準偏差ｍσｔを演算する。
【００８２】
ステップＳ６０２において、重量決定部は、１０−６を所定の閾値として、標準偏差ｍσｔが当該閾値未満か否かを判定する。
【００８３】
ステップＳ６０２において標準偏差ｍσｔが当該閾値未満である場合には、ステップＳ６０３において、重量決定部は、重量成分の平均ｍａｖｅｔを抽出した重量成分ｍｔと決定する。そして安定フラグをオンにする。
【００８４】
ステップＳ６０２において標準偏差ｍσｔが当該閾値未満である場合には、ステップＳ６０４において、重量決定部は、安定フラグがオンとなっているか否かを判定する。
【００８５】
ステップＳ６０４において安定フラグがオンとなっている場合には、ステップＳ６０５において、重量決定部は、前回以前の前記安定フラグをオンにしたときの（例えば一つ前の演算における）重量成分ｍｔ（例えばｍｔ−１）を今回抽出した重量成分ｍｔと決定する。
【００８６】
ステップＳ６０４において安定フラグがオフとなっている場合には、ステップＳ６０６において、重量決定部は、最新の抽出した重量成分ｍｔをそのまま重量成分と決定する。また、ステップＳ６０１において演算回数ｔが所定の回数(n)以下である場合には、重量決定部は、ステップＳ６０６に進み、最新の抽出した重量成分ｍｔをそのまま重量成分と決定する。
【００８７】
なお、ステップＳ６０３において重量成分ｍｔを決定した場合には、プロセスノイズ共分散修正部１０６は、ステップＳ５０２において、ステップＳ５０１で求めた重量成分の標準偏差ｍσｔに基づいてプロセスノイズ共分散Ｑｔ−１を修正してもよい。
【００８８】
ステップＳ６０３において標準偏差からステップＳ６０６において重量成分抽出部１０５が重量成分ｍｔを決定すると、プロセスノイズ共分散修正部１０６は、ステップＳ４０２において当該重量成分を用いてプロセスノイズ共分散Ｑｔ−１を修正し、修正されたプロセスノイズ共分散Ｑｔを演算する。図６のフローによれば、場合に応じて今回または前回以前の演算結果が適宜に利用されることから、このような状態量の推定は、十分に高い精度で迅速に車両の状態量を推定する観点から好適である。
【００８９】
そして、図２のフローのステップＳ２１０において、推定量演算部１０４は、推定状態量ｘａｔおよび推定共分散Ｐｔをデータ格納部１０１に格納し、プロセスノイズ共分散修正部１０６は、修正したプロセスノイズ共分散Ｑｔをデータ格納部１０１に格納する。このように、推定状態量ｘａｔが演算回数ｔにおける状態量ｘａｔとして、推定共分散Ｐｔが演算回数ｔにおける状態共分散Ｐｔとして、また、修正されたプロセスノイズ共分散Ｑｔが演算回数ｔにおけるプロセスノイズ共分散Ｑｔとして、データ格納部１０１に書き込まれる。こうして、状態量、状態共分散およびプロセスノイズ共分散が、例えば所定の間隔で連続してデータ格納部１０１に書き込まれる。データ格納部１０１に書き込まれたこれらのデータは、演算回数が更新されると、新たな演算において前回の演算による結果として、新たな状態量の推定において上記のように用いられる。
【００９０】
上記の一連の処理は、車両が、推定すべき状態量が変わり得る状態になった場合に、終了させてよい。当該処理は、当該状態量は車重成分を含むことから、車両の車重が変わり得る状態になった場合に終了させることができる。車両の車重が変わり得る場合の例には、車両が停止した場合、車両のドアが開閉された場合、および、車両の燃料供給口が開閉された場合、が含まれる。
【００９１】
本実施形態に係る状態量推定装置は、車重を含む状態量の推定にカルマンフィルタを用いる。よって、推定の条件に応じて観測量を破棄する必要がなく、実質的に全ての、獲得した観測量を演算に用い、状態量の推定に利用することが可能である。よって、車両の状態量の推定値を迅速かつ十分に取得することが可能である。
【００９２】
また、本実施形態に係る状態量推定装置では、プロセスノイズ共分散が修正される。たとえば、車両の車重は、種々のセンサ値に応じて推定する場合では、車両の走行状態に応じて変化することがある。しかしながら、車重の推定を所定の期間（例えば車両が走行し始めてから停止するまで等）ごとに行う場合では、車重は実質的に一定である。本実施形態では、プロセスノイズ共分散を修正することで、車両の走行状態に起因して車重を変化させるパラメータを実質的に無視することができ、車両の状態量の推定を安定して実行することが可能となる。よって、車両の状態量を十分に高い精度で迅速に推定することが可能である。
【００９３】
本実施形態に係る状態量推定装置は、車重を含む状態量が推定される。当該状態量推定装置は、車両に用いられる制御装置であって、車重の推定値を用いる制御装置に適用できる。この場合、状態量推定装置は、重量成分を含む状態量を外部に出力する出力部をさらに有することが可能である。
【００９４】
当該出力部は、車重を含む状態量を外部に出力可能な範囲において適宜に選ばれ得る。車重を出力する場合では、出力部は、重量成分ｍを抽出する重量成分抽出部であることが好ましい。このように状態量推定装置を、車重に基づいて車両を制御するための制御装置に適用する場合では、重量成分抽出部１０５は、抽出した重量成分ｍｔを車両の現在の重量として出力してもよい。この場合、車重の推定結果を車両の走行状態の制御に利用する観点から好適である。
【００９５】
本実施形態に係る状態量推定装置は、前述したように、センサ値を破棄することなく車両の状態量の推定に用いることができ、かつ、プロセスノイズ共分散の修正によって、当該状態量の安定した推定が可能である。よって、状態量の推定値は、適度なロバスト性を有し、車両の走行開始後の短期間で適切な値に収束する。この推定値を用いて車両の走行状態を制御すると、車両の走行状態をより好適に制御することが可能となる。
【００９６】
以下、本実施形態に係る状態量推定装置を車両における懸架装置のための制御装置に適用する態様を説明する。
【００９７】
〔実施形態２〕
本発明の他の実施形態について、以下に説明する。なお、説明の便宜上、上記実施形態にて説明した部材と同じ機能を有する部材については、同じ符号を付記し、その説明を繰り返さない。
【００９８】
当該車両は、車両を走行させる各種装備の他、プロセッサ、センサ類、および記憶媒体を有する。当該プロセッサは、前述した予測量演算部、取得部、カルマンゲイン演算部、推定量演算部およびプロセスノイズ共分散修正部を含む。センサ類は、前述した前後Ｇセンサ、車輪軸トルクおよび車速センサを含む。記憶媒体は、前述したデータ格納部に対応しており、例えば磁気ディスクである。
【００９９】
図７は、上記の状態量推定装置を有する車両の構成の一例を模式的に示す図である。図７に示されるように、車両９００は、懸架装置（サスペンション）１５０、車体２００、車輪３００、車速（Ｖ）を検出する車速センサ４５０、エンジン５００およびＥＣＵ（Electronic Control Unit）６００を備えている。ＥＣＵ６００は、前述したプロセッサに該当している。
【０１００】
なお、符号中のアルファベットＡ〜Ｅは、それぞれ、車両９００における位置を表している。Ａは、車両９００の左前の位置を表し、Ｂは、車両９００の右前の位置を表し、Ｃは、車両９００の左後ろを表し、Ｄは、車両９００の右後ろを表し、Ｅは、車両９００の後ろを表している。
【０１０１】
また、車両９００は、車両９００の前後方向の加速度（以下、「前後Ｇ」とも言う）を検出する前後Ｇセンサ３４０（前述の「Ｇｓｅｎｓ」を取得する）、および、エンジン５００が発生させるトルク（車輪トルク（τ））を推定する車輪トルクセンサ５１０を備えている。上記車輪トルクは、前述したように車両９００の車輪に掛けられているトルクであり、車両９００の駆動力に該当する。
【０１０２】
さらに、車両９００は、不図示の燃料タンク中の燃料（ガソリン）の量を検出するフューエルセンサ３６０、車両９００のそれぞれのドアの開閉を検出するドア開閉センサ５５０、および車両９００のギアの接続位置を検出するギアポジションセンサ５４０を備えている。なお、ドア開閉センサ５５０Ｅは、車両９００のトランクのドア（バックドア）の開閉を検出する。
【０１０３】
各種センサの出力値のＥＣＵ６００への供給、および、ＥＣＵ６００から各部への制御信号の伝達は、ＣＡＮ（Controller Area Network）３７０を介して行われる。各センサは、後述の状態量の推定のために新たに設けられてもよいが、コストの面から、車両９００に既存のセンサであることが好ましい。
【０１０４】
記憶媒体には、状態量の推定に要する種々の情報が記憶されている。当該情報の例には、車輪半径、Ｃｄ値（空気抵抗係数）、ならびに、状態量、状態共分散、状態標準偏差およびプロセスノイズの共分散の初期値、が含まれる。初期値は、テスト等により事前に決定された、初期値として適切な数値である。
【０１０５】
たとえば、車両９００の状態量として車重を推定する場合では、当該車重の初期値は、車両９００の仮の総重量であってよい。当該仮の総重量は、例えば、車両９００の乾燥重量と、車両への搭載物の総重量とまたはその見込み量との和であってもよい。また、仮の総重量は、例えば、車両９００の乾燥重量に、フューエルセンサ３６０が検出したガソリン残量を加算した値であってもよい。また、当該仮の総重量は、例えば、車両９００の乾燥重量に、ドア開閉センサ５５０が検出した、開いたフロアドアの数に人間一人分相当の重量を乗じた値と、ドア開閉センサ５５０がバックドアの開閉を検出した場合には所定の重量との和であってよい。また、当該仮の総重量は、車両９００におけるシートベルトの着用の有無をさらに参照した値であってもよい。シートベルトの着用を参照することにより、車両９００に搭乗している人間の数をより高い精度で把握される。搭乗者数の高い精度での把握は、例えば車両９００の重量の下限値の設定に有効であり、上記仮の総重量の信頼性をより一層高める観点から好ましい。
【０１０６】
車両９００の状態量の推定は、実施形態１で前述したように実施される。車両９００における状態量の推定プロセスは、基本的に所定の間隔で連続して行われる。状態量の推定値が安定する場合には、状態量の推定の間隔をより長い間隔で間欠的に行ってもよい。ＥＣＵ６００が推定した状態量は、当該状態量の推定値を必要とする車両９００の他の装置において、制御などの目的で使用される。他の装置は、推定状態量を推定量演算部から受け取ってもよいし、データ格納部から読み込んでもよい。
【０１０７】
たとえば、車重は、一般に、車両の重量と車載物の総重量とによって決まり、車載物の重量は、通常、走行中には実質的には変動しない。したがって、車両が走行し始めてから停止するまでの間、車重は実質的には一定と考えることが可能である。本実施形態では、前述した状態量推定装置によって車両の重量成分を含む状態量が推定されるため、車両が走行し始めると、観測量に基づいて車重の推定値が迅速に得られる。車重は、車重が変化し得るように車両の走行が中断するまで（交差点または目的地に到着することによる車両の停止、停止時におけるドアまたは燃料キャップの開閉、など）は、実質的に一定である。よって、このような走行時に車重の推定値を変化させる観測量および演算結果はプロセスのノイズと見なすことができる。よって、本実施形態におけるＣＰＵによれば、走行状態にある車両の車重は迅速かつ簡易に推定され、車重の推定を安定して実行することが可能となる。
【０１０８】
本実施形態において、車重が空の場合の重量（空虚重量）以下等、走行中の車両の車重として適切でない値が車重の推定値として算出された場合には、推定結果を破棄してもよい。また、前述したように、車両における車重の推定は、通常は連続して行われるが、十分に信頼可能な値を獲得している場合には、車重の推定値を破棄してもよいし、あるいは車重の推定を断続的（間欠的）に実行してもよい。
【０１０９】
前述した実施形態では、状態量推定装置における各種値を演算する前述の演算部または修正部は、算出した値をデータ格納部に書き込むが、本発明では、これらの値のデータのそれぞれをデータ書き込み部に書き込むための書き込み部をＣＰＵがさらに含んでいてもよい。
【０１１０】
上記の状態量推定装置による状態量の推定は、当該状態量が変化する可能性がある車両の状態に応じてリセットすることが可能である。たとえば、車両の車重を推定する場合では、車両のキーポジションがオフの位置となることや、ドアが開かれることによって、推定状態量の初期値をリセットし、状態量の推定を新たに開始してもよい。
【０１１１】
なお、上記の実施形態では、四輪自動車を例に説明したが、本発明において車両はそれに限定されない。たとえば、車両は、バイクであってもよく、また鉄道車両であってもよい。
【０１１２】
本実施形態では、前述したように、車両の走行開始から短時間で車重の安定した推定値を取得することが可能である。よって、車重の推定値に基づいて車両の走行状態を制御する場合に、当該制御をより安定して好適に実行することが可能である。
【０１１３】
〔実施形態３〕
本実施形態は、ＣＰＵが前述した重量成分抽出部をさらに含む。それ以外は、本実施形態は、前述した実施形態２と同様である。
【０１１４】
本実施形態では、重量成分抽出部において車両の重量成分が抽出される。本実施形態によれば、車両における車重の推定値を用いる制御を、より高い精度でより迅速に実行することが可能となる。車重の推定値を必要とする車両の他の装置は、重量成分抽出部で決定された重量成分を車重の推定値として受け取ってもよい。重量成分抽出部は、決定した重量成分を推定車重としてデータ格納部にさらに書き込んでもよい。データ格納部への書き込みは、書き込み部を介して実行されてもよい。
【０１１５】
〔実施形態４〕
本実施形態に係る状態量推定装置で推定された車重は、懸架装置１５０の制御に用いることができる。車両９００における懸架装置１５０の制御装置は、車両９００に作用する状態量を推定して、懸架装置１５０の減衰力を状態量に応じて制御するように構成される。本実施形は、本実施形態に係る状態量推定装置による状態量または車重の推定値を用いる以外は、車重の推定値に応じて懸架装置１５０の減衰力を制御する公知の方法によって実施することが可能である。
【０１１６】
懸架装置１５０は、例えば、車両における車体と車輪との間に介在するアブソーバと、アブソーバのストロークに伴い伸縮するように配置されているスプリングとを有する。アブソーバは、シリンダと、シリンダ内を二室に仕切るとともに摺動可能なピストンと、ピストンに固定されたピストンロッドと、二室を連通する連通路と、連通路を開閉自在なソレノイドバルブとを有している。ピストンで仕切られるいずれの室にも作動油が満たされている。スプリングは、ピストンロッドの外周を囲むように配置されており、シリンダの端部とピストンロッドの端部とで支持されている。ＥＣＵ６００は、車重の推定量に応じて、車重の推定量が大きい程、懸架装置１５０の減衰力を大きくするように、ソレノイドバルブの開度を調整させる。
【０１１７】
なお、懸架装置１５０の構造は限定されない。たとえば、懸架装置１５０におけるソレノイドの位置は限定されない。懸架装置１５０は、例えば、ピストンイン型であってもよいし、外筒接続型であってもよい。また、懸架装置１５０の種類は、特に限定されず、例えばストラット式であってもよいし、ダブルウィッシュボーン式であってもよい。このように種々の構造または種類の懸架装置を懸架装置１５０に採用することができる。また、懸架装置１５０における減衰力の調整方法も限定されず、上記のようにアブソーバにおける減衰力を調整する方法であってもよいし、上記スプリングのスプリングレートを可変させる(自動可変プリロード)方法であってもよいし、これらの両方であってもよい。
【０１１８】
本実施形態によれば、走行開始時に迅速に車両の状態量の安定した推定値が得られる。このように早期に安定して得られる、車重を含む当該状態量の推定値を用いて懸架装置を制御することから、車両の走行状態に応じた懸架装置の制御をより好適に実行することが可能となる。
【０１１９】
〔ソフトウェアによる実現例〕
状態量推定装置の制御ブロック（特に予測量演算部１０２からプロセスノイズ共分散修正部１０６までの各部）は、集積回路（ＩＣチップ）等に形成された論理回路（ハードウェア）によって実現してもよいし、ソフトウェアによって実現してもよい。
【０１２０】
後者の場合、状態量制御装置は、各機能を実現するソフトウェアであるプログラムの命令を実行するコンピュータを備えている。このコンピュータは、例えば１つ以上のプロセッサを備えていると共に、上記プログラムを記憶したコンピュータ読み取り可能な記録媒体を備えている。そして、上記コンピュータにおいて、上記プロセッサが上記プログラムを上記記録媒体から読み取って実行することにより、本発明の目的が達成される。上記プロセッサとしては、例えばＣＰＵ（Central Processing Unit）を用いることができる。上記記録媒体としては、「一時的でない有形の媒体」、例えば、ＲＯＭ（Read Only Memory）等の他、テープ、ディスク、カード、半導体メモリ、プログラマブルな論理回路などを用いることができる。また、上記プログラムを展開するＲＡＭ（Random Access Memory）などをさらに備えていてもよい。また、上記プログラムは、該プログラムを伝送可能な任意の伝送媒体（通信ネットワークや放送波等）を介して上記コンピュータに供給されてもよい。なお、本発明の一態様は、上記プログラムが電子的な伝送によって具現化された、搬送波に埋め込まれたデータ信号の形態でも実現され得る。
【０１２１】
〔まとめ〕
以上の説明から明らかなように、本発明の実施の形態に係る状態量推定装置（ＥＣＵ６００）は、車両（９００）に作用する状態量、状態量の共分散である状態共分散、および、演算ノイズの共分散であるプロセスノイズ共分散、を格納するデータ格納部（１０１）と、状態量から予測状態量を演算し、状態共分散およびプロセスノイズ共分散から予測状態量の共分散である予測共分散を演算する予測量演算部（１０２）と、車両（９００）のセンサ値を取得する取得部（１０７）と、センサ値および予測共分散を用いてカルマンゲインを演算するカルマンゲイン演算部（１０３）と、カルマンゲイン、予測状態量およびセンサ値、からカルマンフィルタを用いて、推定状態量を演算し、カルマンゲイン、予測共分散およびセンサ値、からカルマンフィルタを用いて、推定共分散を演算する推定量演算部（１０４）と、推定状態量を用いてプロセスノイズ共分散を修正するプロセスノイズ共分散修正部（１０６）と、を備える。そして、推定状態量が状態量として、推定共分散が状態共分散として、修正されたプロセスノイズ共分散がプロセスノイズ共分散として、データ格納部（１０１）に書き込まれる。当該構成によれば、カルマンフィルタを用いて演算するため、観測量が条件に合わないものを含んでいる場合でも観測量を破棄する必要がない。よって、センサ類から求められた観測量を演算に利用することが可能である。また、プロセスノイズ共分散を修正するため、種々異なる推定環境においても安定して車両の状態量を精度よく、また迅速に推定することが可能である。
【０１２２】
上記の状態量推定装置は、推定状態量から車両の重量成分を抽出する、重量成分抽出部（１０５）をさらに備えてもよい。このような構成によれば、車重をより精密に推定することが可能である。
【０１２３】
また、上記の状態量推定装置において、重量成分抽出部（１０５）は、重量成分の平均と重量成分の標準偏差とを演算してもよい。このような構成によれば、車重を十分な精度で簡素に推定する観点からより効果的である。
【０１２４】
また、上記の状態量推定装置において、重量成分抽出部（１０５）は、初期状態をオフとする安定フラグと、推定状態量の演算回数、標準偏差の値、および、安定フラグのオンオフ、からなる群から選ばれる一以上を表す情報に応じて重量成分を決定する重量決定部とを含んでいてよい。そして、重量決定部は、推定状態量の演算回数が所定回数以下の場合には、最新の抽出した重量成分を抽出した重量成分と決定し、推定状態量の演算回数が所定回数を超え、かつ、標準偏差が所定の閾値以下の場合には、平均を抽出した重量成分と決定して安定フラグをオンにし、安定フラグがオンであるが、推定状態量の演算回数が所定回数を超え、かつ標準偏差が所定の閾値を超えている場合には、前回以前の安定フラグをオンにしたときの抽出した重量成分を抽出した重量成分と決定し、推定状態量の演算回数が所定回数を超え、標準偏差が所定の閾値を超え、かつ、安定フラグがオフの場合には、最新の抽出した重量成分を抽出した重量成分と決定してよい。このような構成によれば、車重を十分な精度で、また簡素かつ迅速に推定する観点からより一層効果的である。
【０１２５】
また、上記の状態量推定装置において、プロセスノイズ共分散修正部（１０６）は、重量成分抽出部（１０５）が抽出した重量成分を用いてプロセスノイズ共分散を修正してもよい。このような構成によれば、車重を十分な精度で簡素に推定する観点からより効果的である。
【０１２６】
また、上記の状態量推定装置において、プロセスノイズ共分散修正部（１０６）は、重量成分抽出部（１０５）が演算した標準偏差を用いてプロセスノイズ共分散を修正してもよい。このような構成によれば、車重を十分な精度で簡素に推定する観点からより一層効果的である。
【０１２７】
また、上記の状態量推定装置において、重量成分抽出部（１０５）は、抽出した重量成分を車両（９００）の現在の重量として出力してもよい。このような構成によれば、車重の推定結果を車両の走行状態の制御に利用する観点から好適である。
【０１２８】
また、本発明の実施の形態に係る制御装置は、懸架装置（１０５）を有する車両（９００）に作用する状態量を推定して、懸架装置（１５０）の減衰力を状態量に応じて制御する制御装置である。当該制御装置は、車両に作用する状態量、状態量の共分散である状態共分散、および、演算ノイズの共分散であるプロセスノイズ共分散を格納するデータ格納部（１０１）と、状態量から予測状態量を演算し、状態共分散およびプロセスノイズ共分散から予測状態量の共分散である予測共分散を演算する予測量演算部（１０２）と、車両のセンサ値を取得する取得部（１０７）と、センサ値および予測共分散を用いてカルマンゲインを演算するカルマンゲイン演算部（１０３）と、カルマンゲイン、予測状態量およびセンサ値、からカルマンフィルタを用いて、推定状態量を演算し、カルマンゲイン、予測共分散およびセンサ値、からカルマンフィルタを用いて、推定共分散を演算する推定量演算部（１０４）と、推定状態量を用いてプロセスノイズ共分散を修正し、するプロセスノイズ共分散修正部（１０６）とを備え、推定状態量が状態量として、推定共分散が状態共分散として、データ格納部（１０１）に書き込まれる。このような構成によれば、車重を含む状態量を十分な精度で迅速に推定することが可能であるため、車両のより快適な走行を実現するように車両の懸架装置を制御することが可能となる。
【０１２９】
また、本発明の実施の形態に係る状態量推定方法は、車両（９００）に作用する状態量、状態量の共分散である状態共分散、および、演算ノイズの共分散であるプロセスノイズ共分散、をデータ格納部に格納するステップと、状態量から予測状態量を演算し、状態共分散およびプロセスノイズ共分散から予測状態量の共分散である予測共分散を演算するステップと、車両（９００）のセンサ値を取得するステップと、センサ値および予測共分散を用いてカルマンゲインを演算するステップと、カルマンゲイン、予測状態量およびセンサ値、からカルマンフィルタを用いて、推定状態量を演算するステップと、カルマンゲイン、予測共分散およびセンサ値、からカルマンフィルタを用いて、推定共分散を演算するステップと、推定状態量を用いてプロセスノイズ共分散を修正するステップと、推定状態量を状態量とし、推定共分散を状態共分散として、データ格納部に書き込むステップとを含む。このような構成によれば、前述したように、センサ類から求められた観測量を演算に利用することが可能であり、また求めたい状態量に実質的には関係しないプロセスノイズ共分散の変動が求めたい状態量の推定に実質的に影響を及ぼすことを抑制することが可能である。よって、車両の状態量を精度よく、また迅速に推定することが可能である。
【０１３０】
本発明は上述した各実施形態に限定されるものではなく、請求項に示した範囲で種々の変更が可能であり、例えば、Ｕｎｓｃｅｎｔｅｄ  カルマンフィルタのような代替容易に考えられる、既知の技術的手段を適宜組み合わせて得られる実施形態についても本発明の技術的範囲に含まれる。
【符号の説明】
【０１３１】
１００  状態量推定装置
１０１  データ格納部
１０２  予測量演算部
１０３  カルマンゲイン演算部
１０４  推定量演算部
１０５  重量成分抽出部
１０６  プロセスノイズ共分散修正部
１０７  取得部
１０８  センサ類
１５０  懸架装置
２００  車体
３００  車輪
３４０  前後Ｇセンサ
３６０  フューエルセンサ
４５０  車速センサ
５００  エンジン
５１０  車輪トルクセンサ
５４０  ギアポジションセンサ
５５０  ドア開閉センサ
６００  ＥＣＵ（車重推定装置）
９００  車両
【図１】
2020064015.tif 000020
【図２】
2020064015.tif 000021
【図３】
2020064015.tif 000022
【図４】
2020064015.tif 000023
【図５】
2020064015.tif 000024
【図６】
2020064015.tif 000025
【図７】
2020064015.tif 000026
【手続補正書】
【提出日】令和1年7月23日(2019.7.23)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
少なくとも車重を含む車両に作用する状態量、前記状態量の共分散である状態共分散、および、演算ノイズの共分散であるプロセスノイズ共分散、を格納するデータ格納部と、
前記状態量から予測状態量を演算し、前記状態共分散および前記プロセスノイズ共分散から前記予測状態量の共分散である予測共分散を演算する予測量演算部と、
車両のセンサ値を取得する取得部と、
前記センサ値および前記予測共分散を用いてカルマンゲインを演算するカルマンゲイン演算部と、
前記カルマンゲイン、前記予測状態量および前記センサ値、からカルマンフィルタを用いて、推定状態量を演算し、前記カルマンゲイン、前記予測共分散および前記センサ値、からカルマンフィルタを用いて、推定共分散を演算する推定量演算部と、
前記推定状態量から前記車両の重量成分を抽出する重量成分抽出部と、
前記重量成分抽出部が抽出した前記重量成分を用いて前記プロセスノイズ共分散を修正するプロセスノイズ共分散修正部と、を備え、
前記推定状態量が前記状態量として、前記推定共分散が前記状態共分散として、修正された前記プロセスノイズ共分散が前記プロセスノイズ共分散として、前記データ格納部に書き込まれる、状態量推定装置。
【請求項２】
前記重量成分抽出部は、前記重量成分の平均と、前記重量成分の標準偏差と、を演算する、請求項１に記載の状態量推定装置。
【請求項３】
前記状態量は、車両が走行している路面の勾配成分を更に含み、
前記プロセスノイズ共分散は、少なくとも前記車重の推定に係る車重ノイズおよび前記勾配成分の推定に係る勾配成分ノイズを含み、
前記プロセスノイズ共分散修正部は、前記車重の前記標準偏差に、前記車重ノイズが直接的にまたは間接的に比例するように前記プロセスノイズ共分散を修正する、
請求項２に記載の状態量推定装置。
【請求項４】
前記重量成分抽出部は、前記推定状態量の演算回数、および、前記標準偏差の値、からなる群から選ばれる一以上を表す情報に応じて前記重量成分を決定する重量決定部を含む、請求項２または３に記載の状態量推定装置。
【請求項５】
前記重量成分抽出部は、初期状態をオフとする安定フラグをさらに含み、
前記重量決定部の前記群は前記安定フラグのオンオフをさらに含み、
前記重量決定部は、
前記推定状態量の演算回数が所定回数以下の場合には、最新の抽出した前記重量成分を抽出した前記重量成分と決定し、
前記推定状態量の演算回数が前記所定回数を超え、かつ、前記標準偏差が所定の閾値未満の場合には、前記平均を抽出した前記重量成分と決定して前記安定フラグをオンにし、
前記安定フラグがオンであるが、前記推定状態量の演算回数が前記所定回数を超え、かつ前記標準偏差が所定の閾値以上である場合には、前回以前の前記安定フラグをオンにしたときに抽出した前記重量成分を抽出した前記重量成分と決定し、
前記推定状態量の演算回数が前記所定回数を超え、前記標準偏差が所定の閾値以上であり、かつ、前記安定フラグがオフの場合には、最新の抽出した前記重量成分を抽出した前記重量成分と決定する、
請求項４に記載の状態量推定装置。
【請求項６】
前記プロセスノイズ共分散修正部は、前記重量成分抽出部が演算した前記標準偏差を用いて前記プロセスノイズ共分散を修正する、請求項２、４または５に記載の状態量推定装置。
【請求項７】
前記重量成分抽出部は、抽出した前記重量成分を前記車両の現在の重量として出力する、請求項１〜５のいずれか一項に記載の状態量推定装置。
【請求項８】
懸架装置を有する車両に作用する状態量を推定して、前記懸架装置の減衰力を前記状態量に応じて制御する制御装置であって、
少なくとも車重を含む車両に作用する状態量、前記状態量の共分散である状態共分散、および、演算ノイズの共分散であるプロセスノイズ共分散、を格納するデータ格納部と、
前記状態量から予測状態量を演算し、前記状態共分散および前記プロセスノイズ共分散から前記予測状態量の共分散である予測共分散を演算する予測量演算部と、
車両のセンサ値を取得する取得部と、
前記センサ値および前記予測共分散を用いてカルマンゲインを演算するカルマンゲイン演算部と、
前記カルマンゲイン、前記予測状態量および前記センサ値、からカルマンフィルタを用いて、推定状態量を演算し、前記カルマンゲイン、前記予測共分散および前記センサ値、からカルマンフィルタを用いて、推定共分散を演算する推定量演算部と、
前記推定状態量から前記車両の重量成分を抽出する重量成分抽出部と、
前記重量成分抽出部が抽出した前記重量成分を用いて前記プロセスノイズ共分散を修正し、するプロセスノイズ共分散修正部と、を備え、
前記推定状態量が前記状態量として、前記推定共分散が前記状態共分散として、修正された前記プロセスノイズ共分散が前記プロセスノイズ共分散として、前記データ格納部に書き込まれる、制御装置。
【請求項９】
少なくとも車重を含む車両に作用する状態量、前記状態量の共分散である状態共分散、および、演算ノイズの共分散であるプロセスノイズ共分散、をデータ格納部に格納するステップと、
前記状態量から予測状態量を演算し、前記状態共分散および前記プロセスノイズ共分散から前記予測状態量の共分散である予測共分散を演算するステップと、
車両のセンサ値を取得するステップと、
前記センサ値および前記予測共分散を用いてカルマンゲインを演算するステップと、
前記カルマンゲイン、前記予測状態量および前記センサ値、からカルマンフィルタを用いて、推定状態量を演算するステップと、
前記カルマンゲイン、前記予測共分散および前記センサ値、からカルマンフィルタを用いて、推定共分散を演算するステップと、
前記推定状態量から前記車両の重量成分を抽出するステップと、
前記車両の重量成分を用いて前記プロセスノイズ共分散を修正するステップと、
前記推定状態量を前記状態量とし、前記推定共分散を前記状態共分散として、修正された前記プロセスノイズ共分散を前記プロセスノイズ共分散として、前記データ格納部に書き込むステップと、
を含む、状態量推定方法。