JP6075453B2 - 運転状態推定装置 - Google Patents

Description

本発明は、運転状態推定装置に関する。

従来、運転状態推定装置としては、例えば、特許文献１に記載の技術がある。
特許文献１に記載の技術では、ステアリングホイールの操舵角に基づき、操舵操作が滑らかに行われたと仮定した場合の操舵角推定値と実際の操舵角との差（以下、「分布データ」とも呼ぶ）を算出する。続いて、算出した分布データに基づき、普段の操舵特性に対応する相対的に時間的範囲の長い走行状態分布と、現在の操舵特性に対応する相対的に時間的範囲の短い走行状態分布とを更新する。そして、更新した２つの走行状態分布の分布間の相違量に基づき運転者の運転状態を推定する。これにより、特許文献１に記載の技術では、交通環境の違いによらず、運転の不安定状態を精度よく検出可能となっている。
ここで、特許文献１に記載の技術では、２つの走行状態分布や運転者の運転状態の推定等、分布データに基づく演算を分布データが算出するたびに実行していた。

特開２００９−９４９５号公報

しかしながら、上記特許文献１に記載の技術では、２つの走行状態分布や運転者の運転状態の推定等、分布データに基づく演算を分布データが算出するたびに実行していた。それゆえ、運転状態推定装置の演算負荷が増大する可能性があった。そのため、スマートフォン、廉価な車載コントローラ等、演算能力が比較的低い機器で運転状態推定装置を実現することが困難であった。
本発明は、上記のような点に着目したもので、運転状態推定装置の演算負荷を低減可能とすることを目的とする。

上記課題を解決するために、本発明の一態様では、走行状態データに基づき走行状態分布用の分布データを予め定めた第１周期毎に取得する。続いて、取得した分布データに基づき第１周期よりも長い第２周期毎に時間的範囲の異なる複数の走行状態分布を更新する。続いて、更新した複数の走行状態分布に基づき運転者の運転状態を推定する。続いて、運転状態の推定に対して外乱となる運転状況にあると判定すると、下記（ａ）〜（ｅ）の処理の少なくともいずれか１つの処理を実行する。
（ａ）前記走行状態分布更新部による走行状態分布の更新を中断する更新中断処理
（ｂ）前記走行状態分布更新部が更新した走行状態分布を前記外乱となる運転状況にあると判定しているときの分布データを含まない走行状態分布で置き換えるリストア処理
（ｃ）前記走行状態分布更新部が更新した走行状態分布のうち、予め定めた相対的に時間範囲の短い第２走行状態分布を、前記第２走行状態分布よりも時間的範囲が長い第１走行状態分布で置き換えるリセット処理
（ｄ）前記データ一時蓄積部が蓄積している分布データのうち前記外乱となる運転状況にあると判定する前までの分布データで走行状態分布を更新する更新周期変更処理
（ｅ）前記外乱となる運転状況から前記外乱とならない運転状況に変化したと判定すると前記走行状態分布更新部による走行状態分布の更新の中断を解除しその時点を前記第２周期の開始点に変更するレジューム処理

本発明の一態様によれば、予め定めた第１周期毎に分布データを取得するとともに、第１周期よりも長い第２周期毎に複数の走行状態分布を更新する。それゆえ、複数の走行状態分布の更新や運転状態の推定等、分布データに基づく演算の実行頻度を低減できる。これにより、運転状態推定装置の演算負荷を低減できる。

運転状態推定装置を搭載した車両の構成を表す図である。 運転状態推定装置のシステム構成例を表すブロック図である。 運転支援部１００Ａの構成を表すブロック図である。 車両の運転状況を説明するための図である。 走行状態分布更新部１３０の構成を表すブロック図である。 運転不安定度判定処理を表すフローチャートである。 絶対エントロピーＨp1、Ｈp2、相対エントロピーＲＨp算出に用いる記号を説明するための図である。 操舵角予測誤差θeを説明するための図である。 予測誤差区分ｂｉの範囲を説明するための図である。 第２周期の開始点を説明するための図である。 第１走行状態分布および第２走行状態分布の更新方法を説明するための図である。 第１走行状態分布および第２走行状態分布の更新方法を説明するための図である。 第１走行状態分布および第２走行状態分布の特性を説明するための図である。 相対エントロピーＲＨpを説明するための図である。 情報呈示オン例を説明するための図である。 情報呈示オン例を説明するための図である。 運転状態推定装置を搭載した車両の動作を表す図である。 運転状態推定装置を搭載した車両の構成を表す図である。 運転状態推定装置を搭載した車両の構成を表す図である。

まず、本発明に係る実施形態について図面を参照しつつ説明する。
（構成）
図１は、本実施形態の運転状態推定装置を搭載した車両の構成を表す図である。
図１に示すように、車両は、アクセルペダル開度量センサ１、ブレーキペダル操作量センサ２、操舵角センサ３、車輪速センサ４、ウインカ検出センサ５、およびナビゲーション装置６を備える。また、車両は、Ｇセンサ７、シフトセンサ８、前方車両検出装置９、およびコントローラ１００を備える。

アクセルペダル開度量センサ１は、アクセルペダルの開度量を検出する。そして、アクセルペダル開度量センサ１は、検出した開度量をコントローラ１００に出力する。
ブレーキペダル操作量センサ２は、ブレーキペダルの操作量を検出する。そして、ブレーキペダル操作量センサ２は、検出した操作量をコントローラ１００に出力する。
操舵角センサ３は、ステアリングホイール（不図示）の操舵角を検出する。そして、操舵角センサ３は、検出した操舵角をコントローラ１００に出力する。操舵角センサ３としては、例えば、ステアリングコラムの回転角を検出する角度センサを採用できる。

車輪速センサ４は、車輪の回転数（以下、「車輪速」とも呼ぶ）を検出する。続いて、車輪速センサ４は、検出した車輪速に基づいて車速を算出する。そして、車輪速センサ４は、検出した車輪速および算出した車速のそれぞれをコントローラ１００に出力する。
ウインカ検出センサ５は、ウインカレバー（不図示）の操作状態（以下、「ウインカ操作」とも呼ぶ）を検出する。ウインカ操作としては、例えば、操作の有無がある。そして、ウインカ検出センサ５は、検出したウインカ操作をコントローラ１００に出力する。

シフトセンサ８は、シフトレバー（不図示）の操作状態（以下、「シフト操作」とも呼ぶ）を検出する。シフト操作としては、例えば、Ｐ、Ｄ、Ｒ等のシフトレバーの位置がある。そして、シフトセンサ８は、検出したシフト操作をコントローラ１００に出力する。
情報呈示装置は、コントローラ１００が出力した制御信号（後述）に従って、警報その他の情報を運転者に呈示する。呈示方法としては、音声や画像がある。情報呈示装置としては、例えば、ブザー音や音声により運転者への情報提供を行うスピーカ１０、および画像やテキストの表示により運転者への情報提供を行う表示ユニットを採用できる。表示ユニットとしては、例えば、ナビゲーション装置６の表示モニタを流用してもよい。

ナビゲーション装置６は、ＧＰＳ(Global Positioning System)受信機、地図データベース、および表示モニタを備える。そして、ナビゲーション装置６は、ＧＰＳ受信機および地図データベースから車両の現在位置および道路情報を取得する。続いて、ナビゲーション装置６は、取得した車両の現在位置および道路情報に基づいて車両が走行する道路の種別や道路幅員等の各種情報を取得する。続いて、ナビゲーション装置６は、取得した情報に基づいて経路探索の結果および経路案内の結果等を表示モニタに表示する。

Ｇセンサ７は、車両に発生した前後加速度および横加速度を検出する。そして、Ｇセンサ７は、検出した前後加速度および横加速度をコントローラ１００に出力する。
前方車両検出装置９は、車両の進行方向前方に存在する他の車両その他の障害物の情報（例えば、障害物までの距離）を検出する。そして、前方車両検出装置９は、検出した情報をコントローラ１００に出力する。前方車両検出装置９としては、例えば、車両の進行方向前方にレーザー光を出射して反射光を検出するレーザ距離計を採用できる。

コントローラ１００は、ＣＰＵ(Central Processing Unit)、並びにＲＯＭ(Read Only Memory)、ＲＡＭ(Random Access Memory)およびＡ/Ｄ(Analog to Digital)変換回路等のＣＰＵ周辺部品を備える。そして、コントローラ１００（ＣＰＵ、ＣＰＵ周辺部品）は、運転不安定度判定処理を行う運転支援部１００Ａを備える。運転不安定度判定処理では、運転支援部１００Ａは、アクセルペダル開度量センサ１、およびブレーキペダル操作量センサ２等が出力した検出結果に基づいて、運転者が操作可能な運転操作子の操作状態、および車両状態の少なくとも一方を含む走行状態データを取得する。運転操作子としては、例えば、ステアリングホイール、アクセルペダル、およびブレーキペダルがある。車両状態としては、前方車両に対する車間情報がある。本実施形態では、走行状態データとして、操舵角センサ３が出力した操舵角の情報（以下、「操舵角情報」とも呼ぶ）を採用する。

続いて、運転支援部１００Ａは、取得した走行状態データ（操舵角情報）に基づき複数の走行状態分布（第１走行状態分布（後述）、第２走行状態分布（後述））用の分布データ（操舵角予測誤差θe（後述））を予め定めた第１周期（例えば、５０ミリ秒）毎に取得する。続いて、運転支援部１００Ａは、取得した分布データ（操舵角予測誤差θe）に基づき第１周期（例えば、５０ミリ秒）よりも長い第２周期（例えば、５秒）毎に時間的範囲の異なる複数の走行状態分布（第１走行状態分布、第２走行状態分布）を更新する。

続いて、運転支援部１００Ａは、更新した複数の走行状態分布（第１走行状態分布、第２走行状態分布）に基づいて絶対エントロピーＨp1、Ｈp2（後述）を算出する。続いて、続いて、運転支援部１００Ａは、算出した絶対エントロピーＨp1、Ｈp2に基づいて運転者の運転状態（運転の不安定度（後述））を判定する。続いて、運転支援部１００Ａは、更新した複数の走行状態分布（第１走行状態分布、第２走行状態分布）の分布間の相違量（相対エントロピーＲＨp（後述））に基づいて運転者の運転状態（運転の不安定度）を推定する。そして、運転支援部１００Ａは、推定した運転状態（運転の不安定度）に基づいて警報その他の情報（以下、「呈示情報」とも呼ぶ）を運転者に呈示させる制御信号を情報呈示装置に出力する。これにより、運転支援部１００Ａは、運転者に呈示情報を呈示し、運転の不安定度（運転の不安定状態）について運転者の注意を喚起する。

