(57)【特許請求の範囲】
【請求項１】
車両の各車輪の車輪速を検出する車輪速検出手段と、
前記車両の車速を検出する車速検出手段と、
前記車輪速検出手段によって検出した前記各車輪の車輪速に基づく左右輪逆相の車輪速、及び、前記車速検出手段によって検出した前記車速で所定の前後輪のホイールベースを割った値に基づいて、前記車両のばね下状態量及び／又はばね上状態量を推定する状態量推定手段と、を備え、
前記状態量推定手段は、
前記左右輪逆相の車輪速として、前輪位置の左右輪逆相の車輪速及び後輪位置の左右輪逆相の車輪速を表した下記［数１］を、前記車速で所定の前後輪のホイールベースを割った値を用いて下記［数２］に変換し、
【数１】
0006003931.tif 000028
【数２】
0006003931.tif 000029
上記［数２］から、前記ばね下状態量として、下記［数３］に示す前記車両の前輪ばね下左右輪逆相上下変位、及び、下記［数４］に示す前記車両の後輪ばね下左右輪逆相上下変位を推定し、前記ばね上状態量として、下記［数５］に示す前記車両のばね上ロール角を推定することを特徴とする車両状態推定装置。
【数３】
0006003931.tif 000030
【数４】
0006003931.tif 000031
【数５】
0006003931.tif 000032
（上記［数１］から［数５］において、ΔωＦ：前輪位置の左右輪逆相の車輪速、ΔωＲ：後輪位置の左右輪逆相の車輪速、Δω１：ＦＲ輪車輪速、Δω２：ＦＬ輪車輪速、Δω３：ＲＲ輪車輪速、Δω４：ＲＬ輪車輪速、αＸＦ，αＸＲ：サス上下ストロークに対するアクスル前後変位、αθＦ，αθＲ：サス上下ストロークに対するアクスルピッチ角、ｒＦ，ｒＲ：タイヤ半径（前輪、後輪）、ＷＦ，ＷＲ：アクスル位置からボデー重心の左右距離（トレッソの半分）、φＢＧ：ばね上ロール角、ＺＡ１〜ＺＡ４：右前、左前、右後、左後の各車輪のばね下の上下変位変動量、ＺＡＦ：前輪ばね下左右輪逆相上下変位、ＺＡＲ：後輪ばね下左右輪逆相上下変位、Ｌ：ホイールベース、Ｕ：車速、ｓ：ラプラス演算子）
【請求項２】
前記車両状態推定装置は、
前記状態量推定手段によって推定した前記前輪ばね下左右輪逆相上下変位、前記後輪ばね下左右輪逆相上下変位、及び、前記ばね上ロール角に基づいて、ロール共振周波数、及び、前記ばね下状態量に対する前記ばね上状態量の伝達関数最大値を算出し、算出した前記ロール共振周波数及び前記伝達関数最大値に基づいて、前記車両の積載状態を推定する積載状態推定手段、
を更に備える請求項１に記載の車両状態推定装置。
【請求項３】
車両の各車輪の車輪速を検出する車輪速検出手段と、
前記車両の車速を検出する車速検出手段と、
前記車輪速検出手段によって検出した前記各車輪の車輪速に基づく左右輪逆相の車輪速、及び、前記車速検出手段によって検出した前記車速で所定の前後輪のホイールベースを割った値に基づいて、前記車両のばね下状態量及び／又はばね上状態量を推定する状態量推定手段と、
前記状態量推定手段によって推定した前記ばね下状態量及び／又は前記ばね上状態量に基づいて前記車両のサスペンション装置を制御する制御手段と、を備え、
前記状態量推定手段は、
前記左右輪逆相の車輪速として、前輪位置の左右輪逆相の車輪速及び後輪位置の左右輪逆相の車輪速を表した下記［数６］を、前記車速で所定の前後輪のホイールベースを割った値を用いて下記［数７］に変換し、
【数６】
0006003931.tif 000033
【数７】
0006003931.tif 000034
上記［数７］から、前記ばね下状態量として、下記［数８］に示す前記車両の前輪ばね下左右輪逆相上下変位、及び、下記［数９］に示す前記車両の後輪ばね下左右輪逆相上下変位を推定し、前記ばね上状態量として、下記［数１０］に示す前記車両のばね上ロール角を推定することを特徴とする車両制御装置。
【数８】
0006003931.tif 000035
【数９】
0006003931.tif 000036
【数１０】
0006003931.tif 000037
（上記［数６］から［数１０］において、ΔωＦ：前輪位置の左右輪逆相の車輪速、ΔωＲ：後輪位置の左右輪逆相の車輪速、Δω１：ＦＲ輪車輪速、Δω２：ＦＬ輪車輪速、Δω３：ＲＲ輪車輪速、Δω４：ＲＬ輪車輪速、αＸＦ，αＸＲ：サス上下ストロークに対するアクスル前後変位、αθＦ，αθＲ：サス上下ストロークに対するアクスルピッチ角、ｒＦ，ｒＲ：タイヤ半径（前輪、後輪）、ＷＦ，ＷＲ：アクスル位置からボデー重心の左右距離（トレッソの半分）、φＢＧ：ばね上ロール角、ＺＡ１〜ＺＡ４：右前、左前、右後、左後の各車輪のばね下の上下変位変動量、ＺＡＦ：前輪ばね下左右輪逆相上下変位、ＺＡＲ：後輪ばね下左右輪逆相上下変位、Ｌ：ホイールベース、Ｕ：車速、ｓ：ラプラス演算子）
【請求項４】
車両の各車輪の車輪速を検出する車輪速検出手順と、
前記車両の車速を検出する車速検出手順と、
前記車輪速検出手順において検出した前記各車輪の車輪速に基づく左右輪逆相の車輪速、及び、前記車速検出手順において検出した前記車速で所定の前後輪のホイールベースを割った値に基づいて、前記車両のばね下状態量及び／又はばね上状態量を推定する状態量推定手順と、を含み、
前記状態量推定手順は、
前記左右輪逆相の車輪速として、前輪位置の左右輪逆相の車輪速及び後輪位置の左右輪逆相の車輪速を表した下記［数１１］を、前記車速で所定の前後輪のホイールベースを割った値を用いて下記［数１２］に変換し、
【数１１】
0006003931.tif 000038
【数１２】
0006003931.tif 000039
上記［数１２］から、前記ばね下状態量として、下記［数１３］に示す前記車両の前輪ばね下左右輪逆相上下変位、及び、下記［数１４］に示す前記車両の後輪ばね下左右輪逆相上下変位を推定し、前記ばね上状態量として、下記［数１５］に示す前記車両のばね上ロール角を推定することを特徴とする車両状態推定方法。
【数１３】
0006003931.tif 000040
【数１４】
0006003931.tif 000041
【数１５】
0006003931.tif 000042
（上記［数１１］から［数１５］において、ΔωＦ：前輪位置の左右輪逆相の車輪速、ΔωＲ：後輪位置の左右輪逆相の車輪速、Δω１：ＦＲ輪車輪速、Δω２：ＦＬ輪車輪速、Δω３：ＲＲ輪車輪速、Δω４：ＲＬ輪車輪速、αＸＦ，αＸＲ：サス上下ストロークに対するアクスル前後変位、αθＦ，αθＲ：サス上下ストロークに対するアクスルピッチ角、ｒＦ，ｒＲ：タイヤ半径（前輪、後輪）、ＷＦ，ＷＲ：アクスル位置からボデー重心の左右距離（トレッソの半分）、φＢＧ：ばね上ロール角、ＺＡ１〜ＺＡ４：右前、左前、右後、左後の各車輪のばね下の上下変位変動量、ＺＡＦ：前輪ばね下左右輪逆相上下変位、ＺＡＲ：後輪ばね下左右輪逆相上下変位、Ｌ：ホイールベース、Ｕ：車速、ｓ：ラプラス演算子）
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、車両状態推定装置、車両制御装置及び車両状態推定方法に関する。
【背景技術】
【０００２】
従来、車輪速から車両の状態量を推定する技術が報告されている。
【０００３】
例えば、特許文献１には、車両の前後軸回りの回転運動を検出するロール検出装置であって、車両左右の車輪速度を各々検出する第１，第２の車輪速度検出手段と、該検出された左右の車輪速度に基づいて、左右輪各々についてのばね上共振周波数領域の車輪速度の変動量を求める第１，第２の車輪速変動量抽出手段と、該求められた左右輪についての変動量の逆相成分に基づいて、車両前後軸周りの回転運動の大きさを演算するロール演算手段とを備えたロール検出装置が開示されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開平５−３１９０５１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、従来技術では、各入力（路面、操舵、アクセル、ブレーキ等）に対するばね上・ばね下挙動の伝達特性が必要であり、それらの伝達特性は車両状態（積載状態等）により変化する。そのため、従来技術においては、車輪速から車両の状態量を推定する場合、車両状態の変化（積載状態の変化等）によって車両特性（ばね下挙動に対するばね上挙動の伝達特性等）が変化するため、高精度での推定が困難であった。このように、従来技術は車輪速から車両の状態量を推定する場合の精度を向上することについて、なお改良の余地がある。
【０００６】
本発明の目的は、車両特性の変化によらず、車輪速から車両の状態量を高精度に推定することができる車両状態推定装置、車両制御装置及び車両状態推定方法を提供することである。
【課題を解決するための手段】
【０００７】
本発明の車両状態推定装置は、車両の各車輪の車輪速を検出する車輪速検出手段と、前記車両の車速を検出する車速検出手段と、前記車輪速検出手段によって検出した前記各車輪の車輪速に基づく左右輪逆相の車輪速、及び、前記車速検出手段によって検出した前記車速で所定の前後輪のホイールベースを割った値に基づいて、前記車両のばね下状態量及び／又はばね上状態量を推定する状態量推定手段と、を備えることを特徴とする。
【０００８】
上記車両状態推定装置において、前記状態量推定手段は、前記左右輪逆相の車輪速として、前輪位置の左右輪逆相の車輪速及び後輪位置の左右輪逆相の車輪速を表した下記［数１］を、前記車速で所定の前後輪のホイールベースを割った値を用いて下記［数２］に変換し、
【数１】
0006003931.tif 000002
【数２】
0006003931.tif 000003
上記［数２］から、前記ばね下状態量として、下記［数３］に示す前記車両の前輪ばね下左右輪逆相上下変位、及び、下記［数４］に示す前記車両の後輪ばね下左右輪逆相上下変位を推定し、前記ばね上状態量として、下記［数５］に示す前記車両のばね上ロール角を推定することが好ましい。
【数３】
0006003931.tif 000004
【数４】
0006003931.tif 000005
【数５】
0006003931.tif 000006
（上記［数１］から［数５］において、ΔωＦ：前輪位置の左右輪逆相の車輪速、ΔωＲ：後輪位置の左右輪逆相の車輪速、Δω１：ＦＲ輪車輪速、Δω２：ＦＬ輪車輪速、Δω３：ＲＲ輪車輪速、Δω４：ＲＬ輪車輪速、αＸＦ，αＸＲ：サス上下ストロークに対するアクスル前後変位、αθＦ，αθＲ：サス上下ストロークに対するアクスルピッチ角、ｒＦ，ｒＲ：タイヤ半径（前輪、後輪）、ＷＦ，ＷＲ：アクスル位置からボデー重心の左右距離（トレッソの半分）、φＢＧ：ばね上ロール角、ＺＡ１〜ＺＡ４：右前、左前、右後、左後の各車輪のばね下の上下変位変動量、ＺＡＦ：前輪ばね下左右輪逆相上下変位、ＺＡＲ：後輪ばね下左右輪逆相上下変位、Ｌ：ホイールベース、Ｕ：車速、ｓ：ラプラス演算子）
【０００９】
上記車両状態推定装置において、前記車両状態推定装置は、前記状態量推定手段によって推定した前記前輪ばね下左右輪逆相上下変位、前記後輪ばね下左右輪逆相上下変位、及び、前記ばね上ロール角に基づいて、ロール共振周波数、及び、前記ばね下状態量に対する前記ばね上状態量の伝達関数最大値を算出し、算出した前記ロール共振周波数及び前記伝達関数最大値に基づいて、前記車両の積載状態を推定する積載状態推定手段、を更に備えることが好ましい。
【００１０】
本発明の車両制御装置は、車両の各車輪の車輪速を検出する車輪速検出手段と、前記車両の車速を検出する車速検出手段と、前記車輪速検出手段によって検出した前記各車輪の車輪速に基づく左右輪逆相の車輪速、及び、前記車速検出手段によって検出した前記車速で所定の前後輪のホイールベースを割った値に基づいて、前記車両のばね下状態量及び／又はばね上状態量を推定する状態量推定手段と、及び前記状態量推定手段によって推定した前記ばね下状態量及び／又は及び前記ばね上状態量に基づいて前記車両のサスペンション装置を制御する制御手段と、を備えることを特徴する。
【００１１】
本発明の車両状態推定方法は、車両の各車輪の車輪速を検出する車輪速検出手順と、前記車両の車速を検出する車速検出手順と、前記車輪速検出手順において検出した前記各車輪の車輪速に基づく左右輪逆相の車輪速、及び、前記車速検出手順において検出した前記車速で所定の前後輪のホイールベースを割った値に基づいて、前記車両のばね下状態量及び／又はばね上状態量を推定する状態量推定手順と及び、を含むことを特徴とする。
【発明の効果】
【００１２】
本発明に係る車両状態推定装置は、車両の各車輪の車輪速を検出する車輪速検出手段と、前記車両の車速を検出する車速検出手段と、前記車輪速検出手段によって検出した前記各車輪の車輪速に基づく左右輪逆相の車輪速、及び、前記車速検出手段によって検出した前記車速で所定の前後輪のホイールベースを割った値に基づいて、前記車両のばね下状態量及び／又はばね上状態量を推定する状態量推定手段と及び、を備える。本発明に係る車両状態推定装置によれば、車両特性の変化によらず、車輪速から車両の状態量を高精度に推定することができるという効果を奏する。
【図面の簡単な説明】
【００１３】
【図１】図１は、実施形態に係る車両の概略構成図である。
【図２】図２は、実施形態に係る車両状態推定装置のブロック図である。
【図３】図３は、車輪の回転角度の変動量を説明する図である。
【図４】図４は、前後変位変動量の説明図である。
【図５】図５は、ピッチングによる前後変位変動量の説明図である。
【図６】図６は、車両の回転による前後変位変動量の説明図である。
【図７】図７は、サスペンション装置のストロークによる前後変位変動量の説明図である。
【図８】図８は、ピッチングによるばね下ピッチ角の説明図である。
【図９】図９は、サスペンション装置のストロークによるばね下ピッチ角の説明図である。
【図１０】図１０は、図２の車両状態推定装置のブロック図と本発明の実施形態に係る理論式とを対応付けて表した図である。
【図１１】図１１は、積載量、重心高さとロール共振周波数との関係を示す図である。
【図１２】図１２は、積載量、重心高さと伝達関数最大値との関係を示す図である。
【発明を実施するための形態】
【００１４】
以下に、本発明の実施形態に係る車両状態推定装置、車両制御装置及び車両状態推定方法につき図面を参照しつつ詳細に説明する。なお、この実施形態によりこの発明が限定されるものではない。また、下記の実施形態における構成要素には、当業者が容易に想定できるものあるいは実質的に同一のものが含まれる。
【００１５】
［実施形態］
図１から図１２を参照して、実施形態について説明する。本実施形態は、車両状態推定装置、車両制御装置及び車両状態推定方法に関する。図１は、実施形態に係る車両の概略構成図、図２は、本発明の実施形態に係る車両状態推定装置のブロック図、図３は、車輪の回転角度の変動量を説明する図、図４は、前後変位変動量の説明図、図５は、ピッチングによる前後変位変動量の説明図、図６は、車両の回転による前後変位変動量の説明図、図７は、サスペンション装置のストロークによる前後変位変動量の説明図、図８は、ピッチングによるばね下ピッチ角の説明図、図９は、サスペンション装置のストロークによるばね下ピッチ角の説明図、図１０は、図２の車両状態推定装置のブロック図と本発明の実施形態に係る理論式とを対応付けて表した図である。図１１は、積載量、重心高さとロール共振周波数との関係を示す図である。図１２は、積載量、重心高さと伝達関数最大値との関係を示す図である。
【００１６】
本実施形態に係る車両状態推定装置１０１は、左右輪逆相の車輪速とばね上、ばね下挙動の関係が車両特性（ばね下質量等）に影響されない特徴を活かして、左右輪逆相の車輪速から左右輪逆相のばね下・ばね上挙動を推定する。また、本実施形態に係る車両状態推定装置１０１は、推定した左右輪逆相のばね下・ばね上挙動からばね下に対するばね上の伝達関数を算出し、その伝達関数から車両の積載状態等を推定し、サスペンション制御やブレーキ制御のモデル設定に反映させる。以下にその詳細について説明する。
【００１７】
図１を参照して、本実施形態に係る車両１００について説明する。図１に示すように、車両１００は、車両状態推定装置１０１、車輪２（２ＦＲ，２ＦＬ，２ＲＲ，２ＲＬ）、サスペンション装置１０（１０ＦＲ，１０ＦＬ，１０ＲＲ，１０ＲＬ）を含んで構成されている。本実施形態の車両状態推定装置１０１は、ＥＣＵ１（後述する状態量推定手段や積載状態推定手段として機能を有する）、及び、車輪速検出手段４（４ＦＲ，４ＦＬ，４ＲＲ，４ＲＬ）を少なくとも含んで構成されている。また、後述するように、本実施形態のＥＣＵ１は、サスペンション装置１０（１０ＦＲ，１０ＦＬ，１０ＲＲ，１０ＲＬ）を制御する制御手段としての機能も有している。従って、車両１００は、ＥＣＵ１及び車輪速検出手段４を含む車両制御装置１０２を搭載している。
【００１８】
車両１００は、車輪２として、右前輪２ＦＲ、左前輪２ＦＬ、右後輪２ＲＲ及び左後輪２ＲＬを有している。また、車両１００は、サスペンション装置１０として、右前サスペンション装置１０ＦＲ、左前サスペンション装置１０ＦＬ、右後サスペンション装置１０ＲＲ及び左後サスペンション装置１０ＲＬを有している。なお、本明細書では、各構成要素の符号の添字ＦＲは右前輪２ＦＲに係るものを示す。同様にして、符号の添字ＦＬは左前輪２ＦＬ、添字ＲＲは右後輪２ＲＲ、添字ＲＬは左後輪２ＲＬに係るものであることを示す。
【００１９】
サスペンション装置１０は、ばね下とばね上とを接続している。ここで、ばね下は、フロントやリヤのアーム部材や、ナックル等を含むものであり、車両１００のうち、サスペンション装置１０に対して車輪２側に接続されている部分である。ばね上は、車両１００のうち、サスペンション装置１０によって支持される部分であり、例えばボデー３である。サスペンション装置１０は、伸縮することでばね上とばね下との相対変位を許容する。サスペンション装置１０は、例えば、鉛直軸に対して所定の角度で傾斜して設けられている。
【００２０】
サスペンション装置１０は、ショックアブソーバ１１（１１ＦＲ，１１ＦＬ，１１ＲＲ，１１ＲＬ）とサスペンションアクチュエータ１２（１２ＦＲ，１２ＦＬ，１２ＲＲ，１２ＲＬ）を含んで構成されている。ショックアブソーバ１１は、ばね上とばね下との相対運動を減衰させる減衰力を発生させる。サスペンションアクチュエータ１２は、ショックアブソーバ１１が発生させる減衰力（減衰係数）を調節する。サスペンションアクチュエータ１２は、相対的にソフトな減衰特性（減衰力小）から、相対的にハードな減衰特性（減衰力大）まで、ショックアブソーバ１１の減衰特性を任意の特性に変化させることができる。
【００２１】
各車輪２ＦＲ，２ＦＬ，２ＲＲ，２ＲＬには、それぞれの車輪速を検出する車輪速検出手段４（４ＦＲ，４ＦＬ，４ＲＲ，４ＲＬ）が配置されている。右前車輪速検出手段４ＦＲは、右前輪２ＦＲの車輪速を検出する。同様に、左前車輪速検出手段４ＦＬ、右後車輪速検出手段４ＲＲ、左後車輪速検出手段４ＲＬは、それぞれ左前輪２ＦＬ、右後輪２ＲＲ、左後輪２ＲＬの車輪速を検出する。各車輪速検出手段４ＦＲ，４ＦＬ，４ＲＲ，４ＲＬの検出結果を示す信号は、ＥＣＵ１に出力される。
【００２２】
本実施形態のＥＣＵ１は、コンピュータを有する電子制御ユニットである。ＥＣＵ１は、車両状態を推定する各推定手段としての機能を有する。また、ＥＣＵ１は、各車輪２ＦＲ，２ＦＬ，２ＲＲ，２ＲＬのサスペンション装置１０と電気的に接続されており、サスペンション装置１０を制御する。
【００２３】
本実施形態のＥＣＵ１の機能について図２を参照して説明する。図２において、Δω１〜Δω４は、各輪車輪速変動を示す。Δω１は、ＦＲ輪車輪速を示す。Δω２は、ＦＬ輪車輪速を示す。Δω３は、ＲＲ輪車輪速を示す。Δω４は、ＲＬ輪車輪速を示す。ＺＡＦは、前輪ばね下左右輪逆相上下変位を示す。ＺＡＲは、後輪ばね下左右輪逆相上下変位を示す。φＢＧは、ばね上ロール角を示す。ＷＦ，ＷＲは、アクスル位置からボデー重心の左右距離（トレッソの半分）を示す。αＸＦ，αＸＲは、サス上下ストロークに対するアクスル前後変位を示す。αθＦ，αθＲは、サス上下ストロークに対するアクスルピッチ角を示す。ｒＦ，ｒＲは、タイヤ半径（前輪、後輪）を示す。Ｌは、ホイールベースを示す。Ｕは、車速を示す。ｓは、ラプラス演算子を示す。
【００２４】
図２のブロック図に示すように、本実施形態に係るＥＣＵ１は、車輪速検出手段４によって検出した各車輪２ＦＲ，２ＦＬ，２ＲＲ，２ＲＬの車輪速に基づく左右輪逆相の車輪速「１／２・（ΔωＲ−ΔωＬ）」、及び、車両１００の車速を検出する車速検出手段によって検出した車両１００の車速Ｕで所定の前後輪のホイールベースＬを割った値「Ｌ／Ｕ」に基づいて、車両１００のばね下状態量（ＺＡＦ，ＺＡＲ）及び／又はばね上状態量（φＢＧ）を推定する状態量推定手段２１を含んで構成されている。ここで、図２には図示していないが、車両１００の車速を検出する車速検出手段は、ＥＣＵ１内の処理部の一つである。本実施形態において、車速検出手段は、車輪速検出手段４によって検出した各車輪２ＦＲ，２ＦＬ，２ＲＲ，２ＲＬの車輪速に基づいて車両１００の車速を検出してもよいし、車輪速以外の前後加速度の積分値等に基づいて車両１００の車速を検出してもよい。なお、文字の上のドット（・）記号は、微分値を示す。
【００２５】
状態量推定手段２１は、各推定手段２１ａ〜２１ｋを有する。推定手段２１ａは、車輪速検出手段４によって検出した車輪２ＦＲの車輪速（図２において、Δω１ＦＲ輪車輪速）と、車輪速検出手段４によって検出した車輪２ＦＬの車輪速（図２において、Δω２ＦＬ輪車輪速）とを加算する。その際、推定手段２１ａは、車輪２ＦＲの車輪速の値についてはマイナス値（−Δω１）で加算する。推定手段２１ｂは、推定手段２１ａで得られた値「−Δω１＋Δω２」に「１／２」を乗算し、推定手段２１ｃは、推定手段２１ｂで得られた値「（−Δω１＋Δω２）／２」に「１／（ＷＦ（αＸＦ／ｒＦ−αθＦ））」を更に乗算する。
【００２６】
推定手段２１ｄは、車輪速検出手段４によって検出した車輪２ＲＲの車輪速（図２において、Δω３ＲＲ輪車輪速）と、車輪速検出手段４によって検出した車輪２ＲＬの車輪速（図２において、Δω４ＲＬ輪車輪速）とを加算する。その際、推定手段２１ｄは、車輪２ＲＲの車輪速の値についてはマイナス値（−Δω３）で加算する。推定手段２１ｅは、推定手段２１ｄで得られた値「−Δω３＋Δω４」に「１／２」を乗算し、推定手段２１ｆは、推定手段２１ｅで得られた値「（−Δω３＋Δω４）／２」に「１／（ＷＲ（αＸＲ／ｒＲ−αθＲ））」を更に乗算する。
【００２７】
推定手段２１ｇは、推定手段２１ｃで得られた値「｛（−Δω１＋Δω２）／２｝・｛１／（ＷＦ（αＸＦ／ｒＦ−αθＦ））｝」から、推定手段２１ｆで得られた値「｛（−Δω３＋Δω４）／２｝・｛１／（ＷＲ（αＸＲ／ｒＲ−αθＲ））｝」を減算する。そして、推定手段２１ｈは、推定手段２１ｇで得られた値に、「１／｛（１／ＷＦ）−（１／ＷＲ）ｅ＾−Ｌ／Ｕ・ｓ｝」」を乗算して、車両１００のばね下状態量として、車両１００の前輪ばね下左右輪逆相上下変位（ＺＡＦ）を得る。更に、推定手段２１ｉは、推定手段２１ｈで得られた値に、「ｅ＾−Ｌ／Ｕ・ｓ」を乗算して、車両１００のばね下状態量として、車両１００の後輪ばね下左右輪逆相上下変位（ＺＡＲ）を得る。また、推定手段２１ｊは、推定手段２１ｈで得られた値に「１／ＷＦ」を乗算した値から、推定手段２１ｃで得られた値を減算して、車両１００のばね上状態量として、車両１００のばね上ロール角（φＢＧ）を得る。本実施形態のＥＣＵ１は、推定された車両１００のばね下状態量及び／又は車両１００のばね上状態量に基づいて、サスペンション装置１０を制御する。
【００２８】
続いて、本実施形態の状態量推定手段２１による車両１００のばね下状態量（ＺＡＦ，ＺＡＲ）及びばね上状態量（φＢＧ）の推定処理の詳細について説明する。
【００２９】
まず、本実施形態の車両状態量推定装置１０１による車両状態量の推定に係る理論式から詳細を説明する。タイヤ回転速度変動量ωは、下記［数６］によって、算出することができる。なお、タイヤ回転速度変動量ωは、車輪２の回転速度の変動量である。タイヤ回転速度変動量ωは、現在の車速に応じた車輪２の回転速度に対する回転速度の変動量であり、例えば、路面入力やばね上の挙動によって生じる変動量である。言い換えると、タイヤ回転速度変動量ωは、定常的に動いている分に対する変動分（動的な変動分）である。タイヤ回転速度変動量ωは、以下に図３を参照して説明するように、車輪２の半径ｒと、ばね下の前後変位変動量ＸＡによって、［数６］のように表すことができる。
【数６】
0006003931.tif 000007
【００３０】
図３に破線で示した車輪２の位置は、車輪２が定常的に移動した場合の車輪位置、例えば、車速に応じて決まるある所定時刻における車輪位置である。実線で示した車輪２の位置は、所定時刻における実際の車輪位置を示す。路面入力や操作入力により、破線で示す車輪位置と、実線で示す車輪位置とで車両前後方向の位置の変動が生じている。この車両前後方向における車輪位置の変動は、ばね下の前後変位変動量ＸＡに対応している。ばね下の前後変位変動量ＸＡに応じた車輪２の回転角度の変動量βは、図３に示すように、ＸＡ／ｒで近似することができる。タイヤ回転速度変動量ωは、回転角度の変動量βを微分した値であるから、［数６］のように近似的に求めることができる。
【００３１】
ここで、［数６］の右辺における括弧内の第１項について説明する。この第１項は、車軸位置ＴＣ（図４参照）におけるばね上の前後変位変動量ＸＢの微分値である。なお、車軸位置ＴＣとは、各車輪２の中心軸線上でかつ車輪２の幅方向の中央の点とする。車軸位置ＴＣにおけるばね上の前後変位変動量ＸＢは、図４に示すばね上（重心位置ＰＧ）の前後変位変動量ＸＢＧと、図５に示すばね上ピッチ角θＢＧに基づく変動量と、図６に示すばね上ヨー角ΨＢＧに基づく変動量との和として表すことができる。
【００３２】
図４に示すように、ばね上が車両前後方向に位置変動する場合、車軸位置ＴＣで見たばね上の前後変位変動量ＸＢには、ばね上の車両１００の前後方向の位置変動による成分ＸＢａが含まれる。この車両１００の前後方向の位置変動による成分ＸＢａは、下記式（１）で示すように、重心位置ＰＧにおけるばね上の前後変位変動量ＸＢＧとして表される。
ＸＢａ＝ＸＢＧ…（１）
【００３３】
図５に示すように車両１００にピッチングが生じた場合、車軸位置ＴＣで見たばね上の前後変位変動量ＸＢには、ピッチングによる成分ＸＢｂが含まれる。このピッチングによる成分ＸＢｂは、ばね上ピッチ角θＢＧに基づいて、下記式（２）として表される。ここで、Ｈ：車軸位置ＴＣと重心位置ＰＧとの高さ方向における距離、である。
ＸＢｂ＝−θＢＧ×Ｈ…（２）
【００３４】
図６に示すように車両１００にヨー方向の挙動（鉛直軸回りの回転）が生じた場合、車軸位置ＴＣで見たばね上の前後変位変動量ＸＢには、ヨー方向の挙動による成分ＸＢｃが含まれる。ヨー方向の挙動による成分ＸＢｃは、ばね上ヨー角ΨＢＧに基づいて、下記式（３）として表される。ここで、Ｗ：車幅方向における重心位置ＰＧと車軸位置ＴＣとの距離、である。
ＸＢｃ＝ΨＢＧ×Ｗ…（３）
【００３５】
以上の３つの成分ＸＢａ、ＸＢｂ、ＸＢｃの和が、車軸位置ＴＣで見たばね上の前後変位変動量ＸＢである。つまり、下記式（４）が導かれる。
ＸＢ＝−θＢＧ×Ｈ＋ＸＢＧ＋ΨＢＧ×Ｗ…（４）
【００３６】
次に、［数６］の右辺における括弧内の第２項について説明する。この第２項は、ばね上とばね下との相対変位、すなわちサスペンション装置１０の上下ストロークに係る項である。サスペンション装置１０のストロークにより、図７に示すように、車輪２と、ばね上であるボデー３とが上下方向に相対変位する。ここで、サスペンション装置１０が上下方向に対して車両前後方向に傾斜している場合、車輪２と、ボデー３とは車両前後方向においても相対変位する。ばね上とばね下との前後変位変動量の差（ＸＡ−ＸＢ）は、下記式（５）で表すことができる。
ＸＡ−ＸＢ＝αＸ（ＺＡ−ＺＢ）…（５）
ここで、ＺＡ：ばね下の上下変位変動量、ＺＢ：車両前後方向の車軸位置ＴＣにおけるばね上の上下変位変動量、αＸ：サスペンション装置１０の単位ストローク量あたりのばね上とばね下との車両前後方向の相対変位量、である。
【００３７】
上記［数６］及び上記式（１）乃至式（５）から、下記［数７］が導かれる。
【数７】
0006003931.tif 000008
【００３８】
次に、図８及び図９を参照して、ばね下ピッチ角θＡについて説明する。ばね下ピッチ角θＡは、ばね下のピッチ角であり、路面に対するばね下の車両前後方向における傾斜角を示す。車両１００では、サスペンション装置１０の伸縮により、ばね上ピッチ角θＢＧとばね下ピッチ角θＡとが異なる場合がある。言い換えると、ばね下ピッチ角θＡは、ばね上ピッチ角θＢＧに対応する成分θＡａと、サスペンション装置１０の伸縮による成分θＡｂとを含んでいる。
【００３９】
図８に示すように、ばね上ピッチ角θＢＧに対応する成分θＡａは、ばね上ピッチ角θＢＧと等しい。つまり、ばね上ピッチ角θＢＧに対応する成分θＡａは、下記式（６）で表される。
θＡａ＝θＢＧ…（６）
図９には、路面入力等によって車両１００の前輪２ＦＲ，２ＦＬのサスペンション装置１０ＦＲ，１０ＦＬが伸縮し、後輪２ＲＲ，２ＲＬのサスペンション装置１０ＲＲ，１０ＲＬは伸縮していない状態が示されている。サスペンション装置１０の伸縮による成分θＡｂは、図９に示すように、下記式（７）で表される。
θＡｂ＝−αθ（ＺＡ−ＺＢ）…（７）
ここで、αθ：サスペンション装置１０の単位ストローク量あたりのばね下ピッチ角、である。
【００４０】
上記式（６）及び式（７）より、ばね下ピッチ角θＡの角速度変動（ばね下ピッチ角θＡの微分値）は、下記［数８］のように表される。
【数８】
0006003931.tif 000009
【００４１】
次に、車輪速変動量Δωについて説明する。車輪速変動量Δωは、車輪速検出手段４によって検出される車輪速の変動量である。車輪速変動量Δωは、下記［数９］で表される。車輪速検出手段４によって検出される車輪速は、路面に対する前後方向の相対移動による回転速度成分だけでなく、ばね下ピッチ角θＡの変化による回転速度成分を含んでいる。すなわち、車輪速変動量Δωは、［数９］のように、タイヤ回転速度変動量ωと、ばね下ピッチ角θＡの角速度変動との差分となる。ＥＣＵ１は、車輪速検出手段４によって検出された車輪速から現在の車速に対応する車輪速を除いて車輪速変動量Δωを算出する車輪速変動量推定手段を有している。
【数９】
0006003931.tif 000010
【００４２】
上記［数７］、［数８］を［数９］に代入すると、下記［数１０］が導かれる。
【数１０】
0006003931.tif 000011
【００４３】
次に、上下変位変動量について説明する。車軸位置ＴＣにおけるばね上の上下変位変動量ＺＢは、下記式（８）で表すことができる。下記式（８）の右辺第１項は、ばね上の重心位置ＰＧの上下変位変動量である。また、式（８）の右辺第２項は、ばね上のピッチングによる上下変位変動量（近似値）である。また、式（８）の右辺第３項は、ばね上のロールによる上下変位変動量（近似値）である。
ＺＢ＝ＺＢＧ±Ｌ×θＢＧ±Ｗ×φＢＧ…（８）
ここで、Ｌ：車両前後方向における車軸位置ＴＣと重心位置ＰＧとの距離（図４参照）、φＢＧ：ばね上ロール角、である。なお、車軸位置ＴＣと重心位置ＰＧとの距離Ｌについては、前輪２ＦＲ，２ＦＬの場合の当該距離ＬＦと後輪２ＲＲ，２ＲＬの場合の当該距離ＬＲとが異なる値であってもよい。
【００４４】
上記式（８）について、右前輪２ＦＲを例に説明すると、右前輪２ＦＲの車軸位置ＴＣにおけるばね上の上下変位変動量ＺＢ１は、重心位置ＰＧの上下変位変動量ＺＢＧに、ばね上のピッチングによる上下変位変動量とばね上のロールによる上下変位変動量を加算又は減算して求めることができる。例えば、前輪側が沈み込むようにばね上がピッチングする場合、重心位置ＰＧの上下変位変動量ＺＢＧからばね上のピッチングによる上下変位変動量（Ｌ×θＢＧ）が減算される。また、車両右側が沈み込むようにばね上がロールする場合、ばね上のロールによる上下変位変動量（Ｗ×φＢＧ）が減算される。
【００４５】
これとは逆に、前輪側が浮き上がるようにばね上がピッチングする場合、重心位置ＰＧの上下変位変動量ＺＢＧに対して、ばね上のピッチングによる上下変位変動量（Ｌ×θＢＧ）が加算される。また、車両右側が浮き上がるようにばね上がロールする場合、ばね上のロールによる上下変位変動量（Ｗ×φＢＧ）が加算される。他の車輪２ＦＬ，２ＲＲ，２ＲＬについても同様にして車軸位置ＴＣにおけるばね上の上下変位変動量ＺＢが算出される。
【００４６】
上記式（８）から、車軸位置ＴＣにおけるばね上の上下変位変動速度は、下記［数１１］のように表される。
【数１１】
0006003931.tif 000012
【００４７】
上記［数６］から［数１１］を拡張し、４輪分の車輪速変動量をマトリクスで表記すると、下記［数１２］が得られる。なお、マトリクス［Ｄ］は下記［数１３］、マトリクス［Ｇ］は下記［数１４］、マトリクス［Ｅ］は下記［数１５］、マトリクス［Ｆ］は下記［数１６］、マトリクス［Ｃ］は下記［数１７］で示すものである。
【数１２】
0006003931.tif 000013
【数１３】
0006003931.tif 000014
【数１４】
0006003931.tif 000015
【数１５】
0006003931.tif 000016
【数１６】
0006003931.tif 000017
【数１７】
0006003931.tif 000018
【００４８】
ここで、Δω１：右前輪２ＦＲの車輪速変動量、Δω２：左前輪２ＦＬの車輪速変動量、Δω３：右後輪２ＲＲの車輪速変動量、Δω４：左後輪２ＲＬの車輪速変動量である。すなわち、各変数の添字１は右前輪２ＦＲに係る値、添字２は左前輪２ＦＬに係る値、添字３は右後輪２ＲＲに係る値、添字４は左後輪２ＲＬに係る値を示す。
【００４９】
その他の変数は、以下のようである。
ω１，ω２，ω３，ω４：右前、左前、右後、左後の各車輪２のタイヤ回転速度変動量
θＡ１，θＡ２，θＡ３、θＡ４：右前、左前、右後、左後の各車輪２の位置のばね下ピッチ角
ＸＡ１，ＸＡ２，ＸＡ３，ＸＡ４：右前、左前、右後、左後の各車輪２のばね下前後変位変動量
ＸＢ１，ＸＢ２，ＸＢ３，ＸＢ４：右前、左前、右後、左後の各車輪２の車軸位置ＴＣにおけるばね上の前後変位変動量
ＺＡ１，ＺＡ２，ＺＡ３，ＺＡ４：右前、左前、右後、左後の各車輪２のばね下の上下変位変動量
ＺＢ１，ＺＢ２，ＺＢ３，ＺＢ４：右前、左前、右後、左後の各車輪２の車軸位置ＴＣにおけるばね上の上下変位変動量
ＹＢＧ：ばね上の重心位置ＰＧの左右変位変動量
ｒＦ：前輪２ＦＲ，２ＦＬの半径
ｒＲ：後輪２ＲＲ，２ＲＬの半径
ＬＦ：前輪２ＦＲ，２ＦＬの車軸位置ＴＣと重心位置ＰＧとの車両前後方向の距離
ＬＲ：後輪２ＲＲ，２ＲＬの車軸位置ＴＣと重心位置ＰＧとの車両前後方向の距離
ＷＦ：前輪２ＦＲ，２ＦＬの車軸位置ＴＣと重心位置ＰＧとの左右方向（車幅方向）距離
ＷＲ：後輪２ＲＲ，２ＲＬの車軸位置ＴＣと重心位置ＰＧとの左右方向（車幅方向）距離
ＨＦ：前輪２ＦＲ，２ＦＬの車軸位置ＴＣと重心位置ＰＧとの上下方向の距離
ＨＲ：後輪２ＲＲ，２ＲＬの車軸位置ＴＣと重心位置ＰＧとの上下方向の距離
αＸＦ：サスペンション装置１０ＦＲ，ＦＬの単位ストローク量あたりの各前輪２ＦＲ，２ＦＬにおけるばね上とばね下との車両前後方向の相対変位量
αＸＲ：サスペンション装置１０ＲＲ，ＲＬの単位ストローク量あたりの各後輪２ＲＲ，２ＲＬにおけるばね上とばね下との車両前後方向の相対変位量
αθＦ：サスペンション装置１０ＦＲ，ＦＬの単位ストローク量あたりの各前輪２ＦＲ，２ＦＬにおけるばね下ピッチ角
αθＲ：サスペンション装置１０ＲＲ，ＲＬの単位ストローク量あたりの各後輪２ＲＲ，２ＲＬにおけるばね下ピッチ角
【００５０】
次に、上記［数１２］を展開すると、下記［数１８］が得られる。
【数１８】
0006003931.tif 000019
ここで、Δω１乃至Δω４は、其々右前、左前、右後、左後の各車輪２の車輪速を示す。なお、マトリクス［Ｄ］は上記［数１３］、マトリクス［Ｇ］は上記［数１４］、マトリクス［Ｅ］は上記［数１５］、マトリクス［Ｆ］は上記［数１６］、マトリクス［Ｃ］は上記［数１７］で示すものである。ＺＢＧ，θＢＧ，ＸＢＧ，ＹＢＧ，φＢＧ，ΨＢＧは、其々６方向、すなわち上下方向（ＺＢＧ）、ばね上ピッチ角方向（θＢＧ）、車両前後方向（ＸＢＧ）、車両左右方向（ＹＢＧ）、ばね上ロール角方向（φＢＧ）及びばね上ヨー角方向（ΨＢＧ）の各方向についてのばね上の上下変位変動量を示す。ＺＡ１〜ＺＡ４は、其々右前ＦＲ、左前ＦＬ、右後ＲＲ、左後ＲＬの各車輪のばね下の上下変位変動量を示す。
【００５１】
そして、上記［数１８］に示す４輪分の車輪速変動量｛Δω１，Δω２，Δω３，Δω４｝から、前輪位置の左右輪逆相の車輪速「ΔωＦ＝（−Δω１＋Δω２）／２」及び後輪位置の左右輪逆相の車輪速「ΔωＲ＝（−Δω３＋Δω４）／２」がそれぞれ得られる。
【００５２】
前輪位置の左右輪逆相の車輪速「ΔωＦ＝（−Δω１＋Δω２）／２」は、右前ＦＲの車輪速Δω１と左前ＦＬの車輪速Δω２の差を２で割った値であり、具体的には、上記［数１８］の４行のマトリクスのうち上から１行目と２行目の差を２で割った値となる。後輪位置の左右輪逆相の車輪速「ΔωＲ＝（−Δω３＋Δω４）／２」は、右後ＲＲの車輪速Δω３と左後ＲＬの車輪速Δω４の差を２で割った値であり、具体的には、上記［数１８］の４行のマトリクスのうち上から３行目と４行目の差を２で割った値となる。
【００５３】
ここで、左右輪逆相入力に対しては、ＺＢＧ＝θＢＧ＝ＸＢＧ＝０であり、また直進時のΨＢＧ≒０であることを考慮し、上記［数１８］の１行目と２行目、３行目と４行目の差から左右輪逆相成分を演算すると、下記［数１９］が得られる。
【数１９】
0006003931.tif 000020
具体的には、左右輪逆相入力を考えると、上記［数１８］に示すように、ばね上の動きは６個｛ＺＢＧ，θＢＧ，ＸＢＧ，ＹＢＧ，φＢＧ，ΨＢＧ｝あるが、これらのうち３個｛ＺＢＧ，θＢＧ，ＸＢＧ｝は発生しなくなる。これは、左右輪逆相入力を考えると、車両１００の真ん中の位置は上下には動かなくなり、またピッチング方向にも動かなくなり、前後方向にも動かなくなるからである。そのため、上下方向（ＺＢＧ）、ばね上ピッチ角方向（θＢＧ）、車両前後方向（ＸＢＧ）についてのばね上の動きについては０であると考えることができる。よって、逆相入力で発生するのは、車両左右方向（ＹＢＧ）、ばね上ロール角方向（φＢＧ）及びばね上ヨー角方向（ΨＢＧ）の３個のばね上の動きとなる。また、路面入力の４方向の動きについても非常に微小であるため、ばね上ヨー角方向（ΨＢＧ）についても０とみなすことできる。従って、上記［数１８］に示したばね上の動きは６個｛ＺＢＧ，θＢＧ，ＸＢＧ，ＹＢＧ，φＢＧ，ΨＢＧ｝のうち４個｛ＺＢＧ，θＢＧ，ＸＢＧ，ΨＢＧ｝を０とすることができる。そのため、前輪位置の左右輪逆相の車輪速「ΔωＦ＝（−Δω１＋Δω２）／２」及び後輪位置の左右輪逆相の車輪速「ΔωＲ＝（−Δω３＋Δω４）／２」は、其々上記［数１９］のように表すことができる。
【００５４】
上記［数１９］において、未知数は、φＢＧと、｛（−ＺＡ１＋ＺＡ２）／２｝と、｛（−ＺＡ３＋ＺＡ４）／２｝との３つになる。上記［数１９］では、３つの未知数に対し、２つの式が存在するため、未知数を１つ減らして２つにする必要がある。
【００５５】
前輪入力と後輪入力にはホイールベースＬ／車速Ｕの入力時間差があるため、これを考慮すると、上記［数１９］の未知数のうち、｛（−ＺＡ１＋ＺＡ２）／２｝と、｛（−ＺＡ３＋ＺＡ４）／２｝については、下記［数２０］に示すように、其々「ＺＡＦ」と「ＺＡＦｅ＾−（Ｌ／Ｕ・ｓ）」に置き換えることができる。
【数２０】
0006003931.tif 000021
【００５６】
そうすると、上記［数２０］において、未知数は、φＢＧ及びＺＡＦの２つになる。上記［数２０］では、２つの未知数に対して２つの式が存在するため、各未知数を上記［数２０］から求めると、前輪ばね下上下速度の左右輪逆相成分ＺＡＦ、後輪ばね下上下速度の左右輪逆相成分ＺＡＲ、ばね上ロール角速度φＢＧは、其々下記［数２１］，［数２２］，［数２３］で表される。その結果、車両モデルを用いずにばね下、ばね上（ロール）の状態量を推定可能となる。
【数２１】
0006003931.tif 000022
【数２２】
0006003931.tif 000023
【数２３】
0006003931.tif 000024
【００５７】
以上説明したように、状態量推定手段２１は、左右輪逆相の車輪速として、前輪位置の左右輪逆相の車輪速「ΔωＦ＝（−Δω１＋Δω２）／２」及び後輪位置の左右輪逆相の車輪速「ΔωＲ＝（−Δω３＋Δω４）／２」を表した上記［数１９］を、車速Ｕで所定の前後輪のホイールベースＬを割った値「Ｌ／Ｕ」を用いて上記［数２０］に変換する。そして、状態量推定手段２１は、上記［数２０］から、ばね下状態量として、上記［数２１］に示す車両１００の前輪ばね下左右輪逆相上下変位「ＺＡＦ」、及び、上記［数２２］に示す車両１００の後輪ばね下左右輪逆相上下変位「ＺＡＲ」を推定する。また、状態量推定手段２１は、ばね上状態量として、上記［数２３］に示す車両１００のばね上ロール角「φＢＧ」を推定する。
【００５８】
このように本実施形態の状態量推定手段２１は、左右輪逆相の車輪速「１／２・（ΔωＲ−ΔωＬ）」及び前後輪のホイールベース／車速「Ｌ／Ｕ」に基づいて、ばね下・ばね上挙動を推定することができる。この状態量推定手段２１による推定方法は、車両特性により変化するばね上上下変位変動（ＺＢＧ）・ばね上ピッチ角（θＢＧ）・ばね上前後変位変動（ＸＢＧ）と相関のない左右輪逆相の車輪速に基づく推定方法であるため、車両特性の変化によらず、ばね下・ばね上挙動を高精度に推定可能となる。実際には、上記［数２１］により「前輪ばね下左右輪逆相上下変位」を推定し、上記［数２２］により「後輪ばね下左右輪逆相上下変位」を推定し、上記［数２３］により「ばね上ロール角」を推定できる。
【００５９】
次に、本実施形態に係る車両状態推定装置１０１による車両状態の推定方法と上記理論式の対応関係等について説明する。上記理論式［数２１］，［数２２］，［数２３］と、上述の図２のブロック図とを対応付けると、図１０に示すようになる。
【００６０】
図１０中の［数２１］に対応する部分において、状態量推定手段２１の推定手段２１ｇは、状態量推定手段２１の推定手段２１ｃで得られた値「｛（−Δω１＋Δω２）／２｝・｛１／（ＷＦ（αＸＦ／ｒＦ−αθＦ））｝」から、状態量推定手段２１の推定手段２１ｆで得られた値「｛（−Δω３＋Δω４）／２｝・｛１／（ＷＲ（αＸＲ／ｒＲ−αθＲ））｝」を減算する。そして、状態量推定手段２１の推定手段２１ｈは、推定手段２１ｇで得られた値に、「１／｛（１／ＷＦ）−（１／ＷＲ）ｅ＾−Ｌ／Ｕ・ｓ｝」」を乗算して、車両１００のばね下状態量として、車両１００の前輪ばね下左右輪逆相上下変位（ＺＡＦ）を得る。得られた車両１００の前輪ばね下左右輪逆相上下変位（ＺＡＦ）は、上記［数２１］に示すものである。
【００６１】
図１０中の［数２２］に対応する部分において、状態量推定手段２１の推定手段２１ｉは、推定手段２１ｈで得られた値に、「ｅ＾−Ｌ／Ｕ・ｓ」を乗算して、車両１００のばね下状態量として、車両１００の後輪ばね下左右輪逆相上下変位（ＺＡＲ）を得る。得られた車両１００の後輪ばね下左右輪逆相上下変位（ＺＡＲ）は、上記［数２２］に示すものである。
【００６２】
図１０中の［数２３］に対応する部分においては、状態量推定手段２１の推定手段２１ｊは、推定手段２１ｈで得られた値に「１／ＷＦ」を乗算した値から、推定手段２１ｃで得られた値を減算して、車両１００のばね上状態量として、車両１００のばね上ロール角（φＢＧ）を得る。得られた車両１００のばね上ロール角（φＢＧ）は、上記［数２３］に示すものである。
【００６３】
本実施形態において、車両状態推定装置１０１は、上述の状態量推定手段２１によって推定した前輪ばね下左右輪逆相上下変位（ＺＡＦ）、後輪ばね下左右輪逆相上下変位（ＺＡＲ）、及び、ばね上ロール角（φＢＧ）に基づいて、ロール共振周波数（ωｎ）、及び、ばね下状態量に対するばね上状態量の伝達関数最大値（ｐ）を算出し、算出したロール共振周波数（ωｎ）及び伝達関数最大値（ｐ）に基づいて、車両１００の積載状態を推定する積載状態推定手段、を更に備える。本実施形態において、車両１００の積載状態は、積載量（Ｍ）、ばね上重心高（Ｈ）を含む。
【００６４】
具体的には、積載状態推定手段は、下記［数２４］に示すように、ばね下左右輪逆相上下変位（ＺＡＦ，ＺＡＲ）、及び、ばね上ロール角（φＢＧ）に基づいて、ばね下に対するばね上（ロール）の伝達関数を演算する。
【数２４】
0006003931.tif 000025
【００６５】
ここで、上記［数２４］の左辺の上の項（∂φＢＧ／∂ＺＡＦ）は、前輪入力のみによる車両１００のロールの動きを表す。また、上記［数２４］の左辺の下の項（∂φＢＧ／∂ＺＡＲ）は、後輪入力のみによる車両１００のロールの動きを表す。上記［数２４］の左辺の∂は偏微分記号を示す。上記［数２４］に示す伝達関数では、前輪の路面入力でも車両１００はロールし、後輪の路面入力でも車両１００がロールするため、この２つの入力を合成した形で車両１００のロールの動きを重回帰分析の式にて表している。なお、上記［数２４］に示す伝達関数の算出概念は、文献「実走行時の車両の重心地解析技術；吉見剛、香村伸吾、塩谷崇洋著；自動車技術会春季学術講演会前刷集、６９−１３号、２７−３２頁、２０１３年５月２２日発行」に従うものである。
【００６６】
そして、積載状態推定手段は、演算した上記［数２４］に示す伝達関数から決定される伝達関数の共振周波数（ωｎ）や最大値（ｐ）等と、積載状態との関係から、積載状態を推定する。ここで、図１１に示すように、積載量が増加し重心高さが上がると、ロール共振周波数（ωｎ）は低くなり、積載量が減少し重心高さが下がると、ロール共振周波数（ωｎ）は高くなる。また、図１２に示すように、積載量が増加し重心高さが上がると、伝達関数最大値（ｐ）は小さくなり、積載量が減少し重心高さが下がると、伝達関数最大値（ｐ）は大きくなる。図１１及び図１２に示したロール共振周波数（ωｎ）と伝達関数最大値（ｐ）を積載量（Ｍ）と重心高（Ｈ）の関係を数式化すると、下記［数２５］となる。
【数２５】
0006003931.tif 000026
【００６７】
ここで、上記［数２５］において、ωｎ：ロール共振周波数、ｐ：伝達関数最大値、Ｍ：積載量、Ｈ：ばね上重心高（空車時との差）、∂ωｎ／∂Ｍ：積載量に対するロール共振周波数の感度、∂ωｎ／∂Ｈ：重心高に対するロール共振周波数の感度、∂ｐ／∂Ｍ：積載量に対する伝達関数最大値の感度、∂ｐ／∂Ｈ：重心高に対する伝達関数最大値の感度、ωｎ０：空車時ロール共振周波数、ｐ０：空車時伝達関数最大値、を示す。
【００６８】
積載状態推定手段は、上記［数２４］で計算された伝達関数のロール共振周波数（ωｎ）と伝達関数最大値（ｐ）から、上記［数２５］を用いて、下記［数２６］に示すように、車両１００の積載状態として、積載量（Ｍ）とばね上重心高（Ｈ）を逆算することで推定する。
【数２６】
0006003931.tif 000027
【００６９】
このように本実施形態の積載状態推定手段は、前輪ばね下左右輪逆相上下変位（ＺＡＦ）、後輪ばね下左右逆相上下変位（ＺＡＲ）、ばね上ロール角（φＢＧ）から演算したロール共振周波数（ωｎ）とばね下に対するばね上の伝達関数最大値（ｐ）とに基づいて、積載状態（積載量（Ｍ）とばね上重心高（Ｈ））を推定することができる。従って、本実施形態によれば、追加のセンサ等を使用することなく、検出される車輪速と車速から車両１００の積載状態を高精度に推定することができる。
【００７０】
また、本実施形態において、車両制御装置１０２は、状態量推定手段２１によって推定したばね下状態量（前輪ばね下左右輪逆相上下変位（ＺＡＦ）、後輪ばね下左右逆相上下変位（ＺＡＲ））及びばね上状態量（ばね上ロール角（φＢＧ））に基づいて車両１００のサスペンション装置１０を制御する制御手段を備える。具体的には、制御手段は、上述のように積載状態推定手段により推定された車両１００の積載状態（積載量（Ｍ）とばね上重心高（Ｈ））に基づいて車両１００のサスペンション装置１０を制御する。
【００７１】
以上説明したように、本実施形態の車両状態推定装置１０１及び車両制御装置１０２は、以下に示す車両状態推定方法を実行することができる。
【００７２】
車両状態推定装置１０１及び車両制御装置１０２が実行する車両状態推定方法は、車両１００の各車輪の車輪速を検出する車輪速検出手順と、車両１００の車速を検出する車速検出手順と、車輪速検出手順において検出した各車輪の車輪速に基づく左右輪逆相の車輪速、及び、車速検出手順において検出した車両１００の車速で所定の前後輪のホイールベースを割った値に基づいて、車両１００のばね下状態量及び／又はばね上状態量を推定する状態量推定手順と、を含む。また、車両状態推定方法は、状態量推定手順において推定した前輪ばね下左右輪逆相上下変位、後輪ばね下左右輪逆相上下変位、及び、ばね上ロール角に基づいて、ロール共振周波数、及び、ばね下状態量に対するばね上状態量の伝達関数最大値を算出し、算出したロール共振周波数及び伝達関数最大値に基づいて、車両１００の積載状態を推定する積載状態推定手順、を更に含む。本実施形態では、車輪速検出手段４が車輪速検出手順を実行する。状態量推定手段２１は、状態量推定手順を実行する。積載状態推定手段は、積載状態推定手順を実行する。
【００７３】
上記実施形態では、理論式において近似式を用いているものもあるが、近似式に代えてより詳細（高次）の式が用いられてもよい。例えば、車輪２のタイヤの変形等を考慮した高次の式が用いられてもよい。
【００７４】
上記実施形態では、制御手段としてのＥＣＵ１は、ばね下状態量及びばね上状態量に基づいてサスペンション装置１０を制御したが、制御に用いられる状態量は、これには限定されない。制御手段は、ばね下状態量及びばね上状態量の少なくとも何れか一方に基づいてサスペンション装置１０を制御してもよい。
【００７５】
ばね下状態量及び／又はばね上状態量に基づく制御の対象は、サスペンション装置１０以外のものであってもよい。車両１００の挙動に基づいて制御される他の装置、例えばブレーキ装置や加減速装置や操舵装置が制御手段によって制御されてもよい。
【符号の説明】
【００７６】
１  ＥＣＵ
２  車輪
３  ばね上（ボデー）
４  車輪速検出手段
１０  サスペンション装置
１１  ショックアブソーバ
１２  サスペンションアクチュエータ
２１  状態量推定手段
１００  車両
１０１  車両状態推定装置
１０２  車両制御装置
ＰＧ  重心位置
ＴＣ  車軸位置
【図１】
0006003931.tif 000043
【図２】
0006003931.tif 000044
【図３】
0006003931.tif 000045
【図４】
0006003931.tif 000046
【図５】
0006003931.tif 000047
【図６】
0006003931.tif 000048
【図７】
0006003931.tif 000049
【図８】
0006003931.tif 000050
【図９】
0006003931.tif 000051
【図１０】
0006003931.tif 000052
【図１１】
0006003931.tif 000053
【図１２】
0006003931.tif 000054