(57)【特許請求の範囲】
【請求項１】
走行車両のヨーレート、横加速度、及び車体速度、を測定し、
該測定結果と、前記横加速度の測定誤差を補償するために予め規定された遅れ要素Ｈ3とに基づき、次式
△βH＝Ｈ3（△θ）−Ｇｙ／Ｖｂ（但し、△θ：ヨーレート、Ｇｙ：横加速度、Ｖｂ：車体速度）を用いて車体スリップ角速度推定値△βHを算出し、
該車体スリップ角速度推定値△βHを積分演算することにより、車体スリップ角推定値βHを算出する車体スリップ角推定方法であって、
前記遅れ要素Ｈ3は、
前記車両における実ヨーレートとロール角との関係を次式の物理モデル
φ＝Ｈ1・Ｖｂ・Δψ  （但し、φ：ロール角、Ｖｂ：車体速度、Δψ：実ヨーレート、Ｈ1＝Ｋａ／（Ｊｓ2＋Ｄｓ＋Ｋ）、Ｋａ：定数、Ｊ：ロール慣性、Ｄ：ダンパー定数、Ｋ：ばね定数、ｓ：ラプラス変数）で表す場合の前記実ヨーレートの２次遅れ要素Ｈ1と、
前記横加速度と前記ロール角との関係を次式の物理モデル
φ＝Ｈ2・Ｇｙ  （但し、φ：ロール角、Ｇｙ：横加速度、Ｈ2＝Ｋｂ／（ｓ＋Ｋｃ）、Ｋｂ，Ｋｃ：定数、ｓ：ラプラス変数）で表す場合の前記横加速度の１次遅れ要素Ｈ2と、
を用いた次式Ｈ3＝Ｈ1／Ｈ2にて規定されることを特徴とする車体スリップ角推定方法。
【請求項２】
前記横加速度の測定を、前記車両のロールレートを測定した後、該測定結果に基づき、該ロールレートと前記横加速度との関係を記述した物理モデルを用いて前記横加速度を算出することにより行うことを特徴とする請求項１に記載の車体スリップ角推定方法。
【請求項３】
走行車両のヨーレートを測定するヨーレート測定手段と、
該車両の横加速度を測定する横加速度測定手段と、
該車両の車体速度を測定する車体速度測定手段と、
前記ヨーレート、横加速度、及び車体速度の測定結果と、前記横加速度の測定誤差を補償するために予め規定された遅れ要素Ｈ3とに基づき、次式
△βH＝Ｈ3（△θ）−Ｇｙ／Ｖｂ
（但し、△θ：ヨーレート、Ｇｙ：横加速度、Ｖｂ：車体速度）を用いて車体スリップ角速度推定値△βHを算出する車体スリップ角速度推定値算出手段と、
該車体スリップ角速度推定値算出手段にて算出された車体スリップ角速度推定値△βHを積分演算することにより、車体スリップ角推定値βHを算出する車体スリップ角推定値算出手段と、
を備えた車体スリップ角推定装置であって、
前記車体スリップ角速度推定値算出手段において、
前記遅れ要素Ｈ3は、
前記車両における実ヨーレートとロール角との関係を次式の物理モデル
φ＝Ｈ1・Ｖｂ・Δψ  （但し、φ：ロール角、Ｖｂ：車体速度、Δψ：実ヨーレート、Ｈ1＝Ｋａ／（Ｊｓ2＋Ｄｓ＋Ｋ）、Ｋａ：定数、Ｊ：ロール慣性、Ｄ：ダンパー定数、Ｋ：ばね定数、ｓ：ラプラス変数）で表す場合の前記実ヨーレートの２次遅れ要素Ｈ1と、
前記横加速度と前記ロール角との関係を次式の物理モデル
φ＝Ｈ2・Ｇｙ  （但し、φ：ロール角、Ｇｙ：横加速度、Ｈ2＝Ｋｂ／（ｓ＋Ｋｃ）、Ｋｂ，Ｋｃ：定数、ｓ：ラプラス変数）で表す場合の前記横加速度の１次遅れ要素Ｈ2と、
を用いた次式Ｈ3＝Ｈ1／Ｈ2にて規定されることを特徴とする車体スリップ角推定装置。
【請求項４】
前記横加速度測定手段は、
前記車両のロールレートを測定するロールレート測定手段と、
該ロールレート測定手段にて測定されたロールレートに基づき、該ロールレートと前記横加速度との関係を記述した物理モデルを用いて前記横加速度を算出する横加速度算出手段と、
を備えたことを特徴とする請求項３に記載の車体スリップ角推定装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、車両走行時における車体スリップ角を推定する車体スリップ角推定方法及び装置に関する。
【０００２】
【従来の技術】
従来より、車両走行時、特に車両旋回時の操縦安定性を高める技術として、車体スリップ角制御が知られている。
この車体スリップ角制御は、車両における車体前後軸方向と車両の実際の進行方向とが成す角度である車体スリップ角を推定し、この車体スリップ角の推定値に応じて、旋回外輪側の前輪にかかる制動力（ホイールシリンダ圧）を適宜増加させることにより、車体スリップ角を減少させる技術である。
【０００３】
そして、上記車体スリップ角制御においては、車体スリップ角の推定を高精度に行うことが必要である。従来、車体スリップ角βの推定は、車両に取り付けられたヨーレートセンサ等にて検出されるヨーレート△θと、横加速度センサ等にて検出される横加速度Ｇｙと、車輪速度センサ等を用いて検出される車体速度Ｖｂとに基づき、下記式（１）を用いて車体スリップ角速度△βの推定値△βpを算出し、さらに、これを積分演算することにより、行われていた。
【０００４】
△βp＝△θ−Ｇｙ／Ｖｂ    …（１）
【０００５】
【発明が解決しようとする課題】
しかし、上記式（１）に基づいて車体スリップ角速度△βの推定値△βpを算出するのでは、車両走行時に、車両がスピンした場合や、車両のロール角が大きくなった場合に、推定値△βpが、真値△βに対して大きく異なってしまい、その結果、車体スリップ角βを正確に推定することができないという問題があった。
【０００６】
図１０を用いて、上記問題が生じる理由について説明する。
まず、例えば、図１０（ａ）に示すように、旋回前の走行経路（ｘ軸方向に沿った経路）に対してヨー角ψの方向に進行している車両がスピンし、車体スリップ角がβとなった場合（換言すれば、図１０（ａ）中に一点鎖線で示される半径Ｒの円を旋回軌道としてヨーレートΔψで旋回している車両がスピンし、車体スリップ角がβとなった場合）、車両に取り付けられたヨーレートセンサでは、ヨー角ψに対するヨーレート△ψと、車体スリップ角βに対する車体スリップ角速度△βとに基づく下記式（２）で表される値を、ヨーレート△θとして検出する。
【０００７】
△θ＝△ψ＋△β    …（２）
尚、上記式（２）において、△ψは、車両がスピンしない状態で上記のように半径Ｒの円に沿って旋回した場合にヨーレートセンサで検出されるヨーレートであり、以下の説明では実ヨーレートとも記載する。
【０００８】
そして、上記のように車両がスピンした場合（換言すれば、車体スリップ角βの絶対値が０に比べて十分大きい場合）は、車両に取り付けられた横加速度センサにて検出される横加速度Ｇｙが、車両に実際に作用する横加速度Ｇｙ’（真値）と異なった値となる。
【０００９】
つまり、横加速度センサでは、車両に実際に作用する横加速度Ｇｙ’のうち、車体前後軸方向（図１０（ａ）中、「車体の向き」と記した方向）に対して直交する方向の成分、即ち、車両に実際に作用する横加速度Ｇｙ’の方向に対して車体スリップ角βだけ傾いた方向の成分、を横加速度Ｇｙとして検出するため、上記のように車両がスピンして、横加速度が車体前後軸方向に対して直交する方向に作用しなくなった状態では、横加速度を正確に検出できなくなる。
【００１０】
また、例えば、図１０（ｂ）に示すように、車両のロール角φが大きくなった場合も、横加速度センサにて検出される横加速度Ｇｙが、車両に実際に作用する横加速度Ｇｙ’（真値）と異なった値となる。
つまり、横加速度センサでは、車両に実際に作用する横加速度Ｇｙ’のうち、車体上下軸（図１０（ｂ）参照）に対して直交する方向の成分、即ち、車両に実際に作用する横加速度Ｇｙ’の方向に対してロール角φだけ下方を向いた成分を横加速度Ｇｙとして検出するため、上記のようにロール角φが大きくなって、横加速度が車体上下軸に対して直交する方向に作用しなくなった状態では、横加速度を正確に検出できなくなる。
【００１１】
従って、上記式（１）に基づいて車体スリップ角速度△βの推定値△βpを算出するのでは、車両走行時に、車両がスピンした場合や、車両のロール角φが大きくなった場合に横加速度を正確に検出できないため、推定値△βpが、真値△βに対して大きく異なってしまい、その結果、車体スリップ角βを高精度に推定することができなかった。
【００１２】
本発明は上記問題点を解決するためになされたものであり、その目的は、車両がスピンした場合や、車両のロール角が大きくなった場合であっても、車体スリップ角を正確に推定することができる車体スリップ角推定方法及び装置を提供することである。
【００１３】
【課題を解決するための手段及び発明の効果】
かかる目的を達成するためになされた請求項１に記載の車体スリップ角推定方法では、まず、走行車両のヨーレート△θ、横加速度Ｇｙ、及び車体速度Ｖｂを測定する。そして、横加速度Ｇｙの測定誤差を補償するために予め規定された遅れ要素Ｈ3を用いた次式△βH＝Ｈ3（△θ）−Ｇｙ／Ｖｂに上記測定されたヨーレート△θ、横加速度Ｇｙ、及び車体速度Ｖｂを適用して、車体スリップ角速度推定値△βHを算出し、さらに、車体スリップ角速度推定値△βHを積分演算することにより、車体スリップ角推定値βHを算出する。
【００１４】
従って、本発明方法（請求項１）によれば、横加速度Ｇｙの測定誤差が大きい場合、即ち、車両がスピンした場合や、車両のロール角が大きくなった場合であっても、車体スリップ角速度推定値△βHを、実際の車体スリップ角速度Δβ（真値）に近い値として算出することができ、その結果、車体スリップ角βを推定値βHとして高精度に推定することができる。
【００１５】
また、さらに、この車体スリップ角推定値βHに基づき、例えば、車両旋回時において、旋回外輪側前輪にかかる制動力（ホイールシリンダ圧）を適宜増加させ、車体スリップ角を減少させれば、車両走行時、特に車両旋回時における車両の操縦安定性を高めることができる。
【００１６】
そして、上記車体スリップ角速度推定値ΔβH算出用の式及び遅れ要素Ｈ3の具体的態様としては、車体スリップ角速度推定値△βH算出用の式を、ヨーレートと車両におけるロール角との関係をヨーレートの２次遅れ要素Ｈ1を用いて記述した物理モデル、及び横加速度とロール角との関係を横加速度の１次遅れ要素Ｈ2を用いて記述した物理モデルから導出した式とし、遅れ要素Ｈ3を、２次遅れ要素Ｈ1及び１次遅れ要素Ｈ2を用いた次式Ｈ3＝Ｈ1／Ｈ2にて規定される遅れ要素とするものである。
【００１７】
そして、本発明方法（請求項１）の態様における車体スリップ角速度推定値ΔβH算出用の式にて算出された車体スリップ角速度推定値ΔβHを積分演算することにより算出された車体スリップ角推定値βHは、後述実験例に示す実験の結果、車両がスピンした場合や、車両のロール角が大きくなった場合であっても、実際の車体スリップ角βと略一致しており、この推定値βHを、実際の車体スリップ角βを正確に推定した値として適用可能であることが判った。
【００１８】
つまり、本発明方法（請求項１）によれば、上記物理モデルから導出された車体スリップ角速度推定値ΔβH算出用の式で、ヨーレート△θの補正値であるＨ3（△θ）が、横加速度Ｇｙの測定誤差を加味した値として設定されているため、車両がスピンした場合や、車両のロール角が大きくなった場合であっても、車体スリップ角速度推定値△βHを真値△βに近い値として算出することができ、その結果、車体スリップ角βを推定値βHとして高精度に推定することができるのである。
【００１９】
尚、この態様における、ヨーレートとロール角との関係をヨーレートの２次遅れ要素Ｈ1を用いて記述した物理モデル、横加速度とロール角との関係を横加速度の１次遅れ要素Ｈ2を用いて記述した物理モデル、及びこれらの物理モデルと車体スリップ角速度推定値△βHの算出式「△βH＝Ｈ3（△θ）−Ｇｙ／Ｖｂ」との関係（導出手順）については、後述実施例において詳述する。
【００２０】
また、一方、車体スリップ角速度推定値△βH及び車体スリップ角推定値βHを算出する際に使用するヨーレート△θ、横加速度Ｇｙ、及び車体速度Ｖｂの具体的測定方法としては、ヨーレート△θを、例えば、車両に取り付けられたヨーレートセンサにて検出し、横加速度Ｇｙを、例えば、車両に取り付けられた横加速度センサにて検出し、車体速度Ｖｂを、例えば、左右の従動輪（フロントエンジン・リアドライブ方式の車両の場合は、前輪）の回転軸に夫々取り付けられた車輪速度センサより得られる左右の従動輪の回転速度の平均値（平均従動輪速度）として検出するものであっても良い。
【００２１】
また、車体速度Ｖｂの測定方法としては、上記の態様に限らず、後述実施例に示すように、まず、各車輪の回転軸に夫々取り付けられた車輪速度センサより得られる各車輪の回転速度を、他の方法にて測定されたヨーレート△θを用いて補正し、車両が減速状態にあるときは、補正された各車輪の回転速度のうち、現時点で最も車輪速度が大きい車輪についての値を車体速度Ｖｂ（推定値）として検出し、逆に、減速状態にないときは、補正された各車輪の回転速度のうち、現時点で最も車輪速度が小さい車輪についての値を車体速度Ｖｂ（推定値）として検出するものであっても良い。
【００２２】
また、横加速度Ｇｙの測定方法としては、請求項２に記載の車体スリップ角推定方法のように、まず、車両のロールレートを測定し、次いで、ロールレートと横加速度との関係を記述した物理モデルに上記測定されたロールレートを適用することにより、横加速度を算出するものであっても良い。
【００２３】
つまり、本発明方法（請求項２）では、車体スリップ角速度推定値△βH及び車体スリップ角推定値βHの算出に使用する横加速度Ｇｙを測定する際、横加速度Ｇｙを直接測定するのではなく、まず、車両のロールレートを測定し、この測定されたロールレートに基づいて横加速度Ｇｙ（具体的には、横加速度Ｇｙの推定値ＧｙH）を算出する。
【００２４】
従って、本発明方法（請求項２）によれば、横加速度センサを用いなくても、車体スリップ角速度推定値△βH及び車体スリップ角推定値βHの算出に使用する横加速度Ｇｙ（推定値ＧｙH）を測定することができる。尚、ロールレートと横加速度との関係を記述した物理モデルについては、後述実施例において詳述する。
【００２５】
次に、請求項３に記載の発明は、請求項１に記載の発明方法を実現するための構成を備えた車体スリップ角推定装置の発明であり、まず、走行車両のヨーレート△θをヨーレート測定手段にて測定し、横加速度Ｇｙを横加速度測定手段にて測定し、車体速度Ｖｂを車体速度測定手段にて測定する。次に、車体スリップ角速度推定値算出手段では、横加速度Ｇｙの測定誤差を補償するために予め設定された遅れ要素Ｈ3を用いた次式△βH＝Ｈ3（△θ）−Ｇｙ／Ｖｂに上記測定されたヨーレート△θ、横加速度Ｇｙ、及び車体速度Ｖｂを適用して、車体スリップ角速度推定値△βHを算出する。そして、車体スリップ角推定値算出手段では、車体スリップ角速度推定値算出手段にて算出された車体スリップ角速度推定値△βHを積分演算することにより、車体スリップ角推定値βHを算出する。
【００２６】
従って、本発明（請求項３）によれば、車体スリップ角速度推定値算出手段にて、横加速度Ｇｙの測定誤差が大きい場合、即ち、車両がスピンした場合や、車両のロール角が大きくなった場合であっても、車体スリップ角速度推定値△βHを、実際の車体スリップ角速度Δβ（真値）に近い値として算出することができる。そして、その結果、車体スリップ角推定値算出手段にて、実際の車体スリップ角βを正確に推定した値として、車体スリップ角推定値βHを算出することができる。
【００２７】
また、さらに、この車体スリップ角推定値βHに基づき、例えば、車両旋回時において、旋回外輪側前輪にかかる制動力（ホイールシリンダ圧）を適宜増加させ、車体スリップ角を減少させるよう構成すれば、車両走行時、特に車両旋回時における車両の操縦安定性を高めることができる。
【００２８】
そして、本発明（請求項３）は、車体スリップ角速度推定値算出手段にて使用される上記車体スリップ角速度推定値△βH算出用の式を、ヨーレートと車両におけるロール角との関係をヨーレートの２次遅れ要素Ｈ1を用いて記述した物理モデル、及び横加速度とロール角との関係を横加速度の１次遅れ要素Ｈ2を用いて記述した物理モデルから導出した式とし、遅れ要素Ｈ3を、２次遅れ要素Ｈ1及び１次遅れ要素Ｈ2を用いた次式Ｈ3＝Ｈ1／Ｈ2にて規定される遅れ要素とするよう構成したものである。
【００２９】
従って、本発明（請求項３）によれば、車体スリップ角速度推定値算出手段にて使用される、上記物理モデルから導出された車体スリップ角速度推定値ΔβH算出用の式で、ヨーレート△θの補正値であるＨ3（△θ）が、横加速度Ｇｙの測定誤差を加味した値として設定されているため、車両がスピンした場合や、車両のロール角が大きくなった場合であっても、車体スリップ角速度推定値△βHを真値△βに近い値として算出することができる。そして、その結果、車体スリップ角推定値算出手段にて、実際の車体スリップ角βを正確に推定した値として、車体スリップ角推定値βHを算出することができる。
【００３０】
次に、請求項４に記載の発明は、請求項２に記載の発明方法を実現するための構成を備えた車体スリップ角推定装置の発明であり、横加速度測定手段を、ロールレート測定手段にて、車両のロールレートを測定させ、横加速度算出手段にて、ロールレートと横加速度との関係を記述した物理モデルに上記測定されたロールレートを適用させて横加速度を算出させるよう構成したものである。
【００３１】
従って、本発明（請求項４）によれば、横加速度センサを用いなくても、車体スリップ角速度推定値△βH及び車体スリップ角推定値βHの算出に使用する横加速度Ｇｙ（具体的には、横加速度Ｇｙの推定値ＧｙH）を測定することができる。
【００３２】
【発明の実施の形態】
以下に、本発明の実施例を図面と共に説明する。
まず、図１は、本発明の一実施例としての車体スリップ角推定装置が適用された車体スリップ角制御装置の全体的構成を表す概略構成図である。尚、本実施例の車体スリップ角推定装置は、フロントエンジン・リアドライブ（ＦＲ）方式の車両に適用される。
【００３３】
図１に示す如く、この車両では、内燃機関２１から変速機２２を介して出力される車両の駆動力（駆動トルク）が、プロペラシャフト２３及びディファレンシャルギヤ２４を介して左右の後輪（駆動輪）（左後輪２５ＲＬ、右後輪２５ＲＲ）に分配される。
【００３４】
また、車両の各車輪（左前輪２５ＦＬ、右前輪２５ＦＲ、左後輪２５ＲＬ、右後輪２５ＲＲ）には、各車輪２５ＦＬ〜２５ＲＲに制動力を与える油圧式のブレーキ装置（以下、ホイールシリンダとも記す）４１ＦＬ、４１ＦＲ、４１ＲＬ、４１ＲＲが夫々設けられている。
【００３５】
そして、運転者によってブレーキベダル４２が踏み込まれると、マスタシリンダ４３からブレーキ油が圧送され、油圧回路４０を介して、各ホイールシリンダ４１ＦＬ〜４１ＲＲに加わるブレーキ油圧が増圧され、各車輪２５ＦＬ〜２５ＲＲに制動力が加えられるよう構成されている。
【００３６】
また、さらに、各車輪２５ＦＬ〜２５ＲＲには、各車輪２５ＦＬ〜２５ＲＲの回転速度（以下、車輪速度とも記す）を検出するための車輪速度センサ３１ＦＬ、３１ＦＲ、３１ＲＬ、３１ＲＲが夫々設けられている。
そして、車輪速度センサ３１ＦＬ〜３１ＲＲからの検出信号は、車両に設けられた他のセンサである、請求項３のヨーレート測定手段としてのヨーレートセンサ３２や、請求項３の横加速度測定手段としての横加速度センサ３３等からの検出信号と同様に、電子制御装置（ＥＣＵ）３０に入力される。
【００３７】
そして、ＥＣＵ３０は、これら各センサからの入力信号に基づき、運転者によるブレーキペダル４２操作とは別に、油圧回路４０内に設けられた各種アクチュエータ（図示はしない）を駆動して、各ホイールシリンダ４１ＦＬ〜４１ＲＲに加わるブレーキ油圧を調節することにより、各車輪２５ＦＬ〜２５ＲＲに加わる制動力を制御する。
【００３８】
即ち、ＥＣＵ３０は、車両走行（旋回）時に、上記各センサからの入力信号を用いて、車両における車体前後軸方向と車両の実際の進行方向とが成す角度である車体スリップ角を推定し、この車体スリップ角の推定値に応じて（具体的には、この車体スリップ角の推定値と、車体スリップ角の推定値に対する車体スリップ角速度の推定値とに応じて）、車体スリップ角を減少させるよう、左右の前輪２５ＦＬ、２５ＦＲのうち旋回外輪側の前輪に加わる制動力（ホイールシリンダ圧）を適宜増加させる車体スリップ角制御処理等を実行する。
【００３９】
尚、ＥＣＵ３０は、ＣＰＵ、ＲＯＭ、ＲＡＭ等を備えたマイクロコンピュータを中心に構成されており、該ＥＣＵ３０には、ブレーキペダル４２の踏込時にオン（ＯＮ）状態となり、ストップランプ（図示はしない）を点灯させるブレーキスイッチ４４からの検出信号も入力される。
【００４０】
次に、車両走行（旋回）時にＥＣＵ３０にて繰り返し実行される車体スリップ角制御処理について、図２〜図４に示すフローチャートに沿って説明する。
図２に示す如く、車体スリップ角制御処理が開始されると、まず、Ｓ１１０（Ｓはステップを表す）にて、上記各センサ及びブレーキスイッチ４４からの検出信号を読み込む。そして、続くＳ１２０では、車体速度推定処理を実行する。
【００４１】
この車体速度推定処理は、Ｓ１１０にて読み込まれた検出信号に基づき、車両の車体速度Ｖｂを推定（測定）する処理であり、図３に示す如く実行される。
即ち、車体速度推定処理では、まずＳ３１０にて、各車輪２５ＦＬ〜２５ＲＲの車輪速度の補正演算を行う。具体的には、車輪速度センサ３１ＦＬ〜３１ＲＲからの入力信号から検出された各車輪２５ＦＬ〜２５ＲＲの車輪速度ＶＷＦＬ、ＶＷＦＲ、ＶＷＲＬ、ＶＷＲＲを、ヨーレートセンサ３２からの入力信号より検出された車両のヨーレート△θを用いた下記式（３）〜（６）を用いて、車両の重心位置の速度に換算する補正演算を行う。
【００４２】
ＶＷＦＬforＶｂ＝ＶＷＦＬ−Ｌｆ×△θ    …（３）
ＶＷＦＲforＶｂ＝ＶＷＦＲ−Ｌｆ×△θ    …（４）
ＶＷＲＬforＶｂ＝ＶＷＲＬ−Ｌｒ×△θ    …（５）
ＶＷＲＲforＶｂ＝ＶＷＲＲ−Ｌｒ×△θ    …（６）
尚、上記式（３）〜（６）におけるパラメータのうち、ＶＷＦＬforＶｂ、ＶＷＦＲforＶｂ、ＶＷＲＬforＶｂ、ＶＷＲＲforＶｂは、夫々各車輪２５ＦＬ〜２５ＲＲの補正後の車輪速度を示す。また、Ｌｆは、車両の重心からフロント軸（フロントドライブシャフト）までの最短距離を示し、Ｌｒは、車両の重心からリア軸（リアドライブシャフト）までの最短距離を示す。また、ヨーレートセンサ３２からの検出値であるヨーレート△θは、左旋回（自転）時に正、右旋回（自転）時に負の値として出力されるものである。
【００４３】
次に、Ｓ３２０では、現在車両が減速状態にあるか否かを判断する。この判断は、ブレーキスイッチ４４からの入力信号がオン（ＯＮ）状態であるか否かを判断することにより行われる。
そして、Ｓ３２０にて車両が減速状態にあると判断された場合は、Ｓ３３０に移行し、下記式（７）にて車体速度Ｖｂ（推定値）を算出する。
【００４４】
Ｖｂ＝ｍａｘ（ＶＷ＊＊forＶｂ）    …（７）
尚、上記式（７）におけるパラメータのうち、ＶＷ＊＊forＶｂは、上記式（３）〜（６）にて算出された全車輪の補正後の車輪速度を示し、＊＊は、各車輪２５ＦＬ〜２５ＲＲを示す。そして、上記式（７）は、全車輪の補正後の車輪速度のうち、最大車輪速度であるものを車体速度Ｖｂとして算出することを示す。
【００４５】
つまり、車両が減速状態であれば、例えば、ブレーキの各車輪２５ＦＬ〜２５ＲＲへの作用、あるいはエンジンブレーキの各車輪２５ＦＬ〜２５ＲＲへの作用により、車輪が減速スリップに陥ることがあり、減速スリップに陥った車輪の車輪速度（補正後の車輪速度）が車体速度より極端に小さくなる可能性がある。
【００４６】
そこで、Ｓ３３０では、全車輪の補正後の車輪速度のうち、減速スリップに陥っていない車輪（換言すれば、路面にグリップした車輪）の補正後の車輪速度に対応する最大車輪速度を車体速度Ｖｂとして算出するのである。
一方、Ｓ３２０にて車両が減速状態にないと判断された場合は、Ｓ３４０に移行し、下記式（８）にて車体速度Ｖｂ（推定値）を算出する。
【００４７】
Ｖｂ＝ｍｉｎ（ＶＷ＊＊forＶｂ）    …（８）
上記式（８）は、全車輪の補正後の車輪速度のうち、最小車輪速度であるものを車体速度Ｖｂとして算出することを示す。
つまり、車両が減速状態にない場合としては、例えば、車両の加速中が考えられ、この際には加速スリップに陥った車輪の車輪速度（補正後の車輪速度）が車体速度より極端に大きくなる可能性がある。
【００４８】
そこで、Ｓ３４０では、全車輪の補正後の車輪速度のうち、加速スリップに陥っていない車輪（換言すれば、路面にグリップした車輪）の補正後の車輪速度に対応する最小車輪速度を車体速度Ｖｂとして算出するのである。
そして、Ｓ３３０、或いはＳ３４０にて車体速度Ｖｂが算出されると、今度は、Ｓ３５０に移行する。
【００４９】
Ｓ３５０では、上記のように算出された車体速度Ｖｂの変化勾配の制限を行い、車体速度推定処理を完了する。
即ち、Ｓ３５０では、前回フロー時のＳ１３０の処理（後述）にて推定値として算出された車体の前後軸方向の加速度である車体前後加速度△Ｖｂに応じて、前回フロー時に算出された車体速度から今回フロー時に算出された車体速度への変化量を制限することにより、今回フローにて算出された車体速度Ｖｂを補正する。
【００５０】
こうして、車体速度推定処理（Ｓ１２０）により、車体速度Ｖｂが推定（測定）されると、今度は、Ｓ１３０（図２参照）にて、車体前後加速度推定処理を実行する。
この車体前後加速度推定処理は、今回フローにおける車体前後加速度△Ｖｂ（推定値）を算出する処理であり、図４に示す如く実行される。
【００５１】
即ち、車体前後加速度推定処理では、まずＳ４１０にて、今回フローまでの所定時間内に車体速度推定処理（Ｓ１２０）にて算出された車体速度Ｖｂを、ローパスフィルタにてフィルタリングする。具体的には、例えば、１０Ｈｚ以下の周波数のみを通すローパスフィルタを用いて、推定値としての車体速度Ｖｂ中のノイズを除去する。
【００５２】
次に、Ｓ４２０では、Ｓ４１０にてフィルタリングされた車体速度Ｖｂを微分処理する。
そして、Ｓ４３０では、Ｓ４２０にて算出された値、即ち、車体前後加速度△Ｖｂをフィルタリングして、車体前後加速度推定処理を完了する。具体的には、例えば、Ｓ４２０にて算出された車体前後加速度△Ｖｂを、２Ｈｚ以下の周波数のみを通すローパスフィルタにかける。
【００５３】
こうして、車体前後加速度推定処理（Ｓ１３０）により、今回フローにおける車体前後加速度△Ｖｂ（推定値）が算出されると、今度は、Ｓ１４０（図２参照）にて、ヨーレートセンサ３２からの入力信号より検出された車両のヨーレート△θと、横加速度センサ３３からの入力信号より検出された車両の横加速度Ｇｙと、車体速度推定処理（Ｓ１２０）により推定（測定）された車体速度Ｖｂと、後述の式（１３）及び式（１４）にて定義されるＨ1及びＨ2にてＨ3＝Ｈ1／Ｈ2と表されるＨ3と、をパラメータとする次式（９）にて、車体スリップ角速度推定値△βHを算出する。
【００５４】
△βH＝Ｈ3（△θ）−Ｇｙ／Ｖｂ    …（９）
尚、横加速度センサ３３からの検出値である横加速度Ｇｙは、左旋回（公転）時に正、右旋回（公転）時に負の値として出力されるものである。また、車両が極低速で走行している場合（例えば、５km/h以下のような車輪速度センサ３１ＦＬ〜３１ＲＲの分解能の限界以下の極低速で走行している場合）には、上記式（９）においてＶｂを０とした演算を行わないようにするため、Ｖｂに固定値（例えば、５km/h）を代入して車体スリップ角速度推定値△βHを算出する。
【００５５】
そして、続くＳ１５０では、Ｓ１４０にて算出された車体スリップ角速度推定値△βHを積分演算することにより、車体スリップ角推定値βHを算出する。車体スリップ角速度推定値△βHの積分演算は、誤差の蓄積によるオーバーフローを防ぐため、下記式（１０）にて行う。
【００５６】
βH(n)＝（βH(n-1)＋△βH）・ｋ１    …（１０）
尚、上記式（１０）において、ｋ１は定数であり、例えば、ｋ１＝０．９９と設定される。また、ｎは演算回数を示し、βH(n)は、今回フローにて算出された車体スリップ角推定値βHであり、βH(n-1)は、前回フローにて算出された車体スリップ角推定値βHである。また、このように算出される車体スリップ角推定値βHは、車両における車体前後軸方向が車両の実際の進行方向に対して左方向に傾いた（自転した）場合に正、右方向に傾いた（自転した）場合に負の値として出力されるものである。
【００５７】
ここで、Ｓ１４０の処理にて車体スリップ角速度推定値△βHを求める際に使用する上記式（９）は、以下のように設定されている。
即ち、まず、実ヨーレート△ψ（「発明が解決しようとする課題」の項で述べた△ψについての説明参照）と車両におけるロール角φとの関係を記述した物理モデル（ラプラス変換されたモデル）と、横加速度Ｇｙとロール角φとの関係を記述した物理モデル（ラプラス変換されたモデル）とが、夫々、下記式（１１）と式（１２）とのように表される。
【００５８】
【数１】
0004244453.tif 000002
【００５９】
（但し、Ｋａ：定数、Ｊ：ロール慣性、Ｄ：ダンパー定数、Ｋ：ばね定数）
【００６０】
【数２】
0004244453.tif 000003
【００６１】
（但し、Ｋｂ，Ｋｃ：定数）
そして、ヨーレートセンサによる検出値△θと実ヨーレート△ψと車体スリップ角速度△βとの関係を表す式（２）に基づき、上記式（１１）に△ψ＝△θ−△βを代入すれば、検出値であるヨーレート△θとロール角φとの関係式を得ることができ、さらに、上記式（１１）、（１２）に対して、
【００６２】
【数３】
0004244453.tif 000004
【００６３】
【数４】
0004244453.tif 000005
【００６４】
とおき、ロール角φに対するヨーレート△θ（実ヨーレート△ψ）の２次遅れ要素Ｈ1と、ロール角φに対する横加速度Ｇｙの１次遅れ要素Ｈ2とを定義すれば、上記式（１１）と式（１２）とが夫々下記式（１５）と式（１６）とのように変形される。
【００６５】
φ＝Ｈ1（Ｖｂ・（△θ−△β））    …（１５）
φ＝Ｈ2（Ｇｙ）                    …（１６）
そして、上記式（１５）の右辺と式（１６）の右辺とを等号で結び、Ｈ3＝Ｈ1／Ｈ2とすれば、下記式（１７）が得られ、
Ｈ3（△β）＝Ｈ3（△θ）−Ｇｙ／Ｖｂ    …（１７）
さらに、上記式（１７）におけるＨ3（△β）を車体スリップ角速度推定値△βHとすることで上記式（９）が得られる。
【００６６】
尚、上記式（９）にて算出される車体スリップ角速度推定値△βHは、上記式（１７）での表記から、車体スリップ角速度△βではなく、Ｈ3（△β）と表されるものであるが、後述実験例に示す如く、この推定値△βHは、実用上、実際の車体スリップ角速度△βを正確に推定した値として適用可能なものである。
【００６７】
つまり、上記式（９）は、従来、車体スリップ角△βを推定する際に用いられた式（１）においてヨーレート△θをＨ3にて補正した式であり、この補正値Ｈ3（△θ）が、横加速度センサ３３にて検出される横加速度Ｇｙの測定誤差を加味した値として設定されているため、車両がスピンした場合や、車両のロール角φが大きくなった場合であっても、上記式（９）にて、車体スリップ角速度推定値△βHを真値△βに近い値として算出することができるのである。
【００６８】
そして、その結果、この車体スリップ角速度推定値△βHを用いたＳ１５０の処理にて、車両がスピンした場合や、車両のロール角φが大きくなった場合であっても、車体スリップ角推定値βHを、実際の車体スリップ角βと略一致した推定値として算出することができる。つまり、この推定値βHを、実際の車体スリップ角βを正確に推定した値として適用することができるのである。
【００６９】
そして、上記のようにして、Ｓ１４０の処理にて車体スリップ角速度推定値△βHが算出され、Ｓ１５０の処理にて車体スリップ角推定値βHが算出されると、今度は、Ｓ１６０に移行する。
Ｓ１６０では、上記車体スリップ角速度推定値△βH（今回フローにおける算出値）と、車体スリップ角推定値βH（今回フローにおける算出値）とを用いた下記式（１８）にて、車体スリップ角βのＰＩ制御用パラメータｗβを算出する。
【００７０】
ｗβ＝ｋ２・βH＋ｋ３・△βH   …（１８）
（但し、ｋ２，ｋ３：定数）
次に、続くＳ１７０では、ＰＩ制御用パラメータｗβの絶対値が評価係数Ｋｄより大きいか否かの判断を行う。
【００７１】
そして、Ｓ１７０にてＰＩ制御用パラメータｗβの絶対値が評価係数Ｋｄより大きいと判断された場合、つまり、車両の走行（旋回）状態が不安定である（例えば、車両がスピンする状態である）と判断された場合は、Ｓ１８０に移行し、今度は、ＰＩ制御用パラメータｗβが０より大きいか否かの判断を行う。
【００７２】
そして、Ｓ１８０にてＰＩ制御用パラメータｗβが０より大きい、即ち、正であると判断された場合は、Ｓ１９０に移行して、旋回外輪側の前輪を右前輪２５ＦＲであると判断し、油圧回路４０内の各種アクチュエータを駆動して、右前輪２５ＦＲに加わる制動力、つまり、ホイールシリンダ４１ＦＲに加わるブレーキ油圧を適宜増加させて車体スリップ角βを減少させ、当該車体スリップ角制御処理を終了する。
【００７３】
また、Ｓ１８０にてＰＩ制御用パラメータｗβが０より小さい、即ち、負であると判断された場合は、Ｓ２００に移行して、旋回外輪側の前輪を左前輪２５ＦＬであると判断し、油圧回路４０内の各種アクチュエータを駆動して、左前輪２５ＦＬに加わる制動力、つまり、ホイールシリンダ４１ＦＬに加わるブレーキ油圧を適宜増加させて車体スリップ角βを減少させ、当該車体スリップ角制御処理を終了する。
【００７４】
また、一方、Ｓ１７０にてＰＩ制御用パラメータｗβの絶対値が評価係数Ｋｄより大きくはないと判断された場合は、車両の走行（旋回）状態が安定した状態となっており、左右の前輪２５ＦＬ、２５ＦＲのうち一方の制動力を増加させる必要がないと判断して、当該車体スリップ角制御処理を終了する。
【００７５】
尚、Ｓ１２０の処理は、請求項３の車体速度測定手段に相当し、Ｓ１４０の処理は、請求項３の車体スリップ角速度推定値算出手段に相当し、Ｓ１５０の処理は、請求項３の車体スリップ角推定値算出手段に相当する。
以上説明したように、本実施例では、Ｓ１１０及びＳ１２０にて測定されたヨーレート△θ、横加速度Ｇｙ、及び車体速度Ｖｂに基づき、式（９）を用いて車体スリップ角速度推定値△βHを算出し（Ｓ１４０）、さらに、この車体スリップ角速度推定値△βHを式（１０）を用いて積分演算する（Ｓ１５０）ことにより、車体スリップ角推定値βHを算出する。
【００７６】
そして、本実施例では、このように算出された車体スリップ角推定値βHが、車両がスピンした場合や、車両のロール角φが大きくなった場合であっても、実際の車体スリップ角βと略一致しており、この推定値βHを、実際の車体スリップ角βを正確に推定した値として適用することができる。
【００７７】
つまり、本実施例では、式（９）において、ヨーレート△θの補正値であるＨ3（△θ）が、横加速度Ｇｙの測定誤差を加味した値として設定されているため、車両がスピンした場合や、車両のロール角φが大きくなった場合であっても、車体スリップ角速度推定値△βHを真値△βに近い値として算出する（Ｓ１４０）ことができ、その結果、車体スリップ角βを正確に推定した値として、車体スリップ角推定値βHを算出する（Ｓ１５０）ことができるのである。
【００７８】
また、本実施例では、このように算出された車体スリップ角推定値βH（具体的には、車体スリップ角推定値βH及び車体スリップ角速度推定値△βH）に対する評価を行い（Ｓ１７０〜Ｓ１８０）、車両の走行（旋回）状態が不安定である（例えば、車両がスピンする状態である）と判断された場合は、左右の前輪２５ＦＬ、２５ＦＲのうち一方の制動力を適宜増加させ（Ｓ１９０、Ｓ２００）、車体スリップ角βを減少させる。従って、車両走行時、特に車両旋回時における車両の操縦安定性を高めることができる。
【００７９】
以上、本発明の一実施例について説明したが、本発明は、上記実施例に限定されるものではなく、種々の態様を採ることができる。
例えば、上記実施例では、本発明をフロントエンジン・リアドライブ（ＦＲ）方式の車両（換言すれば、後輪駆動車）に適用した場合について説明したが、本発明は、フロントエンジン・フロントドライブ（ＦＦ）方式の車両（換言すれば、前輪駆動車）に適用しても、上記実施例と同様の効果を得ることができる。
【００８０】
また、上記実施例では、車体スリップ角速度推定値△βH及び車体スリップ角推定値βHの算出（Ｓ１４０、Ｓ１５０）に使用する横加速度Ｇｙを測定する（Ｓ１１０）際、横加速度センサ３３にて横加速度Ｇｙを直接測定したが、横加速度センサ３３の代わりにロールレートセンサを設けて、まず、該ロールレートセンサにて、車両のロールレート△φを測定させ、次いで、ＥＣＵ３０にて、該ロールレートセンサからの入力信号より検出されたロールレート△φを、下記式（１９）にて表される、ロールレートと横加速度との関係を記述した物理モデルに適用して、横加速度Ｇｙに対する推定値ＧｙHを算出しても良い。
【００８１】
【数５】
0004244453.tif 000006
【００８２】
（但し、Ｋｅ，Ｋｆ：定数）
そして、このようにすれば、横加速度センサ３３を用いなくても、車体スリップ角速度推定値△βH及び車体スリップ角推定値βHの算出に使用する横加速度Ｇｙ（具体的には、横加速度Ｇｙの推定値ＧｙH）を測定することができる。
【００８３】
尚、この場合のロールレートセンサは請求項４のロールレート測定手段に相当し、上記式（１９）を用いてＥＣＵ３０にて行われる横加速度Ｇｙ（推定値ＧｙH）の算出処理は請求項４の横加速度算出手段に相当する。
また、一方、上記実施例では、車体速度推定処理（Ｓ１２０）において、車輪速度ＶＷＦＬ〜ＶＷＲＲを車両の重心位置の速度に換算する（Ｓ３１０）際に、車両のヨーレート△θを用いた式（３）〜（６）を用いたが、式（３）〜（６）におけるヨーレート△θを、横加速度センサ３３からの入力信号より検出される横加速度Ｇｙを用いた下記式（２０）にて換算値△θHとして算出しても良い。
【００８４】
△θH＝Ｇｙ／Ｖｂ(n-1)   …（２０）
尚、上記式（２０）において、Ｖｂ(n-1)は、前回フローにて算出された車体速度Ｖｂである。また、上記式（２０）を使用する態様では、車両が極低速で走行している場合（例えば、５km/h以下のような車輪速度センサ３１ＦＬ〜３１ＲＲの分解能の限界以下の極低速で走行している場合）には、上記式（２０）においてＶｂ(n-1)を０とした演算を行わないようにするため、Ｖｂ(n-1)に固定値（例えば、５km/h）を代入して△θHを算出すれば良い。
【００８５】
また、上記実施例では、車体速度推定処理（Ｓ１２０）において、車両が減速状態にあるか否かを判断する（Ｓ３２０）際に、ブレーキスイッチ４４からの入力信号を用いて判断を行ったが、例えば、下記３つの態様（▲１▼〜▲３▼）のうち１つを採用しても良い。また、下記３つの態様（▲１▼〜▲３▼）及び上記実施例の態様（即ち、ブレーキスイッチ４４からの入力信号を用いる態様）のうち少なくとも２つを組み合わせた態様を採用しても良い。
【００８６】
▲１▼  車体前後加速度推定処理（Ｓ１３０）にて算出される車体前後加速度△Ｖｂ（前回フローにて算出された値）が正か負かを判断することによって車両が減速状態にあるか否かを判断する。
▲２▼  ＥＣＵ３０にてＩＤＬ（アイドル）信号（例えば、アクセルペダル（図示はしない）が踏み込まれているか否かによって出力が変化する信号）を検出することによって車両が減速状態にあるか否かを判断する。
【００８７】
▲３▼  マスタシリンダ４３内の圧力を検出するマスタシリンダ圧センサ（図示はしない）を設け、該マスタシリンダ圧センサによる検出値（圧力）が基準値以上であって、実質的に車輪制動力を発生している状態か否かを判断することによって車両が減速状態にあるか否かを判断する。
【００８８】
また、一方、車体速度Ｖｂの測定（推定）処理（請求項３の車体速度測定手段に相当）としては、左右の従動輪（フロントエンジン・リアドライブ方式の車両の場合は、前輪２５ＦＬ、２５ＦＲ）の車輪速度を検出する車輪速度センサより得られる左右の従動輪の回転速度の平均値（平均従動輪速度）を車体速度Ｖｂ（推定値）として検出するものであっても良い。
【００８９】
次に、上述した上記実施例の有する効果を裏付ける各種実験例について説明する。
［実験例１］
この実験では、上記実施例と同様の車体スリップ角推定装置が備えられ、且つ、操舵角センサ（図示はしない）が取り付けられた実験車両を走行させ、スピン発生を伴うレーンチェンジを行った場合（速度５５km/h）と、通常旋回（換言すれば、車両がスピンしない旋回）を行った場合（速度４０km/h）とにおいて、式（９）にて算出される車体スリップ角速度推定値△βHが、実際の車体スリップ角速度△βを正確に推定した値として適用可能であるか否かを検証した。
【００９０】
この実験の結果について図５及び図６を用いて説明する。
まず、図５は、上記実験車両がスピン発生を伴うレーンチェンジを行った場合（速度５５km/h）における各検出値（算出値）、即ち、△θ、Ｈ3（△θ）、Ｇｙ／Ｖｂ、△βH、及び操舵角センサ検出値（図中に「ステアリング」として示した検出値）の時間経過に対する変化を示したものである。尚、操舵角センサからの検出値は、左操舵時に正、右操舵時に負の値として出力されるものである。
【００９１】
そして、図５から、式（９）にて算出される車体スリップ角速度推定値△βHが、式（１）にて算出される従来の推定値△βpと異なり、実際の車体スリップ角速度△βに近い値として検出（算出）されることが判った（βとβHの定量的な比較は実験例２で後述）。
【００９２】
つまり、まず、ヨーレート△θ（Ｈ3によって補正されていない、ヨーレートセンサ３２による検出値）は、この実験車両がスピンしたと考えられる時点ｔ１より以前（換言すれば、時点ｔ１以前で、実際の車体スリップ角速度△βが０に近い値になると考えられる時点）の、ステアリングが操舵された時点（換言すれば、時点ｔ１以前で、操舵角センサ検出値が負となった時点）の直後から、Ｇｙ／Ｖｂと離れはじめ、その結果、式（１）にて算出される車体スリップ角速度△βの推定値△βp（図示はしない）が、時点ｔ１より以前において、０と異なった値となり、さらに、時点ｔ１以降においても、実験車両がいつスピンし始めたのかを正確に検出することもできないことが判った。従って、式（１）では、真値△βに近い値を算出することができないことが判った。
【００９３】
一方、ヨーレート△θの補正値Ｈ3（△θ）とＧｙ／Ｖｂとは、時点ｔ１より以前において、略等しい値となり、その結果、式（９）の算出値である車体スリップ角速度推定値△βHが０に近い値となり、この推定値△βHが、真値△βに近い値となることが判った。
【００９４】
そして、この実験車両がスピンすると（時点ｔ１）、今度は、Ｈ3（△θ）とＧｙ／Ｖｂとが離れはじめ、対応する車体スリップ角速度推定値△βHが、０と異なる値に変化し始める。つまり、車体スリップ角速度推定値△βHが、真値△βと同様の変化を示すことが判った。
【００９５】
即ち、この実験車両では、時点ｔ１より以前の、ステアリングが操舵された時点の直後からロール角φが大きくなると考えられ、さらに、時点ｔ１以降ではスピンも発生するため、横加速度センサ３３にて検出される横加速度Ｇｙが、実験車両に実際に作用する横加速度と異なった値となり、式（１）では、実際の車体スリップ角速度△βに近い値を検出（算出）することができないが、式（９）によれば、このような場合であっても、実際の車体スリップ角速度△βに近い値を推定値△βHとして検出（算出）可能であることが判った。
【００９６】
また、図６は、上記実験車両が通常旋回（スピンを伴わない旋回）を行った場合（速度４０km/h）における上記各検出値（算出値）の時間経過に対する変化を示したものである。
そして、図６においても、式（９）にて算出される車体スリップ角速度推定値△βHが、式（１）にて算出される従来の推定値△βpと異なり、実際の車体スリップ角速度△βに近い値として検出（算出）されることが判った。
【００９７】
つまり、この場合は実際の車体スリップ角速度△βが常に０に近い値になると考えられるが、図６では、ヨーレート△θが、旋回時においてＧｙ／Ｖｂと異なった値となり、その結果、式（１）にて算出される従来の推定値△βpが旋回時に０に近い値とならず、式（１）では真値△βに近い値を算出することができないことが判った。
【００９８】
一方、Ｈ3（△θ）とＧｙ／Ｖｂとは、常に略等しい値となり、その結果、車体スリップ角速度推定値△βHが常に０に近い値となり、この推定値△βHが、真値△βに近い値となることが判った。
即ち、この実験車両では、上記通常旋回時においてロール角φが大きくなると考えられ、横加速度センサ３３にて検出される横加速度Ｇｙが、実験車両に実際に作用する横加速度と異なった値となり、式（１）では、実際の車体スリップ角速度△βに近い値を検出（算出）することができないが、式（９）によれば、このような場合であっても、実際の車体スリップ角速度△βに近い値を推定値△βHとして検出（算出）可能であることが判った。
【００９９】
従って、図５及び図６に示した実験結果から、ヨーレート△θの補正値Ｈ3（△θ）が、横加速度センサ３３にて検出される横加速度Ｇｙの測定誤差を加味した値として設定されており、車両がスピンした場合や、車両のロール角φが大きくなった場合であっても、式（９）にて、車体スリップ角速度推定値△βHを真値△βに近い値として算出することができることを確認できた。
【０１００】
つまり、この実験例１の実験結果から、車両がスピンした場合や、車両のロール角φが大きくなった場合であっても、式（９）にて算出される車体スリップ角速度推定値△βHが、実際の車体スリップ角速度△βを正確に推定した値として適用可能なものであることを確認できた。
［実験例２］
この実験では、上記実施例と同様の車体スリップ角推定装置が備えられ、且つ、周知の対地速度計（図示はしない）が車両後部に取り付けられた実験車両を走行させ、式（９）にて算出された車体スリップ角速度推定値△βHを積分演算して得られた車体スリップ角推定値βHが、実際の車体スリップ角βを正確に推定した値として適用可能であるか否かを検証した。
【０１０１】
つまり、この実験では、上記実験車両において、上記実施例における車体スリップ角推定値βHの他に、式（１）にて算出される従来の車体スリップ角速度推定値△βpを積分演算して得られる従来の車体スリップ角推定値βpを検出すると共に、対地速度計による検出値から算出される実際の車体スリップ角βも検出し、これら３つの値を比較することにより、上記検証を行った。
【０１０２】
尚、対地速度計では、実験車両における車体前後軸方向の進行速度Ｖx、及び車体左右軸方向（車体前後軸方向に対して垂直で、例えば車軸（フロントドライブシャフト、リアドライブシャフト）に平行な方向）の進行速度Ｖyを検出した。そして、この実験では、下記式（２１）にて実際の車体スリップ角β（検出値）を算出した。
【０１０３】
β＝ｔａｎ-1（Ｖy／Ｖx）    …（２１）
この実験の結果について図７〜図９を用いて説明する。
図７は、上記実験車両にアスファルト上を走行させた場合について、図８は、上記実験車両に圧雪上を走行させた場合について、図９は、上記実験車両に氷上を走行させた場合について、上記実施例における車体スリップ角推定値βH（図７〜図９に「βH（補正値）」として示した検出値）、従来の車体スリップ角推定値βp（図７〜図９に「無補正」として示した検出値）、及び実際の車体スリップ角βの検出値（図７〜図９に「対地速度計」として示した検出値）の時間経過に対する変化を示したものである。
【０１０４】
そして、図７〜図９いずれの場合であっても、車体スリップ角βが０に比べて十分大きくなるような走行を行った際（つまり、実験車両がスピンしたり、実験車両のロール角φが大きくなるような走行を行った際）に、上記実施例における車体スリップ角推定値βHが、従来の車体スリップ角推定値βpに比べ、実際の車体スリップ角β（検出値）に近い値となることが判った。
【０１０５】
従って、この実験例２の実験結果から、車両がスピンした場合や、車両のロール角φが大きくなった場合であっても、上記実施例における車体スリップ角推定値βHが、実際の車体スリップ角βを正確に推定した値として適用可能なものであることを確認できた。
【図面の簡単な説明】
【図１】  実施例の車体スリップ角推定装置が適用された車体スリップ角制御装置の全体的構成を表す概略構成図である。
【図２】  電子制御装置（ＥＣＵ）において実行される車体スリップ角制御処理を表すフローチャートである。
【図３】  図２のＳ１２０にて実行される車体速度推定処理を表すフローチャートである。
【図４】  図２のＳ１３０にて実行される車体前後加速度推定処理を表すフローチャートである。
【図５】  実験例１による検出結果を表すグラフである。
【図６】  実験例１による検出結果を表すグラフである。
【図７】  実験例２による検出結果を表すグラフである。
【図８】  実験例２による検出結果を表すグラフである。
【図９】  実験例２による検出結果を表すグラフである。
【図１０】  従来の車体スリップ角の推定方法では車体スリップ角を正確に推定できないことがあることを説明する説明図である。
【符号の説明】
３０…電子制御装置（ＥＣＵ）、３１ＦＬ、３１ＦＲ、３１ＲＬ、３１ＲＲ…車輪速度センサ、３２…ヨーレートセンサ、３３…横加速度センサ、４４…ブレーキスイッチ。
【図１】
0004244453.tif 000007
【図２】
0004244453.tif 000008
【図３】
0004244453.tif 000009
【図４】
0004244453.tif 000010
【図５】
0004244453.tif 000011
【図６】
0004244453.tif 000012
【図７】
0004244453.tif 000013
【図８】
0004244453.tif 000014
【図９】
0004244453.tif 000015
【図１０】
0004244453.tif 000016