なお、走行状態データとしては、前方車両に対する車間情報（車間距離、車間時間）や、アクセルペダルやブレーキペダルの操作に基づく加減速情報等を採用してもよい。車間情報（車間距離、車間時間）や、加減速情報等を採用する場合、走行状態分布（第１走行状態分布、第２走行状態分布）および分布間の相違量（相対エントロピーＲＨp）の算出は、例えば、国際公開番号ＷＯ２００９／０１３８１５（特願２００９−５２４３４２号）の公報等に記載しているような公知の方法によって算出すればよい。

コントローラ１００としては、例えば、スマートフォンや、車載コントローラ、クラウドサーバを採用してもよい。クラウドサーバを採用する場合、車両が、アクセルペダル開度量センサ１、およびブレーキペダル操作量センサ２等が出力した検出結果をクラウドサーバに送信する。これにより、クラウドサーバが、車両が送信した検出結果に基づいて運転不安定度判定処理を実行し、運転不安定度判定処理の結果に基づいて呈示情報を運転者に呈示させる制御信号を情報呈示装置に送信する。そして、車両が、クラウドサーバが送信した制御信号に基づいて情報呈示装置で運転者に呈示情報を呈示する。

図２は、本実施形態の運転状態推定装置のシステム構成例を表すブロック図である。
図２に示すように、本実施形態では、情報呈示装置として、視覚情報呈示装置、および聴覚情報呈示装置を例示している。また、視覚情報呈示装置としては、ナビゲーション装置６の表示モニタ、聴覚情報呈示装置としては、スピーカ１０を例示している。
図３は、本実施形態の運転支援部１００Ａの構成を表すブロック図である。
図３に示すように、運転支援部１００Ａは、走行状態データ取得部１１０、運転状況判定部１２０、走行状態分布更新部１３０、運転不安定度判定部１４０、および情報呈示部１５０を備える。
走行状態データ取得部１１０は、操舵角センサ３が出力した検出結果を取得する。そして、走行状態データ取得部１１０は、取得した検出結果を走行状態データとする。

図４は、車両の運転状況を説明するための図である。
運転状況判定部１２０は、アクセルペダル開度量センサ１、およびブレーキペダル操作量センサ２等が出力した検出結果に基づいて、車両の運転状況（第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、第４外乱運転状況、通常運転状況（後述））を判定する。具体的には、運転状況判定部１２０は、アクセルペダル開度量センサ１、およびブレーキペダル操作量センサ２等が出力した検出結果に基づいて、運転者が操作可能な運転操作子の操作状態、走行環境、および車両状態を検出する。続いて、運転状況判定部１２０は、検出した運転操作子の操作状態、走行環境、および車両状態に基づいて、車両の運転状況が第１外乱運転状況（運転状態の推定に対して、外乱となる運転状況）にあるかを判定する。第１外乱運転状況としては、例えば、車両の加減速、ブレーキペダル（不図示）の操作、スイッチ・レバー操作等のイベントで発生する運転状況がある。

車両の加減速は、例えば、アクセルペダル開度量センサ１、ブレーキペダル操作量センサ２、Ｇセンサ７（前後Ｇ）に基づき検出する。ブレーキペダル（不図示）の操作は、例えば、ブレーキペダル操作量センサ２に基づき検出する。スイッチ・レバー操作は、例えば、ワイパー（不図示）の作動状態を検出するワイパーセンサ（不図示）、ライト（不図示）の作動状態を検出するライトセンサ（不図示）、オーディオ装置（不図示）の操作状態を検出するオーディオセンサ（不図示）に基づき検出する。

また、運転状況判定部１２０は、検出した運転操作子の操作状態、走行環境、および車両状態に基づいて、車両の運転状況が第２外乱運転状況（運転状態の推定に対して、外乱となる運転状況）にあるかを判定する。第２外乱運転状況としては、例えば、横滑り、ＶＤＣ(Vehicle Dynamics Control)作動、ＬＤＰ(Lane Departure Preservation)作動、および各種事象（車線変更、車両の右左折、車両の加減速、ブレーキペダル（不図示）の操作、変速操作、スイッチ・レバー操作、トンネル出入口、うねり、路面のジョイント）が設定時間継続等のイベントで発生する運転状況がある。ここで、ＶＤＣ作動は、例えば、ＶＤＣの作動を表すＶＤＣ作動フラグのフラグ値に基づき検出する。ＬＤＰ作動は、例えば、ＬＤＰの作動を表すＬＤＰ作動フラグのフラグ値に基づき検出する。

また、運転状況判定部１２０は、検出した運転操作子の操作状態、走行環境、および車両状態に基づいて、車両の運転状況が第３外乱運転状況（運転状態の推定に対して、外乱となる運転状況）にあるかを判定する。第３外乱運転状況としては、例えば、車線変更、車両の右左折、変速操作、トンネル出入口、およびうねり等のイベントで発生する運転状況がある。車線変更、車両の右左折は、例えば、ウインカ検出センサ５に基づき検出する。変速操作は、例えば、クラッチの操作状態を検出するクラッチセンサ（不図示）、シフトセンサ８に基づき検出する。トンネル出入口は、例えば、ヘッドライト（不図示）の作動状態を検出するライトセンサ（不図示）に基づき検出する。うねりは、例えば、Ｇセンサ７（横Ｇ）に基づき検出する。さらに、運転状況判定部１２０は、検出した運転操作子の操作状態、走行環境、および車両状態に基づいて、車両の運転状況が第４外乱運転状況（運転状態の推定に対して、外乱となる運転状況）にあるかを判定する。第４外乱運転状況としては、例えば、路面のジョイント等のイベントで発生する運転状況がある。路面のジョイントは、例えば、車輪速センサ４（車輪速変化）に基づき検出する。
一方、運転状況判定部１２０は、第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、第４外乱運転状況のいずれでもないと判定した場合には、運転の不安定度の判定に対して外乱となる運転状況にない（以下、「通常運転状況」とも呼ぶ）と判定する。

走行状態分布更新部１３０は、走行状態データ取得部１１０が取得した走行状態データ（操舵角情報）に基づいて走行状態分布（第１走行状態分布、第２走行状態分布）用の分布データ（操舵角予測誤差θe）を取得する。操舵角予測誤差θeの取得周期は、予め定めた第１周期（例えば、５０ミリ秒）とする。続いて、走行状態分布更新部１３０は、取得した分布データ（操舵角予測誤差θe）に基づいて時間的範囲の異なる複数の走行状態分布（第１走行状態分布、第２走行状態分布）を更新する。第１走行状態分布および第２走行状態分布の更新周期は、第１周期（例えば、５０ミリ秒）よりも長い第２周期（例えば、５秒）とする。具体的には、走行状態分布更新部１３０は、取得した分布データ（操舵角予測誤差θe）に基づき、予め定めた相対的に長い時間的範囲（例えば、２１６０秒）で取得した分布データ（操舵角予測誤差θe）の度数分布を第１走行状態分布とする。
また、走行状態分布更新部１３０は、取得した分布データ（操舵角予測誤差θe）に基づき、第１走行状態分布（例えば、２１６０秒）よりも短い時間的範囲（例えば、９０秒）で取得した分布データ（操舵角予測誤差θe）の度数分布を第２走行状態分布とする。第１走行状態分布および第２走行状態分布の算出例については、後述する。

図５は、本実施形態の走行状態分布更新部１３０の構成を表すブロック図である。
図５に示すように、走行状態分布更新部１３０は、分布データ算出部１３０Ａ、データ一時蓄積部１３０Ｂ、第１走行状態分布更新部１３０Ｃ、第２走行状態分布更新部１３０Ｄ、分布蓄積部１３０Ｅ、および分布選択部１３０Ｆを備える。
分布データ算出部１３０Ａは、走行状態データ取得部１１０が取得した走行状態データ（操舵角情報（操舵角θ））に基づき分布データ（操舵角予測誤差θe）を算出する。操舵角予測誤差θeは、予め定めた第１周期（例えば、５０ミリ秒）毎に算出する。
データ一時蓄積部１３０Ｂは、分布データ算出部１３０Ａが算出した分布データ（操舵角予測誤差θe）を取得する。そして、データ一時蓄積部１３０Ｂは、取得した分布データ（操舵角予測誤差θe）を蓄積する。分布データ（操舵角予測誤差θe）は、第１周期（例えば、５０ミリ秒）毎にデータ一時蓄積部１３０Ｂに蓄積する。

なお、データ一時蓄積部１３０Ｂは、運転状況判定部１２０が第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、および第４外乱運転状況のいずれかの運転状況にあると判定すると、分布データ（操舵角予測誤差θe）のデータ一時蓄積部１３０Ｂへの蓄積を中断する。これにより、データ一時蓄積部１３０Ｂは、第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、および第４外乱運転状況にあると判定していない期間、つまり、通常運転状況にあると判定している期間の分布データ（操舵角予測誤差θe）のみを蓄積する。また、データ一時蓄積部１３０Ｂは、運転状況判定部１２０が第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、および第４外乱運転状況のいずれかの運転状況から通常運転状況に変化したと判定すると、分布データ（操舵角予測誤差θe）のデータ一時蓄積部１３０Ｂへの蓄積を再開する。その際、データ一時蓄積部１３０Ｂは、データ一時蓄積部１３０Ｂが蓄積している分布データ（操舵角予測誤差θe）を破棄する。

第１走行状態分布更新部１３０Ｃは、データ一時蓄積部１３０Ｂが蓄積している分布データ（操舵角予測誤差θe）に基づいて第１走行状態分布を更新する。また、第１走行状態分布更新部１３０Ｃは、データ一時蓄積部１３０Ｂが蓄積している分布データ（操舵角予測誤差θe）のうち、運転状況判定部１２０が第４外乱運転状況にあると判定する前までの分布データ（操舵角予測誤差θe）で第１走行状態分布を更新する（以下、「更新周期変更処理」とも呼ぶ）。これにより、第１走行状態分布更新部１３０Ｃは、運転状況判定部１２０が第４外乱運転状態にあると判定すると、通常運転状況にあると判定している期間の分布データ（操舵角予測誤差θe）のみに基づいて第１走行状態分布を更新する。第１走行状態分布は、第２周期（＞第１周期（５０ミリ秒）。例えば、５秒）毎に更新する。
なお、第１走行状態分布更新部１３０Ｃは、運転状況判定部１２０が第１外乱運転状況にあると判定すると、第１走行状態分布更新部１３０Ｃによる第１走行状態分布の更新を中断する（以下、「更新中断処理」とも呼ぶ）。具体的には、第１走行状態分布更新部１３０Ｃは、分布蓄積部１３０Ｅが蓄積している第１走行状態分布のうち、１周期前に更新した第１走行状態分布を現在の第１走行状態分布として設定する。

