JP2008249796A - 運転年齢計測システム及び運転年齢計測プログラム - Google Patents

Abstract

【課題】ドライビングシミュレーションにおいて、操舵円滑度指標をしきい値と比較し、円滑度の低い受験者に対して警告を発することで年齢による運転技術の低下を自覚させる。
【解決手段】本発明の運転年齢計測システムは、操舵角センサ３の操舵角信号を時間微分して角速度データを得、曲り角の入点から出点までの期間内に属する角速度データ群をフーリエ変換し、かつ、変換値を２乗して周波数軸上のエネルギ分布を算出し、このエネルギ分布から操舵円滑度指標Ｒを求め、安全運転指標しきい値と比較し、当該操舵円滑度指標が安全運転指標しきい値を超えるときに警告出力することで、ドライビングシミュレーションによって特に高齢ドライバに顕著な指標としての曲り角での操舵にぎこちなさが現われることを警告し、運転年齢の衰えを自覚させる。
【選択図】図２

Description

本発明は、曲り角の操舵の円滑度の良否によって受験者の運転年齢を推定する運転年齢計測システム及び運転年齢計測プログラムに関する。

従来、ドライバの運転操作状況を監視する技術として操舵角に注目し、操舵角データをフーリエ変換したのちの周波数分析を行い、その分析結果に基づいて運転操作状況を判断する技術が特開平９−２７７８４８号公報（特許文献１）に記載されたもの、特開平５−５８１９２号公報（特許文献２）に記載されたものが知られている。また、特開２００１−２５３２６６号公報（特許文献３）、特開２００３−２４６２２７号公報（特許文献４）にも、ドライバの運転能力を判断する技術が記載されている。しかしながら、これらの従来技術はいずれも、道路の曲り角での操舵の円滑さに着目したものでない。

近年の高齢化社会の到来により、高齢ドライバも増加傾向にある。運動能力の低下に伴い、高齢ドライバの場合、反射神経の鈍りによるブレーキ操作の遅れや、道路交通標識の見落としを起こしがちであるが、それと共に、運転操作の円滑さが低下する傾向にある。この運転操作の円滑さの低下、つまり、「運転操作のぎこちなさ」は、特に曲り角の操舵のぎこちなさに顕著に現われる。

本願発明者等は、３０−４５歳の壮年層のドライバ２０名と、６５歳以上の高齢層のドライバ８０名に対して、各自の常用にしている自動車にて一般道を走行してもらい、そのときのドライバ各人の操舵角データを収集し、分析してみた。そして、曲り角での操舵角の微分値、つまり、角速度をフーリエ変換し、さらに２乗した角速度のエネルギ分布に年齢層によって有意差があることを発見した。

図６は走行コースの地図を示している。この走行コースにおいて、曲り角Ｃ２，Ｃ３，Ｃ６，Ｃ７，Ｃ９での操舵角データをサンプリングし、その時間微分によって角速度データを算出し、これをフーリエ変換し、さらに２乗することで図７、図８に示すような周波数対エネルギ分布を求めた。図７は壮年層の一人ドライバの結果、図８は高齢層のドライバ一人の結果を示している。そして、しきい値を周波数０．５／ｓｅｃに設定して、しきい値０．５／ｓｅｃ以上の高域積分値が全域の積分値に示す割合（ここでは、操舵円滑度指標、またＣ ｉｎｄｅｘと呼ぶことにする。）を求め、ドライバ毎に５つの曲り角でのその円滑度指標の平均を求めた。この結果は、図９に示す分布となった。そして、この分布に対して、ウェルチのＴ検定を実施したところ、高齢層ドライバと壮年層ドライバとの間に有意差（Ｐ＜０．０００１）が認められた。
特開平９−２７７８４８号公報 特開平５−５８１９２号公報 特開２００１−２５３２６６号公報 特開２００３−２４６２２７号公報

