(57)【要約】【課題】多数の車両に囲まれて併走する場合や、渋滞している中での車線変更をする際等、複雑な状況での対応が可能な車両用推奨操作量生成装置を提供する。【解決手段】周囲車両検出手段６と、自車状態検出手段５と、周囲車両挙動予測手段３ａと、自車に対する運転操作の望ましさを算出する評価関数を構成する評価関数構成手段３ｂと、自車にとって望ましい操作を計算する推奨操作量演算手段３ｃとを備えた車両用推奨操作量生成装置であって、周囲車両挙動予測手段３ａは、自車の予測応答を出力とする自車モデル３ａ−１と、周囲車両の予測応答を出力とする他車モデル３ａ−２と、自車を含む車両の情報から自車モデル３ａ−１および他車モデル３ａ−２の計算に必要な情報を算出する車両情報抽出関数群３ｂ−３とを有し、他車モデル３ａ−２と自車モデル３ａ−１とを、車両情報抽出関数群３ａ−３で結合することによって構成される。【選択図】    図２
2004152125.tif 000002
【特許請求の範囲】
【請求項１】
自車の周囲車両を検出する周囲車両検出手段と、
前記自車の状態を検出する自車状態検出手段と、
前記周囲車両の挙動を予測する周囲車両挙動予測手段と、
前記周囲車両検出手段の出力と前記自車状態検出手段の出力とから前記自車に対する運転操作の望ましさを算出する評価関数を構成する評価関数構成手段と、
前記周囲車両挙動予測手段の出力と前記評価関数構成手段の出力とから前記自車にとって望ましい操作を計算する推奨操作量演算手段、とを備えた車両用推奨操作量生成装置であって、
前記周囲車両挙動予測手段は、
前記自車の予測応答を出力とする自車モデルと、
前記周囲車両の予測応答を出力とする他車モデルと、
前記自車を含む車両の情報から前記自車モデルおよび前記他車モデルの計算に必要な情報を算出する車両情報抽出関数群と、を有し、
前記他車モデルと前記自車モデルとを、前記車両情報抽出関数群で結合することによって構成されることを特徴とする車両用推奨操作量生成装置。
【請求項２】
請求項１記載の車両用推奨操作量生成装置において、
前記自車モデルは、運転者の操作を入力とし、
前記他車モデルは、検出された前記周囲車両の中で特定の位置関係を満たす車両に関する情報を入力とし、前記周囲車両の予測応答を出力とし、
前記車両情報抽出関数群は、前記自車を含む全検出車両の情報から前記他車モデルの入力に該当する情報を算出し、
前記車両情報抽出関数群は、検出車両ごとに割り当てられた前記他車モデルと前記自車モデルを結合することを特徴とする車両用推奨操作量生成装置。
【請求項３】
請求項１記載の車両用推奨操作量生成装置において、
前記評価関数構成手段は、
検出された前記周囲車両の中で特定の位置関係を満たす車両との相互関係を評価する項を含み、
前記評価関数構成手段の車両情報抽出関数群は、全検出車両の情報から、評価対象となる位置関係を満たす車両の情報を算出することを特徴とする車両用推奨操作量生成装置。
【請求項４】
請求項１乃至３のいずれか記載の車両用推奨操作量生成装置において、
前記周囲車両挙動予測手段および前記評価関数構成手段の前記車両情報抽出関数群は、
予測および評価の対象となる特定の位置関係の一つごとに定義される一つ以上の車両情報抽出関数から構成され、
各前記車両情報抽出関数は、
注目車両と他の検出車両の縦方向の相対的な位置関係を判別する前後判別関数と、
他の検出車両の走行車線位置を判別する車線判別関数と、
前記注目車両に最も近い車両の情報を抽出する近傍情報抽出関数と
から構成されることを特徴とする車両用推奨操作量生成装置。
【請求項５】
請求項４記載の車両用推奨操作量生成装置において、
前記車線判別関数は、
任意の前記検出車両と前記注目車両との横方向の位置の偏差に基づいて、前記走行車線位置の判別を行なうことを特徴とする車両用推奨操作量生成装置。
【請求項６】
請求項４記載の車両用推奨操作量生成装置において、
前記車線判別関数は、
任意の前記検出車両と前記注目車両との横方向の位置の偏差、および横方向の相対速度に基づいて、前記走行車線位置の判別を行なうことを特徴とする車両用推奨操作量生成装置。
【請求項７】
請求項４記載の車両用推奨操作量生成装置において、
前記車線判別関数は、
検出された車線ごとに基準となる横位置を設定し、
任意の前記検出車両の横位置と前記車線ごとの基準位置の偏差に基づいて、前記走行車線位置の判別を行なうことを特徴とする車両用推奨操作量生成装置。
【請求項８】
請求項４記載の車両用推奨操作量生成装置において、
前記車線判別関数は、
検出された車線ごとに基準となる横位置を設定し、
任意の前記検出車両の横位置と前記車線ごとの基準位置の偏差、および前記検出車両の横速度に基づいて、前記走行車線位置の判別を行なうことを特徴とする車両用推奨操作量生成装置。
【請求項９】
請求項４乃至８のいずれか記載の車両用推奨操作量生成装置において、
前記車線判別関数および前記前後判別関数は、
前記検出車両の各情報を表す変数に関して一階微分可能な関数で構成することを特徴とする車両用推奨操作量生成装置。
【請求項１０】
請求項１乃至９のいずれか記載の車両用推奨操作量生成装置において、
前記周囲車両挙動予測手段および前記評価関数構成手段の前記車両情報抽出関数群は、
注目車両について予測および評価の対象となる特定の位置関係に他の検出車両が存在する場合には、該当する車両の縦方向の位置情報を出力し、
存在しない場合には、前記注目車両の縦方向の位置にあらかじめ定められた値を加算して出力することを特徴とする車両用推奨操作量生成装置。
【請求項１１】
請求項１０記載の車両用推奨操作量生成装置において、
前記近傍情報抽出関数は、
二変数の最大値関数または最小値関数をすべての検出車両の縦位置情報に関して逐次的に適用することで実現することを特徴とする車両用推奨操作量生成装置。
【請求項１２】
請求項１０記載の車両用推奨操作量生成装置において、
前記近傍情報抽出関数は、
二変数の最大値関数または最小値関数を微分可能な関数で近似し、
近似された関数をすべての検出車両の縦位置情報に関して逐次的に適用することで実現することを特徴とする車両用推奨操作量生成装置。
【請求項１３】
請求項１０記載の車両用推奨操作量生成装置において、
前記近傍情報抽出関数は、
適当な実数ｎをとり、検出車両の縦位置情報のｎ乗和のｎ乗根、または縦位置情報の逆数のｎ乗和のｎ乗根の逆数をとることによって近似的に実現することを特徴とする車両用推奨操作量生成装置。
【請求項１４】
請求項１乃至９のいずれか記載の車両用推奨操作量生成装置において、
前記周囲車両挙動予測手段および前記評価関数構成手段の前記車両情報抽出関数群は、
注目車両について予測および評価の対象となる特定の位置関係に他の検出車両が存在する場合には、該当する車両の縦方向の速度、横方向の位置、横方向の速度のいずれかを出力し、
存在しない場合には、前記注目車両の縦方向の速度、横方向の位置、横方向の速度のいずれかにあらかじめ定められた値を加算して出力することを特徴とする車両用推奨操作量生成装置。
【請求項１５】
請求項１４記載の車両用推奨操作量生成装置において、
前記近傍情報抽出関数は、
二車両の縦方向の位置を比較し、値が大きい方の車両の縦方向の速度、横方向の位置、横方向の速度のいずれかを出力する関数を、すべての検出車両に関して逐次的に適用することで実現することを特徴とする車両用推奨操作量生成装置。
【請求項１６】
請求項１４記載の車両用推奨操作量生成装置において、
前記近傍情報抽出関数は、
二車両の縦方向の位置を比較し、値が大きい方の車両の縦方向の速度、横方向の位置、横方向の速度のいずれかを出力する関数を微分可能な関数を用いて近似的に表現し、
近似された関数をすべての検出車両に関して逐次的に適用することで実現することを特徴とする車両用推奨操作量生成装置。
【請求項１７】
請求項１４記載の車両用推奨操作量生成装置において、
前記近傍情報抽出関数は、
適当な実数ｎをとり、検出車両の抽出対象情報で重みづけした縦位置変数またはその逆数のｎ乗和を縦位置変数またはその逆数のｎ乗和で除算した量を計算することによって近似的に実現することを特徴とする車両用推奨操作量生成装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、車両に搭載された車両用推奨操作量生成装置に関する。
【０００２】
【従来の技術】
【特許文献】特開２００１−５２２９７号。
【０００３】
車線変更も含むような自動車の運転を支援する装置に関する発明として、上記特許文献に記載されたものがある。この発明においては、周囲の情報をセンサにより取得できるという前提のもとに、自車の希望する運転動作（例えば車線変更）を指定しておき、指定された希望運転動作が実施できるか否かを周囲状況をもとに判断し、希望動作が実施可能と判断された後で、運転動作指令信号を生成、出力する、という運転支援の方法等が示されている。
【０００４】
【発明が解決しようとする課題】
上記特許文献に示された方法では、あらかじめ運転動作のモデルが決められており、動作モデルの切り替えを周囲情報をもとにした判断ロジックで実行するという構成になっている。言い換えると、判断してから切り替える、という二段階の処理から構成されていると言える。このため、多数の車両に囲まれて併走する場合や、渋滞している中での車線変更をする際等、複雑な状況での対応をすることが難しかった。
本発明の目的は、多数の車両に囲まれて併走する場合や、渋滞している中での車線変更をする際等、複雑な状況での対応が可能な車両用推奨操作量生成装置を提供することにある。
【０００５】
【課題を解決するための手段】
上記課題を解決するため、本発明は、自車の周囲車両検出手段と、自車状態検出手段と、周囲車両挙動予測手段と、前記周囲車両検出手段と前記自車状態検出手段から自車に対する運転操作の望ましさを算出する評価関数を構成する評価関数構成手段と、前記周囲車両挙動予測手段と前記評価関数構成手段から自車にとって望ましい操作を計算する推奨操作量演算手段とを備え、前記周囲車両挙動予測手段は、自車の予測応答を出力とする自車モデルと、周囲車両の予測応答を出力とする他車モデルと、自車を含む車両の情報から前記自車モデルおよび前記他車モデルの計算に必要な情報を算出する車両情報抽出関数群とからなり、前記他車モデルと前記自車モデルは、前記車両情報抽出関数群で結合することによって構成される。
【０００６】
【発明の効果】
本発明によれば、多数の車両に囲まれて併走する場合や、渋滞している中での車線変更をする際等、複雑な状況での対応を行うことができる。
【０００７】
【発明の実施の形態】
以下、図面を用いて本発明の実施の形態について詳細に説明する。なお、以下で説明する図面で、同一機能を有するものは同一符号を付け、その繰り返しの説明は省略する。
《第一の実施の形態》
本発明の第一の実施の形態を図１から図８までの図面に基づいて説明する。
図１は本発明である車両用推奨操作量生成装置を構成するのに必要な第一の実施の形態の一配置図である。
図１において、前方レーダー１ａは車両前面に取り付けられ、自車両前方に位置する複数の車両の位置を測定する。画像センサ１ｂも車両前面の適当な位置に取り付けられ、前方レーダー１ａの測定情報を補完するとともに、道路上に引かれた車線を検出する。後方レーダー１ｃは車両背面に取り付けられ、自車両後方に位置する複数の車両の位置を測定する。側方センサ１ｄは左右の車両側面に一個ずつ取り付けられ、前方レーダー１ａと後方レーダー１ｃの死角となる車両側方に位置する車両の位置を検出する。なお、側方センサ１ｄとしては、レーダーを用いることもできるが、超音波センサや画像センサを用いることができる。
【０００８】
車速センサ２はロータリーエンコーダーをホイールに取り付けることで実現できる。ホイールの回転速度に応じた周期のパルス列を出力し、車速の計測値を得る。
【０００９】
演算部３はマイクロコンピュータとその周辺部品から構成され、各センサからの信号を内蔵メモリに記録されたプログラムに従って処理し、計算結果を表示装置４に送る。
【００１０】
表示装置４は液晶画面などの表示用ディスプレイおよび該ディスプレイに表示する画像を描画するためのマイクロコンピュータとその周辺部品から構成され、演算部３から送られてきた信号を内蔵メモリに記録された描画プログラムに従って処理し、画像をディスプレイに描画することで運転者に情報を提示する。
【００１１】
演算部３はマイクロコンピュータのソフトウェア形態により、図２に示すブロック３ａ〜３ｃを構成する。
【００１２】
図２において、３ａは周囲車両挙動予測手段、３ｂは評価関数構成手段、３ｃは推奨操作量演算手段、５は車速センサ２と車線検出器としての画像センサ１ｂにより構成される自車状態検出手段、６は周囲センサ１ａ〜１ｄにより構成される周囲車両検出手段、７は周囲車両全情報である。周囲車両挙動予測手段３ａは、自車モデル３ａ−１、他車モデル３ａ−２、車両情報抽出関数群３ａ−３を有する。評価関数構成手段３ｂは、相互関係評価項３ｂ−１、その他の評価項３ｂ−２、車両情報抽出関数群３ｂ−３を有する。
【００１３】
以下、図３に示すような道路状況における動作の例に基づいて、各ブロックの具体的な構成方法を説明する。
図３は片側二車線の直線道路を自車（車両０と定義）と３台の他車両（他車１、他車２、他車３と図に示すように定義）が走行している場面である。ここでは、道路の進行方向に沿ってｘ座標をとり、道路の進行方向と垂直にｙ座標をとる。各車両のｘ座標の値を表す変数をｘｉ、座標の値を表す変数をｙｉと表記する。ここで、ｉは０から３までの値をとり、車両を区別するインデックスである。また、検出された車線の中心に車線の基準位置をとり、左車線および右車線の基準位置をそれぞれｙＬ、ｙＲと表記する。
【００１４】
周囲センサ１ａ〜１ｄからは周囲車両と自車両との相対的な位置および相対速度を得ることができる。従って、座標系の原点を適当に定めれば、自車および各車のｘ座標とｙ座標の値を具体的に確定することができる。また、車速センサ２からは自車両の縦方向の絶対速度を得ることができるので、自車速に各車の相対速度を加算することで、各車の絶対速度を得ることができる。以上より、センサから得られた周囲車両の全情報を次のようなベクトルの形でまとめることができる。
【００１５】
【数１】
2004152125.tif 000003
ただし、ｖｉ＝ｘｉ・は各車の縦方向の絶対速度を表す。なお、「・」は時間微分を表す。また、括弧（  ）付きの数式においては「・」は変数の上にあるものとし、上付きの「・」と、変数の上にある「・」は同意とする。
【００１６】
センサから以上のような情報を受け取ると、周囲車両挙動予測手段３ａは各検出車両ごとに予測モデルを生成する。予測モデルには自車の挙動を予測する自車モデル３ａ−１と他車の挙動を予測する他車モデル３ａ−２の二種類が用意される。
【００１７】
自車モデル３ａ−１は、自車に対する操作の結果としての車両の挙動を予測する式であり、例えば以下のような予測式をたてることができる。
【００１８】
【数２】
2004152125.tif 000004
【００１９】
【数３】
2004152125.tif 000005
ただし、ｕｘは自車両に対する加減速指令値、ｕｙは自車の横方向の走行位置に対する指令値、ｗは横方向の運動の速さを規定する時定数の逆数である。
【００２０】
他車モデル３ａ−２は、周囲の車両の挙動から他車両の挙動を予測する式である。一般に、ある車両の挙動はその周囲を走行する車両すべての影響を直接的に受けるとは考えにくく、ある特定の位置関係にある車両だけから直接的な影響を受けると考えられる。典型的な例としては先行車に対する追従操作を挙げることができ、追従操作においては同一車線上を走行する先行車の挙動が直接影響するものと考えることができる。すなわち、この場合、同一車線上の先行車（以下、特に混同の恐れがない場合には、単に先行車と表記する）、という特定の位置関係にある車両の挙動から、他車両の挙動を予測することができると考える。例えば、他車両が車線変更をしないという仮定のもとでは、次のような追従モデルを他車挙動の予測式として用いることができる。
【００２１】
【数４】
2004152125.tif 000006
ここで、車両ｉは固有の目標車間時間（車間時間とは車間距離を追従側の車両の速度で割った値として定義される）と希望走行車速を持っていると仮定し、それらをそれぞれ、ｈｉ、ｖｉｄと表記している。また、ｋ１ｉ、ｋ２ｉ、ｋ３ｉは車両ｉの追従特性を決めるパラメータである。ｘｐ、ｖｐが先行車の位置と速度である。
【００２２】
なお、先行車が存在しない時には、各車は希望走行車速で走行すると予想される。そのようなモデルとしては、例えば、
【００２３】
【数５】
2004152125.tif 000007
のような式が考えられるが、先行車の有無によって予測モデルを切り替えるのは煩雑である。そこで、先行車が存在しない場合には、車両ｉの遠方に車両ｉと同じ速度で走行する仮想車両が存在するとみなし、さらに（４）式におけるパラメータｋ１ｉ、ｋ２ｉ、ｋ３ｉを車間距離ｘｐ−ｘｉに応じて変化するパラメータであるとして、ｘｐ−ｘｉが大きくなるにつれて、
【００２４】
【数６】
2004152125.tif 000008
となるようにパラメータを設計すれば、（４）式は（５）式を含むような構成となるので、予測モデルを切り替える必要がなくなる。
【００２５】
自車モデル３ａ−１と異なり、他車モデル３ａ−２は単独では予測計算を行なうことができない。なぜなら、（４）式には先行車の位置と速度という抽象的に定義された変数が残っており、これらの変数に実際の車両に対応する変数を割り当てなければ具体性を持った式として完成しないからである。予測式を具体的に構成するためには、どの車両がどの車両の先行車になっているかを特定した上で、自車モデルと各検出車両に割り当てられた他車モデルを結合して、車群全体の予測式としてまとめあげるという処理が必要である。
【００２６】
車両ｉの先行車が具体的にはどの車両であるのかを特定することは、車群の全情報ベクトルＸが与えられればそれほど困難なことではない。例えば、図３の場面において、各モデルの結合関係が以下のようになることは明らかである。
【００２７】
【数７】
2004152125.tif 000009
ここで、記号ｊ→ｋは車両ｋのモデルの先行車を表す変数に車両ｊの変数を割り当てることを意味し、（Ｖ）は先行車が存在しない車両に割り当てる仮想車両を表す。
【００２８】
問題は、予測計算を行なう場合には、（７）式のような位置関係がそのまま保存されるとは限らず、車線変更が起きるような予測が成立した場合には、先行車を切り替える必要が出てくることである。例えば、図３の場面でも、自車０が３秒後に右車線（車両２の後方）に車線変更するという予測が成立した場合には、モデルの結合関係は、最初の３秒間は（７）式のままであるが、３秒後からは以下のような関係に切り替えなければならない。
【００２９】
【数８】
2004152125.tif 000010
この場合、具体的には車両３のモデルの先行車を車両２から車両０へと切り替える必要があることを示している。以上の様子を模式的に示したのが、図４である。
【００３０】
簡単のため、ここでは、車両０以外の車両は車線変更しないという仮定を置いているが、この仮定を緩めて他車両が車線変更することも想定した場合には、さらに多様なモデルの結合関係のパターンが発生する。
【００３１】
このように、（４）式のような追従モデルを利用した予測式は、予測計算の精度を向上させる可能性がある一方で、先行車の割り当てと切り替え処理が必要になるという問題点が発生する。切り替え処理を記述することは不可能ではないが、検出車両と車線数が増加するにつれて、また車両挙動に関する仮定を緩めるにつれて、起こり得る事象のパターンは爆発的に増加していくので、すべての切り替えパターンを網羅することは非常に大きな設計上の負担となる。そこで、このような変数の割り当てと切り替え処理を自動的に行なうような関数を構成して、各モデルを結合することができれば、切り替え処理を記述する必要がなくなるので、車両数や車線数、車線変更に関する仮定とは無関係に予測式を構成することができる。そのような機能を実現するのが、車両情報抽出関数群３ａ−３である。
【００３２】
車両情報抽出関数群３ａ−３とは、検出車両の全情報ベクトルＸの中から、ある車両に関して特定の位置関係にある車両の情報を抽出する関数の集合である。この場合、予測モデルを結合するためには先行車の位置と速度が必要なので、各車両ｉに関して、全情報ベクトルＸから先行車の位置と速度を抽出する関数の集合として車両情報抽出関数群３ａ−３が構成されることになる。以下、具体的な構成方法を説明する。
【００３３】
最初に先行車の位置を抽出する関数の構成方法を説明する。車両ｉの先行車の位置を抽出する関数は、以下のような値を出力する関数として定義される。
【００３４】
【数９】
2004152125.tif 000011
ここで、Ｒｄｕｍは先行車が存在しないときに配置される仮想車両である車両ｉからの距離を表す適当な定数である。以上のような関数を実現するためには、まず、先行車を数学的に定義する必要がある。車両ｊが車両ｉの先行車であるための条件として、以下の三条件を考えることができる。
条件Ａ：車両ｊは車両ｉよりも前方を走行している。
条件Ｂ：車両ｊは車両ｉと同一車線上を走行している。
条件Ｃ：車両ｊは条件Ａ、Ｂを満たす車両の中で最も車両ｉに近い。
このうち、条件Ｂの同一車線の判定は必ずしも数学的な定義ではないので、ここでは次のような条件をもって同一車線上の車両か否かを判定する。
条件Ｂ’：車両ｊと車両ｉの横方向位置の偏差が車線幅に比べて十分小さい。
【００３５】
条件Ａを数学的に表現すると、
【００３６】
【数１０】
2004152125.tif 000012
となる。
【００３７】
条件Ｂ’の数学的な表現は、同一車線上の車両とみなすか否かのしきい値をｄとして、
【００３８】
【数１１】
2004152125.tif 000013
となる。
【００３９】
条件Ｃの最も車両ｉに近いという条件を満たす車両の位置を返す関数が近傍情報抽出関数である。この場合、最も車両ｉに近いという条件は、条件Ａ、Ｂを満たすｘｊの中で最も値が小さい、という条件と等価になるので、近傍情報抽出関数は入力引数の中から最小値をとる引数を選び出す最小値関数で表現することができる。以上をまとめると、先行車の位置を抽出する関数は、次のように書ける。
【００４０】
【数１２】
2004152125.tif 000014
ただし、ここでは、Ａ，Ｂは、それぞれ条件Ａ、Ｂ’を満たすインデックスの集合であり、δｉｊは
【００４１】
【数１３】
2004152125.tif 000015
である。ｍｉｎは最小値関数である。
【００４２】
（１２）式は、（９）式の条件を満たす関数になっていることは事実である。しかし、推奨操作量を計算することを考えた場合、（１２）式は必ずしも都合のよい関数ではない。なぜなら、（１２）式は微分可能な関数とはなっていないので、推奨操作量を効率良く計算する手段を与えてくれる最適制御理論の枠組みで扱えないからである。
【００４３】
最適制御理論に基づく効率的な計算方法を利用するためには、（１２）式を微分可能な関数で再構成しなければならない。
【００４４】
簡単のため、ここではすべてのインデックスｊが条件Ａ、Ｂを満たしていると仮定する。このとき、（１２）式はすべてのインデックスに関する最小値関数となる。最小値関数自体が微分不可能な関数なので、まず最小値関数を微分可能な関数で近似する方法を考える。
【００４５】
一般に、Ｎ個の正の値をとる変数ｚｊ、ｊ＝１，２，…，Ｎがあったとして、ｚｊの最大値に対応するインデックスをＭ、最小値に対応するインデックスをｍであらわすとする。すなわち、
【００４６】
【数１４】
2004152125.tif 000016
【００４７】
【数１５】
2004152125.tif 000017
である。このとき、適当な正の定数ｎをとり、
【００４８】
【数１６】
2004152125.tif 000018
という式を計算すると、（１６）式はｎの値を大きくするにつれてｚｊの最大値ｚＭへと収束していく。なぜなら、（１６）式は以下のように変形できるからである。
【００４９】
【数１７】
2004152125.tif 000019
（１７）式の右辺の根号内部は、ｚＭが最大値であるという仮定より、ｚＭの項だけが１として残り、その他の項はｎを大きくするにつれて０に収束していくからである。
【００５０】
以上より、
【００５１】
【数１８】
2004152125.tif 000020
という近似式を構成することができた。（１８）式の左辺は、ｚｊに関して微分可能な関数となっているので、推奨操作量の計算を行なうのに好都合な表現となっている。
【００５２】
最大値関数を構成できれば、それから最小値関数を構成することができる。（１６）式のかわりに
【００５３】
【数１９】
2004152125.tif 000021
を計算すると、（１８）式より、
【００５４】
【数２０】
2004152125.tif 000022
が成り立つので、最小値関数は、
【００５５】
【数２１】
2004152125.tif 000023
で近似できる。従って、（１２）式に表れる最小値関数は（２１）式で置き換えれば、微分可能な関数として再構成することができる。
【００５６】
なお、（２１）式の最小値関数を用いるときには、変数ｚｊはすべて正であるという条件がついていた。これは（１２）式に適用する場合には、車両の位置ｘｊがすべて正であるという条件となる。座標系の原点は任意にとることができるので、この条件を満たすことは容易である。すなわち、後方物体を検出するセンサの検出レンジよりもさらに後方に原点を定めれば、ｘｊがすべて正という条件は自動的に満たされる。
【００５７】
次に、条件Ａ、Ｂ’を（２１）式の中に反映させることを考える。そのためには、条件Ａ、Ｂ’を満たさないインデックスに対応する項を消去、すなわち、０になるような仕組みを導入すればいいことがわかる。
【００５８】
そこで、（１０）式で表される条件Ａを満たすときに値１をとり、満たさないときに値０をとる関数である前後判別関数ｃｘ（ｘｉ，ｘｊ）と、（１１）式で表される条件Ｂ’を満たすときに値１をとり、満たさないときに値０をとる関数である車線判別関数ｃｙ（ｙｉ，ｙｊ）を導入する。
【００５９】
そのような性質を持つ二つの判別関数をｚｊ−ｎの項に作用させれば、条件Ａ，Ｂ’を満たさない項は０となるので、最小値演算の対象から自動的に除外されることになる。
【００６０】
以上をまとめると、（１２）式は以下の関数で近似できることになる。
【００６１】
【数２２】
2004152125.tif 000024
関数ｃｘ（ｘｉ，ｘｊ）、ｃｙ（ｙｉ，ｙｊ）を微分可能な関数で構成すれば、（２２）式全体も微分可能な関数となる。そこで、前後判別関数を次のような微分可能な関数で構成する。
【００６２】
【数２３】
2004152125.tif 000025
と定義する。ここで、ｄｘは関数の連続化のために導入した中間状態の幅に相当する適当な正の定数である。（２３）式をグラフとしてプロットしたものを図５に示す。また、車線判別関数ｃｙ（ｙｉ，ｙｊ）を次のような微分可能な関数で構成する。
【００６３】
【数２４】
2004152125.tif 000026
と定義する。ここで、ｄ１は車両ｉの運転者が車両ｊを同一車線上の車両と認識する横方向の位置偏差の限界値であり、ｄ２は車両ｉの運転者が車両ｊを隣接車線上の車両と認識する横方向の位置偏差の限界値であり、関数の連続化のためにｄ２＞ｄ１となるように設定され、幅ｄ２−ｄ１の中間領域が形成される。（２４）式をグラフとしてプロットしたものを図６に示す。
【００６４】
以上、先行車の縦位置を検出車両の全情報ベクトルＸから抽出する微分可能な関数を構成した。関数（２２）式を用いることで、任意の車両ｉの先行車の縦方向の位置ｘｐを微分可能な関数を用いて特定することができる。
【００６５】
（４）式の予測モデルでは、先行車の位置ｘｐだけでなく、先行車の速度ｖｐの情報も必要になる。そこで、次に先行車の速度を抽出する関数の構成方法について説明する。
【００６６】
車両ｉの先行車の速度を抽出する関数は、以下のような値を出力する関数として定義される。
【００６７】
【数２５】
2004152125.tif 000027
正の値をとる変数ｚｊ、ｊ＝１，２，…，Ｎに対応して、同じく正の値をとる変数ｖｊが定義されているものとする。このとき、ｖｊ、ｊ＝１，２，…，Ｎの中から、ｖＭとｖｍを抽出する関数を構成できればよいことになる。そこで、次のような式を構成する。
【００６８】
【数２６】
2004152125.tif 000028
（２６）式は次のように変形できる。
【００６９】
【数２７】
2004152125.tif 000029
ｚＭは（１６）式で近似できるので、（２７）式の両辺をｚＭで割って、さらに両辺をｎ乗すれば、次式が得られる。
【００７０】
【数２８】
2004152125.tif 000030
同様にして、ｖｍについても次式のような近似式を導出することができる。
【００７１】
【数２９】
2004152125.tif 000031
以上より、先行車の条件Ａ，Ｂ’を考慮すれば、先行車の速度を抽出する微分可能な関数として、次式のような関数を構成することができる。
【００７２】
【数３０】
2004152125.tif 000032
以上、先行車の速度を検出車両の全情報ベクトルＸから抽出する微分可能な関数を構成した。関数（３０）式を用いることで、任意の車両ｉの先行車の縦方向の速度ｖｐを微分可能な関数を用いて特定することができる。
【００７３】
（２２）式、（３０）式を用いることで、図３の状況から、他車の車線変更がないという仮定のもとで、周囲車両の挙動を予測する単一の予測式系を以下のように構成することができる。
【００７４】
【数３１】
2004152125.tif 000033
（３１）式は、自車の車線変更に伴う式の切り替えが不要な上、式全体が微分可能な関数で構成されており、最適制御理論を用いた計算に適した形式となっている。
【００７５】
以上が周囲車両挙動予測手段３ａ（図２）で行なわれる処理の内容である。
【００７６】
次に評価関数構成手段３ｂについて説明する。
評価関数構成は、自車に対する操作の結果として予想される周囲車両の未来の挙動を評価し、自車に対する操作がどれくらい望ましいものであるかを評価する関数であり、一般形として次式のように表現することができる。
【００７７】
【数３２】
2004152125.tif 000034
ここで、ｔは現在の時刻を表す時間変数、Ｔは評価時間の長さを表す正の定数、τは予測計算における時間を表示するための時間変数である。関数Ｌ（Ｘ，ｕｘ，ｕｙ，τ）は望ましい操作を定義する一つ以上の評価項から構成される。評価項として表現される運転操作に対する要請としてはさまざまなものが考えられるが、大きく分類すると、自車の状態だけを評価する項と、他車との位置関係を評価する項の二つから構成されると考えられる。自車の状態だけを評価する項としては、例えば自車に対する加減速を小さくするという要請を表現した
【００７８】
【数３３】
2004152125.tif 000035
といった項や、あるいはなるべく車線の中央を走行するという要請を表現した
【００７９】
【数３４】
2004152125.tif 000036
という項や、あるいはなるべく希望車速に近い速度で走行するという要請を表現した
【００８０】
【数３５】
2004152125.tif 000037
といった評価項が考えられる。
【００８１】
一方、他車との位置関係を評価する項としては、基本的な要請として他車両に近づきすぎないという要請を考慮する必要がある。自動車の運転の場合、車線によって走行方向が規定されている場の中での運動となるので、そのような制約を考慮した上で評価項を構成する必要がある。例えば、縦方向の車間距離が同じであっても、同一車線上の車両と隣接車線上の車両とを同等に評価することは不適切である。また、たとえ同一車線上の車両であったとしても、直前の先行車ともう一台別の車両を挟んだ二台目の先行車とでは、やはり同等に評価することはできない。他にも、前方車両と後方車両に対する評価を異なるものにした方がよい場合もある。このように、評価関数の構成においても、他車両との位置関係を特定することは重要な処理となる。また、予測を立てる場合には、評価対象となる位置関係を占める車両が予測区間の中で切り替わることもあり得るので、効率的に評価関数を構成するためには予測式で用いたような車両情報抽出関数群３ａ−３が必要となってくる。
【００８２】
図７に示すように、図３のような片側二車線道路において考慮する必要がある位置関係というのは、左右両車線の先行車と後続車の全部で四方向であると考えられる。この四方向に存在する他車両と適切な距離が保たれていれば、リスクの大きな追従操作や車線変更操作を避けることができる。そこで、他車との位置関係を評価する評価項として、次式を導入する。
【００８３】
【数３６】
2004152125.tif 000038
ここで、ｘｐＬ、ｘｐＲ、ｘｆＬ、ｘｆＲはそれぞれ、左車線の先行車、右車線の先行車、左車線の後続車、右車線の後続車の縦方向の位置を表し、関数ｃＬ（ｙ０）は自車が左車線を走行している場合に値１をとり、右車線を走行している場合には０になる関数で、この場合、
【００８４】
【数３７】
2004152125.tif 000039
のような関数を利用することができる。関数ｃＲ（ｙ０）はその逆で、左車線で１、右車線で０をとる関数で、
【００８５】
【数３８】
2004152125.tif 000040
のような関数を利用することができる。
【００８６】
ｘｐＬ、ｘｐＲに具体的な検出車両の変数を割り当てるには、予測式の構成の部分で説明した先行車抽出関数をそのまま利用することができる。ただし、この場合、左車線および右車線上の先行車を抽出することになるので、車線判別関数については、比較の対象となるのは自車との横方向位置の偏差ではなく、ｙＬまたはｙＲとの偏差を比較することになる。以上より、ｘｐＬ、ｘｐＲを抽出する関数の具体形として、
【００８７】
【数３９】
2004152125.tif 000041
【００８８】
【数４０】
2004152125.tif 000042
を導くことができる。
【００８９】
ｘｆＬ、ｘｆＲについても考え方はまったく同じであるが、この場合、自車の前方車両ではなく後方車両の中から該当車両を探すことになるため、前後判別関数に代入する変数の順序を逆にする必要があることと、最小値関数ではなく最大値関数を利用する部分を修正する必要がある。式の具体形を書き下すと、
【００９０】
【数４１】
2004152125.tif 000043
【００９１】
【数４２】
2004152125.tif 000044
と書くことができる。
【００９２】
以上に導入した各評価項の線形加重和をとることで評価関数全体を構成する。すなわち、
【００９３】
【数４３】
2004152125.tif 000045
を評価関数（の被積分関数）とする。ここで、ｗＸ、ｗＹ、ｗＶ、ｗＲは各評価項に対する重みパラメータである。
【００９４】
なお、ここでは用いなかったが、他車との位置関係を評価するとき、他車の位置だけでなく、速度も評価式の中に入れることも可能である。例えば、（３６）式の評価項では車間距離の逆数の二乗が評価式として用いられているが、車間距離ではなく車間時間の逆数の二乗を評価式として用いる方法も考えられる。その場合、例えば左車線の先行車に対する評価項は
【００９５】
【数４４】
2004152125.tif 000046
となり、左車線の先行車の速度ｖｐＬを求めるのに予測式の部分で説明した速度抽出関数を変形して利用することができる。ｖｐＬ、ｖｐＲ、ｖｆＬ、ｖｆＲをそれぞれ、左車線の先行車、右車線の先行車、左車線の後続車、右車線の後続車の速度とするなら、速度抽出関数の具体形は次のようになる。
【００９６】
【数４５】
2004152125.tif 000047
【００９７】
【数４６】
2004152125.tif 000048
【００９８】
【数４７】
2004152125.tif 000049
【００９９】
【数４８】
2004152125.tif 000050
以上のように予測式と評価関数を構成すると、推奨操作量を計算する問題は最適制御問題に帰着することができる。すなわち、（３１）式のような車群挙動が予想されている中で、（３２）、（４３）式から定義される評価関数を最小にするような操作量の系列ｕｘ（τ）、ｕｙ（τ）、ｔ≦τ≦ｔ＋Ｔを求めるという最適制御問題となる。最適制御問題の系統的な解法はよく研究されており（例えば、文献１：加藤寛一郎著  工学的最適制御  非線形へのアプローチ、文献２：Ｔ．  Ｏｈｔｓｕｋａ，  Ｃｏｎｔｉｎｕａｔｉｏｎ／ＧＭＲＥＳ  ｍｅｔｈｏｄ  ｆｏｒ  ｆａｓｔ  ａｌｇｏｒｉｔｈｍ  ｏｆ  ｎｏｎｌｉｎｅａｒ  ｒｅｃｅｄｉｎｇｈｏｒｉｚｏｎ  ｃｏｎｔｒｏｌ  Ｐｒｏｃ．  ３９ｔｈ  ＩＥＥＥ  Ｃｏｎｆｅｒｅｎｃｅ  ｏｎ  Ｄｅｃｉｓｉｏｎ  ａｎｄ  Ｃｏｎｔｒｏｌ，  ｐｐ．７６６−７７１，  ２０００．  など）、効率的に解を計算する手法が数多く提案されている。推奨操作量演算手段はそのような計算を実行することによって、実際に推奨操作量を算出する処理を行なう。推奨操作量を計算して得られる結果の例を図８に示す。
【０１００】
図８の例では、加速しながら他車３の前方の右車線に車線変更した後、緩く減速しながら他車１の前方の左車線に車線変更する操作が生成されている。このように車線変更が二度起きるような例では、本来は予測式と評価関数の切り替え処理を最適化計算の最中に行なわなければいけないが、車両情報抽出関数群３ａ−３を用いることによって切り替え処理が不要になり、どのような動作に対しても単一の予測式と評価関数で最適制御問題を定義、計算できるようになる。生成された操作量系列は表示装置４を通して運転者に知らされる。
【０１０１】
まとめとして、全体としての処理の流れを図９のフローチャートに沿って説明する。
ステップ１では、周囲センサ１ａ〜１ｄおよび自車状態検出手段（車速センサ２、車線検出器としての画像センサ１ｂ）のセンサ信号を読み込む。この時点で各車の速度および自車と周囲車両との相対的な位置を算出することができる。
【０１０２】
ステップ２では、縦方向に適当な原点をとって座標系を導入し、自車と周囲車両の縦方向の位置を表す変数ｘｉ、ｉ＝０，…，Ｎ（Ｎは検出した他車両の台数）の値を具体的に定める。このとき、（２２）式をはじめとする車両情報抽出関数群３ａ−３は、ｘｉ＞０でないと正しく機能しないので、すべての座標値が正となるように座標系の原点を定める。具体的には、例えば最も後方に検出した車両のさらに１０ｍ後方に原点をとる、といったようなルールを定めればよい。
【０１０３】
ステップ３では、横方向に適当な原点をとって座標系を導入し、自車と周囲車両の横方向の位置を表す変数ｙｉ、ｉ＝０，…，Ｎの値を具体的に定める。座標系の原点はどこにとってもよいが、例えば左車線の左側の白線上などにとることができる。このとき、各車両の横方向の座標だけでなく、左車線と右車線の基準位置ｙＬ、ｙＲの値も確定する。以上の処理により、周囲車両の全情報をまとめた（１）式で表されるベクトルＸが構成される。
【０１０４】
ステップ４では、検出した車両ごとに、必要な車両情報抽出関数群３ａ−３が構成される。この第一の実施の形態の場合、他車については周囲車両群挙動予測式を構成するために、（２２）、（３０）式の先行車抽出関数が各車ごとに割り当てられる。自車については評価関数を構成するために、（３９）〜（４３）式の車両情報抽出関数群３ａ−３が割り当てられる。
【０１０５】
ステップ５では、自車モデルが生成される。具体的には（２）、（３）式がメモリー上に読み出される。
【０１０６】
ステップ６では、検出車両１台ごとに、他車モデルが生成される。具体的には、（４）式のモデルが適当に割り当てられたパラメータとともにメモリー上に生成される。
【０１０７】
ステップ７では、ステップ５とステップ６で生成された自車モデルと他車モデルをステップ４で生成した車両情報抽出関数群３ａ−３で結合する。具体的には、（３１）式で表されるような微分方程式系としてモデルをまとめあげることになる。
【０１０８】
ステップ８では、評価関数を設定する。具体的にはあらかじめ設定されている（４３）式のような評価関数を読み出して、ステップ４で構成した車両情報抽出関数群３ａ−３を適用することで、評価関数を現在の周囲状況に適合する形に再構成する。
【０１０９】
ステップ９では、ステップ７で設定された車両情報抽出関数群３ａ−３と、ステップ８で設定された評価関数を、推奨操作量演算手段に入力し、推奨操作量を生成する。
【０１１０】
ステップ１０では、生成された推奨操作量を表示装置４に転送して処理を終える。
【０１１１】
上記が、第一の実施の形態における車両用推奨操作量生成装置の処理内容である。
以上が本発明の第一の実施の形態である。
【０１１２】
以上説明したように、本実施の形態では、自車の周囲車両を検出する周囲車両検出手段６と、自車の状態を検出する自車状態検出手段５と、周囲車両の挙動を予測する周囲車両挙動予測手段３ａと、周囲車両検出手段６の出力と自車状態検出手段５の出力とから自車に対する運転操作の望ましさを算出する評価関数を構成する評価関数構成手段３ｂと、周囲車両挙動予測手段３ａの出力と評価関数構成手段３ｂの出力とから自車にとって望ましい操作を計算する推奨操作量演算手段３ｃとを備えた車両用推奨操作量生成装置であって、周囲車両挙動予測手段３ａは、自車の予測応答を出力とする自車モデル３ａ−１と、周囲車両の予測応答を出力とする他車モデル３ａ−２と、自車を含む車両の情報から自車モデル３ａ−１および他車モデル３ａ−２の計算に必要な情報を算出する車両情報抽出関数群３ａ−３とを有し、他車モデル３ａ−２と自車モデル３ａ−１とを、車両情報抽出関数群３ａ−３で結合することによって構成されることを特徴とする。
このような構成により、予測の結果として起こる車線変更などによって車群の配置状態に変化が生じても、車両情報抽出関数群３ａ−３によって各車両モデルの間の依存関係が自動的に切り替わるので、事前に起こり得る事象をすべて列挙して予測式の切り替え規則を作る必要がなくなり、単一の予測式系で複雑な車両挙動の予測計算ができるようになる。この結果、多数の車両に囲まれて併走する場合や、渋滞している中での車線変更をする際等、複雑な状況での対応を行うことができる。すなわち、周囲車両挙動予測手段３ａの設計の手間を大幅に軽減することができる。
【０１１３】
また、自車モデル３ａ−１は、運転者の操作を入力とし、他車モデル３ａ−２は、検出された周囲車両の中で特定の位置関係を満たす車両に関する情報を入力とし、周囲車両の予測応答を出力とし、車両情報抽出関数群３ｂ−３は、自車を含む全検出車両の情報から他車モデル３ａ−２の入力に該当する情報を算出し、車両情報抽出関数群３ａ−３は、検出車両ごとに割り当てられた他車モデル３ａ−２と自車モデル３ａ−１を結合することを特徴とする。
このような構成により、予測の結果として起こる車線変更などによって車群の配置状態に変化が生じても、車両情報抽出関数群３ａ−３によって各車両モデルの間の依存関係が自動的に切り替わるので、事前に起こり得る事象をすべて列挙して予測式の切り替え規則を作る必要がなくなり、単一の予測式系で複雑な車両挙動の予測計算ができるようになる。
【０１１４】
また、評価関数構成手段３ｂは、検出された周囲車両の中で特定の位置関係を満たす車両との相互関係を評価する項（相互関係評価項３ｂ−１）を含み、評価関数構成手段３ｂの車両情報抽出関数群３ｂ−３は、全検出車両の情報から、評価対象となる位置関係を満たす車両の情報を算出することを特徴とする。
このような構成により、特定の位置関係にある車両（例えば先行車など）だけを評価する評価関数を構成する際に、評価の対象となる車両が予測の途中で切り替わっても、車両情報抽出関数群３ｂ−３によって車両の切り替えが自動的に行なわれるので、事前に起こり得る事象をすべて列挙して評価関数の切り替え規則を作る必要がなくなり、単一の評価関数でそのような評価を表現することが可能になる。すなわち、評価関数構成手段３ｂの設計の手間を大幅に軽減することができる。
【０１１５】
また、周囲車両挙動予測手段３ａおよび評価関数構成手段３ｂの車両情報抽出関数群３ａ−３、３ｂ−３は、予測および評価の対象となる特定の位置関係の一つごとに定義される一つ以上の車両情報抽出関数から構成され、各車両情報抽出関数は、注目車両と他の検出車両の縦方向の相対的な位置関係を判別する前後判別関数と、他の検出車両の走行車線位置を判別する車線判別関数と、  注目車両に最も近い車両の情報を抽出する近傍情報抽出関数とから構成されることを特徴とする。
このような構成により、前後関係の判別、車線の判別、最近傍車両情報の抽出、という三つの要素から関数を構成することにより、予測と評価の対象となる特定の位置関係として各車線ごとの先行車および後続車を表現することができる。人間の運転行動における他車との相互作用も他車との前後関係と走行車線位置の関係によって影響されていると考えることができるので、実際の運転者の感覚に近い予測式と評価関数を構成することができる。
【０１１６】
また、車線判別関数は、任意の検出車両と注目車両との横方向の位置の偏差に基づいて、走行車線位置の判別を行なうことを特徴とする。
このような構成により、横方向の相対的な位置の偏差に基づいて他車両の車線位置を判別しているので、予測と評価の対象となる特定の位置関係として、例えば右方向の隣接車線上を走行する車両、といったような相対的な車線位置に基づく位置関係を規定することができる。
【０１１７】
また、車線判別関数は、検出された車線ごとに基準となる横位置を設定し、任意の検出車両の横位置と車線ごとの基準位置の偏差に基づいて、走行車線位置の判別を行なうことを特徴とする。
このような構成により、車線を基準とした絶対的な位置に基づいて他車両の車線位置を判別しているので、予測と評価の対象となる特定の位置関係として、例えば左から２番目の車線上を走行する車両、といったような絶対的な車線位置に基づく位置関係を規定することができる。
【０１１８】
また、車線判別関数および前後判別関数は、検出車両の各情報を表す変数に関して一階微分可能な関数で構成することを特徴とする。
このような構成により、車線判別関数を微分可能な関数で構成することにより、推奨操作量を計算する場合に偏微分計算に基づく最適制御理論を活用することができるので、効率的な計算を行なうことができる。
【０１１９】
また、周囲車両挙動予測手段３ａおよび評価関数構成手段３ｂの車両情報抽出関数群３ａ−３、３ｂ−３は、注目車両について予測および評価の対象となる特定の位置関係に他の検出車両が存在する場合には、該当する車両の縦方向の位置情報を出力し、存在しない場合には、注目車両の縦方向の位置にあらかじめ定められた値を加算して出力することを特徴とする。
このような構成により、車両の情報として縦方向の位置を出力することで、縦方向の車間距離を入力とする他車モデル３ａ−２や縦方向の車間距離を評価する評価関数を用いることができるようになり、追従動作の予測や追突回避といった基本的な予測と評価を表現することができる。また、該当車両が存在しない場合に自車の遠方に仮想的な車両が存在するものとして、その仮想的な車両の縦位置を出力する構成とし、遠方車両の影響が小さくなるような予測式と評価関数を設定することで、該当車両の有無によって予測式や評価関数を切り替える必要がなくなる。
【０１２０】
また、近傍情報抽出関数は、適当な実数ｎをとり、検出車両の縦位置情報のｎ乗和のｎ乗根、または縦位置情報の逆数のｎ乗和のｎ乗根の逆数をとることによって近似的に実現することを特徴とする。
このような構成により、本来は不連続関数である近傍情報抽出関数を微分可能な関数で近似的に実現することによって、車両の縦位置を抽出する関数全体を微分可能な関数で構成することが可能になり、推奨操作量を計算する場合に偏微分計算に基づく最適制御理論を活用することができるので、効率的な計算を行なうことができる。また、近傍情報抽出関数は、二変数の最大値関数または最小値関数をすべての検出車両の縦位置情報に関して逐次的に適用することで実現する場合と異なり、逐次的な構成方法によらずに関数が構成できるので、関数の表現形式がより簡潔なものとなり、最適化計算のアルゴリズムを容易に設計することができる。
【０１２１】
また、近傍情報抽出関数は、適当な実数ｎをとり、検出車両の抽出対象情報で重みづけした縦位置変数またはその逆数のｎ乗和を縦位置変数またはその逆数のｎ乗和で除算した量を計算することによって近似的に実現することを特徴とする。
このような構成により、本来は不連続関数である近傍情報抽出関数を微分可能な関数で近似的に実現することによって、車両情報抽出関数全体を微分可能な関数で構成することが可能になり、推奨操作量を計算する場合に偏微分計算に基づく最適制御理論を活用することができるので、効率的な計算を行なうことができる。また、車両情報抽出関数群３ａ−３が、注目車両について予測および評価の対象となる特定の位置関係に他の検出車両が存在する場合には、該当する車両の縦方向の速度、横方向の位置、横方向の速度のいずれかを出力し、存在しない場合には、注目車両の縦方向の速度、横方向の位置、横方向の速度のいずれかにあらかじめ定められた値を加算して出力する場合と異なり、逐次的な構成方法によらずに関数が構成できるので、関数の表現形式がより簡潔なものとなり、最適化計算のアルゴリズムを容易に設計することができる。
【０１２２】
《第二の実施の形態》
次に、本発明の第二の実施の形態を説明する。
本発明の第二の実施の形態における装置の基本的な構成は、第一の実施の形態と同様であり、図１に示した配置図と図２に示した構成図は同一である。
【０１２３】
第一の実施の形態では、車両情報抽出関数群３ａ−３を構成するのに必要な最大値関数あるいは最小値関数を（１８）式あるいは（２１）式を用いて構成していた。（１８）式あるいは（２１）式の方法は、べき乗関数による収束を利用して最大値あるいは最小値関数を近似する方法であった。この方法は比較的簡潔な数式表現で実現できることができる一方、べき乗関数の指数ｎの値に注意が必要である。ｎが小さすぎると収束が十分でなくなり、抽出される値と実際の値との間に生じる誤差が大きくなる恐れがある。逆にｎが大きすぎると、今度は数値計算自体が不安定になる恐れがある。適当な指数ｎが見つかれば実際の問題にはならないが、本実施の形態では、べき乗関数の収束性を利用しないで、最大値または最小値関数を微分可能な関数で近似する方法を示す。
【０１２４】
ここで、再びＮ個の変数ｚｊ、ｊ＝１，２，…，Ｎを考え、ｚｊの最大値ｚＭを抽出する方法を考える。もし、Ｎ＝２、すなわち、二つの変数のうちで値が大きい方の変数を抽出する問題を考えると、例えば次のような最大値関数の近似関数を構成することができる。
【０１２５】
【数４９】
2004152125.tif 000051
ここで、ｋは適当な正のパラメータである。（４９）式は、ｚ１とｚ２との値の差が十分に大きければ、最大値関数とほぼ同じ値を出力する。
【０１２６】
二変数の最大値関数が構成できれば、三変数の最大値関数は次のようにすれば構成できる。
【０１２７】
【数５０】
2004152125.tif 000052
ｚ１とｚ２の大きい方の値とｚ３とを比較し、大きい方の値が出力されれば、出力された値は三つの変数の中で最も大きい値になっているはずであり、（５０）はそれを利用したものである。変数がＮ個ある場合には、（５０）式と同様の比較を繰り返していけばよい。すなわち、
【０１２８】
【数５１】
2004152125.tif 000053
でＮ変数の最大値関数を構成できる。（５１）式の表現は再帰的な表現であり、陽な表現形式とはなっていないが、微分可能な関数ｆｍａｘの合成関数として定義されているので、全体も微分可能な関数になっている。最小値関数については、関数ｆｍａｘの替わりに以下の関数を用いればよい。
【０１２９】
【数５２】
2004152125.tif 000054
Ｎ変数関数については、（５１）式と同様、以下のように構成できる。
【０１３０】
【数５３】
2004152125.tif 000055
（２２）式に対応する先行車抽出関数は、（５３）式を利用すれば、
【０１３１】
【数５４】
2004152125.tif 000056
と構成することができる。ただし、（５３）式が正しく機能するためには、座標系の原点を、ｘｉ＋Ｒｄｕｍ＜０となるようにとらなければならない。すなわち、座標系の原点を周囲センサの前方検出レンジよりもさらに前方に設定する必要がある。
【０１３２】
次に、（３０）式に対応する先行車の速度を抽出する関数を以上のような逐次比較の方法を用いて構成する。今、（ｚｊ，ｖｊ）、ｊ＝１，…，Ｎというペアを考え、ｚｊが最大になるインデックスＭに対応するｖＭを抽出するという問題を考える。ここでも、まず、Ｎ＝２の問題を考えると、次のような抽出関数を構成することができる。
【０１３３】
【数５５】
2004152125.tif 000057
二変数間の抽出操作を繰り返せば、変数の数がいくつになってもｖＭを抽出することができる。すなわち、
【０１３４】
【数５６】
2004152125.tif 000058
で、Ｎ変数の抽出関数になる。式は複雑であるが、微分可能な関数の合成関数から構成されているので、全体の式も微分可能になっている。
【０１３５】
最小値関数については、
【０１３６】
【数５７】
2004152125.tif 000059
を定義して、
【０１３７】
【数５８】
2004152125.tif 000060
とすれば構成できる。
【０１３８】
（３０）式に相当する先行車速度抽出関数は、
【０１３９】
【数５９】
2004152125.tif 000061
で構成される。
【０１４０】
第一の実施の形態における他の抽出関数（３９）〜（４２）式、および（４４）〜（４７）式についても、（５４）式、（５９）式に準じて構成することができる。実施の形態の他の要素については、第一の実施の形態の場合とまったく同じである。
以上が本発明の第二の実施の形態である。
【０１４１】
上記のように、本実施の形態では、近傍情報抽出関数は、二変数の最大値関数または最小値関数をすべての検出車両の縦位置情報に関して逐次的に適用することで実現することを特徴とする。
このような構成により、二変数の逐次比較によって最大値関数または最小値関数が構成されるので、前後判定関数と車線判別関数と組み合わせることにより、注目車両に最も近い先行車または後続車の縦方向の位置を系統的に抽出することができる。
【０１４２】
また、近傍情報抽出関数は、二変数の最大値関数または最小値関数を微分可能な関数で近似し、近似された関数をすべての検出車両の縦位置情報に関して逐次的に適用することで実現することを特徴とする。
このような構成により、本来は不連続関数である近傍情報抽出関数を微分可能な２変数比較関数を逐次的に適用することで近似的に実現することによって、車両の縦位置を抽出する関数全体を微分可能な関数で構成することが可能になり、推奨操作量を計算する場合に偏微分計算に基づく最適制御理論を活用することができるので、効率的な計算を行なうことができる。
【０１４３】
《第三の実施の形態》
本発明の第三の実施の形態を図１０〜図１２の図面に基づいて説明する。
本発明の第二の実施の形態における装置の基本的な構成は、第一の実施の形態と同様であり、図１に示した配置図と図２に示した構成図は同一である。
【０１４４】
第一の実施の形態では、自車だけが車線変更し、他車は車線変更しないという仮定を置いていた。第三の実施の形態では、逆に自車は車線変更せず、他車の車線変更が予想される場面に基づいて説明を行なう。場面としては第一の実施の形態と同様、図３のような場面を考える。ただし、本実施の形態においては、他車３が前方を走行している他車２を追い越すために、左車線に車線変更する意志をもっているものとする。
【０１４５】
本実施の形態では、周囲センサは周囲車両の横方向の位置だけでなく横方向の速度も算出する構成を仮定する。各車の横方向の速度をｗｉであらわすと、周囲車両の全情報ベクトルＸは、
【０１４６】
【数６０】
2004152125.tif 000062
となる。
【０１４７】
他車１から３は、縦方向の追従モデルは（４）式が割り当てられているものとする。一方横方向の運動については、次のようなダイナミクスを割り当てる。
【０１４８】
【数６１】
2004152125.tif 000063
ここで、ｕｉｙは横方向の速度に対する指令値であり、後述する車線変更の判断ロジックに基づいて決定される。
【０１４９】
各車の挙動を予測する予測モデルの構成において、車線変更の動作を予測するために、各車にはダイナミクス以外に車線変更するか否かを判断するロジックも割り当てられる。車線変更のロジックは、車線変更の意志を決定するロジックと車線変更の可能性を判断するロジックに分割される。
【０１５０】
車線変更の意志を決定するロジックは、各車の希望走行車速と実際の走行車速とが比較され、実際の走行車速が希望車速よりも低い状態がある所定の水準以上に継続した場合に、車線変更の意志を固めるというロジックを設定する。具体的には、車線変更の意志を反映する変数ｚｉ（ｔ）が定義され、次式が計算される。
【０１５１】
【数６２】
2004152125.tif 000064
ただし、ｔ０は車両ｉが最後に車線変更した時刻かあるいは車両ｉが検出された時刻とする。ここで定義したｚｉ（ｔ）の値が、あるしきい値ｚｉｔｈｒよりも大きくなった場合に、車線変更の意志を固めて車線変更の可能性を判断するロジックを起動する。
【０１５２】
車線変更の判断は、車線変更先の先行車および後続車との車間時間をチェックし、車間時間が両方とも許容値よりも長い場合に車線変更を開始する。車線変更開始の条件が成立したら、横方向の速度指令値ｕｉｙにあらかじめ定められた典型的な車線変更パターンの時系列信号が与えられ、（６１）式に従って車線変更運動の計算が実行される。以上の車線変更のロジックを模式的に説明したものが図１１である。
【０１５３】
このような車線変更のロジックが加えられる以外は、車両情報抽出関数群３ａ−３を用いて予測モデルを結合して予測式を構成するという手順は同じである。ただ、他車両の横方向の速度が測定できる場合には、車両情報抽出関数群３ａ−３に改良を加えることができる。
【０１５４】
車両情報抽出関数群３ａ−３を構成する際に定義した車線判別関数（２４）式は、先行車であるための条件Ｂを、車両間の横方向の偏差が所定値以下であるという数学的な条件Ｂ’に置き換えられるものとして構成された関数であった。しかし、実際の運転においては、横方向の偏差だけでなく、横方向の速度も考慮して他車両の車線位置を判断している場合もあると考えられる。例えば、車線変更の動作が開始され、横方向の速度がある程度大きくなった場合には、車両自体がまだ車線内に進入してくるより前に、該当車両が車線変更してくることを前提とした運転に切り替えることがあり得る。従って、条件Ｂの数学的な表現は必ずしも条件Ｂ’だけに限定されるものではなく、以下に示すような条件Ｂ’’による表現も考えることができる。
【０１５５】
条件Ｂ’’：車両ｊと車両ｉのｌ秒後に予想される横方向位置の偏差が車線幅に比べて十分小さい。
【０１５６】
車両ｊと車両ｉの横方向の速度ｗｊとｗｉがｌ秒間にわたって一定に保たれると仮定すれば、車線判別関数は次のように再構成される。
【０１５７】
【数６３】
2004152125.tif 000065
（６３）式をプロットした三次元グラフを図１０に示す。
【０１５８】
このような車線判別関数を採用すると、車線変更してくる車両をより早く同一車線上の車両と判断して、迅速に割り込みに対応する操作を算出することが可能になる。
【０１５９】
（６３）式の車線判別関数を用いた場合、（２２）式の先行車抽出関数は（２４）式が（６３）式に置き換えられて、
【０１６０】
【数６４】
2004152125.tif 000066
となる。他の関数についても同様の置き換えによって関数を再構成することが可能である。
【０１６１】
第一の実施の形態の評価関数は、縦方向の位置だけを評価する評価項を構成したが、横方向の位置を評価する評価項を構成することも可能である。例えば、次のような評価関数を構成することが可能である。
【０１６２】
【数６５】
2004152125.tif 000067
ＬＸ、ＬＶは第一の実施の形態と同じく、（３３）、（３５）式を利用する。ＬＹについては、自車両が車線変更しないという仮定を置いているので、
【０１６３】
【数６６】
2004152125.tif 000068
と変更する。ＬＲＸは次のように定義する。
【０１６４】
【数６７】
2004152125.tif 000069
抽出関数ｘｐＬ（Ｘ）、ｘｆＬ（Ｘ）は、
【０１６５】
【数６８】
2004152125.tif 000070
【０１６６】
【数６９】
2004152125.tif 000071
と構成される。ＬＲＹは次のように定義する。
【０１６７】
【数７０】
2004152125.tif 000072
λは適当な正の定数である。抽出関数ｘｐＲ（Ｘ）、ｘｆＲ（Ｘ）は、（６８）、（６９）式における車線判別関数の引数ｙＬをｙＲに置き換えることで、それぞれ構成することができる。
【０１６８】
抽出関数ｘｐＲ（Ｘ）、ｘｆＲ（Ｘ）は、
【０１６９】
【数７１】
2004152125.tif 000073
【０１７０】
【数７２】
2004152125.tif 000074
と構成される。ただし、（７１）、（７２）式では、該当車両が存在しない場合の仮想車両は、車線の基準位置に配置されるものとして構成した。
【０１７１】
（７０）式のかわりに、次式のようなｌ秒後に予想される横方向の位置を評価する評価項も考えることが可能である。
【０１７２】
【数７３】
2004152125.tif 000075
この場合、抽出関数ＷｐＲ（Ｘ）、ＷｆＲ（Ｘ）は、
【０１７３】
【数７４】
2004152125.tif 000076
【０１７４】
【数７５】
2004152125.tif 000077
と構成される。ここでは、該当車両が存在しない場合の仮想車両は横方向の速度が０になるものとして構成した。
【０１７５】
推奨操作量と周囲車両の挙動予測を計算して得られる結果の例を図１２に示す。図１２の例では、他車３の接近に伴って自車両０は若干左方向に移動する操作量となり、ついで他車３が自車０を追い抜いて左車線に車線変更を始めると、他車３が左車線に進入する前から減速を指示する操作量が生成される。
【０１７６】
以上が本発明の第三の実施の形態である。
【０１７７】
上記のように、本実施の形態では、車線判別関数は、任意の検出車両と注目車両との横方向の位置の偏差、および横方向の相対速度に基づいて、走行車線位置の判別を行なうことを特徴とする。
このような構成により、横方向の相対的な位置の偏差に加えて横方向の相対速度の情報も用いて他車両の車線位置を判別しているので、他車両が車線変更してくる場合に、物理的に車線位置が変わるよりも早い段階で車線変更を検知して予測と評価に反映させることができるようになり、より実際の運転者の感覚に近い予測と評価を実現することができる。
【０１７８】
また、車線判別関数は、検出された車線ごとに基準となる横位置を設定し、任意の検出車両の横位置と車線ごとの基準位置の偏差、および検出車両の横速度に基づいて、走行車線位置の判別を行なうことを特徴とする。
このような構成により、横方向の絶対的な位置に加えて横方向の速度の情報も用いて他車両の車線位置を判別しているので、他車両が車線変更してくる場合に、物理的に車線位置が変わるよりも早い段階で車線変更を検知して予測と評価に反映させることができるようになり、より実際の運転者の感覚に近い予測と評価を実現することができる。
【０１７９】
また、周囲車両挙動予測手段３ａおよび評価関数構成手段３ｂの車両情報抽出関数群３ａ−３、３ｂ−３は、注目車両について予測および評価の対象となる特定の位置関係に他の検出車両が存在する場合には、該当する車両の縦方向の速度、横方向の位置、横方向の速度のいずれかを出力し、存在しない場合には、注目車両の縦方向の速度、横方向の位置、横方向の速度のいずれかにあらかじめ定められた値を加算して出力することを特徴とする。
このような構成により、車両の情報として縦方向の速度、横方向の位置、横方向の速度を出力することで、それらの情報を入力とする他車モデル３ａ−２やそれらの情報を評価する評価関数を用いることができるようになり、縦位置だけでは表現できない多様な予測と評価の方法を表現することができる。また、該当車両が存在しない場合には、自車の遠方に仮想的な車両が存在するものとし、その仮想的な車両が予測と評価に与える影響を小さくするように情報を割り当てることで、該当車両の有無によって予測式や評価関数を切り替える必要がなくなる。
【０１８０】
また、近傍情報抽出関数は、二車両の縦方向の位置を比較し、値が大きい方の車両の縦方向の速度、横方向の位置、横方向の速度のいずれかを出力する関数を、すべての検出車両に関して逐次的に適用することで実現することを特徴とする。
このような構成により、二車両の情報の比較を繰り返すことで、注目車両に最も近い先行車あるいは後続車に付随する縦方向の速度、横方向の位置、横方向の速度を系統的に抽出することができる。
【０１８１】
また、近傍情報抽出関数は、二車両の縦方向の位置を比較し、値が大きい方の車両の縦方向の速度、横方向の位置、横方向の速度のいずれかを出力する関数を微分可能な関数を用いて近似的に表現し、近似された関数をすべての検出車両に関して逐次的に適用することで実現することを特徴とする。
このような構成により、本来は不連続関数である近傍情報抽出関数を微分可能な２変数比較関数を逐次的に適用することで近似的に実現することによって、車両情報抽出関数全体を微分可能な関数で構成することが可能になり、推奨操作量を計算する場合に偏微分計算に基づく最適制御理論を活用することができるので、効率的な計算を行なうことができる。
【０１８２】
《第四の実施の形態》
本発明の第四の実施の形態を説明する。
本発明の第四の実施の形態における装置の基本的な構成は、第一の実施の形態と同様であり、図１に示した配置図と図２に示した構成図は同一である。
【０１８３】
本実施の形態では、自車モデルは、（２）、（３）式のように自車の操作量を入力とするモデルではなく、次のようなモデルを考える。
【０１８４】
【数７６】
2004152125.tif 000078
【０１８５】
【数７７】
2004152125.tif 000079
ここで、ｙｃｍｄはｙＬまたはｙＲのどちらかの値をとる変数であり、その遷移規則が次のように定められているとする。
【０１８６】
【数７８】
2004152125.tif 000080
また、他車の挙動は時系列信号として予測が与えられるものとする。すなわち、任意の未来の時刻ｔに関して、ｘｉ（ｔ）とｖｉ（ｔ）（ｉ＝１，２，３）の予測値を得ることができるとする。例えば、
【０１８７】
【数７９】
2004152125.tif 000081
のような時系列信号が得られているとする。ただし、ｔ０は現在時刻である。
【０１８８】
評価関数としては、（３５）式の速度評価項と（３６）式の他車評価項だけを用いて、
【０１８９】
【数８０】
2004152125.tif 000082
のように構成した被積分関数を用いて、
【０１９０】
【数８１】
2004152125.tif 000083
と構成する。
【０１９１】
推奨操作量は、自車モデル（７６）、（７７）、（７８）に含まれるパラメータを変更して、評価関数値が最も小さくなるパラメータから生成される操作量を推奨操作量として提示するものとする。ここでは、例えばｈ０とｖｄ０を可変のパラメータとして、考えられるパラメータの組をあらかじめメモリに記録しておく。すなわち、（ｈ０（ｋ）  ｖｄ０（ｋ））ｋ＝１，２，…，ｎｐというｎｐ組のパラメータを用意する。ｈ０（ｋ）とｖｄ０（ｋ）をパラメータとして設定したときの評価関数（８１）式の値をＪｋとする。
【０１９２】
【数８２】
2004152125.tif 000084
として最適なｋ＊が求まれば、ｈ０（ｋ＊）とｖｄ０（ｋ＊）をパラメータとしたときに、周囲車両群挙動予測式（７６）〜（７９）を計算して得られた軌道が推奨操作量になっている。具体的な操作量としては、（７６）の第二式の右辺を時系列信号としてまとめたものが推奨加減速パターン、変数ｙｃｍｄの時系列信号が推奨車線変更パターンとなっている。
【０１９３】
以上本発明を実施の形態に基づいて具体的に説明したが、本発明は上記実施の形態に限定されるものではなく、その要旨を逸脱しない範囲において種々変更可能であることは勿論である。
【図面の簡単な説明】
【図１】本発明の第一の実施の形態の一配置図である。
【図２】本発明の第一の実施の形態における演算部のソフトウェアブロックを示す図である。
【図３】本発明の第一の実施の形態の適用場面の一例を示す図である。
【図４】本発明の第一の実施の形態における予測式の構成方法を説明した図である。
【図５】本発明の第一の実施の形態における前後判別関数をプロットした図である。
【図６】本発明の第一の実施の形態における車線判別関数をプロットした図である。
【図７】本発明の第一の実施の形態における評価対象となる位置関係を示した図である。
【図８】本発明の第一の実施の形態における推奨操作量計算の結果の例を示した図である。
【図９】本発明の第一の実施の形態におけるフローチャートである。
【図１０】本発明の第三の実施の形態における車線判別関数をプロットした図である。
【図１１】本発明の第三の実施の形態における車線変更のロジックを説明した図である。
【図１２】本発明の第三の実施の形態における推奨操作量計算の結果の例を示した図である。
【符号の説明】
１ａ…前方レーダー
１ｂ…画像センサ
１ｃ…後方レーダー
１ｄ…側方センサ
２…車速センサ
３…演算部
３ａ…周囲車両挙動予測手段
３ａ−１…自車モデル
３ａ−２…他車モデル
３ａ−３…車両情報抽出関数群
３ｂ…評価関数構成手段
３ｂ−１…相互関係評価項
３ｂ−２…その他の評価項
３ｂ−３…車両情報抽出関数群
３ｃ…推奨操作量演算手段
４…表示装置
５…自車状態検出手段
６…周囲車両検出手段
７…周囲車両全情報
【図１】
2004152125.tif 000085
【図２】
2004152125.tif 000086
【図３】
2004152125.tif 000087
【図４】
2004152125.tif 000088
【図５】
2004152125.tif 000089
【図６】
2004152125.tif 000090
【図７】
2004152125.tif 000091
【図８】
2004152125.tif 000092
【図９】
2004152125.tif 000093
【図１０】
2004152125.tif 000094
【図１１】
2004152125.tif 000095
【図１２】
2004152125.tif 000096