第２走行状態分布更新部１３０Ｄは、データ一時蓄積部１３０Ｂが蓄積している分布データ（操舵角予測誤差θe）に基づいて第２走行状態分布を更新する。また、第２走行状態分布更新部１３０Ｄは、データ一時蓄積部１３０Ｂが蓄積している分布データ（操舵角予測誤差θe）のうち、運転状況判定部１２０が第４外乱運転状況にあると判定する前までの分布データ（操舵角予測誤差θe）で第２走行状態分布を更新する（更新周期変更処理）。これにより、第２走行状態分布更新部１３０Ｄは、運転状況判定部１２０が第４外乱運転状態にあると判定すると、通常運転状況にあると判定している期間の分布データ（操舵角予測誤差θe）のみに基づいて第２走行状態分布を更新する。第２走行状態分布は、第２周期（＞第１周期（５０ミリ秒）。例えば、５秒）毎に更新する。

なお、第２走行状態分布更新部１３０Ｄは、運転状況判定部１２０が第１外乱運転状況にあると判定すると、第２走行状態分布更新部１３０Ｄによる第２走行状態分布の更新を中断する（更新中断処理）。具体的には、第２走行状態分布更新部１３０Ｄは、分布蓄積部１３０Ｅが蓄積している第２走行状態分布のうち、１周期前に更新した第２走行状態分布を現在の第２走行状態分布として設定する。
分布蓄積部１３０Ｅは、第１走行状態分布更新部１３０Ｃおよび第２走行状態分布更新部１３０Ｄが更新した第１走行状態分布および第２走行状態分布を取得する。そして、分布蓄積部１３０Ｅは、取得した第１走行状態分布および第２走行状態分布を蓄積する。

分布選択部１３０Ｆは、運転状況判定部１２０が運転状態の推定（運転の不安定度の判定）に対して外乱となる運転状況（第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、第４外乱運転状況）にあると判定すると、第１走行状態分布更新部１３０Ｃおよび第２走行状態分布更新部１３０Ｄが更新する走行状態分布（第１走行状態分布、第２走行状態分布）を変更する。具体的には、分布選択部１３０Ｆは、運転状況判定部１２０の判定結果に基づいて、車両の運転状況が第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、および第４外乱運転状況のいずれに該当するのかを判定する。そして、分布選択部１３０Ｆは、第２外乱運転状況に該当すると判定した場合には、第２走行状態分布更新部１３０Ｄが更新した第２走行状態分布を、第１走行状態分布更新部１３０Ｃが更新した第１走行状態分布で置き換える（以下、「リセット処理」とも呼ぶ）。

一方、分布選択部１３０Ｆは、第３外乱運転状況に該当すると判定した場合には、第１走行状態分布更新部１３０Ｃおよび第２走行状態分布更新部１３０Ｄが更新した第１走行状態分布および第２走行状態分布を、第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、および第４外乱運転状況のいずれかにあると判定しているときの分布データ（操舵角予測誤差θe）を含まない第１走行状態分布および第２走行状態分布で置き換える（以下、「リストア処理」とも呼ぶ）。具体的には、分布選択部１３０Ｆは、分布蓄積部１３０Ｅが記憶している第１走行状態分布および第２走行状態分布のうち、少なくとも２周期以上前に更新した第１走行状態分布および第２走行状態分布で置き換える。一方、分布選択部１３０Ｆは、第１外乱運転状況および第４外乱運転状況のいずれかに該当すると判定した場合には、リセット処理およびリストア処理を行わず、第１走行状態分布および第２走行状態分布の置き換えを行わない。

運転不安定度判定部１４０は、走行状態分布更新部１３０が更新した第１走行状態分布および第２走行状態分布（置き換えた場合には、置き換え後の第１走行状態分布、第２走行状態分布）に基づいて運転者の運転状態（運転の不安定度）を推定する。
情報呈示部１５０は、運転不安定度判定部１４０が推定した運転者の運転状態（運転の不安定度）に基づいて運転者に呈示情報を呈示する処理（以下、「情報呈示処理」とも呼ぶ）を行う。情報呈示処理では、情報呈示部１５０は、呈示情報、つまり、運転者に呈示する警報その他の情報を運転者に呈示させる制御信号を情報呈示装置に出力する。

（運転不安定度判定処理）
次に、運転支援部１００Ａが実行する運転不安定度判定処理について説明する。運転不安定度判定処理は、予め定めた制御周期毎に実施する。
図６は、運転不安定度判定処理を表すフローチャートである。
図６に示すように、まず、ステップＳ１０１では、運転支援部１００Ａ（走行状態データ取得部１１０、運転状況判定部１２０）は、車両情報を取得する。車両情報としては、例えば、走行状態データ（操舵角情報）、および運転操作子の操作状態の情報がある。
続いてステップＳ１０２に移行して、運転支援部１００Ａ（運転状況判定部１２０）は、交通環境情報を取得する。交通環境情報としては、例えば、走行環境の情報がある。

続いてステップＳ１０３に移行して、運転支援部１００Ａ（運転状況判定部１２０）は、ステップＳ１０１で取得した車両情報、およびステップＳ１０２で取得した交通環境情報に基づいて、車両の運転状況（第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、第４外乱運転状況、および通常運転状況）を判定する。具体的には、運転支援部１００Ａ（運転状況判定部１２０）は、ステップＳ１０１で取得した車両情報、およびステップＳ１０２で取得した交通環境情報に基づいて、運転操作子の操作状態、走行環境、または車両状態を検出する。続いて、運転支援部１００Ａ（運転状況判定部１２０）は、検出した運転操作子の操作状態、走行環境、または車両状態に基づいて、車両の運転状況が第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、第４外乱運転状況、および通常運転状況のいずれの状況にあるか判定する。

続いてステップＳ１０４に移行して、走行状態分布更新部１３０（分布データ算出部１３０Ａ）は、ステップＳ１０１で取得した走行状態データ（操舵角情報）に基づいて操舵角予測誤差θeを算出する。ここで、図７に、絶対エントロピーＨp1、Ｈp2、相対エントロピーＲＨpを算出するために用いる特殊記号および当該特殊記号の名称を示す。操舵角円滑値θn-tildeは、量子化ノイズの影響を低減した操舵角θである。また、操舵角の推定値θn-hatは、ステアリングホイールを滑らかに操作したと仮定してサンプリング時点における操舵角θを推定した値である。操舵角推定値θn-hatは、以下の（式１）に示すように、操舵角円滑値θn-tildeに対して二次のテイラー展開を施して得る。

（式１）において、ｔnは操舵角θnのサンプリング時刻である。
操舵角円滑値θn-tildeは、量子化ノイズの影響を低減するために、３個の隣接した操舵角θnの平均値として以下の（式２）から算出する。

（式２）において、ｌは、操舵角円滑値θn-tildeの算出時間間隔を１５０ミリ秒、すなわち、手動操作において人間が断続的に操作可能な最小時間間隔とした場合に、１５０ミリ秒内に含まれる操舵角θnのサンプル数を表す。
操舵角θnのサンプリング間隔をＴsとすると、サンプル数ｌは、以下の（式３）で表する。
ｌ＝round（０．１５／Ｔs） ・・・（式３）
（式３）において、ｋ＝１、２、３の値をとり、（ｋ＊１）により１５０ミリ秒間隔の操舵角とそれに隣接する合計３個の操舵角θnに基づいて、円滑値θn-tildeを求めることができる。したがって、このような円滑値θn-tildeに基づいて算出する推定値θn-hatは、実質的に１５０ミリ秒間隔で得た操舵角θにより算出したことになる。

図８は、操舵角予測誤差θeを説明するための図である。
サンプリング時点における操舵角予測誤差θeは、図８に示すように、ステアリングホイールを滑らかに操作したと仮定した場合のサンプリング時点における操舵角推定値θn-hatと実際の操舵角θnとの差として、以下の（式４）から算出できる。

ただし、操舵角予測誤差θeは、手動操作において人間が断続的に操作可能な最小時間間隔、すなわち、１５０ミリ秒毎の操舵角θnに対してのみ算出するものとする。
以下に、操舵角予測誤差θeの具体的な算出方法を説明する。なお、操舵角θのサンプリング間隔Ｔsは、例えば、５０ミリ秒とする。まず、１５０ミリ秒間隔の隣接する３個の操舵角θnを用いて、上記（式２）から３個の操舵角円滑値θn-tildeを算出する。３個の操舵角円滑値θn-tildeは、以下の（式５）で表す。

次に、算出した３個の操舵角円滑値θn-tildeを用いて、上記（式１）から操舵角の推定値θn-hatを算出する。推定値θn-hatは、以下の（式６）で表す。

そして、算出した操舵角推定値θn-hatと実際の操舵角θnとを用いて、上記（式４）から操舵角予測誤差θeを算出する。

図９は、予測誤差区分ｂｉの範囲を説明するための図である。
続いてステップＳ１０５に移行して、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）は、ステップＳ１０４で算出した操舵角予測誤差θeをデータ一時蓄積部１３０Ｂに蓄積する。具体的には、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）は、ステップＳ１０４で算出した操舵角予測誤差θeを９個の予測誤差区分ｂｉ（＝ｂ１、ｂ２、ｂ３、ｂ４、ｂ５、ｂ６、ｂ７、ｂ８、ｂ９）に分類する。

予測誤差区分ｂｉの範囲は、図９に示すように、ステアリングエントロピーの算出に用いるα値に基づいて設定する。α値としては、例えば、操舵角θの時系列データに基づいて一定時間内の操舵角予測誤差θe、すなわち、ステアリングホイールを滑らかに操作したと仮定した場合の操舵角推定値θn-hatと実際の操舵角θnとの差を求め、操舵角予測誤差θeの分布（ばらつき）を測定して９０パーセントタイル値（操舵角予測誤差θeの９０％が含まれる分布の範囲）を算出したものがある。すなわち、α値は、操舵角予測誤差θeの９０％が区間[−α、α]の中に含まれるように設定する。

具体的には、予測誤差区分ｂ１は５α未満とし、予測誤差区分ｂ２は−５α以上で且つ−２．５α未満とし、予測誤差区分ｂ３は−２．５α以上で且つ−α未満とし、予測誤差区分ｂ４は−α以上で且つ−０．５α未満とし、予測誤差区分ｂ５は−０．５α以上で且つ０．５α未満とする。また、予測誤差区分ｂ６は０．５α以上で且つα未満とし、予測誤差区分ｂ７はα以上で且つ２．５α未満とし、予測誤差区分ｂ８は２．５α以上で且つ５α未満とし、予測誤差区分ｂ９は５α以上とする。