本発明は、上記実験結果に基づいてなされたもので、ドライビングシミュレーションにおいて、高齢層ドライバと若壮年層ドライバとの間で有意差が認められる操舵円滑度指標を予め設定したしきい値と比較し、円滑度の低い受験者に対して警告を発することで年齢による運転技術の低下を自覚させる運転年齢計測システム及び運転年齢計測プログラムを提供することを目的とする。

本発明の１つの特徴は、少なくとも模擬ステアリングホイールとドライブシミュレーション映像の再生表示装置を備えたドライビングシュミレータと、前記ドライブシミュレーション映像上の曲り角を検出する曲り角検出手段と、前記模擬ステアリングホイールの操舵角を検出し、第１の所定周期で出力する操舵角センサと、前記操舵角センサの操舵角信号を時間微分して角速度データを得る角速度演算手段と、前記曲り角検出手段が検出した曲り角の入点から出点までの期間において、前記角速度演算手段が算出した角速度データをフーリエ変換し、かつ、変換値を２乗して周波数軸上のエネルギ分布を算出するフーリエ変換手段と、前記フーリエ変換手段の算出した前記エネルギ分布のうち、所定の周波数しきい値よりも高い周波数領域の値の積分値（高域積分値）と、該角速度フーリエ変換値の全周波数領域の値の積分値（全域積分値）とを求める積分手段と、前記積分手段の求めた前記高域積分値が前記全域積分値に占める割合を操舵円滑度指標として求める操舵円滑度指標演算手段と、前記操舵円滑度指標演算手段の求めた操舵円滑度指標を安全運転指標しきい値と比較し、当該操舵円滑度指標が安全運転指標しきい値を超えるときに警告出力する運転年齢計測手段とを備えた運転年齢計測システムである。

本発明の別の特徴は、少なくとも模擬ステアリングホイールとドライブシミュレーション映像の再生表示装置を備えたドライビングシュミレータと、前記ドライブシミュレーション映像上の曲り角を検出する曲り角検出手段と、前記模擬ステアリングホイールの操舵角を検出し、第１の所定周期で出力する操舵角センサと、前記操舵角センサの操舵角信号を時間微分して角速度データを得る角速度演算手段と、前記曲り角検出手段が検出した曲り角の入点から出点までの期間において、前記角速度演算手段が算出した角速度データをフーリエ変換し、かつ、変換値を２乗して周波数軸上のエネルギ分布を算出するフーリエ変換手段と、前記フーリエ変換手段の算出した前記エネルギ分布のうち、所定の周波数しきい値よりも高い周波数領域の値の積分値（高域積分値）と、該角速度フーリエ変換値の全周波数領域の値の積分値（全域積分値）とを求める積分手段と、前記積分手段の求めた前記高域積分値と前記全域積分値との比率を操舵円滑度指標として求める操舵円滑度指標演算手段と、操舵円滑度指標と運転年齢指標との対応表手段と、前記対応表手段を参照し、前記操舵円滑度指標演算手段の求めた操舵円滑度指標を運転年齢指標に変換して出力する運転年齢計測手段とを備えた運転年齢計測システムである。

本発明のまた別の特徴は、ドライブシミュレーション映像を再生表示するステップと、前記ドライブシミュレーション映像上の曲り角の入点と出点を検出する曲り角検出ステップと、模擬ステアリングホイールに装着された操舵角センサの操舵角信号を時間微分して角速度データを得る角速度演算ステップと、前記曲り角検出ステップで検出した曲り角の入点から出点までの期間に属する角速度データ群をフーリエ変換し、かつ、変換値を２乗して周波数軸上のエネルギ分布を算出するフーリエ変換ステップと、前記フーリエ変換ステップで算出した前記エネルギ分布のうち、所定の周波数しきい値よりも高い周波数領域の値の積分値（高域積分値）と、該角速度フーリエ変換値の全周波数領域の値の積分値（全域積分値）とを求める積分ステップと、前記積分ステップで求めた前記高域積分値が前記全域積分値に占める割合を操舵円滑度指標として求める操舵円滑度指標演算ステップと、前記操舵円滑度指標演算ステップで求めた操舵円滑度指標を安全運転指標しきい値と比較し、当該操舵円滑度指標が安全運転指標しきい値を超えるときに警告出力する運転年齢計測ステップとをコンピュータに実行させる運転年齢計測プログラムである。

