(57)【要約】【課題】ドライバの運転適否状態を精度よく推定することができるようにする。【解決手段】アクセル開度センサ１２、車速センサ１４、及び車間距離計測部１６によって、自車両の走行状態及び操作状態を検出する。そして、運転信号生成部２４によって、複数の運転適否状態の各々について、運転適否状態における走行状態及び操作状態の時間変化を複数の運転モードの各々に対する線形モデルで表わした区分線形モデル、及び検出された走行状態及び操作状態に基づいて、次の時刻の複数の運転適否状態の各々に対する走行状態又は操作状態を推定する。そして、ドライバ状態推定部２６によって、運転信号生成部２４によって推定された複数の運転適否状態の各々に対する走行状態及び操作状態と、検出された走行状態及び操作状態とを比較し、比較結果に基づいて、ドライバの運転適否状態を推定する。【選択図】図１
2009154675.tif 000002
【特許請求の範囲】
【請求項１】
自車両の走行状態及びドライバが自車両を操作したときの操作状態の少なくとも一方の状態を検出する状態検出手段と、
前記状態検出手段によって検出された前記状態に基づいて、該状態が検出されたタイミングより後のタイミングにおける予め定められた複数のドライバの運転適否状態の各々に対する前記状態を推定する状態推定手段と、
前記状態推定手段によって推定された前記複数の運転適否状態に対する前記状態の各々と、前記状態検出手段によって該タイミングに検出された前記状態とを比較し、比較結果に基づいて、前記ドライバの運転適否状態を推定するドライバ状態推定手段と、
を含むドライバ状態推定装置。
【請求項２】
前記状態推定手段は、前記複数の運転適否状態の各々について、前記運転適否状態における前記状態の時間変化を複数の運転モードの各々に対する線形モデルで表わした区分線形モデル、及び前記状態検出手段によって検出された前記状態に基づいて、該タイミングにおける前記複数の運転適否状態の各々に対する前記状態を推定する請求項１記載のドライバ状態推定装置。
【請求項３】
前記状態推定手段は、前記複数の運転適否状態の各々について、前記状態検出手段によって検出された前記状態に対して、前記運転モードを同定し、前記同定された運転モードに対する線形モデルに基づいて、前記状態検出手段によって検出された前記状態から、該タイミングにおける前記状態を推定する請求項２記載のドライバ状態推定装置。
【請求項４】
前記状態検出手段によって検出される前記状態とは異なり、かつ、自車両の走行状態及びドライバが自車両を操作したときの操作状態の少なくとも一方の第２の状態を検出する第２状態検出手段を更に含み、
前記状態推定手段は、前記状態検出手段によって検出された前記状態、及び前記第２状態検出手段によって検出された前記第２の状態に基づいて、該状態が検出されたタイミングより後のタイミングにおける前記複数の運転適否状態の各々に対する前記状態を推定する請求項１記載のドライバ状態推定装置。
【請求項５】
前記状態推定手段は、前記複数の運転適否状態の各々について、前記運転適否状態における前記状態及び前記第２の状態の時間変化を複数の運転モードの各々に対する線形モデルで表わした区分線形モデル、前記状態検出手段によって検出された前記状態、及び前記第２状態検出手段によって検出された前記第２の状態に基づいて、該タイミングにおける前記複数の運転適否状態の各々に対する前記状態を推定する請求項４記載のドライバ状態推定装置。
【請求項６】
前記状態推定手段は、前記複数の運転適否状態の各々について、前記状態検出手段によって検出された前記状態及び前記第２状態検出手段によって検出された前記第２の状態に対して、前記運転モードを同定し、前記同定された運転モードに対する線形モデルに基づいて、前記状態検出手段によって検出された前記状態、及び前記第２状態検出手段によって検出された前記第２の状態から、該タイミングにおける前記状態を推定する請求項５記載のドライバ状態推定装置。
【請求項７】
前記状態は、ドライバが自車両を操作したときの操作状態である請求項１〜請求項６の何れか１項記載のドライバ状態推定装置。
【請求項８】
コンピュータを、
自車両の走行状態及びドライバが自車両を操作したときの操作状態の少なくとも一方の状態を検出する状態検出手段によって検出された前記状態に基づいて、該状態が検出されたタイミングより後のタイミングにおける予め定められた複数のドライバの運転適否状態の各々に対する前記状態を推定する状態推定手段、及び
前記状態推定手段によって推定された前記複数の運転適否状態に対する前記状態の各々と、前記状態検出手段によって該タイミングに検出された前記状態とを比較し、比較結果に基づいて、前記ドライバの運転適否状態を推定するドライバ状態推定手段
として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ドライバ状態推定装置及びプログラムに係り、特に、車両の走行状態やドライバの操作状態に基づいて、ドライバの運転適否状態を推定するドライバ状態推定装置及びプログラムに関する。
【背景技術】
【０００２】
従来より、加減速の滑らかさを表す車両のジャークを算出し、ジャークに関する正常運転時との異なりに基づいて、ドライバ状態を推定する運転状態推定装置が知られている（例えば、特許文献１）。
【０００３】
また、走行レーンの曲率から基準となるハンドル操作量を算出し、現在のハンドル操作量と基準操作量とのずれの標準偏差に基づいて、ドライバの覚醒度を推定する車両用覚醒度検出装置が知られている（例えば、特許文献２）。
【特許文献１】特開２００６−３３５２７７号公報
【特許文献２】特開平７−９８７９号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、上記の特許文献１、２に記載の技術では、現在の車両信号と、基準となる車両信号との異なりに基づいて、ドライバの状態を推定しており、基準とどれだけ異なっていれば異常と推定することが可能か、という点と、信号が基準と異なっている場合は必ず異常状態であるか、という点とに関して、明確な基準が与えられていないため、誤推定が生じてしまう、という問題があった。
【０００５】
本発明は、上記の問題点を解決するためになされたもので、ドライバの運転適否状態を精度よく推定することができるドライバ状態推定装置及びプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００６】
上記の目的を達成するために本発明に係るドライバ状態推定装置は、自車両の走行状態及びドライバが自車両を操作したときの操作状態の少なくとも一方の状態を検出する状態検出手段と、前記状態検出手段によって検出された前記状態に基づいて、該状態が検出されたタイミングより後のタイミングにおける予め定められた複数のドライバの運転適否状態の各々に対する前記状態を推定する状態推定手段と、前記状態推定手段によって推定された前記複数の運転適否状態に対する前記状態の各々と、前記状態検出手段によって該タイミングに検出された前記状態とを比較し、比較結果に基づいて、前記ドライバの運転適否状態を推定するドライバ状態推定手段とを含んで構成されている。
【０００７】
本発明に係るプログラムは、コンピュータを、自車両の走行状態及びドライバが自車両を操作したときの操作状態の少なくとも一方の状態を検出する状態検出手段によって検出された前記状態に基づいて、該状態が検出されたタイミングより後のタイミングにおける予め定められた複数のドライバの運転適否状態の各々に対する前記状態を推定する状態推定手段、及び前記状態推定手段によって推定された前記複数の運転適否状態に対する前記状態の各々と、前記状態検出手段によって該タイミングに検出された前記状態とを比較し、比較結果に基づいて、前記ドライバの運転適否状態を推定するドライバ状態推定手段として機能させるためのプログラムである。
【０００８】
本発明によれば、状態検出手段によって、自車両の走行状態及びドライバが自車両を操作したときの操作状態の少なくとも一方の状態を検出する。そして、状態推定手段によって、状態検出手段によって検出された状態に基づいて、該状態が検出されたタイミングより後のタイミングにおける予め定められた複数のドライバの運転適否状態の各々に対する状態を推定する。
【０００９】
そして、ドライバ状態推定手段によって、状態推定手段によって推定された複数の運転適否状態に対する状態の各々と、状態検出手段によって該タイミングに検出された状態とを比較し、比較結果に基づいて、ドライバの運転適否状態を推定する。
【００１０】
このように、複数の運転適否状態に対して推定された状態の各々と、検出された状態とを比較して、ドライバの運転適否状態を推定することにより、ドライバの運転適否状態を精度よく推定することができる。
【００１１】
本発明に係る状態推定手段は、複数の運転適否状態の各々について、運転適否状態における状態の時間変化を複数の運転モードの各々に対する線形モデルで表わした区分線形モデル、及び状態検出手段によって検出された状態に基づいて、該タイミングにおける複数の運転適否状態の各々に対する状態を推定することができる。これによって、複数の運転モードに対する線形モデルによって状態の時間変化を表わし、複数の運転適否状態の各々に対する状態を精度よく推定することができるため、ドライバの運転適否状態を精度よく推定することができる。
【００１２】
また、上記の状態推定手段は、複数の運転適否状態の各々について、状態検出手段によって検出された状態に対して、運転モードを同定し、同定された運転モードの線形モデルに基づいて、状態検出手段によって検出された状態から、該タイミングにおける状態を推定することができる。これによって、同定された運転モードの線形モデルに基づいて、後のタイミングにおける状態を精度よく推定することができる。
【００１３】
本発明に係るドライバ状態推定装置は、状態検出手段によって検出される状態とは異なり、かつ、自車両の走行状態及びドライバが自車両を操作したときの操作状態の少なくとも一方の第２の状態を検出する第２状態検出手段を更に含み、状態推定手段は、状態検出手段によって検出された状態、及び第２状態検出手段によって検出された第２の状態に基づいて、該状態が検出されたタイミングより後のタイミングにおける複数の運転適否状態の各々に対する状態を推定することができる。これによって、複数種類の状態を検出して、後のタイミングにおける複数の運転適否状態の各々に対する状態を精度よく推定することができる。
【００１４】
上記の第２状態検出手段を含む発明に係る状態推定手段は、複数の運転適否状態の各々について、運転適否状態における状態及び第２の状態の時間変化を複数の運転モードの各々に対する線形モデルで表わした区分線形モデル、状態検出手段によって検出された状態、及び第２状態検出手段によって検出された第２の状態に基づいて、該タイミングにおける複数の運転適否状態の各々に対する状態を推定することができる。これによって、区分線形モデルによって状態及び第２の状態の時間変化を表わし、複数の運転適否状態の各々に対する状態を精度よく推定することができるため、ドライバの運転適否状態を精度よく推定することができる。
【００１５】
また、上記の状態推定手段は、複数の運転適否状態の各々について、状態検出手段によって検出された状態及び第２状態検出手段によって検出された第２の状態に対して、運転モードを同定し、同定された運転モードに対する線形モデルに基づいて、状態検出手段によって検出された状態、及び第２状態検出手段によって検出された第２の状態から、該タイミングにおける状態を推定することができる。これによって、同定された運転モードの線形モデルに基づいて、後のタイミングにおける状態を精度よく推定することができる。
【００１６】
上記の状態検出手段によって検出する状態を、ドライバが自車両を操作したときの操作状態とすることが好ましい。
【発明の効果】
【００１７】
以上説明したように、本発明のドライバ状態推定装置及びプログラムによれば、複数の運転適否状態に対して推定された状態の各々と、検出された状態とを比較して、ドライバの運転適否状態を推定することにより、ドライバの運転適否状態を精度よく推定することができる、という効果が得られる。
【発明を実施するための最良の形態】
【００１８】
以下、図面を参照して本発明の実施の形態を詳細に説明する。なお、車両に搭載されたドライバ状態推定装置に本発明を適用した場合を例に説明する。
【００１９】
図１に示すように、第１の実施の形態に係るドライバ状態推定装置１０は、ドライバが自車両を操作したときの操作状態としてのアクセル開度を検出するアクセル開度センサ１２と、自車両の走行状態としての車速を検出する車速センサ１４と、走行状態としての車間距離を計測する車間距離計測部１６と、アクセル開度センサ１２、車速センサ１４、及び車間距離計測部１６からの出力に基づいて、ドライバの運転適否状態が正常運転状態及び異常運転状態の何れであるかを推定し、推定されたドライバの運転適否状態に応じて、表示装置４０に警告メッセージを表示させるコンピュータ１８とを備えている。
【００２０】
アクセル開度センサ１２は、検出したアクセル開度を示すアクセル開度信号を出力し、車速センサ１４は、検出した自車両の車速を示す車速信号を出力し、また、車間距離計測部１６は、計測した前方車両との車間距離を示す前方車間距離信号を出力する。また、アクセル開度センサ１２、車速センサ１４、及び車間距離計測部１６は、ある一定間隔で（例えば、１００ｍｓ毎）で、各信号をコンピュータ１８に入力する。
【００２１】
コンピュータ１８は、ＣＰＵと、ＲＡＭと、後述するドライバ状態推定処理ルーチンを実行するためのプログラムを記憶したＲＯＭとを備え、機能的には次に示すように構成されている。コンピュータ１８は、ドライバが正常運転状態であるときの運転行動をモデル化した行動モデルに関する情報を記憶した正常運転時行動モデル記憶部２０と、ドライバが異常運転状態として例えば居眠り状態であるときの運転行動をモデル化した行動モデルに関する情報を記憶した異常運転時行動モデル記憶部２２と、アクセル開度センサ１２から入力された時刻ｔのアクセル開度信号、車速センサ１４から入力された時刻ｔの車速信号、及び車間距離計測部１６から入力された時刻ｔの前方車間距離信号、正常運転状態時の行動モデル、及び異常運転状態時の行動モデルに基づいて、次の時刻ｔ＋１における正常運転状態時のアクセル開度信号と時刻ｔ＋１における異常運転状態時のアクセル開度信号とを生成する運転信号生成部２４と、アクセル開度センサ１２から入力されたアクセル開度信号、及び運転信号生成部２４によって生成されたアクセル開度信号とを比較して、ドライバの運転適否状態を推定するドライバ状態推定部２６と、ドライバ状態推定部２６によって推定されたドライバの運転適否状態に応じて、表示装置４０に警告メッセージを表示させる表示制御部２８とを備えている。
【００２２】
運転信号生成部２４は、以下に説明するように、次時刻のアクセル開度を算出して、アクセル開度信号を生成する。まず、入力された車速信号及び前方車間距離信号から得られる、走行状態としての前方車の面積変化率（ＫｄＢ）を求める。ＫｄＢは、非特許文献（津留等、“前後方向の接近に伴う危険状態評価に関する研究（第５報）”、自動車技術会２００７春季大会）に記載された接近離間状態を評価する指標であり、本実施の形態では、以下の（１）式により導出する。
【００２３】
【数１】
2009154675.tif 000003
【００２４】
次に、入力された車速信号が示す車速、入力された前方車間距離信号が示す車間距離、算出された面積変化率、及び入力されたアクセル開度信号が示すアクセル開度を、バッファ（図示省略）に格納する。ここで、バッファサイズをＴとし、バッファインデックスｉのバッファの値は、現在からＴ−ｉ＋１時刻前の各信号を表すものとする。
【００２５】
そして、バッファ中に含まれる信号から、正常運転状態時の行動モデルのパラメータ及び異常運転状態時の行動モデルのパラメータに基づいて、次時刻のアクセル開度信号を繰り返し生成し、生成されたアクセル開度信号系列の情報を正常時算出用バッファ（図示省略）及び異常時算出用バッファ（図示省略）に格納する。この生成処理では、バッファのインデックスｉの各信号の情報から、時刻ｉの次の時刻ｉ＋１のアクセル開度信号を求める処理を、時刻１〜Ｔ−１に対応するバッファのインデックス１〜Ｔ−１の各々に対して繰り返して行うことで、時刻２〜Ｔまでのアクセル開度信号系列を得て、アクセル開度信号系列の情報を正常時算出用バッファ及び異常時算出用バッファのインデックス２〜Ｔに格納する。
【００２６】
ここで、インデックスｉ＋１のアクセル開度信号を生成する処理は、バッファのインデックスｉに格納された信号の情報の組み合わせに対するモード同定と、同定された運転モードに応じた信号生成との２つのステップに分けられる。
【００２７】
モード同定は、後述する複数の運転モードに関する予め計算されたモード境界に基づいて、同定対象となる信号の組み合わせが、どの運転モードの領域に含まれるかを判定して、運転モードを同定する。具体的には、ｏｎｅ−ｖｓ−ｒｅｓｔ法に基づいたＳＶＭ識別器を設計し、多クラス識別問題を解くことによって、モード同定を行なう。すなわち信号の組み合わせを表わすベクトルｘに対して、以下の（２）式を最大にする運転モードｍを、同定対象の信号の組み合わせが属する運転モードと判定する。
【００２８】
【数２】
2009154675.tif 000004
【００２９】
ここで、Ｎは、サポートベクトルの個数、ｘｉ、ｍは、運転モードｍにおけるサポートベクトル、ｙｉ、ｍは、運転モードｍにおけるサポートベクトルの符号、αｉ、ｍは、運転モードｍにおける重み係数のベクトルである。サポートベクトルの符号は、ｏｎｅ−ｖｓ−ｒｅｓｔ法に従い、対象データｘが運転モードｍに属するとき＋１とし、運転モードｍに属しないときに−１とする。また、ｙｉ、ｍ、αｉ、ｍ、ｂは、後述するモード境界の推定によって求められる。
【００３０】
また、運転モードについて説明する。ドライバの運転行動は、複数の連続的な動作、及び入力情報に基づいた判断によってそれらの動作を切り替えることで実現されていると考えられる。そこで、本実施の形態では、ドライバの運転適否状態の各々について、運転適否状態における運転行動を、各信号の時間変化を複数の運転モードの各々に対する線形モデルで表わした区分線形モデルにおける複数の運転モードと、その運転モードの切り替えとによってモデル化する。
【００３１】
ここで、運転モードの分類について説明する。予め学習により、各信号の組み合わせの時間変化を線形モデルで表わし、線形モデルに関するパラメータを、クラスタリングによって分類し、分類に対応して複数の運転モードを定める。
【００３２】
以下に、クラスタリングによる運転モードの分類方法について説明する。まず、車速信号、前方車間距離信号、アクセル開度信号、及び面積変化率の組み合わせの時系列データＰが得られた場合、各データＰｉ（ｉ＝１、２、・・・、ｎ）に対し、距離の近いｃ個のデータを集め、局所集合ＬＤｓｉを作り、特徴量ξｉを得る。ここで、特徴量ξｉは、ＬＤｓｉ中のデータに対して、以下の（３）式で表される線形モデルを想定し、最小二乗推定を施すことにより得られる線形モデルに関するパラメータａｖ、ａｄ、ａｋ、ａｘと、ＬＤｓｉ中のデータの平均値から構成される。
ｘ（ｔ＋１）＝
ａｖ＊Ｖ（ｔ）＋ａｄ＊ｄ（ｔ）＋ａｋ＊ｋｄｂ（ｔ）＋ａｘ＊ｘ（ｔ）
・・・（３）
【００３３】
ここで、Ｖ（ｔ）が時刻ｔにおける車速、ｄ（ｔ）が時刻ｔにおける前方車間距離、ｋｄｂ（ｔ）が時刻ｔにおける面積変化率、ｘ（ｔ）が時刻ｔにおけるアクセル開度である。
【００３４】
そして、時系列データＰ中の各データＰｉ（ｉ＝１、２、・・・、ｎ）を、特徴量ξｉに変換することで、特徴空間における分布を得る。そして、特徴空間上のデータに対し、Ｋ−ｍｅａｎｓ法に基づいて、クラスタリングを行って、得られた複数のクラスタを、複数の運転モードとして決定する。
【００３５】
例えば、図２（Ａ）に示すように、ドライバの運転適否状態が正常運転状態である場合には、一般に運転者は、車速の上昇に合わせて車間距離が広くなるように調整して運転する。車速は、法廷速度に従い上限が定められ、この上限速度付近では、車速に対する車間距離のばらつきが大きくなる。
【００３６】
よって、おおまかには、上限速度付近で分類される２つの運転モードに従ってドライバは運転すると考えられるため、学習により、車速信号が示す車速及び前方車間距離信号が示す車間距離の組み合わせの時系列データの各データを、線形モデルで表わし、線形モデルに関するパラメータを、クラスタリングによって分類することにより、２つの運転モードに分類される。
【００３７】
また、図２（Ｂ）に示すように、ドライバが異常運転状態である場合については、正常運転状態である場合と同様に、上限速度付近において運転モードの変化が存在すると考えられる。学習により、車速信号が示す車速及び前方車間距離信号が示す車間距離の組み合わせの時系列データの各データを、線形モデルで表わし、線形モデルに関するパラメータを、クラスタリングによって分類することにより、上限速度付近で、２つの運転モードに分類される。なお、各運転モードに分類された線形モデルに関するパラメータが、正常運転状態と異なる。例えば、飲酒状態である場合には、車速が高いにもかかわらず車間距離をつめて走行すると考えられるため、各運転モードにおける車速に対する車間距離が小さくなる。
【００３８】
また、上記のように、クラスタリング（例えばｋ−ｍｅａｎｓクラスタリング）により分類された線形モデルに関するパラメータに基づいて、最小二乗法を用いて、各運転モードに対して、線形モデルに関するパラメータを求め、以下の（４）式から（６）式のように、区分線形モデルの形式で、運転行動をモデル化する。
運転モード１：
ｘ（ｔ＋１）＝
ａｖｍ１＊Ｖ（ｔ）＋ａｄｍ１＊ｄ（ｔ）＋ａｋｍ１＊ｋｄｂ（ｔ）＋ａｘｍ１＊ｘ（ｔ）
・・・（４）
運転モード２：
ｘ（ｔ＋１）＝
ａｖｍ２＊Ｖ（ｔ）＋ａｄｍ２＊ｄ（ｔ）＋ａｋｍ２＊ｋｄｂ（ｔ）＋ａｘｍ２＊ｘ（ｔ）
・・・（５）
運転モード３：
ｘ（ｔ＋１）＝
ａｖｍ３＊Ｖ（ｔ）＋ａｄｍ３＊ｄ（ｔ）＋ａｋｍ３＊ｋｄｂ（ｔ）＋ａｘｍ３＊ｘ（ｔ）
・・・（６）
【００３９】
ここで、ａｖｍｉ、ａｄｍｉ、ａｋｍｉ、及びａｘｍｉの各々は、運転モードｉに対する区分線形モデルに関するパラメータである。
【００４０】
また、運転モード間の境界の推定方法について説明する。本実施の形態では、ＳＶＭ（サポートベクターマシン）による境界の推定を行う。運転モード間の境界を簡略に表現するため、非線形式としての識別関数（各運転モード間の境界）は考えず、線形式としての識別関数を用いることとして、運転モード間の境界面を求める。以下に線形ＳＶＭの概要を示す。ＳＶＭは基本的に２クラスの識別器を構成するための学習法である。
【００４１】
例えば、運転モード１と運転モード２との間の境界を求める場合、境界を飛び越えて存在する学習データの数（例えば、運転モード２側に存在する運転モード１の学習データの数）が、最小となる境界面を求める。
【００４２】
これは、以下の（７）式、（８）式を制約条件とした、以下の（９）式に示すような２次計画問題として定式化することができる
αｉ，ｍＴｘｉ，ｍ−ｂ≧１  ｉｆ  ｙｉ，ｍ＝１    ・・・（７）
αｉ，ｍＴｘｉ，ｍ−ｂ≧−１  ｉｆ  ｙｉ，ｍ＝−１    ・・・（８）
【００４３】
【数３】
2009154675.tif 000005
【００４４】
ここで、ｘｉ，ｍは、運転モードｍに該当する学習データ、ｙｉ，ｍは、運転モードｍにおけるサポートベクトルの符号、αｉ，ｍは境界面の重み係数のベクトル、ｂは境界面の切片である。
【００４５】
次に、同定された運転モードにおける信号生成の方法について説明する。先に述べたように、運転行動を、区分線形の形式でモデル化するため、同定された運転モードに対する区分線形モデルに関するパラメータを用いて、上記（３）式に基づいて、時刻ｔの車速信号、時刻ｔの前方車間距離、時刻ｔにおける面積変化率、及び時刻ｔにおけるアクセル開度から、次の時刻ｔ＋１におけるアクセル開度信号を生成する。
【００４６】
上記の区分線形モデルのパラメータ及びモード同定に用いたＳＶＭのパラメータは、正常運転の行動モデル及び異常運転の行動モデルの各々について予め求められ、正常運転の行動モデルに対するパラメータは、正常運転時行動モデル記憶部２０に記憶され、異常運転の行動モデルに対するパラメータは、異常運転時行動モデル記憶部２２に記憶される。
【００４７】
ドライバ状態推定部２６は、運転信号生成部２４によって生成された正常時算出用バッファ中のインデックス２〜Ｔの正常時のアクセル開度信号（ｓｉｇ）と、バッファ（ｂｕｆ）中のインデックス２〜Ｔのアクセル開度信号との差の総和ｄｎを、以下の（１０）式によって計算する。
【００４８】
【数４】
2009154675.tif 000006
【００４９】
同様に、ドライバ状態推定部２６は、異常時算出用バッファ中のインデックス２〜Ｔの異常時のアクセル開度信号（ｓｉｇ）と、バッファ（ｂｕｆ）中のインデックス２〜Ｔのアクセル開度信号との差の総和ｄａを求める。そして、求めた差の総和ｄｎ、ｄａの大小比較から、ドライバの運転適否状態が正常運転状態および異常運転状態の何れであるかを推定する。
【００５０】
表示制御部２８は、ドライバ状態推定部２６によってドライバの運転適否状態が異常運転状態であると判定された場合に、異常運転状態であることを示す警告メッセージを表示装置４０に表示させる。
【００５１】
次に、第１の実施の形態に係るドライバ状態推定装置１０の作用について説明する。まず、以下に説明するような学習処理が行われる。正常運転に関する学習処理では、予め用意した正常運転時の学習データ（アクセル開度、車速、車間距離、及び面積変化率（ＫｄＢ）の組み合わせの学習データ）の時系列データから、各データについて、線形モデルに関するパラメータを求め、線形モデルに関するパラメータに対して、クラスタリングを行い、線形モデルに関するパラメータを複数の運転モードに分類する。
【００５２】
そして、分類された運転モードの各々について、運転モードに属する線形モデルに関するパラメータに基づいて、区分線形モデルに関するパラメータを求め、正常運転時行動モデル記憶部２０に記憶する。また、ＳＶＭによる境界推定を行い、モード同定に用いるＳＶＭのパラメータを求め、正常運転時行動モデル記憶部２０に記憶する。
【００５３】
また、異常運転に関する学習処理では、予め用意した異常運転時の学習データ（アクセル開度、車速、車間距離、及び面積変化率（ＫｄＢ）の組み合わせの学習データ）の時系列データに対して、各データについて、線形モデルに関するパラメータを求め、線形モデルに関するパラメータに対して、クラスタリングを行い、線形モデルに関するパラメータを複数の運転モードに分類する。
【００５４】
そして、分類された運転モードの各々について、運転モードに属する線形モデルに関するパラメータに基づいて、区分線形モデルに関するパラメータを求め、異常運転時行動モデル記憶部２２に記憶する。また、ＳＶＭによる境界推定を行い、モード同定に用いるＳＶＭのパラメータを求め、異常運転時行動モデル記憶部２２に記憶する。
【００５５】
そして、ドライバ状態推定装置１０を搭載した車両の走行中に、コンピュータ１８において、図３に示すドライバ状態推定処理ルーチンが実行される。
【００５６】
ステップ１０２において、アクセル開度センサ１２、車速センサ１４、及び車間距離計測部１６の各々から出力された信号を取得する。
【００５７】
そして、ステップ１０４では、上記ステップ１０２で取得した車速信号及び車間距離信号に基づいて、前方車の面積変化率（ＫｄＢ）を算出する。次のステップ１０６では、上記ステップ１０２で取得した車速信号、車間距離信号、及びアクセル開度信号と、上記ステップ１０４で算出した面積変化率とを、ＦＩＦＯバッファに格納する。
【００５８】
ステップ１１２では、変数ｔを初期値１に設定し、ステップ１１４で、バッファのインデックスｔから、アクセル開度信号、車速信号、車間距離信号、及び面積変化率を取得し、ステップ１１６において、時刻ｔ＋１における正常運転時のアクセル開度信号を生成する。上記ステップ１１６では、正常運転時行動モデル記憶部２０に記憶されたＳＶＭのパラメータを用いて、上記ステップ１１４で取得したアクセル開度信号、車速信号、車間距離信号、及び面積変化率に基づいて、上記（２）式を最大にする運転モードを決定して、モード同定を行う。そして、正常運転時行動モデル記憶部２０から、同定された運転モードに対応する区分線形モデルに関するパラメータを取得し、取得したパラメータを用いて、上記ステップ１１４で取得したアクセル開度信号、車速信号、車間距離信号、及び面積変化率に基づいて、上記（３）式に従って、次の時刻のアクセル開度信号を算出する。
【００５９】
次のステップ１１８において、上記ステップ１１６で生成されたアクセル開度信号を、正常時算出用バッファのインデックスｔ＋１に格納する。
【００６０】
そして、ステップ１２０において、時刻ｔ＋１における異常運転時のアクセル開度信号を算出して生成する。上記ステップ１２０では、異常運転時行動モデル記憶部２２に記憶されたＳＶＭのパラメータを用いて、上記ステップ１１４で取得したアクセル開度信号、車速信号、車間距離信号、及び面積変化率に基づいて、上記（２）式を最大にする運転モードを決定して、モード同定を行う。そして、異常運転時行動モデル記憶部２２から、同定された運転モードに対応する区分線形モデルに関するパラメータを取得し、取得したパラメータを用いて、上記ステップ１１４で取得したアクセル開度信号、車速信号、車間距離信号、及び面積変化率に基づいて、上記（３）式に従って、次の時刻のアクセル開度信号を算出する。
【００６１】
次のステップ１２２において、上記ステップ１２０で生成されたアクセル開度信号を、異常時算出用バッファのインデックスｔ＋１に格納する。
【００６２】
そして、ステップ１２４において、変数ｔが定数Ｔ−１以上であるか否かを判定し、変数ｔがＴ−１未満である場合には、ステップ１２６で、変数ｔをインクリメントして、ステップ１１４へ戻る。一方、上記ステップ１２４において、変数ｔがＴ−１に達した場合には、ステップ１２８へ移行する。
【００６３】
ステップ１２８では、上記（４）式に従って、正常時算出用バッファの値とバッファの値との差の総和ｄｎを算出し、ステップ１３０において、異常時算出用バッファの値とバッファの値との差の総和ｄａを算出する。
【００６４】
そして、ステップ１３２において、上記ステップ１２８で算出した差の総和ｄｎが、上記ステップ１３０で算出した差の総和ｄａより大きいか否かを判定し、ｄｎがｄａより大きい場合には、ドライバの運転適否状態が異常運転状態であると推定されるため、ステップ１３４で、ドライバの運転適否状態が異常運転状態であることを示す警告メッセージを表示装置４０に表示させて、ドライバ状態推定処理ルーチンを終了する。一方、上記ステップ１３２において、ｄｎがｄａ以下である場合には、ドライバの運転適否状態が正常運転状態であると推定されるため、ステップ１３６において、ドライバの運転適否状態が正常運転状態であることを示すメッセージを表示装置４０に表示させて、ドライバ状態推定処理ルーチンを終了する。
【００６５】
以上説明したように、第１の実施の形態に係るドライバ状態推定装置によれば、複数の運転適否状態に対して推定されたアクセル開度の各々と、検出されたアクセル開度とを比較して、ドライバの運転適否状態を推定することにより、ドライバの運転適否状態を精度よく推定することができる。
【００６６】
また、複数の運転モードに対する線形モデルによって、車速、アクセル開度、車間距離、及び面積変化率の時間変化を表わし、同定された運転モードの線形モデルに基づいて、複数の運転適否状態の各々における次時刻のアクセル開度を精度よく推定することができるため、ドライバの運転適否状態を精度よく推定することができる。
【００６７】
また、複数種類の操作状態又は走行状態として、車速、アクセル開度、車間距離、及び面積変化率を検出して、車速、アクセル開度、車間距離、及び面積変化率の組み合わせに対して運転モードを同定することにより、運転モードに対する線形モデルに従って、次時刻における複数の運転適否状態の各々に対するアクセル開度を精度よく推定することができる。
【００６８】
また、正常運転時および異常運転時における運転行動を、区分線形モデルによってモデル化する。そして、その区分線形モデルに関するパラメータから生成される正常運転状態におけるアクセル開度信号及び異常運転状態におけるアクセル開度信号の何れかが、実際の運転で検出されたアクセル開度信号に近いかを計算する。ドライバの運転適否状態が正常運転状態であればその時のアクセル開度信号は、正常運転状態に対する区分線形モデルから生成されたアクセル開度信号との差が小さくなることが期待されるため、生成されたアクセル開度信号と走行中に検出されたアクセル開度信号との差から、ドライバの運転適否状態を精度よく推定することが可能となる。
【００６９】
また、運転適否状態として、居眠り状態、飲酒状態、漫然状態等、正常とは異なるドライバの状態を、検出される走行状態や操作状態によって推定することが可能となる。
また、正常運転時の信号との差が大きかったとしても、異常運転時の信号との差が、それ以上に大きければ、ドライバの運転適否状態が正常運転状態であると推定することができるため、誤推定が生じることを防止することができる。
【００７０】
従来技術では、異常運転状態を、対象信号（例えばハンドル操作量）の標準偏差に基づいて判定していた。これは、時々刻々変化する多様な運転環境を、対象信号の変化としてのみ捕らえているに過ぎない。また、直線、カーブといったおおまかな走行路の区別だけでなく、前方車両が接近している場合、アクセルを強く踏み込んでいる場合などにも、ハンドル操作量の特性が変化することが予想されるが、ハンドル操作量の標準偏差では、そのような変化を捉えることはできない。一方、本実施の形態では、運転行動モデルのモデル化を、複数の連続的な動作と、入力情報に基づいた判断による、それらの動作の切り替えとで実現するため、ドライバの運転行動を的確にモデル化することができる。
【００７１】
なお、上記の実施の形態では、検出する操作状態又は走行状態の特徴量として、車速、車間距離、及びアクセル開度を採用した場合を例に説明したが、これに限定されるものではなく、他の特徴量の組み合わせであってもよい。例えば、アクセル開度、ブレーキ操作量、及びハンドル操作量などの操作状態の特徴量の組み合わせを検出するようにしてもよい。また、走行状態を表わす特徴量のみを検出するようにしてもよい。
【００７２】
また、１種類の特徴量だけを検出するようにしてもよい。例えば、アクセル開度を検出し、検出されたアクセル開度から、次の時刻のアクセル開度を推定するようにしてもよい。
【００７３】
また、１時刻の検出信号から次の時刻のアクセル開度を推定する場合を例に説明したが、検出信号の時系列データから、次の時刻のアクセル開度を推定するようにしても良い。
【００７４】
また、ドライバの運転適否状態が、正常運転状態及び異常運転状態の何れであるかを推定する場合を例に説明したが、ドライバの運転適否状態の異常運転状態として、居眠り状態や飲酒状態など複数種類の状態の何れかであるかを推定するようにしてもよい。この場合には、複数種類の異常運転状態の各々における運転行動を、区分線形モデルでモデル化すればよい。
【００７５】
次に、第２の実施の形態について説明する。なお、第１の実施の形態と同様の構成の部分については、同一符号を付して説明を省略する。
【００７６】
第２の実施の形態では、時刻Ｔまで先に予測してアクセル開度信号を生成している点が第１の実施の形態と異なっている。
【００７７】
図４に示すように、第２の実施の形態に係るドライバ状態推定装置２１０は、アクセル開度センサ１２と、車速センサ１４と、自車両の加速度を検出する加速度センサ２１６と、アクセル開度センサ１２、車速センサ１４、及び加速度センサ２１６からの出力に基づいて、ドライバの状態を推定し、推定されたドライバの運転適否状態に応じて、表示装置４０に警告メッセージを表示させるコンピュータ２１８とを備えている。
【００７８】
コンピュータ２１８は、正常運転時行動モデル記憶部２０と、異常運転時行動モデル記憶部２２と、アクセル開度センサ１２から入力された時刻ｔのアクセル開度信号、車速センサ１４から入力された時刻ｔの車速信号、及び加速度センサ２１６のから入力された時刻ｔの加速度信号、正常運転時の行動モデル、及び異常運転時の行動モデルに基づいて、時刻ｔ＋１以降における正常運転時のアクセル開度信号と時刻ｔ＋１以降における異常運転時のアクセル開度信号とを生成する運転信号生成部２２４と、アクセル開度センサ１２から入力されたアクセル開度信号、及び運転信号生成部２２４によって生成されたアクセル開度信号とを比較して、ドライバの運転適否状態を推定するドライバ状態推定部２２６と、表示制御部２８とを備えている。
【００７９】
運転信号生成部２２４は、以下に説明するように、次時刻以降のアクセル開度を予測して、アクセル開度信号を生成する。
【００８０】
まず、入力された車速信号が示す車速、入力された加速度信号が示す加速度、及び入力されたアクセル開度信号が示すアクセル開度を、正常時予測用バッファ及び異常時予測用バッファの各々の先頭のインデックス（インデックス１）に格納する。そして、インデックスｉの各信号に基づいて、次の時刻のアクセル開度信号を生成して、インデックスｉ＋１に格納することを繰り返して、時刻Ｔまで先に運転行動を予測する。
【００８１】
上記の生成する処理では、正常時予測用バッファ及び異常時予測用バッファの各々に対して、インデックスｉの各信号の情報から、時刻ｉ＋１のアクセル開度信号を求める処理を行い、時刻ｉ＋１のアクセル開度信号に基づいて、時刻ｉ＋１の車速信号及び加速度信号を生成して、インデックスｉ＋１に格納する。このように、予測された次時刻のアクセル開度信号に基づいて、同じ時刻の自車速や自車加速度を求めることができる。
【００８２】
ここで、インデックスｉ＋１のアクセル開度信号を生成する処理は、上記第１の実施の形態と同様に、インデックスｉに格納された各信号の情報のモード同定と、同定されたモードにおける信号生成との２つのステップに分けられる。
【００８３】
モード同定は、上記の第１の実施の形態と同様に、同定対象となる信号の組み合わせに対して、上記（２）式を最大にする運転モードｍを求めて、どの運転モードの領域に含まれるかを判定し、運転モードを同定する。
【００８４】
次に、同定された運転モードにおける信号生成では、当該運転モードに対応して定められた区分線形モデルに関するパラメータを用いて、次時刻のアクセル開度信号を生成する。すなわち、当該時刻ｔのアクセル開度信号が示すアクセル開度ｘ（ｔ）、車速信号が示す車速ｖ（ｔ）、加速度信号が示す加速度ａ（ｔ）に対して、以下の（１１）式により、次時刻のアクセル開度ｘ（ｔ＋１）を算出する。
ｘ（ｔ＋１）＝ａｖｍ＊ｖ（ｔ）＋ａａｍ＊ａ（ｔ）＋ａｘｍ＊ｘ（ｔ）
・・・（１１）
ここでａｖｍ、ａｄｍ、ａｘｍの各パラメータについては、例えばｋ−ｍｅａｎｓクラスタリングなどで得られた運転モード別の学習データ（車速、加速度、及びアクセル開度に関する学習データ）に対する線形モデルに関するパラメータから、最小二乗法を用いてあらかじめ求めておく。
【００８５】
上記の区分線形モデルに関するパラメータ及びモード同定に用いたＳＶＭのパラメータは、正常運転の行動モデル及び異常運転の行動モデルの各々について予め求められ、正常運転の行動モデルに対するパラメータは、正常運転時行動モデル記憶部２０に記憶され、異常運転の行動モデルに対するパラメータは、異常運転時行動モデル記憶部２２に記憶される。
【００８６】
ドライバ状態推定部２２６は、運転信号生成部２２４によって生成された正常時予測用バッファ中のインデックス２〜Ｔの正常時のアクセル開度信号と、アクセル開度センサ１２からインデックス２〜Ｔに対応するタイミングに入力されたアクセル開度信号との差の総和を、計算する。同様に、ドライバ状態推定部２２６は、異常時予測用バッファ中のインデックス２〜Ｔの異常時のアクセル開度信号と、アクセル開度センサ１２からインデックス２〜Ｔに対応するタイミングに入力されたアクセル開度信号との差の総和を、計算する。そして、求めた差の総和の大小比較から、ドライバの運転適否状態が正常運転状態および異常運転状態の何れであるかを推定する。
【００８７】
次に、第２の実施の形態に係るドライバ状態推定装置２１０の作用について説明する。まず、以下に説明するような学習処理が行われる。正常運転に関する学習処理では、予め用意した正常運転時の学習データ（アクセル開度、車速、及び加速度の組み合わせの学習データ）の時系列データから、各データについて、線形モデルに関するパラメータを求め、線形モデルに関するパラメータに対して、クラスタリングを行い、線形モデルに関するパラメータを複数の運転モードに分類する。
【００８８】
そして、分類された運転モードの各々について、運転モードに属する線形モデルに関するパラメータに基づいて、区分線形モデルに関するパラメータを求め、正常運転時行動モデル記憶部２０に記憶する。また、ＳＶＭによる境界推定を行い、モード同定に用いるＳＶＭのパラメータを求め、正常運転時行動モデル記憶部２０に記憶する。
【００８９】
また、異常運転に関する学習処理では、予め用意した異常運転時の学習データ（アクセル開度、車速、及び加速度の各々の学習データ）の組み合わせの時系列データに対して、各データについて、線形モデルに関するパラメータを求め、線形モデルに関するパラメータに対して、クラスタリングを行い、線形モデルに関するパラメータを分類して、分類に対応して複数の運転モードを決定する。
【００９０】
そして、分類された運転モードの各々について、クラスタリングされたパラメータに基づいて、区分線形モデルに関するパラメータを求め、異常運転時行動モデル記憶部２２に記憶する。また、ＳＶＭによる境界推定を行い、モード同定に用いるＳＶＭのパラメータを求め、異常運転時行動モデル記憶部２２に記憶する。
【００９１】
そして、ドライバ状態推定装置２１０を搭載した車両の走行中に、コンピュータ２１８において、図５に示すドライバ状態推定処理ルーチンが実行される。なお、第１の実施の形態と同様の処理については、同一符号を付して詳細な説明を省略する。
【００９２】
まず、ステップ２５０において、アクセル開度センサ１２、車速センサ１４、及び加速度センサ２１６の各々から出力された信号を取得し、ステップ２５２において、バッファのインデックス及び検出タイミングを示す変数ｔを初期値である１に設定する。
【００９３】
そして、ステップ２５４において、正常時予測用バッファ及び異常時予測用バッファの各々のインデックスｔ（初期値である１）に、上記ステップ２５０で取得した各信号を格納する。
【００９４】
次のステップ２５６では、正常時予測用バッファのインデックスｔから、アクセル開度信号、車速信号、及び加速度信号を取得し、ステップ２５８において、時刻ｔ＋１における正常運転時のアクセル開度信号を生成する。上記ステップ２５８では、正常運転時行動モデル記憶部２０に記憶されたＳＶＭのパラメータを用いて、上記ステップ２５６で取得したアクセル開度信号、車速信号、及び加速度信号に基づいて、上記（２）式を最大にする運転モードを決定して、モード同定を行う。そして、正常運転時行動モデル記憶部２０から、同定された運転モードに対応する区分線形モデルに関するパラメータを取得し、取得したパラメータを用いて、上記ステップ２５６で取得したアクセル開度信号、車速信号、及び加速度信号に基づいて、上記（１１）式に従って、次の時刻のアクセル開度信号を生成する。
【００９５】
そして、ステップ２６０では、上記ステップ２５６で取得した車速信号、加速度信号、及び上記ステップ２５８で生成された次時刻のアクセル開度信号に基づいて、次の時刻の車速信号及び加速度信号を生成し、ステップ２６２において、上記ステップ２５８で生成された次時刻のアクセル開度信号、上記ステップ２６０で生成された車速信号、及び加速度信号を、正常時予測用バッファのインデックスｔ＋１に格納する。
【００９６】
そして、ステップ２６４では、異常時予測用バッファのインデックスｔから、アクセル開度信号、車速信号、及び加速度信号を取得し、ステップ２６６において、時刻ｔ＋１における異常運転時のアクセル開度信号を生成する。上記ステップ２６６では、異常運転時行動モデル記憶部２２に記憶されたＳＶＭのパラメータを用いて、上記ステップ２６４で取得したアクセル開度信号、車速信号、及び加速度信号に基づいて、上記（２）式を最大にする運転モードを決定して、モード同定を行う。そして、正常運転時行動モデル記憶部２０から、同定された運転モードに対応する区分線形モデルに関するパラメータを取得し、取得したパラメータを用いて、上記ステップ２６４で取得したアクセル開度信号、車速信号、及び加速度信号に基づいて、上記（５）式に従って、次の時刻のアクセル開度信号を生成する。
【００９７】
そして、ステップ２６８では、上記ステップ２６４で取得した車速信号、加速度信号、及び上記ステップ２６６で生成された次時刻のアクセル開度信号に基づいて、次の時刻の車速信号及び加速度信号を生成し、ステップ２７０において、上記ステップ２６６で生成された次時刻のアクセル開度信号、上記ステップ２６８で生成された車速信号、及び加速度信号を、異常時予測用バッファのインデックスｔ＋１に格納する。
【００９８】
そして、ステップ２７２において、変数ｔが定数Ｔ−１以上であるか否かを判定し、変数ｔがＴ−１未満である場合には、ステップ２７４で、変数ｔをインクリメントして、ステップ２５６へ戻る。一方、上記ステップ２７２において、変数ｔがＴ−１に達した場合には、ステップ２７８へ移行する。
【００９９】
ステップ２７８において、アクセル開度センサ１２、車速センサ１４、及び加速度センサ２１６の各々から出力された信号を取得する。
【０１００】
そして、ステップ２８０において、上記ステップ２７８で取得した車速信号、加速度信号、及びアクセル開度信号を、ＦＩＦＯバッファに格納する。
【０１０１】
ステップ２８２では、正常時予測用バッファの値とバッファの値との差の総和ｄｎを算出し、ステップ２８４において、異常時予測用バッファの値とバッファの値との差の総和ｄａを算出する。
【０１０２】
そして、ステップ１３２において、上記ステップ２８２で算出した差の総和ｄｎが、上記ステップ２８４で算出した差の総和ｄａより大きいか否かを判定し、ｄｎがｄａより大きい場合には、ステップ１３４で、ドライバの運転適否状態が異常運転状態であることを示す警告メッセージを表示装置４０に表示させて、ドライバ状態推定処理ルーチンを終了する。一方、ｄｎがｄａ以下である場合には、ステップ１３６において、ドライバの運転適否状態が正常運転状態であることを示すメッセージを表示装置４０に表示させて、ドライバ状態推定処理ルーチンを終了する。
【０１０３】
以上説明したように、第２の実施の形態に係るドライバ状態推定装置によれば、複数の運転適否状態の各々に対して予測されたアクセル開度の時系列データと、検出されたアクセル開度の時系列データとを比較して、ドライバの運転適否状態を推定することにより、ドライバの運転適否状態を精度よく推定することができる。
【図面の簡単な説明】
【０１０４】
【図１】本発明の第１の実施の形態に係るドライバ状態推定装置の構成を示す概略図である。
【図２】（Ａ）正常運転状態時の車速及び車間距離の時間変化を示すグラフ、及び（Ｂ）異常運転状態時の車速及び車間距離の時間変化を示すグラフである。
【図３】本発明の第１の実施の形態に係るドライバ状態推定装置におけるドライバ状態推定処理ルーチンの内容を示すフローチャートである。
【図４】本発明の第２の実施の形態に係るドライバ状態推定装置の構成を示す概略図である。
【図５】本発明の第２の実施の形態に係るドライバ状態推定装置におけるドライバ状態推定処理ルーチンの内容を示すフローチャートである。
【符号の説明】
【０１０５】
１０、２１０  ドライバ状態推定装置
１２   アクセル開度センサ
１４   車速センサ
１６   車間距離計測部
１８、２１８  コンピュータ
２０   正常運転時行動モデル記憶部
２２   異常運転時行動モデル記憶部
２４、２２４  運転信号生成部
２６、２２６  ドライバ状態推定部
２１６  加速度センサ
【図１】
2009154675.tif 000007
【図２】
2009154675.tif 000008
【図３】
2009154675.tif 000009
【図４】
2009154675.tif 000010
【図５】
2009154675.tif 000011