続いて、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）は、分類結果に基づいて区分毎標本数Ｎi（＝Ｎ１、Ｎ２、Ｎ３、Ｎ４、Ｎ５、Ｎ６、Ｎ７、Ｎ８、Ｎ９）を設定する。具体的には、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）は、操舵角予測誤差θeを予測誤差区分ｂ１に分類した場合には区分毎標本数Ｎ１に１を加算し、以下同様に、操舵角予測誤差θeを予測誤差区分ｂｊ（ｊは２〜９のいずれか）に分類した場合には区分毎標本数Ｎｊに１を加算する。区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９）の初期値は０とする。続いて、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）は、区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９）をデータ一時蓄積部１３０Ｂに蓄積する。これにより、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）は、ステップＳ１０４で算出した操舵角予測誤差θeを９個の予測誤差区分ｂｉ（＝ｂ１〜ｂ９）に分類し、分類した操舵角予測誤差θeを予測誤差区分ｂｉ（＝ｂ１〜ｂ９）毎の度数として蓄積する。

なお、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）は、この運転不安定度判定処理の開始時にデータ一時蓄積部１３０Ｂが蓄積しているデータをすべて破棄し、データ一時蓄積部１３０Ｂが蓄積している区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９）を初期化する。また、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）は、ステップＳ１０３で第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、および第４外乱運転状況のいずれかの運転状況にあると判定すると、ステップＳ１０４で算出した操舵角予測誤差θeのデータ一時蓄積部１３０Ｂへの蓄積を中断する。すなわち、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）は、操舵角予測誤差θeの分類や、区分毎標本数Ｎi（＝Ｎ１〜Ｎ９）の算出・蓄積を中断する。これにより、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）は、ステップＳ１０３で第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、および第４外乱運転状況にあると判定していない期間、つまり、通常運転状況にあると判定している期間の分布データ（操舵角予測誤差θe）のみを蓄積する。

また、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）は、ステップＳ１０３で第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、および第４外乱運転状況のいずれかの運転状況から通常運転状況に変化したと判定すると、ステップＳ１０４で算出した操舵角予測誤差θeのデータ一時蓄積部１３０Ｂへの蓄積を再開する。すなわち、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）は、操舵角予測誤差θeの分類や、区分毎標本数Ｎi（＝Ｎ１〜Ｎ９）の算出・蓄積を再開する。その際、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）は、データ一時蓄積部１３０Ｂが蓄積しているデータをすべて破棄し、データ一時蓄積部１３０Ｂが蓄積している区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９）を初期化する。

続いてステップＳ１０６に移行して、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ、第１走行状態分布更新部１３０Ｃ、第２走行状態分布更新部１３０Ｄ）は、この運転不安定度判定処理を開始してから予め定めた設定時間（例えば、５秒）経過したか否かを判定する。具体的には、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ、第１走行状態分布更新部１３０Ｃ、第２走行状態分布更新部１３０Ｄ）は、この運転不安定度判定処理を開始してからの経過時間をカウントするタイマー（不図示）のタイマー値Ｔが設定時間（例えば、５秒）以上であるか否かを判定する。そして、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ、第１走行状態分布更新部１３０Ｃ、第２走行状態分布更新部１３０Ｄ）は、タイマー値Ｔが設定時間（例えば、５秒）以上であると判定した場合には（Ｙｅｓ）、ステップＳ１０７に移行する。

一方、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ、第１走行状態分布更新部１３０Ｃ、第２走行状態分布更新部１３０Ｄ）は、タイマー値Ｔが設定時間（例えば、５秒）未満であると判定した場合には（Ｎｏ）、ステップＳ１０１に移行する。これにより、運転支援部１００Ａは、この運転不安定度判定処理を開始してから設定時間（例えば、５秒）が経過するまでは、ステップＳ１０１〜Ｓ１０６のフローを繰り返し実行する。すなわち、運転支援部１００Ａは、設定時間（例えば、５秒）が経過するまではステップＳ１０７以降のフローは実行しない。したがって、運転支援部１００Ａは、後述するステップＳ１０７〜Ｓ１１７、つまり、第１走行状態分布、第２走行状態分布の更新等を第２周期（例えば、５秒）毎に実行する。
ここで、本実施形態では、運転支援部１００Ａは、ステップＳ１０１〜Ｓ１０６のフローを予め定めた第１周期（例えば、５０ミリ秒毎）で実行する。

図１０は、第２周期の開始点を説明するための図である。
なお、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ、第１走行状態分布更新部１３０Ｃ、第２走行状態分布更新部１３０Ｄ）は、図１０に示すように、ステップＳ１０３で第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、および第４外乱運転状況のいずれかの運転状況から通常運転状況に変化したと判定するとタイマー値Ｔを０とする。これにより、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ、第１走行状態分布更新部１３０Ｃ、第２走行状態分布更新部１３０Ｄ）は、第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、または第４外乱運転状況のいずれかの運転状況から通常運転状況に変化したときに第１走行状態分布更新部１３０Ｃ、第２走行状態分布更新部１３０Ｄが更新する第１走行状態分布、第２走行状態分布の更新の中断を解除しその時点を第２周期（例えば、５秒）の開始点に変更する（以下、「レジューム処理ともよぶ）。

また、本実施形態では、タイマー値Ｔが設定時間（例えば、５秒）以上であるか否かを判定する例を示したが、他の構成を採用することもできる。例えば、データ一時蓄積部１３０Ｂが予め定めた設定数の分布データ（操舵角予測誤差θe）を蓄積したか否かを判定する。具体的には、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ、第１走行状態分布更新部１３０Ｃ、第２走行状態分布更新部１３０Ｄ）が、データ一時蓄積部１３０Ｂが第２周期（例えば、５秒）分の分布データ（操舵角予測誤差θe）を蓄積したか否かを判定する。第２周期（例えば、５秒）分の分布データ（操舵角予測誤差θe）としては、例えば、１００（＝５秒／５０ミリ秒）個の分布データ（操舵角予測誤差θe）がある。

そして、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ、第１走行状態分布更新部１３０Ｃ、第２走行状態分布更新部１３０Ｄ）が、データ一時蓄積部１３０Ｂが第２周期（例えば、５秒）分の分布データ（操舵角予測誤差θe）を蓄積したと判定した場合には（Ｙｅｓ）、ステップＳ１０７に移行する。一方、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ、第１走行状態分布更新部１３０Ｃ、第２走行状態分布更新部１３０Ｄ）が、データ一時蓄積部１３０Ｂが第２周期（例えば、５秒）分の分布データ（操舵角予測誤差θe）を蓄積していないと判定した場合には（Ｎｏ）、ステップＳ１０１に移行する。

図１１、図１２は、第１走行状態分布および第２走行状態分布の更新方法を説明するための図である。
ステップＳ１０７では、運転支援部１００Ａ（第１走行状態分布更新部１３０Ｃ）は、データ一時蓄積部１３０Ｂが蓄積している操舵角予測誤差θe（区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９））に基づいて第１走行状態分布を更新する。具体的には、運転支援部１００Ａ（第１走行状態分布更新部１３０Ｃ）は、図１１に示すように、データ一時蓄積部１３０Ｂが蓄積している区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９）に基づき、以下の（式７）から各予測誤差区分ｂｉに含まれる操舵角予測誤差θeの度数の全度数に対する確率ｐi（＝ｐ１、ｐ２、ｐ３、ｐ４、ｐ５、ｐ６、ｐ７、ｐ８、ｐ９）を求める。そして、運転支援部１００Ａ（第１走行状態分布更新部１３０Ｃ）は、更新した第１走行状態分布を分布蓄積部１３０Ｅに蓄積する。

（式７）において、ｐiは、６４ビット実数（double型）を用いて算出する。また、ｐioldはこの運転不安定度判定処理の前回の実行時に算出したｐiであり、Ｋwindowは第１走行状態分布の更新に用いる標本数（操舵角予測誤差θeの数）、Ｎallは区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９）の合計値である。本実施形態では、（式７）のＫwindowは、４３２０００（＝２１６０秒／５０ミリ秒/回）になる。また、Ｎallは、ステップＳ１０３で運転状況が通常運転状況にあると判定すると、１００（＝５秒／５０ミリ秒/回）になる。なお、Ｎall、つまり、区分毎標本数Ｎ１〜Ｎ９の合計値は、ステップＳ１０３で運転状況が第１外乱運転状況、第２外乱運転状況、第３外乱運転状況および第４外乱運転状況のいずれかであると判定すると１００より小さい数になる。これにより、運転支援部１００Ａ（第１走行状態分布更新部１３０Ｃ）は、ステップＳ１０３で運転状況が第４外乱運転状況にあると判定すると、ステップＳ１０３で運転状況が通常運転状況にあると判定している期間の操舵角予測誤差θe（区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９））のみに基づいて第１走行状態分布を更新する（更新周期変更処理）。
なお、運転支援部１００Ａ（第１走行状態分布更新部１３０Ｃ）は、ステップＳ１０３で運転状況が第１外乱運転状況にあると判定すると、第１走行状態分布の更新を中断する。すなわち、運転支援部１００Ａ（第１走行状態分布更新部１３０Ｃ）は、分布蓄積部１３０Ｅが蓄積している第１走行状態分布のうち、１周期前に更新した第１走行状態分布を現在の第１走行状態分布とする（更新中断処理）。

続いてステップＳ１０８に移行して、運転支援部１００Ａ（第２走行状態分布更新部１３０Ｄ）は、データ一時蓄積部１３０Ｂが蓄積している操舵角予測誤差θe（区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９））に基づいて第２走行状態分布を更新する。具体的には、運転支援部１００Ａ（第２走行状態分布更新部１３０Ｄ）は、データ一時蓄積部１３０Ｂが蓄積している区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９）に基づき、以下の（式８）から各予測誤差区分ｂｉに含まれる操舵角予測誤差θeの度数の全度数に対する確率ｑi（＝ｑ１、ｑ２、ｑ３、ｑ４、ｑ５、ｑ６、ｑ７、ｑ８、ｑ９）を求める。そして、運転支援部１００Ａ（第２走行状態分布更新部１３０Ｄ）は、算出した第２走行状態分布を分布蓄積部１３０Ｅに蓄積する。