本発明のさらに別の特徴は、ドライブシミュレーション映像を再生表示するステップと、前記ドライブシミュレーション映像上の曲り角の入点と出点を検出する曲り角検出ステップと、模擬ステアリングホイールに装着された操舵角センサの操舵角信号を時間微分して角速度データを得る角速度演算ステップと、前記曲り角検出ステップで検出した曲り角の入点から出点までの期間に属する角速度データ群をフーリエ変換し、かつ、変換値を２乗して周波数軸上のエネルギ分布を算出するフーリエ変換ステップと、前記フーリエ変換ステップで算出した前記エネルギ分布のうち、所定の周波数しきい値よりも高い周波数領域の値の積分値（高域積分値）と、該角速度フーリエ変換値の全周波数領域の値の積分値（全域積分値）とを求める積分ステップと、前記積分ステップで求めた前記高域積分値が前記全域積分値に占める割合を操舵円滑度指標として求める操舵円滑度指標演算ステップと、操舵円滑度指標と運転年齢指標との対応表を参照し、前記操舵円滑度指標演算ステップで求めた操舵円滑度指標を運転年齢指標に変換して出力する運転年齢計測ステップとをコンピュータに実行させる運転年齢計測プログラムである。

上記の各発明において、ドライブシミュレーション映像中に曲り角の入点と出点を示すマークが付しておき、曲り角検出は、ドライブシミュレーション中に曲り角の入点と出点のマークを検出すればその間を曲り角期間とするものとすることができる。

また、上記の各発明において、高域積分値と判定する所定のしきい値は、０．５／ｓｅｃに設定することができる。

さらに、上記の各発明において、安全運転指標しきい値は、０．１０〜０．２５の範囲内のいずれかの値に設定することができる。

本発明によれば、操舵角センサの操舵角信号を時間微分して角速度データを得、曲り角の入点から出点までの期間内に属する角速度データ群をフーリエ変換し、かつ、変換値を２乗して周波数軸上のエネルギ分布を算出し、このエネルギ分布から操舵円滑度指標を求め、安全運転指標しきい値と比較し、当該操舵円滑度指標が安全運転指標しきい値を超えるときに警告出力するので、ドライビングシミュレーションによって特に高齢ドライバに顕著な指標としての曲り角での操舵にぎこちなさが現われることを警告し、運転年齢の衰えを自覚させることができる。

また、本発明によれば、操舵角センサの操舵角信号を時間微分して角速度データを得、曲り角の入点から出点までの期間内に属する角速度データ群をフーリエ変換し、かつ、変換値を２乗して周波数軸上のエネルギ分布を算出し、このエネルギ分布から操舵円滑度指標を求め、対応表を参照して操舵円滑度指標を運転年齢指標に変換して出力するので、ドライビングシミュレーションによって受験者に自分の運転年齢を提示することができ、特に高齢ドライバに運転年齢の衰えを自覚させることができる。

以下、本発明の実施の形態を図に基づいて詳説する。

（第１の実施の形態）本発明の第１の実施の形態の運転年齢計測システムについて説明する。図１に示すように、本実施の形態の運転年齢計測システムは、一般的なドライビングシミュレータ４に備えられている模擬ステアリングホイール１のホイール軸２に取り付けられ、模擬ステアリングホイール１の操舵角瞬時値を検出する操舵角センサ３、操舵角センサ３の操舵角検出データを定周期にサンプリングして時系列に保存するデータストレージ５、このデータストレージ５の保存されたデータを後述する手法にて解析し、操舵円滑度を判定し、警告を出力する運転年齢計測部６、警告を出力するディスプレイ、ブザー、プリンタのような出力装置７から構成されている。尚、データストレージ５と運転年齢計測部６と出力装置７は後述する演算処理を実行するプログラムをインストールしたパーソナルコンピュータとそれに接続される周辺装置にて構成することができる。また、ドライビングシミュレータ４のディスプレイを出力装置として利用することもできる。

