前後輪のタイヤモデルと車両モデルとに基づいて路面摩擦係数と車体すべり角との少なくとも1つを推定する走行状態推定手段を有する走行状態推定装置において、車速を検出する車速検出手段と、前輪がある地点を通過してから後輪が同地点を通過するまでの遅延時間を車速が高くなるほど短くなるように算出する遅延時間算出手段と、を備え、前記走行状態推定手段は、後輪タイヤモデルで用いられる後輪路面摩擦係数を、前記遅延時間分遅れた前輪タイヤモデルで用いられる前輪路面摩擦係数の値を用いることを特徴とする走行状態推定装置。請求項1に記載の走行状態推定装置において、前記走行状態推定手段は、前輪路面摩擦係数の急変を推定もしくは検出したときには、前輪路面摩擦係数が急変してから前記遅延時間の間、路面摩擦係数の推定速度を速くすることを特徴とする走行状態推定装置。請求項1に記載の走行状態推定装置において、前記走行状態推定手段は、前輪路面摩擦係数の急変を推定もしくは検出したときには、前輪路面摩擦係数が急変してから前記遅延時間の後に、車体すべり角の推定速度を速くすることを特徴とする走行状態推定装置。請求項1ないし請求項3のいずれか1項に記載の走行状態推定装置において、ドライバの操作入力をドライバ入力検出値として検出するドライバ入力検出手段と、前記ドライバ入力検出値と、前記前輪路面摩擦係数推定値と、前記遅延時間とから、後輪路面摩擦係数が変化するであろう時刻以降の車両状態を予測する車両状態予測手段と、この予測した前記車両状態に基づいて、車両運動制御の開始を予測する制御開始予測手段と、を有することを特徴とする走行状態推定装置。請求項4に記載の走行状態推定装置において、ブレーキもしくはクラッチをアクチュエータとする車両運動制御装置を有し、前記車両運動制御の開始を予測したときには、制御開始時刻までに、前記ブレーキもしくは前記クラッチのクリアランスを予め小さくしておくことを特徴とする走行状態推定装置。請求項1ないし請求項5のいずれか1項に記載の走行状態推定装置において、前記前輪路面摩擦係数の変化を検出した時刻と、前記後輪路面摩擦係数が変化すると予測される時刻との間で、後輪の制駆動力制限値を、変化前の路面摩擦係数相当の制限値から変化後の路面摩擦係数相当の制限値へ滑らかに変化させることを特徴とする走行状態推定装置。請求項1ないし請求項6のいずれか1項に記載の走行状態推定装置において、前記走行状態推定手段は、前輪を左右独立に路面摩擦係数を推定し、左後輪路面摩擦係数推定値を、前記遅延時間分遅れた左前輪路面摩擦係数推定値とし、右後輪路面摩擦係数推定値を、前記遅延時間分遅れた右前輪路面摩擦係数推定値とすることを特徴とする走行状態推定装置。請求項1ないし請求項7のいずれか1項に記載の走行状態推定装置において、前輪の横力もしくはセルフアライニングトルクを検出する前輪状態検出部と、車両状態を検出する車両状態検出部と、車両へのドライバの入力を検出するドライバ入力検出部と、を備え、前記走行状態推定手段は、前記前輪状態と、前記車両状態と、前記ドライバ入力とから、前後輪のタイヤモデルと車両モデルとに基づいて、少なくとも前記前輪路面摩擦係数を推定することを特徴とする走行状態推定装置。請求項8に記載の走行状態推定装置において、前記走行状態推定手段は、前記前輪状態と、前記車両状態と、前記ドライバ入力とから、前後輪のタイヤモデルと車両モデルとに基づいて、前記前輪路面摩擦係数と前記車体すべり角とを推定することを特徴とする走行状態推定装置。請求項9に記載の走行状態推定装置において、前輪の路面摩擦係数を推定もしくは検出する前輪路面摩擦係数推定部と、前記車両状態と、前記ドライバ入力と、前記前輪摩擦係数推定部とから、前後輪のタイヤモデルと車両モデルとに基づいて、少なくとも前記車体すべり角を推定する車体すべり角推定部と、を有することを特徴とする走行状態推定装置。しかしながら上記従来技術にあっては、摩擦特性モデルは各車輪とも同一としているため、摩擦係数が急変する路面を走行中には前輪の路面摩擦係数と後輪の路面摩擦係数とが大きく異なり、これによるステア特性の変化をモデルに反映することができず、横滑り角の推定精度が悪化するおそれがあった。本発明は、上記問題に着目されたもので、その目的とするところは、路面摩擦係数が急変する路面を走行しているときであっても、車両挙動の推定精度の悪化を抑制することができる走行状態推定装置を提供することである。上述の目的を達成するため、本発明では、前輪がある地点を通過してから後輪が同地点を通過するまでの遅延時間を車速が高くなるほど短くなるように算出し、後輪タイヤモデルで用いられる後輪路面摩擦係数を、遅延時間分遅れた前輪タイヤモデルで用いられる前輪路面摩擦係数の値を用いて、路面摩擦係数と車体すべり角との少なくとも1つを推定するようにした。よって本発明においては、路面摩擦係数が急変する路面を走行しているときであっても車両挙動の推定精度の悪化を抑制することができる走行状態推定装置を提供することである。実施例1図1は実施例1の走行状態推定装置を搭載した車両14を示す図である。この車両14は駆動系として、後輪4RL4RRに制駆動力を発生させる駆動モータ1と、リチウムイオンバッテリ3から駆動モータ1に電力を供給する駆動回路2とを有している。車両14の操舵系はステアバイワイヤであって、ステアリングホイール13の回転軸に取り付けられて操舵反力を発生させる反力モータ10と、前輪4FL4FRを転舵させる転舵モータ9L9Rとを有している。車両14は制御系として、駆動回路2、反力モータ10、転舵モータ11を制御する統合コントローラ8と、アクセル開度APOを検出するアクセル開度センサ12と、各車輪4FL4FR4RL4RRの車輪速Vwを検出する車輪速センサ5FL5FR5RL5RRと、車両重心のヨーレートγを検出するヨーレートセンサ6と、前後加速度Gxおよび横加速度Gyを検出する加速度センサ7と、ステアリングホイール13の回転軸に取り付けられて操舵角δを検出する操舵角センサ11と、転舵モータ9L9Rの回転角を検出する回転角センサ15L15Rとを有している。各センサの検出値は統合コントローラ8に入力されている。統合コントローラ8は、駆動モータ1の制駆動トルクを演算して、駆動回路2にトルク指令値を出力する。また統合コントローラ8は、転舵時に前輪4FL4FRに発生するセルフアライニングトルクを推定し、推定したセルフアライニングトルク相当の反力トルクを反力モータ10により発生させ、操舵反力としてドライバに伝える。また統合コントローラ8は、検出された操舵角δに応じて転舵モータ9L9Rによりはステアリングラックを車幅方向に変位させて、左右前輪4FL4FRを同相に転舵させる。[統合コントローラの構成]図2は統合コントローラ8の制御ブロック図である。図2に示すように、統合コントローラ8は、車速検出器81と、前後輪遅延時間算出器82と、前輪状態検出器83と、車両状態検出器84と、ドライバ入力検出器85と、状態推定器86と、車両制御器87とを有している。車速検出器車速検出器81では、車速車体速Vtを検出する。例えば、従動輪である前輪4FR4FLの車輪速Vwを基にして検出しても良いし、横滑り防止装置やGPS全地球測位システムなど、既存の車載装置において検出された値を用いても良い。前後輪遅延時間算出器前後輪遅延時間算出器82では、車速検出値Vtsensから前後輪遅延時間tdelayを求める。この前後輪遅延時間tdelayは、前輪4FL4FRがある地点を通過してから後輪4RL4RRが同地点を通過するまでの時間であり、状態推定器86において後輪路面摩擦係数後輪路面摩擦係数は、前後輪遅延時間tdelay前の前輪路面摩擦係数前輪路面摩擦係数を用いるようにしている。前後輪遅延時間tdelayは、例えば次の式を用いて計算できる。ここで、lwはホイールベース長である。また、旋回半径が小さい場合には、旋回時の軌跡から幾何的に導いた次の式から前後輪遅延時間tdelayを算出してもよい。ここでRtは旋回半径である。なお、この上記式では前後輪遅延時間tdelayは車速に反比例したものとなっているが、特にこれに限定せずに、単純な比例関係や、低車速でもより短い時間に設定したり、逆に高車速でも長い時間に設定したり、車速に関する変化の感度を適宜調整することが可能である。図3は前後輪遅延時間tdelayの設定例を示す図である。前後輪遅延時間tdelayは車速が高くなるほど短くなるという関係を満たすものであれば良く、例えば図3aに示すように単純比例一次関数の関係のものでも良いし、図3bに示すように上側に膨らむ二次関数的なものでも良いし、図3cに示すように下側に膨らむ二次関数的なものでも良いし、図3dに示すように階段状に減少するものであっても良い。前輪状態検出器前輪状態検出器83では、セルフアライニングトルクMzを検出する。例えば、転舵モータ9L9Rの指令電流からモータトルク定数や、サスペンションリンクの構造を考慮して、セルフアライニングトルクMzを検出する。前輪状態はセルフアライニングトルクMzに限られるものではなく、前輪タイヤ横力、前輪スリップ率など、前輪路面摩擦係数に応じて高感度に影響を受ける物理量ならば、同様に用いることができる。その場合、以下実施例中のセルフアライニングトルクMzを対応する物理量に読み替える。車両状態検出器車両状態検出器84では、車両重心回りのヨーレートγを検出する。例えば、前述のヨーレートセンサ6の検出値を用いても良いし、既存の横すべり防止装置に内蔵されたセンサで計測された値を用いても良い。またヨーレートγよりも横加速度Gyのほうが精度良く検出できる場合は、横加速度Gyで代用しても良い。ドライバ入力検出器ドライバ入力検出器85では、操舵角δを検出する。例えば、前述の操舵角センサ11の検出値を用いても良いし、既存の横すべり防止装置に内蔵されたセンサで計測された値を用いても良い。状態推定器状態推定器86は、前後輪遅延時間tdelayと、セルフアライニングトルク検出値Mzsensと、ヨーレート検出値γsensと、操舵角検出値δsensとを入力し、車体すべり角推定値βestと、前輪路面摩擦係数推定値μfestとを求める。状態推定器86の詳しい構成については、制御ブロック図を用いて後で述べる。車両制御器車両制御器87は、前後輪遅延時間tdelayと、前輪路面摩擦係数推定値μfestと、車体すべり角推定値βestとを入力し、後述する制動器87hを制御する。車両制御器87の詳しい構成については、制御ブロック図を用いて後で述べる。[状態推定器の構成]図4は状態推定器86の制御ブロック図である。図4に示すように状態推定器86は、車両ダイナミクスモデル86aと、前輪路面摩擦係数ダイナミクスモデル86bと、後輪路面摩擦係数モデル86cと、セルフアライニングトルク算出器86dと、偏差算出器86eと、補正量算出器86fと、路面μ補正ゲイン増加判定器86gと、車両状態補正ゲイン判定器86hとを有している。車両ダイナミクスモデル車両ダイナミクスモデル86aは、操舵角検出値δsensと、後述する前輪路面摩擦係数ダイナミクスモデル86bで推定された前輪路面摩擦係数推定値μfestと、後述する後輪路面摩擦係数モデル86cで推定された後輪路面摩擦係数推定値μrestと、後述する補正量算出器86fで求められた車体すべり角補正量Cβと、ヨーレート補正量Cγとを入力し、車体すべり角推定値βestと、前輪路面摩擦係数推定値μfestとを求める。2輪モデル車両の運動と制御第2版、山海堂、安部正人著を離散時間系に変換することで導出した次の式に基づいて、車体すべり角βおよびヨーレートγを推定する。Δtは1計算ステップあたりの時間である。Izはヨー方向の車体慣性モーメント、mは車両の質量、lfは車両重心点と前輪車軸間の距離、lrは車両重心点と後輪車軸間の距離である。添え字kは現在の計算ステップにおける状態であることを表し、添え字k1は次回計算ステップにおける状態であることをあらわす。上の式を用いて、現在計算ステップの状態から次回計算ステップの状態を逐次的に計算できる。式においてFyfは前輪タイヤ横力であり、前輪タイヤ横力Fyfは、例えば次の式のフィアラモデルを用いて求められる。Fzf0は静的な前輪輪荷重である。またφfは次の式によって求められる。μfestは前輪路面摩擦係数ダイナミクスモデル86bにより推定された前輪路面摩擦係数推定値である。Cαfは前輪コーナリングパワーである。αfは前輪すべり角であり、次の式により求められる。δsensはドライバ入力検出器85により検出された操舵角検出値である。式においてFyrは後輪タイヤ横力であり、後輪タイヤ横力Fyrは、例えば次の式のフィアラモデルを用いて求められる。Fzr0は静的な後輪輪荷重である。またφrは次の式によって求められる。μrestは後輪路面摩擦係数モデル86cにより推定された後輪路面摩擦係数推定値である。Cαrは後輪コーナリングパワーである。αrは後輪滑り角であり、次の式により求められる。この車両ダイナミクスモデル86aにおいて求められた車体すべり角推定値βestは、状態推定器86の出力値とされる。前輪路面摩擦係数ダイナミクスモデル前輪路面摩擦係数ダイナミクスモデル86bは、前輪路面摩擦係数補正量Cμfを入力し、前輪路面摩擦係数ダイナミクスモデルに基づいて前輪路面摩擦係数推定値μfestを求める。前輪路面摩擦係数推定値μfestは、次の式により求められる。前輪路面摩擦係数補正量Cμfは前輪路面摩擦係数に関する補正量であり、補正量算出器86fにより求められる。この前輪路面摩擦係数ダイナミクスモデル86bにおいて求められた前輪路面摩擦係数μfestは、状態推定器86の出力値とされる。後輪路面摩擦係数モデル後輪路面摩擦係数モデル86cは、前輪路面摩擦係数推定値μfestと、前後輪遅延時間tdelayとを入力し、後輪路面摩擦係数μrestを求める。まず、次の式より前後輪遅延時間tdelayから遅延ステップ数Ndelayを求める。tdelayは、前後輪遅延時間算出器82で求められた前後輪遅延時間である。floorは引数が小数の場合は、小数点以下を切り捨てて整数に変換する演算子を表わす。後輪路面摩擦係数推定値μrestは、現在の前輪路面摩擦係数μfestに対して遅延ステップ数Ndelay前の前輪路面摩擦係数μfestの値に設定される。セルフアライニングトルク算出器セルフアライニングトルク算出器86dは、前輪摩擦係数推定値μfestと、車体すべり角推定値βestと、ヨーレート推定値γestを入力し、セルフアライニングトルク推定値Mzestを求める。セルフアライニングトルク推定値Mzestは、例えば次の式のフィアラモデルを用いて求められる。ここで、lcはタイヤ接地長であり、例えばタイヤ空気圧とタイヤ空気圧の関係を実験的に取得したマップを用いて、タイヤ空気圧から推定される。またφfは次の式によって求められる。μfestは前輪路面摩擦係数ダイナミクスモデル86bにより推定された前輪路面摩擦係数推定値である。Cαfは前輪コーナリングパワーである。αfは前輪すべり角であり、次の式により求められる。偏差算出器偏差算出器86eは、車両ダイナミクスモデル86aで求められたヨーレート推定値γestと、車両状態検出器84で検出されたヨーレート検出値γsensとの偏差eγ、およびセルフアライニングトルク算出器86dで求められたセルフアライニングトルク推定値Mzestと、前輪状態検出器83で検出されたセルフアライニングトルク検出値Mzsensとの偏差eMzを求める。ヨーレート推定値γestとヨーレート検出値γsensとの偏差eγは次の式により求められる。セルフアライニングトルク推定値Mzestとセルフアライニングトルク検出値Mzsensとの偏差eMzは次の式により求められる。補正量算出器補正量算出器86fは、偏差算出器86eで求めた偏差eγeMzを入力し、車体すべり角補正量Cβと、ヨーレート補正量Cγと、前輪路面摩擦係数補正量Cμfとを求める。車体すべり角補正量Cβ、ヨーレート補正量Cγ、前輪路面摩擦係数補正量Cμfは、それぞれ次の式~により求められる。補正ゲインKcβ1Kcβ2Kcγ1Kcγ2Kcμf1Kcμf2は、設計者により適当に定められる。補正ゲインは、各々の絶対値が大きいほど補正量の絶対値が大きくなり、車体すべり角β、ヨーレートγ、前輪路面摩擦係数μfが急に変化した場合でも各々の推定値をすばやく追従させることでき高い応答性が実現できる。一方、補正ゲインの絶対値が大きいときには、車体すべり角推定値βest、ヨーレート推定値γest、前輪路面摩擦係数推定値μfestが、ヨーレート検出値γsensおよびセルフアライニングトルク検出値Mzsensに含まれるノイズの影響を受け振動的になる。したがって、応答性の設計値と、検出値に含まれるノイズ強度とのバランスを考慮して、補正ゲインを適当に定めることとする。なお、走行状況に応じた各補正ゲインの調節は、路面摩擦係数補正ゲイン増加判定器86gおよび車両状態補正ゲイン増加判定器86hの信号に基づいて行う。路面摩擦係数補正ゲイン増加判定器路面摩擦係数補正ゲイン増加判定器86gは、操舵角検出値δsensとセルフアライニングトルク検出値Mzsensとを入力し、路面摩擦係数補正ゲイン増加判定信号を出力するか否かを判定する。図5は路面摩擦係数補正ゲイン増加判定器86gの制御ブロック図である。路面摩擦係数補正ゲイン増加判定器86gは、操舵角速度推定器86g1と、セルフアライニングトルク変化率推定器86g2と、比計算器86g3と、セルフアライニングトルク変化判定器86g4とを有している。操舵角速度推定器86g1は、操舵角検出値δsensを入力し、操舵角速度推定値rδestを求める。操舵角速度推定値rδestは次の式で求められる。ここで、ΔTはサンプル時間、kは現在の計算時間ステップを表わす。操舵角検出値δsensにノイズが含まれるときには、操舵角速度推定値rδestを適当な特性を有するローパスフィルターで処理した後に出力しても良い。セルフアライニングトルク変化率推定器86g2は、セルフアライニングトルク検出値Mzsensを入力し、セルフアライニングトルク推定変化率rMzestを求める。セルフアライニングトルク推定変化率rMzestは次の式で求められる。比計算器86g3は、操舵角速度推定値rδestと、セルフアライニングトルク推定変化率rMzestとを入力し、これらの比λを求める。比λは次の式で求められる。セルフアライニングトルク変化判定器86g4は、比λを入力し、路面摩擦係数増加信号を出力するか否かを判定する。比λの絶対値が、予め定めた閾値よりも大きいときには、セルフアライニングトルクMzの変化が操舵角δの変化ではなく、路面摩擦係数μが比較的大きく変化したと判断する。そして、今後も路面摩擦係数μが急速に変化することが見込まれるとして、前輪路面摩擦係数推定値μfestの補正に用いる路面摩擦係数補正ゲインKcμf1Kcμf2の絶対値を大きくするように路面摩擦補正ゲイン増加信号を出力する。なお、路面摩擦補正ゲイン増加信号は、比λの絶対値が予め定めた閾値よりも大きくなったと判定してから前後輪遅延時間tdelayの間に出力される。路面摩擦係数補正ゲインの増加判定の方法は、上記の方法に限られるものではない。例えば、路面を記録した画像に基づいて路面摩擦係数μの変化を判定し、前輪路面摩擦係数推定値μfestの補正ゲインの増加を判定してもよい。すなわち、前方のバンパーに備えたCCDカメラで路面を記録し、これを画像処理することにより、湿潤した路面か、乾燥した路面かを判断する。ここで、乾燥した路面から湿潤した路面に変化したこと、もしくは、湿潤した路面から乾燥した路面に変化したことが検知されたときには、前輪路面摩擦係数μfも変化することが見込まれるため、路面摩擦係数補正ゲイン増加信号を出力する。車両状態補正ゲイン増加判定器車両状態補正ゲイン増加判定器86hは、前輪路面摩擦係数推定値μfestと前後輪遅延時間tdelayとを入力し、車両状態補正ゲイン増加信号を出力するか否かを判定する。所定の時間内に、前輪路面摩擦係数推定値μfestが予め定められた変化量より大きく変化したときには、この変化を検知した時間から前後輪遅延時間tdelay後に、車体すべり角補正ゲインKcβと、ヨーレート補正ゲインKcγの値を大きくするように車両状態補正ゲイン増加信号を出力する。[車両制御器の構成]図6は車両制御器87の制御ブロック図である。図6に示すように車両制御器87は、前輪路面摩擦係数変化時刻判定器87aと、後輪路面摩擦係数変化予測時間算出器87bと、後輪路面摩擦係数モデル87cと、ドライバ入力予測器87dと、車速予測器87eと、車両状態予測器87fと、閾値比較器87fgと、制動器87hとを有している。前輪路面摩擦係数変化時刻判定器前輪路面摩擦係数変化時刻判定器87aは、前輪路面摩擦係数推定値μfestと前後輪遅延時間tdelayとを入力して、所定の時間内に所定の値よりも大きな範囲で前輪路面摩擦係数推定値μfestが変化したときに、前輪路面摩擦係数推定値μfestが変化を開始した時刻を、前輪路面摩擦係数変化時刻tμfとして出力する。後輪路面摩擦係数変化予測時間算出器後輪路面摩擦係数変化予測時間算出器87bは、前輪路面摩擦係数変化時刻tμfと、前後輪遅延時間tdelayと、現在時刻tとから後輪路面摩擦係数変化予測時間Δtchangeを算出する。後輪路面摩擦係数変化予測時間Δtchangeは、現在時刻tから何秒後に後輪路面摩擦係数μrが変化するかを示す。後輪路面摩擦係数変化予測時間Δtchangeは、次の式により求められる。後輪路面摩擦係数モデル後輪路面摩擦係数モデル87cは、前輪路面摩擦係数推定値μfestと前後輪遅延時間tdelayを入力し、後輪路面摩擦係数μrestを求める。後輪路面摩擦係数μrestの求め方は、状態推定器86の後輪路面摩擦係数モデル86cと同様である。ドライバ入力予測器ドライバ入力予測器87dは、後輪路面摩擦係数変化予測時間Δtchangeと、操舵角検出値δsensとを入力し、操舵角予測値δpreを求める。図7は操舵角予測値δpreのマップである。現在からNサンプル前までの操舵角検出値δsensを直線でフィッティングし、この直線を延長して後輪路面摩擦係数変化予測時間Δtchange後の操舵角δを予測する。フィッティングする対象の操舵角検出値の数Nは、設計者により適当に定められる。フィッティングのアルゴリズムには、例えば最小2乗法を用いればよい。ドライバ入力の予測は上記の方法に限定されるものではなく、例えば現在の操舵角δが将来も引き続き保持されるとして、現在の操舵角δを用いて車両状態を予測することとしても良い。車速予測器車速予測器87eは、後輪路面摩擦係数変化予測時間Δtchangeと、車速検出値Vtsensとを入力し、車速予測値Vtpreを求める。図8は車速予測値Vtpreのマップである。現在からNサンプル前までの車速検出値Vtsensを直線でフィッティングし、この直線を延長して後輪路面摩擦係数変化予測時間Δtchange後の車速Vtを予測する。フィッティングする対象の操舵角検出値の数Nは、設計者により適当に定められる。フィッティングのアルゴリズムには、例えば最小2乗法を用いればよい。車速の予測は上記の方法に限定されるものではなく、例えば現在の車速Vtが将来も引き続き保持されるとして、現在の車速Vtを用いて車両状態を予測することとしても良い。車両状態予測器車両状態予測器87fでは、操舵角予測値δpreと、車速予測値Vtpreと、後輪路面摩擦係数変化予測時間Δtchangeと、前輪路面摩擦係数推定値μfestと、後輪路面摩擦係数推定値μrestと、車体すべり角推定値βestと、ヨーレート検出値γsensとを入力し、車体すべり角予測値βpreと、ヨーレート予測値γpreとを求める。2輪モデル車両の運動と制御第2版、山海堂、安部正人著を離散時間系に変換することで導出した次の式に基づいて、車体すべり角およびヨーレートを推定する。Δtは1計算ステップあたりの時間である。Izはヨー方向の車体慣性モーメント、mは車両の質量、lfは車両重心点と前輪車軸間の距離、lrは車両重心点と後輪車軸間の距離である。添え字kは現在の計算ステップにおける状態であることを表し、添え字k1は次回計算ステップにおける状態であることをあらわす。式の求めた車体すべり角予測値βpreとヨーレート予測値γpreの初期値は、状態推定器86の車両ダイナミクスモデル86aにおける式で求めた車体すべり角予測値βestとヨーレート予測値γestを用いることとする。式においてFyfは前輪タイヤ横力であり、前輪タイヤ横力Fyfは、例えば次の式のフィアラモデルを用いて求められる。Fzf0は前輪輪荷重である。またφfは次の式によって求められる。Cαfは前輪コーナリングパワーである。αfは前輪すべり角であり、次の式により求められる。式においてFyrは後輪タイヤ横力であり、後輪タイヤ横力Fyrは、例えば次の式のフィアラモデルを用いて求められる。Fzr0は後輪輪荷重である。またφrは次の式によって求められる。Cαrは後輪コーナリングパワーである。αrは後輪滑り角であり、次の式により求められる。式の反復計算を、少なくとも後輪路面摩擦係数変化予測時間Δtchangeよりも後の時間まで所定回数繰り返して算出した値を車体すべり角予測値βpreとして出力する。また、式の反復計算を同様に繰り返して算出した値をヨーレート予測値γpreとして出力する。閾値比較器閾値比較器87gは、車体すべり角予測値βpreと、ヨーレート予測値γpreとを入力し、予め定めた車体すべり角閾値βthresと、ヨーレート閾値γthresとを比較した結果に基づいて制御開始信号を出力する。例えば、車体すべり角予測値βpreの絶対値が車体すべり角閾値βthresよりも大きいとき、あるいはヨーレート予測値γpreの絶対値がヨーレート閾値γthresよりも大きいときには、操舵角等のドライバ入力の傾向が維持されれば、後輪路面摩擦係数変化予測時間Δtchange後に車両が不安定となる可能性が高いとして、制御開始信号を出力する。すなわち、または、であるときに制御開始信号を出力する。制動器制動器87hは、制御開始信号に基づいてブレーキ与圧を付加する。ブレーキの与圧を付加する輪は、制動力により旋回方向と逆方向にヨーモーメントを発生することができる輪を選択する。例えば左旋回時には、右前輪と右後輪に備えられたブレーキに与圧を付加し、右旋回時には、左前輪と左後輪に備えられたブレーキに与圧を付加する。ブレーキ与圧の強度は、リターンスプリングによる圧力ないしは、これよりもやや大きな圧力を付加することとする。なお、旋回方向は、ヨーレート検出値の符号により判定する。[作用]車両挙動推定の高精度化従来は、前輪と後輪の路面摩擦係数μは同一であるとして、車両挙動を推定していた。しかし、例えば乾燥路面から湿潤路面へ進入するときなどのように、路面摩擦係数が急変するときには、車両挙動の推定精度が悪化してしまうことがあった。そこで実施例1では前後輪遅延時間算出器82において、前輪4FL4FRがある地点を通過してから後輪4RL4RRが同地点を通過するまでの時間である前後輪遅延時間tdelayを演算するようにした式または参照。また、後輪路面摩擦係数モデル86cにおいて、後輪路面摩擦係数推定値μrestは、現在の前輪路面摩擦係数μfestに対して前後輪遅延時間tdelay遅延ステップ数Ndelay前の前輪路面摩擦係数μfestの値に設定するようにした式参照。車両ダイナミクスモデル86aでは、この後輪路面摩擦係数モデル86cを用いて、車体すべり角推定値βestを推定するようにした式~参照。図9は車体すべり角βと前輪路面摩擦係数μfのシミュレーション結果を示す図である。図9aは後輪路面摩擦係数μrとして前後輪遅延時間tdelay分遅れた前輪路面摩擦係数μfの値を用いたときのシミュレーション結果である。図9bは後輪路面摩擦係数μrとして現在の前輪路面摩擦係数μfの値を用いたときのシミュレーション結果である。図9aに示すように、後輪路面摩擦係数μrとして前後輪遅延時間tdelay分遅れた前輪路面摩擦係数μfの値を用いた場合は、路面摩擦係数が急変するときであっても、車体すべり角推定値βestは真値に追従している。一方、後輪路面摩擦係数μrとして現在の前輪路面摩擦係数μfの値を用いた場合は、路面摩擦係数が急変するときには車体すべり角推定値βestは真値に対して誤差が発生している。すなわち、後輪路面摩擦係数μrとして前後輪遅延時間tdelay分遅れた前輪路面摩擦係数μfの値を用いることにより、路面摩擦係数が急変するときであっても、車両挙動の推定精度を高く保つことができる。前輪路面摩擦係数推定の高応答性また実施例1では路面摩擦係数補正ゲイン増加判定器86gにおいて、操舵角速度推定値rδestに対するセルフアライニングトルク推定変化率rMzestの比λが予め設定した閾値よりも大きくなったと判定してから前後輪遅延時間tdelayの間、前輪路面摩擦係数推定値μfestの補正に用いる補正ゲインKcμf1Kcμf2の絶対値を大きくするように路面摩擦補正ゲイン増加信号を出力するようにした。路面摩擦係数補正ゲインKcμf1Kcμf2の絶対値が大きいほど前輪路面摩擦係数補正量Cμfの絶対値が大きくなり式参照、前輪路面摩擦係数推定値μfestをすばやく追従させることでき式参照、高い応答性を実現できる。車両状態推定の高応答性また実施例1では車両状態補正ゲイン増加判定器86hにおいて、所定の時間内に、前輪路面摩擦係数推定値μfestが予め定められた変化量より大きく変化したときには、この変化を検知した時間から前後輪遅延時間tdelay後に、車体すべり角補正ゲインKcβ1Kcβ2とヨーレート補正ゲインKcγ1Kcγ2の値を大きくするように車両状態補正ゲイン増加信号を出力するようにした。図10は、旋回時に路面摩擦係数が高い路面から低い路面に進入した場合の車体すべり角βとヨーレートγのタイムチャートをである。区間1は前輪路面摩擦係数μfも後輪路面摩擦係数μrも高い区間、区間2は前輪路面摩擦係数μfは低く後輪路面摩擦係数μrは高い区間、区間3は前輪路面摩擦係数μfも後輪路面摩擦係数μrも低い区間とする。区間3に変化した直後、前輪路面摩擦係数μfだけでなく後輪路面摩擦係数μrも低い状態となり、ステア特性がオーバーステアに変化したことで、車体すべり角βおよびヨーレートγが急激に変化を始めることが分かる。一方、前輪路面摩擦係数μfの変化が起こったことは区間2の初期に検知が可能であり、この時点から前後輪遅延時間tdelay経過した後に、車体すべり角βおよびヨーレートγが急激に変化を始める。したがって、前輪路面摩擦係数推定値μfestが予め定められた変化量より大きく変化したときには、この変化を検知した時間から前後輪遅延時間tdelay後に、車両状態補正ゲイン増加信号を出力する。これにより、車体すべり角補正ゲインKcβ1Kcβ2の絶対値が大きいほど車体すべり角補正量Cβの絶対値が大きくなり式参照、車体すべり角推定値βestをすばやく追従させることでき式参照、高い応答性が実現できる。同じく、ヨーレート補正ゲインKcγ1Kcγ2の絶対値が大きいほどヨーレート補正量Cγの絶対値が大きくなり式参照、車体すべり角推定値βestをすばやく追従させることでき式参照、高い応答性を実現できる。車両運動制御の開始予測また実施例1では車両状態予測器87fにおいて、後輪路面摩擦係数変化予測時間Δtchange後の車体すべり角予測値βpreおよびヨーレート予測値γpreを予測するようにした式~参照。また閾値比較器87gにおいて、車体すべり角予測値βpreの絶対値が車体すべり角閾値βthresよりも大きいとき、あるいはヨーレート予測値γpreの絶対値がヨーレート閾値γthresよりも大きいときには、操舵角等のドライバ入力の傾向が維持されれば、後輪路面摩擦係数変化予測時間Δtchange後に車両が不安定となる可能性が高いとし制御開始信号を出力するようにした。よって、この予測に基づいて適切な車両制御を実施することができる。ブレーキ予圧また実施例1では制動器87hにおいて、制御開始信号が入力されるとてブレーキ与圧を付加して、ブレーキのクリアランスを予め小さくするようにした。よって、予めブレーキ与圧を付加することで、車両挙動が実際に不安定化したときに、より素早く制動力を生じさせることができ、高応答に車両運動を安定化することができる。図11は前輪路面摩擦係数μfが低くなり制御開始信号が出力されたとしたときのブレーキ圧のタイムチャートである。前輪路面摩擦係数μfが低くなった後に制御開始信号が出力されるとブレーキ与圧が付加される。さらに、後輪路面摩擦係数μrが低くなると車両運動を安定化するヨーモーメントを発生するのに十分なブレーキ圧を付加する。前輪路面摩擦係数推定の高精度化また実施例1の状態推定器86では、セルフアライニングトルク検出値Mzsens、ヨーレート検出値γsens、操舵角検出値δsensから前輪路面摩擦係数推定値μfestを推定することとした。よって、前輪路面摩擦係数μfに応じて高感度に変化するセルフアライニングトルクMzを用いるため、前輪路面摩擦係数μfの推定精度を高くすることができる。前輪路面摩擦係数と車体すべり角を同時推定また実施例1の状態推定器86では、セルフアライニングトルク検出値Mzsens、ヨーレート検出値γsens、操舵角検出値δsensから前輪路面摩擦係数推定値μfestおよび車体すべり角推定値βestを推定することとした。よって、互いに関連性がある前輪路面摩擦係数μfと車体すべり角βを同時に推定でき、前輪路面摩擦係数推定値μfestと車体すべり角推定値βestの推定精度を高くすることができる。車体すべり角推定の高精度化また実施例1の状態推定器86では、前輪路面摩擦係数推定値μfestを推定し、操舵角検出値δsens、前輪路面摩擦係数推定値μfestから車体すべり角推定値βestを推定することとした。よって、前輪路面摩擦係数推定値μfestから精度良く車体すべり角推定値βestを推定できる。[効果]実施例1の効果について以下に列記する。1前後輪のタイヤモデルと車両モデルとに基づいて路面摩擦係数μと車体すべり角βとの少なくとも1つを推定する状態推定器86走行状態推定手段を有する走行状態推定装置において、車速を検出する車速検出器81車速検出手段と、前輪がある地点を通過してから後輪が同地点を通過するまでの前後輪遅延時間tdelayを車速Vtが高くなるほど短くなるように算出する前後輪遅延時間算出器82遅延時間算出手段と、を備え、状態推定器86は、後輪路面摩擦係数モデル86cで用いられる後輪路面摩擦係数μrを、前後輪遅延時間tdelay分遅れた前輪路面摩擦係数ダイナミクスモデル86bで用いられる前輪路面摩擦係数推定値μfestの値を用いるようにした。よって、路面摩擦係数が急変するときであっても、車両挙動の推定精度を高く保つことができる。2状態推定器86は、前輪路面摩擦係数μfestの急変を推定もしくは検出したときには、前輪路面摩擦係数μfが急変してから前後輪遅延時間tdelayの間、前輪路面摩擦係数μfの推定速度を速くするようにした。よって、前輪路面摩擦係数推定値μfestをすばやく追従させることでき高い応答性を実現できる。3状態推定器86は、前輪路面摩擦係数μfの急変を推定もしくは検出したときには、前輪路面摩擦係数μfが急変してから前後輪遅延時間tdelayの後に、車体すべり角βの推定速度を速くするようにした。よって、車体すべり角推定値βestをすばやく追従させることでき高い応答性を実現できる。4ドライバの操作入力を操舵角検出値δsensドライバ入力検出値として検出するドライバ入力検出器85ドライバ入力検出手段と、操舵角検出値δsensと、前輪路面摩擦係数推定値μfestと、前後輪遅延時間tdelayとから、後輪路面摩擦係数が変化するであろう時刻以降の車両状態車体すべり角β、ヨーレートγを予測する車両状態予測器87f車両状態予測手段と、この予測した車両状態に基づいて、車両運動制御の開始を予測する閾値比較器制御開始予測手段とを設けた。よって、この予測に基づいて適切な車両制御を実施することができる。5ブレーキをアクチュエータとする制動器87h車両運動制御装置を有し、車両運動制御の開始を予測したときには、制御開始時刻までに、ブレーキのクリアランスを予め小さくしておくようにした。よって、予めブレーキ与圧を付加することで、車両挙動が実際に不安定化したときに、より素早く制動力を生じさせることができ、高応答に車両運動を安定化することができる。6セルフアライニングトルク検出値Mzsensを検出する前輪状態検出器83前輪状態検出部と、ヨーレート検出値γ車両状態を検出する車両状態検出器84車両状態検出部と、操舵角検出値δsens車両へのドライバの入力を検出するドライバ入力検出器85ドライバ入力検出部と、を備え、状態推定器86は、セルフアライニングトルク検出値Mzsensと、ヨーレート検出値γと、操舵角検出値δsensとから、前後輪のタイヤモデルと車両モデルとに基づいて、少なくとも前輪路面摩擦係数推定値μfestを求めるようにした。よって、前輪路面摩擦係数μfに応じて高感度に変化するセルフアライニングトルクMzを用いるため、前輪路面摩擦係数μfの推定精度を高くすることができる。7状態推定器86は、セルフアライニングトルク検出値Mzsensと、ヨーレート検出値γと、操舵角検出値δsensとから、前後輪のタイヤモデルと車両モデルとに基づいて、前輪路面摩擦係数推定値μfestと車体すべり角βestとを求めるようにした。よって、互いに関連性がある前輪路面摩擦係数μfと車体すべり角βを同時に推定でき、前輪路面摩擦係数推定値μfestと車体すべり角推定値βestの推定精度を高くすることができる。8前輪路面摩擦係数推定値μfestを求める前輪路面摩擦係数ダイナミクスモデル86b前輪路面摩擦係数推定部と、ヨーレート検出値γと、操舵角検出値δsensと前輪路面摩擦係数推定値μfestとから、前後輪のタイヤモデルと車両モデルとに基づいて、少なくとも車体すべり角推定値βestを求め車体ダイナミクスモデル86a車体すべり角推定部とを設けた。よって、前輪路面摩擦係数推定値μfestから精度良く車体すべり角推定値βestを推定できる。実施例2実施例1では、車両制御器87の閾値比較器87gから制御開始信号が出力されると、制動器87hはブレーキ予圧を負荷するようにしていた。実施例2では、車両制御器87に制御開始判定器87mを設け、この制御開始判定器87mから制御開始信号が出力されると、滑らかに後輪タイヤ縦力Fxrが変化するように後輪制駆動トルク指令値Tdを出力するようにした。実施例2は、状態推定器86と車両制御器87の構成が実施例1と異なる。以下、実施例2について説明するが、実施例1と同じ構成については同一の符号を付して説明を省略する。[状態推定器の構成]図12は、状態推定器86の制御ブロック図である。図12にしめすように、状態推定器86は、車両ダイナミクスモデル86aと、後輪路面摩擦係数モデル86cと、前輪路面摩擦係数推定器86iと、偏差算出器86jと、車両状態補正量算出器86kとを有している。車両ダイナミクスモデル86aと後輪路面摩擦係数モデル86cは実施例1と同じものである。前輪路面摩擦係数推定器前輪路面摩擦係数推定器86iは、前輪路面摩擦係数推定値μfestを求める。例えば、前輪の近くに備えられた路面を撮影するカメラから得られる画像を解析することで、アスファルトや雪道などの路面の表面構成と、路面表面がどの程度湿潤しているかを判定し、前輪路面摩擦係数μfを推定する。また、ナビゲーションシステムの情報やワイパーの作動頻度を参照して、前輪路面摩擦係数μfを推定するようにしても良い。この前輪路面摩擦係数推定器86iにより求められた前輪路面摩擦係数推定値μfestは、状態推定器86の出力値とされる。偏差算出器偏差算出器86jは、車両ダイナミクスモデル86aで求められたヨーレート推定値γestと、車両状態検出器84で検出されたヨーレート検出値γsensとの偏差eγを求める。ヨーレート推定値γestとヨーレート検出値γsensとの偏差eγは次の式により求められる。車両状態補正量算出器車両状態補正量算出器86kは、偏差算出器86jで求めた偏差eγを入力し、車体すべり角補正量Cβと、ヨーレート補正量Cγとを求める。車体すべり角補正量Cβ、ヨーレート補正量Cγは、それぞれ次の式により求められる。補正ゲインKcβKcγは、設計者により適当に定められる。[車両制御器の構成]図13は車両制御器87の制御ブロック図である。図13に示すように車両制御器87は、後輪制駆動トルク検出器87iと、後輪タイヤ縦力算出器87jと、後輪スリップ率推定器87kと、後輪タイヤ横力推定器87lと、制御開始判定器87mと、目標後輪タイヤ縦力応答生成器87nと、後輪タイヤ縦力偏差算出器87pと、後輪制駆動トルク指令値算出器87qとを有している。後輪制駆動トルク検出器後輪制駆動トルク検出器87iは、後輪4RL4RRの駆動モータ1から後輪制駆動トルク検出値Trsensを検出する。例えば、検出された駆動モータ1の電流にトルク定数を乗算することで、回生制動時および駆動時の制駆動トルク検出値Trsensを検出する。後輪タイヤ縦力算出器後輪タイヤ縦力算出器87jは、後輪制駆動トルク検出値Trsensを入力し、後輪タイヤ縦力検出値Fxrsensを求める。後輪制駆動トルク検出値Trsensにタイヤ有効半径を乗算することでタイヤ周での相当する力に変換する。また、後輪制駆動トルクTrが変化してから後輪タイヤ縦力Fxrが変化するまでにはタイヤの変形に起因する遅れが存在するため、この遅れ特性を考慮して後輪タイヤ縦力検出値Fxrsensを求める。後輪スリップ率推定器後輪スリップ率推定器87kは、例えば4輪の車輪速Vwから推定された車体速Vtと、後輪の車輪速Vwに基づいて、後輪推定スリップ率Srestを求める。後輪タイヤ横力推定器後輪タイヤ横力推定器87lは、車体すべり角推定値βestと、ヨーレート検出値γsensと、後輪推定スリップ率Srestとを入力し、後輪タイヤ横力推定値Fyrestを求める。まず、後輪すべり角αrを次の式により求める。この後輪すべり角αrと、後輪推定スリップ率Srestとから、例えばブラッシュモデル車両の運動と制御第2版、山海堂、安部正人著等のタイヤモデルを用いて、後輪タイヤ横力Fyrestを求める。タイヤモデルはこれに限られるものではなく、実験的に同定したマジックフォーミュラ等を用いても良い。制御開始判定器制御開始判定器87mは、後輪タイヤ横力推定値Fyrestと、前輪路面摩擦係数推定値μfestと、前後輪遅延時間tdelayとを入力し、制駆動力制限制御を開始するか否か判断する。そして、開始するべきであると判断したときには制御開始信号を出力する。例えば、次の条件を満たすときに制御開始信号を出力する。条件1:前輪路面摩擦係数推定値μfestが、所定の時間内に所定変化幅以上に低下すること。条件2:後輪タイヤ横力推定値fyrestが、前輪路面摩擦係数推定値μfestから算出される摩擦円半径よりも小さいこと。すなわち、次の式を満たしていること。ここで、Fzr0は静的な後輪輪荷重である。条件3:タイヤ力が、前輪路面摩擦係数推定値μfestから算出される摩擦円半径よりも大きいこと。すなわち、次の式を満たしていること。以上の条件1から条件3を全て満たすときに制御開始信号を出力する。以下の目標後輪タイヤ縦力応答生成器87n、後輪タイヤ縦力偏差算出器87p、後輪制駆動トルク指令値算出器87qにおける演算は、制御開始信号が出力されたときに実行する。目標後輪タイヤ縦力応答生成器目標後輪タイヤ縦力応答生成器87nは、前輪路面摩擦推定値μfestと、後輪タイヤ縦力検出値Fxrsensと、前後輪遅延時間tdelayとから、目標後輪タイヤ縦力応答Fxrtarを生成する。まず、前輪路面摩擦推定値μfestから算出される許容される後輪タイヤ力許容値は、前輪路面摩擦係数推定値μfestと後輪輪荷重Fzr0との積であたえられる。制御の結果、前後輪遅延時間tdelay後に実現したい後輪タイヤ力は、この後輪タイヤ力許容値よりも小さい値でなければならない。したがって、後輪タイヤ縦力許容値Fxrtolを前後輪遅延時間tdelay後に次の式で求める値とすれば良い。ここで、Sはマージンを設定する係数であり、0より大きく1より小さい値に設定する。図14は、目標後輪タイヤ縦力応答Fxrtarのマップである。目標後輪タイヤ縦力応答Fxrtarは、例えば図14に示すように、現在時刻と現在時刻における後輪タイヤ縦力とでプロットされる図中A点と、現在時刻から前後輪遅延時間後の時刻と許容後輪タイヤ縦力とでプロットされる図中B点とを結ぶ線分を目標後輪タイヤ縦力応答Fxrtarとすれば良い。後輪縦力偏差算出器後輪タイヤ縦力偏差算出器87pでは、後輪タイヤ縦力検出値Fxrsensと目標後輪タイヤ縦力応答Fxrtarとの偏差を算出し、後輪タイヤ縦力偏差eFxrとして出力する。後輪制駆動トルク指令値算出器後輪制駆動トルク指令値算出器87qでは、後輪タイヤ縦力偏差eFxrから後輪制駆動トルク指令値Tdを算出する。後輪制駆動トルク指令値Tdの算出には、例えばPI制御器を用いる。すなわち、入出力関係を伝達関数で表現すると、次の式に示すようになる。ここで、Kpは比例ゲイン、Kiは積分ゲインである。[作用]実施例2では、前輪路面摩擦係数μfの変化を検出した時間から前後輪遅延時間tdelayの間、目標後輪タイヤ縦力応答Fxrtarを、前輪路面摩擦係数μfの変化を検出した時間の前の後輪タイヤ縦力検出値Fxrsensから、変化後の後輪タイヤ縦力許容値Fxrtolへ滑らかに変化させるようにした。制御を開始してから、前後輪遅延時間tdelayのうちに後輪4RLRRがタイヤ限界を超えないような後輪タイヤ縦力Fxrを実現するように滑らかに後輪制駆動トルクTrを制御することで、制御時にドライバが感じる違和感を低減しつつ車両挙動を安定化することができる。[効果]実施例2の効果について以下に述べる9前輪路面摩擦係数μfの変化を検出した時刻と、後輪路面摩擦係数μrが変化すると予測される時刻との間で、目標後輪タイヤ縦力応答Fxtar後輪制駆動力制限値を、変化前の後輪タイヤ縦力検出値Fxrsens路面摩擦係数相当の制限値から変化後の後輪タイヤ縦力許容値Fxrtol路面摩擦係数相当の制限値へ滑らかに変化させるようにした。よって、制御時にドライバが感じる違和感を低減しつつ車両挙動を安定化することができる。実施例3実施例1では、車両制御器87の閾値比較器87gから制御開始信号が出力されると、制動器87hはブレーキ予圧を負荷するようにしていた。実施例3では、車両制御器87に制御開始判定器87mを設け、この制御開始判定器87mから制御開始信号が出力されると、滑らかに横加速度Gyが変化するように操舵角指令値δdを出力するようにした。実施例3は、車両制御器87の構成が実施例1と異なる。以下、実施例3について説明するが、実施例1と同じ構成については同一の符号を付して説明を省略する。[車両制御器の構成]図15は車両制御器87の制御ブロック図である。図15に示すように車両制御器87は、横加速度検出器87rと、制御開始判定器87sと、目標横加速度応答生成器87tと、横加速度偏差算出器87uと、操舵角指令値算出器87vとを有している。横加速度検出器横加速度検出器87rは横加速度Gyを検出する。ここで、横加速度検出値Gysensは、例えば既存の横滑り防止装置に内蔵されたセンサで計測された量を用いればよい。制御開始判定器制御開始判定器87sは、前輪路面摩擦係数推定値μfestと、横加速度検出値Gysensとを入力し、横加速度制御を開始するか否か判断する。そして、開始するべきであると判断したときには制御開始信号を出力する。例えば、次の条件を満たすときに制御開始信号を出力する。条件1:所定時間内に前輪路面摩擦係数推定値μfestが所定値より大きく変動すること。条件2:現在の横加速度検出値Gysensが前輪路面摩擦係数推定値μfestと重力加速度gで与えられる加速度よりも大きいこと。以上の条件1および条件2を満たすときに制御開始信号を出力する。以下の目標横加速度応答生成器87t、横加速度偏差算出器87u、操舵角指令値算出器87vにおける演算は制御開始信号が出力されたときに実行する。目標横加速度応答生成器目標横加速度応答生成器87tは、前輪路面摩擦係数推定値μfestと、横加速度検出値Gysensと、前後輪遅延時間tdelayとを入力し、目標横加速度応答Gytarを生成する。まず、前輪路面摩擦係数推定値μfestから算出される許容される横加速度許容値は、前輪路面摩擦係数推定値μfestと重力加速度gとの積であたえられる。制御の結果、前後輪遅延時間tdelay後に実現したい横加速度である目標最終横加速度Gyfinは、この横加速度許容値よりも小さい値でなければならない。つまり、目標最終横加速度は、次の式で求める値とすれば良い。ここで、Sはマージンを設定する係数であり、0より大きく1より小さい値に設定する。図16は、目標横加速度応答Gytarのマップである。目標横加速度応答Gytarは、例えば、図16に示すように、現在時刻と現在時刻における横加速度検出値Gysensとでプロットされる図中A点と、現在時刻から前後輪遅延時間後の時刻と目標最終横加速度Gyfinとでプロットされる図中B点とを結ぶ線分を目標横加速度応答Gytarとすれば良い。横加速度偏差算出器横加速度偏差算出器87uは、横加速度検出値Gysensと目標横加速度応答Gytarとの偏差を算出し、横加速度偏差eGyとして出力する。操舵角指令値算出器操舵角指令値算出器87vは、横加速度偏差eGyを小さくするように操舵角指令値δdを算出する。操舵角指令値δdの算出には、例えばPI制御器を用いる。すなわち、入出力関係を伝達関数で表現すると、次の式に示すようになる。ここで、Kpは比例ゲイン、Kiは積分ゲインである。[作用]実施例3では、前輪路面摩擦係数μfの変化を検出した時間から前後輪遅延時間tdelayの間、目標横加速度応答Gytarを、前輪路面摩擦係数μfの変化を検出した時間の前の横加速度検出値Gysensから、変化後の目標最終横加速度Gyfinへ滑らかに変化させるようにした。制御を開始してから、前後輪遅延時間tdelayのうちに滑らかに目標横加速度を制御することで、制御時にドライバが感じる違和感を低減しつつ、後輪路面摩擦係数が変化したときのときの不安定挙動を予め防止することができる。[効果]実施例3の効果について以下に述べる。10前輪路面摩擦係数μfの変化を検出した時刻と、後輪路面摩擦係数μrが変化すると予測される時刻との間で、目標横加速度応答Gytarを、変化前の横加速度検出値Gytarから変化後の目標最終横加速度Gyfinへ滑らかに変化させるようにした。よって、制御時にドライバが感じる違和感を低減しつつ車両挙動を安定化することができる。実施例4実施例1では、前輪路面摩擦係数推定値μfestを求め、その推定値から後輪路面摩擦係数推定値μrestを求めるようにしていた。実施例4では、左右前輪独立して、左前輪路面摩擦係数推定値μflestおよび右前輪路面摩擦係数推定値μfrestを求め、その推定値から左後輪路面摩擦係数推定値μrlestおよび右後輪路面摩擦係数推定値μrrestを求めるようにした。実施例4の統合コントローラ8の制御ブロック自体の構成は実施例1と同様であり状態推定器86の各構成の処理の内容が一部異なる。よって、実施例1の図4を用いて説明する。なお以下では、実施例1と処理が異なる構成のみ説明する。[状態推定器]状態推定器86は、図4に示すように状態推定器86は、車両ダイナミクスモデル86aと、前輪路面摩擦係数ダイナミクスモデル86bと、後輪路面摩擦係数モデル86cと、セルフアライニングトルク算出器86dと、偏差算出器86eと、補正量算出器86fと、路面μ補正ゲイン増加判定器86gと、車両状態補正ゲイン判定器86hとを有している。前輪路面摩擦係数ダイナミクスモデル前輪路面摩擦係数ダイナミクスモデル86bは、左前輪路面摩擦係数補正量Cμflと右前輪路面摩擦係数補正量Cμfrとを入力し、左前輪路面摩擦係数推定値μflestと右前輪路面摩擦係数推定値μfrestとを求める。左前輪路面摩擦係数推定値μflestは次の式で、右前輪路面摩擦係数推定値μfrestは次の式で求められる。左前輪路面摩擦係数補正量Cμflは左前輪路面摩擦係数に関する補正量であり、右前輪路面摩擦係数補正量Cμfrは右前輪路面摩擦係数に関する補正量であり、補正量算出器86fにより求められる。後輪路面摩擦係数モデル後輪路面摩擦係数モデル86cは、左前輪路面摩擦係数推定値μflestと、右前輪路面摩擦係数推定値μfrestと、前後輪遅延時間tdelayとを入力し、左後輪路面摩擦係数μrestと、右後輪路面摩擦係数μrrestとを求める。まず、次の式より前後輪遅延時間tdelayから遅延ステップ数Ndelayを求める。tdelayは、前後輪遅延時間算出器82で求められた前後輪遅延時間である。floorは引数が小数の場合は、小数点以下を切り捨てて整数に変換する演算子を表わす。左後輪路面摩擦係数推定値μrlestは、現在の左前輪路面摩擦係数μflestに対して遅延ステップ数Ndelay前の左前輪路面摩擦係数μflestの値に設定される。右後輪路面摩擦係数推定値μrrestは、現在の右前輪路面摩擦係数μfrestに対して遅延ステップ数Ndelay前の右前輪路面摩擦係数μfrestの値に設定される。セルフアライニングトルク算出器セルフアライニングトルク算出器86dは、左前輪路面摩擦係数推定値μflestと、右前輪路面摩擦係数推定値μfrestと、車体すべり角推定値βestと、ヨーレート推定値γestを入力し、セルフアライニングトルク推定値Mzestを求める。セルフアライニングトルク推定値Mzestは、例えば実施例1に記載のフィアラモデルを用いて、左前輪セルフアライニングトルクMzl、右前輪セルフアライニングトルクMzrを算出する。偏差算出器偏差算出器86eは、車両ダイナミクスモデル86aで求められたヨーレート推定値γestと、車両状態検出器84で検出されたヨーレート検出値γsensとの偏差eγ、およびセルフアライニングトルク算出器86dで求められた左前輪セルフアライニングトルク推定値Mzlestと、前輪状態検出器83で検出された左前輪セルフアライニングトルク検出値Mzlsensとの偏差eMzl、およびセルフアライニングトルク算出器86dで求められた右前輪セルフアライニングトルク推定値Mzrestと、前輪状態検出器83で検出された右前輪セルフアライニングトルク検出値Mzrsensとの偏差eMzrを求める。ヨーレート推定値γestとヨーレート検出値γsensとの偏差eγは次の式により求められる。左前輪セルフアライニングトルク推定値Mzlesと左前輪セルフアライニングトルク検出値Mzlsensとの偏差eMzlは次の式により求められる。左前輪セルフアライニングトルク推定値Mzresと左前輪セルフアライニングトルク検出値Mzrsensとの偏差eMzrは次の式により求められる。補正量算出器補正量算出器86fは、偏差算出器86eで求めた偏差eγeMzleMzrを入力し、車体すべり角補正量Cβと、ヨーレート補正量Cγと、左前輪路面摩擦係数補正量Cμflと、右前輪路面摩擦係数補正量Cμfrとを求める。車体すべり角補正量Cβ、ヨーレート補正量Cγ、左前輪路面摩擦係数補正量Cμfl、右前輪路面摩擦係数補正量Cμfrは、それぞれ次の式~により求められる。補正ゲインKcβ1Kcβ2Kcβ3Kcγ1Kcγ2Kcγ3Kcμf1Kcμf2Kcμf3は、設計者により適当に定められる。[作用]実施例4では、左前輪路面摩擦係数推定値μflestと右前輪路面摩擦係数推定値μfrestとを独立して推定するようにした。そして、左後輪路面摩擦係数推定値μrlestを、現在の左前輪路面摩擦係数推定値μflestに対して前後輪遅延時間tdelay前の左前輪路面摩擦係数推定値μflestに設定し、右後輪路面摩擦係数推定値μrrestを、現在の右前輪路面摩擦係数推定値μfrestに対して前後輪遅延時間tdelay前の右前輪路面摩擦係数推定値μfrestに設定するようにした。これにより、左右輪で路面摩擦係数μが異なる路面を跨ぐ所謂スプリットμ時に、より精度良く、車体すべり角βおよび路面摩擦係数を推定することができる。[効果]11状態推定器86は、前輪4FL4FRを左右独立に路面摩擦係数を推定し、左後輪路面摩擦係数推定値μrlestを、遅延時間tdelay分遅れた左前輪路面摩擦係数推定値μflestとし、右後輪路面摩擦係数推定値μrrestを、遅延時間tdelay分遅れた右前輪路面摩擦係数推定値μfrestとするようにした。よって、左右輪で路面摩擦係数μが異なる路面を跨ぐ所謂スプリットμ時に、より精度良く、車体すべり角βおよび路面摩擦係数を推定することができる。他の実施例以上、本発明を実施するための最良の形態を、図面に基づく実施例により説明したが、本発明の具体的な構成は、実施例に限定されるものではなく、発明の要旨を逸脱しない範囲の設計変更等があっても本発明に含まれる。例えば、実施例1では制動器87hをブレーキ予圧しているが、車両挙動制御に関連する他の機器のブレーキやクラッチのクリアランスを予め小さくするものであっても良い。