（式８）において、ｑiは、６４ビット実数（double型）を用いて算出する。また、ｑioldはこの運転不安定度判定処理の前回の実行時に算出したｑiであり、Ｋwindowは第２走行状態分布の更新に用いる標本数（操舵角予測誤差θeの数）、Ｎallは区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９）の合計値である。本実施形態では、（式８）のＫwindowは、１８０００（＝９０秒／５０ミリ秒/回）になる。また、Ｎallは、ステップＳ１０３で運転状況が通常運転状況にあると判定すると、１００（＝５秒／５０ミリ秒/回）になる。なお、Ｎall、つまり、区分毎標本数Ｎ１〜Ｎ９の合計値は、ステップＳ１０３で運転状況が第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、および第４外乱運転状況のいずれかであると判定すると１００より小さい数になる。これにより、運転支援部１００Ａ（第１走行状態分布更新部１３０Ｃ）は、ステップＳ１０３で運転状況が第４外乱運転状況にあると判定すると、運転状況が通常運転状況にあると判定している期間の操舵角予測誤差θe（区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９））のみに基づいて第２走行状態分布を更新する（更新周期変更処理）。

これにより、第２走行状態分布は、図１３に示すように、運転が不安定状態にある場合（例えば、ふらつき時）には、第１走行状態分布（通常時）に比べ、中央部付近（予測誤差区分ｂ５付近）が低減し、端部付近（予測誤差区分ｂ１、ｂ９付近）が増大する。
なお、運転支援部１００Ａ（第２走行状態分布更新部１３０Ｄ）は、ステップＳ１０３で運転状況が第１外乱運転状況にあると判定すると、第２走行状態分布の更新を中断する。すなわち、運転支援部１００Ａ（第２走行状態分布更新部１３０Ｄ）は、分布蓄積部１３０Ｅが蓄積している第２走行状態分布のうち、１周期前に更新した第２走行状態分布を現在の第２走行状態分布として設定する（更新中断処理）。

続いてステップＳ１０９に移行して、運転支援部１００Ａ（分布選択部１３０Ｆ）は、ステップＳ１０３の判定結果に基づいて、車両の運転状況が第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、第４外乱運転状況、および通常運転状況のいずれに該当するのかを判定する。そして、運転支援部１００Ａ（分布選択部１３０Ｆ）は、第２外乱運転状況に該当すると判定した場合には、ステップＳ１１１に移行する。一方、運転支援部１００Ａ（分布選択部１３０Ｆ）は、第３外乱運転状況に該当すると判定した場合には、ステップＳ１１２に移行する。一方、運転支援部１００Ａ（分布選択部１３０Ｆ）は、第１外乱運転状況、第４外乱運転状況、および通常運転状況のいずれかに該当すると判定した場合には、ステップＳ１１０に移行する。

ステップＳ１１０では、運転支援部１００Ａ（分布選択部１３０Ｆ）は、リセット処理、リストア処理を行わず、ステップＳ１０７、Ｓ１０８で更新した第１走行状態分布および第２走行状態分布の置き換えを行わずに、ステップＳ１１３に移行する。
一方、ステップＳ１１１では、運転支援部１００Ａ（分布選択部１３０Ｆ）は、リセット処理を行った後、ステップＳ１１３に移行する。リセット処理では、運転支援部１００Ａ（分布選択部１３０Ｆ）は、ステップＳ１０８で更新した第２走行状態分布を、ステップＳ１０７で更新した第１走行状態分布で置き換える。

一方、ステップＳ１１２では、運転支援部１００Ａ（分布選択部１３０Ｆ）は、リストア処理を行った後、ステップＳ１１３に移行する。リストア処理では、ステップＳ１０７、Ｓ１０８で更新した第１走行状態分布および第２走行状態分布を、第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、および第４外乱運転状況のいずれかにあると判定しているときの分布データ（操舵角予測誤差θe）を含まない第１走行状態分布および第２走行状態分布で置き換える。
具体的には、運転支援部１００Ａ（分布選択部１３０Ｆ）は、分布蓄積部１３０Ｅが記憶している第１走行状態分布および第２走行状態分布のうち、少なくとも２周期以上前に更新した第１走行状態分布および第２走行状態分布で置き換える。例えば、運転支援部１００Ａ（分布選択部１３０Ｆ）は、運転操作子の操作状態、走行環境、および車両状態に基づいて、変速操作およびうねりのいずれかが発生しているか否かを判定する。

そして、運転支援部１００Ａ（分布選択部１３０Ｆ）は、変速操作およびうねりのいずれかが発生していると判定した場合には、ステップＳ１０７で更新した第１走行状態分布を２周期前に更新した第１走行状態分布で置き換えるとともに、ステップＳ１０８で更新した第２走行状態分布を２周期前に更新した第２走行状態分布で置き換える。
一方、運転支援部１００Ａ（分布選択部１３０Ｆ）は、変速操作およびうねりのいずれも発生していないと判定した場合には、ステップＳ１０７で更新した第１走行状態分布を３周期前に更新した第１走行状態分布とするとともに、ステップＳ１０８で更新した第２走行状態分布を３周期前に更新した第２走行状態分布で置き換える。

続いてステップＳ１１３に移行して、運転支援部１００Ａ（運転不安定度判定部１４０）は、ステップＳ１０７、Ｓ１０８で更新した第１走行状態分布および第２走行状態分布（置き換えた場合には、置き換え後の第１走行状態分布、第２走行状態分布）に基づいて絶対エントロピーＨp1、Ｈp2を算出する。具体的には、走行状態分布更新部１３０（運転不安定度判定部１４０）は、ステップＳ１０７で更新した第１走行状態分布に基づき、以下の（式９）から絶対エントロピーＨp1を算出する。また、走行状態分布更新部１３０（運転不安定度判定部１４０）は、ステップＳ１０８で更新した第２走行状態分布に基づき、以下の（式１０）から絶対エントロピーＨp2を算出する。

上記（式９）（式１０）より、絶対エントロピーＨp1、Ｈp2が小さいほど、第１走行状態分布および第２走行状態分（操舵角予測誤差θeの分布）の峻険度が増大し、操舵角予測誤差θeの分布が一定の範囲に収まる。すなわち、運転操作が滑らかに行われ、運転が安定状態にあることを示す。一方、絶対エントロピーＨp1、Ｈp2が大きいほど、操舵角予測誤差θeの分布の峻険度が低減し、操舵角予測誤差θeの分布がばらつく。すなわち、運転操作が滑らかに行われず、運転が不安定状態にあることを示す。

続いてステップＳ１１４に移行して、運転支援部１００Ａ（運転不安定度判定部１４０）は、Ｒ³-ステアリングエントロピー法によって、ステップＳ１０７、Ｓ１０８で更新した第１走行状態分布および第２走行状態分布（置き換えた場合には、置き換え後の第１走行状態分布、第２走行状態分布）の分布間の相違量（相対エントロピーＲＨp）を算出する。具体的には、運転支援部１００Ａ（運転不安定度判定部１４０）は、ステップＳ１０７で算出した確率ｐiおよびステップＳ１０８で算出した確率ｑiに基づき、以下の（式１１）から相対エントロピーＲＨpを算出する。

図１４は、相対エントロピーＲＨpを説明するための図である。
上記（式１１）より、相対エントロピーＲＨpが小さいほど、第１走行状態分布の確率ｐi（＝ｐ１〜ｐ９）と第２走行状態分布の確率ｑi（＝ｑ１〜ｑ９）とのずれが小さくなる。すなわち、図１４に示すように、運転者の現在の運転操作が普段の運転操作と同様に滑らかに行われ、運転が安定状態にあることを示す。一方、相対エントロピーＲＨpが大きいほど、第１走行状態分布の確率ｐi（＝ｐ１〜ｐ９）と第２走行状態分布の確率ｑi（＝ｑ１〜ｑ９）とのずれが大きくなる。すなわち、運転者の現在の運転操作が普段の運転操作と比べて滑らかに行われず、運転が不安定状態にあることを示す。

続いてステップＳ１１５に移行して、運転支援部１００Ａ（運転不安定度判定部１４０）は、ステップＳ１１３で算出した絶対エントロピーＨp1、Ｈp2に基づいて運転者の運転状態を推定する（運転が不安定状態にあるか否かを判定する）。具体的には、運転支援部１００Ａ（運転不安定度判定部１４０）は、ステップＳ１１３で算出した絶対エントロピーＨp1、Ｈp2の差（Ｈp2−Ｈp1）が予め定めた判定閾値よりも大きいか否かを判定する。そして、運転支援部１００Ａ（運転不安定度判定部１４０）は、絶対エントロピーＨp1、Ｈp2の差（Ｈp2−Ｈp1）が判定閾値よりも大きいと判定した場合には、運転が不安定状態にあると判定する。一方、運転支援部１００Ａ（運転不安定度判定部１４０）は、絶対エントロピーＨp1、Ｈp2の差（Ｈp2−Ｈp1）が判定閾値以下であると判定した場合には、運転が安定状態にあると判定する。

続いてステップＳ１１６に移行して、運転支援部１００Ａ（運転不安定度判定部１４０）は、ステップＳ１１４で算出した相対エントロピーＲＨpに基づいて運転者の運転状態を推定する（運転が不安定状態にあるか否かを判定する）。具体的には、運転支援部１００Ａ（運転不安定度判定部１４０）は、ステップＳ１１４で算出した相対エントロピーＲＨpが予め定めた判定閾値よりも大きいか否かを判定する。そして、運転支援部１００Ａ（運転不安定度判定部１４０）は、相対エントロピーＲＨpが判定閾値よりも大きいと判定した場合には、運転が不安定状態にあると判定する。一方、運転支援部１００Ａ（運転不安定度判定部１４０）は、相対エントロピーＲＨpが判定閾値以下であると判定した場合には、運転が安定状態にあると判定する。

続いてステップＳ１１７に移行して、運転支援部１００Ａ（情報呈示部１５０）は、ステップＳ１１５、Ｓ１１６で推定（判定）した運転状態に基づいて運転者に呈示情報を呈示する処理（情報呈示処理）を行う。具体的には、運転支援部１００Ａ（情報呈示部１５０）は、ステップＳ１１５およびＳ１１６の両方で不安定状態と判定した状態が予め定めた設定時間（例えば、５秒）以上継続したか否かを判定する。そして、運転支援部１００Ａ（情報呈示部１５０）は、ステップＳ１１５およびＳ１１６の両方で不安定状態と判定した状態が設定時間（例えば、５秒）以上継続したと判定した場合には、情報呈示処理を行う。一方、運転支援部１００Ａ（情報呈示部１５０）は、ステップＳ１１５およびＳ１１６の両方で不安定状態と判定した状態が設定時間（例えば、５秒）以上継続していないと判定した場合には、情報呈示処理を行わない。