ドライビングシミュレータ４は市販されている一般的なものを採用することができる。ただし、そして該ドライビングシミュレータ４が再生表示するドライブシミュレーション映像データには曲り角の入点と出点を示すマークを付す。これはチャプターマークであってもよいし、デジタルデータであるの場合にはデジタル的なマークとすることもできる。

運転年齢計測部６の実行する演算処理方法、演算処理プログラムについて、図２のブロック図、図３のフローチャートを用いて説明する。ドライビングシミュレータ４にてテスト道路の映像を再生表示させ、受験者に模擬ステアリングホイール１を握らせて模擬運転を行わせる。この模擬運転中に、受験者が模擬ステアリングホイール１を操作すれば、操舵角センサ３はホイール軸２の回転角度を操舵角として検出し、データストレージ５に出力し、データストレージ５は所定のサンプリング周期で操舵角データを時系列に蓄積していく。

模擬運転が完了すると、運転年齢計測部６は、データストレージ５に蓄積されている時系列の操舵角データと地点データを用いて曲り角での操舵円滑度指標を算定する。これは、次の手順による。

ステップＳ１：模擬ステアリングホイール１のホイール軸２の部分に装着した操舵角センサ３にて、１／２０ｓｅｃ周期でデジタル値を出力する。

ステップＳ２，Ｓ３：移動平均演算部１１にて３個程度のサンプル数で移動平均を求め、これをデータストレージ５に時系列に保存していく。これは、人間の意図しない成分を排除することを目的としている。システムの簡素化が必要であれば、このステップＳ２，Ｓ３の処理機能は省略することも可能である。次に、リサンプリング部１１にて１／１０ｓｅｃなどの周期でリサンプリングする。

ステップＳ４：データストレージ５から取り込んだ操舵角データに対して、システムの始動時に、操舵角センサ３のゼロ点をキャリブレートする必要がないよう（ハンドルを正置させる必要がないよう）、減衰器１３にて時間と共に出力値が適度にダンピングされる。この機能をソフトウェアとして実装させている。これによりハンドルの取り付けのメカニカルな誤差や、角度センサの取り付け位置のずれを許容できるようになる。ただし減衰係数が大きすぎるとハンドルの戻しのタイミングの測定に誤差が出るなどの悪影響があるので、係数は「弱め」に設定する。この処理ステップＳ４もシステムの簡素が必要な場合には省略することが可能である。

ステップＳ５：差分器１４にてサンプリングタイミング毎に、１つ前のサンプリングタイミングでの操舵角データとの差分を取ることで角速度データに変換し、データストレージ５に保存していく。

ステップＳ６，Ｓ７：曲り角検出部１０はテスト道路のドライブシミュレーション映像データ８に付された曲り角の入点と出点のマークを検出する。そして、データ切出し部１５は、この曲り角期間内に属する角速度データ群を切り出す。

ステップＳ８，Ｓ９：ＦＦＴ演算部１６にて、切り出された角速度データ群に対して、折り返し効果を低減するためにハミング窓関数を乗じる。そして、離散フーリエ変換ＦＦＴを行う。この結果は複素数のベクトルである。

ステップＳ１０：２乗器１７にてＦＦＴの結果を２乗して、時間軸上のエネルギ分布を表すベクトル（実数のベクトル）に変換する。

