(57)【要約】【課題】車両の運動状態の推定を精度高く行なうことができる状態推定装置、サスペンション制御装置及びサスペンションシステムを提供する。【解決手段】車両２からのばね上速度（観測出力）と実車両モデル状態量推定用オブザーバ２３の車両近似モデルからの推定ばね上速度（推定観測出力）との差に相当する出力偏差から、実減衰力推定用オブザーバ２１の第１オブザーバゲインが動特性補償信号を算出し、実車両モデル状態量推定用オブザーバ２３の第２オブザーバゲインが車両モデル補償信号を算出する。動特性補償信号は、実車両モデル状態量推定用オブザーバ２３の動特性付与手段に入力され、動特性付与手段の設定内容の調整に用いられる。このため、制御上の時間遅れを発生が抑制され、良好な精度をもって振動抑制を行うことができる。【選択図】図１
2010058541.tif 000002
【特許請求の範囲】
【請求項１】
ばね上とばね下との間に制御力を発生する制御サスペンションを有して第１運動状態量を出力する車両に用いられ、前記車両の第２運動状態量を推定する状態推定装置であって、
前記制御サスペンションへの入力から前記制御サスペンションの静特性に基づいて推定静的制御力を算出する制御力推定手段と、
前記推定静的制御力を入力とするオブザーバとからなり、
前記オブザーバは、
前記推定静的制御力から推定動的制御力を算出する動特性付与手段と、
前記推定動的制御力から前記車両の推定第１運動状態量及び推定第２運動状態量を算出する車両モデルと、
前記車両からの前記第１運動状態量と前記車両モデルからの前記推定第１運動状態量との差に相当する出力偏差から動特性補償信号を算出する動特性ゲイン手段と、
前記出力偏差から車両モデル補償信号を算出する車両モデルゲイン手段とを有し、
前記動特性補償信号を前記動特性付与手段に、前記車両モデル補償信号を車両モデルに、それぞれフィードバックすることを特徴とする状態推定装置。
【請求項２】
請求項１に記載の状態推定装置において、制御サスペンションが、セミアクティブダンパであることを特徴とする状態推定装置。
【請求項３】
請求項２に記載の状態推定装置において、制御力は減衰力であり、制御力推定手段は、指令信号及び推定ピストン速度を入力信号として受けることにより推定静的制御力に相当する推定静的減衰力を算出することを特徴とする状態推定装置。
【請求項４】
請求項２に記載の状態推定装置において、第１運動状態量が、ばね上速度及びばね上・ばね下相対速度のうち一方であり、第２運動状態量が他方であることを特徴とする状態推定装置。
【請求項５】
請求項１に記載の状態推定装置において、制御サスペンションが、アクティブサスペンションであることを特徴とする状態推定装置。
【請求項６】
請求項１に記載の状態推定装置において、制御サスペンションが、ロール剛性可変装置であることを特徴とする状態推定装置。
【請求項７】
請求項１に記載の状態推定装置において、動特性をｎ次遅れ系で構成したことを特徴とする状態推定装置。
【請求項８】
請求項７に記載の状態推定装置において、ｎ次遅れ系は１次遅れ系であることを特徴とする状態推定装置。
【請求項９】
制御力が可変の制御サスペンションに指令信号を送るサスペンション制御装置であって、
車両が出力する当該車両の運動状態量に相当する第１運動状態量以外の当該車両の第２運動状態量を少なくとも１つ推定する第２運動状態量推定手段と、
前記第１運動状態量および第２運動状態量から前記指令信号を算出する車両振動制御器とからなり、
前記第２運動状態量推定手段は、
前記制御サスペンションの静特性に基づいて推定静的制御力を算出する制御力推定手段と、
前記推定静的制御力を入力変数とするオブザーバとからなり、
前記オブザーバは、
前記推定静的制御力から推定動的制御力を算出する動特性付与手段と、
前記推定動的制御力から前記車両の推定第１運動状態量及び推定第２運動状態量を算出する車両モデルと、
前記車両からの前記第１運動状態量と前記車両モデルからの前記推定第１運動状態量との差に相当する出力偏差から動特性補償信号を算出する動特性ゲイン手段と、
前記出力偏差から車両モデル補償信号を算出する車両モデルゲイン手段とを有し、
前記動特性補償信号を前記動特性付与手段に、前記車両モデル補償信号を車両モデルに、それぞれフィードバックすることを特徴とするサスペンション制御装置。
【請求項１０】
請求項９に記載のサスペンション制御装置において、車両振動制御器がスカイフック制御器又はＨ∞制御器であることを特徴とするサスペンション制御装置。
【請求項１１】
請求項９に記載のサスペンション制御装置において、制御サスペンションが、セミアクティブダンパであることを特徴とするサスペンション制御装置。
【請求項１２】
請求項９に記載のサスペンション制御装置において、制御サスペンションが、アクティブサスペンションであることを特徴とするサスペンション制御装置。
【請求項１３】
請求項９に記載のサスペンション制御装置において、制御サスペンションが、ロール剛性可変装置あることを特徴とするサスペンション制御装置。
【請求項１４】
請求項９に記載のサスペンション制御装置において、動特性をｎ次遅れ系で構成した特徴とするサスペンション制御装置。
【請求項１５】
制御力が可変の制御サスペンションと、
前記制御サスペンションに指令信号を送るサスペンション制御装置とからなるサスペンションシステムであって、
前記サスペンション制御装置は、車両が出力する当該車両の運動状態量に相当する第１運動状態量以外の当該車両の第２運動状態量を少なくとも１つ推定する第２運動状態量推定手段と、
前記第１運動状態量および第２運動状態量から前記指令信号を算出する車両振動制御器とからなり、
前記第２運動状態量推定手段は、
前記制御サスペンションの静特性に基づいて推定静的制御力を算出する制御力推定手段と、
前記推定静的制御力を入力変数とするオブザーバとからなり、
前記オブザーバは、
前記推定静的制御力から推定動的制御力を算出する動特性付与手段と、
前記推定動的制御力から、前記車両の推定第１運動状態量及び推定第２運動状態量を算出する車両モデルと、
前記車両からの前記第１運動状態量と前記車両モデルからの前記推定第１運動状態量との差に相当する出力偏差から動特性補償信号を算出する動特性ゲイン手段と、
前記出力偏差から車両モデル補償信号を算出する車両モデルゲイン手段とを有し、
前記動特性補償信号を前記動特性付与手段に、前記車両モデル補償信号を車両モデルに、それぞれフィードバックすることを特徴とするサスペンションシステム。
【請求項１６】
請求項１５に記載のサスペンションシステムにおいて、車両振動制御器がスカイフック制御器又はＨ∞制御器であることを特徴とするサスペンションシステム。
【請求項１７】
請求項１５に記載のサスペンションシステムにおいて、制御サスペンションが、セミアクティブダンパであることを特徴とするサスペンションシステム。
【請求項１８】
請求項１５に記載のサスペンションシステムにおいて、制御サスペンションが、アクティブサスペンションであることを特徴とするサスペンションシステム。
【請求項１９】
請求項１５に記載のサスペンションシステムにおいて、制御サスペンションが、ロール剛性可変装置あることを特徴とするサスペンションシステム。
【請求項２０】
請求項１５に記載のサスペンションシステムにおいて、動特性をｎ次遅れ系で構成した特徴とするサスペンションシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、車両の振動制御に用いられる状態推定装置、サスペンション制御装置及びサスペンションシステムに関する。
【背景技術】
【０００２】
車両の振動制御に用いられるサスペンション制御装置及びサスペンションシステムなどに、車両の運動状態推定などのためにオブザーバを用いることがある。そして、このように、オブザーバにより車両の運動状態を推定するためには、一般に、車両の線形近似モデルを作成する必要がある。また、このモデルは可検出である必要があり、さらに、車両近似モデルに入力する推定制御力が実際の車両に入力する実制御力と同一である必要がある。これらの条件を満たさない場合、推定精度が低下、もしくは、推定不能となる場合もある。
特に、サスペンション制御装置及びサスペンションシステムなどに用いられるセミアクティブダンパやアクティブサスペンションなどの制御力が可変の制御サスペンションは、非線形特性を有するため、推定制御力を実制御力と同一にすることは難しい。
【０００３】
この解決策として、特許文献１に示される装置がある。特許文献１に示される装置は、予め制御サスペンションの制御力〔例えばセミアクティブダンパが発生する減衰力及びアクティブサスペンションが発生する力（抵抗力に相当する減衰力に比して言わば能動的な力と称せられる力）など〕の非線形特性を実験により把握し、実験データに基づいて推定制御力を算出するマップ、もしくは、近似関数を作成して、これを制御力推定手段として用いて推定制御力を実制御力に近づけることで、車両の運動状態を推定している。この特許文献１に示される装置では、制御力推定手段は、静的システムである。
【特許文献１】特開２００１−２８７５２８号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
ところで、制御サスペンションは静的システムではなく、動的システムである。このため、制御力推定手段が静的システムである前記特許文献１に示される装置においては、実制御力は時間遅れを生じる。そして、このように時間遅れが生じる分、推定制御力は実制御力と異なり、その影響で、車両の運動状態の推定精度は低下し、所望の制御効果を発揮することが難しい。
制御効果を向上するためには、車両の運動状態の推定精度を向上する必要がある。そのためには、制御サスペンションの非線形特性だけでなく、動特性(時間遅れ)も考慮して実制御力を推定する必要がある。
【０００５】
本発明は、車両の運動状態の推定を精度高く行なうことができる状態推定装置、サスペンション制御装置及びサスペンションシステムを提供することを目的とする
【課題を解決するための手段】
【０００６】
本発明は、ばね上とばね下との間に制御力を発生する制御サスペンションを有して第１運動状態量を出力する車両に用いられ、前記車両の第２運動状態量を推定する状態推定装置であって、前記制御サスペンションの静特性に基づいて推定静的制御力を算出する制御力推定手段と、前記推定静的制御力を入力変数とするオブザーバとからなり、前記オブザーバは、前記推定静的制御力から推定動的制御力を算出する動特性付与手段と、前記推定動的制御力から、前記車両の推定第１運動状態量及び推定第２運動状態量を算出する車両モデルと、前記車両からの前記第１運動状態量と前記車両モデルからの前記推定第１運動状態量との差に相当する出力偏差から動特性補償信号を算出する動特性ゲイン手段と、前記出力偏差から車両モデル補償信号を算出する車両モデルゲイン手段とを有し、前記動特性補償信号を前記動特性付与手段に、前記車両モデル補償信号を車両モデルに、それぞれフィードバックすることを特徴とする。
【発明の効果】
【０００７】
本発明によれば、動的システムを有する車両の運動状態の推定を精度高く行なうことができる。
【発明を実施するための最良の形態】
【０００８】
以下、本発明の第１実施の形態に係るセミアクティブサスペンションシステム（サスペンションシステム）を図面に基づいて説明する。
図１〜図５において、本発明の第１実施の形態に係るセミアクティブサスペンションシステム１は、路面外乱及び実減衰力（実制御力）の入力を受ける車両（以下、実車両モデルともいう。）２に搭載して用いられている。セミアクティブサスペンションシステム１は、車両２の車体３（ばね上）と車輪４（ばね下）との間に介装されて減衰力（制御力）を発生する減衰力（制御力）可変のセミアクティブダンパ５（制御サスペンション）と、セミアクティブダンパ５に指令信号を送るサスペンション制御装置６と、から大略構成されている。本実施形態では、車両２は、サスペンション制御装置６の構成要素の一部を成し第１運動状態量検出手段の一例であるばね上加速度センサ７を有し、ばね上加速度センサ７が検出する当該車両２の上下方向の加速度を積分してばね上速度（車両２の第１運動状態量）を出力する。また、本実施形態では、セミアクティブダンパ５は、図５に示すように、その本体（ダンパ本体）８内部にピストン９が収納され、ピストン９に連結したピストンロッド１０が車体３に保持され、セミアクティブダンパ５が車輪４側に保持されている。
【０００９】
サスペンション制御装置６は、前記ピストン９の速度（以下、ピストン速度ともいう。車両の第２運動状態量）を推定して推定ピストン速度を得る状態推定装置１２（第２運動状態量推定手段）と、ばね上速度（第１運動状態量）及びピストン速度（第２運動状態量）から前記指令信号を算出する車両振動制御器１３と、から大略構成されている。
車両振動制御器１３は、フィードバック制御器１４及び指令信号算出手段１５からなっている。
【００１０】
状態推定装置１２は、セミアクティブダンパ５（制御サスペンション）の静特性に基づいて推定静的減衰力を算出する推定静的減衰力算出用変換マップ〔減衰力（制御力）推定手段〕１７を含むセミアクティブダンパモデル静的システム１８と、推定静的減衰力算出用変換マップ１７からの推定静的減衰力を入力変数とするオブザーバ２０と、から大略構成されている。
【００１１】
前記オブザーバ２０は、図１〜図４に示すように、前記推定静的減衰力から推定動的制御力を算出する実減衰力推定用オブザーバ２１と、前記推定動的制御力から、推定ばね上速度（車両２の推定第１運動状態量）及び推定ピストン速度（推定第２運動状態量）を算出する車両モデル（以下、車両近似モデルともいう。）２２を含むオブザーバ（以下、実車両モデル状態量推定用オブザーバという。）２３と、
前記ばね上加速度センサ７（車両２）からのばね上速度（第１運動状態量）と前記実車両モデル状態量推定用オブザーバ２３の車両近似モデル２２からの推定ばね上速度（推定観測出力とも言う。車両２の推定第１運動状態量）との差に相当する出力偏差から動特性補償信号Ｌ１を算出するオブザーバゲイン（以下、第１オブザーバゲインという。）２５と、前記出力偏差から車両モデル補償信号Ｌ２を算出するオブザーバゲイン（以下、第２オブザーバゲインという。）２６と、出力偏差の入力を受けて路面外乱を推定して推定外乱として出力する外乱推定手段２７と、を有している。
本実施形態では、第１オブザーバゲイン２５が動特性ゲイン手段を構成し、第２オブザーバゲイン２６が車両モデルゲイン手段を構成している。
動特性補償信号Ｌ１は、後述する動特性付与手段３０に入力され、その設定内容の調整に用いられる。車両モデル補償信号Ｌ２は、車両近似モデル２２にフィードバックしてその設定内容の調整に用いられるようになっている。
【００１２】
実車両モデルから出力される観測出力（ばね上速度ひいては減衰力）はフィードバック制御器１４に入力され、観測出力（ばね上速度ひいては減衰力）の入力を受けたフィードバック制御器１４は、目標減衰力を出力する。ここで、フィードバック制御器１４としては、スカイフック制御もしくはＨ∞制御等の制御を行える、振動制御に適した制御器を用いている。
指令信号算出手段１５は、フィードバック制御器１４からの目標減衰力と実車両モデル状態量推定用オブザーバ２３（車両近似モデル２２）により得られる推定ピストン速度（推定第２運動状態量）とから指令信号を算出し、これをセミアクティブダンパ５及び減衰力推定手段に出力する。
【００１３】
セミアクティブダンパ５は、実ピストン速度と指令信号とを入力としてこれらに依存した実減衰力を出力し、それを実車両モデルに伝える。
推定ピストン速度は、推定静的減衰力算出用変換マップ１７（制御力推定手段）が算出する推定静的減衰力と前記出力偏差とを基にして、オブザーバ２０（実車両モデル状態量推定用オブザーバ２３）により算出される。
【００１４】
実減衰力推定用オブザーバ２１は、図２に示すように、前記動特性付与手段３０と、実減衰力推定用オブザーバゲインＬｓを有する前記第１オブザーバゲイン２５と、からなっている。
動特性付与手段３０は、減衰力の動特性（時間遅れ）の実験データを基に一次遅れ系や二次遅れ系などで近似される動的システムとなっている。動特性付与手段３０では、出力偏差が生じる場合、実減衰力推定用オブザーバゲインＬｓ（第１オブザーバゲイン２５）により推定減衰力は修正される。
実車両モデル状態量推定用オブザーバ２３は、図１及び図３に示すように、車両近似モデル２２を含んで構成されており、実減衰力推定用オブザーバ２１で算出された推定動的減衰力、出力偏差及び推定外乱を基に、車両近似モデル２２を介して、推定観測出力と推定ピストン速度を算出する。
【００１５】
前記オブザーバ２０（実減衰力推定用オブザーバ２１、実車両モデル状態量推定用オブザーバ２３）は、本実施形態では、以下に述べるように設計されている。本実施形態では、図５に示す単輪モデルに対してオブザーバ２０を４輪分設計し、車両２の運動状態推定に用いる。
図５の単輪モデルから運動方程式を導出し、車両近似モデル２２の状態空間表現を次のように表す。
ここで、ｍ１はばね上質量、ｍ２はばね下質量、ｃ２はばね上質量とばね下質量の間の不変減衰係数、ｃｕは可変減衰係数、ｚ０は路面変位、ｚ１はばね下質量変位、ｚ２はばね上質量変位、ｋ１はタイヤのばね定数、ｋ２はサスペンションスプリングのばね定数である。
【００１６】
【数１】
2010058541.tif 000003
ここで、
【数２】
2010058541.tif 000004
は状態ベクトルで、Ａｍ，Ｂｍｕ，Ｂｍｚ，Ｃｍは
【数３】
2010058541.tif 000005
ｙｍは観測出力、Ｆｕは実減衰力、
【数４】
2010058541.tif 000006
は路面速度外乱である。
【００１７】
次に、セミアクティブダンパ５のモデル（セミアクティブダンパモデル）の動的システムを
【数５】
2010058541.tif 000007
で表す。ここで、ｆｕは減衰力の状態ベクトル、Ｆｓは静的減衰力である。
【００１８】
実施形態では、時定数Ｔの一次遅れ系で近似すると、式（３）のＡｓ、Ｂｓ、Ｃｓは、
Ａｓ＝−１／Ｔ、Ｂｓ＝１／Ｔ、Ｃｓ＝１    （４）
となり、実施形態では、ｆｕ＝Ｆｕとなる。ここで、式（１）と式（２）から、拡大系を構築すると
【数６】
2010058541.tif 000008
となる。ここで、
【数７】
2010058541.tif 000009
とすると、
【数８】
2010058541.tif 000010
となる。
【００１９】
式（６）より、この系のオブザーバ２０は、
【数９】
2010058541.tif 000011
として構成される。
ここで、Ｌはオブザーバゲインであり、実施形態では、定常カルマンフィルタによる方法を用いて求める。
ここで、オブザーバゲインＬは、
【数１０】
2010058541.tif 000012
である。
【数１１】
2010058541.tif 000013
は
【数１２】
2010058541.tif 000014
としたときの推定値である。今、誤差を
【数１３】
2010058541.tif 000015
とし、
【数１４】
2010058541.tif 000016
とすると、
【数１５】
2010058541.tif 000017
となる。まず、オブザーバゲインＬを定常カルマンフィルタによる方法を用いて求める。いま、外乱を無視した系〔式（４）〕がガウス性白色雑音ｗ（ｔ）、ｖ（ｔ）を受けるとすると、
【数１６】
2010058541.tif 000018
となる。また、ｗ（ｔ）、ｖ（ｔ）の平均値と共分散行列Ｑ、Ｒを
Ｅ〔ｗ（ｔ）〕＝０、
Ｅ〔ｗ（ｔ）ｗ（ｔ）Ｔ〕＝Ｑδ（ｔ−τ）        （１１）
Ｅ〔ｖ（ｔ）〕＝０、
Ｅ〔ｖ（ｔ）ｖ（ｔ）Ｔ〕＝Ｒδ（ｔ−τ）        （１２）
とする。ここで、Ｅ〔・〕は期待値であり、δはデルタ関数である。
オブザーバゲインＬはリカッチ方程式
ＡＰ+ＰＡＴ−ＰＣＴＲ−１ＣＰ＋Ｑ＝０          （１３）
の非負定対象行列解Ｐより、
Ｌ＝ＰＣＴＲ−１            （１４）
と決定される。ここで、ｙは車体３の上下速度で、加速度センサ７によって得られた値を積分して得る。
【００２０】
次に、ある小さな値αを用いて、
【数１７】
2010058541.tif 000019
とすると
【数１８】
2010058541.tif 000020
はＡ０＝Ａ−ＬＣとして、
【数１９】
2010058541.tif 000021
のように修正される。この修正により、ある行列Ｓが存在し、（ＳＣ，Ａ０，Ｂ´２）が強正実となるので、
【数２０】
2010058541.tif 000022
を満足する正定対称な行列Ｇ，Ｈが存在することになる。
【００２１】
ｅについて正定なリアプノフ関数を
【数２１】
2010058541.tif 000023
のように選ぶ。ここでΓは正の定数であり、
【数２２】
2010058541.tif 000024
である。
Ｖの時間微分をとると、
【数２３】
2010058541.tif 000025
となる。ここでＳは定数であり、αが十分小さいことを考慮すると、
【数２４】
2010058541.tif 000026
を
【数２５】
2010058541.tif 000027
のようにきめればよいが、
【数２６】
2010058541.tif 000028
の有界性を確保するためにσ修正法を用いて
【数２７】
2010058541.tif 000029
のように修正する。
以上のようにしてオブザーバ２０が設計されている。
【００２２】
本実施形態では、ばね上加速度センサ７（車両２）から検出した上下加速度を積分して得たばね上速度（第１運動状態量）と実車両モデル状態量推定用オブザーバ２３の車両近似モデル２２からの推定ばね上速度（車両の推定第１運動状態量）との差に相当する出力偏差から、第１オブザーバゲイン２５（動特性ゲイン手段）が動特性補償信号Ｌ１を算出し、また、第２オブザーバゲイン２６（車両モデルゲイン手段）が車両モデル補償信号Ｌ２を算出する。
そして、動特性補償信号Ｌ１は、動特性付与手段３０に入力され、動特性付与手段３０の設定内容の調整に用いられる。このため、制御上の時間遅れの発生が抑制され、良好な精度をもって振動抑制を行うことができる。さらに、車両モデル補償信号Ｌ２は、車両近似モデル２２にフィードバックしてその設定内容の調整に用いられ、出力偏差に応じた車両近似モデル２２の内容調整ができ、この分、振動抑制の精度向上を図ることができる。
【００２３】
また、本実施形態では、セミアクティブダンパ５の動特性(時間遅れ)モデルと車両モデルから成る拡大系に対してオブザーバ２０を設計しているので，車両の運動の推定状態量（推定運動状態量）だけでなく推定動的制御力も実値に近づくように修正できる。
推定動的制御力が実値に近づくことで、車両運動の推定状態量の推定精度は従来法に比べ高くなる。
さらに、車両２の運動の推定状態量の推定精度を向上できるので、ばね上加速度の応答を低減できる。
【００２４】
本実施形態では、上記のように設計したオブザーバ２０を用いて、実減衰力推定用オブザーバ２１を含む状態推定装置１２、サスペンション制御装置６ひいてはセミアクティブサスペンションシステム１が構成されている。
本発明は、実減衰力推定用オブザーバ２１が導入されているところに、その特徴があり、この特徴を有することにより、出力偏差が大きく生じても、推定減衰力の修正を適性に行い、ひいては良好な振動抑制効果を奏するようにしている。
【００２５】
本願発明者等は、上記実施形態に係る装置６（実減衰力推定用オブザーバ２１を備えている）が奏する上述した効果を確認するべく、実施形態に係る装置６に比して、図６に示すように、実減衰力推定用オブザーバ２１を備えていない装置（以下、比較対象装置５０という）を設け、実施形態に係る装置６と比較対象装置５０（実減衰力推定用オブザーバ２１を備えていない）とを対象にして、比較検証を行った。この比較検証では、０．５Ｈｚ〜２０Ｈｚの周波数成分を含むランダム波加振シミュレーションを行い、このシミュレーションを行ったときのばね上加速度のパワースペクトル密度（ＰＳＤ）及び推定ピストン速度と実ピストン速度との誤差のＰＳＤについて計測し、その結果、夫々図７及び図８に示す結果を得た。
【００２６】
比較対象装置５０では、実減衰力推定用オブザーバ２１を備えておらず、推定静的減衰力がそのまま推定減衰力として実車両モデル状態量推定用オブザーバ２３に入力され、出力偏差が大きく生じても、推定減衰力は、修正されないこととなり、振動抑制効果が劣るものになる。
これに対して、実施形態に係る装置６によれば、図７及び図８に示されるように、推定ピストン速度と実際のピストン速度との推定誤差を低減することで，ばね上加速度のパワースペクトル密度を低減でき、ひいては振動抑制を良好に行えることがわかる。
【００２７】
上記実施形態では、制御サスペンションとしてセミアクティブダンパ５を用いた場合を例にしたが、本発明はこれに限らず、アクティブサスペンションやロール剛性可変装置（アクティブスタビライザなど）など他の制御サスペンションを用い、セミアクティブダンパ５が発生する減衰力（制御力）の制御を行うのに代えて、アクティブサスペンションやロール剛性可変装置などが発生する言わばアクティブな力（制御力）の制御を行うように構成してもよい。
例えば、アクティブサスペンションの場合であれば、制御サスペンションであるアクティブサスペンションへの入力は指令信号であり、オブザーバからの出力は推定したい任意の運動状態量（各車輪のばね上速度など）となる。また、ロール剛性可変装置の場合であれば、制御サスペンションであるロール剛性可変装置への入力は指令信号であり、オブザーバ出力は推定したい任意の運動状態量（ロールの大きさなど）となる。
また、上記実施形態では、オブザーバ２０からの出力をばね上速度とピストン速度とし、ばね上速度を実測値と比較したが、車高センサなどを活用することによりピストン速度を実測できる場合は、ピストン速度を実測値と比較してもよい。また、オブザーバ２０からの出力についてもこれらの値に限られるものではなく、推定したい任意の運動状態量とすることができる。
【図面の簡単な説明】
【００２８】
【図１】本発明の一実施形態に係るサスペンションシステムを模式的に示すブロック図である。
【図２】図１の実減衰力推定用オブザーバを模式的に示すブロック図である。
【図３】図１の実車両モデル状態量推定用オブザーバを模式的に示すブロック図である。
【図４】図１の外乱推定手段を模式的に示すブロック図である。
【図５】車両の単輪モデルを示す図である。
【図６】比較対象装置を模式的に示すブロック図である。
【図７】実施形態及び比較対象装置による比較検証により得た、ばね上加速度のパワースペクトル密度（ＰＳＤ）を示す図である。
【図８】実施形態及び比較対象装置による比較検証により得た、ピストン速度推定誤差のパワースペクトル密度（ＰＳＤ）を示す図である。
【符号の説明】
【００２９】
１…セミアクティブサスペンションシステム（スペンションシステム）、５…セミアクティブダンパ（制御サスペンション）、６…サスペンション制御装置、１２…状態推定装置（第２運動状態量推定手段）、１７…推定静的減衰力算出用変換マップ〔制御力推定手段〕、２０…オブザーバ、２１…実減衰力推定用オブザーバ、２２…車両近似モデル（車両モデル）、２３…実車両モデル状態量推定用オブザーバ、２５…第１オブザーバゲイン（動特性ゲイン手段）、２６…第２オブザーバゲイン（車両モデルゲイン手段）、３０…動特性付与手段。
【図１】
2010058541.tif 000030
【図２】
2010058541.tif 000031
【図３】
2010058541.tif 000032
【図４】
2010058541.tif 000033
【図５】
2010058541.tif 000034
【図６】
2010058541.tif 000035
【図７】
2010058541.tif 000036
【図８】
2010058541.tif 000037