情報呈示処理の例を、図１５に示す。この例では、警告表示を行うとともに、「ピー！！そろそろ休憩しませんか。」等と音声で警告の呈示を行う。
なお、本実施形態では、不安定状態と判定した状態が設定時間（例えば、５秒）以上継続したと判定した場合に、呈示情報を呈示する処理（情報呈示処理）を行う例を示したが、他の構成を採用することもできる。例えば、図１６に示すように、レベルゲージによって複数段階の呈示レベルで呈示するとともに、各呈示レベルに対応した聴覚情報を呈示する構成としてもよい。この場合、運転状態として運転の不安定度が高くなるほど、レベルゲージの呈示レベルを高くする。図１６では、８段階の呈示レベルとし、左側の表示ほど呈示レベルが高い状態（運転の不安定度が高い状態）を示している。

（動作その他）
次に、本実施形態の運転状態推定装置を搭載した車両の動作について説明する。
車両の走行中、運転支援部１００Ａが、運転不安定度判定処理を実行したとする。すると、運転支援部１００Ａ（走行状態データ取得部１１０、運転状況判定部１２０）が、走行状態データ（操舵角情報）、車両情報および交通環境情報を取得する（図６のステップＳ１０１、Ｓ１０２）。続いて、運転支援部１００Ａ（運転状況判定部１２０）が、取得した車両情報および交通環境情報に基づき、車両の運転状況（第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、第４外乱運転状況、通常運転状況）を判定する（図６のステップＳ１０３）。続いて、走行状態分布更新部１３０（分布データ算出部１３０Ａ）が、取得した走行状態データ（操舵角情）に基づいて分布データ（操舵角予測誤差θe）を算出する（ステップＳ１０４）。

ここで、運転支援部１００Ａ（運転状況判定部１２０）が、ステップＳ１０３で運転状況が通常運転状況にあると判定したとする。すると、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）が、算出した分布データ（操舵角予測誤差θe）を９個の予測誤差区分ｂｉ（＝ｂ１〜ｂ９）に分類し、分類結果に基づいて区分毎標本数Ｎi（＝Ｎ１〜Ｎ９）を設定する（図６のステップＳ１０５）。続いて、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）が、設定した区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９。分布データ（操舵角予測誤差θe））をデータ一時蓄積部１３０Ｂに蓄積する（図６のステップＳ１０５）。
続いて、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ、第１走行状態分布更新部１３０Ｃ、第２走行状態分布更新部１３０Ｄ）が、この運転不安定度判定処理を開始してから設定時間（例えば、５秒）経過していないと判定する（図６のステップＳ１０６「Ｎｏ」）。そして、運転支援部１００Ａが、ステップＳ１０１〜Ｓ１０６のフローを繰り返し、第１周期（例えば、５０ミリ秒）毎に区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９。分布データ（操舵角予測誤差θe））を算出・蓄積する。

また、上記ステップＳ１０１〜Ｓ１０６のフローを繰り返すうちに、分布データ（操舵角予測誤差θe）の算出・蓄積を１００回実行し、この運転不安定度判定処理を開始してから５秒経過したとする。すると、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ、第１走行状態分布更新部１３０Ｃ、第２走行状態分布更新部１３０Ｄ）が、この運転不安定度判定処理を開始してから設定時間（例えば、５秒）経過したと判定する（図６のステップＳ１０６「Ｙｅｓ」）。続いて、運転支援部１００Ａ（第１走行状態分布更新部１３０Ｃ）が、データ一時蓄積部１３０Ｂが蓄積している操舵角予測誤差θe（区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９））に基づいて第１走行状態分布および第２走行状態分布を更新する（図６のステップＳ１０７、Ｓ１０８）。第１走行状態分布および第２走行状態分布の更新には、６４ビット実数（double型）、つまり、精度が高く比較的演算負荷の高い変数を用いる。続いて、運転支援部１００Ａ（分布選択部１３０Ｆ）が、ステップＳ１０３の判定結果に基づき、車両の運転状況が通常運転状況に該当すると判定する（図６のステップＳ１０９、Ｓ１１０）。続いて、運転支援部１００Ａ（運転不安定度判定部１４０）が、更新した第１走行状態分布および第２走行状態分布に基づいて絶対エントロピーＨp1、Ｈp2を算出する（図６のステップＳ１１３）。絶対エントロピーＨp1、Ｈp2の算出では、（式９）（式１０）に示すように、比較的演算負荷の高いlog計算を行う。

続いて、運転支援部１００Ａ（運転不安定度判定部１４０）が、Ｒ3-ステアリングエントロピー法によって、更新した第１走行状態分布および第２走行状態分布の分布間の相違量（相対エントロピーＲＨp）を算出する（図６のステップＳ１１４）。相対エントロピーＲＨpの算出では、（式１１）に示すように、比較的演算負荷の高いlog計算を行う。続いて、運転支援部１００Ａ（運転不安定度判定部１４０）が、算出した絶対エントロピーＨp1、Ｈp2に基づき運転者の運転状態を推定する（運転者の運転状態が不安定状態にあるか否かを判定する）（図６のステップＳ１１５）。続いて、運転支援部１００Ａ（運転不安定度判定部１４０）が、算出した相対エントロピーＲＨpに基づいて運転者の運転状態を推定する（運転者の運転状態が不安定状態にあるか否かを判定する）（図６のステップＳ１１６）。続いて、運転支援部１００Ａ（情報呈示部１５０）が、推定した運転状態に基づいて運転者に呈示情報を呈示する処理を行う（図６のステップＳ１１７）。これにより、運転支援部１００Ａが、第２周期（例えば、５秒）毎に第１走行状態分布、第２走行状態分布、絶対エントロピーＨp1、Ｈp2、相対エントロピーＲＨpを算出する。

このように、本実施形態では、運転支援部１００Ａは、操舵角情報に基づいて第１走行状態分布および第２走行状態分布用の分布データ（操舵角予測誤差θe）を予め定めた第１周期（例えば、５０ミリ秒）毎に取得する。そして、運転支援部１００Ａは、取得した分布データ（操舵角予測誤差θe）に基づいて第１周期（例えば、５０ミリ秒）よりも長い第２周期（例えば、５秒）毎に第１走行状態分布および第２走行状態分布を更新する。それゆえ、第１走行状態分布および第２走行状態分布の更新（６４ビット実数による演算）や、絶対エントロピーＨp1、Ｈp2、相対エントロピーＲＨpの算出（log計算）に基づく演算の実行頻度、つまり、比較的演算負荷の高い演算の実行頻度を低減できる。これにより、運転状態推定装置の演算負荷を低減できる。それゆえ、スマートフォン、廉価な車載コントローラ等、演算能力が比較的低い機器で運転状態推定装置を実現できる。

図１７は、本実施形態の運転状態推定装置を搭載した車両の動作を表す図である。図１７では、第１周期を５０ミリ秒、第２周期を２．５秒とした場合の動作を表す。図１７において、ＲＨpnewは、分布データ（操舵角予測誤差θe）を取得する周期（第１周期。５０ミリ秒）よりも長い周期（第２周期。２．５秒）で第１走行状態分布等を算出する場合、つまり、本実施形態の相対エントロピーである。また、ＲＨpoldは、分布データ（操舵角予測誤差θe）を取得する周期と第１走行状態分布等を算出する周期とが同じ（５０ミリ秒）である場合（以下、「改良前」とも呼ぶ）の相対エントロピーである。
図１７（ａ）に示すように、本実施形態の相対エントロピーＲＨpnewは、改良前の相対エントロピーＲＨoldと２．５秒毎に重なる。また、図１７（ｂ）に示すように、車両のトリップ中、本実施形態の相対エントロピーＲＨpnewは、改良前の相対エントロピーＲＨoldとほぼ重なっており、実用上問題がないものと考える。

一方、車両が、加速を行ったとする。すると、運転支援部１００Ａ（運転状況判定部１２０）が、取得した車両情報および交通環境情報に基づき、車両の運転状況が第１外乱運転状況、つまり、第１走行状態分布および第２走行状態分布の更新を行わない運転状況にあると判定する（図６のステップＳ１０３）。続いて、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）が、算出した分布データ（操舵角予測誤差θe）のデータ一時蓄積部１３０Ｂへの蓄積を中断する（図６のステップＳ１０５）。すなわち、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）は、操舵角予測誤差θeの分類や、区分毎標本数Ｎi（＝Ｎ１〜Ｎ９）の算出を中断する。

そして、運転支援部１００Ａが、上記ステップＳ１０１〜Ｓ１０６のフローを繰り返すうちに、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ、第１走行状態分布更新部１３０Ｃ、第２走行状態分布更新部１３０Ｄ）が、この運転不安定度判定処理を開始してから設定時間（例えば、５秒）経過したと判定したとする（図６のステップＳ１０６「Ｙｅｓ」）。すると、運転支援部１００Ａ（第１走行状態分布更新部１３０Ｃ、第２走行状態分布更新部１３０Ｄ）が、第１走行状態分布、第２走行状態分布の更新を中断すると判定する（図６のステップＳ１０７、Ｓ１０８）。すなわち、分布蓄積部１３０Ｅが蓄積している第１走行状態分布、第２走行状態分布のうち、１周期前に更新した第１走行状態分布を現在の第１走行状態分布とするとともに、１周期前に更新した第２走行状態分布を現在の第２走行状態分布とする（更新中断処理）。

このように、本実施形態では、運転支援部１００Ａは、第１外乱運転状況にあると判定すると、更新した走行状態データ（操舵角情報）を用いず、１周期前に更新した第１走行状態分布、第２走行状態分布を現在の第１走行状態分布、第２走行状態分布とする（図６のステップＳ１０７、Ｓ１０８）。これにより、本実施形態は、運転の不安定度の判定（運転の不安定状態の判定）において、不安定状態にあると誤検知することを防止できる。
また、車両が、加速を終了したとする。すると、運転支援部１００Ａ（運転状況判定部１２０）が、取得した車両情報および交通環境情報に基づき、車両の運転状況が第１外乱運転状況から通常運転状況に変化したと判定する（図６のステップＳ１０３）。続いて、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）が、データ一時蓄積部１３０Ｂが蓄積しているデータすべてを破棄し、データ一時蓄積部１３０Ｂが蓄積している区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９）を初期化する（図６のステップＳ１０５）。続いて、運転支援部１００Ａは、タイマー値Ｔを０とする（図６のステップＳ１０６）。これにより、運転支援部１００Ａが、車両の運転状況が通常運転状況に変化したときに、第１走行状態分布、第２走行状態分布の更新周期（第２周期）の計測を開始する（レジューム処理）。