ステップＳ１１：比率（操舵円滑度指標）算出部１８にて、実数のベクトルを全周波数域で積分して全域積分値Ｓ１を求め、また、しきい値Ｒｅｆ１（ここでは、最適値として０．５／ｓｅｃに設定している。）以上の周波数域でも積分して高域積分値Ｓ２を求め、この高域積分値Ｓ２の全域積分値Ｓ１に対する比率Ｒ（＝Ｓ２／Ｓ１）を操舵円滑度指標として算出する。尚、Ｒｅｆ１は、可変設定若しくは微調整可能にすることができる。

ステップＳ１２：操舵円滑度判定部１９にて、ステップＳ１０で求めた操舵円滑度指標Ｒを運転年齢判定指標Ｒｅｆ２（ここでは、０．１０〜０．２５の範囲で可変設定できる。）と比較し、このＲｅｆ２を超える場合に運転能力が低下している高齢ドライバと判定し、警告出力部２０にて、例えば、「運転能力が顕著に低下してきています。十分に注意して運転なさって下さい。」とか、「運転能力が極度に低下してきています。運転免許証を返上をされることをお勧めします。」といった警告を出力させる。この警告出力は、表示装置７に画面表示し、また、プリントアウトして受験者本人に手渡しできるようにする。

これにより、本実施の形態の運転年齢計測システム、これにて実行される運転年齢計測プログラムによれば、運転能力が低下している高齢ドライバに運転能力の低下を自覚させることができる。

尚、ステップＳ５〜Ｓ７の角速度演算、曲り角検出及びデータの切り出しについては、次の手順であってもよい。まず、ステップＳ５で曲り角検出を先にし、ステップＳ６で操舵角データのうち曲り角走行期間内に属するデータ群を切り出し、ステップＳ７で切り出した操舵角データ群に対して差分を実施して曲り角走行期間内に属する角速度データ群を得、これをデータストレージ５に保存し、ステップＳ８のＦＦＴをこの角速度データ群に対して実行することになる。結果は同じである。

（第２の実施の形態）本発明の第２の実施の形態の運転年齢計測システム、これにて実行する運転年齢計測プログラムについて説明する。

本実施の形態の特徴は、受験者の操舵円滑度指標Ｒをランク分けし、ランク毎の運転年齢を推定して受験者に提示する点にある。本実施の形態の運転年齢計測システムのハードウェア構成は、第１の実施の形態と同様に図１に示すものである。そして、運転年齢計測部６の機能構成は図４に示してあり、図２に示した第１の実施の形態における警告出力部２０に代えて、運転年齢ランキング部２２を備えている。尚、図４において、第１の実施の形態と共通する要素には同一の符号を用いて示してある。

本実施の形態によれば、図５のフローチャートにおけるステップＳ１２’で、次のような処理を行う。すなわち、操舵円滑度判定部１９のＲｅｆ２を複数段に設定し、操舵円滑度指標算出部１８が求めた操舵円滑度指標Ｒがいずれのランクに属するかを判定する。そして運転年齢ランキング部２２は、そのランクに対応して、Ｒ＝０．１を超えていれば「運転年齢６０歳を超えています。運転能力低下傾向にありますので、十分慎重に運転なさって下さい。」、さらにＲ＝０．１５を超えていれば「運転年齢７０歳を超えています。運転能力が低下しています。運転は中止された方が良いです。」、そしてさらにＲ＝０．２０を超えていれば「運転年齢８０歳を超えています。運転は無理です。免許証の返上をお勧めします。」といった運転年齢のランクに応じた警告を出力する。

尚、この図５のフローチャートにおいて、図３に示したフローチャートと共通するステップ番号は同様の演算処理を行う。また第１の実施の形態と同様に、ステップＳ５〜Ｓ７の処理は入れ替えることができる。

これにより、本実施の形態によれば、曲り角での操舵円滑度指標のランクに応じて運転年齢を提示すると共に警告も出力するようにしたので、運転能力が低下している高齢ドライバに運転能力の低下を自覚させることができる。

