(57)【要約】  車両の現在の荷重を推定するための方法、装置、及びコンピュータ・プログラム製品であって、フィルタ・バンクが異なる重量区分に対するフィルタを含み、各フィルタが、車両の現在の質量を推定するための車両モデルを実装する。車両の現在の走行状況を示す車両データと、それぞれの重量区分に固有のフィルタ・パラメータとに基づいて、各フィルタは、車両の現在の荷重のフィルタ固有の推定として荷重推定値を提供する。
2011505566.tif 000002
【特許請求の範囲】
【請求項１】
車両の現在の荷重を推定する方法であって、
少なくとも２つの重量区分を定義する段階と、
前記少なくとも２つの重量区分の各々に対して、各々が前記車両の現在の質量を推定するための車両モデルを実装するフィルタを含むフィルタ・バンクを設ける段階と、
各フィルタに対して、前記車両の現在の走行状況を示す少なくとも１つの車両データ、及び前記それぞれの重量区分に固有の少なくとも１つの現在のフィルタ・パラメータを供給する段階と、
各フィルタを用いて、前記車両の現在の荷重のフィルタ固有の推定として荷重推定値を特定する段階と、
各荷重推定値に対して、前記車両の現在の荷重が前記それぞれの重量区分に属すると仮定できるか否かの現在指標値を特定する段階と、
全ての荷重推定値の中で、前記現在の走行状況における前記車両の現在の荷重の推定として最良の現在指標値を有する荷重推定値を選択する段階、及び／又は前記重量区分における前記現在指標値を重み付けして、前記現在の走行状況における前記車両の現在の荷重についての荷重推定として用いられる全体荷重推定値を前記フィルタ・バンクにおいて得る段階と、
を含む方法。
【請求項２】
各フィルタに対して少なくとも１つの車両データ及び少なくとも１つの現在のフィルタ・パラメータを供給する前記段階、各フィルタを用いて荷重推定値を特定する前記段階、前記各荷重推定値について現在指標値を特定する前記段階、及び全ての荷重推定値の中で前記最良の現在指標値が特定された前記質量推定を選択する前記段階が繰り返される、
請求項１に記載の方法。
【請求項３】
前記段階の繰り返しは、所定の時間間隔で及び／又は前記車両の実際の質量が推定されることになる各時点で実施される、
請求項２に記載の方法。
【請求項４】
各フィルタが、モデルベース・フィルタ、カルマン・フィルタ、拡張カルマン・フィルタ、無香カルマン・フィルタ、制約条件付き拡張カルマン・フィルタ、粒子フィルタ、シグマ点フィルタ、点質量フィルタ、格子ベース・フィルタのうちの少なくとも１つを含む、
請求項１から請求項３のうちの一項に記載の方法。
【請求項５】
前記車両モデルが、
ギアボックス効率（ηtm）パラメータと、
エンジン・トルク（Ｔe）パラメータと、
ギアボックス減衰定数（ζ）パラメータと、
エンジン角速度（ωe）パラメータと、
エンジン慣性（Ｊe）パラメータと、
ホイール角速度（ωw）パラメータと、
ホイール角加速度（ωwドット）パラメータと、
ギア比（ｉg）パラメータと、
最終駆動ギア比（ｉgf）パラメータと、
ホイール慣性（Ｊe）パラメータと、
空気密度（ρ）パラメータと、
空力抗力係数（Ｃd）パラメータと、
車両空力前面面積（Ａf）パラメータと、
車速（ｖ）パラメータと、
前記それぞれのフィルタからの過去に特定された荷重推定値（ｍヤマ）パラメータと、
前記車両の現在の荷重が前記それぞれのフィルタからのそれぞれの重量区分に属すると仮定できるか否かについて過去に特定された指標値パラメータと、
重力定数（ｇ）パラメータと、
転がり抵抗係数（Ｃr）パラメータと、
ホイール半径（ｒ）パラメータと、
うちのの少なくとも１つを含む、
請求項１から請求項４のうちの一項に記載の方法。
【請求項６】
前記現在の走行状況を示す少なくとも１つの車両データを供給する前記段階が、
現在の車両加速度を示す加速度センサ値（ｙ）を供給する段階と、
前記車両のエンジンによって発生した現在のトルクを示すエンジン・トルク値（Ｔe）を供給する段階と、
前記車両の少なくとも１つのホイールの現在の角速度を示す少なくとも１つのホイール角速度値（ｗw）を供給する段階と、
前記車両のエンジンの現在の角速度を示すエンジン角速度値（ωe）を供給する段階と、
現在の最終駆動比を示す最終駆動比値（ｉfd）を供給する段階と、
のうちの少なくとも１つを含む、
請求項１から請求項５のうちの一項に記載の方法。
【請求項７】
前記それぞれの重量区分に固有の少なくとも１つの現在のフィルタ・パラメータを供給する前記段階が、
前記それぞれのフィルタの過去の荷重推定値パラメータと、
前記それぞれのフィルタの過去の平均値パラメータと、
前記それぞれのフィルタの過去の共分散パラメータと、
前記それぞれのフィルタの過去の確率パラメータと、
のうちの少なくとも１つを供給する段階を含む、
請求項１から請求項６のうちの一項に記載の方法。
【請求項８】
前記過去の共分散が、前記それぞれの車両モデルのモデル誤差の共分散、及び前記少なくとも１つの車両データに関連する測定ノイズの共分散のうちの少なくとも一方を含む、
請求項７に記載の方法。
【請求項９】
前記過去の平均値が、前記それぞれの車両モデルに対して仮定された初期平均値であり、又は以前に実施された荷重推定値を特定する段階から結果として得られたそれぞれの車両モデルの平均値であり、及び／又は、
前記過去の共分散が、前記それぞれの車両モデルに対して仮定された初期共分散であり、又は以前に実施された荷重推定値を特定する段階から結果として得られたそれぞれの車両モデルの共分散であり、及び／又は
前記過去の確率が、前記それぞれの車両モデルに対して仮定された初期確率であり、又は以前に実施された現在の確率を特定する段階から結果として得られたそれぞれの車両モデルの確率である、
請求項７又は請求項８に記載の方法。
【請求項１０】
前記車両の実際の質量が前記それぞれの重量区分に属することに関する現在の確率を特定する前記段階が、ベイズの定理を用いる段階を含む、
請求項１から請求項９のうちの一項に記載の方法。
【請求項１１】
前記車両の実際の質量が前記それぞれの重量区分に属することに関する現在の確率を特定する前記段階が、前記それぞれの荷重推定値の尤度及び以前に算出されたそれぞれの重量区分の確率に基づいて前記現在の確率を算出する段階を含む、
請求項１から請求項１０のうちの一項に記載の方法。
【請求項１２】
前記それぞれの荷重推定値の尤度が、荷重のない車両の前記実際の質量が既知である場合に、前記車両の少なくとも２つの異なる荷重で前記それぞれの車両モデルを励振することによる車両の加速度の推定に基づいて特定される、
請求項１１に記載の方法。
【請求項１３】
荷重推定値を特定する前記段階が、前記それぞれの車両モデルについて前記車両の現在の荷重と、前記車両の現在の荷重で励振されたそれぞれの車両モデルから得られる荷重推定値との間の残差の重量関連関数を特定する段階を含む、
請求項１から請求項１２のうちの一項に記載の方法。
【請求項１４】
荷重推定値を特定する前記段階が、各車両モデルについて前記それぞれの荷重推定値の不確実性を特定する段階を含む、
請求項１から請求項１３のうちの一項に記載の方法。
【請求項１５】
前記不確実性が、
前記車両のギアボックスの動作モードと、
加速度の大きさと、
エンジン・トルクの大きさと、
のうちの少なくとも１つによって特定される、
請求項１から請求項１４のうちの一項に記載の方法。
【請求項１６】
前記車両モデルが、前記車両の加速度についてのモデル、及び車両のセンサ手段によって供給される車両の現在の加速度を示す車両加速度値についてのモデルを含む、
請求項１から請求項１５のうちの一項に記載の方法。
【請求項１７】
前記車両の加速度についての前記モデルが、次式、すなわち、
【数１】
2011505566.tif 000041
ここで、ηtmがギアボックス効率、Ｔeが現在のエンジン・トルク、ζがギアボックス減衰定数、ωeが現在のエンジン角速度、ｉtが現在のギア比と現在の最終駆動ギア比との積、ρが現在の空気密度、Ｃdが空力抗力係数、Ａfが前記車両の空力前面面積、ｖが現在の車速、ｍヤマが過去の推定車両荷重、γmが前記車両の加速度に対する質量因子、ｇが重力定数、Ｃrが転がり抵抗係数、ｒが現在のホイール半径である、
に基づいている請求項１６に記載の方法。
【請求項１８】
前記車両の加速度値についての前記モデルが、次式、すなわち、
【数２】
2011505566.tif 000042
ここで、ｖドットが車両加速度、ｇが重力定数、αが前記センサ手段の測定軸と地球の重力場との間の角度、Ｃrが転がり抵抗係数、γmが前記車両の加速度に対する質量因子である、
に基づいている、請求項１６又は請求項１７に記載の方法。
【請求項１９】
前記センサ手段の出力についての前記モデルが、次式、すなわち、
【数３】
2011505566.tif 000043
ここで、ｙaccが前記車両の測定加速度、ｇが重力定数、Ｃrが転がり抵抗係数、γmが前記車両の加速度に対する質量因子である、
に基づいている、請求項１６又は請求項１７に記載の方法。
【請求項２０】
前記質量因子γmが、下式、すなわち、
【数４】
2011505566.tif 000044
ここで、ｍが前記車両の現在の質量、ｊdlが駆動系慣性、Ｊeがエンジン慣性、ｉtが現在のギア比と現在の最終駆動ギア比との積である、
によって表される、請求項１６から請求項１９のうちの一項に記載の方法。
【請求項２１】
車両の現在の荷重の推定のための装置であって、
各々が、前記車両の質量を推定するための車両モデルを実装し且つ少なくとも２つの重量区分のうちの１つに関連付けられた少なくとも２つのフィルタを含むフィルタ・バンクと、
各フィルタに対して、前記車両の現在の走行状況を示す少なくとも１つの車両データ及び前記それぞれの重量区分に固有の少なくとも１つの現在のフィルタ・パラメータを供給するための入力と、
を備え、
前記各フィルタが、前記車両の現在の荷重のフィルタ固有の推定として荷重推定値を特定するように適合されており、前記装置が更に、
各荷重推定値に対して、前記車両の現在の荷重が前記それぞれの重量区分に属すると仮定できるか否かの現在指標値を特定するように適合された特定手段と、
全ての荷重推定値の中で、前記現在の走行状況における前記車両の現在の荷重の推定として最良の現在指標値が特定される前記荷重推定値を選択するように適合されており、及び／又は前記重量区分における前記現在指標値を重み付けして、前記現在の走行状況における車両の現在の荷重についての荷重推定として用いられる全体荷重推定値を前記フィルタ・バンクにおいて得るように適合された手段と、
を備える装置。
【請求項２２】
前記各フィルタが、モデルベース・フィルタ、カルマン・フィルタ、拡張カルマン・フィルタ、無香カルマン・フィルタ、制約条件付き拡張カルマン・フィルタ、粒子フィルタ、シグマ点フィルタ、点質量フィルタ、格子ベース・フィルタのうちの少なくとも１つを含む、
請求項２１に記載の装置。
【請求項２３】
少なくとも１つの車両データを供給するための前記入力が、
車両センサによって供給される現在の車両加速度を示す加速度センサ値（ｙ）を受け取るための入力と、
前記車両のエンジンによって発生した現在のトルクを示すエンジン・トルク値（Ｔe）を受け取るための入力と、
前記車両の少なくとも１つのホイールの現在の角速度を示す少なくとも１つのホイール角速度値（ｗw）を受け取るための入力と、
前記車両のエンジンの現在の角速度を示すエンジン角速度値（ωe）を受け取るための入力と、
現在の最終駆動比を示す最終駆動比値（ｉfd）を受け取るための入力と、
のうちの少なくとも１つを備える、
請求項２１又は請求項２２に記載の装置。
【請求項２４】
処理システム上で実行されたときに、車両の現在の荷重を推定する方法を実施するためのプログラム・コードを含むコンピュータ・プログラム製品であって、前記方法が、
少なくとも２つの重量区分を定義する段階と、
前記少なくとも２つの重量区分の各々に対して、各々が前記車両の現在の質量を推定するための車両モデルを実装するフィルタを含むフィルタ・バンクを設ける段階と、
各フィルタに対して、前記車両の現在の走行状況を示す少なくとも１つの車両データ、及び前記それぞれの重量区分に固有の少なくとも１つの現在のフィルタ・パラメータを供給する段階と、
各フィルタを用いて、前記車両の現在の荷重のフィルタ固有の推定として荷重推定値を特定する段階と、
各荷重推定値に対して、前記車両の現在の荷重が前記それぞれの重量区分に属すると仮定できるか否かの現在指標値を特定する段階と、
全ての荷重推定値の中で、前記現在の走行状況における前記車両の現在の荷重の推定として最良の現在指標値を有する荷重推定値を選択する段階、及び／又は前記重量区分における前記現在指標値を重み付けして、前記現在の走行状況における前記車両の現在の荷重についての荷重推定として用いられる全体荷重推定値を前記フィルタ・バンクにおいて得る段階と、
を含む、
コンピュータ・プログラム製品。
【請求項２５】
処理システム上で実行されたときに、前記請求項２から請求項２０のうちの少なくとも一項に定義される段階を実施するためのプログラム・コードを含む、
請求項２４に記載のコンピュータ・プログラム製品。
【請求項２６】
コンピュータ読み取り可能記憶媒体上又は処理システムのストレージ内に記憶された、請求項２４又は請求項２５に記載のコンピュータ・プログラム製品。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、一般的に車両の荷重の推定に関し、より具体的にはフィルタ装置によって処理された車両関連情報に基づく車両の現在質量の推定に関する。
【背景技術】
【０００２】
車両の全重量（すなわち、車両の自重、並びに存在する場合には荷重の重量）を認識することは、多くの自動車用途において大きな関心事である。このパラメータを正確に推定することができると、多くの制御判定及び故障診断システムを改善することができる。例えば、間接的タイヤ空気圧監視システムは、低タイヤ空気圧の指標値としてホイールのロール半径を用いる。ロール半径は車両の質量に相関付けられるので、車両の荷重の認識は極めて重要である。
【０００３】
米国特許第５，９７３，２７３号は、車両上に追加のセンサとして垂直加速度計を装着し、車両の垂直方向変位の周波数挙動の変化を求めることを開示している。この場合、追加のセンサが必要となることは不利であり、また、既存の車両において修正することなくこの手法を用いることはできない。
【０００４】
米国特許公開第２００２／００３８１９３  Ａ１号は、エアサスペンション・システムの圧力に関する情報に基づいて大型トラックの荷重を推定することを開示している。この手法は、エアサスペンションを有する車両に限定されるので、自動車分野において一般的に適用することは可能ではない。
【０００５】
欧州特許第１  ８２９  ７１４  Ａ１号は、直接タイヤ空気圧センサからの情報及びホイールのロール半径と共にタイヤモデルを用いて車両の荷重を算出する。この手法は、タイヤ空気圧及び／又はホイール半径が特定される用途では有用ではない。
【０００６】
更なる手法は、米国特許公開第２００５／００１０３５６  Ａ１号に開示されるような、車両の縦方向動特性に基づいた車両質量の算出、及び米国特許第６，１６７，３５７号に開示された再帰最小二乗（ＲＬＳ）フィルタリングを用いた車両質量推定を含む。
【０００７】
車両の縦方向動特性に基づく質量算出は、例えば、推定において共に用いられる車両の加速度及び／又はエンジンの出力トルクが低い場合には、信頼性の高い質量推定が得られないことになる。再帰最小二乗（ＲＬＳ）フィルタリングを用いると、有用な結果を得るためにかなりの時間が必要とされ（例えば１０分程度）、走行状況の突然の変化に対する影響を受け易い。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】米国特許第５，９７３，２７３号公報
【特許文献３】米国特許出願第２００２／００３８１９３  Ａ１号公報
【特許文献４】欧州特許第１  ８２９  ７１４  Ａ１号公報
【特許文献５】米国特許出願第２００５／００１０３５６  Ａ１号公報
【特許文献６】米国特許第６，１６７，３５７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００９】
本発明の目的は、特定の車両のタイプに限定されることなく使用可能な、車両の現在の荷重を高い信頼性で迅速に推定する解決策を提供することである。
【課題を解決するための手段】
【００１０】
上記の目的を解決するために、本発明は、各々が独立請求項に定義された車両の現在の荷重を推定するための方法、装置、及びコンピュータ・プログラム製品を提供する。
【００１１】
本発明の１つの態様によれば、車両の現在の荷重を推定する方法が提供され、本方法は、
−少なくとも２つの重量区分を定義する段階と、
−少なくとも２つの重量区分の各々に対して、各々が車両の現在の質量を推定するための車両モデルを実装するフィルタを含むフィルタ・バンクを準備する段階と、
−各フィルタに対して、車両の現在の走行状況を示す少なくとも１つの車両データ、及びそれぞれの重量区分に固有の少なくとも１つの現在のフィルタ・パラメータを供給する段階と、
−各フィルタを用いて、車両の現在の荷重のフィルタ固有の推定として荷重推定値を特定する段階と、
−各荷重推定値に対して、車両の現在の荷重がそれぞれの重量区分に属すると仮定できるか否かの現在指標値を特定する段階と、
−全ての荷重推定値の中で、現在の走行状況における車両の現在の荷重の推定として最良の現在指標値を有する荷重推定値を選択する段階、及び／又は重量区分における現在指標値を重み付けして、現在の走行状況における車両の現在の荷重についての荷重推定として用いられる全体荷重推定値をフィルタ・バンクにおいて得る段階と、
を含む。
【００１２】
本発明の別の態様によれば、車両の現在の荷重を推定するための装置が提供され、本装置は、
−各々が、車両の質量を推定するための車両モデルを実装し且つ少なくとも２つの重量区分のうちの１つに関連付けられた少なくとも２つのフィルタを含むフィルタ・バンクと、
−各フィルタに対して、車両の現在の走行状況を示す少なくとも１つの車両データ、及びそれぞれの重量区分に固有の少なくとも１つの現在のフィルタ・パラメータを供給するための入力と、
を備え、
−各フィルタが、車両の現在の荷重のフィルタ固有の推定として荷重推定値を特定するように適合されており、本装置が更に、
−各荷重推定値に対して、車両の現在の荷重がそれぞれの重量区分に属すると仮定できるか否かの現在指標値を特定するように適合された特定手段と、
−全ての荷重推定値の中で、現在の走行状況における車両の現在の荷重の推定として最良の現在指標値が特定される荷重推定値を選択するように適合されており、及び／又は重量区分における現在指標値を重み付けして、現在の走行状況における車両の現在の荷重についての荷重推定として用いられる全体荷重推定値をフィルタ・バンクにおいて得るように適合された手段と、
を備える。
【００１３】
本発明の別の態様によれば、処理システム上で実行されたときに、車両の現在の荷重を推定する方法を実施するためのプログラム・コードを含むコンピュータ・プログラム製品が提供され、当該方法が、
−少なくとも２つの重量区分を定義する段階と、
−少なくとも２つの重量区分の各々に対して、各々が車両の現在の質量を推定するための車両モデルを実装するフィルタを含むフィルタ・バンクを準備する段階と、
−各フィルタに対して、車両の現在の駆動状況を示す少なくとも１つの車両データ、及びそれぞれの重量区分に固有の少なくとも１つの現在のフィルタ・パラメータを供給する段階と、
−各フィルタを用いて、車両の現在の荷重のフィルタ固有の推定として荷重推定値を特定する段階と、
−各荷重推定値に対して、車両の現在の荷重がそれぞれの重量区分に属すると仮定できるか否かの現在指標値を特定する段階と、
−全ての荷重推定値の中で、現在の走行状況における車両の現在の荷重の推定として最良の現在指標値を有する荷重推定値を選択する段階、及び／又は重量区分における現在指標値を重み付けして、現在の走行状況における車両の現在の荷重についての荷重推定として用いられる全体荷重推定値をフィルタ・バンクにおいて得る段階と、
を含む。
【００１４】
本発明の更なる態様、特徴、及び利点は、以下の説明、添付図面、及び附属の請求項から明らかとなるであろう。
【図面の簡単な説明】
【００１５】
【図１】本発明の実施可能な用途を説明するための変形したホイールの図である。
【図２】斜面上を走行している車両に作用する力を示す図である。
【図３】車両の駆動伝達系を示す図である。
【図４】本発明によるフィルタ・バンクを示す図である。
【図５】本発明による、異なる重量区分におけるモデルの残差分布を例示するグラフである。
【図６】本発明の好ましい実施形態における尤度対推定車両質量及び荷重をそれぞれ示す図である。
【図７】本発明の好ましい実施形態を用いた荷重／質量推定の結果を示す図である。
【図８】本発明による荷重／質量推定の結果をそれぞれの実際の車両質量と比較した図である。
【図９】本発明の好ましい実施形態のフロー線図である。
【発明を実施するための形態】
【００１６】
次に、本発明の実施形態を例証として添付図面を参照しながら説明する。
【００１７】
しかしながら、図面の説明に入る前に、本発明の更なる態様に対する幾つかの更なる考察を示す。本発明の方法に関する態様の詳細な考察は、明示しない場合でも、本発明の装置及びコンピュータ・プログラムに関する態様にも適用される。
【００１８】
本発明は、車両の現在の荷重の推定に関する。「空載の」車両の重量に基づいて発明の進歩性のある荷重を導出することができ、その逆もまた同様であるので、上記のことは、潜在的に車両の現在の重量（すなわち、「空載の」車両の重量、並びに存在する場合には荷重の重量）の推定にも関する。
【００１９】
少なくとも２つの重量区分が異なる荷重範囲を示すことができ、これらの範囲は重複する場合があり、又は重複しない場合もある。これらの範囲は、フィルタによる固有荷重推定においてそれぞれのフィルタ及び車両モデルがそれぞれ用いることができる。
【００２０】
これらの少なくとも２つの重量区分は、総車両重量又は全車両重量の異なる範囲を示すこともでき、これらの範囲は重複する場合があり、又は重複しない場合もある。これらの範囲は、フィルタ固有荷重推定において、それぞれのフィルタ及び車両モデルがそれぞれ用いることができる。
【００２１】
現在指標値を特定する段階は、確率、尤度、共分散、及び他の何れかの統計値、範囲、又は閾値のうちの少なくとも１つを表す情報及び／又はデータを特定する段階を含み、これに基づいて車両の現在の荷重がそれぞれの重量区分に属するか否かを仮定することができる。
【００２２】
最良の現在指標値を有する荷重推定値を特定する段階は、最も高い確率を有する荷重推定値、最も小さい共分散を有する荷重推定値、及びその荷重推定値が残りの荷重推定値と比較して現在の車両荷重のより良好な推定であることを示している別の統計値を有する荷重推定値のうちの少なくとも１つを特定する段階を含むことができる。
【００２３】
上記の供給段階、荷重推定値の特定段階、現在指標値の特定段階、及び選択段階は、繰り返すことができる。このことによって、例えば、異なる走行状況に対して固有の荷重推定を得るために、荷重推定を繰り返すことが可能になる。この繰り返しを用いて、個々のフィルタ推定及び／又はフィルタ・バンクの全体的な推定を改善するために、フィルタを再帰的に動作させることもできる。例えば、段階の繰り返しは、所定の時間間隔で、及び／又は車両の実際の質量を推定する度に実施することができる。
【００２４】
各フィルタは、モデルベース・フィルタ、カルマン（Ｋａｌｍａｎ）フィルタ、拡張カルマン・フィルタ、無香カルマン・フィルタ、制約条件付き拡張カルマン・フィルタ、粒子フィルタ、シグマ点フィルタ、点質量フィルタ、格子ベース・フィルタのうちの少なくとも１つを含むことができる。
【００２５】
車両モデルは、以下のパラメータのうちの少なくとも１つを含むことができる。
−ギアボックス効率（ηtm）
−エンジン・トルク（Ｔe）
−ギアボックス減衰定数（ζ）
−エンジン角速度（ωe）
−エンジン慣性（Ｊe）
−ホイール角速度（ωw）
−ホイール角加速度（ωwドット）
−ギア比（ｉg）
−最終駆動ギア比（ｉgf）
−ホイール慣性（Ｊe）
−空気密度（ρ）
−空力抗力係数（Ｃd）
−車両空力前面面積（Ａf）
−車速（ｖ）
−それぞれのフィルタから過去に特定された荷重推定値（ｍヤマ）
−それぞれのフィルタから車両の現在の荷重がそれぞれの重量区分に属すると仮定できるか否かについて過去に特定された指標値
−重力定数（ｇ）
−転がり抵抗係数（Ｃr）
−ホイール半径（ｒ）
【００２６】
【数１】
2011505566.tif 000003
【００２７】
現在の走行状況を示す少なくとも１つの車両データを供給する段階は、以下のうちの少なくとも１つを含むことができる。
−現在の車両加速度を示す加速度センサ値（ｙ）を供給する段階
−車両のエンジンによって発生した現在のトルクを示すエンジン・トルク値（Ｔe）を供給する段階
−車両の少なくとも１つのホイールの現在の角速度を示す少なくとも１つのホイール角速度値（ｗw）を供給する段階
−車両のエンジンの現在の角速度を示すエンジン角速度値（ωe）を供給する段階
−現在の最終駆動比を示す最終駆動比値（ｉfd）を供給する段階
【００２８】
それぞれの重量区分に固有の少なくとも１つの現在のフィルタ・パラメータを供給する段階は、以下のフィルタ・パラメータのうちの少なくとも１つを供給する段階を含むことができる。
−それぞれのフィルタの過去の荷重推定値
−それぞれのフィルタの過去の平均値
−それぞれのフィルタの過去の共分散
−それぞれのフィルタの過去の確率
【００２９】
該当する場合、共分散は、それぞれの車両モデルのモデル誤差の共分散、及び少なくとも１つの車両データに関連する測定ノイズの共分散のうちの少なくとも一方を含むことができる。
【００３０】
更に該当する場合には、過去の平均値は、それぞれの車両モデルに対して仮定された初期平均値とすることができ、又は以前に実施された荷重推定値を特定する段階から結果として得られたそれぞれの車両モデルの平均値であり、及び／又は過去の共分散は、それぞれの車両モデルに対して仮定された初期共分散であり、又は以前に実施された荷重推定値を特定する段階から結果として得られたそれぞれの車両モデルの共分散であり、及び／又は過去の確率は、それぞれの車両モデルに対して仮定された初期確率であり、又は以前に実施された現在の確率を特定する段階から結果として得られたそれぞれの車両モデルの確率である。
【００３１】
例えば、初期平均値の場合には、それぞれの重量区分の中央質量値を用いることができる。初期共分散の場合には、低いプロセスノイズを仮定することができる。
【００３２】
車両の実際の質量がそれぞれの重量区分に属することに関する現在の確率を特定する段階は、ベイズ（Ｂａｙｅｓ）の定理を用いる段階を含むことができる。
【００３３】
車両の実際の質量がそれぞれの重量区分に属することに関する現在の確率を特定する段階は、それぞれの荷重推定値の尤度と以前に算出されたそれぞれの重量区分の確率とに基づいて現在の確率を算出する段階を含むことができる。
【００３４】
それぞれの荷重推定値の尤度は、荷重のない車両の実際の質量が既知である場合に、車両の少なくとも２つの異なる荷重でそれぞれの車両モデルを励振することによる、車両の加速度の推定に基づいて特定することができる。
【００３５】
荷重推定値を特定する段階は、それぞれの車両モデルについて、車両の現在の荷重と、車両の現在の荷重で励振されるそれぞれの車両モデルから得られた荷重推定値との間の残差の重量関連関数を特定する段階を含むことができる。
【００３６】
荷重推定値を特定する段階は、各車両モデルについて、それぞれの荷重推定値の不確実性を特定する段階を含むことができる。
【００３７】
不確実性は、以下のうちの少なくとも１つによって特定することができる。
−車両のギアボックスの動作モード
−加速度の大きさ
−エンジン・トルクの大きさ
【００３８】
車両モデルは、車両の加速度についてのモデル、及び車両のセンサ手段によって供給される、車両の現在の加速度を示す車両加速度値についてのモデルを含むことができる。
【００３９】
車両の加速度のモデルは、次式に基づくことができる。
【００４０】
【数２】
2011505566.tif 000004
【００４１】
ここで、
ηtmは、ギアボックス効率、
Ｔeは、現在のエンジン・トルク、
ζは、ギアボックス減衰定数、
ωeは、現在のエンジン角速度、
ｉtは、現在のギア比と現在の最終駆動ギア比との積、
ρは、現在の空気密度、
Ｃdは、空力抗力係数、
Ａfは、車両の空力前面面積、
ｖは、現在の車速、
ｍヤマは、過去の推定車両荷重、
γmは、車両の加速度に対する質量因子、
ｇは、重力定数、
Ｃrは、転がり抵抗係数、
ｒは、現在のホイール半径である。
【００４２】
車両加速度値についてのモデルは、次式に基づくことができる。
【００４３】
【数３】
2011505566.tif 000005
【００４４】
ここで、
ｖドットは、車両加速度、
ｇは、重力定数、
αは、センサ手段の測定軸と地球の重力場との間の角度、
Ｃrは、転がり抵抗係数、
γmは、車両の加速度に対する質量因子である。
【００４５】
【数４】
2011505566.tif 000006
【００４６】
センサ手段の出力についてのモデルは、次式に基づくことができる。
【００４７】
【数５】
2011505566.tif 000007
【００４８】
ここで、
ｙaccは、車両の測定加速度、
ｇは、重力定数、
Ｃrは、転がり抵抗係数、
γmは、車両の加速度に対する質量因子である。
【００４９】
上記の質量因子γmは、次式で表すことができる。
【００５０】
【数６】
2011505566.tif 000008
【００５１】
ｍは、車両の現在の質量、
ｊdlは、駆動系慣性、
Ｊeは、エンジン慣性、
ｉtは、現在のギア比と現在の最終駆動ギア比との積である。
【００５２】
本発明の装置では、各フィルタは、モデルベース・フィルタ、カルマン・フィルタ、拡張カルマン・フィルタ、無香カルマン・フィルタ、制約条件付き拡張カルマン・フィルタ、粒子フィルタ、シグマ点フィルタ、点質量フィルタ、格子ベース・フィルタのうちの少なくとも１つを含むことができる。
【００５３】
少なくとも１つの車両データを供給するための入力は、以下のうちの少なくとも１つを含むことができる。
−車両センサによって供給される現在の車両加速度を示す加速度センサ値（ｙ）を受け取るための入力
−車両のエンジンによって発生した現在のトルクを示すエンジン・トルク値（Ｔe）を受け取るための入力
−車両の少なくとも１つのホイールの現在の角速度を示す少なくとも１つのホイール角速度値（ｗw）を受け取るための入力
−車両のエンジンの現在の角速度を示すエンジン角速度値（ωe）を受け取るための入力
−現在の最終駆動比を示す最終駆動比値（ｉfd）を受け取るための入力
【００５４】
本発明のコンピュータ・プログラム製品は、処理システム上で実行されたときに、上記で言及した方法段階を実施するためのプログラム・コードを含むことができる。
【００５５】
本発明のコンピュータ・プログラム製品は、コンピュータ読み取り可能記憶媒体上に記憶することができ、又は処理システムのストレージ内に記憶される。
【００５６】
車両の質量又は荷重が認識された場合、車両の多くの制御判定及び故障診断システムを改善することができる。このことは、例えば、間接的タイヤ空気圧監視を含む用途にも当てはまる。間接的タイヤ空気圧監視では、タイヤのロール半径の減少は、タイヤの空気が抜けたことの主な指標である。しかしながら、これは、車両の荷重によって（更に）生じる場合がある。従って、車両の質量又は荷重を認識又は推定することは、例えば、タイヤが不十分なタイヤ空気圧及び／又は車両荷重によって変形されているかどうかを特定するのに用いることができる。
【００５７】
図１に例示しているホイール状態は、車両の荷重状況を認識していなければ、タイヤ空気圧が低いという仮定を導く可能性がある。一方、車両が高荷重であることがわかった場合には、これはタイヤの変形が原因であるとみなされる可能性が高い。
【００５８】
モデルベース・パラメータ推定を用いて車両の質量を推定する場合には、質量が重要なパラメータである系のモデルが好ましい。質量は、多くの車両動特性に影響を及ぼす可能性がある。例えば、トランク内の大きな荷重は、重心が後方に移動すること及び側方スリップを小さくするタイヤの法線力が大きいことに起因して、車両のアンダーステア・グラディエントを変化させる。
【００５９】
本発明では、車両に作用する縦方向力の関係を用いることが好ましい。車両の加速中の縦方向動特性は、ニュートンの第２法則に基づいて表すことができる。更に縦方向動特性についての情報は、一般に、既に設置されている車両センサからのデータ及び／又は信号の形態で利用可能であり、従って、付加的なハードウェアを追加する必要はない。
【００６０】
車両質量及び該当する場合にはその荷重（すなわち、実際の車両質量からの重量偏差）の良好な推定を得るためには、走行状況と車両質量との間の関係を記述するモデルは、可能な限り正確でなければならず、及び／又は提供可能（例えば、計算処理の観点から）でなければならない。好ましい実施形態の説明では、縦方向動特性モデルを用いる。より具体的には、本明細書で説明する好ましい実施形態は、加速中の車両縦方向動特性を用いる。
【００６１】
当該目的に好適なモデルは、入力信号と測定値との間の物理的関係を記述する。これは、縦方向動特性モデルにおいては、一方では「入力信号」としてのエンジン・トルクと、他方では「測定値」としての車速及び加速度との間の関係である。
【００６２】
車両を剛体と仮定すると、ある一定の加速度を生じる力は、運動中の物体についてのニュートンの第２法則に従って表すことができる。
【００６３】
【数７】
2011505566.tif 000009
【００６４】
車両を加速する牽引力は、トランスミッション及びホイールを通じて道路にまで作用し、ホイールと道路との間の摩擦を用いて車両に前進力を与えるエンジン・トルクであり、並びに、車両が坂を下っている場合には重力の縦方向成分ｍ＊ｇ＊ｓｉｎ（α）であり、ここで、道路の傾斜角であるαは、加速力を生じる負であると考えられる。
【００６５】
車両を減速する力は、空気抵抗及び転がり抵抗を含み、更に、車両が坂を上っている場合には重力の縦方向成分ｍ＊ｇ＊ｓｉｎ（α）を含み、ここで道路の傾斜角であるαは、減速力を生じる正であると考えられる。
【００６６】
図２は、このような加速力及び減速力を例示しており、ここでＦは牽引力を示し、Ｆaは空気抵抗から生じる減速力（すなわち空気抵抗力）を示し、Ｆrは転がり抵抗を示し、Ｆgは車両に作用する重力から生じる力を示し、ａは結果として生じる車両加速度を示す。ニュートンの第２法則により次式が得られる。
【００６７】
【数８】
2011505566.tif 000010
【００６８】
図２では、車両に作用する加速度はまた、道路と比較した車両本体の加速度ｖドットと、地球重力の車両の縦方向に作用する成分ｇ＊ｓｉｎαとに分割される。一定速度での登坂走行時には、車両は、平坦な道路上での加速と同様に地球重力場の影響を受け、相違点は、回転部品、ホイール、駆動系及びエンジンに角加速度がないことである。
【００６９】
上記の力の各１つについての式は、図３に例示している車両の駆動伝達系を下記に説明するようにモデル化することによる牽引力Ｆに基づいて、更にまた下記に更に説明する外力に基づいて導出することができる。
【００７０】
エンジンのモデル化では、エンジンから伝達されるトルクが、トランスミッションを通って地面までの途中で損失を生じることを考慮に入れることができる。加速時には、車両を加速する目的のエンジン力は、車両の質量を加速するために通常は大部分が用いられることになる。しかしながら、エンジン力はまた、ホイールへの駆動伝達系の伝搬力の一部によって消費されることになり、例えば、力は、エンジン、プロペラシャフト、及びホイールのような回転部品の角加速度によって消費することができる。このような部品における慣性をＪe（エンジン慣性）、Ｊps（プロペラシャフト慣性）、及びＪw（ホイール慣性）で表す。平坦な道路上での加速時には、トルクは、エンジンの角加速度ωeドット（エンジン角加速度）及びホイールの角加速度ωｗドット（ホイール角加速度）によって失われる。回転加速におけるニュートンの第２法則によると、質量が角加速度ωドットにより与えられるような、慣性Ｊを有する質量（本体）を加速するためのトルクＴは次式である。
【００７１】
【数９】
2011505566.tif 000011
【００７２】
【数１０】
2011505566.tif 000012
【００７３】
従ってエンジンからギアボックスに伝達されるトルクは次式である。
【００７４】
【数１１】
2011505566.tif 000013
【００７５】
ここでＴgbはギアボックスに伝達されるトルクを表し、Ｔeは、エンジンによって発生するトルクを表し、Ｊeはエンジンの慣性を表し、ωeドットはエンジンの角加速度を表す。
【００７６】
プロペラシャフトに対して作用するトルクは、ギアボックスへの入力トルクＴgbにギア比を乗じたものである。しかしながら、ギアボックスでは、例えば、ギア歯間の摩擦、並びにベアリング及びシールにおける摩擦に起因した損失も存在する。この損失は、入力トルクＴgbに比例した一定の効率損失としてモデル化することができる。更に、ギアボックス及び最終駆動部における粘性損失も存在する可能性がある。一般にこれらの粘性損失は、ギアボックスの回転速度と共に変化し、従って、回転速度に比例する粘性減衰としてモデル化される。このことによって、ギアボックスからのトルク出力はギア比ｉgを用いてモデル化され、ギアボックス効率は、比例効率ηtmとして及び減衰係数ζを用いて粘性減衰としてモデル化される。これは次式で表すことができる。
【００７７】
【数１２】
2011505566.tif 000014
【００７８】
プロペラ及びドライブシャフトの慣性Ｊps及びホイールの慣性Ｊw、並びにホイールのトルク出力が最終駆動ギア比ｉgfに依存することを考慮すると、ホイール・トルクＴwは次式で表すことができる。
【００７９】
【数１３】
2011505566.tif 000015
【００８０】
ここで、ωpsドットはプロペラ及びドライブシャフトの角加速度を表し、ωeドットはホイールの角加速度を表す。
【００８１】
【数１４】
2011505566.tif 000016
【００８２】
最終駆動ギア比ｉgfが一定であると仮定すると、プロペラ及びドライブシャフトとホイールとの全体の慣性は次式で表すことができる。
【００８３】
【数１５】
2011505566.tif 000017
【００８４】
一方、最終駆動ギア比ｉgfとギア比ｉgとの積は、一般にｉtで表すことができる。
【００８５】
次に、結果として生じる牽引力Ｆは、結果として生じるホイール・トルクＴwをホイール半径ｒで除算することで表すことができる。
【００８６】
【数１６】
2011505566.tif 000018
【００８７】
牽引力Ｆの上記の式では、いかなる縦方向スリップも仮定していないが、通常の走行条件下ではスリップは通常２％を下回るので、これは妥当である点に留意されたい。
【００８８】
空力抵抗Ｆaは、速度ｖの二乗に比例し、次式で表すことができる。
【００８９】
【数１７】
2011505566.tif 000019
【００９０】
ここでＡfは車両の空力前面面積を表し、Ｃdは転がり抵抗係数を表し、ρは空気密度を表す。
【００９１】
転がり抵抗Ｃdは、様々な多少は複雑な手法でモデル化することができる。ここでは、転がり抵抗力が比例定数Ｃrで法線力に比例すると仮定している。道路の傾き角αで表すと、次式となる。
Ｆr＝Ｃrｍｇ
Ｆg＝ｍｇ・ｓｉｎ（α）
【００９２】
【数１８】
2011505566.tif 000020
【００９３】
【数１９】
2011505566.tif 000021
【００９４】
【数２０】
2011505566.tif 000022
【００９５】
質量因子は質量に依存し、これをパラメータとして用いて車両の質量を推定するときには質量依存性が失われる。
【００９６】
また項ｇ＊ｓｉｎαもγmで除算され、このことは、この補償に道路勾配情報が必要とされることを意味する。この手法は、γmの質量依存性及び道路勾配を無視できる平坦な道路に対して、加速度計においてどのような偏り誤差もないと仮定できる場合に良好に機能し、ωｗドット及びωeドットを得る上で、ホイール速度を微分する必要がない利点を有し、すなわち次式となる。
【００９７】
【数２１】
2011505566.tif 000023
【００９８】
通常、車両には、とりわけ縦方向加速度計を含むセンサ手段が装備される。この加速度計は、車両の実加速度／減速度を記録すると共に、センサが地球の重力場に対して垂直ではないとき、すなわち、車両が降坂又は登坂しており、センサが加速度の重力成分を記録するとき、或いは車両が加速に起因するピッチトルク又はトランクの荷重によって影響を受けてサスペンションが圧縮／収縮され、車両にピッチ角をもたらすときには、センサは、重力場の成分ｇ＊ｓｉｎαを記録する。第１の場合の角度をａslopeと呼び、他方の場合の角度をａpitchと呼ぶことによってこれらの異なる効果を区別する。更に加速度計が、ｂで表す偏り誤差と、センサノイズν（ｔ）とを有すると仮定する。縦方向加速度計のモデルは次式となる。
【００９９】
【数２２】
2011505566.tif 000024
【０１００】
この式は、角加速度を明示的に算出することなく、次式で再度定式化することができる。
【０１０１】
【数２３】
2011505566.tif 000025
【０１０２】
上記のモデル、少なくともその一部を用いて、本発明を実装することができる。以下では、各々がある特定の（真の）パラメータ・セット（例えば異なる質量仮定）に調整された異なるモデル／フィルタで複数のモデル又はフィルタ・バンク方式を用いる実施形態を参照する。より多くの観測値、例えば加速度計信号が得られるように、各フィルタにおける確率を再帰的に計算することができる。図４は、各々がフィルタ固有の車両モデルを実装する幾つかのフィルタを含むこのようなフィルタ・バンクを概略的に示している。各フィルタ及びその車両モデルはそれぞれ、再帰的に更新され、荷重推定及び／又はフィルタ・バンクによって全体として提供される全体推定を改善することになる。より具体的には、図４は、拡張カルマン・フィルタＥＫＦを有する例証的なフィルタ・バンクを示し、この各々は、測定値としてのセンサ加速度値ｙを現在の走行状況を示す車両データとして利用し、更に以前に特定されたフィルタ平均値ｘ、共分散Ｐ、及び確率γをフィルタ・パラメータとして用いる。これらのフィルタは、各々が独立して点推定を生成するように再帰的に用いることができ、又は互いに相互作用して組み合わせることができる。
【０１０３】
ある重量区分に属する真の質量の確率を算出するために、本明細書では例示的な実施形態としてベイズの定理を用いる。観測値が与えられると、ベイズの定理は、ある概念が真である確率を算出できるようにする。この算出は、この概念が真であり、且つ概念の確率の認識が事前に与えられた場合、観測値が真である確率を認識することによって行われる。
【０１０４】
【数２４】
2011505566.tif 000026
【０１０５】
質量推定の場合には、この確率は、真の質量がテスト対象の重量区分に属する確率として解釈することができる。この確率は、現在の重量区分に対する現在の観測値の尤度、及びこのサンプルの前に算出された重量区分の確率である事前の認識値から算出することができる。分母は、確率を合計で１にするためのスケール因子である。
【０１０６】
【数２５】
2011505566.tif 000027
【０１０７】
この場合、全ての重量区分についての確率が算出された後、この結果がその和を用いて正規化される。この目的のためには、観測値についての尤度が必要とされる。１つの方法は、荷重が存在しないと仮定されるモデルを用いて推定加速度を算出することとすることができる。この方法は、真の重量が既知である場合の様々な荷重を用いた異なるテストケースについて行うことができる。それぞれの残差の分布は、各重量区分における観測値についての確率関数と解釈することができる。
【０１０８】
このような種類の例示的な分布が図５に示されており、ここではテストケースを３つの重量区分に分割し、モデルをシミュレートする際に、モデル残差の分布を同じ質量（荷重なし）を用いて算出した。このことによって高い荷重では、分布はより左に歪んでいる。残差分布は、重量をゼロと仮定したときに、各重量区分において特定のモデル誤差が発生する確率を与える。
【０１０９】
上記の手法は、各１つが、発生する尤度が高い異なる系の状態をモデル化する、複数のモデルが存在しているかのように見える。この場合、真の系に対応していると考えられるものがこれらのうちの何れであるかを特定しなければならない。図５でわかるように、重量が幾つかの重量区分に離散化されている場合、モデル誤差の全体分布は、複数の幅狭なガウス分布の和で近似することができる。現在の重量を推定するために、測定値の残差を異なるガウス隆起に対して重み付けすることができる。この重み付けにより次式のガウス和フィルタが得られる。
【０１１０】
【数２６】
2011505566.tif 000028
【０１１１】
好ましい実施形態では、モデルベース・フィルタとして、拡張カルマン・フィルタ（ＥＫＦ）を用いることができる。非線形離散時間モデルが与えられると、この種のフィルタは、次式で表すことができる。
【０１１２】
【数２７】
2011505566.tif 000029
【０１１３】
Ｑ及びＲは、それぞれプロセスノイズ及び測定ノイズの共分散行列である。線形化は、ｆ及びｈのヤコビアンを算出することによって行われる。一般的にはこれらは次式で与えられる。
【０１１４】
【数２８】
2011505566.tif 000030
【０１１５】
【数２９】
2011505566.tif 000031
【０１１６】
【数３０】
2011505566.tif 000032
【０１１７】
【数３１】
2011505566.tif 000033
【０１１８】
更に、好ましい実施形態によれば、領域パラメータ化拡張カルマン・フィルタ（ＲＰＥＫＦ）と比較することができる手法を用いることができる。特に、１つのカルマン・フィルタが各重量区分に対して用いられ、初期重量はその重量区分の中央値に設定され、プロセスノイズを低く仮定される。各フィルタは、低速で突発的なモデル誤差の影響を受け難くなるが、アルゴリズムは、最も近いフィルタが最初から選ばれたときには極めて高速にすることができる。別の利点は、フィルタをある特定の基準に対して機能するように最適化された異なるものにすることができる点である。
【０１１９】
荷重がｎ個の重量区分に分割され、ｎ個の拡張カルマン・フィルタが稼働している場合には、各フィルタが最良のものである確率は、プロセスノイズがガウス分布であるという仮定に基づいて算出される。平均値μ及び標準偏差σを有するガウス分布関数は次式である。
【０１２０】
【数３２】
2011505566.tif 000034
【０１２１】
モデルが真の質量を用いてシミュレートされるときには、モデル誤差は、平均値ａ＝０のガウス分布を有すると仮定される。このことは、推定質量が真の質量に等しい場合に、残差が、平均値μ＝０のガウス分布を有さなければならないことを意味する。予測推定値の標準偏差は、共分散行列Ｐで記述される。残差の標準偏差を得るためには、Ｐを測定値空間にまで投影する必要があり、測定の不確実性を追加する必要があり、このことは、（２．６）の行列Ｓによって記述されたものそのものである。上記によって、フィルタ・バンク内の各フィルタについての残差の尤度関数が得られる。
【０１２２】
【数３３】
2011505566.tif 000035
【０１２３】
重量区分の尤度を算出するために、ベイズの定理によって次式が得られる。
【０１２４】
【数３４】
2011505566.tif 000036
【０１２５】
全フィルタの推定値及び共分散を得るために、次式に従って確率を共に重み付けすることができる。
【０１２６】
【数３５】
2011505566.tif 000037
【０１２７】
これはまた、異なる重量区分に関連付けられたフィルタを有するフィルタ・バンクの使用を示す図６から導くことができる。結果として得られた、各フィルタが真である尤度は、測定ノイズがガウス分布であるという仮定に基づいて特定される。次いで、尤度を共に重み付けし、推定車両荷重又は質量それぞれを特定し、全体的なフィルタ共分散を特定する。
【０１２８】
更なる実施形態では、最尤フィルタを推定値として用いることができ、これによってこのフィルタは、全てのフィルタが等確率である最初の時点からより速くなる。
【０１２９】
尤度と共に作動させる場合の別の利点は、モデルの不確実性を算入させるのが容易である点である。これは、サンプルがモデルに有利な状況の間に取得される場合には、モデルがモデル誤差を有することが多い状態で収集されたサンプルよりも可能性が高い推定をフィルタに保持させることができる。実施例としてギアボックスモデルにおける不確実性を考える。ギアボックスがロックアップモードに入ると、モデルは比較的良好になるが、市街地交通においては、ギアボックスは、滅多にロックアップモードには入らない。このことによって、ロックアップモードが発生するときにのみフィルタを稼働させることは、僅かなサンプルしか収集されず、統計の確実性を低くするので困難なことになる。その一方でフィルタは、ギアボックスがロックアップモードに入った場合のサンプルに対しては高い信頼性を置くことができる。またフィルタは、トルク入力が大幅に増大しているとき又はモデルの励振が低いときには、収集されたサンプルに対して極めて疑わしい可能性がある。
【０１３０】
この場合、式（２．８）においてモデル不確実性を追加することができる。測定ノイズ共分散Ｒ又はモデル誤差共分散Ｑは、モデル状態に依存するものとする。モデルが不確実である可能性の高い状態に対して共分散をより大きくした場合、Ｓは増大し、確率は小さくなる。この手法を用いると、多くのサンプルが収集され、推定値が得られることになるが、フィルタは、好ましい状態が発生して適正なフィルタに向かって高速に移行するまでは、比較的低速で不確実になる。モデルの状態は、３つのレベル、すなわち低いモデル精度、良好なモデル精度、及び最良のモデル精度に区分されている。この区分は、トランスミッションのロックアップ状態、加速度の大きさ、並びにトルクの大きさ及び微分値等の基準でモデルの状態確率を区分することによって行われる。この区分は、アルゴリズムのロバスト性を改善しており、モデルの不確実性対モデルの複雑性に関して問題の一部を解決し、多くのサンプルが破棄される場合にモデルの不確実性と統計の不確実性との間で妥協点を見いだすようにする。
【０１３１】
式（１．１２）を、式（２．４）と同様の離散状態空間モデル形式にすると次式が得られる。
【０１３２】
【数３６】
2011505566.tif 000038
【０１３３】
【数３７】
2011505566.tif 000039
【０１３４】
図７は、フィルタ・バンクアルゴリズムを用いた質量推定の結果を示している。アルゴリズムの始動時にどのようにして異なるフィルタが選ばれるかが分かる。これらの結果は、５個のフィルタに基づいて得られたものである。最初は、最大荷重を有するフィルタが最高の尤度を有するが、アルゴリズムは、より低い重量を有するフィルタをより尤度が高いものとして選び始め、フィルタは、約２０８０ｋｇの最低重量区分にまで段階的に降下する。次に、この場合での真の質量は最低重量区分よりも若干低かったので、フィルタは、真値に向かって緩慢に移行する。しかしながら、適正なフィルタが選ばれたときには、推定誤差は２５ｋｇを上回ることはない。長い時間の後、全てのフィルタは、緩慢にしか移行しないものの、真値に向かって移行することになる。
【０１３５】
ＥＫＦフィルタ・バンクを用いたテストケースによる結果は、アルゴリズムが、この章の冒頭に記載した全ての基準をどのように満たしているかを示している。アルゴリズムは、適正な値に向かって迅速に移行し、最良のモデル確率区分に対する全ての基準が満たされた場合には極めて迅速に移行する。更にこのアルゴリズムは、各ＥＫＦフィルタが極めて低速であり、大きいモデル誤差がある間は全てのフィルタの確率が低いので、突然のモデル誤差に対する影響を受け難い。アルゴリズムの性能は、何らかの形で、フィルタの収束時間の間に自動車がどのように走行されるかに基づいている。より好ましいモデル状態が始めから満たされた場合、フィルタは、適正なフィルタを即座に割り出すことが多い。アルゴリズムがどれほど迅速であるかについての尺度は、開始時点からフィルタが小さな変分しかなく、その後ある一定の値を有するまでの時間として設定される。
【０１３６】
【数３８】
2011505566.tif 000040
【０１３７】
変分ｖが５０ｋｇに設定された場合、アルゴリズムは、テストケースのほとんどにおいて、５分未満の時間内に適正な値から１００ｋｇ以内の質量を求め、すなわち、５分未満の時間内に適正な値を求めるほど十分に迅速である。このことはまた、例えば、市街地交通及び高速道路の速度を含む通常走行条件下での２５個のテストケースの結果を示す図７にも示されている。
【０１３８】
図９は、本発明の方法の例示的な実施形態、及び／又は装置の例示的な動作を示すフロー線図を示している。
【０１３９】
段階（１００）では、本方法及び／又は装置が初期化され、フィルタ確率は、各フィルタについての１次及び２次統計量と共に所与の事前統計量によって初期化することができる。
【０１４０】
段階（１０１）では、測定値及びフィルタ情報を用いて、ギア比、車速、その他等の物理信号を計算する。
【０１４１】
段階（１０２）では、荷重推定を更新すべきか否かが判定される。例えば、現在の走行状況及び／又は車両状態が信頼性の高い荷重推定を可能にする場合には、更新を実施することができる。例えば、車両の加速度、エンジン・トルク、速度、ヨーレート及び／又は横加速度、及び／又は制動操作及び／又は変速操作が、最適以下の又は信頼性が高くない荷重推定を生じると想定される場合には、いかなる更新も実施しないようにすることができる。
【０１４２】
段階（１０３）では、フィルタ（すなわちフィルタ・バンク）を用いて、車両からの情報（例えば、加速度センサ情報、ホイール速度／加速度情報、速度情報、エンジン・トルク情報、ギア（比）情報、その他）、及びフィルタ固有情報（例えば、確率及び／又は共分散）を利用してこれらのフィルタの荷重推定を行う。段階（１０３）では更に、例えば再帰的に動作させることによってフィルタ（すなわちフィルタ・バンク）を更新することができる。
【０１４３】
段階（１０４）では、工程フローをリセットすることができる。例えば、車両のドアが開けられた場合、及び／又は車両停止の場合、及び／又は荷重変化の可能性がある何れかの状況において、リセットを行うことができる。
【符号の説明】
【０１４４】
Ｆ  牽引力；  Ｆa  減速力（空気抵抗力）；  Ｆr  転がり抵抗；  
Ｆg  車両に作用する重力から生じる力；  ａ  車両加速度。
【図１】
2011505566.tif 000045
【図２】
2011505566.tif 000046
【図３】
2011505566.tif 000047
【図４】
2011505566.tif 000048
【図５】
2011505566.tif 000049
【図６】
2011505566.tif 000050
【図７】
2011505566.tif 000051
【図８】
2011505566.tif 000052
【図９】
2011505566.tif 000053
【国際調査報告】
2011505566.tif 000054
2011505566.tif 000055