一方、車両が、横滑りを行ったとする。すると、運転支援部１００Ａ（運転状況判定部１２０）が、取得した車両情報および交通環境情報に基づき、車両の運転状況が第２外乱運転状況、つまり、リセット処理を行う運転状況にあると判定する（図６のステップＳ１０３）。そして、ステップＳ１０４〜Ｓ１０８を経て、運転支援部１００Ａが、ステップＳ１０３の判定結果に基づき、車両の運転状況が第２外乱運転状況に該当すると判定する（図６のステップＳ１０９）。続いて、運転支援部１００Ａ（分布選択部１３０Ｆ）が、ステップＳ１０８で更新した第２走行状態分布を、ステップＳ１０７で更新した第１走行状態分布で置き換える（リセット処理。図６のステップＳ１１１）。
このように、本実施形態では、運転支援部１００Ａは、第２外乱運転状況にあると判定すると、リセット処理を行い、第２走行状態分布を第１走行状態分布で置き換える（図６のステップＳ１１１）。これにより、本実施形態は、運転の不安定度の判定（運転の不安定状態の判定）において、不安定状態にあると誤検知することを防止できる。

一方、運転者が、変速操作を行ったとする。すると、運転支援部１００Ａ（運転状況判定部１２０）が、取得した車両情報および交通環境情報に基づき、車両の運転状況が第３外乱運転状況、つまり、リストア処理を行う運転状況にあると判定する（図６のステップＳ１０３）。そして、ステップＳ１０４〜Ｓ１０８を経て、運転支援部１００Ａ（分布選択部１３０Ｆ）が、車両の運転状況が第３外乱運転状況に該当すると判定する（図６のステップＳ１０９）。続いて、運転支援部１００Ａ（分布選択部１３０Ｆ）が、分布蓄積部１３０Ｅが蓄積している第１走行状態分布、第２走行状態分布のうち、２周期前に更新した第１走行状態分布を現在の第１走行状態分布とするとともに、２周期前に更新した第２走行状態分布を現在の第２走行状態分布とする（図６のステップＳ１１２）。
このように、本実施形態では、運転支援部１００Ａは、変速操作を基に第３外乱運転状況にある判定すると、２周期前に更新した第１走行状態分布を現在の第１走行状態分布とし、２周期前に更新した第２走行状態分布を現在の第２走行状態分布とする（図６のステップＳ１１２）。これにより、本実施形態は、運転の不安定度の判定（運転の不安定状態の判定）において、不安定状態にあると誤検知することを防止できる。

一方、運転者が、車線変更を行ったとする。すると、運転支援部１００Ａ（運転状況判定部１２０）が、取得した車両情報および交通環境情報に基づき、車両の運転状況が第３外乱運転状況、つまり、リストア処理を行う運転状況にあると判定する（図６のステップＳ１０３）。そして、ステップＳ１０４〜Ｓ１０８を経て、運転支援部１００Ａ（分布選択部１３０Ｆ）が、車両の運転状況が第３外乱運転状況に該当すると判定する（図６のステップＳ１０９）。続いて、運転支援部１００Ａ（分布選択部１３０Ｆ）が、分布蓄積部１３０Ｅが蓄積している第１走行状態分布、第２走行状態分布のうち、３周期前に更新した第１走行状態分布を現在の第１走行状態分布とするとともに、３周期前に更新した第２走行状態分布を現在の第２走行状態分布とする（図６のステップＳ１１２）。
このように、本実施形態では、運転支援部１００Ａは、車線変更を基に第３外乱運転状況にあると判定すると、３周期前に更新した第１走行状態分布を現在の第１走行状態分布とし、３周期前に更新した第２走行状態分布を現在の第２走行状態分布とする（図６のステップＳ１１２）。これにより、本実施形態は、運転の不安定度の判定（運転の不安定状態の判定）において、不安定状態にあると誤検知することを防止できる。

一方、車両が、路面のジョイントを走行したとする。すると、運転支援部１００Ａ（運転状況判定部１２０）が、取得した車両情報および交通環境情報に基づき、車両の運転状況が第４外乱運転状況、つまり、更新周期変更処理を行う運転状況にあると判定する（図６のステップＳ１０３）。続いて、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）が、算出した分布データ（操舵角予測誤差θe）のデータ一時蓄積部１３０Ｂへの蓄積を中断する（図６のステップＳ１０５）。すなわち、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）は、操舵角予測誤差θeの分類や、区分毎標本数Ｎi（＝Ｎ１〜Ｎ９）の算出を中断する。これにより、データ一時蓄積部１３０Ｂは、通常運転状況にあると判定している期間の操舵角予測誤差θe（区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９））のみを蓄積する。そして、運転支援部１００Ａが、ステップＳ１０１〜Ｓ１０６のフローを繰り返す（図６のステップＳ１０６「Ｎｏ」）。

また、上記ステップＳ１０１〜Ｓ１０６のフローを繰り返すうちに、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ、第１走行状態分布更新部１３０Ｃ、第２走行状態分布更新部１３０Ｄ）が、この運転不安定度判定処理を開始してから設定時間（例えば、５秒）経過したと判定したとする（図６のステップＳ１０６「Ｙｅｓ」）。すると、運転支援部１００Ａ（第１走行状態分布更新部１３０Ｃ）が、データ一時蓄積部１３０Ｂが蓄積している操舵角予測誤差θe（区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９））に基づいて第１走行状態分布および第２走行状態分布を更新する（図６のステップＳ１０７、Ｓ１０８）。これにより、運転支援部１００Ａ（データ一時蓄積部１３０Ｂ）は、通常運転状況にあると判定している期間、つまり、第４外乱運転状況にあると判定していない期間の操舵角予測誤差θe（区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９））に基づいて第１走行状態分布および第２走行状態分布を更新する（更新周期変更処理）。

このように、本実施形態では、運転支援部１００Ａは、第４外乱運転状況を判定すると、通常運転状況にあると判定している期間、つまり、第４外乱運転状況にあると判定していない期間の操舵角予測誤差θe（区分毎標本数Ｎｉ（＝Ｎ１〜Ｎ９））に基づいて第１走行状態分布および第２走行状態分布を更新する（図６のステップＳ１０７、Ｓ１０８）。これにより、本実施形態は、運転の不安定度の判定（運転の不安定状態の判定）において、不安定状態にあると誤検知することを防止できる。

本実施形態では、図３の走行状態データ取得部１１０、図６のステップＳ１０１が走行状態データ取得部を構成する。以下同様に、図３の走行状態分布更新部１３０、図５の分布データ算出部１３０Ａ、図６のステップＳ１０４が分布データ取得部を構成する。また、図３の走行状態分布更新部１３０、図５の第１走行状態分布更新部１３０Ｃ、第２走行状態分布更新部１３０Ｄ、図６のステップＳ１０７、Ｓ１０８が走行状態分布更新部を構成する。さらに、図３の運転不安定度判定部１４０、図６のステップＳ１１４、Ｓ１１６が運転状態推定部を構成する。さらに、図３の走行状態分布更新部１３０、図５のデータ一時蓄積部１３０Ｂがデータ一時蓄積部を構成する。また、図３の運転状況判定部１２０、図６のステップＳ１０３が運転状況判定部を構成する。さらに、図３の走行状態分布更新部１３０、図５の分布蓄積部１３０Ｅ、図６のステップＳ１０７、Ｓ１０８が分布蓄積部を構成する。

（本実施形態の効果）
本実施形態は、次のような効果を奏する。
（１）運転支援部１００Ａは、走行状態データ（操舵角情報）に基づき走行状態分布（第１走行状態分布、第２走行状態分布）用の分布データ（操舵角予測誤差θe）を予め定めた第１周期（例えば、５０ミリ秒）毎に取得する。続いて、運転支援部１００Ａは、取得した分布データ（操舵角予測誤差θe）に基づき第１周期（例えば、５０ミリ秒）よりも長い第２周期（例えば、５秒）毎に複数の走行状態分布（第１走行状態分布、第２走行状態分布）を更新する。続いて、運転支援部１００Ａは、更新した複数の走行状態分布（第１走行状態分布、第２走行状態分布）に基づき運転者の運転状態を推定する。
このような構成により、予め定めた第１周期（例えば、５０ミリ秒）毎に分布データ（操舵角予測誤差θe）を取得するとともに、第１周期（例えば、５０ミリ秒）よりも長い第２周期（例えば、５秒）毎に複数の走行状態分布（例えば、第１走行状態分布、第２走行状態分布）を更新する。それゆえ、複数の走行状態分布（第１走行状態分布、第２走行状態分布）の更新や運転状態の推定等、分布データ（操舵角予測誤差θe）に基づく演算の実行頻度を低減できる。これにより、運転状態推定装置の演算負荷を低減できる。

（２）運転支援部１００Ａは、取得した分布データ（操舵角予測誤差θe）を蓄積するデータ一時蓄積部１３０Ｂを備える。そして、運転支援部１００Ａは、データ一時蓄積部１３０Ｂが蓄積している第２周期（例えば、５秒）分の分布データ（操舵角予測誤差θe）で各走行状態分布（第１走行状態分布、第２走行状態分布）を更新する。
このような構成により、第２周期（例えば、５秒）分の分布データ（操舵角予測誤差θe）で各走行状態分布（第１走行状態分布、第２走行状態分布）を更新する。それゆえ、走行状態分布（第１走行状態分布、第２走行状態分布）をより精度よく更新できる。

（３）運転支援部１００Ａは、データ一時蓄積部１３０Ｂが第２周期（例えば、５秒）分の分布データ（例えば、１００個の操舵角予測誤差θe）を蓄積したと判定すると、各走行状態分布（第１走行状態分布、第２走行状態分布）を更新する。
このような構成により、第２周期（例えば、５秒）分の分布データ（操舵角予測誤差θe）に基づいて走行状態分布（第１走行状態分布、第２走行状態分布）を算出できる。
（４）運転支援部１００Ａは、分布データ（操舵角予測誤差θe）を９個の予測誤差区分ｂｉ（＝ｂ１〜ｂ９）に分類し、分類した分布データ（操舵角予測誤差θe）を予測誤差区分ｂｉ（＝ｂ１〜ｂ９）毎の度数としてデータ一時蓄積部１３０Ｂに蓄積する。
このような構成により、９個の予測誤差区分ｂｉ（＝ｂ１〜ｂ９）毎に分布データ（操舵角予測誤差θe）の度数をデータ一時蓄積部１３０Ｂに蓄積する。それゆえ、９個の予測誤差区分ｂｉ（＝ｂ１〜ｂ９）毎の分布データ（操舵角予測誤差θe）の度数に基づいて、各走行状態分布（第１走行状態分布、第２走行状態分布）を更新できる。

（５）運転支援部１００Ａは、分布データとして、ステアリングホイールを滑らかに操作したと仮定した場合の操舵角推定値と実際の操舵角との差（操舵角予測誤差θe）を取得する。
このような構成により、運転者の操舵操作に関する運転状態を推定できる。
（６）運転支援部１００Ａは、運転状態（運転の不安定度）の推定に対して外乱となる運転状況（第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、第４外乱運転状況）にあると判定すると、運転支援部１００Ａが更新する走行状態分布（第１走行状態分布、第２走行状態分布）を変更する。
このような構成により、運転状態の推定に対して外乱となる運転状況（第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、第４外乱運転状況）が発生した場合にも、走行状態分布（第１走行状態分布、第２走行状態分布）をより適切に更新できる。