本発明の第１の実施の形態の運転年齢計測システムのブロック図。 本発明の第１の実施の形態における運転年齢計測部の機能ブロック図。 本発明の第１の実施の形態による運転年齢計測処理のフローチャート。 本発明の第２の実施の形態における運転年齢計測部の機能ブロック図。 本発明の第２の実施の形態による運転年齢計測処理のフローチャート。 本発明の実用性の立証のための設定したテストコースの地図。 壮年層ドライバの操舵角速度のエネルギ分布のグラフ。 高齢ドライバの操舵角速度のエネルギ分布のグラフ。 壮年層ドライバと高齢ドライバとの操舵円滑度指標の分布を示すグラフ。

符号の説明

１…模擬ステアリングホイール
２…ステアリングホイール軸
３…操舵角センサ
４…ドライビングシミュレータ
５…データストレージ
６…運転年齢計測部
７…出力部
８…道路地図データ
１０…曲り角検出部
１４…差分器
１６…ＦＦＴ演算器
１７…２乗器
１８…操舵円滑度指標算出部
１９…操舵円滑度判定部
２０…警告出力部
２２…運転年齢ランキング部

Claims (10)

1. 少なくとも模擬ステアリングホイールとドライブシミュレーション映像の再生表示装置を備えたドライビングシュミレータと、
   前記ドライブシミュレーション映像上の曲り角を検出する曲り角検出手段と、
   前記模擬ステアリングホイールの操舵角を検出し、第１の所定周期で出力する操舵角センサと、
   前記操舵角センサの操舵角信号を時間微分して角速度データを得る角速度演算手段と、
   前記曲り角検出手段が検出した曲り角の入点から出点までの期間において、前記角速度演算手段が算出した角速度データをフーリエ変換し、かつ、変換値を２乗して周波数軸上のエネルギ分布を算出するフーリエ変換手段と、
   前記フーリエ変換手段の算出した前記エネルギ分布のうち、所定の周波数しきい値よりも高い周波数領域の値の積分値（高域積分値）と、該角速度フーリエ変換値の全周波数領域の値の積分値（全域積分値）とを求める積分手段と、
   前記積分手段の求めた前記高域積分値が前記全域積分値に占める割合を操舵円滑度指標として求める操舵円滑度指標演算手段と、
   前記操舵円滑度指標演算手段の求めた操舵円滑度指標を安全運転指標しきい値と比較し、当該操舵円滑度指標が安全運転指標しきい値を超えるときに警告出力する運転年齢計測手段とを備えた運転年齢計測システム。
2. 少なくとも模擬ステアリングホイールとドライブシミュレーション映像の再生表示装置を備えたドライビングシュミレータと、
   前記ドライブシミュレーション映像上の曲り角を検出する曲り角検出手段と、
   前記模擬ステアリングホイールの操舵角を検出し、第１の所定周期で出力する操舵角センサと、
   
   前記操舵角センサの操舵角信号を時間微分して角速度データを得る角速度演算手段と、
   前記曲り角検出手段が検出した曲り角の入点から出点までの期間において、前記角速度演算手段が算出した角速度データをフーリエ変換し、かつ、変換値を２乗して周波数軸上のエネルギ分布を算出するフーリエ変換手段と、
   前記フーリエ変換手段の算出した前記エネルギ分布のうち、所定の周波数しきい値よりも高い周波数領域の値の積分値（高域積分値）と、該角速度フーリエ変換値の全周波数領域の値の積分値（全域積分値）とを求める積分手段と、
   前記積分手段の求めた前記高域積分値と前記全域積分値との比率を操舵円滑度指標として求める操舵円滑度指標演算手段と、
   操舵円滑度指標と運転年齢指標との対応表手段と、
   前記対応表手段を参照し、前記操舵円滑度指標演算手段の求めた操舵円滑度指標を運転年齢指標に変換して出力する運転年齢計測手段とを備えた運転年齢計測システム。
