(57)【要約】【課題】 路面状態を精度よく推定する。【解決手段】 車両が走行することによりタイヤから発生するタイヤ発生音を、タイヤ発生音検出センサ１０が検出する。車輪速を車輪速検出センサ１２が検出する。前処理器１４は、検出されたタイヤ発生音から、路面状態及びタイヤ発生音の発生源に対応する特徴量を計算する。特徴ベクトル生成器１６は、前処理器１４での特徴量と、検出された車輪速と、を成分とする特徴ベクトルを生成する。判定処理器２０は、特徴ベクトル生成器１６で生成された特徴ベクトルと、判定マップメモリ２２に記憶された・複数の路面状態の各々毎に対応して記憶された路面状態及びタイヤ発生音の発生源に対応する特徴ベクトルと、に基づいて、車両が走行する路面の状態を推定する。【選択図】 図１
2004168286.tif 000002
【特許請求の範囲】
【請求項１】
車両が走行することによりタイヤから発生するタイヤ発生音を検出するタイヤ発生音検出手段と、
前記タイヤ発生音検出手段により検出されたタイヤ発生音から、路面状態及びタイヤ発生音の発生源に対応する少なくとも１つの特徴量を計算する計算手段と、
複数の路面状態の各々毎に、路面状態及びタイヤ発生音の発生源に対応する特徴量を記憶する記憶手段と、
前記計算手段により計算された特徴量と、前記記憶手段に複数の路面状態の各々毎に対応して記憶された特徴量と、に基づいて、車両が走行する路面の状態を推定する推定手段と、
を備えた路面状態推定装置。
【請求項２】
車輪速を検出する車輪速検出手段を更に備え、
前記記憶手段は、複数の路面状態の各々及び複数の車輪速の各々毎に対応して特徴量を記憶し、
前記推定手段は、前記計算手段により計算された特徴量、前記車輪速検出手段により検出された車輪速、及び、前記記憶手段に複数の路面状態の各々及び複数の車輪速の各々毎に対応して記憶された路面状態に基づいて、車両が走行する路面の状態を推定する、
ことを特徴とする請求項１記載の路面状態推定装置。
【請求項３】
前記計算手段は、前記タイヤ発生音検出手段により検出されたタイヤ発生音の波形と少なくとも１つの予め定められた長さの基本波形との相関値を求めることにより、路面状態及びタイヤ発生音の発生源に対応する少なくとも１つの特徴量を計算することを特徴とする請求項１又は請求項２記載の路面状態推定装置。
【請求項４】
前記計算手段は、ウェーブレット変換を用いることを特徴とする請求項１乃至請求項３の何れか１項に記載の路面状態推定装置。
【請求項５】
前記計算手段は、前記タイヤ発生音検出手段により検出されたタイヤ発生音の波形と少なくとも１つのウェーブレットの基本波形との相関値を求めることを特徴とする請求項１乃至請求項４の何れか１項に記載の路面状態推定装置。
【請求項６】
前記基本波形は、各々振幅が第１の範囲内の値である２つの第１の領域と、該２つの第１の領域に挟まれ、振幅が該第１の範囲より大きくかつ該第１の範囲より広い第２の範囲内の値である第２の領域と、を備えたことを特徴とする請求項３乃至請求項５の何れか１項に記載の路面状態推定装置。
【請求項７】
前記計算手段は、前記タイヤ発生音検出手段により検出されたタイヤ発生音の波形と少なくとも１つの路面状態に対応する基本波形との相関値を求めることを特徴とする請求項３記載の路面状態推定装置。
【請求項８】
前記計算手段は、予め定められた原波形を路面状態に対応するように変形することにより、前記基本波形を設定することを特徴とする請求項７記載の路面状態推定装置。
【請求項９】
路面状態を表す特徴量を検出する検出手段と、
各々異なる路面状態を表す複数の特徴量に対応して路面とタイヤとの間の摩擦状態である路面摩擦状態を記憶する記憶手段と、
前記検出手段により検出された路面状態を表す特徴量と、前記記憶手段に各々異なる路面状態を表す複数の特徴量に対応して記憶された路面摩擦状態と、に基づいて、前記検出手段により前記特徴量が検出される対象となる路面の路面摩擦状態を推定する推定手段と、
前記検出手段により前記特徴量が検出される対象となる路面の予め求められた路面摩擦状態と、前記推定手段により推定された路面摩擦状態と、に基づいて、前記記憶手段に記憶された路面摩擦状態を修正する修正手段と、
を備えた路面摩擦状態推定装置。
【請求項１０】
前記検出手段は、前方路面の路面状態を表す特徴量を検出することを特徴とする請求項９に記載の路面摩擦状態推定装置。
【請求項１１】
前記検出手段は、現在走行している路面の路面状態を表す特徴量を検出することを特徴とする請求項９に記載の路面摩擦状態推定装置。
【請求項１２】
前記検出手段は、請求項１に記載の路面状態推定装置であることを特徴とする請求項１１に記載の路面摩擦状態推定装置。
【請求項１３】
現在走行している路面の路面状態を表す特徴量を検出する第１の検出手段と、
前方路面の路面状態を表す特徴量を検出する第２の検出手段と、
前記第２の検出手段により検出された前方路面の路面状態を表す特徴量を、前記第２の検出手段により前記特徴量が検出される対象となった路面の前記第１の検出手段により検出された前記特徴量に同期させる同期手段と、
各々異なる路面状態を表す複数の特徴量に対応して路面とタイヤとの間の摩擦状態である路面摩擦状態を記憶する記憶手段と、
前記第１の検出手段により検出された現在走行している路面状態を表す特徴量、前記同期手段により同期された前記第２の検出手段により検出された前方路面の路面状態を表す特徴量、及び前記記憶手段に各々異なる路面状態を表す複数の特徴量に対応して記憶された路面摩擦状態と、に基づいて、前記第１の検出手段により前記特徴量が検出される対象となる路面の路面摩擦状態を推定する推定手段と、
を備えた路面摩擦状態推定装置。
【請求項１４】
前記第１の検出手段により前記特徴量が検出される対象となる路面の予め求められた路面摩擦状態と、前記推定手段により推定された路面摩擦状態と、に基づいて、前記記憶手段に記憶された路面摩擦状態を修正する修正手段を備えた請求項１３記載の路面摩擦状態推定装置。
【請求項１５】
車両が走行予定の前方路面の輝度情報を検出する検出手段と、
前記検出手段により検出された前方路面の輝度情報と、路面状態を表す路面状態物理量を演算するためのパラメータと、に基づいて、路面状態物理量を演算する演算手段と、
を備えた路面状態物理量演算装置であって、
前記演算手段は、前記検出手段により検出された前方路面の輝度情報に基づいて前記パラメータを修正し、該修正したパラメータと、前記検出手段により検出された前方路面の輝度情報と、に基づいて、路面状態を表す路面状態物理量を演算することを特徴とする路面状態物理量演算装置。
【請求項１６】
前記検出手段は、前方路面側の領域を複数の領域に分割しかつ該分割した複数の領域各々の輝度情報を検出し、
前記パラメータは、前記検出された複数の領域各々の輝度情報を各路面状態毎に分類するための値であり、
前記演算手段は、前記検出された複数の領域各々の輝度情報を前記パラメータにより各路面状態毎に分類することにより、前記路面状態物理量を演算する、
ことを特徴とする請求項１５記載の路面状態物理量演算装置。
【請求項１７】
前記検出手段は、前方路面側の領域を複数の領域に分割しかつ該分割した複数の領域各々の輝度情報を検出し、
前記パラメータは、前記検出された複数の領域各々の輝度情報が得られる確率を最大にする値であり、
前記演算手段は、前記パラメータにより最大となった、前記検出された複数の領域各々の輝度情報が得られる確率を演算することにより、前記路面状態物理量を演算する、
ことを特徴とする請求項１５記載の路面状態物理量演算装置。
【請求項１８】
複数の路面状態の各々毎に、輝度情報に基づいて定まる各路面の路面状態を表す路面状態物理量を記憶する記憶手段と、
前方路面の輝度情報を検出する検出手段と、
前記検出手段により検出された前方路面の輝度情報と、前記記憶手段により複数の路面状態の各々毎に記憶された各路面の路面状態を表す路面状態物理量と、に基づいて、前方路面の路面状態を推定する推定手段と、
を備えた路面状態推定装置。
【請求項１９】
前記推定手段は、前記検出手段により検出された前方路面の輝度情報から路面状態物理量を求め、該求めた路面状態物理量と、前記記憶手段により複数の路面状態の各々毎に記憶された各路面の路面状態物理量と、に基づいて、前方路面の路面状態を推定する請求項１８記載の路面状態推定装置。
【請求項２０】
請求項１記載の路面状態推定装置と、
前記路面状態推定装置により推定された路面状態と、前記路面状態推定装置により路面状態が推定される際に前記検出手段により検出された輝度情報と、に基づいて、前記記憶手段に記憶されかつ該路面状態推定装置により推定された路面状態に対応する路面状態物理量を更新する更新手段と、
を更に備えた請求項１８記載の路面状態推定装置。
【請求項２１】
車両が現在走行している現在走行路面の路面状態を検出する現在路面状態検出手段と、
前方路面の路面状態を検出する前方路面状態検出手段と、
前記現在路面状態検出手段により検出された現在走行路面の路面状態と、前記前方路面状態検出手段により検出された前方路面の路面状態と、を比較し、何れか一方の路面状態が変化した場合に、該変化した路面状態を報知する報知手段と、
を備えた路面状態報知装置。
【請求項２２】
前記報知手段は、前記前方路面の路面状態が前記現在路面の路面状態から変化した後、前記現在路面の路面状態が変化した場合、該変化した現在路面の路面状態を報知する請求項２１記載の路面状態報知装置。
【請求項２３】
現在走行している現在走行路面の路面状態を検出する現在路面状態検出手段と、
前方路面の路面状態を検出する前方路面状態検出手段と、
前記前方路面状態検出手段により検出された前方路面の路面状態を、前記現在路面状態検出手段により現在路面の路面状態が検出される対象となった現在路面に同期させる同期手段と、
前記現在路面状態検出手段により検出された現在路面の路面状態と、前記同期手段により、前記現在路面状態検出手段により現在路面の路面状態が検出される対象となった現在路面に同期された前記前方路面状態検出手段により検出された前方路面の路面状態と、に基づいて、現在路面の路面状態を報知する報知手段と、
を備えた路面状態報知装置。
【請求項２４】
車両が走行する路面の路面状態を検出する路面状態検出手段と、
前記路面と前記車両のタイヤとの間の摩擦状態である路面摩擦状態を推定する路面摩擦状態推定手段と、
前記路面状態検出手段により検出された路面状態と、前記路面摩擦状態推定手段により推定された路面摩擦状態と、に基づいて、最終の路面摩擦状態を推定する最終路面摩擦状態推定手段と、
を備えた路面摩擦状態推定装置。
【請求項２５】
複数の路面状態の各々毎に、各路面状態に対応する路面摩擦状態の範囲を予め定めておき、
前記各路面状態に対応する路面摩擦状態の範囲、前記路面状態検出手段により検出された路面状態、及び前記路面摩擦状態推定手段により検出された路面摩擦状態に基づいて、前記路面摩擦状態検出手段により検出された路面摩擦状態が、前記路面状態検出手段により検出された路面状態に対応する路面摩擦状態の範囲内か否かを判断する判断手段を更に備え、
前記最終路面摩擦状態推定手段は、前記判断手段により、前記路面摩擦状態検出手段により検出された路面摩擦状態が、前記路面状態推定手段により検出された路面状態に対応する路面摩擦状態の範囲内であると判断された場合に、前記路面状態検出手段により検出された路面状態と、前記路面摩擦状態推定手段により検出されかつ該路面状態に対応する路面摩擦状態の範囲内であると判断された路面摩擦状態と、に基づいて、最終の路面摩擦状態を推定する、
ことを特徴とする請求項２４路面摩擦状態推定装置。
【請求項２６】
車両が走行することによりタイヤから発生するタイヤ発生音を検出するタイヤ発生音検出手段と、
車輪速を検出する車輪速検出手段と、
複数の路面状態の各々毎に、タイヤ発生音の発生源に対応する特徴量と車輪速との関係を記憶する記憶手段と、
前記タイヤ発生音検出手段により検出されたタイヤ発生音、前記車輪速検出手段により検出された車輪速、及び前記記憶手段に前記複数の路面状態の各々毎に記憶された関係に基づいて、路面状態を推定する推定手段と、
を備えた路面状態推定装置。
【請求項２７】
前記タイヤ発生音検出手段により検出された車両が走行することによりタイヤから発生するタイヤ発生音と、前記車輪速検出手段により検出された車輪速と、に基づいて、前記記憶手段に記憶された関係を更新する更新手段を更に備え、
前記推定手段は、前記タイヤ発生音検出手段により検出されたタイヤ発生音、前記車輪速検出手段により検出された車輪速、及び前記記憶手段に前記複数の路面状態の各々毎に記憶されかつ前記更新手段により更新された関係に基づいて、路面状態を推定する、
ことを特徴とする請求項２６記載の路面状態推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、路面状態推定装置、路面摩擦状態推定装置、路面状態物理量演算装置、及び路面状態報知装置にかかり、特に、車両が走行することによりタイヤから発生するタイヤ発生音や路面の輝度情報などから路面状態を推定する路面状態推定装置、及び、この路面状態推定装置に関連する路面摩擦状態推定装置、路面状態物理量演算装置、及び路面状態報知装置に関する。
【背景技術】
【０００２】
タイヤから発生するタイヤ発生音を用いて路面状態を推定する装置が提案されている（特許文献１、２）。これらの装置では、タイヤから発生するタイヤ発生音をＦＦＴ解析し、タイヤ発生音の各周波数成分を入力とし、ニューラルネットワークを利用して路面状態を推定する。
【０００３】
また、路面をカメラによって撮影して得られた画像の水平、垂直偏向画像のパワースペクトル分布等の路面状態を表す特徴量を検出し、各々異なる路面状態を表す複数の特徴量に対応して路面とタイヤとの間の摩擦係数である路面摩擦係数を記憶しておき、検出された路面状態を表す特徴量と、各々異なる路面状態を表す複数の特徴量に対応して記憶された路面摩擦係数と、に基づいて、路面の路面摩擦係数を推定する装置が提案されている（特許文献３）。
【特許文献１】特開平６-１３８０１８号公報
【特許文献２】特開平７-１５６７８２号公報
【特許文献３】特開平８-３２７５３０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、上記装置では、タイヤから発生するタイヤ発生音をＦＦＴ解析し、タイヤ発生音の各周波数成分を入力とし、広範囲の周波数成分を入力するので、計算時間が多く必要となり、誤差も多くなる。これらの要因により、実際に車両に搭載しても、要求される水準を満たすことは困難である。
【０００５】
また、上記特許文献３の装置では、各々異なる路面状態を表す複数の特徴量に対応して路面とタイヤとの間の摩擦係数である路面摩擦係数を記憶しておき、検出された路面状態を表す特徴量と、各々異なる路面状態を表す複数の特徴量に対応して記憶された路面摩擦係数と、に基づいて、路面摩擦係数を推定するので、記憶した路面状態と路面摩擦係数の範囲内でしか、路面摩擦係数を推定することができない。よって、記憶した路面状態と路面摩擦係数が正しくない場合、精度よく路面摩擦係数を推定することができない。
【０００６】
本発明は、上記事実に鑑み成されたもので、路面状態を精度よく推定することの可能な路面状態推定装置、及び、この路面状態推定装置に関連する路面摩擦状態推定装置、路面状態物理量演算装置、及び路面状態報知装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記目的達成のため請求項１記載の路面状態推定装置は、車両が走行することによりタイヤから発生するタイヤ発生音を検出するタイヤ発生音検出手段と、前記タイヤ発生音検出手段により検出されたタイヤ発生音から、路面状態及びタイヤ発生音の発生源に対応する少なくとも１つの特徴量を計算する計算手段と、複数の路面状態の各々毎に、路面状態及びタイヤ発生音の発生源に対応する特徴量を記憶する記憶手段と、前記計算手段により計算された特徴量と、前記記憶手段に複数の路面状態の各々毎に対応して記憶された特徴量と、に基づいて、車両が走行する路面の状態を推定する推定手段と、を備えている。
【０００８】
タイヤ発生音検出手段は、車両が走行することによりタイヤから発生するタイヤ発生音を検出する。
【０００９】
計算手段は、前記タイヤ発生音検出手段により検出されたタイヤ発生音から、路面状態及びタイヤ発生音の発生源に対応する少なくとも１つの特徴量を計算する。
【００１０】
記憶手段は、複数の路面状態の各々毎に対応して、路面状態及びタイヤ発生音の発生源に対応する特徴量を記憶する。
【００１１】
推定手段は、前記計算手段により計算された特徴量と、前記記憶手段に複数の路面状態の各々毎に対応して記憶された特徴量と、に基づいて、車両が走行する路面の状態を推定する。
【００１２】
上記タイヤ発生音の発生源としては、少なくともタイヤへの物体の相対的な衝突を含むようにしてもよい。以下これを例にとり、本発明の原理を説明する。なお、以下は例示であり、本発明は、以下に限定されるものではない。
【００１３】
タイヤが回転（車両が走行）し、タイヤに、路面上の小石、路面自体の突起部、路面上の水等が、相対的に衝突することにより、このような衝突を音源としたタイヤ発生音が発生する。湿った路面や氷盤路面上には水が存在するので、タイヤが回転（車両が走行）し、タイヤに路面上の水が相対的に衝突すると、タイヤ発生音の発生源としての、少なくともタイヤへの物体の相対的な衝突に特有な特徴量が得られる。
【００１４】
そして、計算手段により計算された特徴量が、タイヤ発生音の発生源としての、少なくともタイヤへの物体の相対的な衝突に特有な特徴量に一致又は近似していると判断できれば、湿った路面や氷盤路面上、即ち、少なくとも湿った路面と推定することができる。
【００１５】
このように、タイヤ発生音自体を扱うのではなく、タイヤ発生音から、路面状態及びタイヤ発生音の発生源に対応する少なくとも１つの特徴量を計算しているので、路面の状態を推定するのに必要な情報を少なくすることができ、計算された特徴量と、複数の路面状態の各々毎に対応して記憶された路面状態及びタイヤ発生音の発生源に対応する特徴量と、に基づいて、車両が走行する路面の状態を推定するので、路面の状態を精度よく推定することができる。
【００１６】
ところで、上記特徴量は、車輪速によって値が異なる。
【００１７】
そこで、請求項２記載の発明は、車輪速を検出する車輪速検出手段を更に備え、前記記憶手段は、複数の路面状態の各々及び複数の車輪速の各々毎に対応して特徴量を記憶し、前記推定手段は、前記計算手段により計算された特徴量、前記車輪速検出手段により検出された車輪速、及び、前記記憶手段に複数の路面状態の各々及び複数の車輪速の各々毎に対応して記憶された路面状態に基づいて、車両が走行する路面の状態を推定する。
【００１８】
このように、複数の路面状態の各々及び複数の車輪速の各々毎に対応して特徴量を記憶し、計算された特徴量、検出された車輪速、及び、複数の路面状態の各々及び複数の車輪速の各々毎に対応して記憶された路面状態に基づいて、車両が走行する路面の状態を推定するので、速度を加味して路面の状態を推定することだでき、より路面状態を精度よく推定することができる。
【００１９】
ここで、計算手段は、請求項３のように、タイヤ発生音検出手段により検出されたタイヤ発生音の波形と少なくとも１つの予め定められた長さの基本波形との相関値を求めることにより、路面状態及びタイヤ発生音の発生源に対応する少なくとも１つの特徴量を計算するようにしてもよい。なお、請求項４のように、計算手段はウェーブレット変換を用いるようにしてもよい。
【００２０】
また、請求項５のように、計算手段は、タイヤ発生音検出手段により検出されたタイヤ発生音の波形と少なくとも１つのウェーブレットの基本波形との相関値を求めるようにしてもよい。
【００２１】
また、基本波形として、請求項６記載のように、各々振幅が第１の範囲内の値である２つの第１の領域と、該２つの第１の領域に挟まれ、振幅が該第１の範囲より大きくかつ該第１の範囲より広い第２の範囲内の値である第２の領域と、を備えた波形を用いるようにしてもよい。
【００２２】
更に、基本波形は、ドベシイ、シムレット、コアレット、及びメイエ等のウェーブレット関数であるとしてもよい。
【００２３】
また、請求項７のように、計算手段は、タイヤ発生音検出手段により検出されたタイヤ発生音の波形と少なくとも１つの路面状態に対応する基本波形との相関値を求めるようにしてもよい。この場合、請求項８のように、計算手段は、予め定められた原波形を路面状態に対応するように変形することにより、基本波形を設定する。
【００２４】
請求項９記載の発明の路面摩擦状態推定装置は、路面状態を表す特徴量を検出する検出手段と、各々異なる路面状態を表す複数の特徴量に対応して路面とタイヤとの間の摩擦状態である路面摩擦状態を記憶する記憶手段と、前記検出手段により検出された路面状態を表す特徴量と、前記記憶手段に各々異なる路面状態を表す複数の特徴量に対応して記憶された路面摩擦状態と、に基づいて、前記検出手段により前記特徴量が検出される対象となる路面の路面摩擦状態を推定する推定手段と、前記検出手段により前記特徴量が検出される対象となる路面の予め求められた路面摩擦状態と、前記推定手段により推定された路面摩擦状態と、に基づいて、前記記憶手段に記憶された路面摩擦状態を修正する修正手段と、を備えている。
【００２５】
本発明にかかる検出手段は、路面状態を表す特徴量を検出する。記憶手段は、各々異なる路面状態を表す複数の特徴量に対応して路面とタイヤとの間の摩擦状態である路面摩擦状態（路面摩擦状態を表す物理量、例えば、路面摩擦係数（以下同様））を記憶する。推定手段は、検出手段により検出された路面状態を表す特徴量と、記憶手段に各々異なる路面状態を表す複数の特徴量に対応して記憶された路面摩擦状態と、に基づいて、検出手段により特徴量が検出される対象となる路面の路面摩擦状態を推定する。
【００２６】
そして、本発明の修正手段は、検出手段により特徴量が検出される対象となる路面の予め求められた路面摩擦状態と、推定手段により推定された路面摩擦状態と、に基づいて、記憶手段に記憶された路面摩擦状態を修正する。
【００２７】
このように、検出手段により特徴量が検出される対象となる路面の予め求められた路面摩擦状態と、推定手段により推定された路面摩擦状態と、に基づいて、記憶手段に記憶された路面摩擦状態を修正するので、記憶手段に記憶された路面摩擦状態を適正な値にすることができ、路面摩擦状態を精度よく推定することができる。
【００２８】
ここで、請求項１０のように、前記検出手段は、前方路面の路面状態を表す特徴量を検出するようにしてもよく、請求項１１のように、検出手段は、現在走行している路面の路面状態を表す特徴量を検出するようにしてもよい。また、請求項１２のように、検出手段は請求項１に記載の路面状態推定装置でもよい。なお、検出手段は上記何れの路面状態推定装置でもよい。
【００２９】
請求項１３記載の路面摩擦状態推定装置は、現在走行している路面の路面状態を表す特徴量を検出する第１の検出手段と、前方路面の路面状態を表す特徴量を検出する第２の検出手段と、前記第２の検出手段により検出された前方路面の路面状態を表す特徴量を、前記第２の検出手段により前記特徴量が検出される対象となった路面の前記第１の検出手段により検出された前記特徴量に同期させる同期手段と、各々異なる路面状態を表す複数の特徴量に対応して路面とタイヤとの間の摩擦状態である路面摩擦状態を記憶する記憶手段と、前記第１の検出手段により検出された現在走行している路面状態を表す特徴量、前記同期手段により同期された前記第２の検出手段により検出された前方路面の路面状態を表す特徴量、及び前記記憶手段に各々異なる路面状態を表す複数の特徴量に対応して記憶された路面摩擦状態と、に基づいて、前記第１の検出手段により前記特徴量が検出される対象となる路面の路面摩擦状態を推定する推定手段と、を備えている。
【００３０】
本発明の第１の検出手段は、現在走行している路面の路面状態を表す特徴量を検出する。なお、第１の検出手段は請求項１乃至請求項４の何れか１項に記載の路面状態推定装置でもよい。第２の検出手段は、前方路面の路面の路面状態を表す特徴量を検出する。このように、特徴量を検出する対象は、第１の検出手段では、現在走行している路面であるのに対し、第２の検出手段では、前方路面である。従って、各々の検出手段が検出する特徴量は、時間的にずれている。そこで、同期手段は、第２の検出手段により検出された前方路面の路面状態を表す特徴量を、第２の検出手段により特徴量が検出される対象となった路面の第１の検出手段により検出された特徴量に同期させる。例えば、第２の検出手段により検出された前方路面の路面状態を表す特徴量を、該前方路面を走行するまで遅延させるようにしてもよい。
【００３１】
記憶手段は、各々異なる路面状態を表す複数の特徴量に対応して路面とタイヤとの間の摩擦状態である路面摩擦状態を記憶する。
【００３２】
推定手段は、第１の検出手段により検出された現在走行している路面状態を表す特徴量、同期手段により同期された第２の検出手段により検出された前方路面の路面状態を表す特徴量、及び記憶手段に各々異なる路面状態を表す複数の特徴量に対応して記憶された路面摩擦状態と、に基づいて、第１の検出手段により特徴量が検出される対象となる路面の路面摩擦状態を推定する。
【００３３】
このように、現在走行している路面の路面摩擦状態を推定するのに、現在走行している路面状態を表す特徴量の他に、前方路面の路面状態を表す特徴量を用いているので、多くの情報を得て路面摩擦状態を推定することができ、現在走行している路面の路面摩擦状態を精度よく推定することができる。
【００３４】
請求項１４記載に発明は、前記第１の検出手段により前記特徴量が検出される対象となる路面の予め求められた路面摩擦状態と、前記推定手段により推定された路面摩擦状態と、に基づいて、前記記憶手段に記憶された路面摩擦状態を修正する修正手段を備えている。
【００３５】
このように、前記第１の検出手段により前記特徴量が検出される対象となる路面の予め求められた路面摩擦状態と、推定手段により推定された路面摩擦状態と、に基づいて、記憶手段に記憶された路面摩擦状態を修正するので、記憶手段に記憶された路面摩擦状態を適正な値にすることができ、路面摩擦係数を精度よく推定することができる。
【００３６】
請求項１５記載の路面状態物理量演算装置は、車両が走行予定の前方路面の輝度情報を検出する検出手段と、前記検出手段により検出された前方路面の輝度情報と、路面状態を表す路面状態物理量を演算するためのパラメータと、に基づいて、路面状態物理量を演算する演算手段と、を備えた路面状態物理量演算装置であって、前記演算手段は、前記検出手段により検出された前方路面の輝度情報に基づいて前記パラメータを修正し、該修正したパラメータと、前記検出手段により検出された前方路面の輝度情報と、に基づいて、路面状態を表す路面状態物理量を演算することを特徴とする。
【００３７】
即ち、本発明は、車両が走行予定の前方路面の輝度情報を検出する検出手段と、前記検出手段により検出された前方路面の輝度情報と、路面状態を表す路面状態物理量を演算するためのパラメータと、に基づいて、路面状態物理量を演算する演算手段と、を備えた路面状態物理量演算装置である。
【００３８】
このように、前方路面の輝度情報と、路面状態を表す路面状態物理量を演算するためのパラメータと、に基づいて、路面状態物理量を演算すると、環境条件が変わると、前方路面の輝度情報の値が変わり、精度良く路面状態物理量を得られない。
【００３９】
そこで、本発明では、演算手段は、検出手段により検出された前方路面の輝度情報に基づいてパラメータを、環境条件に応じた路面状態物理量が適正に演算されるように、修正し、該修正したパラメータと、前記検出手段により検出された前方路面の輝度情報と、に基づいて、路面状態を表す路面状態物理量を演算するようにしている。よって、環境条件が変わってもパラメータが修正され、精度良く路面状態物理量を得ることができる。
【００４０】
この場合、請求項１６のように、検出手段は、前方路面側の領域を複数の領域に分割しかつ該分割した複数の領域各々の輝度情報を検出し、上記パラメータは、検出された複数の領域各々の輝度情報を各路面状態毎に分類するための値であり、演算手段は、検出された複数の領域各々の輝度情報を前記パラメータにより各路面状態毎に分類し、路面状態物理量を演算するようにしてもよい。また、請求項１７のように、検出手段は、路面側の領域を複数の領域に分割しかつ該分割した複数の領域各々の輝度情報を検出し、上記パラメータは、検出された複数の領域各々の輝度情報が得られる確率を最大にする値であり、演算手段は、パラメータにより最大となった、上記検出された複数の領域各々の輝度情報が得られる確率を演算することにより、路面状態物理量を演算するようにしてもよい。
【００４１】
請求項１８記載の路面状態推定装置は、複数の路面状態の各々毎に、輝度情報に基づいて定まる各路面の路面状態を表す路面状態物理量を記憶する記憶手段と、前方路面の輝度情報を検出する検出手段と、前記検出手段により検出された前方路面の輝度情報と、前記記憶手段により複数の路面状態の各々毎に記憶された各路面の路面状態を表す路面状態物理量と、に基づいて、前方路面の路面状態を推定する推定手段と、を備えている。
【００４２】
即ち、記憶手段は、複数の路面状態の各々毎に、輝度情報に基づいて定まる各路面の路面状態を表す路面状態物理量を記憶する。なお、輝度情報に基づいて定まる路面状態物理量には、例えば、輝度分布等である。
【００４３】
検出手段は、前方路面の輝度情報を検出する。
【００４４】
そして、推定手段は、検出手段により検出された前方路面の輝度情報と、記憶手段により複数の路面状態の各々毎に記憶された各路面の路面状態を表す路面状態物理量と、に基づいて、前方路面の路面状態を推定する。即ち、推定手段は、例えば、検出手段により検出された前方路面の輝度情報から路面状態物理量を求め、該求めた路面状態物理量と、記憶手段により複数の路面状態の各々毎に記憶された各路面の路面状態物理量と、に基づいて、前方路面の路面状態を推定する。即ち、記憶手段により複数の路面状態の各々毎に記憶された各路面の路面状態物理量の内、検出手段により検出された前方路面の輝度情報から求められた路面状態物理量に最も近い路面状態物理量を求め、求めた路面状態物理量に対応する路面状態を前方路面の路面状態として推定する。
【００４５】
このように、複数の路面状態の各々毎に路面状態物理量を記憶し、検出された前方路面の輝度情報と、複数の路面状態の各々毎に記憶された各路面の路面状態物理量と、に基づいて、前方路面の路面状態を推定するので、前方路面の路面状態を精度よく推定することができる。
【００４６】
なお、推定手段は、検出手段により検出された前方路面の輝度情報と、路面状態物理量を演算するためのパラメータと、に基づいて、路面状態物理量を演算し、該演算した路面状態物理量と、前記記憶手段により複数の路面状態の各々毎に記憶された各路面の路面状態を表す路面状態物理量と、に基づいて、前方路面の路面状態を推定するようにしてもよい。
【００４７】
この場合、推定手段は、検出手段により検出された前方路面の輝度情報に基づいて上記パラメータを修正し、該修正したパラメータと、検出手段により検出された前方路面の輝度情報と、に基づいて、路面状態を表す路面状態物理量を演算するようにしてもよい。
【００４８】
また、検出手段は、前方路面側の領域を複数の領域に分割しかつ該分割した複数の領域各々の輝度情報を検出し、パラメータは、検出された複数の領域各々の輝度情報が得られる確率を最大にする値であり、推定手段は、パラメータにより最大となった、上記検出された複数の領域各々の輝度情報が得られる確率を演算することにより、路面状態物理量を演算するようにしてもよい。
【００４９】
請求項２０記載の路面状態推定装置は、請求項１８記載の発明において、請求項１記載の路面状態推定装置と、前記路面状態推定装置により推定された路面状態と、前記路面状態推定装置により路面状態が推定される際に前記検出手段により検出された輝度情報と、に基づいて、前記記憶手段に記憶されかつ該路面状態推定装置により推定された路面状態に対応する路面状態物理量を更新する更新手段と、を更に備えている。
【００５０】
即ち、本発明は、路面状態推定装置により推定された路面状態と、路面状態推定装置により路面状態が推定される際に検出手段により検出された輝度情報と、に基づいて、記憶手段に記憶されかつ該路面状態推定装置により推定された路面状態に対応する路面状態物理量を更新するので、前方路面の路面状態をより精度よく推定することができる。
【００５１】
請求項２１記載の路面状態報知装置は、車両が現在走行している現在走行路面の路面状態を検出する現在路面状態検出手段、前方路面の路面状態を検出する前方路面状態検出手段、及び、前記現在路面状態検出手段により検出された現在走行路面の路面状態と、前記前方路面状態検出手段により検出された前方路面の路面状態と、を比較し、何れか一方の路面状態が変化した場合に、該変化した路面状態を報知する報知手段と、を備えている。
【００５２】
現在路面状態検出手段は、車両が現在走行している現在走行路面の路面状態を検出し、前方路面状態検出手段は、車両が走行予定の前方路面の路面状態を検出する。
【００５３】
報知手段は、現在路面状態検出手段により検出された現在走行路面の路面状態と、前方路面状態検出手段により検出された前方路面の路面状態と、を比較し、何れか一方の路面状態が変化した場合に、該変化した路面状態を報知する。
【００５４】
即ち、例えば、前方路面の路面状態が現在路面の路面状態から変化した場合、例えば、現在路面の路面状態が乾燥路面で、前方路面の路面状態が乾燥路面から圧雪路面に変化したとき、報知手段は、該変化した路面状態である、圧雪路面を報知する。更に、現在路面の路面状態が、乾燥路面から圧雪路面に変化した場合、報知手段は、該変化した路面状態である、圧雪路面を報知する。
【００５５】
このように、現在路面状態検出手段により検出された現在走行路面の路面状態と、前方路面状態検出手段により検出された前方路面の路面状態と、を比較し、何れか一方の路面状態が変化した場合に、該変化した路面状態を報知するので、路面状態の変化に対する路面状態の報知の応答性を向上させることができる。
【００５６】
請求項２３記載の発明にかかる路面状態報知装置は、現在走行している現在走行路面の路面状態を検出する現在路面状態検出手段と、前方路面の路面状態を検出する前方路面状態検出手段と、前記前方路面状態検出手段により検出された前方路面の路面状態を、前記現在路面状態検出手段により現在路面の路面状態が検出される対象となった現在路面に同期させる同期手段と、前記現在路面状態検出手段により検出された現在路面の路面状態と、前記同期手段により、前記現在路面状態検出手段により現在路面の路面状態が検出される対象となった現在路面に同期された前記前方路面状態検出手段により検出された前方路面の路面状態と、に基づいて、現在路面の路面状態を報知する報知手段と、を備えている。
【００５７】
即ち、現在路面状態検出手段は、車両が現在走行している現在走行路面の路面状態を検出し、前方路面状態検出手段は、前方路面の路面状態を検出する。
【００５８】
同期手段は、前方路面状態検出手段により検出された前方路面の路面状態を、現在路面状態検出手段により現在路面の路面状態が検出される対象となった現在路面に同期させる。
【００５９】
報知手段は、現在路面状態検出手段により検出された現在路面の路面状態と、同期手段により、現在路面状態検出手段により現在路面の路面状態が検出される対象となった現在路面に同期された、前方路面状態検出手段により検出された前方路面の路面状態と、に基づいて、現在路面の路面状態を報知する。
【００６０】
このように、検出された現在路面の路面状態と検出かつ同期された前方路面の路面状態とに基づいて、現在路面の路面状態を報知するので、路面状態を精度よく報知することができる。
【００６１】
請求項２４記載の発明にかかる路面摩擦状態推定装置は、車両が走行する路面の路面状態を検出する路面状態検出手段と、前記路面と前記車両のタイヤとの間の摩擦状態である路面摩擦状態を推定する路面摩擦状態推定手段と、前記路面状態検出手段により検出された路面状態と、前記路面摩擦状態推定手段により検出された路面摩擦状態と、に基づいて、最終の路面摩擦状態を推定する最終路面摩擦状態推定手段と、を備えている。
【００６２】
即ち、路面状態検出手段は、車両が走行する路面の路面状態を検出し、路面摩擦状態推定手段は、該路面と該車両のタイヤとの間の摩擦状態である路面摩擦状態を推定する。
【００６３】
最終路面摩擦状態推定手段は、路面状態検出手段により検出された路面状態と、路面摩擦状態推定手段により検出された路面摩擦状態と、に基づいて、最終の路面摩擦状態を推定する。即ち、最終路面摩擦状態推定手段は、路面状態検出手段により検出された路面状態と、路面摩擦状態推定手段により検出された路面摩擦状態と、に基づいて、検出された路面状態が変わらない又は変化量が所定値以内の場合には、路面摩擦状態推定手段により検出された路面摩擦状態の平均値（移動平均等）を、最終の路面摩擦状態として求める。
【００６４】
このように、路面状態検出手段により検出された路面状態と、路面摩擦状態推定手段により検出された路面摩擦状態と、に基づいて、最終の路面摩擦状態を推定するので、路面摩擦状態を路面状態の変化に応じて処理することができ、路面摩擦状態をより精度よく推定することができる。
【００６５】
なお、路面摩擦状態を路面状態の変化に応じて処理することができるとしても、何らかの影響により、検出された路面摩擦状態が正確なものでない場合には、最終的に推定された路面摩擦状態が、正確でなくなる。
【００６６】
そこで、請求項２５記載の発明では、請求項２４において、複数の路面状態の各々毎に、各路面状態に対応する路面摩擦状態の範囲を予め定めておき、前記各路面状態に対応する路面摩擦状態の範囲、前記路面状態検出手段により検出された路面状態、及び前記路面摩擦状態推定手段により検出された路面摩擦状態に基づいて、前記路面摩擦状態検出手段により検出された路面摩擦状態が、前記路面状態検出手段により検出された路面状態に対応する路面摩擦状態の範囲内か否かを判断する判断手段を更に備え、前記最終路面摩擦状態推定手段は、前記判断手段により、前記路面摩擦状態検出手段により検出された路面摩擦状態が、前記路面状態推定手段により検出された路面状態に対応する路面摩擦状態の範囲内であると判断された場合に、前記路面状態検出手段により検出された路面状態と、前記路面摩擦状態推定手段により検出されかつ該路面状態に対応する路面摩擦状態の範囲内であると判断された路面摩擦状態と、に基づいて、最終の路面摩擦状態を推定するようにしている。
【００６７】
即ち、複数の路面状態の各々毎に、各路面状態に対応する路面摩擦状態の範囲が予め定められている。
【００６８】
判断手段は、各路面状態に対応する路面摩擦状態の範囲、路面状態検出手段により検出された路面状態、及び路面摩擦状態推定手段により検出された路面摩擦状態に基づいて、路面摩擦状態検出手段により検出された路面摩擦状態が、路面状態検出手段により検出された路面状態に対応する路面摩擦状態の範囲内か否かを判断する。即ち、例えば、検出された路面摩擦状態が、路面状態検出手段により検出された路面状態に対応する路面摩擦状態の範囲外と判断された場合には、検出された路面摩擦状態は正確でないと判断することができる。よって、この正確でない路面摩擦状態を用いて最終の路面摩擦状態を推定すると、路面摩擦状態の推定値は正確でなくなる。
【００６９】
そこで、最終路面摩擦状態推定手段は、判断手段により、路面摩擦状態検出手段により検出された路面摩擦状態が、路面状態推定手段により検出された路面状態に対応する路面摩擦状態の範囲内であると判断された場合に、路面状態検出手段により検出された路面状態と、路面摩擦状態推定手段により検出されかつ該路面状態に対応する路面摩擦状態の範囲内であると判断された路面摩擦状態と、に基づいて、最終の路面摩擦状態を推定するようにしている。よって、路面摩擦状態を精度よく推定することができる。
【００７０】
請求項２６記載の路面状態推定装置は、車両が走行することによりタイヤから発生するタイヤ発生音を検出するタイヤ発生音検出手段と、車輪速を検出する車輪速検出手段と、複数の路面状態の各々毎に、タイヤ発生音の発生源に対応する特徴量と車輪速との関係を記憶する記憶手段と、前記タイヤ発生音検出手段により検出されたタイヤ発生音、前記車輪速検出手段により検出された車輪速、及び前記記憶手段に前記複数の路面状態の各々毎に記憶された係数に基づいて、路面状態を推定する推定手段と、を備えている。
【００７１】
タイヤ発生音検出手段は、車両が走行することによりタイヤから発生するタイヤ発生音を検出し、車輪速検出手段は、車輪速を検出する。
【００７２】
記憶手段は、タイヤ発生音の発生源に対応する特徴量と車輪速との関係を記憶する。なお、当該関係を表す式を定めるための係数を記憶することにより、当該関係を記憶するようにしてもよい。なお、該式としては例えば２次式でもよい。
【００７３】
推定手段は、タイヤ発生音検出手段により検出されたタイヤ発生音、車輪速検出手段により検出された車輪速、及び記憶手段に複数の路面状態の各々毎に記憶された関係に基づいて、路面状態を推定する。
【００７４】
このように、路面状態を推定するために、複数の路面状態の各々毎にタイヤ発生音と車輪速との関係を記憶しているので、各路面状態毎にタイヤ発生音及び車輪速を記憶する場合よりも記憶量を少なくすることができる。
【００７５】
そして、請求項２７記載の発明のように、前記タイヤ発生音検出手段により検出された車両が走行することによりタイヤから発生するタイヤ発生音と、前記車輪速検出手段により検出された車輪速と、に基づいて、前記記憶手段に記憶された関係を更新する更新手段を更に備え、前記推定手段は、前記タイヤ発生音検出手段により検出されたタイヤ発生音、前記車輪速検出手段により検出された車輪速、及び前記記憶手段に前記複数の路面状態の各々毎に記憶されかつ前記更新手段により更新された関係に基づいて、路面状態を推定するようにしてもよい。
【００７６】
なお、以上説明した種々の発明においては、タイヤ発生音検出手段は、タイヤ発生音に、該タイヤ発生音の発生源に基づくタイヤ発生音以外の外乱の影響が小さくなる位置に配置する。例えば、タイヤ近傍のバンパー内部等である。
【００７７】
また、路面の輝度情報を検出する手段は、車両の一部及び前方路面を含む領域の輝度情報を検出し、該手段により検出された輝度情報に基づいて、車両の一部の輝度情報と前方路面の輝度情報との関係を演算かつ出力する演算出力手段を更に設けるようにしても良い。
【発明の効果】
【００７８】
以上説明したように請求項１乃至請求項８に記載の発明は、タイヤ発生音自体を扱うのではなく、タイヤ発生音から、路面状態及びタイヤ発生音の発生源に対応する少なくとも１つの特徴量を計算しているので、路面の状態を推定するのに必要な情報を少なくすることができ、計算された特徴量と、複数の路面状態の各々毎に対応して記憶された路面状態及びタイヤ発生音の発生源に対応する特徴量と、に基づいて、車両が走行する路面の状態を推定するので、路面の状態を精度よく推定することができる、という効果を有する。
【００７９】
請求項９乃至請求項１２に記載の発明は、検出手段により特徴量が検出される対象となる路面の予め求められた路面摩擦状態と、推定手段により推定された路面摩擦状態と、に基づいて、記憶手段に記憶された路面摩擦状態を修正するので、記憶手段に記憶された路面摩擦状態を適正な値にすることができ、路面摩擦状態を精度よく推定することができる、という効果を有する。
【００８０】
請求項１３及び請求項１４記載の発明は、現在走行している路面の路面摩擦状態を推定するのに、現在走行している路面状態を表す特徴量の他に、前方路面の路面状態を表す特徴量を用いているので、多くの情報を得て路面摩擦状態を推定することができ、現在走行している路面の路面摩擦状態を精度よく推定することができる、という効果を有する。
【００８１】
請求項１５乃至請求項１７記載の発明は、検出された前方路面の輝度情報に基づいてパラメータを修正し、該修正したパラメータと、検出された前方路面の輝度情報と、に基づいて、路面状態を表す路面状態物理量を演算するので、環境条件が変わってもパラメータが修正され、精度良く路面状態物理量を得ることができる、という効果を有する。
【００８２】
請求項１８乃至請求項２０記載の発明は、複数の路面状態の各々毎に路面状態物理量を記憶し、検出された前方路面の輝度情報と、複数の路面状態の各々毎に記憶された各路面の路面状態物理量と、に基づいて、前方路面の路面状態を推定するので、前方路面の路面状態を精度よく推定することができる、という効果を有する。
【００８３】
請求項２１及び請求項２２記載の発明は、現在路面状態検出手段により検出された現在走行路面の路面状態と、前方路面状態検出手段により検出された前方路面の路面状態と、を比較し、何れか一方の路面状態が変化した場合に、該変化した路面状態を報知するので、路面状態の変化に対する路面状態の報知の応答性を向上させることができる、という効果を有する。
【００８４】
請求項２３記載の発明は、検出された現在路面の路面状態と検出かつ同期された前方路面の路面状態とに基づいて、現在路面の路面状態を報知するので、路面状態を精度よく報知することができる、という効果を有する。
【００８５】
請求項２４及び請求項２５記載の発明は、路面状態検出手段により検出された路面状態と、路面摩擦状態推定手段により検出された路面摩擦状態と、に基づいて、最終の路面摩擦状態を推定するので、路面摩擦状態を路面状態の変化に応じて処理することができ、路面摩擦状態をより精度よく推定することができる、という効果を有する。
【００８６】
請求項２６及び請求項２７記載の発明は、路面状態を推定するために、複数の路面状態の各々毎にタイヤ発生音と車輪速との関係を記憶しているので、各路面状態毎にタイヤ発生音及び車輪速を記憶する場合よりも記憶量を少なくすることができる、という効果を有する。
【発明を実施するための最良の形態】
【００８７】
以下、図面を参照して、本発明の第１の実施の形態を詳細に説明する。
【００８８】
図１に示すように、本実施の形態に係る路面状態推定装置は、車両が走行することによりタイヤから発生するタイヤ発生音を検出するタイヤ発生音検出手段としてのタイヤ発生音検出センサ１０と、車輪速を検出する車輪速検出手段としての車輪速検出センサ１２と、路面状態推定器１００と、を備えている。
【００８９】
ここで、タイヤ発生音検出センサ１０は、タイヤ近傍で、エンジン音及び風切り音などの外部の影響を受け難く、さらに、水、石等の障害物が極力当たらない位置に配置されている。具体的には、図示しないタイヤ近傍に車両振動音の影響を受けないように配置している。例えば、エンジンが車両の前側に位置する場合には、タイヤ発生音検出センサ１０を車両後部のタイヤ近傍のボデイに配置する。なお、タイヤ発生音検出センサ１０としてはマイクを用いることができる。また、車輪速検出センサ１２は、車輪速として車輪速の回転角速度を検出するセンサである。
【００９０】
また、路面状態推定器１００は、タイヤ発生音検出センサ１０により検出されたタイヤ発生音から、路面状態及びタイヤ発生音の発生源に対応する少なくとも１つの特徴量を計算する計算手段としてのＩＣ等で構成された前処理器１４と、前処理器１４により計算された特徴量と、車輪速検出センサ１２により検出された車輪速と、を成分とする特徴ベクトルを生成するＩＣ等で構成された特徴ベクトル生成器１６と、推定装置１８と、を備えている。
【００９１】
推定装置１８は、複数の路面状態の各々及び複数の車輪速の各々毎に対応して特徴量をマップとして記憶する判定マップメモリ２２と、特徴ベクトル生成器１６により生成された特徴ベクトルと、判定マップメモリ２２に複数の路面状態の各々及び複数の車輪速の各々毎に対応してマップとして記憶された特徴量と、に基づいて、車両が走行する路面の状態を推定するＩＣ等で構成された判定処理器２０と、を備えている。
【００９２】
ここで、判定マップメモリ２２に記憶されている特徴量、即ち、複数の路面状態の各々及び複数の車輪速の各々毎に対応してマップとして記憶された特徴量を説明するが、最初に、タイヤと路面間で生じる音の発生要因について説明する。
【００９３】
タイヤと路面間で生じる音の発生源について、例えば、横浜ゴム株式会社編「自動車用タイヤの研究」山海堂、酒井秀男著「タイヤ工学」グランプリ出版に記載されており、その主な発生源は、
ポンピングノイズ
ゴムの衝突ノイズ
タイヤボディの振動ノイズ
スキールノイズ
である。
【００９４】
ポンピングノイズは、トレッド構内の空気が圧縮、放出のポンピング作用による気柱共鳴で発生し、その基本波周波数は１．５ｋＨｚ〜２．５ｋＨｚである。ポンピングノイズは、路面が平滑だと大きくなる。
【００９５】
ゴムの衝突ノイズは、トレッド表面と路面表面の小石や水等とが相対的に衝突したときに発生し、不等ピッチパターンでホワイトノイズ化され、その周波数は１ｋＨｚ以上である。
【００９６】
タイヤボディの振動ノイズは、転動（タイヤの回転）時にトレッドとサイドウォールが加振され、その振動が空気伝播してノイズとなるもので、加振力としてトレッドゴムブロックによる第１の振動ノイズと、路面凹凸による第２の振動ノイズと、がある。これらの周波数は数百〜数ｋＨｚである。
【００９７】
スキールノイズは急ブレーキ、急コーナリング時におけるトレッドゴムのスティックスリップにより発生し、周波数は１ｋＨｚ〜２ｋＨｚとなる。
【００９８】
以上説明した発生源から分かるように、タイヤ発生音は、路面状態及び発生源、例えば、タイヤの種類、トレッドパターン、走行速度に影響を受ける。
【００９９】
図２は、スタッドレスタイヤで異なる路面を所定速度で走行した時の時間−周波数解析結果を示すものである。音の発生源からこの周波数分布を考察すると、２０００Ｈｚ以下の部分Ａがタイヤボディの振動ノイズ、基本波が２０００Ｈｚの部分Ｂがポンピングノイズにそれぞれ相当すると考えられる。また１０００Ｈｚ以上の部分Ｃは、ホワイトノイズ成分であり路面表面上の水等とトレッド表面とが衝突したり、該水等が跳ね上げられたりする際に発生する成分と考えられる。
【０１００】
従って、湿潤アスファルト路面は、図２（Ｂ）に示すように、１ｋＨｚ以上のホワイトノイズ成分（部分Ｃ参照）が大きく、逆に乾燥アスファルト路面ではこれが小さい。
【０１０１】
氷盤路面は平坦なため、図２（Ｃ）に示すように、基本波が２０００Ｈｚのポンピングノイズ成分（部分Ｂ参照）が大きく、逆に路面凹凸によるタイヤボディ振動ノイズ（部分Ａ参照）が小さくなっている。
【０１０２】
圧雪路面は、図２（Ｄ）に示すように、ポンピングノイズ、ホワイトノイズ成分、トレッドゴムブロック加振による成分とも小さくなっている。
【０１０３】
このように路面状態により音の発生源に依存する音の成分の大きさが異なる為、これら成分を抽出する事により路面状態を判別する事が可能となる。
【０１０４】
図３は、乾燥アスファルト路及び湿潤アスファルト路、図４は、永盤路及び圧雪路の、それぞれ予め把握された各路面のタイヤ発生音を、ドベシィのウェーブレット（Ｎ＝３）を用いて、ウェーブレット変換処理したものである。また図５は、各ウェーブレット成分を時間−周波数解析したものである。
【０１０５】
図３〜図５から理解されるように、特にｄ１の成分は２０００Ｈｚのポンピングノイズと２０００Ｈｚ以上のゴム衝突ノイズを表しており、ｄ４の成分は２ｋＨｚ以下のトレッドゴムブロック加振によるタイヤボディの第１の振動ノイズ、ｄ５〜ｄ９の成分は１ｋＨｚ以下の路面凹凸によるタイヤボディの第２の振動ノイズを表している事が分かる。したがって特にこの成分が各路面で異なっている事が分かる。
【０１０６】
以上は、車輪速がある所定値の場合であるが、車輪速が変化すると、タイヤ発生音をウェーブレット変換処理した値も変化する。
【０１０７】
本実施形態では、複数の車輪速のそれぞれ毎及び複数の路面状態のそれぞれ毎に、路面状態及びタイヤ発生音の発生源に対応する特徴量としての・タイヤ発生音をウェーブレット変換した値を、判別マップメモリ２２に記憶している。即ち、予め複数の路面の各々毎の代表的なタイヤ発生音と車速とから、前処理器１４と特徴ベクトル生成器１６で得られた特徴ベクトル（特徴量（タイヤ発生音をウェーブレット変換処理した値）と車速との組み合わせ）を、各路面状態毎にクラス分けし、プロトタイプとして判定マップメモリ２２に記憶している。
【０１０８】
次に、本実施の形態の作用を説明する。
【０１０９】
車両が走行することによりタイヤから発生するタイヤ発生音を、タイヤ発生音検出センサ１０が検出する。また、このときの車輪速を車輪速検出センサ１２が検出する。
【０１１０】
前処理器１４は、タイヤ発生音検出センサ１０により検出されたタイヤ発生音から、路面状態及びタイヤ発生音の発生源に対応する特徴量を計算する。即ち、前処理器１４は、タイヤ発生音検出センサ１０により検出されたタイヤ発生音を信号処理解析することにより、タイヤ発生音を路面状態に依存する各発生源の成分に分離する。ここで信号処理解析として、タイヤ発生音の波形と少なくとも１つの予め定められた長さの基本波形との相関値、即ち、例えば、タイヤ発生音の信号とウェーブレット基本波との内積をとって、時間領域の相似相関を検出するウェーブレット解析を用いることができる。なお、これに代えて、各種フィルタ処理を用いるようにしてもよい。
【０１１１】
図６に、ウェーブレット変換を用いた前処理器１４のアルゴリズムを示す。
【０１１２】
これは、Mallatによって１９８９年に発見された高速ウェーブレット変換アルゴリズムである（A．Mallat：“A Theory for Multiresolution Signal Decomposition：the Wavelet Representation，”IEEE Trans．PAMI，Vol. 11，No．7，ppp．674−693，1989．参照）。
【０１１３】
ステップ３２で、タイヤ発生音をサンプリングして得られる数列をＳk(j)とする。ステップ３４で、レベルjのスケーリング係数から、１レベル精度の低いレベルj＋１のウェーブレット展開係数及びスケーリング係数を導出する。
【０１１４】
【数１】
2004168286.tif 000003
【０１１５】
ここで、ＳK(j)はスケーリング係数、ＷK(j)はウェーブレット展開係数、ｐはスケーリング関数の数列、ｑはウェーブレットの数列(基本波形)を表し、ｋはシフト量、ｊはスケールを表す。
【０１１６】
特徴ベクトル生成器１６は、前処理器１４で分離抽出した各成分と、車輪速検出センサ１２により検出された車輪速と、の組み合わせ、即ち、これらを成分とする特徴ベクトルを生成する。そして、特徴ベクトル生成器１６は、このようにして求められた特徴ベクトルを判定処理器２０に供給する。
【０１１７】
判定処理器２０は、特徴ベクトル生成器１６で生成された特徴ベクトルと、判定マップメモリ２２に記憶された・複数の路面状態の各々毎に対応して記憶された路面状態及びタイヤ発生音の発生源に対応する特徴ベクトルと、に基づいて、車両が走行する路面の状態を推定する。即ち、判定処理器２０では、入力された特徴ベクトルが、特徴空間（特徴量と速度を成分する空間）上で判定マップメモリ２２に記憶されているプロトタイプと比較され、最も距離の近いプロトタイプの属するクラスを路面状態の判別結果として出力する。距離としてはユークリッド距離を用いる。
【０１１８】
なお、識別処理は、最近傍決定則（ＮＮ法）、ベイズ決定則、線形識別関数、ニューラルネットワークなどを適用することができる。識別手法としてニューラルネットワークを用いる場合は、ニューラルネットワークの重み係数が識別辞書のデータとなる。
【０１１９】
以上説明したように本実施の形態では、広範囲の周波数成分を有するタイヤ発生音自体を扱うのではなく、タイヤ発生音から、路面状態及びタイヤ発生音の発生源に対応する少なくとも１つの特徴量を計算しているので、路面の状態を推定するのに必要な情報を少なくすることができ、計算された特徴量と、複数の路面状態の各々毎に対応して記憶された路面状態及びタイヤ発生音の発生源に対応する特徴量と、に基づいて、車両が走行する路面の状態を推定するので、路面の状態を精度よく推定することができる。
【０１２０】
即ち、図３及び図４に示すように、例えば、湿潤路面を走行したときのタイヤ発生音をウエーブレットで分解した波形から理解できるように、タイヤ発生音は、低周波数から高周波数（a12,d12,d11,…d1）に渡る各信号で構成されている。各信号の波形は、周期的ではなく、突発的は変化をする特徴を有する。各信号の周期も異なる。フーリエ変換は、波形が、周期的ではなく、突発的は変化をする特徴を有する信号の解析としては不十分である。よって、フーリエ変換を用いる場合では、路面の状態を精度よく推定することができない。
【０１２１】
上記のようにウエーブレット変換に代えて、タイヤ発生音を、複数のタイヤ発生音の発生源各々に応じて定まる帯域フィルタを通過させることにより、分解し、解析するようにしてもよい。即ち、図１３には、４つの異なる路面を走行したときのタイヤ発生音を、帯域フィルタを通過させることにより、分解した波形が示されている。この分解された波形を具体的に解析する場合には、各波形の実効値を求める。図１４（Ａ）及び図１４（Ｂ）は、４つの異なる路面を走行したときのタイヤ発生音を、帯域フィルタを通過させることにより分解した波形から求めた特徴量（波形の実効値）を示した図である。
【０１２２】
ただし、図１４（Ａ）及び図１４（Ｂ）に示すように、４つの異なる路面を走行したときのタイヤ発生音の特徴量は比較的狭い領域内に分布している。従って、核路面の特徴量を明確に区別することは難しい。よって、タイヤ発生音を、帯域フィルタを通過させることにより、分解した波形から求めた特徴量からでは、タイヤ発生音をウエーブレットで分解する場合に比較すると、路面状態を正確に判別することができない。
【０１２３】
一方、上記実施の形態では、複数の路面状態、例えば、４つの異なる路面を走行したときのタイヤ発生音をウエーブレットで分解した波形から特徴量を求めており、各特徴量は広い範囲に散らばって分布しているので、特徴量から、路面状態を正確に判別することができる。以下、これを更に詳細に説明する。
【０１２４】
タイヤ発生音をウエーブレットで分解する際に使用する基本波形は、各々振幅が第１の範囲内の値である２つの第１の領域と、該２つの第１の領域に挟まれ、振幅が該第１の範囲より大きくかつ該第１の範囲より広い第２の範囲内の値である第２の領域と、を有する。即ち、例えば、図１５に示すドベシイのウェーブレットの基本波形を説明する。図１５に示すように、ドベシイのウェーブレットの基本波形は、各々振幅が第１の範囲Ｌ１内の値である２つの第１の領域Ｒ１（Ｒ１１、Ｒ１２）と、該２つの第１の領域Ｒ１（Ｒ１１、Ｒ１２）に挟まれ、振幅が該第１の範囲Ｌ１より大きくかつ該第１の範囲Ｌ１より広い第２の範囲Ｌ２内の値である第２の領域Ｒ２と、を有する。なお、以下、同様に、図１６〜図１８に示すように、シムレット、コアレット、及びメイエのウェーブレットの基本波形も同様である。
【０１２５】
複数の路面、例えば、乾燥、湿潤、氷盤、及び圧雪の４つの異なる路面を走行したときのタイヤ発生音を、上記ドベシイのウェーブレット（図１５参照）で分解した波形は、図３及び図４に示すようになり、このように分解した波形から求めた特徴量（相関値）は、図１９（Ａ）及び図１９（Ｂ）に示すようになる。なお、図中では、乾燥、湿潤、氷盤、及び圧雪をそれぞれ、Dry, Wet, Ice, snowと記載している。
【０１２６】
上記図１４（Ａ）及び図１４（Ｂ）と図１９（Ａ）及び図１９（Ｂ）と比較すると、図１９（Ａ）及び図１９（Ｂ）に示す、上記ドベシイのウェーブレットで分解した波形から求めた特徴量（相関値）の各路面毎の値は、上記図１４（Ａ）及び図１４（Ｂ）に示す特徴量の分布領域より、広い。このように、各路面の特徴量は広い範囲に散らばって分布しているので、特徴量から、路面状態を正確に判別することができる。
【０１２７】
なお、シムレット、コアレット、及びメイエのウェーブレットで分解した波形から求めた特徴量（相関値）の各路面毎の値も、広い範囲に散らばって分布しているので、特徴量から、路面状態を正確に判別することができる。
【０１２８】
次に、上記前処理器１４の具体的内容を図２０を参照して説明する。
【０１２９】
図２０に示すように、上記前処理器１４は、タイヤ発生源波形波長データを記憶するメモリ１４Ａ、マザー・ウェーブレット（ウエーブレット関数）を発生するウェーブレット関数発生器１４Ｂ、メモリ１４Ａに記憶されたタイヤ発生源波形波長データと、ウェーブレット関数発生器１４Ｂにより発生されたマザー・ウェーブレットと、から、タイヤ発生源の波形を解析するためのウェーブレット基底を演算するウェーブレット基底演算器１４Ｃ、及び、タイヤ発生音検出センサ１０により検出されたタイヤ発生音とウェーブレット基底演算器１４Ｃにより演算されたウェーブレット基底との内積を演算して相似相関値を検出する相似相関演算器１４Ｃを備えている。
【０１３０】
ウェーブレット関数発生器１４Ｂは、マザー・ウェーブレットを発生する。このマザー・ウェーブレットは、以下の式を満たす２乗可積分関数に属する。
【０１３１】
【数２】
2004168286.tif 000004
【０１３２】
さらに、マザー・ウェーブレットは、
【０１３３】
【数３】
2004168286.tif 000005
【０１３４】
となる関数の総称である。数３における左の式（第1の条件）から、ψ（ｔ）は、振動的で、台（関数が０とならない区間の閉包）が有界となるものである。そして、数３における右の式（第2の条件）から、ψ（ｔ）は、その変動が有限で１に標準化されている。マザー・ウェーブレットは、これら２つの条件をみたす限り任意のものでよい。
【０１３５】
ウェーブレット基底演算器１４Ｃは、メモリ１４Ａに記憶されたタイヤ発生源波形波長データと、ウェーブレット関数発生器１４Ｂにより発生されたマザー・ウェーブレットと、から、タイヤ発生源の波形を解析するためのウェーブレット基底を演算する。即ち、ウェーブレット基底演算器１４Ｃは、タイヤ発生音の発生源となる各発生源波形を検出するための、ウェーブレット基底を演算する。
【０１３６】
ウェーブレット基底は次式で表される。
【０１３７】
【数４】
2004168286.tif 000006
【０１３８】
aは、波長の拡大・縮小を制御するスケール・パラメータ、ｂは時間を調整するシフト・パラメータである。
【０１３９】
ウェーブレット基底演算器１４Ｃは、メモリ１４Ａに記憶されたタイヤ発生源波形波長データからウエーブレット基底の波長がタイヤ発生源波長と等しくなるようスケール・パラメータaを設定する。即ち、複数の路面状態の各々の路面を予め走行することにより、ウエーブレット基底の波長がタイヤ発生源波長と等しくなるようスケール・パラメータaを設定する。
【０１４０】
相似相関演算器１４Ｃは、タイヤ発生音検出センサ１０により検出されたタイヤ発生音とウェーブレット基底演算器１４Ｃにより演算された各ウェーブレット基底との内積を演算して相似相関値を検出する。即ち、相似相関演算器１４Ｃは、次式のように、設定したウェーブレット基底とタイヤ発生音ｘ（ｔ）との内積を演算する。
【０１４１】
【数５】
2004168286.tif 000007
【０１４２】
そして、相似相関演算器１４Ｃは、シフト・パラメータｂ変化させて、各時点での相似相関値を演算して、タイヤ発生音から、その発生源となる音の成分を検出する。
【０１４３】
このように、一般に間欠音との相似性の強いウェーブレット関数を用いて、さらにタイヤ発生源となる各波形の波長データから、そのタイヤ発生源波形と相似なウェーブレット基底を求め、そのウェーブレット基底とタイヤ発生音の内積を演算するので、タイヤ発生音からその発生源となる成分のみを精度良く検出できる。
【０１４４】
ここで、
【０１４５】
【数６】
2004168286.tif 000008
【０１４６】
とおき、変形すると、数５の式は、
【０１４７】
【数７】
2004168286.tif 000009
【０１４８】
となり、インパルス応答がｇ（ｔ）であるフィルタにｘ（ｔ）を入力したときの出力と等価である。
【０１４９】
即ち、相関を利用した信号処理法に整合フィルタがある。これはある波形と一致した波形が入力されたときに大きな出力が得られるように設計されたフィルタである。たとえばｚ（ｔ）という波形に一致したときに出力が大きくなるようなフィルタを考える。フィルタのインパルス応答をｇ（ｔ）＝ｚ（−ｔ）のようにする。入力をｕ（ｔ）、出力をｙ（ｔ）とすると、
【０１５０】
【数８】
2004168286.tif 000010
【０１５１】
となる。すなわち入力にｚ（ｔ）と同じ波形が入ると、出力はｚ（ｔ）の自己相関関数をｔだけずらしたものになる。一方で、ｚ（ｔ）と相関のない入力に対しては出力は小さくなる。
【０１５２】
ウェーブレット変換はスケーリング付の整合フィルタと考えることが出来る。すなわち、ウェーブレット変換はインパルス応答がｇ（ｔ）であるフィルタにｘ（ｔ）を入力したときの出力と等価である。
【０１５３】
このフィルタを整合フィルタと考えると、
【０１５４】
【数９】
2004168286.tif 000011
【０１５５】
即ち、スケーリングしたアナライジングウェーブレットと一致した波形が入力されたときに大きな出力が得られる整合フィルタということになる。
（変形例）
次に、第１の実施の形態の変形例を説明する。以下に示す２つの変形例は何れも、上記前処理器１４の部分が異なる。よって、主として異なる部分のみを説明する。
【０１５６】
図２１に示すように、第１の変形例にかかる前処理器１４Ｎ１は、タイヤ発生音の各発生源波形に対し、相関のある各相似波形を記憶する相似波形メモリ１４Ｇ、及び、タイヤ発生音と、該記憶した各相似波形の相関を演算して、タイヤ発生音からその各発生源波形を検出する相関演算器１４Ｈを備えている。
【０１５７】
相似波形メモリ１４Ｇは、タイヤ発生音の各発生源波形に対し、相関のある各相似波形を記憶する。具体的には、図２３に示すように、タイヤ発生音と相関をとるために、複数の路面、例えば、上記のように、乾燥路面、湿潤路面、圧雪路面、及び氷盤路面などの各路面毎に対応して基本波形を記憶する。即ち、各路面を予め走行することにより検出した各路面の基本波形を記憶する。なお、図２３（Ａ）は、乾燥路面の基本波形、図２３（Ｂ）は、湿潤路面の基本波形を示している。
【０１５８】
相関演算器１４Ｈは、タイヤ発生音と、上記記憶した各基本波形の相関を演算して、タイヤ発生音からその各発生源波形を検出する。
【０１５９】
具体的には、タイヤ発生音ｘ（ｉ）と基本波形ｓ（ｉ）との相関を次式より演算する。
【０１６０】
【数１０】
2004168286.tif 000012
【０１６１】
即ち、図２４に示すように、一定期間内のタイヤ発生音ｘ（ｉ）（図２４（Ａ）参照）に対し、例えば、乾燥路面の基本波形ｓ（ｉ）（図２３（Ａ）、図２４（Ｂ）参照）を用いて、相関をとると、図２４（Ｃ）に示すように相関値Ｒｘｓ（ｉ）が得られる。図２４（Ａ）乃至図２４（Ｃ）に示すように、一定期間内のタイヤ発生音ｘ（ｉ）の中で、基本波形ｓ（ｉ）に近似している部分について大きな相関値Ｒｘｓ（ｉ）が得られる。
【０１６２】
次に、第２の実施の形態の変形例を説明する。図２２に示すように、第２の変形例にかかる前処理器１４Ｎ２は、タイヤ発生音の各発生源波形に対し、相関のある各相似波形を設定する相似波形設定器１４Ｉ、及び、タイヤ発生音と、該設定された各相似波形の相関を演算して、タイヤ発生音からその各発生源波形を検出する相関演算器１４Ｊを備えている。
【０１６３】
相似波形設定器１４Ｉは、各路面状態のタイヤ発生源波形と相似な局在波形（原波形）を用いて、その波長を設定するすることにより、基本波形を設定する。図２５は、次式による局在波系の一例である。
【０１６４】
【数１１】
2004168286.tif 000013
【０１６５】
相似波形設定器１４Ｊは、上式のパラメータａを、各タイヤ発生源波形を最も良く検出できるよう設定する。即ち、各路面を予め走行することにより、各路面の波形に一致又は近似するように、上式のパラメータａを設定する。例えば、図２５に示すように、局在波形（a＝１）Ｗ１を、走行した路面の波形に一致又は近似するように、パラメータａを１より小さくして、基本波形Ｗ２を設定する。
【０１６６】
なお、相関演算器１４Ｊは、タイヤ発生音と、上記設定された各基本波形の相関を演算して、タイヤ発生音からその各発生源波形を検出する。なお、詳細は上記第１の変形例と同様であるので、その説明を省略する。
【０１６７】
なお、第２の変形例では、各路面の基本波形を、上記数式のパラメータを変更することにより、設定しているので、各路面の基本波形事態を記憶する第１の変形例に比べると、記憶量を少なくすることができる。
【０１６８】
次に、本発明の第２の実施の形態を説明する。本実施の形態は、路面と車輪そとの間の最大摩擦係数である最大路面摩擦係数μsを推定する路面摩擦係数推定装置（μs推定装置）５０に、路面状態推定器を備えたものである。即ち、μs推定装置５０は、図７に示すように、第１の実施の形態の路面状態推定器１００と同構成の路面状態推定器５２と、各路面状態毎に対応して最大路面摩擦係数μsを記憶する路面μ推定器５４と、から構成されている。第１の実施の形態と比較すると、路面状態から最大路面摩擦係数μsを推定する機能を付加した点相違する。
即ち、路面μ推定器５４は、路面状態推定器５２により推定された路面状態と、路面μ推定器５４が各路面状態毎に対応して記憶する最大路面摩擦係数μsと、から、現在走行している路面（タイヤ発生音が検出された路面）の最大路面摩擦係数μsを推定するものである。
【０１６９】
最大路面摩擦係数μsは路面状態と走行速度、水膜圧、タイヤ種類に依存することが、横浜ゴム株式会社編「自動車用タイヤの研究」山海堂に示されている。
【０１７０】
路面μ推定器５４では、路面状態判別部で判別した路面状態から最大路面摩擦係数μsを関連付けて出力する。
【０１７１】
次に、本発明の第３の実施の形態を説明する。図８に示すように、本実施の形態は、前述した第２の実施の形態（図７参照）と略同様の構成であるが、車輪速が路面μ推定器６４にも入力されるように構成されている点で相違する。
【０１７２】
即ち、本実施の形態に係るμs推定装置６０は、第１の実施の形態の路面状態推定器１００と同構成の路面状態推定器６２と、各路面状態毎及び車輪速毎に対応して最大路面摩擦係数μsを記憶する路面μ推定器６４と、から構成されている。
【０１７３】
路面μ推定器６４は、車輪速及び路面状態推定器６２で判別された路面状態と、
路面μ推定器６４が各路面状態毎及び車輪速毎に対応して記憶している最大路面摩擦係数μsと、から現在走行している路面（タイヤ発生音が検出された路面）の最大路面摩擦係数μsを推定し、これを出力する。
【０１７４】
第１の実施の形態と比較すると、路面状態と車速から最大路面摩擦係数を推定する機能を付加した点相違する。
【０１７５】
次に、本発明の第４の実施の形態を説明する。
【０１７６】
図９に示すように、本実施の形態に係るμs推定装置７０は、第１の実施の形態に係る路面状態推定器１００と同構成の路面状態判別器７２、路面μ推定器７８、μ勾配推定器７４、及びタイヤ種別判定器７６から構成されている。
【０１７７】
本実施の形態に係るμs推定装置７０では、路面状態推定器７２で路面状態を推定する。
【０１７８】
μ勾配推定器７４で車輪と路面との間の摩擦係数のスリップ速度に対する勾配であるμ勾配を推定する。即ち、μ勾配推定器７４は、例えば、特願２００２−０４０７２２号公報に記載の方法により前輪及び後輪の車輪速差と前後加速度とを用いてμ勾配を推定する。
【０１７９】
タイヤ種別判定器７６では、μ勾配に基づいてタイヤ種別を判定する。タイヤ種別判定器７６は、例えば、特願２００２−０４０７２２号公報に記載の方法により、μ勾配推定器７４で推定されたμ勾配を用いて、サマータイヤ（磨耗したスタッドレスタイヤを含む。）であるか、スタッドレスタイヤであるかのタイヤの種別を判定する。
【０１８０】
そして、路面μ推定器７８は、路面状態、車輪速、及びタイヤ種から、路面μ推定器７８で現在走行している路面（タイヤ発生音が検出された路面）の最大路面摩擦係数μsを推定する。第１の実施の形態と比較すると、路面状態とタイヤ種別と車輪速から最大路面摩擦係数μsを推定する機能を付加した点相違する。
【０１８１】
次に、本発明の第５の実施の形態を説明する。
【０１８２】
本実施の形態に係るμs推定装置は、図１０に示すように、定常走行時の最大路面摩擦係数μsを推定するμs推定器８０と、加減速走行時の最大路面摩擦係数μｇ１を推定するμｇ１推定器８２と、旋回走行時の最大路面摩擦係数μｇ２を推定するμｇ２推定器８４と、最大路面摩擦係数μｇ１、μｇ２のいずれか一方を選択するμｇ選択器８６と、重み係数Ｋを決定する重み係数決定器８８と、総合的な最大路面摩擦係数μｍａｘを推定するμｍａｘ推定器９０と、を備えている。
【０１８３】
μs推定器８０には、前述したタイヤ発生音及び車輪速が入力される。μｇ１推定器８２には、図示しないセンサから前後加速度（加減速度）が入力される。μｇ２推定器８４には、図示しないセンサから横加速度が入力される。
【０１８４】
以下、各器を説明する。
【０１８５】
μs推定器８０としては、図７、図８、図９のμs推定器５０、６０、７０の何れも採用することができる。
【０１８６】
μｇ１推定器８２は、例えば、特願２００２−０４０７２２号公報に記載の方法により、車輪速及び前後加速度に基づいてμ勾配を推定し、μ勾配に基づいて加減速走行時のグリップ度Ｇｒ１及び最大路面摩擦係数μｇ１を推定する。推定された最大路面摩擦係数μｇ１はμｇ選択器８６に供給される。
【０１８７】
μｇ２推定器８４は、例えば、特願２００２−０４０７２２号公報に記載の方法により、車輪速からμ勾配を推定し、μ勾配に基づいて旋回走行時のグリップ度Ｇｒ２及び細大路面摩擦係数μｇ２を推定する。推定された最大路面摩擦係数μｇ２はμｇ選択器８６に供給される。
【０１８８】
μｇ選択器８６は、μｇ１推定器８２で推定された最大路面摩擦係数μｇ１と、μｇ２推定器８４で推定された最大路面摩擦係数μｇ２と、のうちの大きな値の方を選択し、選択したものを最大路面摩擦係数μｇとしてμｍａｘ推定器９０に供給する。
【０１８９】
なお、μｇ選択器８６は、グリップ度Ｇｒ１，Ｇｒ２のうち小さい方に対応する最大路面摩擦係数を選択してもよい。また、μｇ選択器８６は、操舵角センサで検出された操舵角に基づいて旋回走行時であるかを判定し、旋回走行時であるときは最大路面摩擦係数μｇ２を選択し、旋回走行時でないときは最大路面摩擦係数μｇ１を選択してもよい。
【０１９０】
重み係数決定器８８は、μｇ選択器８６で選択された最大路面摩擦係数μｇに対するグリップ度Ｇｒに基づいて、重み係数Ｋを決定する。具体的には、等該グリップ度Ｇｒが１に近づくに従って重み係数Ｋを１に近づけ、グリップ度Ｇｒがゼロに近づくに従って重み係数Ｋをゼロに近づける。即ち、車両が定常走行している場合は、グリップ度Ｇｒは大きな値（１に近い値）になるので、μｇよりもμsの信頼度を高くするために、重み係数Ｋを１に近い値にする。
【０１９１】
また、車両が加減速走行又は旋回走行している場合は、グリップ度Ｇｒは小さな値になるので、μsよりもμｇの信頼度を高くするために、重み係数Ｋをゼロに近い値にする。
【０１９２】
なお、重み係数決定器８８は、トラクション・コントロール・システム（TRC：Traction Control System）,アンチロック・ブレーキ・システム（ABS：Antilock Brake System）等の車両安定性制御システム（VSC：Vehicle Stability Control System）の制御が行われている場合は、重み係数Ｋをゼロに近づけてもよい。
【０１９３】
μｍａｘ推定器９０は、例えば、特願２００２−０４０７２２号公報に記載の方法により、定常走行時と、加減速走行時又は旋回走行時とを総合的に考慮した最大路面摩擦係数μｍａｘを推定する。
【０１９４】
以上のように、定常走行時の最大路面摩擦係数μsと、加減速走行時時又は旋回走行時の最大路面摩擦係数μｇとに基づいて、総合的な最大路面摩擦係数μｍａｘを推定することができる。
【０１９５】
このとき、加減速走行時又は旋回走行時のグリップ度Ｇｒに基づいて重み付け係数Ｋを求めるので、定常走行時、加減速走行時又は旋回走行時の状態をそれぞれ考慮することができ、制度良く最大路面摩擦係数μｍａｘを推定することができる。
【０１９６】
次に、本発明の第６の実施の形態を説明する。
【０１９７】
図１１に示すように、本実施の形態に係る路面摩擦係数推定装置は、前方路面特徴量検出器１０２、現在路面特徴量検出器１０４、及び路面摩擦係数真値記憶装置１０７を備えている。
【０１９８】
ここで、前方路面特徴量検出器１０２は、前方路面特徴量を検出する装置として、次の装置を適用することができる。即ち、例えば、カメラで前方路面を撮影して得られた画像の水平、垂直偏光画像のパワースペクトル分布を検出する装置（特開平成８年−３２７５３０号報参照）、カメラで前方路面を撮影して得られた画像の色相、彩度、明度それぞれの平均、及び分散値を検出する装置（特開２００２−１４０７８９号報参照）等である。また、インフラ側の前方路面センサや前方車両を前方路面特徴量検出器１０２として採用することができ、これらにより計測された路面μ値等を前方路面の特徴量として検出するようにしてもよい。
【０１９９】
また、現在路面特徴量検出器１０４は、現在路面特徴量を検出する装置として、次の装置を適用することができる。即ち、例えば、ロードノイズの周波数スペクトルを検出する装置（特開平６−１３８０１８号公報、特開平７−１５６７８２号公報等参照）や、超音波の路面反射特性を検出する装置（特開平９−５４４９号公報参照）、光の路面反射特性を検出する装置（特開平７−４９２２３号公報参照）や、前述した第１の実施の形態乃至第５の実施の形態等である。
【０２００】
更に、路面摩擦係数真値記載装置１０７は、車両状態量から測定される路面摩擦係数で、車両挙動の限界付近の路面摩擦係数、あるいは比較的長い期間の統計をとったより現在路面及び前方路面の確からしい路面摩擦係数を記憶している。なお、具体的内容は、特願２００２−０４０７２２号の出願明細書に記載されている。
【０２０１】
前方路面特徴量検出器１０２には、特徴ベクトル生成器１０４が接続されている。特徴ベクトル生成器１０４には、μ判定処理器１１０及びμ判定マップメモリ１１２を備えた前方路面μ推定器１０８が接続されている。μ判定処理器１１０の出力側には、遅れ時間処理器１１４が接続されている。
【０２０２】
現在路面特徴量検出器１０４には、特徴ベクトル生成器１２０が接続されている。特徴ベクトル生成器１２０には、μ判定処理器１２４及びμ判定マップメモリ１２６を備えた前方路面μ推定器１２２が接続されている。μ判定処理器１２４の出力側には、遅れ時間処理器１２８が接続されている。なお、特徴ベクトル生成器１２０には、前方路面特徴量検出器１０２の前方路面特徴量を遅れ時間処理して特徴ベクトル生成器１２０に入力する遅れ時間処理器１０６が接続されている。
【０２０３】
遅れ時間処理器１１４と路面摩擦係数真値記憶装置１０７とには減算器１１６が接続されていると共に、遅れ時間処理器１２８と路面摩擦係数真値記憶装置１０７とには減算器１３０が接続されている。減算器１１６にはμ判定マップメモリ１１２が記憶している値を修正するマップ修正処理器１１８が接続されていると共に、減算器１３０にはμ判定マップメモリ１２６が記憶している値を修正するマップ修正処理器１３２が接続されている。
【０２０４】
次に、本実施の形態の作用を説明する。本実施の形態は、前方路面μと現在路面μをそれぞれその路面状態とμ判定マップの路面μとを対応させることにより、路面μを推定する装置で、過去の車両挙動等から測定された路面μの真値とその真値の時点の前方路面μの値と現在路面μの値との比較を行い、真値の路面μに近づくようにμ判定マップを学習する。以下、詳細に説明する。
【０２０５】
前方路面特徴量検出器１０２は、上記前方路面の特徴量を検出する。特徴ベクトル生成器１０４は、前方路面特徴量を１つのベクトルとし、前方路面μ判定器１０８に入力する。
【０２０６】
前方路面μ判定器１０８は、前方路面特徴量を入力し、対応する前方路面の摩擦状態（摩擦状態を表す物理量）、例えば、摩擦係数を出力するための識別処理を行う。識別処理は、近傍決定則（ＮＮ法）、ベイズ決定則、線形識別関数、ニューラルネットワークなどを適用することができる。識別処理として近傍決定則（ＮＮ法）を用いる場合、あらかじめ各路面状態の前方路面データをそれに対応する路面摩擦係数でクラス分けし、プロトタイプとしてμ判定マップメモリ１１２に記憶しておく、μ判定処理器１１０では、判定時に入力される前方路面データがμ判定マップメモリ１１２に記憶されているプロトタイプと比較され、最も距離の近いプロトタイプの属する路面摩擦係数のクラスを前方路面摩擦係数の判定結果として出力する。距離としてはユークリッド距離が使われるときが多い。識別手法としてニューラルネットを用いる場合は、ニューラルネットの重み係数がμ判定マップメモリ１１２のデータとなる。
【０２０７】
現在路面特徴量検出器１０４は、上記現在路面の特徴量を検出する。遅れ時間処理処理器１０６は、前方路面特徴量を現在路面特徴量と同期を取るための処理を行う。車速の情報および前方路面データと現在路面データとの相関演算して、常に同期をとるため遅れ時間を制御する。前方路面特徴量はメモリに記憶しておき、遅れ時間後の特徴量をメモリから特徴ベクトル生成器１２０に出力する。
【０２０８】
特徴ベクトル生成器１２０は、遅れ時間処理された前方路面特徴量と、現在路面特徴量を１つのベクトルとし、現在路面μ判定器１２２へ入力する。
【０２０９】
現在路面μ推定器１２２は、遅れ時間処理した前方路面特徴量と現在路面特徴量から得られた特徴ベクトルを入力して、対応する現在路面の摩擦状態（摩擦状態を表す物理量）、例えば、摩擦係数を出力するための識別処理を行う。処理方法は前方路面μ判定器１０８と同様な処理を行う。
【０２１０】
ところで、路面摩擦係数真値記載装置１０７は、前述したように現在路面及び前方路面の確からしい路面摩擦係数を記憶している。
【０２１１】
遅れ時間処理器１１４と遅れ時間処理器１２５は、路面摩擦係数真値が測定された時点の前方路面摩擦係数、現在路面摩擦係数を出力する。前方路面摩擦係数と、現在路面摩擦係数をメモリに記憶しておき、路面摩擦係数真値が測定された時点の値をメモリから出力する。
【０２１２】
マップ修正処理器１１８は、路面摩擦係数真値と遅れ時間処理された前方路面摩擦係数との偏差に応じて、μ判定マップ１に記憶されているプロトタイプの摩擦係数の値を増減させ、路面摩擦係数真値に近づくようμ判定マップメモリ１１２内のμ判定マップを学習する。マップ修正処理器１３２も同様の処理を行い、μ判定マップメモリ１２６内のμ判定マップを学習する。
【０２１３】
以下、マップ修処理器１１８、１３２のマップ修正処理の具体的計算方法を２例示す。
具体例１
マップ修処理器１１８では、μ判定マップメモリ１１２に予め次の代表的な前方路面特徴量の特徴ベクトルを記憶しておく。
ｘ1＝（ｘ11，ｘ12，…ｘ1d）
ｘ2＝（ｘ21，ｘ22，…ｘ2d）
・
・
・
ｘn＝（ｘn1，ｘn2，…ｘnd）
ここで各要素ｘn1，ｘn2，…ｘnd等はそれぞれ前方路面特徴量を示す。
【０２１４】
さらに、上記各々の特徴ベクトルに下記の様に１対１で対応する前方路面摩擦係数μを記憶しておく。
ｘ1  →  μ1
ｘ2  →  μ2
・
・
・
ｘn  →  μn
今現在の前方路面特徴ベクトルをｘaとすると、μ判定処理器１１０では最初に、ｘaとμ判定マップ１に記憶されている特徴ベクトルｘ1，ｘ2，…ｘnとのｄ次元の特徴空間での距離を演算し、最も距離の短いマップの特徴ベクトルを選択する。
【０２１５】
ここで例えばマップの特徴ベクトルｘと対応する前方路面摩擦係数μを以下の様に設定したとする。
【０２１６】
乾燥アスファルト路面に相等  ｘ1  →  μ1  ＝  ０．８
湿潤アスファルト路面に相等  ｘ2  →  μ2  ＝  ０．６
圧雪路面                    ｘ3  →  μ3  ＝  ０．２
氷盤路面                    ｘ4  →  μ4  ＝  ０．１
上記演算の結果、ｘaに最も距離の近い特徴ベクトルがｘ1とすると、ｘ1のポインタが指示するアドレスに記憶されているμ1＝０．８を今現在の前方路面摩擦係数として出力する。
【０２１７】
こうして推定された摩擦係数μsは一旦記憶され、過去の路面摩擦係数真値μriと比較するために、遅れ時間処理し、その時点の推定摩擦係数と次式により比較し、その誤差を求める。
【０２１８】
一旦記憶される推定摩擦係数μsを時系列順に
μs1，μs2…μsi…μsnとする。
過去の路面摩擦係数真値をμri、その時点の推定摩擦係数をμsiとすると、
ｅi＝μri−μsi
μri＝０．９、μsi＝０．８だったとすると
ｅi＝０．１
従ってマップ修正処理器１１８では、次式により、乾燥アスファルト路面に相当するμ値を更新する。
μ1new＝μ1＋ｅi
＝０．８＋０．１
でｘ1に対するμ1を０．９に修正する。
【０２１９】
なお、マップ修正処理器１３２でも、現在路面摩擦推定に対し上記と同様の処理を行う。
【０２２０】
具体例２
ニューラルネットワークを用いて誤差逆伝播法によりマップを修正する計算例を示す。
【０２２１】
μ判定器処理器１１０、１２４を図１２に示すニューラルネットワークで構成する。
【０２２２】
ここで、入力層Ｋへ入力されるｘ0，…，ｘ0'，…，ｘ1dは特徴ベクトルの各要素を示し、出力層Ｌのμ1，…，μl，…，μ1dはその特徴ベクトルに対応する路面摩擦係数の値を示すベクトルの各要素を示す。
【０２２３】
μ判定器処理器１１０、１２４では、予め代表的な各種路面の特徴ベクトル
ｘ1＝（ｘ11，ｘ12，…ｘ1d）
ｘ2＝（ｘ21，ｘ22，…ｘ2d）
・
・
・
ｘn＝（ｘn1，ｘn2，…ｘnd）
をニューラルネットワークの学習データとし、それに１対１に対応する路面摩擦係数μ1，μ2，…，μnを教師データとして、学習データの特徴ベクトルを入力し、それに対応する教師データが出力されるよう、ニューラルネットワークを学習させておく。そして、学習で得られたニューラルネットワークの重み係数をμ判定マップメモリ１１２、１２６に記憶させておく。ここで出力層で出力される路面摩擦係数の値を示すベクトルについて説明する。
【０２２４】
出力層の路面摩擦係数ベクトルＭ＝（μ1，…，μl，…，μc）は例えば次のように摩擦係数値μsと対応させる。
Ｍ＝（１，…０，…０）    →  μs＝０．１
Ｍ＝（０，１…０，…０）  →  μs＝０．２
・
・
・
Ｍ＝（０，０…０，…１）  →  μs＝１．０
従って、特徴ベクトルｘ1＝（ｘ11，ｘ12，…ｘ1d）に対応する教師データμ1＝０．１の場合、ニューラルネットにｘ1を入力すると出力層ｙμs＝０．１に相当するＭ＝（１，…０，…０）になるようにニューラルネットを学習させる。
【０２２５】
このように、予め記憶させたμ判定マップによる摩擦係数の推定値の精度を上げるためには現在走行中の路面状況をより反映したμ判定マップに修正する必要がある。
【０２２６】
そのため走行中に得られた特徴ベクトルを記憶させておき、これを学習データとし、同時点での路面摩擦係数真値を教師データとして再学習しμ判定マップを修正する。
【０２２７】
再学習のタイミングは、推定摩擦係数と路面摩擦係数真値の誤差がある許容範囲を超えたときに行う。
【０２２８】
さて、この場合におけるμ判定マップの修正は、ニューラルネットワークの重み係数を誤差逆伝播法により行うことができる（誤差逆伝播法は「パターン認識」オーム社Ｐ４２〜Ｐ４８に詳細に記載されている）。
【０２２９】
以上説明したように、本実施の形態では、前方路面μと現在路面μのμ判定に用いるμ判定マップを車両の挙動等から得られた過去の真のμ値を用いて学習するため、種々のどんな路面状態においても正確な前方路面μと現在路面μが得られる。
【０２３０】
また、本実施の形態では、現在路面の路面摩擦係数を推定するのに、現在走行している路面状態を表す特徴量の他に、前方路面の路面状態を表す特徴量を用いているので、多くの情報を得て路面摩擦係数を推定することができ、現在走行している路面の路面摩擦係数を精度よく推定することができる。
【０２３１】
更に、本実施の形態では、現在路面の予め求められた路面摩擦係数と、推定された路面摩擦係数と、に基づいて、記憶された路面摩擦係数を修正するので、記憶された路面摩擦係数を適正な値にすることができ、路面摩擦係数を精度よく推定することができる。
【０２３２】
更に種々の実施の形態を説明する。なお、各実施の形態では、同一の構成部分を有するので、同一の構成部分には同一の符号を付して、その構成及び作用の詳細な説明を省略する。
【０２３３】
第７の実施の形態を説明する。
【０２３４】
図２６に示すように、本実施の形態にかかる路面状態物理量演算装置としての路面状態推定装置は、車両が走行予定の前方路面の輝度情報を検出する検出手段としての、前方の路面画像を検出する路面画像検出センサ１０Ａと、前記検出手段により検出された前方路面の路面画像と、路面状態を表す路面状態物理量を演算するためのパラメータと、に基づいて、路面状態物理量を演算し、演算した路面状態物理量から前方路面状態を推定する路面状態推定器１００と、を備えている。なお、路面状態推定器１００は、演算手段を構成する。
【０２３５】
より詳細には、路面画像検出センサ１０Ａは、前方路面側の領域を複数の領域に分割しかつ該分割した複数の領域各々の輝度情報を検出する例えば、ＣＣＤカメラ等により構成することができる。
【０２３６】
路面状態推定器１００は、後述するように、パラメータとしての閾値ｋを修正するとともに、路面画像検出センサ１０Ａにより検出された前方の路面画像の画像データの各画素の輝度データを、該修正された閾値ｋにより、輝度の低い分布のクラスＣ１と、輝度の高い分布のクラスＣ２と、の２つのクラスに分ける輝度分布演算器１１、及び、輝度分布演算器１１により前方の路面画像の各画素の輝度データが分けられ、分けられた輝度データの個数（画素数；路面状態物理量に対応する）を演算し、輝度データが輝度の高いクラスＣ１に、より多く分けられた場合には、圧雪か氷盤の第１の路面状態と判定し、輝度データが輝度の低いクラスＣ２に、より多く分けられた場合には、乾燥か湿潤路面の第２の路面状態と判定する判定処理器２０Ａを備えている。
【０２３７】
次に、本実施の形態の作用を説明する。
【０２３８】
最初に、輝度分布演算器１１Ａの演算の一手法について説明する（大津、粟田、関田著「パターン認識」朝倉書店参照）。
【０２３９】
路面画像検出センサ１０Ａで検出された最大輝度がＬレベルで表され、０からＬまでの間のレベルｉのデータ数をｎｉ、全データ数をＮとすると正規化ヒストグラムは、
【０２４０】
【数１２】
2004168286.tif 000014
【０２４１】
で表され、これは、輝度の確率分布とみなすことができる。輝度の全平均レベル、全分散は、それぞれ
【０２４２】
【数１３】
2004168286.tif 000015
【０２４３】
で与えられる。
【０２４４】
０からＬまでの間のレベルｋを閾値とし、ｋ以下の区間Ｓ１＝［１、・・・、ｋ］のデータをクラスＣ１、ｋより高い輝度の区間Ｓ２＝［ｋ＋１、・・・、Ｌ］のデータをクラスＣ２に分類するものとする。また、正規化ヒストグラムに対して、次の二つの累計量
【０２４５】
【数１４】
2004168286.tif 000016
【０２４６】
を定義する。
【０２４７】
レベルｋを閾値としたときの各クラスの生起確率は、この累計量を用いて、
【０２４８】
【数１５】
2004168286.tif 000017
【０２４９】
で与えられる。また、各クラスの平均レベルは
【０２５０】
【数１６】
2004168286.tif 000018
【０２５１】
で与えられる。
【０２５２】
ヒストグラムを閾値ｋにより二つのクラスに分けたとき、クラスの分離度がよければその閾値はよい閾値であることが期待できる。閾値ｋのよさを評価するために、判別基準
【０２５３】
【数１７】
2004168286.tif 000019
【０２５４】
を考える。ここで、
【０２５５】
【数１８】
2004168286.tif 000020
【０２５６】
は、それぞれ、クラス内分散およびクラス間分散である。クラス内分散、クラス間分散、および全分散の間には、
【０２５７】
【数１９】
2004168286.tif 000021
【０２５８】
の基本的関係がｋによらず常に成立する。
【０２５９】
σ2Wとσ2Bとはともにｋの関数であるが、前者が２次の統計量まで必要とするのに対して、後者は１次の統計量のみに基づいている。つまり、ηがｋに関して最も簡単な評価基準である。したがって、最適閾値として２クラスＣ１、Ｃ２の最もよい分離を与えるｋを探索する場合には、閾値ｋのよさの評価としてηを採用するのが最も効率的である。これは、σ2bを最大とするｋを最適閾値とすることと同値となる。よって最適なｋ*は、
【０２６０】
【数２０】
2004168286.tif 000022
【０２６１】
から、
【０２６２】
【数２１】
2004168286.tif 000023
【０２６３】
のように逐次ｋを進めて探索することができる。ここで、ｋは実質的にはｗ１ｗ２＞０、つまり、０＜ｗ（ｋ）＜１の範囲を動かせばよい。
【０２６４】
以上説明したように、輝度分布演算器１１Ａは、路面画像検出センサ１０Ａにより得られた前方路面の輝度情報から、輝度分布が適切に環境条件（日照条件等）に反映するように、閾値ｋを修正している。即ち、快晴の昼間と夕方とでは、同じ路面状態の路面を走行していても、輝度分布は異なる。そこで、本実施の形態では、このように、例えば、快晴の昼間と夕方とでも、同じ路面状態であれば同じ路面状態と判定されるように、閾値ｋを修正している。
【０２６５】
そして、輝度分布演算器１１は、路面画像検出センサ１０Ａにより検出された前方の路面画像の画像データの各画素の輝度データを、上記のように走行中オンラインで常に修正される閾値ｋにより、輝度の低い分布のクラスＣ１と、輝度の高い分布のクラスＣ２と、の２つのクラスに分ける。
【０２６６】
判定処理器２０Ａは、輝度分布演算器１１により前方の路面画像の各画素の輝度データが分けられた輝度データの個数（画素数；路面状態物理量に対応する）を演算し、輝度データが輝度の高いクラスＣ１に、より多く分けられた場合には、圧雪か氷盤の第１の路面状態と判定し、輝度データが輝度の低いクラスＣ２に、より多く分けられた場合には、乾燥か湿潤路面の第２の路面状態と判定する。
【０２６７】
以上説明したように第７の実施の形態では、路面状態を判定するための輝度分布の閾値は推定中一定ではなく、走行中オンラインで上記の演算をし、各路面状態の輝度分布を常に学習して、日照条件等の走行環境の変化に適応しているため、走行環境にロバストな推定が可能となる。
【０２６８】
次に、第８の実施の形態を説明する。
【０２６９】
図２７に示すように、第８の実施の形態の路面状態推定装置は、上記第７の実施の形態にかかる路面画像検出センサと同様の路面画像検出センサ１０Ａ、前方路面状態を推定する路面状態推定器１１０Ａ、各々上記第１の実施の形態のタイヤ発生音検出センサ及び路面状態推定器と同様の、タイヤ発生音を検出するタイヤ発生音検出センサ２０Ａ、及び、タイヤ発生音から走行中の路面状態を推定する路面状態推定器２００を備えている。
【０２７０】
上記路面状態１１０Ａは、路面状態推定器２００で得られた走行中の路面状態の路面状態物理量としての輝度分布を演算する輝度分布演算器１３Ａ、及び、輝度分布演算器１３Ａで得られた輝度分布により、路面画像検出センサ１０で得られた前方画像の輝度分布がどの路面状態の輝度分布かを判定する判定処理器２０Ａを備えている。なお、輝度分布演算器１３Ａは、複数の路面状態の各々毎に、輝度分布を記憶する記憶手段としてのメモリを備えている。
【０２７１】
次に、本実施の形態の作用を説明する。なお、タイヤ発生音検出センサ２０Ａ及び路面状態推定器２００の作用を前述した第１の実施の形態と同様であるので、その説明を省略する。
【０２７２】
路面状態推定器２００は走行中常に車両が走行している路面の路面状態を推定する。
【０２７３】
路面画像検出センサ１０Ａは前方路面側の画像の輝度情報を検出する。
【０２７４】
輝度分布演算器１３Ａは、前方路面が現在路面になるまで、即ち、上記前方路面を車両が走行するまで、路面画像検出センサ１０Ａにより検出された前方路面側の画像の輝度情報を遅れ時間処理することにより、路面状態推定器２００で推定された現在走行中の路面状態と同期をとる。
【０２７５】
輝度分布演算器１３は、路面状態の輝度分布を求める。以下、輝度分布演算器１３の演算を最尤推定法で求める一手法について説明する。
【０２７６】
まず、現在走行しているある路面の路面状態の輝度のデータの集合をＸ＝｛ｘ１、ｘ２、・・・、ｘｎ｝とし、推定すべき輝度の確率密度関数をｐ（ｘ；θ）で表す。ここでθはパラメータの組を表すベクトルで、パラメータベクトルと呼ばれる。
【０２７７】
ここで、１路面状態の輝度のデータの集合Ｘをもたらしたθとしてはいろいろなものが想定できるが、その中でどのθがいちばん尤もらしいかを考える。データの集合Ｘに含まれるデータは、確率ｐ（ｘ；θ）に従って独立に生起したものと考えられるから、このようなデータ集合が得られる確率ｐ（ｘ；θ）は次式であらわされる。
【０２７８】
【数２２】
2004168286.tif 000024
【０２７９】
したがって、いちばん尤もらしいθは、式（２・１）を最大にするθである。このようなθを
【０２８０】
【数２３】
2004168286.tif 000025
【０２８１】
と置き、これを推定値として用いることにすれば、
【０２８２】
【数２４】
2004168286.tif 000026
【０２８３】
であり、これを求めるには
【０２８４】
【数２５】
2004168286.tif 000027
【０２８５】
あるいは対数をとって
【０２８６】
【数２６】
2004168286.tif 000028
【０２８７】
を解けばよい。
【０２８８】
輝度の確率密度関数をｐ（ｘ；θ）が正規分布であると仮定すると、
【０２８９】
【数２７】
2004168286.tif 000029
【０２９０】
ここでｍは平均ベクトル、Σは共分散行列である。
【０２９１】
式（２・５）と式（２・４）より、ｍ、Σに対して、それぞれ次のような推定値
【０２９２】
【数２８】
2004168286.tif 000030
【０２９３】
が得られる。
【０２９４】
同様に、車両が他の路面を走行した場合には、当該他の路面状態の輝度分布も上記のように演算する。よって、複数の路面状態、例えば、乾燥、湿潤、圧雪、氷盤、の各路面状態に対応して輝度分布を走行中常に記憶（更新）する。
【０２９５】
判定処理器２０Ａは、路面画像検出センサ１０Ａで今回検出された前方路面の輝度分布がどの路面状態のものかを判定する。即ち、上判定処理器１２Ａは、上記のように、複数の路面状態、例えば、乾燥、湿潤、圧雪、氷盤、の各路面状態に対応して輝度分布を記憶している。そして、判定処理器２０Ａは、路面画像検出センサ１０Ａで今回検出された前方路面の輝度分布が、上記のように複数の路面状態の各々毎に記憶された輝度分布の何れにもっとも近いのかを判定し、判定された輝度分布に対応する路面状態を前方路面の路面状態の推定結果として出力する。
【０２９６】
以上説明したように第８の実施の形態では、第７の実施の形態同様、路面状態を判定するための輝度分布の閾値は推定中一定ではなく、走行中オンラインで上記の演算をし、各路面状態の輝度分布を常に学習して、日照条件等の走行環境の変化に適応しているため、走行環境にロバストな推定が可能となる。さらに、第８の実施の形態では、常にタイヤ発生音から推定された現在走行中の路面状態の輝度の分布をオンラインで学習して求め、実際の走行路面状態と推定の演算で用いる路面状態を照合しているため、より正確な推定結果を得ることが出来る。
【０２９７】
なお、輝度分布演算器１３は、第７の実施の形態の輝度分布演算器１１を用いて、輝度分布演算器１１でクラス分けした路面状態を、路面状態推定器２００で推定された路面状態推定結果で修正するようにしても良い。
【０２９８】
次に、第９の実施の形態を説明する。
【０２９９】
図２８に示すように、第９の実施の形態にかかる路面状態推定装置は、路面画像検出センサ１０Ａ、前方路面状態を推定する路面状態推定器１００Ａ（なお、上記路面状態推定器１１０Ａでもよい）、タイヤ発生音検出センサ２０Ａ、現在路面状態を推定する路面状態推定器２００、及び、現在走行路面の路面状態と前方路面の路面状態とを比較し、何れか一方の路面状態が変化した場合に、該変化した路面状態を報知する報知手段としての判定処理器２０を備えている。
【０３００】
なお、路面画像検出センサ１０Ａ、及び路面状態推定器１００Ａ（なお、上記路面状態推定器１１０Ａでもよい）は、車両が走行予定の前方路面の路面状態を検出する前方路面状態検出手段を構成し、タイヤ発生音検出センサ２０Ａ及び路面状態推定器２００は、車両が現在走行している現在走行路面の路面状態を検出する現在路面状態検出手段を構成する。
【０３０１】
次に、本実施の形態の作用を説明する。
【０３０２】
路面状態推定器１００Ａ（あるいは路面状態推定器１１０Ａ）は、第７の実施の形態（あるいは第８の実施の形態）で説明したように、前方路面の輝度情報から、前方の路面が圧雪あるいは氷盤の路面かその他の路面かの２状態を推定する。
【０３０３】
一方、路面状態推定器２００では、現在走行中の路面を乾燥、湿潤、圧雪、氷盤、の４状態で推定する。
【０３０４】
判定処理器３０では、路面状態推定器１００Ａ（あるいは路面状態推定器１１０Ａ）の推定結果と、路面状態推定器２００の推定結果と、から、前方路面の状態が現在路面の状態から変化したとき、例えば、路面状態推定器２００で推定された現在路面状態が乾燥路面で、路面状態推定器１００あるいは路面状態推定器１１０で推定された前方路面状態が乾燥か湿潤の一状態から圧雪か氷盤の一路面状態に変化したとき、判定処理器３０は、路面状態推定器１００あるいは路面状態推定器１１０で推定された結果をすみやかに出力して最終的な路面状態推定結果とする。
【０３０５】
さらに、その後、路面状態推定器１００あるいは路面状態推定器１１０で推定された結果に変化がなく、また路面状態推定器２００で推定された現在路面状態が乾燥路面から例えば圧雪路面に変化した場合は、最終的な推定結果を圧雪か氷盤の一路面状態から、そのうちの圧雪路面と、より詳細な路面状態を出力する。
【０３０６】
以上説明したように第９の実施の形態では、路面状態推定器１００Ａ（あるいは路面状態推定器１１０Ａ）で、前方の路面状態変化をとらえ、さらに、路面状態推定器１００あるいは路面状態推定器１１０では推定できない、より詳細な路面状態を、路面状態推定器２００で推定するため、従来の路面状態推定装置より、推定の応答性と精度を向上できる。
【０３０７】
次に、第１０の実施の形態を説明する。
【０３０８】
図２９に示すように、第１０の実施の形態にかかる路面状態推定装置は、路面画像検出センサ１０Ａ、路面状態推定器１００Ａ（あるいは路面状態推定器１１０Ａ）、遅れ時間処理器３０、タイヤ発生音検出センサ２０、及び路面状態推定器２１０を備えている。
【０３０９】
なお、路面画像検出センサ１０Ａ、及び路面状態推定器１００Ａ（あるいは路面状態推定器１１０Ａ）は、前方路面の路面状態を検出する前方路面状態検出手段を構成する。タイヤ発生音検出センサ２０、及び路面状態推定器２１０は、現在走行している現在走行路面の路面状態を検出する現在路面状態検出手段を構成する。遅れ時間処理器３０は、同期手段を構成する。そして、路面状態推定器２１０は、報知手段を構成する。
【０３１０】
次に、本実施の形態の作用を説明する。
【０３１１】
上記路面状態推定器２１０は、現在走行中の路面の路面状態を推定するのに、タイヤ発生音の特徴量から推定するほかに、路面状態推定器１００Ａ（あるいは路面状態推定器１１０Ａ）で推定された前方路面の路面状態の推定結果を、遅れ時間処理器３０で現在走行中の路面の時刻に同期された推定結果も参考にして、現在走行中の路面の路面状態を推定する。
【０３１２】
その具体的な一方法は、路面状態推定器１００Ａ（あるいは路面状態推定器１１０Ａ）で推定された前方路面の路面状態の推定結果が遅れ時間処理器３０で現在走行中の路面の時刻に同期されたその推定結果が、圧雪あるいは氷盤の第１の路面状態で、タイヤ発生音の特徴ベクトルから推定された現在路面の路面状態が前回の圧雪または氷盤の路面の第１の路面状態から、乾燥または湿潤路面の第２の路面状態に変化しても、乾燥あるいは湿潤路面状態の確率は低いとして、最終的に路面状態推定器２１０は圧雪または氷盤の第１の路面状態を推定結果とする。
【０３１３】
以上説明したように第８の実施の形態では、現在走行中の路面状態を推定する際、前方路面状態の推定結果を参考にすることにより、ノイズ等の誤判定からくる確率の低い現在路面状態推定結果を排除することにより、推定の精度を向上できる。
【０３１４】
次に、第１１の実施の形態を説明する。
【０３１５】
図３０に示すように、本実施の形態にかかる路面状態推定装置は、タイヤ発生音検出センサ２０Ａ、車両が走行する路面の状態を推定する路面状態推定器２００（あるいは路面状態推定器１００Ａ，１１０Ａ、２１０）、路面と車両のタイヤとの間の摩擦状態である路面摩擦状態としての路面摩擦係数を推定する路面摩擦状態推定手段としての路面摩擦係数推定器３００、及び、推定された路面状態と、推定された路面摩擦係数と、に基づいて、最終の路面摩擦係数を推定する最終路面摩擦状態推定手段としての路面摩擦係数学習器４０Ａを備えている。
【０３１６】
なお、タイヤ発生音検出センサ２０Ａ、及び路面状態推定器２００（あるいは路面状態推定器１００Ａ，１１０Ａ、２１０）は、路面状態検出手段を構成する。
【０３１７】
次に、本実施の形態を説明する。
【０３１８】
本実施の形態にかかる路面摩擦係数推定器３００は、車両の車輪速、前後加速度、横加速度から走行中の路面の摩擦係数を推定する。なお、本実施の形態にかかる路面摩擦係数推定器３００は、例えば、特願２００２−０４０７２２号に記載の装置が適応可能である。
【０３１９】
路面摩擦係数学習器４０Ａは、路面状態推定器２００（あるいは路面状態推定器１００Ａ，１１０Ａ、２１０）で推定された各路面状態の、路面摩擦係数推定器３００で推定された摩擦係数を記憶しておき、走行中の各路面状態の最も確からしい摩擦係数を学習して求める。一例として最も簡単な確からしい摩擦係数を求める方法は、その路面状態に対応する路面摩擦係数推定器３００により推定された摩擦係数の移動平均値をその路面状態の摩擦係数として出力することである。以上のように、路面摩擦係数学習器４０は、路面状態推定器２００（あるいは路面状態推定器１００Ａ，１１０Ａ、２１０）で推定された路面状態における、最も確からしい摩擦係数を推定することが出来る。
【０３２０】
以上説明したように第５の発明では、走行中の路面状態の摩擦係数を、常にオンラインで学習により求めているので、地域、季節等による路面状態の変化にも適応して、適切な路面摩擦係数を推定することが出来る。
【０３２１】
以上説明した実施の形態では、路面摩擦係数学習器４０は、路面摩擦係数推定器３００により推定された全ての摩擦係数を平均しているが、本発明はこれに限定されず、図３１に示すように、路面摩擦係数推定器３００と、路面摩擦係数学習器４０と、の間に、路面摩擦係数選別器５０Ａを備えるようにしてもよい。
【０３２２】
即ち、路面摩擦係数推定器３００により推定される路面摩擦係数は常に正しい値でない場合があり、何らかの影響により、現在の路面状態からかけ離れた値を推定される場合もある。このように現在の路面状態からかけ離れた値が推定されても、この値を含めて路面摩擦係数推定器３００により推定された全ての摩擦係数を平均してしまったのでは、最終敵な路面摩擦係数の推定値が正確な値とならない。
【０３２３】
そこで、路面摩擦係数選別器５０Ａでは、複数の路面状態の各々毎に、各路面状態に対応する路面摩擦係数の値の範囲を予め定めておき、各路面状態に対応する路面摩擦状態の範囲、路面状態推定器２００等により推定された路面状態、及び路面摩擦係数推定器３００により検出された路面摩擦係数に基づいて、路面摩擦係数推定器３００により検出された路面摩擦係数が、路面状態推定器２００により推定された路面状態に対応する路面摩擦係数の値の範囲内か否かを判断し、路面摩擦係数推定器３００により検出された路面摩擦係数が、路面状態推定器２００により推定された路面状態に対応する路面摩擦係数の値の範囲内と判断された場合に、路面摩擦係数推定器３００により検出された該路面摩擦係数を、路面摩擦係数学習器４０に出力するようにする。
【０３２４】
このように、路面摩擦係数推定器３００により検出された路面摩擦係数が、路面状態推定器２００により推定された路面状態に対応する路面摩擦係数の値の範囲内と判断された場合に、路面摩擦係数推定器３００により検出された該路面摩擦係数を、路面摩擦係数学習器４０に出力するので、現在の路面状態からかけ離れた値を用いて、路面摩擦係数学習器４０が最終敵な路面摩擦係数を推定することを防止できる。よって、最終敵な路面摩擦係数の推定値をより正確な値とすることができる。
【０３２５】
次に、第１２の実施の形態を説明する。
【０３２６】
図３２に示すように、本実施の形態にかかる路面状態推定装置は、図示しないが、上記第１の実施の形態における特徴ベクトル生成器１６（図１参照）により生成されたタイヤ発生音特徴量と、車輪速センサ１２（図１参照）により検出された車輪速と、から、現在走行中の路面状態を推定する路面状態推定器２２０を備えている。そして、路面状態推定器２２０は判定処理器２１Ａとプロトタイプ演算器２２Ａとからなる。
【０３２７】
ここで、タイヤ発生音特徴量と車輪速の関係を説明する。タイヤ発生音特徴量の値は車輪速が高くなると、それに応じて大きくなる。図３６は、車輪速に対するタイヤ発生音量の関係を示したものである。タイヤ発生音特徴量の大きさは車輪速の２次関数で近似できる。
【０３２８】
したがって、プロトタイプ演算器２２Ａでは、走行中検出されたタイヤ発生音特徴ベクトルと照合するためのプロトタイプの値を、車輪速の２次関数から演算して求める。判定処理器２１では、走行中検出されたタイヤ発生音特徴ベクトルと、プロトタイプ演算器２２で演算されたプロトタイプの値を照合して、現在走行中の路面状態を推定する。
【０３２９】
以上説明したように第１２の実施の形態では、複数の車輪速毎に対応して記憶された路面状態のプロトタイプをすべて持つのではなく、２次関数の係数のみ記憶すればよいので、記憶容量を削減できる。
【０３３０】
次に、第１３の実施の形態を説明する。
【０３３１】
図３３に示すように、本実施の形態にかかる路面状態推定装置は、タイヤ発生音特徴量と車輪速から現在走行中の路面状態を推定する路面状態推定器２３０を備えている。そして、路面状態推定器２３０は判定処理器２１Ａと、プロトタイプ学習器２３Ａと、プロトタイプ演算器２４Ａとからなる。
【０３３２】
プロトタイプ学習器２３Ａは、走行中オンラインで代表点の車輪速に対するタイヤ発生音特徴量を記憶し、記憶した値をもとにタイヤ発生音特徴量の車輪速に関する２次関数の係数を演算し更新する。プロトタイプ演算器２４Ａでは、プロトタイプ学習器２３Ａで求まる２次関数の係数をもとに、プロトタイプ演算器２４Ａでは、走行中検出されたタイヤ発生音特徴ベクトルと照合するためのプロトタイプの値を、車輪速の２次関数から演算して求める。判定処理器２１Ａでは、走行中検出されたタイヤ発生音特徴ベクトルと、プロトタイプ演算器２４Ａで演算されたプロトタイプの値を照合して、現在走行中の路面状態を推定する。
【０３３３】
以上説明したように第１３の実施の形態では、タイヤ発生音特徴ベクトルと照合するためのプロトタイプを演算するためのプロトタイプの２次関数の係数を、走行中常にオンラインで学習して求めるため、地域、季節等による路面状態の変化にも、たとえば同じ乾燥路面でも路面状態の性質が異なり、タイヤ発生音特徴量が変化しても、それに適応して、プロトタイプ点を適切に対応できるため、変化に適応して路面摩擦係数を精度良く推定することが出来る。
【０３３４】
以上説明した種々の実施の形態では、タイヤ発生音検出センサ２０Ａを、タイヤ発生音に、該タイヤ発生音の発生源に基づくタイヤ発生音以外の外乱の影響が小さくなる位置、例えば、タイヤ最近傍のバンパー内部に設置する。
【０３３５】
タイヤ近傍でタイヤ発生音を検出する場合、タイヤで跳ね上げられる、泥、水、などの障害物、また風圧などの音の外乱を回避する必要がある。タイヤ最近傍のバンパー内部では泥、水、などの障害物、また風圧などの音の外乱を回避でき、しかもタイヤ発生音検出センサ２０を比較的容易に設置できる。
【０３３６】
以上説明したように、タイヤ発生音検出センサ２０を、タイヤ最近傍のバンパー内部に設置することにより、泥、水、などの障害物、また風圧などの音の外乱を回避でき、ノイズの少ないタイヤ発生音を検出できる。またその設置が容易なことより、低コストで実現できる。
【０３３７】
また、路面側の画像を検出する上記種々の実施の形態では、路面画像検出センサ１０Ａにより検出された画像データをそのまま出力するようにしているが、本発明はこれに限定されるものではなく、次の各素子を備えるようにしてもよい。
【０３３８】
即ち、路面画像検出センサ１０Ａ、路面輝度検出器１４Ａ、ボンネット輝度検出器１５Ａ、及び輝度比演算器１６Ａを備えるようにしてもよい。
【０３３９】
路面画像検出センサ１０Ａは、車両の一部及び前方路面を含む領域を撮像する。路面輝度検出器１４Ａは、路面画像検出センサ１０Ａで検出された前方路面画像から一部の前方路面の輝度を検出する、ボンネット輝度検出器１５Ａは、路面画像検出センサ１０Ａで検出された前方路面画像から車両の一部のボンネット部の輝度を検出する、輝度比演算器１６Ａは、路面輝度検出器１４Ａで得られた前方路面の輝度と、ボンネット輝度検出器１５Ａで得られた車両のボンネット部の輝度の比を演算して、路面状態推定器１００（あるいは路面状態推定器１１０）に出力する。
【０３４０】
このように、前方路面状態を推定するための前方路面輝度情報として、車両のボンネット部の輝度との比を用いているため、車両の向きとともに日のあたる角度が変化することによる輝度の変化の影響を低減できる。また前方画像全体でなく、ある部分の平均輝度を用いることにより、演算負荷を低減できる。
【図面の簡単な説明】
【０３４１】
【図１】第１の実施の形態に係る路面状態推定装置のブロック図である。
【図２】スタッドレスタイヤで異なる路面を所定速度で走行した時の時間−周波数解析結果を示す図である。
【図３】（Ａ）は乾燥アスファルト路、（Ｂ）は湿潤アスファルト路の、それぞれ予め把握された各路面のタイヤ発生音を、ドベシィのウェーブレット（Ｎ＝３）を用いて、ウェーブレット変換処理した波形を示す図である。
【図４】（Ａ）は永盤路、（Ｂ）は圧雪路の、それぞれ予め把握された各路面のタイヤ発生音を、ドベシィのウェーブレット（Ｎ＝３）を用いて、ウェーブレット変換処理した図である。
【図５】各ウェーブレット成分を時間−周波数解析した図である。
【図６】ウェーブレット変換を用いた前処理器１４のアルゴリズムを示すフローチャートである。
【図７】第２の実施の形態にかかる、路面摩擦係数を推定する路面摩擦係数推定装置のブロック図である。
【図８】第３の実施の形態にかかる、路面摩擦係数を推定する路面摩擦係数推定装置のブロック図である。
【図９】第４の実施の形態にかかる、路面摩擦係数を推定する路面摩擦係数推定装置のブロック図である。
【図１０】第５の実施の形態にかかる、路面摩擦係数を推定する路面摩擦係数推定装置のブロック図である。
【図１１】第６の実施の形態にかかる、路面摩擦係数を推定する路面摩擦係数推定装置のブロック図である。
【図１２】ニューラルネットワークで構成した場合のニュー判定処理器の構成の概念図である。
【図１３】４つの異なる路面を走行したときのタイヤ発生音を、帯域フィルタを通過させることにより、分解した波形を示した図である。
【図１４】（Ａ）及び（Ｂ）は、４つの異なる路面を走行したときのタイヤ発生音を、帯域フィルタを通過させることにより分解した波形から求めた特徴量（波形の実効値）を示した図である。
【図１５】ドベシイのウェーブレットの基本波形を示した図である。
【図１６】シムレットのウェーブレットの基本波形を示した図である。
【図１７】コアレットのウェーブレットの基本波形を示した図である。
【図１８】メイエのウェーブレットの基本波形を示した図である。
【図１９】（Ａ）及び（Ｂ）は、４つの異なる路面を走行したときのタイヤ発生音を、ドベシイのウェーブレットの基本波形で分解した波形から求めた特徴量（相関値）を示した図である。
【図２０】前処理器の詳細を示す第１の実施の形態に係る路面状態推定装置のブロック図である。
【図２１】第１の実施の形態に係る路面状態推定装置の第１の変形例のブロック図である。
【図２２】第１の実施の形態に係る路面状態推定装置の第２の変形例のブロック図である。
【図２３】タイヤ発生音と相関をとるために路面毎に対応して予め記憶された基本波形であり、（Ａ）は、乾燥路面の基本波形を、（Ｂ）は、湿潤路面の基本波形を示した図である。
【図２４】タイヤ発生音と基本波形との相関を説明する説明図である。
【図２５】基本波形の設定の内容を説明する説明図である。
【図２６】第７の実施の形態に係る路面状態推定装置のブロック図である。
【図２７】第８の実施の形態に係る路面状態推定装置のブロック図である。
【図２８】第９の実施の形態に係る路面状態推定装置のブロック図である。
【図２９】第１０の実施の形態に係る路面状態推定装置のブロック図である。
【図３０】第１１の実施の形態に係る路面摩擦状態推定装置のブロック図である。
【図３１】第１１の実施の形態に係る路面摩擦状態推定装置の変形例のブロック図である。
【図３２】第１２の実施の形態に係る路面状態推定装置のブロック図である。
【図３３】第１３の実施の形態に係る路面状態推定装置のブロック図である。
【図３４】タイヤ発生音検出センサの配置位置を示した図である。
【図３５】路面輝度情報出力部のブロックである。
【図３６】車輪速とタイヤ発生音特徴量との関係を示したグラフである。
【符号の説明】
【０３４２】
１０、１０Ａ  タイヤ発生音検出センサ（タイヤ発生音検出手段）
１２、１２Ａ  車輪速検出センサ（車輪速検出手段）
１４、１４Ｎ１、１４Ｎ２  前処理器（計算手段）
１６  特徴ベクトル生成器
２０、２０Ａ  判定処理器（推定手段）
２２  判定マップメモリ（記憶手段）
１０２    前方路面特徴量検出器（検出手段、第２の検出手段）
１１２    μ判定マップメモリ（記憶手段）
１０４    現在路面特徴量検出器（検出手段、第１の検出手段）
１２６    μ判定マップメモリ（記憶手段）
１０７    路面摩擦係数真値記憶装置（記憶手段）
１０６    遅れ時間処理器（同期手段）
１１０    μ判定処理器（推定手段）
１２４    μ判定処理器（推定手段）
１１８    マップ修正処理器（修正手段）
１３２    マップ修正処理器（修正手段）
【図１】
2004168286.tif 000031
【図２】
2004168286.tif 000032
【図３】
2004168286.tif 000033
【図４】
2004168286.tif 000034
【図５】
2004168286.tif 000035
【図６】
2004168286.tif 000036
【図７】
2004168286.tif 000037
【図８】
2004168286.tif 000038
【図９】
2004168286.tif 000039
【図１０】
2004168286.tif 000040
【図１１】
2004168286.tif 000041
【図１２】
2004168286.tif 000042
【図１３】
2004168286.tif 000043
【図１４】
2004168286.tif 000044
【図１５】
2004168286.tif 000045
【図１６】
2004168286.tif 000046
【図１７】
2004168286.tif 000047
【図１８】
2004168286.tif 000048
【図１９】
2004168286.tif 000049
【図２０】
2004168286.tif 000050
【図２１】
2004168286.tif 000051
【図２２】
2004168286.tif 000052
【図２３】
2004168286.tif 000053
【図２４】
2004168286.tif 000054
【図２５】
2004168286.tif 000055
【図２６】
2004168286.tif 000056
【図２７】
2004168286.tif 000057
【図２８】
2004168286.tif 000058
【図２９】
2004168286.tif 000059
【図３０】
2004168286.tif 000060
【図３１】
2004168286.tif 000061
【図３２】
2004168286.tif 000062
【図３３】
2004168286.tif 000063
【図３４】
2004168286.tif 000064
【図３５】
2004168286.tif 000065
【図３６】
2004168286.tif 000066