(57)【特許請求の範囲】
【請求項１】
車両のタイヤ（１〜４）と路面における摩擦力に基づく変動であって前記タイヤの回転と同期した変動成分を含む車輪速度信号を出力する車輪速度センサ（１１〜１４）と、
前記車輪速度信号から車輪速度信号成分とピッチ成分とを分離することで前記ピッチ成分として表される前記変動成分を抽出する手段とを備え、
前記抽出された変動成分に基づき路面状態を判別することを特徴とする路面状態判別装置。
【請求項２】
前記信号をケプストラム解析して前記変動成分の抽出を行うことを特徴とする請求項１に記載の路面状態判別装置。
【請求項３】
車両のタイヤ（１〜４）と路面における摩擦力に基づく変動であって前記タイヤの回転と同期した変動成分を含む信号を出力する手段（１１〜１４）を備え、
ケプストラムの原理に基づき前記信号から前記変動成分を抽出し、この変動成分に基づき路面状態を判別することを特徴とする路面状態判別装置。
【請求項４】
前記信号を出力する手段は、車輪速センサ（１１〜１４）であることを特徴とする請求項３に記載の路面状態判別装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、走行中の路面状態を判別する路面状態判別装置に関し、ドライバーへの危険情報の通知を目的とするシステムに関するもの、アンチスキッドブレーキ装置、トラクションコントロール装置等の走行系の制御へ路面情報を伝達することにより車両走行の安全性を向上するシステムに適用することもできる。
【０００２】
【従来の技術】
従来、路面における摩擦係数μ（以下、μという）の判別方法として、車両の減速・加速状態におけるスリップ率に基づいて判別を行う方法や、旋回中の横加速度、回転モーメントに基づいて判別を行う方法が知られている。
また、検出された摩擦係数μに基づいて制御が行われる装置として、アンチスキッドブレーキ装置やトラクションコントロール装置がある。
【０００３】
【発明が解決しようとする課題】
しかしながら、これらの路面のμの判別方法によると、車両に加速度が発生していない一定速度状態であるとき、或いはステアリングによる操作が行われていないときには、路面のμの判別を行うことができない。
従って、上記したような状態においては、路面のμが低下したという情報等を運転者に伝達して注意を促したり、アンチスキッドブレーキ装置等の制御装置に早いタイミングで伝達することができなかった。
【０００４】
本発明は上記点に鑑みて、車両に加速度が発生していない一定速度状態であるとき、或いはステアリングによる操作が行われていないときにおいても路面のμを判別可能にすることを目的とする。
【０００５】
【課題を解決するための手段】
円形のタイヤは、接地面で一度フラットな状態に変形するが、この時タイヤトレッド面に働く前後力は路面のμやタイヤ周上の各部の材料の強弱によって変化する。図１（ａ）は接地面でのゴムブロックが前後方向に受ける力を図示したものであるが、転動しているタイヤ１００のゴムブロックは最初の接地部で前方向に引っ張られたのち、上下に圧縮され地面から離れようとするときに後部に引っ張られる。この力をΔＦとする。
【０００６】
前後に引っ張られる力ΔＦは路面のμ特性及びタイヤ１００のゴムブロックの剛性により異なる。一般には、路面のμまたはタイヤ剛性が高いほど、力ΔＦの変動は大きくなる。また、路面のμが高く、剛性が高いほど或いはタイヤの剛性が高いほど力ΔＦの変化率（勾配）も大きい。
ところで、一般にタイヤの各部の材料は回転中の共振を減らすため、材料の剛性に強弱をつけて製作される。このため、接地部にてタイヤ各部が異なる状態で変形してタイヤ回転毎に変動が起きる。また、タイヤ１００のトレッド部及びカーカス部に振動や強弱の波を形成しながら転動する。このため、路面μが同様である場合にはタイヤ１００の同一箇所で、つまりタイヤ１００の回転と同一又はその整数倍の周波数毎に同様な波を形成する。
【０００７】
このため、路面のμの大小及び路面の均一性は、上記タイヤ１００の振動や剛性の強弱によってタイヤ１００の回転やサスペンション系に発生するタイヤ１００の回転と同期する変動（以下、タイヤ回転ピッチという）に基づいて計測することができる。つまり、タイヤ１００の同一箇所に発生する変動を物差しとしてその大小、継続性を評価することで路面状態を判別することができる。
【０００８】
具体的には、人工的に形成されたアスファルトのような路面では路面状態が安定して継続する。このため、タイヤ回転ピッチの情報が乱されない。また、路面μが大きいために力ΔＦの変動が大きく、図１（ｂ）に示すようにタイヤ１００の前後歪みが大きくなり、タイヤ回転ピッチも大きくなる。
しかし、自然的に形成される低μ路面では路面のμが小さく、自然路面であるためバラツキがあり、また凹凸もアスファルトより不均一である。このため、タイヤ回転ピッチの情報が乱される。また路面μが小さいために力ΔＦの変動が小さく、図１（ｃ）に示すようにタイヤ１００の前後歪みが小さくなり、タイヤ回転ピッチ成分も小さくなる。
【０００９】
上記理由により車両の走行中、タイヤ回転ピッチを計測することにより路面のμの状態を判別することができる。即ち、アスファルト路のようにグリップがあり一定した路面と、雪や氷のようにμが小さくバラツキがある路面をタイヤ回転ピッチを計測することにより区別できる。
上記目的を達成するために、以下に示す手段を採用する。
【００１０】
請求項１に記載の発明においては、車両のタイヤ（１〜４）と路面の間の摩擦力に基づく変動であってタイヤの回転と同期した変動成分を含む車輪速度信号を出力する車輪速度センサ（１１〜１４）と、この車輪速度信号から速度信号成分とピッチ成分とを分離することでピッチ成分として表される変動成分を抽出する手段とを備え、この変動成分に基づき路面状態を判別することを特徴とする。このように、車輪速センサからの車輪速度信号から車輪速度信号成分とタイヤ回転ピッチ成分とを分離して、分離されたタイヤ回転ピッチ成分から路面のμの判別を行う。
【００１１】
上述したように、タイヤ回転ピッチ成分から路面のμを判別することができる。従って、路面判別をタイヤ回転ピッチで行うことにより車両にスリップが発生していない一定速度状態であるとき、或いはステアリングによる操作が行われていないときにおいても、路面のμを判別することができる。
具体的には、請求項２に示すように、タイヤ回転ピッチの抽出は、ケプストラム解析によって行うことができる。この場合、タイヤ回転ピッチ成分に基づき、路面のμの判別を行うことができる。
【００１２】
また、請求項３に記載の発明では、車両のタイヤ（１〜４）と路面における摩擦力に基づく変動であってタイヤの回転と同期した変動成分を含む信号を出力する手段（１１〜１４）を備え、その信号からケプストラムの原理に基づき変動成分を抽出し、この変動成分に基づき路面状態を判別することを特徴としている。
このように、タイヤの回転と同期した変動成分を含む信号をケプストラム解析することで変動成分を抽出できる。そして、この変動成分に基づいて路面状態を判別することが可能である。なお、変動成分を含む信号としては、車輪加速度、ホイールトルク、サスペンション歪み又はホイールサスペンション振動等がある。
【００１３】
【発明の実施の形態】
以下、本発明を図に示す実施形態について説明する。
図２に、本実施形態における構成図を示す。図２に示すように各車輪１〜４に対応して車輪速センサ１１〜１４が取りつけられており、これらの車輪速センサ１１〜１４からの信号が路面状態判別用の電子制御装置（以下、ＥＣＵという）５に入力される。また、ＥＣＵ５には、ステアリング６に取りつけられたステアリングセンサ７により検出されるステアリング角も入力され、これらに基づき路面状態判別が行われる。なお、ＥＣＵ５によって判別された路面状態は、ウォーニングランプ８にて運転者に警告される。
【００１４】
具体的には、上記路面状態判別は、ケプストラム解析を使用してタイヤ回転ピッチ成分を検出することにより行う。
ここで、ケプストラム解析とは、音声認識等で普及している技術であり（参考文献：斉藤  収三著  音声情報処理の基礎  オーム社）、音声においてケプストラムの原理に基づき各ピッチ成分を検出している。このケプストラムの原理に基づき、タイヤ回転数ピッチの成分を検出する。
【００１５】
即ち、信号ｓ（ｔ）を周期Ｔのピッチ周期の信号とみなした場合、数式１のように、信号ｓ（ｔ）は信号成分ｓ０（ｔ）とピッチ成分Σδの積算値の畳み込み積分で表せる。
【００１６】
【数１】
0003684726.tif 000002
【００１７】
また、数式１をフーリエ変換（以下、ＦＦＴという）することにより、ｓ（ｔ）におけるパワースペクトル値：Ｓ（ω）が信号成分Ｓ０（ω）とピッチ成分の積で表せ、数式２のようになる。
【００１８】
【数２】
0003684726.tif 000003
【００１９】
さらに、数式２について対数変換を行うと信号成分とピッチ成分が分離でき、数式３で表される。
【００２０】
【数３】
0003684726.tif 000004
【００２１】
但し、Ｓ（ｔ）：入力信号、ｔ：時間、ｎ：整数、Ｔ：ピッチ周期、Ｓ０（ｔ）：ピッチ以外の信号成分、＊：畳み込み積分
この数式３における第１項はｓ０（ｔ）のパワースペクトル値であり、第２項の括弧内は線スペクトルを示す。ここで、このｓ０（ｔ）のパワースペクトル値はωに対して比較的ゆるやかな変化を示し、線スペクトルはωに対して急激な変化を示す。
【００２２】
このため、数式３を変数ωの信号列と見なしたとき、ピッチ成分は信号成分に比べてωに対して急激な変化を示す。また、一般には、ｓ０（ω）はωに対して急な変化は少ないため数式３を逆フーリエ変換するとピッチ成分を高周波（高ケフレンシー）側に、信号成分を低周波（低ケフレンシー）に分離することができる。
【００２３】
つまり、Ｓ０（ω）は、低周波成分が主体であるので高周波領域ではあまり波形として現れず、この領域ではタイヤピッチ成分を主体とした波形になるため、逆ＦＦＴを行うことにより、低周波の車輪速度信号成分と高周波のタイヤピッチ成分とを分離し、この波形からタイヤ回転ピッチ成分のみを抽出することができる。
【００２４】
従って、車両においてはケプストラムにて車輪速度信号からタイヤ回転ピッチ成分及び車輪速度信号成分に分離して、このうちのタイヤ回転ピッチ成分から路面状態判別を行うことができる。
また、路面がドライ状態、圧雪状態、氷上状態である場合において、検出される波形の例を図３（ａ）〜（ｃ）に示す。図３（ａ）に示されるように、ドライ状態においてはタイヤ回転ピッチ成分はタイヤ１回転周期となる波形として現れ、タイヤ回転ピッチ成分は大きく現れる。図３（ｂ）に示されるように、圧雪状態においては凹凸や路面μが不均一であるため振動がランダムになり無秩序な波形となる。図３（ｃ）に示されるように、氷上状態においてはタイヤ１回転周期となる波形として現れるが、タイヤ回転ピッチ成分は小さくしか現れない。このように、各状態で特徴ある波形となり、これらの波形に基づきケプストラム解析を行いタイヤ回転ピッチ成分を検出する。
【００２５】
上述したケプストラム解析を使用した路面判別検出における各処理を、図４に示す。また、図５にこれらの処理に基づき演算される波形を示す。
まず、ステップ２０１では、回転角度センサからローターの歯車毎（１パルス毎）の車輪速度パルスを検出して、この車輪速度パルスから車輪速度パルス間隔を求めて入力する。なお、この車輪速度パルスの検出には汎用に普及している電磁ピックアップ方式又はフォトカプラ方式等が利用できる。
【００２６】
ステップ２０２では、サンプルした車輪速度パルス間隔から、サンプル間隔が一定時間となる車輪速度Ｖｗ（ｔ）を求める。具体的には、図６（ａ）に示すように、ローター５０の歯車５０ａ毎に発生する車輪速度パルス毎に、図６（ｂ）に示すような車輪速度Ｖｗ（ｎ）を求め、さらに車輪速度Ｖｗ（ｎ）から２次スプライン関数で補間して、図６（ｃ）に示すサンプリング間隔が一定時間となる車輪速度Ｖｗ（ｔ）を求める。
【００２７】
続くステップ２０３では、ＦＦＴを実行する際の誤差を減らすためにハイパスフィルター（以下、ＨＰＦという）によって速度の低周波成分を除いている。
ステップ２０４〜ステップ２０８は、ケプストラム演算における処理である。上記ケプストラムの原理に基づき、これらのステップについて説明する。
ステップ２０４でウィンドウ処理を行い、対象となるウィンドウエリアを決める。このエリアが狭いと計算精度が下がり、広いと応答性が下がるため適切な広さのエリアに設定する必要がある。なお、計算範囲をタイヤ回転速度に応じて切り替えることにより低速から高速までカバーできる。
【００２８】
本実施形態では、このウィンドウエリアをＴｓ＝２πｒ／Ｖｂ＊Ｎｒｏｔ（但し、Ｔｓ：ウィンドウ時間、ｒ：タイヤ半径、Ｖｂ：車体速度、Ｎｒｏｔ：タイヤ回転数とする）に基づき決定している。基本的には計測エリアをタイヤが何回転するかで決定するものであり、Ｎｒｏｔは例えば約２０回転ほどで実施する。
なお、図５（ａ）に示す波形は、検出された車輪速度Ｖｗ（ｔ）を示す波形であり、図５（ｂ）に示す波形は、ＨＰＦの処理後にウィンドウ処理によって決定されたウィンドウエリア内の車輪速度サンプルを拡大した波形である。そして、この図５（ｂ）における波形が、数式１における信号ｓ（ｔ）によって表された波形としており、これにより信号ｓ（ｔ）が決定する。
【００２９】
ステップ２０５では、ウィンドウエリア内において、上記した数式２に示すようなＦＦＴ演算を行いＳ（ω）を算出する。ステップ２０６ではＦＦＴ演算の結果からパワースペクトル値を演算する。続くステップ２０７では、数式３に示すように、このパワースペクトル値の対数値ｌｏｇＶｗ（ω）を演算する。なお、この時のｌｏｇＶｗ（ω）をプロットすると図５（ｃ）の波形にて表される。そして、ステップ２０８では、この対数値を逆ＦＦＴ演算する。
【００３０】
ステップ２０９及びステップ２１０では、タイヤ回転ピッチ成分と関係した成分のみを取り出すために、ケフレンシー軸上で車体速度と関係したウィンドウをかける処理を行う。まず、ステップ２０９では、ＨＰＦによってノイズ等で発生したオフセット成分を取り除く。このとき、オフセット成分を取り除いた後の波形は図５（ｄ）の波形にて表される。
【００３１】
そして、ステップ２１０では、車輪速度信号成分を取り除きタイヤの回転ピッチのみを検出するためのウィンドウ処理を行う（図５（ｅ）の波形に該当する）。ウィンドウはタイヤの回転周期にあわせて一定の間隔でｎ個並べられ、その配置間隔ΔＴはΔＴ＝２πｒ／Ｖｂによって算出される。また、ｎはケプストラムグラフ上に発生するピッチの次数成分で決められるが、通常のタイヤでは４次程であるため、本実施形態においてはｎ＝４で処理を行っている。また、ウィンドウにはハミングウィンドウ：Ｈｗ（ｉ，ｑ）を使用しており、数式４のような関数として表すことができる。
【００３２】
また、一般に高μ路面の時には高次成分がでやすいため、この場合には数式５に示すように、ウィンドウゲイン：Ｇ（ｉ）を高く設定する。
なお、所定の定数はハミングウィンドウに用いられる所定の定数であり、ｑはケフレンシー、Ｍはウィンドウ幅である。
【００３３】
【数４】
0003684726.tif 000005
【００３４】
【数５】
0003684726.tif 000006
【００３５】
また、上記数式４、５に基づきウィンドウ処理後のケプストラム値：Ｃｗ（ｑ）は数式６により求められる（図５（ｆ）の波形に該当する）。なお、Ｃ（ｑ）はオフセットをＨＰＦで取り除いた後のケプストラム値である。
【００３６】
【数６】
0003684726.tif 000007
【００３７】
そして、ステップ２１１では、ウィンドウ領域における２乗和を算出して、タイヤ回転ピッチの各次数成分の総和（タイヤ回転ピッチパラメータ）：Ｃｐを求める。
【００３８】
【数７】
0003684726.tif 000008
【００３９】
なお、実験的に設定された各種路面状態について、図７（ａ）に示すような速度で走行した場合における各々のケプストラム解析結果Ｃｐを図７（ｂ）に示す。このように、路面の種類に応じてケプストラム解析の結果が異なる。このため、ケプストラム解析によって路面状態の判別を行うことができる。
ステップ２１２では、先程求めたタイヤ回転ピッチ成分と車体加速度及びステアリング角をもとに路面判別を実施する。つまり、タイヤ回転ピッチ成分はステアリング角や車体加速度によって変化するため、より完全に路面状態を判別するためである。
【００４０】
具体的には、この判別は路面の状態に基づき設定された図８に示す路面状態選択マップに基づいて行い、タイヤ回転ピッチパラメータＣｐ、ステアリング角及び車体加速度により表される座標Ａ（タイヤ回転ピッチパターン，ステアリング角，車体加速度）が図８に示される座標軸上のどこに位置するかによって路面判別を行う。なお、車体加速度は、４輪の車輪速度に基づいて算出している。
【００４１】
具体的には、図８に示されるようにアスファルトやコンクリート等の高μ路と雪道や氷上等の低μ路とを区別するしきい面Ｂがあり、座標Ａがこのしきい面Ｂの上側にあれば高μ路と判別され、逆に下側にあれば低μ路と判別される。
このとき、運転者がステアリングを操作していない場合にはステアリング角は零であり、また車速が一定である場合には車体加速度は零である。このような場合には、座標Ａは（Ｃｐ，０，０）となる。
【００４２】
ステップ２１３では、判別結果に基づき運転者への通知を行う。通知には、路面μ低下を示すウォーニングランプを使用し、例えば高μ路から低μ路に移行した場合に数秒間ウォーニングランプが点灯させる等して通知する。また、アンチロックブレーキ装置やトラクションコントロール装置等に対しても予め低μ路である旨の信号を送ることにより、これらの制御初期におけるスリップを有効に抑えることができる。
【００４３】
このように、車両に加速度が発生していない一定速度状態であるとき、或いはステアリングによる操作が行われていないときにおいても路面の摩擦係数μを検出することができる。
（第２の実施形態）
本実施形態においては、第１実施形態と異なるサンプル間隔でケプストラム解析を行う。
【００４４】
ケプストラム解析は音声処理に使用する都合上通常は時間軸上で解析を行うため、第１実施形態においてはこの時間軸上で上記解析を行ったが、本実施形態においては距離軸上で解析を行う。この場合、タイヤ回転ピッチは、タイヤの回転する距離に比例してピッチが発生するため、距離軸上で解析を行った方が直接的な解析を行うことができる。
【００４５】
具体的には、図５（ｂ）に示すように、第１実施形態においては車輪速度におけるサンプリング間隔を一定時間の周期毎にしているため、車輪速度Ｖｗ（ｔ）をローター５０の歯車５０ａ毎に取り出された車輪速度Ｖｗ（ｎ）から補間する必要がある。しかし、図５（ｃ）に示すように、車輪速度におけるサンプリング間隔を一定距離の周期毎にすると、距離はローター５０の歯車５０ａの移動に比例しているため、サンプリング間隔をローター５０の歯車５０ａ毎（１パルス毎）にすることができ、上記補間の必要がなくなる。
【００４６】
このように距離に比例した解析を行うことによりケプストラム値もタイヤ回転に同期した距離の位置に並ぶ。ケプストラム解析を使用した路面判別検出における各処理を図９に示す。以下、図９に基づきこれらの処理を説明する。
ステップ３０１にて、１個分のパルスが伝達されるのに所要時間を１周期として、この周期毎の車輪速度に対応した車輪速度パルスを検出する。次に、ステップ３０２では、車輪速度パルスに基づき、Ｖ＝ｋ×Δｎ／ΔＴの関係から車輪速度Ｖｗ（ｎ）を演算する。
【００４７】
また、ステップ３０３〜ステップ３０９については、第１実施形態におけるステップ２０３〜ステップ２０９と同様の処理を行う（図１０（ａ）の波形に該当する）。
そして、ステップ３１０では、タイヤ回転ピッチを検出するためのウィンドウ処理を行う（図１０（ｂ）の波形に該当する）。但し、この際には算出されたケフレンシーの単位がパルス数となるため、タイヤ１回転の検出パルス数（歯数）をＮｓとすると、このパルス数Ｎｓの定数倍でウィンドウ位置を決定する。
【００４８】
この後、ステップ３１１〜ステップ３１３について、第１実施形態におけるステップ２１１〜ステップ２１３と同様の処理を行い処理を終了する。
なお、このように距離軸上でケプストラム解析を行っているため回転と同期していないノイズを除くという効果が得られる。つまり、距離軸上でケプストラム解析を行う場合、サンプリング時に取り出される信号は、回転と同期した成分を持っている。このため、その信号における最大周波数は、回転周期に比例する周波数となる。従って、サンプル周波数を回転に比例した値とすれば、サンプリング時に取り出される信号と一定の比例関係を保つことができ、車速によらずエリアスが発生しにくくなる。
【００４９】
（他の実施形態）
第１、第２実施形態においては、タイヤ回転ピッチ成分の検出を、車輪速センサからの信号（車輪速度信号）に基づき行ったが、タイヤと路面の摩擦力に基づく変動を信号出力できるもの、例えば上下加速度センサやホイールトルクセンサ、サスペンション歪み検出センサ又はホイールサスペンション振動検出センサ等からの信号に基づき行うこともできる。
【００５０】
また、ケプストラム解析によって、タイヤ回転ピッチ成分を求めたが、ＬＰＣケプストラムの係数によって求めることもできる。
また、スタッドレスタイヤ、サマータイヤ等タイヤの種類に応じてもタイヤ回転ピッチが変化するため、タイヤの種類別に路面状態の判別を行うこともできる。例えば、タイヤの種類別に図６に示すような路面状態判別用マップを設定しておき、車両走行時には予めこのタイヤの種類を示す記号などを入力することで、タイヤの種類に応じた路面状態判別用マップに補正すればよい。
【００５１】
さらに、これと同様にアスファルト等の路面が均一な高μ路を車両が走行中に、状態設定用のスイッチ等を押すことにより、初期状態として入力してその時のタイヤ回転ピッチパラメータＣｐを保存し、この値を基準値として路面状態判別を行うこともできる。
【図面の簡単な説明】
【図１】（ａ）はタイヤが路面からうける力を示す模式図であり、（ｂ）及び（ｃ）は路面状態に対するタイヤ回転ピッチの大きさを示す模式図である。
【図２】車両各センサにおける回路的な構成を示す模式図である。
【図３】路面状態に応じてサンプルとして検出される波形を示した説明図である。
【図４】第１実施形態における路面状態判別処理のフローチャートである。
【図５】図４における処理に対応した波形を示す特性図である。
【図６】サンプル間隔として、時間又は距離を採用した場合における比較図である。
【図７】各種の路面状態に応じたケプストラム解析の結果を示す説明図である。
【図８】路面状態判別を行う基準とする選択マップを示す説明図である。
【図９】第２実施形態における路面状態判別処理のフローチャートである。
【図１０】図９における処理に対応した波形を示す特性図である。
【符号の説明】
１〜４…車輪、５…電子制御装置、６…ステアリングセンサ、
７…ウォーニングランプ、１１〜１４…車輪速センサ、５０…ローター。
【図１】
0003684726.tif 000009
【図２】
0003684726.tif 000010
【図３】
0003684726.tif 000011
【図４】
0003684726.tif 000012
【図５】
0003684726.tif 000013
【図６】
0003684726.tif 000014
【図７】
0003684726.tif 000015
【図８】
0003684726.tif 000016
【図９】
0003684726.tif 000017
【図１０】
0003684726.tif 000018