(57)【特許請求の範囲】
【請求項１】
車輪速度信号を検出する検出手段と、
前記検出手段により検出された車輪速度信号の内、少なくとも１つの共振点又は少なくとも１つの反共振点を含む周波数が低周波数領域より大きい所定範囲の車輪速度信号を抽出する抽出手段と、
前記抽出手段により抽出された車輪速度信号に基づいて、車輪の振動の大きさを表す物理量を演算する演算手段と、
前記検出手段により検出された車輪速度信号に基づいて、車輪のすべり易さを表す物理量として、路面μ勾配、制動トルク勾配または駆動トルク勾配のいずれか１つを推定する推定手段と、
前記演算手段により演算された車輪の振動の大きさを表す物理量と前記推定手段により推定された車輪のすべり易さを表す物理量とに基づいて、車輪の路面に対する摩擦状態を推定する路面摩擦状態推定手段と、
を備えた路面摩擦状態推定装置。
【請求項２】
前記車輪の振動の大きさを表す物理量は、前記抽出手段により抽出された車輪速度信号の振動レベル又は振幅であることを特徴とする請求項１記載の路面摩擦状態推定装置。
【請求項３】
車輪速度信号を検出する検出手段と、
前記検出手段により検出された車輪速度信号に基づいて、車輪のすべり易さを表す物理量として、路面μ勾配、制動トルク勾配または駆動トルク勾配のいずれか１つを推定する推定手段と、
前記推定手段により推定された車輪のすべり易さを表す物理量の時間的変化率を演算する演算手段と、
前記推定手段により推定された車輪のすべり易さを表す物理量と前記演算手段により演算された車輪のすべり易さを表す物理量の時間的変化率とに基づいて、車輪の路面に対する摩擦状態を推定する路面摩擦状態推定手段と、
を備えた路面摩擦状態推定装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、路面摩擦状態推定装置に係り、より詳しくは、車輪と路面との間の摩擦係数の滑り速度に対する勾配である路面μ勾配などの車輪の滑り易さを表す物理量に基づいて、車輪の路面に対する摩擦状態を推定する路面摩擦状態推定装置に関する。
【０００２】
【従来の技術及び発明が解決しようとする課題】
従来より、車輪速度信号を検出し、検出した車輪速度信号に基づいて路面μ勾配を推定する装置が提案されている。この装置では、路面μ勾配を推定するために車輪速度信号を検出しているため、突起の乗り越し時や、圧雪などの荒れた路面では、路面μ勾配の推定値が小さくなる。よって、路面μ勾配の推定精度が悪い。
【０００３】
例えば、図１０に示すように、アスファルト路を走行中に突起を乗り越した場合には、前述したように、車輪速度信号が小さくなり、この時に推定された路面μ勾配の推定値が小さくなり、低μ路における路面μ勾配の推定値とほぼ同様な値となり、突起乗り越し時に低μ路を走行していると判定する場合がある。
【０００４】
本発明は、上記事実に鑑み成されたもので、車輪の路面に対する摩擦状態を精度良く推定することの可能な路面摩擦状態推定装置を提供することを目的とする。
【０００５】
【課題を解決するための手段】
上記目的達成のため請求項１記載の発明は、車輪速度信号を検出する検出手段と、前記検出手段により検出された車輪速度信号の内、少なくとも１つの共振点又は少なくとも１つの反共振点を含む周波数が低周波数領域より大きい所定範囲の車輪速度信号を抽出する抽出手段と、前記抽出手段により抽出された車輪速度信号に基づいて、車輪の振動の大きさを表す物理量を演算する演算手段と、前記検出手段により検出された車輪速度信号に基づいて、車輪のすべり易さを表す物理量として、路面μ勾配、制動トルク勾配または駆動トルク勾配のいずれか１つを推定する推定手段と、前記演算手段により演算された車輪の振動の大きさを表す物理量と前記推定手段により推定された車輪のすべり易さを表す物理量とに基づいて、車輪の路面に対する摩擦状態を推定する路面摩擦状態推定手段と、を備えている。
【０００６】
即ち、本発明の検出手段は、車輪速度信号を検出し、抽出手段は、検出手段により検出された車輪速度信号の内、少なくとも１つの共振点又は少なくとも１つの反共振点を含む周波数が低周波数領域より大きい所定範囲の車輪速度信号を抽出する。
【０００７】
算出手段は、抽出手段により抽出された車輪速度信号に基づいて、車輪の振動の大きさを表す物理量を演算する。なお、車輪の振動の大きさを表す物理量は、抽出手段により抽出された車輪速度信号の振動レベル又は振幅としてもよい。
【０００８】
推定手段は、検出手段により検出された車輪速度信号に基づいて、車輪のすべり易さを表す物理量を推定する。ここで、車輪の滑り易さを表す物理量は、スリップ速度に対する車輪と路面との間の摩擦係数μの勾配である路面μ勾配や、スリップ速度に対する制動トルクの勾配である制動トルク勾配、及び駆動トルクの勾配である駆動トルク勾配のいずれか１つである。
【０００９】
ところで、車輪の振動の大きさを表す物理量と車輪の滑り易さを表す物理量との関係は路面状態毎に定まる。
【００１０】
そこで、推定手段は、演算手段により演算された車輪の振動の大きさを表す物理量と推定手段により推定された車輪のすべり易さを表す物理量とに基づいて、車輪の路面に対する摩擦状態を推定する。
【００１１】
このように、車輪の振動の大きさを表す物理量と車輪のすべり易さを表す物理量とに基づいて、車輪の路面に対する摩擦状態を推定するので、車輪の路面に対する摩擦状態を精度よく判定することができる。
【００１２】
ここで、車輪の振動の大きさを表す物理量は、車輪のすべり易さを表す物理量の逆数に比例して大きくなる。即ち、逆転現象が発生する。このため、低周波数領域の車輪速度信号を抽出して求めた車輪の振動の大きさを表す物理量及び車輪のすべり易さを表す物理量に基づいて、車輪の路面に対する摩擦状態を推定しても、車輪の路面に対する摩擦状態を精度よく判定することができない。
【００１３】
そこで、本発明では、上記摩擦状態を推定するために、低周波数領域より大きい所定範囲の車輪速度信号を抽出して、車輪の路面に対する摩擦状態を精度よく判定している。
【００１４】
請求項３記載の発明は、車輪速度信号を検出する検出手段と、前記検出手段により検出された車輪速度信号に基づいて、車輪のすべり易さを表す物理量として、路面μ勾配、制動トルク勾配または駆動トルク勾配のいずれか１つを推定する推定手段と、前記推定手段により推定された車輪のすべり易さを表す物理量の時間的変化率を演算する演算手段と、前記推定手段により推定された車輪のすべり易さを表す物理量と前記演算手段により演算された車輪のすべり易さを表す物理量の時間的変化率とに基づいて、車輪の路面に対する摩擦状態を推定する路面摩擦状態推定手段と、を備えている。
【００１５】
上記車輪の滑り易さを表す物理量と車輪のすべり易さを表す物理量の変化率との関係は路面状態毎に定まる。
【００１６】
そこで、路面摩擦状態推定手段は、推定手段により推定された車輪のすべり易さを表す物理量と前記演算手段により演算された車輪のすべり易さを表す物理量の変化率とに基づいて、車輪の路面に対する摩擦状態を推定する。
【００１７】
このように、車輪のすべり易さを表す物理量と車輪のすべり易さを表す物理量の変化率とに基づいて、車輪の路面に対する摩擦状態を推定するので、車輪の路面に対する摩擦状態を精度よく判定することができる。
【００１８】
ここで、本発明においても、上記摩擦状態を推定するために、低周波数領域より大きい所定範囲の車輪速度信号を抽出しているため、車輪の路面に対する摩擦状態を精度よく判定するこたができる。
【００１９】
【発明の実施の形態】
以下、図面を参照して本発明の第１の実施の形態を説明する。
【００２０】
図１には、第１の実施の形態に係る路面摩擦推定装置の構成が示されている。同図に示すように、本路面摩擦推定装置は、各車輪の車輪速度信号ω1を検出する車輪速検出手段１と、検出された各車輪の車輪速度信号ω1から路面外乱ΔＴdを受けた車輪共振系の応答出力としての各車輪の車輪速振動Δω1を検出する前処理フィルタ２と、検出された車輪速振動Δω1を満足するような各車輪の伝達関数を最小自乗法を用いて同定する伝達関数同定手段３と、同定された伝達関数に基づいてタイヤと路面との間の摩擦係数μの勾配を各車輪毎に演算するμ勾配演算手段４と、を備えている。
【００２１】
図１において、車輪速検出手段１は、車輪速度に応じたセンサ出力信号を出力するいわゆる車輪速センサと、該センサ出力信号から各車輪の実際の回転速度信号を演算する演算手段と、から構成することができる。
【００２２】
また、前処理フィルタ２は、本車輪共振系の共振周波数と予想される周波数を中心として一定の帯域の周波数成分のみを通過させるバンドパスフィルタや、該共振周波数成分を含む高帯域の周波数成分のみを通過させるハイパスフィルタなどで構成することができる。なお、本実施の形態では、このバンドパスフィルタ或いはハイパスフィルタの周波数特性を規定するパラメータを一定値に固定したものである。
【００２３】
なお、この前処理フィルタ２の出力は、直流成分を除去したものとする。すなわち、車輪速度信号ω1の回りの車輪速振動Δω1のみが抽出される。
【００２４】
いまここで、前処理フィルタ２の伝達関数Ｆ（ｓ）を、
【００２５】
【数１】
0003769435.tif 000002
【００２６】
とする。ただし、ｃiはフィルタ伝達関数の係数、ｓはラプラス演算子である。
【００２７】
次に、伝達関数同定手段３が依拠する演算式を導出しておく。なお、本実施の形態では、前処理フィルタ２の演算を、伝達関数同定手段３の演算に含めて実施する。
【００２８】
まず、第１の実施の形態で同定すべき伝達関数は、路面外乱ΔＴdを加振入力として、このとき前処理フィルタ２により検出された車輪速振動Δω1を応答出力とする２次のモデルとする。すなわち、
【００２９】
【数２】
0003769435.tif 000003
【００３０】
の振動モデルを仮定する。ここに、ｖは車輪速信号を観測するときに含まれる観測雑音である。（２）式を変形すると、次式を得る。
【００３１】
【数３】
0003769435.tif 000004
【００３２】
まず、（３）式に（１）式の前処理フィルタを掛けて得られた式を離散化する。このとき、Δω1、ΔＴd、ｖは、サンプリング周期Ｔｓ  毎にサンプリングされた離散化データΔω1  （ｋ）、ΔＴｄ  （ｋ）、ｖ（ｋ）（ｋはサンプリング番号：ｋ＝１，２，３，．．．．  ）として表される。また、ラプラス演算子ｓは、所定の離散化手法を用いて離散化することができる。本実施の形態では、１例として、次の双一次変換により離散化するものとする。なお、ｄは１サンプル遅延演算子である。
【００３３】
【数４】
0003769435.tif 000005
【００３４】
また、前処理フィルタの次数ｍは、２以上が望ましいので、本実施の形態では、演算時間も考慮してｍ＝２とし、これによって次式を得る。
【００３５】
【数５】
0003769435.tif 000006
【００３６】
また、最小自乗法に基づいて、車輪速振動Δω1の各データから伝達関数を同定するために、（４）式を、同定すべきパラメータに関して一次関数の形式となるように、次式のように変形する。なお、”T”を行列の転置とする。
【００３７】
【数６】
0003769435.tif 000007
【００３８】
である。上式において、θが同定すべき伝達関数のパラメータとなる。
【００３９】
また、本実施の形態の路面摩擦状態推定装置は、車輪速検出手段によって検出された車輪速信号のうち、少なくとも１つの共振点または少なくとも１つの反共振点を含む周波数が低周波数領域より大きい所定範囲の車輪速度信号を抽出するバンドパスフィルタ５を備えている。バンドパスフィルタ５には、バンドパスフィルタ５より抽出された車輪速信号に基づいて、車輪の振動の大きさを表す物理量、本実施の形態では、振動レベルを演算する振動レベル演算手段６Ａが接続されている。そして、本実施の形態の路面摩擦状態推定装置は、μ勾配演算手段４により演算された路面μ勾配と、振動レベル演算手段６Ａにより演算された振動レベルとに基づいて車輪の路面に対する摩擦状態を推定する路面摩擦状態推定手段７を備えている。
【００４０】
次に、本実施の形態の作用を説明する。
【００４１】
伝達関数同定手段３では、検出された車輪速振動Δω1の離散化データを(２２)式に順次当てはめた各データに対し、最小自乗法を適用することによって、未知パラメータθを推定し、これにより伝達関数を同定する。
【００４２】
具体的には、検出された車輪速振動Δω1を離散化データΔω（ｋ）（ｋ＝１，２，３，．．．）に変換し、該データをＮ点サンプルし、次式の最小自乗法の演算式を用いて、伝達関数のパラメータθを推定する。
【００４３】
【数７】
0003769435.tif 000008
【００４４】
ここに、記号”＾”の冠した量をその推定値と定義することにする。
【００４５】
また、上記最小自乗法は、次の漸化式によってパラメータθを求める逐次型最小自乗法として演算してもよい。
【００４６】
【数８】
0003769435.tif 000009
【００４７】
ここに、ρは、いわゆる忘却係数で、通常は０．９５〜０．９９の値に設定する。このとき、初期値は、
【００４８】
【数９】
0003769435.tif 000010
【００４９】
とすればよい。
【００５０】
また、上記最小自乗法の推定誤差を低減する方法として、種々の修正最小自乗法を用いてもよい。本実施の形態では、補助変数を導入した最小自乗法である補助変数法を用いた例を説明する。該方法によれば、（９）式の関係が得られた段階でｍ（ｋ）を補助変数として、次式を用いて伝達関数のパラメータを推定する。
【００５１】
【数１０】
0003769435.tif 000011
【００５２】
また、逐次演算は、以下のようになる。
【００５３】
【数１１】
0003769435.tif 000012
【００５４】
補助変数法の原理は、以下の通りである。（２８）式に（２２）式を代入すると、
【００５５】
【数１２】
0003769435.tif 000013
【００５６】
となるので、（１９）式の右辺第２項が零となるように補助変数を選べばθの推定値は、θの真値に一致する。そこで、本実施の形態では、補助変数として、ζ（ｋ）＝［−ξy1（ｋ）−ξy2（ｋ）］Tを式誤差ｒ（ｋ）と相関を持たないほどに遅らせたものを利用する。すなわち、
ｍ（ｋ）＝［−ξy1（ｋ−Ｌ）−ξy2（ｋ−Ｌ）］T        （２０）
とする。ただし、Ｌは遅延時間である。
【００５７】
上記のようにして伝達関数を同定した後、μ勾配演算手段４において、路面μ勾配Ｄ０  に関係する物理量を、
【００５８】
【数１３】
0003769435.tif 000014
【００５９】
と演算する。このように（２１）式により路面μ勾配Ｄ０  に関係する物理量を演算できると、例えば、該物理量が小さいとき、タイヤと路面との間の摩擦特性が飽和状態であると容易に判定できる。
【００６０】
ところで、車輪速検出手段１により、検出された車輪速信号は、周波数に応じて分解すると、図２に示すように、２つの共振点と１つの反共振点を有する。２つの共振点のうち周波数が小さい側の共振点は、タイヤ慣性などに基づく前後共振点であり、周波数はｆ1（１５〜２０）Ｈｚである。また、２つの共振点のうち周波数が高い側の共振点は、タイヤの空気圧やタイヤゴム弾性に基づく捩れ共振点であり、周波数はｆ3（３５〜４０）Ｈｚである。そして、車輪速信号には、色々な信号に対して不感帯の反共振点を有し、周波数ｆ2（２０〜２５）Ｈｚである。本実施の形態に係るバンドパスフィルタ５は、車輪速検出手段１により検出された車輪速信号のうち、捩れ共振点（周波数ｆ3）含む所定範囲Δｆの車輪速信号を抽出している。なお、バンドパスフィルタ５は、捩れ共振点以外の、前後共振点や反共振点を含む所定範囲Δｆの車輪速信号を抽出するようにしてもよい。
【００６１】
振動レベル演算手段６Ａは次式より定まる車輪速信号の振動レベルＧ（Ｎ）を演算する。なお、バンドパスフィルタ５の出力はω（ｋ）である。
【００６２】
【数１４】
0003769435.tif 000015
【００６３】
なお、振動レベル演算手段６Ａは、実際には、演算タイミング毎に、
【００６４】
【数１５】
0003769435.tif 000016
【００６５】
の漸化式を逐次的に演算する。このように、振動レベル演算手段６Ａで振動レベルを演算するのは次の理由からである。即ち、前述したように路面μ勾配は、突起乗り越し等により、適正に推定されない（図１０参照）。
【００６６】
一方、このようにして得られた振動レベルと路面μ勾配の推定値とを、アスファルト路面、低μ路路面においてプロットすると、図３に示すようにアスファルト路面と低μ路路面とでは明確に区別されて認識されるとともに、アスファルト路面中に突起を乗り越したとしても、低μ路における領域と明確に区別されて認識される。上記事実に鑑み本発明者は、種々の路面状態毎に振動レベルと路面μ勾配との推定値との関係を求めたところ、図４に示すように、振動レベルと路面μ勾配の推定値とは各路面状態、例えば低μ路、アスファルト路、石畳路、圧雪路、砂利・非整地路毎に領域毎に明確に区別されて認識出来ることが実験の結果得られた。即ち、例えば、圧雪路は路面μ勾配の推定値としては高μ路（アスファルト路、石畳路）より若干小さい値となっているが、振動レベルの違いによって領域に差が表れている。砂利道等は更に振動レベルが高い領域にあるが、路面μ勾配の推定値は高μ路より低い値となっている。
【００６７】
そこで、路面摩擦状態推定手段７は、振動レベル演算手段６Ａにより演算された振動レベル及びμ勾配演算手段４により演算された路面μ勾配の推定値と、上記路面状態毎の両者の関係（図４参照）と、基づいて、車輪の路面に対する摩擦状態を推定する。
【００６８】
ところで、本実施の形態では、バンドパスフィルタ５は、車輪速信号のうち、２つの共振点の何れか一方の共振点の車輪速信号または、反共振点の車輪速信号を抽出する。即ち、図２に示すように、周波数が低周波数領域より大きい車輪速信号を検出するようにしている。このように周波数が低周波数領域より大きい車輪速信号を検出するようにしているのは、図２に示すように周波数が低周波数領域より小さい車輪速信号は誤差が大きいため影響が大きくなり、路面摩擦状態を精度よく推定することが出来ないためである。
【００６９】
即ち、図５には、アスファルト路、低μ路を走行したときの車輪速度信号の周波数特性が示されている。なお、低μ路は、アスファルト路に比較して凹凸が非常に小さくなっている。
【００７０】
４０Ｈｚ付近にある共振特性を見ると、低μ路はアスファルト路に比較して振動レベルが小さくなっており、路面の凹凸をよく反映していることが理解できる。また、共振の強さも小さくなっており、本手法によって低μ路であることが推定できる。
【００７１】
低周波数領域(５Ｈｚ以下）の振動成分をみると、この領域では、振動レベルは路面μ勾配の逆数に比例して大きくなるため、低μ路の低周波数領域の車輪速度信号の振動レベルは、アスファルト路の低周波数領域の車輪速度信号の振動レベルより大きくなっている。
【００７２】
よって、低周波数領域の車輪速度信号に基づいて振動レベル及び路面μ勾配を演算して路面状態を推定すると、次のような問題が生ずる。
【００７３】
即ち、低周波数領域の振動成分は、路面μ勾配に依存し、実路面の凹凸状態を反映していない。このため、路面判定用のマップ作成には多くの適合工数を必要となる。
【００７４】
また、低周波数領域の車輪速度信号には、車両のロール、ピッチング運動による振動成分も含まれているため、走行条件によっては路面状態を精度よく判定することができない場合がある。
【００７５】
更に、上記理由から、判定マップは、タイヤや車両の物理パラメータに応じて用意しなければならないが、タイヤ交換等により車両の緒元を変更したときは対応が困難であり、場合によっては路面状態の判定ができない虞がある。
【００７６】
一方、本実施の形態では、低周波数領域より大きい領域の車輪速度信号を対象としているため、上記問題が発生せず、路面状態を精度よく判定することができる。
【００７７】
次に本発明の第２の実施の形態を説明する。本実施の形態の構成は、前述した第１の実施の形態と同様の構成部分を有するので、同一部分は同一の符合を付してその説明を省略し、異なる部分のみを説明する。
【００７８】
即ち、本実施の形態に係る路面摩擦状態推定装置は、第１の実施の形態に係る振動レベル演算手段６Ａにかえて、振幅演算手段６Ｂを備えている。
【００７９】
次に、本実施の形態に係る作用を説明する。なお、本実施の形態に係る作用は前述した第１の実施の形態に係る作用と同様の作用部分があるので、同一の作用部分はその説明を省略し、異なる作用部分についてのみ説明する。
【００８０】
本実施の形態に係る振幅演算手段６Ｂは、バンドパスフィルタ５により抽出された所定周波数の車輪速信号の振幅を演算する。具体的には、次の通りである。
【００８１】
即ち、所定周波数ｆの振幅を求めるには、後述するｓｉｎ関数の値Ｓｉとｃｏｓ関数の値Ｃｉを用いて求める。
【００８２】
即ち、サンプル時間をＴＳとし、
【００８３】
【数１６】
0003769435.tif 000017
【００８４】
を求め、振幅を求めたい車輪速信号をｙ、現在のサンプリング時の車輪速信号の値をｙ（Ｎ）とすると、
【００８５】
【数１７】
0003769435.tif 000018
【００８６】
となり、振幅Ａが求められる。
【００８７】
ところで、本発明者は、上記のように演算された振幅と路面μ勾配の推定値との関係は、路面状態毎に領域に分かれて明確に区別されて認識出来ることを発見した。
そこで、本実施の形態に係る路面摩擦状態推定手段７は、μ勾配演算手段４により演算された路面μ勾配及び振幅演算手段６Ｂにより演算された振幅と上記関係（図７参照）とに基づいて、路面摩擦状態を推定する。
【００８８】
次に、本発明の第３の実施の形態を説明する。本実施の形態の構成は、前述した第１の実施の形態の構成と同一部分を有するので、同一部分には同一の符号を付してその説明を省略し、異なる部分のみを説明する。
【００８９】
即ち、本実施の形態に係る路面状態推定装置は、車輪速検出手段１、前処理フィルタ２、伝達関数道程手段３、及びμ勾配演算手段４を備えている。更に本実施の形態に係る路面摩擦状態推定装置は、μ勾配演算手段４により演算された路面μ勾配の変化率を演算するμ勾配変化率演算手段６Ｃを備えている。
【００９０】
本発明者は、路面μ勾配の変化率と路面μ勾配の推定値とが、図９に示すように、路面状態に応じて領域毎に分割されて明確に区別されて認識出来ることを発見した。
【００９１】
そこで、本実施の形態に係る路面摩擦状態推定手段７は、μ勾配演算手段４により演算された路面μ勾配及びμ勾配変化率演算手段６Ｃにより演算された路面μ勾配の変化率と、上記関係（図９参照）と、に基づいて車輪の路面に対する摩擦状態を推定する。
【００９２】
前述した第１の実施の形態乃至第３の実施の形態では、路面μ勾配を演算しているが、本発明はこれに限定されるものでなく、路面μ勾配と等価な制動トルク勾配を演算するようにしても良い。
【００９３】
【発明の効果】
以上説明したように本発明は、車輪の路面に対する摩擦状態を推定するために、低周波数領域より大きい所定範囲の車輪速度信号を抽出するので、車輪の路面に対する摩擦状態を精度よく判定することができる、という効果を有する。
【図面の簡単な説明】
【図１】第２の実施の形態のブロック図である。
【図２】車輪速信号の周波数と振幅との関係を示したグラフである。
【図３】車輪速度信号の振動レベルと路面μ勾配との、アスファルト路及び低μ路の関係を示したグラフである。
【図４】車輪速度信号の振動レベルと路面μ勾配との路面状態毎の関係を示したグラフである。
【図５】車輪速度信号の周波数特性を示した図である。
【図６】第２の実施の形態のブロック図である。
【図７】車輪速度信号の振幅と路面μ勾配との路面状態毎の関係を示したグラフである。
【図８】第３の実施の形態のブロック図である。
【図９】路面μ勾配と路面μ勾配の変化率との路面状態毎の関係を示したグラフである。
【図１０】低μ路を走行する前に突起を乗り越した時の路面μ勾配推定値の変化を表すグラフである。
【符号の説明】
１      車輪速検出手段(検出手段）
４      μ勾配演算手段(推定手段）
５      バンドパスフィルタ(抽出手段）
６Ａ    振動レベル演算手段（演算手段）
７      路面摩擦状態推定手段
６Ｂ    振幅演算手段（演算手段）
６Ｃ    μ勾配変化率演算手段（演算手段）
【図１】
0003769435.tif 000019
【図２】
0003769435.tif 000020
【図３】
0003769435.tif 000021
【図４】
0003769435.tif 000022
【図５】
0003769435.tif 000023
【図６】
0003769435.tif 000024
【図７】
0003769435.tif 000025
【図８】
0003769435.tif 000026
【図９】
0003769435.tif 000027
【図１０】
0003769435.tif 000028