（７）運転支援部１００Ａは、運転状態（運転の不安定度）の推定に対して外乱となる運転状況（第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、第４外乱運転状況）にあると判定すると、判定した運転状況（第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、第４外乱運転状況）に基づき、更新中断処理、リストア処理、リセット処理、更新周期変更処理およびレジューム処理の少なくともいずれか１つの処理を実行する。
このような構成により、運転状態の推定に対して外乱となる運転状況（第１外乱運転状況、第２外乱運転状況、第３外乱運転状況、第４外乱運転状況）が発生した場合に、更新中断処理、リストア処理、リセット処理、更新周期変更処理およびレジューム処理の少なくともいずれか１つの処理を実行できる。これにより、実行した処理を基に、走行状態分布（第１走行状態分布、第２走行状態分布）をより適切に更新できる。

（８）運転支援部１００Ａは、走行状態分布（第１走行状態分布、第２走行状態分布）を蓄積する分布蓄積部１３０Ｅを備える。そして、運転支援部１００Ａは、リストア処理として、分布蓄積部１３０Ｅが蓄積している走行状態分布（第１走行状態分布、第２走行状態分布）のうち、少なくとも２周期以上前の走行状態分布（第１走行状態分布、第２走行状態分布）を更新後の走行状態分布（第１走行状態分布、第２走行状態分布）とする。
このような構成により、リストア処理を比較的容易に実行できる。
（９）運転支援部１００Ａは、リストア処理として、運転状況（車線変更、車両の右左折、トンネル出入口、うねり）に基づき、少なくとも２周期以上前の走行状態分布（第１走行状態分布、第２走行状態分布）のうちから走行状態分布（第１走行状態分布、第２走行状態分布）を選択し、選択した走行状態分布（第１走行状態分布、第２走行状態分布）を更新後の走行状態分布（第１走行状態分布、第２走行状態分布）とする。
このような構成により、リストア処理をより適切に実行できる。
（１０）運転支援部１００Ａは、運転者が操作可能な運転操作子の操作状態、走行環境、および車両状態の少なくともいずれかに基づき運転状況を判定する。
このような構成により、運転状況をより適切に判定できる。

（変形例１）
なお、本実施形態では、分布データとして、操舵角予測誤差θeを用いる例を示したが、他の構成を採用することもできる。例えば、ヨーレートを用いてもよい。
図１８は、本実施形態の運転状態推定装置を搭載した車両の構成を表す図である。
具体的には、図１８に示すように、車両は、ヨーレートセンサ１１を備える。
ヨーレートセンサ１１は、車両のヨーレートを検出する。続いて、ヨーレートセンサ１１は、検出したヨーレートをコントローラ１００に出力する。
そして、運転支援部１００Ａは、操舵角予測誤差θeの代わりに、ヨーレートセンサ１１が検出したヨーレートを用いて、絶対エントロピーＨp1、Ｈp2および相対エントロピーＲＨpを算出する（図６のステップＳ１０５〜Ｓ１１４）。
（本変形例の効果）
（１）運転支援部１００Ａは、分布データとして、車両のヨーレートを取得する。
このような構成により、運転者の横方向の運転に関する運転状態を推定できる。

（変形例２）
また、例えば、分布データとして、車両の車線内横位置を用いてもよい。
図１９は、本実施形態の運転状態推定装置を搭載した車両の構成を表す図である。
具体的には、図１９に示すように、車両は、前方カメラ１２を備える。
前方カメラ１２は、車両の前方の道路の画像を撮影する。続いて、前方カメラ１２は、撮影した画像をコントローラ１００に出力する。
そして、運転支援部１００Ａは、操舵角予測誤差θeの代わりに、前方カメラ１２が撮影した画像から車両の車線内横位置を算出し、算出した車線内横位置を用いて、絶対エントロピーＨp1、Ｈp2および相対エントロピーＲＨpを算出する。
（本変形例の効果）
（１）運転支援部１００Ａは、分布データとして、車両の車線内横位置を取得する。
このような構成により、運転者の横方向の運転に関する運転状態を推定できる。
以上、本願が優先権を主張する日本国特許出願２０１３−１５０７２６（２０１３年７月１９日出願）の全内容は、参照により本開示の一部をなす。
ここでは、限られた数の実施形態を参照しながら説明したが、権利範囲はそれらに限定されるものではなく、上記の開示に基づく各実施形態の改変は当業者にとって自明なことである。

１１０ 走行状態データ取得部（走行状態データ取得部）
１２０ 運転状況判定部（運転状況判定部）
１３０ 走行状態分布更新部（分布データ取得部、走行状態分布更新部、データ一時蓄積部、分布蓄積部）
１３０Ａ 分布データ算出部（分布データ取得部）
１３０Ｂ データ一時蓄積部（データ一時蓄積部）
１３０Ｃ 第１走行状態分布更新部（走行状態分布更新部）
１３０Ｄ 第２走行状態分布更新部（走行状態分布更新部）
１３０Ｅ 分布蓄積部（分布蓄積部）
１４０ 運転不安定度判定部（運転状態推定部）
ステップＳ１０１（走行状態データ取得部）
ステップＳ１０３（運転状況判定部）
ステップＳ１０４（分布データ取得部）
ステップＳ１０７、Ｓ１０８（走行状態分布更新部、分布蓄積部）
ステップＳ１１４、Ｓ１１６（運転状態推定部）

Claims (6)

1. 運転者が操作可能な運転操作子の操作状態および車両状態の少なくとも一方を含む走行状態データを取得する走行状態データ取得部と、
   前記走行状態データ取得部が取得した走行状態データに基づき走行状態分布用の分布データを予め定めた第１周期毎に取得する分布データ取得部と、
   前記分布データ取得部が取得した分布データを蓄積するデータ一時蓄積部と、
   前記データ一時蓄積部が蓄積している分布データに基づき前記第１周期よりも長い第２周期毎に時間的範囲の異なる複数の走行状態分布を更新する走行状態分布更新部と、
   前記走行状態分布更新部が更新した複数の走行状態分布に基づき前記運転者の運転状態を推定する運転状態推定部と、
   運転状況を判定する運転状況判定部と、
   前記運転状況判定部が運転状態の推定に対して外乱となる運転状況にあると判定すると、前記走行状態分布更新部が更新する走行状態分布を変更する設定変更部と、を備え、
   前記設定変更部は、前記運転状況判定部が運転状態の推定に対して外乱となる運転状況にあると判定すると、下記（ａ）及び（ｂ）の処理の少なくともいずれか１つの処理を実行することを特徴とする運転状態推定装置。
   （ａ）前記走行状態分布更新部が更新した走行状態分布のうち、予め定めた相対的に時間範囲の短い第２走行状態分布を、前記第２走行状態分布よりも時間的範囲が長い第１走行状態分布で置き換えるリセット処理
   （ｂ）前記外乱となる運転状況から前記外乱とならない運転状況に変化したと判定すると前記走行状態分布更新部による走行状態分布の更新の中断を解除しその時点を前記第２周期の開始点に変更するレジューム処理
2. 前記走行状態分布更新部は、前記データ一時蓄積部が前記第２周期分の分布データを蓄積したと判定すると、各走行状態分布を更新することを特徴とする請求項１に記載の運転状態推定装置。
3. 前記データ一時蓄積部は、前記分布データ取得部が取得した分布データを予め定めた設定数の区分に分類し、分類した分布データを区分毎の度数として蓄積することを特徴とする請求項１または２に記載の運転状態推定装置。
4. 運転者が操作可能な運転操作子の操作状態および車両状態の少なくとも一方を含む走行状態データを取得する走行状態データ取得部と、
   前記走行状態データ取得部が取得した走行状態データに基づき走行状態分布用の分布データを予め定めた第１周期毎に取得する分布データ取得部と、
   前記分布データ取得部が取得した分布データを蓄積するデータ一時蓄積部と、
   前記データ一時蓄積部が蓄積している分布データに基づき前記第１周期よりも長い第２周期毎に時間的範囲の異なる複数の走行状態分布を更新する走行状態分布更新部と、
   前記走行状態分布更新部が更新した複数の走行状態分布に基づき前記運転者の運転状態を推定する運転状態推定部と、
   運転状況を判定する運転状況判定部と、
   前記運転状況判定部が運転状態の推定に対して外乱となる運転状況にあると判定すると、前記走行状態分布更新部が更新する走行状態分布を変更する設定変更部と、
   前記走行状態分布更新部が更新した走行状態分布を蓄積する分布蓄積部と、を備え、
   を備え、
   前記設定変更部は、前記運転状況判定部が運転状態の推定に対して外乱となる運転状況にあると判定すると、（ａ）前記走行状態分布更新部による走行状態分布の更新を中断する更新中断処理、（ｂ）前記走行状態分布更新部が更新した走行状態分布を前記外乱となる運転状況にあると判定しているときの分布データを含まない走行状態分布で置き換えるリストア処理、（ｃ）前記走行状態分布更新部が更新した走行状態分布のうち、予め定めた相対的に時間範囲の短い第２走行状態分布を、前記第２走行状態分布よりも時間的範囲が長い第１走行状態分布で置き換えるリセット処理、（ｄ）前記データ一時蓄積部が蓄積している分布データのうち前記外乱となる運転状況にあると判定する前までの分布データで走行状態分布を更新する更新周期変更処理、及び（ｅ）前記外乱となる運転状況から前記外乱とならない運転状況に変化したと判定すると前記走行状態分布更新部による走行状態分布の更新の中断を解除しその時点を前記第２周期の開始点に変更するレジューム処理、のうち少なくともいずれか１つの処理を実行し、
   前記設定変更部は、前記リストア処理として、前記分布蓄積部が蓄積している走行状態分布のうち、少なくとも２周期以上前の走行状態分布を更新後の走行状態分布とすることを特徴とする運転状態推定装置。
5. 前記設定変更部は、前記リストア処理として、前記運転状況判定部が判定した運転状況に基づき、少なくとも２周期以上前の走行状態分布のうちから走行状態分布を選択し、選択した走行状態分布を更新後の走行状態分布とする走行状態分布を選択することを特徴とする請求項４に記載の運転状態推定装置。
6. 前記運転状況判定部は、前記運転者が操作可能な運転操作子の操作状態、走行環境、および車両状態の少なくともいずれかに基づき運転状況を判定することを特徴とする請求項１から５のいずれか１項に記載の運転状態推定装置。

Images