3. 前記ドライブシミュレーション映像中に曲り角の入点と出点を示すマークが付されており、
   前記曲り角検出手段は、前記曲り角の入点と出点のマークを検出することで曲り角期間を検出することを特徴とする請求項１又は２に記載の運転年齢計測システム。
4. 前記積分手段が前記高域積分値と判定する所定のしきい値は、０．５／ｓｅｃにしたことを特徴とする請求項１〜３のいずれかに記載の運転年齢計測システム。
5. 前記安全運転指標しきい値は、０．１０〜０．２５の範囲内のいずれかの値に設定したことを特徴とする請求項１に記載の運転年齢計測システム。
6. ドライブシミュレーション映像を再生表示するステップと、
   前記ドライブシミュレーション映像上の曲り角の入点と出点を検出する曲り角検出ステップと、
   模擬ステアリングホイールに装着された操舵角センサの操舵角信号を時間微分して角速度データを得る角速度演算ステップと、
   前記曲り角検出ステップで検出した曲り角の入点から出点までの期間に属する角速度データ群をフーリエ変換し、かつ、変換値を２乗して周波数軸上のエネルギ分布を算出するフーリエ変換ステップと、
   前記フーリエ変換ステップで算出した前記エネルギ分布のうち、所定の周波数しきい値よりも高い周波数領域の値の積分値（高域積分値）と、該角速度フーリエ変換値の全周波数領域の値の積分値（全域積分値）とを求める積分ステップと、
   前記積分ステップで求めた前記高域積分値が前記全域積分値に占める割合を操舵円滑度指標として求める操舵円滑度指標演算ステップと、
   前記操舵円滑度指標演算ステップで求めた操舵円滑度指標を安全運転指標しきい値と比較し、当該操舵円滑度指標が安全運転指標しきい値を超えるときに警告出力する運転年齢計測ステップとをコンピュータに実行させることを特徴とする運転年齢計測プログラム。
7. ドライブシミュレーション映像を再生表示するステップと、
   前記ドライブシミュレーション映像上の曲り角の入点と出点を検出する曲り角検出ステップと、
   模擬ステアリングホイールに装着された操舵角センサの操舵角信号を時間微分して角速度データを得る角速度演算ステップと、
   前記曲り角検出ステップで検出した曲り角の入点から出点までの期間に属する角速度データ群をフーリエ変換し、かつ、変換値を２乗して周波数軸上のエネルギ分布を算出するフーリエ変換ステップと、
   前記フーリエ変換ステップで算出した前記エネルギ分布のうち、所定の周波数しきい値よりも高い周波数領域の値の積分値（高域積分値）と、該角速度フーリエ変換値の全周波数領域の値の積分値（全域積分値）とを求める積分ステップと、
   前記積分ステップで求めた前記高域積分値が前記全域積分値に占める割合を操舵円滑度指標として求める操舵円滑度指標演算ステップと、
   操舵円滑度指標と運転年齢指標との対応表を参照し、前記操舵円滑度指標演算ステップで求めた操舵円滑度指標を運転年齢指標に変換して出力する運転年齢計測ステップとをコンピュータに実行させることを特徴とする運転年齢計測プログラム。
8. 前記ドライブシミュレーション映像中に曲り角の入点と出点を示すマークが付されており、
   前記曲り角検出ステップては、前記ドライブシミュレーション映像中に付されている曲り角の入点と出点のマークを検出することで曲り角期間を検出することを特徴とする請求項６又は７に記載の運転年齢計測プログラム。
9. 前記積分ステップでは、前記高域積分値と判定する所定のしきい値を０．５／ｓｅｃにしたことを特徴とする請求項６〜８のいずれかに記載の運転年齢計測プログラム。
10. 前記安全運転指標しきい値は、０．１０〜０．２５の範囲内のいずれかの値に設定したことを特徴とする請求項９に記載の運転年齢計測プログラム。

Images