(57)【特許請求の範囲】
【請求項１】
車両挙動（１０）に影響を及ぼす少なくとも１つの変数（β；ｖｘ，ｖｙ，Ψ，μ）の予測のための方法であって、
前記車両（１０）の動的変数（ＭＱ）を移動中に測定することと、
リアルタイムで前記変数（β；ｖｘ，ｖｙ，Ψ）の予測
0006573978.tif 000169
を、前記測定された動的変数（ＭＱ）に基づいて算出することと、
を含み、
前記少なくとも１つの変数（β；ｖｘ，ｖｙ，Ψ，μ）の前記予測を、予測手順（ＤＶＳβ；ＤＶＳβｖ；ＤＶＳβｖμ）により算出すること（２３０）であって、
前記車両（１０）の移動中に測定される動的変数（ＭＱ）のセットを、各々の時間間隔（ｎｙ，ｎｗ，ｎΨ，ｎｘ，ｎα）において考慮すること、および、
前記測定される動的変数（ＭＱ）のセットにおいて、前記変数（β；ｖｘ，ｖｙ，Ψ）に対して算出される少なくとも１つの最適な非線形回帰関数
0006573978.tif 000170
を適用して、前記変数（β；ｖｘ，ｖｙ，Ψ，μ）の前記予測を取得するために予測することであって、
前記最適な非線形回帰関数
0006573978.tif 000171
は、
基準データ（Ｄｄ）の取得されたセットおよび前記車両（１０）の移動中に測定される前記動的変数（ＭＱ）のセットに基づいて、所望の精度レベル（ε）に対して、動作条件（ＯＣ）の所与のセットにおける前記所望の精度レベル（ε）以下の予測エラーを与える回帰関数
0006573978.tif 000172
を見つけること、
を含む、最適な算出手順（２２０）により取得され、
前記基準データ（Ｄｄ）の取得されたセットは、前記動作条件（ＯＣ）の所与のセットにおいて、前記車両（１０）の前記測定された動的変数（ＭＱ）および前記車両（１０）の横速度（ｖｙ）および縦速度（ｖｘ）に対応する変数を含む変数の基準データ（Ｄｄ）のセットを取得すること（２１０）により取得される、
予測すること、
を備える、算出すること、
により特徴づけられる、方法。
【請求項２】
前記測定される動的変数（ＭＱ）は、ステアリング角（α）、横加速度（ａｙ）、４つの車輪速度（ｗｓ）、ヨーレート
0006573978.tif 000173
、縦加速度（ａｘ）を含み、前記動作条件（ＯＣ）の所与のセットにおいて、基準データ（Ｄｄ）のセットを取得する（２１０）前記動作は、前記測定される動的変数（ＭＱ）および横速度（ｖｙ）および縦速度（ｖｘ）に関するデータを取得することを含むことを特徴とする、請求項１に記載の方法。
【請求項３】
前記基準データ（Ｄｄ）のセットを、テスト車両（１０’）においてテストすることにより、および／または、前記車両（１０）のシミュレータ（１０’’）により、取得すること（２１０）により特徴づけられる、請求項１または請求項２に記載の方法。
【請求項４】
最適な非線形回帰関数
0006573978.tif 000174
を取得する前記動作（２２０）は、
所与のフェージングメモリを伴う関数のクラス
0006573978.tif 000175
に属する関数
0006573978.tif 000176
を考慮して、任意の時間間隔ｔに対する、および、動作条件（ＯＣ）の所与のセット全体に対する、前記変数（β；ｖｘ，ｖｙ，μ）に対する最大予測エラー
0006573978.tif 000177
を最小にする前記クラス
0006573978.tif 000178
における前記関数を見つけるステップと、
前記前のステップにおいて最適な非線形回帰
0006573978.tif 000179
関数として見つかった前記関数を割り当てて、前記車両挙動に影響を及ぼす前記変数の前記予測
0006573978.tif 000180
を取得する（２２０）ステップと、
を含むことを特徴とする、請求項１から３のうちのいずれか１項に記載の方法。
【請求項５】
予測される車両挙動（１０）に影響を及ぼす前記少なくとも１つの変数（（β；ｖｘ，ｖｙ，Ψ，μ）は、前記横滑り角
0006573978.tif 000181
であることを特徴とする、請求項１から４のうちのいずれか１項に記載の方法。
【請求項６】
車両挙動（１０）に影響を及ぼす前記少なくとも１つの変数（β；ｖｘ，ｖｙ，Ψ，μ）は、前記横速度（ｖｘ）および前記縦速度（ｖｙ）を含み、
縦速度の予測
0006573978.tif 000182
は、前記縦速度（ｖｘ）に関連する前記測定された動的変数（ＭＱ）のセットのサブセット（ｒｘ）に基づいて前記最適な算出により算出される最適な非線形回帰関数
0006573978.tif 000183
として取得され、
前記サブセット（ｒｘ）は、ステアリング角（α）、４つの車輪速度（ｗｓ）、および、縦加速度（ａｘ）と関連するデータを含み、
横速度の予測
0006573978.tif 000184
は、横加速度（ａｙ）およびヨーレート
0006573978.tif 000185
を含む前記測定された動的変数（ｒβ）のセットの追加的なサブセット（ｒｙ）および先に取得された縦速度の前記予測
0006573978.tif 000186
に基づいて前記最適な算出により算出される最適な非線形回帰関数
0006573978.tif 000187
として取得されることを特徴とする、請求項１から４のうちのいずれか１項に記載の方法。
【請求項７】
前記横滑り角
0006573978.tif 000188
の前記予測を、前記横速度の予測
0006573978.tif 000189
および前記縦速度の予測
0006573978.tif 000190
の関数として、特に、それらの比率の逆正接として算出することを含むことを特徴とする、請求項６に記載の方法。
【請求項８】
前記基準データセット（Ｄｄ）を複数（Ｌ）の基準サブセット（Ｄｄ１，Ｄｄ２，．．，ＤｄＬ）に、前記動作条件（ＯＣ）の所与のセットにおける対象の動作条件（（μ））の値（μ１，μ２）にしたがって、特に、前記道路タイヤ摩擦係数（μ）にしたがって分割することと、
前記基準サブセット（Ｄｄ１，Ｄｄ２，．．，ＤｄＬ）に対応して、さらに値がリアルタイムで前記車両（１０）において測定され
0006573978.tif 000191
、値が前記対象の条件（μ）に依存する、前記車両挙動に影響を及ぼす変数
0006573978.tif 000192
の複数の予測
0006573978.tif 000193
を、最適な非線形回帰関数
0006573978.tif 000194
を取得する（２２０）前記動作により取得することと、
前記測定された変数
0006573978.tif 000195
と、係数のベクトル
0006573978.tif 000196
により判定される前記係数を有する前記車両挙動に影響を及ぼす前記変数
の予測
0006573978.tif 000197
の線形的な組み合わせとの間の差を最小にする最適化問題の係数解決策のベクトル
0006573978.tif 000198
を取得することと、
を含むことを特徴とする、請求項１から５のうちのいずれか１項に記載の方法。
【請求項９】
前記対象の動作条件の前記実時間値（μ１，μ２）を、係数の前記ベクトル
0006573978.tif 000199
に基づいて検出することを含み、前記対象の動作条件は、特に、前記道路タイヤ摩擦係数であることを特徴とする、請求項８に記載の方法。
【請求項１０】
前記基準サブセット（Ｄｄ１，Ｄｄ２，．．，ＤｄＬ）に対応する前記横滑り角の予測
0006573978.tif 000200
を、最適な非線形回帰関数
0006573978.tif 000201
を取得する前記動作（２２０）により取得することと、
前記横滑り角（β）の前記予測
0006573978.tif 000202
を、前記線形パラメータ
0006573978.tif 000203
の前記横滑り角予測の線形的な組み合わせとして算出することと、
を含むことを特徴とする、請求項８または請求項９に記載の方法。
【請求項１１】
車両挙動（１０）に影響を及ぼす少なくとも１つの変数（β；ｖｘ，ｖｙ，Ψ，μ）の予測のための仮想センサであって、処理モジュール（１１）に実装され、前記車両（１０）の移動中に測定される前記動的変数（ＭＱ）のセットを利用して、前記測定された動的変数（ＭＱ）のセットに少なくとも１つの最適な非線形回帰関数
0006573978.tif 000204
を適用する前記車両（１０）の移動力学を示す少なくとも１つの変数（β；ｖｘ，ｖｙ，Ψ，μ）の前記予測を算出して（２３０）前記予測を取得するよう構成され、前記最適な非線形回帰関数
0006573978.tif 000205
は、請求項１から請求項１０のうちのいずれか１項に記載の方法にしたがって取得される、仮想センサ。
【請求項１２】
前記処理モジュール（１１）は、電子制御ユニットまたは前記車両（１０）の電子制御ボードに備えられる、請求項１１に記載のセンサ。
【請求項１３】
前記センサは、前記車両（１０）の移動中に測定される前記動的変数（ＭＱ）のセットを、ＥＳＣ（電子安定制御）システムを備える前記車両（１０）の動的変数（ＭＱ）を測定するよう構成されるモジュール（１２）から受信することを含む、請求項１１または請求項１２に記載のセンサ。
【請求項１４】
車両挙動（１０）に影響を及ぼす少なくとも１つの変数（β；ｖｘ，ｖｙ，Ψ，μ）の予測のためのシステムに統合され、前記車両（１０）の動的変数（ＭＱ）を移動中に測定するよう構成される前記モジュール（１２）を含み、
前記処理モジュール（１１）は、リアルタイムで、車両挙動（１０）に影響を及ぼす少なくとも１つの変数（β；ｖｘ，ｖｙ，Ψ，μ）の予測
0006573978.tif 000206
を、前記測定される動的変数（ＭＱ）に基づいて算出するよう構成されることを特徴とする、請求項１１から１３のうちのいずれか１項に記載のセンサ。
【請求項１５】
請求項１１から１４のうちのいずれか１項に記載のセンサを備える、特に自動車である、車両。
【請求項１６】
少なくとも１つのコンピュータのメモリ内にロードされることができ、製品が少なくとも１つのコンピュータ上で実行されると、請求項１から１０のうちのいずれか１項の方法のステップを実行可能なソフトウェアコードの部分を備える、コンピュータプログラム製品。
【発明の詳細な説明】
【技術分野】
【０００１】
本記載は、車両挙動に影響を及ぼす少なくとも１つの変数の仮想センサを策定する技術に関し、リアルタイムで対象の変数の信頼できる予測を、生産車および大型車両において典型的に利用可能な電子安定制御システムまたは他のデバイスから測定されるデータを使用して提供することが可能である。
【０００２】
これらの技術は、特に、車両の横滑り角および／または横および縦速度および／またはタイヤ−道路摩擦係数を予測することに向けられている。リアルタイムで生産車において、これらの変数の信頼できる予測を取得することにより、現在のシステムよりも大幅な性能向上を達成する安定性、トラクション、および、ブレーキング制御システムの策定が可能となり得る。
【背景技術】
【０００３】
多くの方法が、この数年で本発明の技術的な問題の焦点に関して、すなわち、車両挙動に影響を及ぼす変数（例えば、横滑り角）の予測に関して提案されており、その直接的な測定は、生産車での使用には複雑および／または高価すぎるセンサの使用を必要とする。これらの方法の全ては、２段階の手順に基づく。最初に、車両の適切なモデル（例えば、所与の自由度を伴う動的または運動学方程式を含むモデル）が導出される；その後、適切な予測アルゴリズム（とりわけ、カルマンフィルタ、または、スライディングモードオブザーバ、または、移動水平予測器、または、粒子フィルタに基づき得る）は、導出されたモデルに基づいて策定され、リアルタイムで対象の変数（例えば、横滑り角）を、ＥＳＣシステムまたは他の車両デバイスから入手可能なデータ（例えば、ステアリング角、車輪速度、ヨーレート、横加速度などに関する信号）を使用して予測するために電子ボードに実装される。
【０００４】
これらの２段階の方法は、重大な欠点に悩まされ、以下のように合成的に記載される。予測アルゴリズムは、実際の車両の動的挙動の近似的記載に過ぎない特定モデルで動作する。非常に正確なモデルが取得され得る場合でさえ、最適な予測（すなわち、最小変数）を見つけることは、特定モデルが非線形の場合は計算的に扱いにくく、計算的に扱い易いが必然的に近似の方法が使用される。モデル化および予測段階での近似値に起因して、どのくらい２段階の方法が正確であり得るかを評価するための方法は存在しない。予測エラーの有界性でさえ、複雑なシステムでは容易に達成されない。より関連性が強くても、車両モデルは、リアルタイムの値が通常の生産車で検出されない異なる動作条件（道路の乾湿、タイヤ摩耗状態、車両負荷など）にしたがって変化し得るパラメータに依存している。これらの問題（モデルおよび予測近似、可変動作条件）に起因して、いずれの上述された方法も、受け入れ可能な予測精度を達成する、横滑り角、縦および横速度、タイヤ−道路摩擦係数に対する予測アルゴリズムを策定する能力を達成したようには見えない。これは、さらに、これらの方法で策定されたこれらの変数の予測器が、商業車で利用可能となっていないことにより証明される。
【発明の概要】
【発明が解決しようとする課題】
【０００５】
１つ以上の実施形態の目的は、先行技術で達成可能な解決策に固有の制限を克服することである。
【０００６】
１つ以上の実施形態によると、その目的は、請求項１に指定される特徴を有する車両挙動に影響を及ぼす変数を予測する方法により、達成される。１つ以上の実施形態は、対応するシステム、仮想センサモジュール、そのようなシステムまたは仮想センサモジュールを備える車両に関するものであり得、さらに、少なくとも１つのコンピュータのメモリ内へロードされることができ、製品が少なくとも１つのコンピュータ上で実行されると、方法のステップを実行可能であるソフトウェアコードの部分を備える、コンピュータプログラム製品に関するものであり得る。本明細書において使用される場合、そのようなコンピュータプログラム製品という記載は、実施形態による方法の実装を統合するために処理システムを制御する命令を包含するコンピュータ可読手段という記載と同等であると、理解される。「少なくとも１つのコンピュータ」という記載は、モジュール形式および／または分散形式で実装される本実施形態の可能性を強調することが明白に意図される。少なくとも１つのコンピュータは、例えば、車両の電子制御ボード、または、上述の方法を実装するソフトウェアコードの部分を備える、いわゆる電子制御ユニット（ＥＣＵ）のレベルであり得る。
【０００７】
請求項は、様々な実施形態と関連して本明細書において提供される技術的な教示の不可欠な部分を形成する。
【課題を解決するための手段】
【０００８】
本明細書に記載の解決手段によれば、本方法は、少なくとも１つの変数の予測を予測手順により算出することを含む。
予測手順は、
各々の時間間隔における車両の移動中に測定される動的変数のセットを利用すること、および、
測定された動的変数のセットにおいて、変数に対して算出される少なくとも１つの最適な非線形回帰関数を適用して、変数の予測を取得するために予測することを備える。
最適な非線形回帰関数は、オフラインの最適な算出手順により取得される。
オフラインの最適な算出手順は、基準データの取得されたセットおよび車両の移動中に測定される動的変数のセットに基づいて、所望の精度レベルに対して、所望の精度レベル以下の予測エラーを動作条件の所与のセットにおいて提供する回帰関数を見つけることを含む。
基準データの取得されたセットは、動作条件の所与のセットにおいて、車両の測定された動的変数および車両の横および縦速度に対応する変数を含む変数の基準データのセットを取得することにより得られる。
【０００９】
様々な実施形態において、テスト車両および／または車両シミュレータから取得される測定変数は、ステアリング角、横加速度、４つの車輪速度、ヨーレート、縦加速度を含み、基準データのセットを取得する動作は、測定される変数および横速度および縦速度に関するデータを取得することを含み、後者は、テスト車両で動作する特定のセンサおよび／または車両シミュレータから取得される。
【００１０】
様々な実施形態において、基準データのセットは、テスト車両でテストすることにより、および／または、車両のシミュレータにより、取得される。
【００１１】
様々な実施形態において、非線形回帰関数を取得する動作は、所与のフェージングメモリを伴う関数のクラスに属する関数を考慮すると、横速度の縦速度に対する比率の逆正接を伴って、前回の時間間隔で測定された変数のデータのベクトルにおいて算出された、当該関数の異なるモジュールの測定の時間間隔における最大値を最小にするクラスにおいて関数を見つけることを含み、見出された関数を非線形回帰関数として割り当てて予測を取得する。
【００１２】
様々な実施形態において、横滑り角は、横速度および縦速度の予測の関数として、特にそれらの比率の逆正接として算出され、縦速度の予測は、縦速度に関連する測定された変数の当該セットのサブセットに基づいて算出される最適な非線形回帰関数として取得され、サブセットはステアリング角、４つの車輪速度および縦加速度に関するデータを含み、横速度の予測は、横加速度およびヨーレートを含む測定された変数のセットのさらなるサブセット、および、事前に取得された縦速度の予測に基づいて算出される最適な非線形回帰関数として取得される。
【００１３】
様々な実施形態において、車両挙動に影響を及ぼす少なくとも１つの変数の予測のためのセンサが、処理モジュールに実装され、車両の移動力学を表す少なくとも１つの変数の予測を、車両の移動中に測定される動的変数のセットを利用して、測定された動的変数のセットへ少なくとも１つの最適な非線形回帰関数を適用して算出し、当該予測を取得するよう構成され、最適な非線形回帰関数は、前述の実施形態のうちの任意の実施形態の方法にしたがって取得される。
【００１４】
様々な実施形態において、センサは、電子制御ユニットまたは車両の電子制御ボードに備えられる処理モジュールを有する。
【００１５】
様々な実施形態において、センサは、車両の移動中に測定される動的変数のセットを、ＥＳＣ（電子安定制御）システムを備える車両の動的変数を測定するよう構成されるモジュールから受信する。
【００１６】
様々な実施形態において、センサは、車両挙動に影響を及ぼす少なくとも１つの変数の予測のためのシステムに統合され、その移動中に車両の動的変数を測定するよう構成される当該モジュールと、リアルタイムで、車両挙動に影響を及ぼす少なくとも１つの変数の予測を、そのような測定された動的変数に基づいて算出するよう構成される処理モジュールと、を含む。
【００１７】
以下、実施形態を、付属の図を参照して限定されない単なる例として説明する。
【図面の簡単な説明】
【００１８】
【図１】図１は、リアルタイムで横滑り角を予測するために、オンラインで車両において１つの実施形態による方法を実施する仮想センサシステムのブロック図を示す。
【図２】図２は、仮想センサを策定するために、当該方法によりオフラインで行われる動作のブロック図を提示する。
【発明を実施するための形態】
【００１９】
次の記載は、実施形態の詳細な理解を目的とする様々な特定の詳細を示す。実施形態は、特定の詳細の１つ以上を伴わずに、または、他の方法、コンポーネント、材料などを伴って、実施されてよい。他の事例において、実施形態の様々な態様が不明瞭にされないように、既知の構造、材料、または、動作は、詳細に説明または記載されていない。
【００２０】
本記載のフレームワークにおける「実施形態」または「１つの実施形態」という記載は、実施形態と関連して記載される特定の構成、構造、または、特徴が、少なくとも１つの実施形態に包含されることを示すように意図されている。同様に、本記載の様々な部分に存在し得る、「実施形態において」または「１つの実施形態において」などの語句は、必ずしも１つの同じ実施形態を指していない。さらに、特定の形態、構造、または、特徴は、１つ以上の実施形態に適切に組み合されてよい。
【００２１】
本明細書において使用される記載は、単に利便性を目的としており、したがって、実施形態の保護の領分または範囲を規定しない。
【００２２】
図１において、車両の横滑り角のリアルタイム予測のための仮想センサの実施形態を概略的に提示するブロック図が示される。
【００２３】
数的な参照１０を伴って、測定される横滑り角βが関連付けられる車両を提示するブロックが示される。ＥＳＣシステム１２は、または、車両１０の安定性を車両１０の移動中における変数の値の測定に基づいて制御する車両の別のシステムは、本質的に既知の手法において、車両１０の以下の測定変数ＭＱ：ステアリング角α、横加速度ａｙ、４つの車輪速度ｗｓ、ヨーレート
0006573978.tif 000002
、縦加速度ａｘを移動中に測定する。上述されたように、これらは基本的にＥＳＣシステムにより測定される変数であり、本明細書に記載される方法は、リアルタイムで、好ましくは、ＥＳＣ１２により測定される変数にのみ基づいて動作する。そのような変数ＭＱは、横滑り角の仮想センサモジュール１１に提供され、この仮想センサモジュール１１は、そのような変数に基づいて、リアルタイムで、車両１０の横滑り角βの予測
0006573978.tif 000003
を供給する。
【００２４】
仮想センサ１１は、車両の移動に影響を及ぼす変数、特にこの事例においては、横滑り角の予測のための手順ＤＶＳβを実施し、図２に概略が記載される方法により策定される。方法は、テスト車両１０’からの測定により、または、車両シミュレータ１０’’を介するシミュレーションにより、または、データの一部を測定してデータの他の部分をシミュレートすることにより、基準データセットＤｄを取得するステップ２１０を含む。基準データの測定による取得に関連して、テスト車両１０’は、センサ１１が実装される車両１０の種類と同じでなければならない。
【００２５】
その後、基準データセットＤｄは、仮想センサＤＶＳβを策定する動作２２０へ渡され、仮想センサＤＶＳβは、その後、ステップ２３０において、ソフトウェアモジュール１１で実施される。このモジュールは、車両１０で利用可能なＥＣＵに埋め込まれ、リアルタイムで車両電子安定制御１２により測定される測定変数ＭＱに関する信号、すなわち、
0006573978.tif 000004
のみを受信し、リアルタイムで、横滑り角βの予測
0006573978.tif 000005
を供給する。
【００２６】
本明細書に記載される実施形態において、基準サンプルデータＤｄのセットは、以下である。
【数１】
0006573978.tif 000006
（１）
【００２７】
ｔが取得の時間間隔を示し、ΔｔからＮ＊Δｔへ変化し、ここで、Δｔはサンプリング時間であり、Ｎは取得サンプルの数である。以下において表記を簡素化するために、サンプリング時間Δｔ＝１と見なす。この基準データセットＤｄは、テスト車両で、特に自動車または大型車両で、典型的に測定される、考慮中の車両の動的性能を評価するための実験データを含んでよく、特に車両モデルを構築およびテストする既知の２段階の方法により使用される変数のデータを含んでよい。しかしながら、そのような基準データは、または、この基準データの一部は、さらに、直接的な測定により取得されるだけでなく、シミュレーションにより取得され得る：車両の信頼できるシミュレータが利用可能である場合、方程式（１）の基準データセットＤｄのデータは、そのようなシミュレータにより生成され得る。いずれの事例においても、基準データセットＤｄのそのようなデータは、実行またはシミュレートされた、速度、ブレーキング、道路−タイヤ摩擦状態、車の負荷、運転スタイル、および、他のパラメータに対する所与の範囲内で行われる運転テストと関連しなければならず、車両１０の異なる動作条件ＯＣを規定し、特に対象の動作条件を規定する。
【００２８】
方程式（１）の通り、基準データセットＤｄは、車両１０において車両の移動中に測定された変数ＭＱに対応し、ＥＳＣシステム１２により仮想センサ１１へ供給される取得基準データ
0006573978.tif 000007
すなわち、ステアリング角α、横加速度ａｙ、４つの車輪速度ｗｓ、ヨーレート
0006573978.tif 000008
、縦加速度ａｘを含む。しかしながら、基準データセットＤｄは、さらに、縦速度ｖｘおよび横速度ｖｙと関連する、さらに取得された基準データを含む。これらの縦速度ｖｘおよび横速度ｖｙデータは、好ましくは、車両シミュレータにより、仮想センサ１１のオフライン策定に必要な測定装置を最小に維持するために提供され得る。しかしながら、基準データセットＤｄがテスト車両における測定から取得される事例において、これらのデータは、ＧＰＳセンサを加えた実験用の光または慣性センサにより、本質的に既知の手法で提供され得る。縦速度ｖｘおよび横速度ｖｙデータは、仮想センサ１１により実施される予測手順ＤＶＳβの既定または策定においてのみ使用され、実時間予測には必要とされない。
【００２９】
本明細書において、一般的に、車両挙動に影響を及ぼす変数の予測に対する解決策が記載される。以下に記載される実施形態は、主に横滑り角の予測を取得することに向けられているが、これらの実施形態の一部において、本明細書に記載される方法は、車両の移動を表す際に対象となる他の変数を予測するために、必ずしも横滑り角をさらに取得する必要なく活用され得る。
【００３０】
本明細書に記載される解決策の実施形態は、一般的に、以下の形式の離散時間型の非線形回帰方程式として、仮想センサＤＶＳβの実現に基づく。
【数２】
0006573978.tif 000009
0006573978.tif 000010
すなわち、変数ＭＱが測定される時間ｔを考慮すると、次の時間間隔ｔ＋１で予測される横滑り角
0006573978.tif 000011
は、所与の整数である各々の時間間隔ｎｙ、ｎｗ、ｎΨ、ｎｘ、および、ｎαにおいて測定される、車両１０のＥＳＣシステム１２から取得される変数のベクトルである引数ベクトル
0006573978.tif 000012
で評価される非線形回帰関数
0006573978.tif 000013
の値により取得される。そのような時間間隔ｎｙ、ｎｗ、ｎΨ、ｎｘ、および、ｎαは、仮想センサＤＶＳβのメモリを規定し、他とは異なるものとして設定され得る。
【００３１】
仮想センサＤＶＳβを策定する動作２２０は、以下のε−ロバスト設計の問題の解決策である関数
0006573978.tif 000014
を見つけることにより取得される：基準データＤｄのセットに包含されるデータおよび測定精度の情報を利用すること、すなわち、各々の測定値は、その測定精度またはエラーを伴って考慮され、所望の精度レベルεに対して、
0006573978.tif 000015
として有界の予測エラーを与える回帰関数
0006573978.tif 000016
を見つけること、すなわち、横滑り角βとその予測
0006573978.tif 000017
との間の差のモジュールは、任意の時間ｔに対する、対象の動作条件ＯＣの範囲全体（動作条件ＯＣは、とりわけ、道路−タイヤ摩擦、車の負荷、タイヤの状態を含み得る）に対して所望の精度レベルε以下であり、テスト車両またはシミュレータは、そのような動作条件ＯＣを受けて、データセットＤｄをステップ２１０で取得する。
【００３２】
３つの実施形態の例が、本明細書において提示される。
【００３３】
（実施の形態１）
第１の実施形態において、ＤＶＳβとして示される仮想センサに至り、オフライン動作２２０は、以下の最適化問題の回帰関数
0006573978.tif 000018
解決策を見つけることと関連する。
【数３】
0006573978.tif 000019
（３）
ここで、基準データセットＤｄに包含されるデータは、（３）の右端側を計算するために使用される。
【００３４】
0006573978.tif 000020
は、以下のように規定される所与のフェージングメモリを伴う関数のクラスである。
【数４】
0006573978.tif 000021
すなわち、クラス
0006573978.tif 000022
の関数
0006573978.tif 000023
は、各時間ｔに対するリプシッツ条件を順守する。リプシッツ条件の定数は、ｋ番目の出力に対するパラメータγ、ρの生成物であり、ここで、０≦γ＜∞、０≦ρ＜∞であり、ｋは１からｍへ変化し、ここで、ｍ＝ｍａｘ［ｎｙ，ｎｗ，ｎΨ，ｎｘ，ｎα］であり、すなわち、変数ＭＱの測定の異なる間隔の中の最大間隔である。
【００３５】
規定により、横滑り角が、方程式（２）を参照して、
0006573978.tif 000024
および
0006573978.tif 000025
であることを考えると、（３）の回帰関数
0006573978.tif 000026
の解決策は、任意の時間間隔ｔに対して、および、基準データセットＤｄを取得するために使用される実験的な条件を構成する対象の動作条件ＯＣ（道路の乾湿または道路−タイヤ摩擦係数、車の負荷、タイヤの状態など）の範囲全体に対して、最大予測エラー
0006573978.tif 000027
を最小にするフェージングメモリ関数の当該クラス
0006573978.tif 000028
における関数
0006573978.tif 000029
として見つけられる。
【００３６】
方程式（３）および（４）におけるρおよびｍの値は、予測手順ＤＶＳβの設計パラメータである。すなわち、大きい値が選択されるほど、達成され得る予測エラーεは低くなるが、予測
0006573978.tif 000030
の過渡応答時間は大きくなる。
【００３７】
基準データＤｄにおいて動作するパラメータｙの値は、Ｃ．Ｎｏｖａｒａ、Ｆ．Ｒｕｉｚ、Ｍ．Ｍｉｌａｎｅｓｅによる論文「Ｄｉｒｅｃｔ  Ｆｉｌｔｅｒｉｎｇ：  Ａ  Ｎｅｗ  Ａｐｐｒｏａｃｈ  ｔｏ  Ｏｐｔｉｍａｌ  Ｆｉｌｔｅｒ  Ｄｅｓｉｇｎ  ｆｏｒ  Ｎｏｎｌｉｎｅａｒ  Ｓｙｓｔｅｍｓ（直接フィルタリング：非線形システムの最適なフィルタ設計に対する新しい手法）」（ＩＥＥＥ  Ｔｒａｎｓ．  ｏｎ  Ａｕｔｏｍａｔｉｃ  Ｃｏｎｔｒｏｌ，  ５８，  ｐｐ．８９−９９，  ２０１３）の項Ｄに記載される手順のステップ５ｂにしたがって選択され得る。
【００３８】
同じ論文の項ＩＩ．Ｄにおいて、最適化問題（３）を解決する方法が記載されている。
【００３９】
最終的に、ε＊が以下のように計算される。
【数５】
0006573978.tif 000031
ここで、基準データセットＤｄに包含されるデータは、（５）の右端側を計算するために使用される。
【００４０】
関数
0006573978.tif 000032
により記載される予測手順ＤＶＳβのエラーは、全ての時間ｔおよびデータセットＤｄに含まれる動作条件ＯＣの範囲全体に対して、
0006573978.tif 000033
として境界を示される。その後、
0006573978.tif 000034
の場合、導き出された仮想センサＤＶＳβは、ε−ロバスト設計問題に対する解決策である。
【００４１】
（３）から（５）の計算の全ては、オフライン動作２２０により、基準データセットＤｄに包含されるデータを使用して行われる。
【００４２】
この実施形態をリアルタイムで車両において実施する仮想センサは、図１の１１により提示され、目下の時間ｔに横滑り角β（ｔ）の予測β＊（ｔ）を以下のように計算する。
【数６】
0006573978.tif 000035
0006573978.tif 000036
ここで、要求時間における
0006573978.tif 000037
の値はＥＳＣモジュール１２から取得され、
0006573978.tif 000038
は、オフラインで（３）の解決策として計算される回帰関数である。
【００４３】
したがって、第１の実施形態は、実施形態において横滑り角βである、車両１０の挙動に影響を及ぼす少なくとも１つの変数の予測のための方法に関し、方法は、車両１０の変数ＭＱを移動中に測定すること、リアルタイムで横滑り角βにより提示される当該変数の予測
0006573978.tif 000039
を、そのような測定変数ＭＱに基づいて算出すること、車両１０の移動中に各々の時間間隔ｎｙ，ｎｗ，ｎΨ，ｎｘ，ｎαにおいて測定される変数ＭＱのセット、すなわち、ＥＳＣシステム１２から取得される変数のベクトル
0006573978.tif 000040
を利用すること、および、そのような測定変数のセットにおいて、そのような横滑り角βに対して算出された最適な非線形回帰関数
0006573978.tif 000041
を適用して、そのような予測を取得するために予測することを備える横滑り角βの予測を予測手順ＤＶＳβにより算出するステップ２３０を行うことを含み、当該の最適な非線形回帰関数
0006573978.tif 000042
は、以下を含む、オフラインの最適な算出手順により、すなわち、動作２２０により取得される：基準データＤｄの取得されたセットに基づいて、所望の精度レベルεに対して、当該の所望の予測レベルε以下の予測エラーを動作条件ＯＣの所与のセットにおいて与える回帰関数ｆβを見つけることであって、取得された基準データＤｄのセットは、動作条件ＯＣの当該所与のセットにおいて、車両１０の当該の測定変数ＭＱおよび車両１０の横速度ｖｙおよび縦速度ｖｘに対応する変数を含む変数の基準データＤｄのセットを取得する動作２１０により取得される、見つけること。
【００４４】
（実施の形態２）
車両の横滑り角の仮想センサに至り、ＤＶＳβｖとして示される第２の実施形態は、それぞれＤＶＳｖｘおよびＤＶＳｖｙと示される、縦速度ｖｘおよび横速度ｖｙの２つの仮想センサの策定に基づく。
【００４５】
これらの２つの仮想センサは、オフラインの動作２２０により策定され、以下のように基準データＤｄを利用する。
−仮想センサＤＶＳｖｘは以下の式で取得される。
【数７】
0006573978.tif 000043
0006573978.tif 000044
【００４６】
ｒｘは、測定変数ＭＱのサブセット、すなわち、４つの車輪速度ｗ、縦加速度ａｘ、および、ステアリング角αからなるベクトルである。様々な実施形態において、ベクトルｒｘを構成するサブセットに対する他の選択が可能である。
【００４７】
0006573978.tif 000045
は、最適化問題の解決策である。
【数８】
0006573978.tif 000046
ここで、基準データセットＤｄに包含されるデータは、（８）の右端側を計算するために使用される。
【００４８】
上記の最適化問題の
0006573978.tif 000047
解決策は、任意のｔに対する、および、基準データセットＤｄを取得するために使用される実験条件を構成する対象の動作条件ＯＣ（道路の乾湿、車の負荷、タイヤの状態など）の範囲全体に対する、最大予測エラー
0006573978.tif 000048
を最小化する縦速度ｖｘ（ｔ）の予測
0006573978.tif 000049
を与える。
−仮想センサＤＶＳｖｙは、以下の式で取得される。
【数９】
0006573978.tif 000050
0006573978.tif 000051
【００４９】
0006573978.tif 000052
は、測定変数ＭＱのサブセット、すなわち、横加速度ａｙ、ヨーレート
0006573978.tif 000053
からなる、（７）および（８）から取得される縦速度
0006573978.tif 000054
の予測を加えたベクトルである。様々な実施形態において、ベクトルｒｙを構成するサブセットに対する他の選択が可能である。
【００５０】
0006573978.tif 000055
は、最適化問題の解決策である。
【数１０】
0006573978.tif 000056
ここで、基準データセットＤｄに包含されるデータは、（１０）の右端側を計算するために使用される。
【００５１】
上記の最適化問題の
0006573978.tif 000057
解決策は、任意のｔに対する、および、基準データセットＤｄを取得するために使用される実験条件を構成する対象の動作条件ＯＣ（道路の乾湿、タイヤの状態など）の範囲全体に対する、最大予測エラー
0006573978.tif 000058
を最小化する縦速度
0006573978.tif 000059
の予測
0006573978.tif 000060
を与える。
−この第２の実施形態にしたがって策定される仮想センサＤＶＳβｖに提供される横滑り角βの予測
0006573978.tif 000061
は、以下のように規定される。
【数１１】
0006573978.tif 000062
すなわち、縦速度ｖｘの最適な予測
0006573978.tif 000063
に対する横速度ｖｙの最適な予測
0006573978.tif 000064
の比率の逆正接として、規定される。
−量
0006573978.tif 000065
は、以下のように計算される。
【数１２】
0006573978.tif 000066
ここで、基準データセットＤｄに包含されるデータは、（１１）の右端側を計算するために使用される。
【００５２】
この量は、この第２の実施形態にしたがって策定される仮想センサＤＶＳβｖの予測エラーにおいて境界を提供し、すなわち、全ての時間ｔに対する、および、データセットＤｄに含まれる動作条件ＯＣの範囲全体に対する
0006573978.tif 000067
を提供する。その後、
0006573978.tif 000068
の場合、導き出される予測手順ＤＶＳβｖは、ε−ロバスト設計問題に対する解決策である。最適化問題（８）および（１０）の
0006573978.tif 000069
および
0006573978.tif 000070
解決策は、先に引用されたＮｏｖａｒａ−Ｒｕｉｚ−Ｍｉｌａｎｅｓｅの論文に記載されているアルゴリズムを使用して取得され得る。
【００５３】
（７）から（１１）までの全ての計算は、オフラインで動作２２０により、上述されたように、基準データセットＤｄに包含されるデータを使用して行われる。
【００５４】
この実施形態をリアルタイムで車両において実施する仮想センサＤＶＳβｖは、図１において１１により提示され、目下の時間ｔに、以下の動作を行う横滑り角
0006573978.tif 000071
の予測
0006573978.tif 000072
を計算する。
−目下の時間ｔにおける縦速度
0006573978.tif 000073
の予測
0006573978.tif 000074
は、以下のように計算される。
【数１３】
0006573978.tif 000075
ここで、要求時間におけるｗｓ，ａｘ，αの値は、オンラインでＥＳＣモジュール１２から取得され、
0006573978.tif 000076
は、オフラインで（８）の解決策として計算される回帰関数である。
−目下の時間ｔにおける横速度
0006573978.tif 000077
の予測
0006573978.tif 000078
は、以下のように計算される。
【数１４】
0006573978.tif 000079
ここで、要求時間における
0006573978.tif 000080
の値は、オンラインでＥＳＣモジュール１２から取得され、要求時間における
0006573978.tif 000081
の値は、オンラインで（１２）から計算される縦速度の最適な予測であり、
0006573978.tif 000082
は、オフラインで（１０）の解決策として計算される回帰関数である。
−仮想センサＤＶＳβｖにより提供される目下の時間ｔにおける横滑り角
0006573978.tif 000083
の予測
0006573978.tif 000084
は、以下のように計算される。
【数１５】
0006573978.tif 000085
ここで、
0006573978.tif 000086
および
0006573978.tif 000087
は、（１２）および（１３）を介して事前に計算された最適な予測である。
【００５５】
この第２の実施形態は、第１の実施形態により与えられる仮想センサＤＶＳβより精度が高くあり得る横滑り角βの仮想センサＤＶＳβｖを与えることに加えて、さらに、仮想センサＤＶＳｖｘおよびＤＶＳｖｙを提供し、それぞれが車両の縦速度および横速度の予測を与える。これらの変数の信頼できる予測を取得することは、本質的に関連する技術の態様を提示し、トラクションおよびブレーキング制御システム、衝突回避システムなどを最適化するのに良好な値である
0006573978.tif 000088
および
0006573978.tif 000089
の知識であることに、留意されたい。
【００５６】
横速度
0006573978.tif 000090
および縦速度
0006573978.tif 000091
の予測を与える（１２）および（１３）の計算は、横滑り角
0006573978.tif 000092
を予測する方程式（１４）の算出の前に起こるので、本明細書に記載される第２の実施形態により、横速度および／または縦速度のみを予測する仮想センサＤＶＳｖｘおよびＤＶＳｖｙを、必ずしも横滑り角をさらに予測する必要なく策定できる。
【００５７】
（実施の形態３）
方法の第３の実施形態は、ＤＶＳβμとして示される横滑り角の仮想センサの策定に至り、ステップ２１０で、基準データセットＤｄをＬ個の基準サブセットＤｄ１，Ｄｄ２，．．，ＤｄＬに分割するように動作することを想定する。動作条件ＯＣの当該の所与のセットにおける対象の判定された動作条件により特定される当該サブセットの各々は、対象の動作条件の同じ値（または、値の範囲）に対して取得されるデータを包含する。本明細書に記載される実施形態において、分割は車両挙動に影響を及ぼす最も関連性の高い動作条件ＯＣであると見なされ得る道路−タイヤ摩擦係数μ１，μ２，．．，μＬにしたがって行われると見なされる。説明の簡易化を目的として、Ｌ＝２サブセットにおける分割が考慮されるが、方法は、以下に示されるように、より大きいＬの値へ容易に拡張され得る。基準データセットＤｄにおいてｔ＝１，．．，Ｍに対するデータが摩擦係数μ≒μ１を伴って取得され、ｔ＝Ｍ＋１，．．，Ｎに対するデータがμ≒μ２を伴って取得されると仮定すると、２つの基準サブセットは以下である。
【数１６】
0006573978.tif 000093
【００５８】
そのような分割ステップの後、オフラインの動作２２０において、以下の予測器が、ε−ロバスト設計問題の解決策である関数を見つけることにより評価される。すなわち、
−横滑り角の第１の予測器は、以下の式で評価される。
【数１７】
0006573978.tif 000094
ここで、ｒβ（ｔ）は方程式（２）において与えられ、
0006573978.tif 000095
は以下の最適化問題の解決策である。
【数１８】
0006573978.tif 000096
ここで、第１の基準データセット
0006573978.tif 000097
に包含されるデータは、（１６）の右端側を計算するために使用される。
−横滑り角の第２の予測器は、以下の式で評価される。
【数１９】
0006573978.tif 000098
ここで、
0006573978.tif 000099
は、以下の最適化問題の解決策である。
【数２０】
0006573978.tif 000100
ここで、第２の基準データセット
0006573978.tif 000101
に包含されるデータは、（１７）の右端側を計算するために使用される。
−動作２２０により、以下の式のヨーレートの第１の予測器が評価される。
【数２１】
0006573978.tif 000102
ここで、
0006573978.tif 000103
は、測定変数ＭＱのサブセット、すなわち、４つの車輪速度ｗｓ、横加速度ａｙ、および、ステアリング角αからなるベクトルであり、
0006573978.tif 000104
は、以下の最適化問題の解決策である。
【数２２】
0006573978.tif 000105
ここで、第１の基準データセット
0006573978.tif 000106
に包含されるデータは、（１８）の右端側を計算するために使用される。
−ヨーレートの第２の予測器は、以下の式により評価される。
【数２３】
0006573978.tif 000107
ここで、
0006573978.tif 000108
は、以下の最適化問題の解決策である。
【数２４】
0006573978.tif 000109
ここで、第２の基準データセット
0006573978.tif 000110
に包含されるデータは、（１９）の右端側を計算するために使用される。
【００５９】
（１６）から（１９）の全ての計算は、オフラインの動作２２０により、上述されたように、基準データセットＤｄ１およびＤｄ２に包含されるデータを使用して行われる。最適化問題（１６）（１７）（１８）（１９）の
0006573978.tif 000111
解決策は、先に引用されたＮｏｖａｒａ−Ｒｕｉｚ−Ｍｉｌａｎｅｓｅの論文に記載されるアルゴリズムを使用して取得され得る。
【００６０】
この実施形態をリアルタイムで車両において実施する仮想センサは、図１において１１により提示され、目下の時間ｔに、以下の動作を行う。
−目下の時間ｔにおけるヨーレート
0006573978.tif 000112
の２つの予測
0006573978.tif 000113
および
0006573978.tif 000114
は、以下のように計算される。
【数２５】
0006573978.tif 000115
0006573978.tif 000116
【数２６】
0006573978.tif 000117
0006573978.tif 000118
ここで、要求時間における
0006573978.tif 000119
の値はＥＳＣモジュール１２から取得され、
0006573978.tif 000120
および
0006573978.tif 000121
は、オフラインで（１８）および（１９）において計算される回帰関数であり、実際の動作条件が、それぞれμ＝μ１またはμ＝μ２であると見なす。
−以下の最適化問題の解決策
0006573978.tif 000122
が計算される。
【数２７】
0006573978.tif 000123
すなわち、リアルタイムでＥＳＣモジュール１２において測定される、目下の時間ｔにおけるヨーレート
0006573978.tif 000124
と、パラメータλの係数関数を有する第１のヨーレート予測
0006573978.tif 000125
および第２のヨーレート予測
0006573978.tif 000126
の線形的な組み合わせとの間の差の絶対値を考慮して、０から１の中に包含されるパラメータλの値を見つけ、そのような差を最小化する。この最適化問題は、オンラインで既知の線形検索方法を使用して効率的に解決され得る。
−目下の時間ｔにおける横滑り角
0006573978.tif 000127
の２つの予測
0006573978.tif 000128
および
0006573978.tif 000129
は、以下のように計算される。
【数２８】
0006573978.tif 000130
0006573978.tif 000131
【数２９】
0006573978.tif 000132
0006573978.tif 000133
ここで、要求時間における
0006573978.tif 000134
の値は、ＥＳＣモジュール１２２から取得され、
0006573978.tif 000135
および
0006573978.tif 000136
は、オフラインで（１６）および（１７）において計算される回帰関数であり、実際の動作条件は、それぞれμ＝μ１またはμ＝μ２であると見なす。
−最終的に、目下の時間ｔでの横滑り角
0006573978.tif 000137
の予測
0006573978.tif 000138
は、以下のように計算される。
【数３０】
0006573978.tif 000139
すなわち、最適なパラメータ
0006573978.tif 000140
の係数関数として有する横滑り角の第１および第２の予測の線形的な組み合わせとして、要するに、最適なパラメータ
0006573978.tif 000141
、および、それを補完する負の値として計算される。
【００６１】
したがって、先ほど記載された第３の実施形態は、横滑り角βの予測
0006573978.tif 000142
を、ヨーレートの予測を介して取得する。この実施形態の論拠は、以下の通りである。
【００６２】
方程式（２０）および（２１）は、リアルタイムでヨーレートの予測を与える２つの仮想センサを規定する。
【００６３】
しかしながら、ヨーレートは、リアルタイムでＥＳＣシステム１２により実際に測定されるため、予測される必要はない。ヨーレート予測は、２つのヨーレート予測をＥＳＣシステム１２により実際に測定される値と比較することにより最適化問題（２２）の解決策として取得される、最適なパラメータ
0006573978.tif 000143
を使用して、タイヤ−道路摩擦係数μの値をリアルタイムで検出するために実際に利用される。
【００６４】
車両１０が、例えば、道路上で実際の摩擦係数μ１を伴って動作している場合、第１のヨーレート予測
0006573978.tif 000144
は、（１８）を参照すると、同じ摩擦係数μ１の条件で測定されるデータから策定される
0006573978.tif 000145
であり、（１９）を参照すると、実際のものとは異なる摩擦係数μ２に対して策定される
0006573978.tif 000146
である第２のヨーレート予測
0006573978.tif 000147
より精度の高い予測を与える。したがって、最適化問題（２２）の解決策は、
0006573978.tif 000148
に至り、したがって、実際に車両１０は摩擦係数μ１で動作していることを検出する。その後、方程式（２５）から、横滑り予測
0006573978.tif 000149
が提供される。方程式（１６）および（２５）から確認され得るように、横滑り角の予測
0006573978.tif 000150
は、動作条件μ１に対する予測エラーを最小化する回帰関数
0006573978.tif 000151
により記載される仮想センサにより取得される。したがって、この予測は、基準データの分割をμの値にしたがって利用し、リアルタイムでタイヤ−道路摩擦係数の値を検出せずに、データセットＤｄが構成するμの値の範囲に対する予測性能のバランスを取る必要がある第１の実施形態の仮想センサＤＶＳβおよび第２の実施形態の仮想センサＤＶＳβｖにより達成可能な精度より良好な（または、最大限でも同等の）精度を達成する。
【００６５】
仮想センサＤＶＳβおよびＤＶＳβｖと比べて上述された予測精度の向上とは別に、この第３の仮想センサＤＶＳβμの追加的な興味深い特性は、この第３の実施形態が、ＤＶＳμと示される、リアルタイムで目下の時間ｔにおけるタイヤ−道路摩擦係数
0006573978.tif 000152
の予測
0006573978.tif 000153
を提供する仮想センサの取得を可能とすることであり、以下のように計算される。
【数３１】
0006573978.tif 000154
【００６６】
リアルタイムでのタイヤ−道路摩擦係数の予測は、本質的に関連する技術態様を提示する。この情報は、実際に、多くの車両挙動の問題に対する、例えば、トラクションおよびブレーキング制御、車両挙動制御、衝突回避などに対する関連値であるが、現段階で、通常の生産車において利用可能となっていない。タイヤ−道路摩擦係数のみを予測する仮想センサＤＶＳμの策定は、最適化問題（１８）および（１９）のみのオフライン解決策、および、動作（２０）（２１）（２２）および（２６）のみのオンラインの実行を必要とする。
【００６７】
上述されたように、Ｌ＝２サブセットへの分割に対して先ほど記載された動作は、より大きいＬ個の基準サブセットの値へ容易に拡張され得る。オフラインの動作２２０において、Ｌ個の関数
0006573978.tif 000155
およびＬ個の
0006573978.tif 000156
は、方程式（１６）（１７）および（１８）（１９）と類似のＬ個の方程式から計算され、各々が対応する基準サブセットＤｄ１，Ｄｄ２，．．，ＤｄＬに基づく。リアルタイムで車両において動作する仮想センサ１１は、目下の時間ｔにおいて、方程式（２０）（２１）と類似するＬ個の方程式により取得されるヨーレート
0006573978.tif 000157
のＬ個の予測
0006573978.tif 000158
および、方程式（２３）（２４）と類似するＬ個の方程式により取得される、横滑り角
0006573978.tif 000159
のＬ個の予測
0006573978.tif 000160
を計算する。ヨーレートの予測は、方程式（２２）のベクトルのバージョンにおいて使用され、ヨーレート
0006573978.tif 000161
と、そのような予測
0006573978.tif 000162
の線形的な組み合わせとの間の差を最小化する合計を有する正の係数のベクトル
0006573978.tif 000163
の値を判定し、係数として当該ベクトル
0006573978.tif 000164
のコンポーネントを有する。最終的に、横滑り角は、方程式（２５）のベクトルのバージョンを使用して、すなわち、横滑り予測の線形的な組み合わせ
0006573978.tif 000165
を使用して取得され、係数として当該ベクトル
0006573978.tif 000166
のコンポーネントを有する。
【００６８】
目下の実施形態は、値がリアルタイムで車両（１０）において測定され、予測される対象の条件に、例えば、道路−摩擦係数μに依存する変数の特定の例としてヨーレート（
0006573978.tif 000167
）に関して示されている。
【００６９】
しかしながら、第３の実施形態に関して記載される動作は、より一般的に、対象の動作条件ＯＣの実時間値を係数の当該ベクトル
0006573978.tif 000168
に基づいて検出するために使用され得ることが明白であり、以下のために使用され得ることが明白である。
基準データセットＤｄを複数のＬ個の基準サブセットＤｄ１，Ｄｄ２，．．，ＤｄＬに、動作条件ＯＣの所与のセットにおける対象の動作条件の値にしたがって分割すること、
当該の基準サブセット（Ｄｄ１，Ｄｄ２，．．，ＤｄＬ）に対応して、その値がリアルタイムで車両１０において測定され対象の当該条件ＯＣに依存する、車両挙動に影響を及ぼす変数の複数の予測を、最適な非線形回帰関数を取得する当該動作２２０により取得すること、
測定変数と、係数の当該ベクトルにより判定される係数を有する車両挙動に影響を及ぼす変数の該当予測の線形的な組み合わせとの間の差を最小化する最適化問題の係数解決策のベクトルΛ＊を取得すること。
【産業上の利用可能性】
【００７０】
本明細書に記載される様々な実施形態による解決策により、以下の利点を取得することが可能となる。
【００７１】
本明細書に記載される様々な実施形態による方法および仮想センサにより、従来技術の方法と異なり、通常の生産車に実装可能な対象の変数のリアルタイム予測が取得可能となる。特に、方法は、計算的に扱い易く、リアルタイムの値が生産車において検出されない生産車が動作しなければならない異なる動作条件（例えば、道路−タイヤ摩擦係数、負荷、タイヤの状態）において達成され得る予測精度を保証することが可能となる。
【００７２】
もちろん、実施形態の原理を害することなく、構築の詳細および実施形態は、次の請求項で規定される本実施形態の範囲から逸脱することなく、単なる例として本明細書におい記載および示されるものに対して幅広く変化してよい。
【００７３】
解決策はＥＳＣシステムにより測定される変数に関して記載されているが、解決策は、さらにリアルタイムで車両において利用可能な他の測定を使用して実施され得る。
【００７４】
本明細書において記載および特許請求される解決策は、車両の安定性、ステアリング、トラクションおよびブレーキングの監視および制御に関する自動車の安全システムの関連変数である車両の横滑り角、縦速度および横速度、タイヤ−道路摩擦係数の予測への特別な注意を伴って開発された。しかしながら、これらの特定の変数の記載は、いかなる手法においても、実施形態を限定するものとして理解されず、例えば、鉛直および回転変数など、サスペンション制御システムに関する他の車両挙動変数；例えば、航空機および海上船舶における姿勢予測など、自動車とは異なる分野に関する変数にも適用可能である。
【図１】
0006573978.tif 000207
【図２】
0006573978.tif 000208