[go: up one dir, main page]

JP2010198120A - Driving condition determination device and method - Google Patents

Driving condition determination device and method Download PDF

Info

Publication number
JP2010198120A
JP2010198120A JP2009039631A JP2009039631A JP2010198120A JP 2010198120 A JP2010198120 A JP 2010198120A JP 2009039631 A JP2009039631 A JP 2009039631A JP 2009039631 A JP2009039631 A JP 2009039631A JP 2010198120 A JP2010198120 A JP 2010198120A
Authority
JP
Japan
Prior art keywords
driving
characteristic
information
model
state
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2009039631A
Other languages
Japanese (ja)
Other versions
JP5267203B2 (en
Inventor
Takayuki Kondo
崇之 近藤
Toshihiro Yamamura
智弘 山村
Kazuto Sato
一人 佐藤
Hirobumi Hashiguchi
博文 橋口
Masachika Kojima
正親 小島
Nobutomo Hisaie
伸友 久家
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nissan Motor Co Ltd
Original Assignee
Nissan Motor Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nissan Motor Co Ltd filed Critical Nissan Motor Co Ltd
Priority to JP2009039631A priority Critical patent/JP5267203B2/en
Publication of JP2010198120A publication Critical patent/JP2010198120A/en
Application granted granted Critical
Publication of JP5267203B2 publication Critical patent/JP5267203B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Traffic Control Systems (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To provide a technique for determining the driving conditions of a driver with a few data to be acquired in a short period of time. <P>SOLUTION: A driving condition determination device is configured to calculate the featured quantities of the driving characteristics of a driver as the information of the predicted error distribution of a driving operation based on driving operation information, and to acquire a driving characteristic model after learning based on a driving characteristic model comprising the calculated featured quantities of the driving conditions and the information of the predicted error distribution of the driving conditions becoming a preset reference, and to determine whether or not the current driving conditions are unstable driving conditions different from normal driving characteristics based on the acquired driving characteristic model after learning and the featured quantities of the current driving characteristics. <P>COPYRIGHT: (C)2010,JPO&INPIT

Description

本発明は、運転者の運転支援のために、運転者の運転状態を判定する技術に関する。   The present invention relates to a technique for determining a driving state of a driver for driving support of the driver.

運転者の運転を支援する技術としては、例えば特許文献1に記載の技術がある。この技術は、運転者毎に記憶した走行パターンと最新の走行状態とを比較して、走行状態が正常か異常か判定する。異常と判別した場合には警報を出す。ここで、逐次取得する運転者の操作と走行情報とに基づき補正値を更新する。そして、逐次更新する補正値を基本マップに加えてなる走行マップを、上記判定用の走行パターンとして使用する。   As a technique for supporting the driving of the driver, there is a technique described in Patent Document 1, for example. In this technology, the traveling pattern stored for each driver is compared with the latest traveling state to determine whether the traveling state is normal or abnormal. When it is determined that there is an abnormality, an alarm is issued. Here, the correction value is updated based on the driver's operation and travel information acquired sequentially. Then, a travel map in which correction values that are sequentially updated are added to the basic map is used as the travel pattern for determination.

特開2005−301832号公報(段落番号0019、0020参照)Japanese Patent Laying-Open No. 2005-301832 (see paragraphs 0019 and 0020)

上記のような従来技術では、実用可能な精度以上の走行パターンや運転者特性を取得するまでには、所要の時間の学習が必要となる。これは、短時間のデータでは走行環境などに大きく影響し、ノイズの多い情報となるためである。このため、運転状態の判定が可能となるまでに時間が掛かる。
本発明は、上記のような点に着目してなされたもので、短時間に得られる少ないデータで、運転者の運転状態の判定が可能な運転状態判定の技術を提示することを課題とする。
In the conventional technology as described above, it is necessary to learn a required time before acquiring a driving pattern or driver characteristics that are more than practical accuracy. This is because short-time data greatly affects the driving environment and becomes information with a lot of noise. For this reason, it takes time until the operating state can be determined.
The present invention has been made paying attention to the above points, and it is an object of the present invention to provide a driving state determination technique capable of determining a driving state of a driver with a small amount of data obtained in a short time. .

上記課題を解決するために、本発明は、運転操作情報に基づき、運転者の運転特性の特徴量を運転操作の予測誤差分布の情報として算出する。算出した運転特性の特徴量、及び、予め設定した基準となる運転操作の予測誤差分布の情報からなる運転特性モデルに基づき、学習後の運転特性モデルを取得する。そして、取得した学習後の運転特性モデルと、現在の運転特性の特徴量とに基づき、現在の運転状態が普段の運転特性とは異なる不安定運転状態か否かを判定する。   In order to solve the above-described problem, the present invention calculates a feature amount of a driving characteristic of a driver as information on a prediction error distribution of the driving operation based on the driving operation information. A learned driving characteristic model is acquired based on the driving characteristic model including the calculated characteristic value of the driving characteristic and information on the prediction error distribution of the driving operation as a reference set in advance. Then, based on the acquired driving characteristic model after learning and the characteristic amount of the current driving characteristic, it is determined whether or not the current driving state is an unstable driving state different from the normal driving characteristic.

本発明によれば、運転者の運転特性の特徴量として運転操作の予測誤差分布の情報を使用する。このため、学習時間を短く出来る。また、予め基準となる運転操作の予測誤差分布の情報からなる運転特性モデルを備えることで、更に学習時間を短く出来る。
以上の事から、短時間に得られる少ないデータで、運転者の運転状態の判定が可能となる。
According to the present invention, the information on the prediction error distribution of the driving operation is used as the feature amount of the driving characteristic of the driver. For this reason, learning time can be shortened. In addition, the learning time can be further shortened by providing a driving characteristic model including information on the prediction error distribution of the driving operation as a reference in advance.
From the above, it is possible to determine the driving state of the driver with a small amount of data obtained in a short time.

本発明に基づく第1実施形態に係るシステム構成を示す図である。It is a figure which shows the system configuration | structure which concerns on 1st Embodiment based on this invention. 本発明に基づく第1実施形態に係る情報提示コントローラの構成を示す図である。It is a figure which shows the structure of the information presentation controller which concerns on 1st Embodiment based on this invention. 運転特性モデルを示す図である。It is a figure which shows a driving | running characteristic model. 運転特性モデルを示す図である。It is a figure which shows a driving | running characteristic model. 本発明に基づく第1実施形態に係る情報提示コントローラの処理を説明する図である。It is a figure explaining the process of the information presentation controller which concerns on 1st Embodiment based on this invention. 本発明に基づく第1実施形態に係る不安定指標の算出処理を説明する図である。It is a figure explaining the calculation process of the unstable index which concerns on 1st Embodiment based on this invention. 本発明に基づく第1実施形態に係る操作誤差分布の算出処理を説明する図である。It is a figure explaining the calculation process of the operation error distribution which concerns on 1st Embodiment based on this invention. 本発明に基づく第1実施形態に係る運転特性特徴量の算出処理を説明する図である。It is a figure explaining the calculation process of the driving characteristic feature-value based on 1st Embodiment based on this invention. 予測誤差分布を説明する図である。It is a figure explaining prediction error distribution. 予測誤差分布を説明する図である。It is a figure explaining prediction error distribution. 予測誤差の頻度分布を示す図である。It is a figure which shows frequency distribution of prediction error. 本発明に基づく第2実施形態に係る情報提示コントローラの処理を説明する図である。It is a figure explaining the process of the information presentation controller which concerns on 2nd Embodiment based on this invention. 本発明に基づく第3実施形態に係る運転特性特徴量の算出処理を説明する図である。It is a figure explaining the calculation process of the driving characteristic feature-value based on 3rd Embodiment based on this invention. 補正係数を説明する図である。It is a figure explaining a correction coefficient. 補正係数を説明する図である。It is a figure explaining a correction coefficient.

次に、本発明の実施形態を図面を参照しつつ説明する。
図1は、本実施形態に係る運転状態判定装置を含む情報提示装置を備えた車両のシステム構成図である。
本実施形態では、運転操作としてステアリング操作を対象とする場合を例に挙げて説明する。そして、運転操作の情報として操舵角度を使用する。
Next, embodiments of the present invention will be described with reference to the drawings.
FIG. 1 is a system configuration diagram of a vehicle including an information presentation device including a driving state determination device according to the present embodiment.
In the present embodiment, a case where a steering operation is targeted as a driving operation will be described as an example. Then, the steering angle is used as driving operation information.

(構成)
本実施形態の車両は、普段の車両と同様に、運転者の操作に応じた状態に車両挙動を制御する。すなわち、運転者が操作するアクセルペダルのアクセル開度に応じた駆動力を発生するように、エンジンやモータなどの駆動輪を駆動する駆動源を制御する。また、運転者が操作するブレーキペダルのブレーキ開度に応じた制動力が発生するように、各輪の制動装置を制御する。また、運転者が操作するステアリングホイールの操舵角に応じて、ステアリング装置が操向輪を転舵する。
(Constitution)
The vehicle according to the present embodiment controls the vehicle behavior in a state corresponding to the driver's operation, like the ordinary vehicle. That is, a drive source for driving drive wheels such as an engine and a motor is controlled so as to generate a driving force corresponding to the accelerator opening of the accelerator pedal operated by the driver. Further, the braking device for each wheel is controlled so that a braking force according to the brake opening degree of the brake pedal operated by the driver is generated. Further, the steering device steers the steered wheels according to the steering angle of the steering wheel operated by the driver.

上記のような車両に対し、運転状態判定装置を含む情報提示装置を設ける。以下、この情報提示装置に関する説明を行う。
操舵角センサ1は、ステアリングホイールの操舵角を検出する。操舵角センサ1は、検出した操舵角を情報提示コントローラ2に出力する。
情報提示コントローラ2は、操舵角に基づく操作誤差分布の情報によって、運転者の普段の運転状態を学習する。学習が完了したら、学習した普段の運転状態に基づき、現在の運転状態が不安定運転状態か否かを判定する。そして、情報提示コントローラ2は、現在の運転状態が不安定運転状態と判定すると、運転者にその旨の報知を行う。すなわち、上記情報提示コントローラ2は、短時間の運転行動の情報から、運転者の(普段の)運転特性を推定する。そして、この推定結果を学習後の初期値とすることで、早期に不安定状態を判定し、情報提示する。また、モデルのパラメータの同定後(学習完了後)も、逐次学習を行う。
An information presentation device including a driving state determination device is provided for the vehicle as described above. Hereinafter, the information presentation apparatus will be described.
The steering angle sensor 1 detects the steering angle of the steering wheel. The steering angle sensor 1 outputs the detected steering angle to the information presentation controller 2.
The information presentation controller 2 learns the driver's normal driving state based on the operation error distribution information based on the steering angle. When the learning is completed, it is determined whether or not the current driving state is an unstable driving state based on the learned normal driving state. And if the information presentation controller 2 determines with the present driving | running state being an unstable driving | running state, it will alert | report that to a driver | operator. That is, the information presentation controller 2 estimates the driver's (ordinary) driving characteristics from information on driving behavior in a short time. And by making this estimation result into the initial value after learning, an unstable state is determined early and information is presented. Further, sequential learning is performed after the identification of model parameters (after learning is completed).

上記情報提示コントローラ2は、運転状態判定装置2A及び情報発信手段2Bを備える。運転状態判定装置2Aは、基準となる運転特性モデル2Aa、運転特性特徴量算出手段2Ab、モデル取得手段2Ac、及び不安定運転状態判定手段2Adを備える。
基準となる運転情報モデルは、予め設定した基準となる運転操作の予測誤差分布の情報からなる。この情報は、実験などによって取得した複数人の運転操作情報に基づき作成すれば良い。図3及び図4にその例を示す。
The information presentation controller 2 includes a driving state determination device 2A and information transmission means 2B. The driving state determination device 2A includes a driving characteristic model 2Aa serving as a reference, a driving characteristic feature amount calculation unit 2Ab, a model acquisition unit 2Ac, and an unstable driving state determination unit 2Ad.
The reference driving information model includes information on a prediction error distribution of a driving operation that is set in advance as a reference. This information may be created based on driving operation information of a plurality of persons acquired through experiments or the like. Examples thereof are shown in FIGS.

運転特性特徴量算出手段2Abは、運転操作情報を取得する手段が取得する運転操作情報に基づき、運転者の運転特性の特徴量を運転操作の予測誤差分布の情報として算出する。
モデル取得手段2Acは、上記運転特性特徴量算出手段2Abが算出した運転特性の特徴量、及び、予め設定した上記基準となる運転特性モデル2Aaに基づき、学習後の運転特性モデルMを取得する。
The driving characteristic feature amount calculation means 2Ab calculates the characteristic amount of the driving characteristic of the driver as information on the prediction error distribution of the driving operation based on the driving operation information acquired by the means for acquiring the driving operation information.
The model acquisition unit 2Ac acquires the learned driving characteristic model M based on the driving characteristic feature amount calculated by the driving characteristic feature amount calculating unit 2Ab and the preset driving characteristic model 2Aa.

不安定運転状態判定手段2Adは、モデル取得手段2Acが取得した学習後の運転特性モデルMと、上記運転特性特徴量算出手段2Abが算出した現在の特徴量とに基づき、現在の運転状態が普段の運転特性とは異なる不安定運転状態か否かを判定する。
情報発信手段2Bは、不安定運転状態判定手段2Adの判定に応じた運転状態情報を情報提示手段3に出力する。
そして、情報提示手段3は、情報発信手段2Bからの運転状態情報に応じた情報を運転者に提示する。情報提示手段3は、例えば、ナビ等の表示部や音声発生装置、座席を振動させる振動装置等によって構成する。
The unstable driving state determination means 2Ad is based on the learned driving characteristic model M acquired by the model acquisition means 2Ac and the current characteristic amount calculated by the driving characteristic feature amount calculation means 2Ab. It is determined whether or not the unstable operation state is different from the operation characteristic of.
The information transmission unit 2B outputs the driving state information corresponding to the determination of the unstable driving state determination unit 2Ad to the information presentation unit 3.
Then, the information presenting means 3 presents information corresponding to the driving state information from the information transmitting means 2B to the driver. The information presenting means 3 is configured by, for example, a display unit such as a navigation, a sound generation device, a vibration device that vibrates the seat, or the like.

次に、上記情報提示コントローラ2の処理を図5を参照して説明する。
この処理は、所定サンプリング周期毎に作動する。
まず、ステップS100では、運転操作の情報を習得する。本実施形態では、操舵角センサ1からの入力に基づき、操舵角を取得する。
ステップS200では、取得した操舵角を運転行動として記録する。この記録は、運転者の状態を正確に推定するため、過去のデータの統計値・分布や、直近の運転行動から現在の運転行動・運転状態を推定するために用いる。
ステップS300では、不安定運転指標の算出を行う。本実施形態では、ステアリングエントロピー法を用いて、不安定運転指標を算出する。不安定運転指標の算出には後述する。不安定運転指標は、学習用の普段の不安定運転指標と、現在(このとき)の不安定運転指標を個別に算出する。
Next, the process of the information presentation controller 2 will be described with reference to FIG.
This process operates every predetermined sampling period.
First, in step S100, driving operation information is acquired. In the present embodiment, the steering angle is acquired based on the input from the steering angle sensor 1.
In step S200, the acquired steering angle is recorded as a driving action. This record is used to estimate the current driving behavior / driving state from the statistical values / distribution of past data and the latest driving behavior in order to accurately estimate the driver's state.
In step S300, an unstable operation index is calculated. In this embodiment, an unstable driving index is calculated using a steering entropy method. The calculation of the unstable operation index will be described later. As the unstable driving index, a normal unstable driving index for learning and a current (at this time) unstable driving index are calculated separately.

ここで、ステアリングエントロピー法について説明する。
一般的に、運転者の注意が運転に集中していない状態では、操舵が行われない時間が運転に集中した正常運転時よりも長くなり、大きな操舵角の誤差が蓄積される。したがって、運転者の注意が運転に戻ったときの修正操舵量が大きくなる。ステアリングエントロピー法は、この特性に着目したものであり、特性値としてα値と、α値を基準に算出した操舵角エントロピーHpを用いる。基準となる操舵角エントロピーと、計測された操舵角に基づいて算出した操舵角エントロピーとを比較することにより、基準に対する運転操作の不安定な状態を検出する。
Here, the steering entropy method will be described.
In general, in a state where the driver's attention is not concentrated on driving, the time during which steering is not performed becomes longer than in normal driving where the driving is concentrated, and a large steering angle error is accumulated. Therefore, the corrected steering amount when the driver's attention returns to driving increases. The steering entropy method focuses on this characteristic, and uses an α value as a characteristic value and a steering angle entropy Hp calculated based on the α value. By comparing the reference steering angle entropy with the steering angle entropy calculated based on the measured steering angle, an unstable state of the driving operation with respect to the reference is detected.

なお、特性値としてのα値は、操舵角の時系列データに基づいて一定時間内の操舵誤差を求め、その操舵誤差の分布(ばらつき)を測定して90パーセントタイル値(操舵誤差の90%が含まれる分布の範囲)を算出したものである。上記操舵誤差は、ステアリングを滑らかに操作したと仮定した場合の操舵角の推定値と、実際の操舵角との差である。
ステアリングエントロピー値としてのHp値は、操舵誤差分布のあいまいさ(不確実性)を表す。Hp値は、α値と同様に、ステアリング操作が滑らかで安定している場合は小さくなり、ガクガクと不安定な場合は大きくなる。Hp値はα値によって補正することで、運転者の技量や癖により影響を受けにくい運転者不安定度として用いることができる。
Note that the α value as the characteristic value is a 90% tile value (90% of the steering error) by calculating a steering error within a predetermined time based on the time-series data of the steering angle and measuring the distribution (variation) of the steering error. The range of the distribution in which is included. The steering error is a difference between an estimated value of the steering angle when it is assumed that the steering is operated smoothly and an actual steering angle.
The Hp value as the steering entropy value represents the ambiguity (uncertainty) of the steering error distribution. Similar to the α value, the Hp value decreases when the steering operation is smooth and stable, and increases when it is unstable and unstable. By correcting the Hp value by the α value, it can be used as a driver instability that is not easily affected by the skill or habit of the driver.

なお、操舵誤差は、道路線形や運転者の負荷の状態に影響を受けるため、運転者が無負荷の状態で、基準となる操舵誤差分布を計測する必要がある。そこで、長時間計測した操舵角データに基づいて操舵誤差分布を求める。
次にステップS400では、学習用の普段の不安定運転指標に基づき運転特性特徴量の算出を行う。この運転特性特徴量の算出については後述する。
Since the steering error is affected by the road alignment and the load state of the driver, it is necessary to measure a reference steering error distribution when the driver is in an unloaded state. Therefore, a steering error distribution is obtained based on the steering angle data measured for a long time.
Next, in step S400, the driving characteristic feature amount is calculated based on the usual unstable driving index for learning. The calculation of the driving characteristic feature amount will be described later.

次にステップS500では、学習完了か否かの判定を行う。すなわち、パラメータの同定が完了しているか否かを判定する。運転者特性の推定処理において、モデルを適合し初期値が入っている場合は、学習完了としてステップS1410へ進む。一方で適合していない場合は、処理を終了して復帰する。学習完了か否かは、例えば、上記運転特性特徴量の算出の処理を何回実施したかをカウントし、そのカウントが所定値を越えると、学習が完了と判定する。「普段」のための学習時間を2160秒とし、サンプリング数を5Hzとすると、データ個数Kは10800となる。   Next, in step S500, it is determined whether learning is completed. That is, it is determined whether or not parameter identification has been completed. In the driver characteristic estimation process, if the model is adapted and an initial value is entered, the learning is completed and the process proceeds to step S1410. On the other hand, if it does not match, the process is terminated and the process returns. Whether or not the learning is completed is, for example, counting how many times the calculation process of the driving characteristic feature amount is performed, and when the count exceeds a predetermined value, it is determined that the learning is completed. If the learning time for “ordinary” is 2160 seconds and the sampling number is 5 Hz, the number of data K is 10800.

ステップS600では、相対エントロピーHpを算出する。
学習した運転特性モデルの予測誤差の分布(pl)と、現在の予測誤差の分布(ps)とに基づき、下記式によって相対エントロピーR_Hplsを算出する。ここで不安定状態の検出は、普段に対する現在このとき(直近)の比較である。
In step S600, the relative entropy Hp is calculated.
Based on the learned prediction error distribution (pl i ) and the current prediction error distribution (ps i ), the relative entropy R_Hpls is calculated by the following equation. Here, the detection of the unstable state is a comparison at this time (most recent) with respect to the usual.

Figure 2010198120
Figure 2010198120

ステップS700では、不安定状態かを判定する。ここでは、相対Hp(R−Hp)が所定の閾値(例えば閾値=0.5)を超えていたら、ステップS800へ進む。一方で、超えていない場合は、処理を終了して復帰する。
ステップS800では、不安定状態の情報提示を行う。情報提示は、情報提示手段3を用いて、視覚、音声などから情報提示を行う。
次に、上記ステップS300にて行う、不安定運転指標の算出処理について、図6を参照して説明する。
先ずステップS310で、予測誤差PEの算出を行う。
その予測誤差PEは、下記式のように、予測操舵角θpと実際の操舵角θの差として求める。
PE= θ −θp
In step S700, it is determined whether the state is unstable. Here, if the relative Hp (R−Hp) exceeds a predetermined threshold (for example, threshold = 0.5), the process proceeds to step S800. On the other hand, when it does not exceed, the process is terminated and the process returns.
In step S800, unstable state information is presented. For information presentation, information presentation means 3 is used to present information from vision, voice, or the like.
Next, the unstable operation index calculation process performed in step S300 will be described with reference to FIG.
First, in step S310, the prediction error PE is calculated.
The prediction error PE is obtained as a difference between the predicted steering angle θp and the actual steering angle θ as in the following equation.
PE = θ−θp

予測操舵角θpは、過去の操舵角に基づき、滑らかに操舵したら得られるであろうと推定する操舵角である。例えば、現在の時点をnとした場合に、の直近3点(n−3、n−2、n−1)の過去の操舵角を用いてn−1時点を中心とする2次テイラー展開により算出する。具体的には、過去3回の150msごと(150ms前と300ms前と450ms前)の操舵角θa1、θa2、θa3から、下記式によって、予測操舵角θpを算出する。
θp=θa1+(θa1−θa2)
+(1/2)(((θa1−θa2)−(θa2−θa3))
次にステップS320では、算出した予測操舵角に基づき、2つの予測誤差分布を更新する。2つの予測誤差分布は、普段の予測誤差分布と、現在(このとき)の予測誤差分布である。ここで予測誤差分布は、頻度分布で表す。
The predicted steering angle θp is a steering angle that is estimated to be obtained by smooth steering based on the past steering angle. For example, when the current time point is n, secondary Taylor expansion around the n-1 time point using the past three steering angles of the last three points (n-3, n-2, n-1). calculate. Specifically, the predicted steering angle θp is calculated from the steering angles θa1, θa2, and θa3 every 150 ms (150 ms before, 300 ms before, and 450 ms before) by the following formula.
θp = θa1 + (θa1-θa2)
+ (1/2) (((θa1-θa2)-(θa2-θa3))
Next, in step S320, two prediction error distributions are updated based on the calculated predicted steering angle. The two prediction error distributions are a normal prediction error distribution and a current (at this time) prediction error distribution. Here, the prediction error distribution is represented by a frequency distribution.

この2つの予測誤差分布の各分布割合pl、psの更新について、図7を参照して説明する。
ここで、2つの予測誤差分布は、下記のように予め設定する。
2つの予測誤差分布を更新するためのそれぞれのタイムウインドウを設定する。普段用のタイムウインドウは、Nlとして2160秒(所定値)を設定する。また、「現在用のタイムウインドウは、Nsとして60秒(所定値)を設定する。また、予測誤差分布のビンの間隔を決定する値αを、所定値0.64に設定する。
Each distribution ratio pl i of the two prediction error distribution for updating ps i, be described with reference to FIG.
Here, the two prediction error distributions are preset as follows.
Each time window for updating the two prediction error distributions is set. The usual time window is set to 2160 seconds (predetermined value) as Nl. “The current time window is set to 60 seconds (predetermined value) as Ns. Further, the value α for determining the bin interval of the prediction error distribution is set to the predetermined value 0.64.

そして、まずステップS321で、αに基づき、2つの予測誤差分布の各ビンの幅(B)をそれぞれ設定する。ビンの幅は、下記式のように、bin1〜bin9にαを乗算することで算出する。bin1〜bin9は、それぞれ[−5、−2.5、−1、−0.5、0.5、1、2.5、5]に設定する。
B1 =α× bin1
・・・
B9 =α× bin9
In step S321, the width (B) of each bin of the two prediction error distributions is set based on α. The bin width is calculated by multiplying bin1 to bin9 by α as in the following equation. bin1 to bin9 are set to [-5, -2.5, -1, -0.5, 0.5, 1, 2.5, 5], respectively.
B1 = α × bin1
...
B9 = α × bin9

次に、ステップS322では、ステップS310で算出した現時点の予測誤差PEが、どこのビンの範囲内に該当するかどうか判定する。たとえば、αが0.64の場合、今回算出した操舵角予測誤差PEが0.2とすると、p5(−0.32〜0.32)のビンが範囲内と判定する。
そして、ステップS322の判定において、範囲外と判定した場合は、ステップS323に移行する。一方、範囲内と判定した場合はステップ324に移行する。そして、以降のステップS323、S324の処理で、現在の操舵角予測誤差分布の分布割合ps{i=1〜9}、普段の操舵角予測誤差分布の分布割合pl{i=1〜9}をそれぞれ更新する。ps、plは予測誤差分布のビンの割合を示す。
Next, in step S322, it is determined in which bin range the current prediction error PE calculated in step S310 falls. For example, when α is 0.64, if the currently calculated steering angle prediction error PE is 0.2, the bin of p5 (−0.32 to 0.32) is determined to be within the range.
If it is determined in step S322 that it is out of range, the process proceeds to step S323. On the other hand, if it is determined to be within the range, the process proceeds to step 324. Then, in the process of subsequent steps S323, S324, the distribution ratio of the current steering angle prediction error distribution ps i {i = 1~9}, the distribution ratio of the usual steering angle prediction error distribution pl i {i = 1~9 } Is updated respectively. ps i and p i i indicate the bin ratio of the prediction error distribution.

ここで、以下の処理では、実際には分布割合ps、plについて個別に更新処理を行うが、同一の処理であるため、代表として分布割合pi(:ps、pl)として説明する。
ステップS323では、下記式によって分布割合piを求める。すなわち、ビンの値を(1+1/N)で割る。つまり確率密度が減少する。その後ステップS325に移行する。
=p/(1+1/N)
ここで、Nは、予測誤差分布のタイムウインドウに応じた標本相当数(Ns、Nl)である。pは、予測誤差分布のタイムウインドウに応じた分布割合(ps{i=1〜9}、pl{i=1〜9})である。
Here, in the following process, in fact distribution ratio ps i, is performed individually update processing for pl i, since the same process, a representative as a distribution ratio pi (: ps i, pl i ) will be described as .
In step S323, the distribution ratio pi is obtained by the following equation. That is, the bin value is divided by (1 + 1 / N). That is, the probability density decreases. Thereafter, the process proceeds to step S325.
p i = p i / (1 + 1 / N)
Here, N is the sample equivalent number (Ns, Nl) corresponding to the time window of the prediction error distribution. p i is a distribution ratio (ps i {i = 1 to 9}, pl i {i = 1 to 9}) corresponding to the time window of the prediction error distribution.

また、ステップS324では、下記式によって分布割合piを求める。すなわち、ビンの値に(1/N)を足し、(1+1/N)で割る。つまり確率密度が増加する。その後ステップS325に移行する。
=(p+1/N)/(1+1/N)
ここで、Nは、予測誤差分布のタイムウインドウに応じた標本相当数(Ns、Nl)である。pは、予測誤差分布のタイムウインドウに応じた分布割合(ps{i=1〜9}、pl{i=1〜9})である。
In step S324, the distribution ratio pi is obtained by the following equation. That is, (1 / N) is added to the bin value and divided by (1 + 1 / N). That is, the probability density increases. Thereafter, the process proceeds to step S325.
p i = (p i + 1 / N) / (1 + 1 / N)
Here, N is the sample equivalent number (Ns, Nl) corresponding to the time window of the prediction error distribution. p i is a distribution ratio (ps i {i = 1 to 9}, pl i {i = 1 to 9}) corresponding to the time window of the prediction error distribution.

ステップS325では、ビンをひとつ進める。すなわち、iをインクリメントする。
ステップS326では、9つのビンが全て更新完了した(i=10)か否かを判定し、
9つのビンが全て更新完了すれば(i=10)、処理を終了して復帰する。一方、全てのビンを更新していない場合に、ステップS322に戻る。
以上で、不安定指標算出部の処理が終了する。
In step S325, the bin is advanced by one. That is, i is incremented.
In step S326, it is determined whether all nine bins have been updated (i = 10).
If all nine bins have been updated (i = 10), the process ends and returns. On the other hand, when all the bins have not been updated, the process returns to step S322.
This completes the process of the unstable index calculation unit.

次に、運転特性特徴量の算出処理について、図8を参照して説明する。
先ず、ステップS410で、普段の予測誤差分布における所定範囲内の予測誤差PEの大きさ(以下、PEdr)を求める。
所定範囲とは、ある分部内(百分位(%tile値)や正規化した範囲の所定値内である。例えば90%タイルを採用する。
次にステップS420では、運転者特性の推定するためのデータ量を満たしているか判定する。データ量というのは、短時間内で運転者特性の推定できるデータ量のことを指す。データ量が満たしていない場合は、パラメータの同定する必要がないため、運転者特性推定の処理を終了して、復帰する。一方、データ量が満たしている場合は、パラメータを同定するため、ステップS430に移行する。
Next, driving characteristic feature amount calculation processing will be described with reference to FIG.
First, in step S410, the size of the prediction error PE (hereinafter referred to as PEdr) within a predetermined range in the normal prediction error distribution is obtained.
The predetermined range is within a certain part (percentile (% tile value) or within a predetermined value of a normalized range. For example, 90% tile is adopted.
Next, in step S420, it is determined whether the data amount for estimating the driver characteristics is satisfied. The amount of data refers to the amount of data for which driver characteristics can be estimated within a short time. If the amount of data is not satisfied, it is not necessary to identify the parameter, so the driver characteristic estimation process is terminated and the process returns. On the other hand, when the data amount is satisfied, the process proceeds to step S430 in order to identify the parameter.

次にステップS430では、運転特性モデルの更新を行う。
ここで、上述のように基準とする基本運転特性モデルを予め設定している。この基準とする基本運転特性モデルは、予測誤差分布の状態で設定してあり、実験などによって、予測誤差分布の度合いを設定しておく。具体的には、基本運転特性モデルの予測誤差分布の情報として、上記所定範囲内のPEの大きさ(以下、PEbasic)と、予測誤差分布の横軸(BIN{i =1〜9})と縦軸の割合(DIST{i=1〜9})とからなる。なお、ビン個数はここでは9個設定しているが、この値に限定するものではない。
そして、ステップS440では、この基本運転特性モデルのPEbasicと、ステップS410で算出したPEdrとから、この両者の比を求める。更に、その比に予測誤差分布の横軸BINに掛け合わせて、下記式のように、学習した運転特性モデルの予測誤差分布pl=BIN_MODEL{i =1〜9}を求める。
すなわち、
BIN_MODEL =(PEdr/PEbasic)×BIN
このステップにより、学習した運転特性モデルの予測誤差分布plは、
[x、y]=[BIN_MODEL、 DIST
と表現することが出来る。
Next, in step S430, the driving characteristic model is updated.
Here, the basic operation characteristic model used as a reference is set in advance as described above. The basic operating characteristic model used as a reference is set in the state of the prediction error distribution, and the degree of the prediction error distribution is set through experiments or the like. Specifically, as information on the prediction error distribution of the basic operating characteristic model, the PE size within the predetermined range (hereinafter referred to as PE basic) and the horizontal axis of the prediction error distribution (BIN i {i = 1 to 9}) And the ratio of the vertical axis (DIST i {i = 1 to 9}). Here, the number of bins is set to nine here, but is not limited to this value.
In step S440, a ratio between the basic operating characteristic model PEbasic and the PEdr calculated in step S410 is obtained. Further, by multiplying the horizontal axis BIN prediction error distribution in the ratio, as shown in Equation obtains a prediction error distribution of the learned driving characteristics model pl i = BIN_MODEL i {i = 1~9}.
That is,
BIN_MODE i = (PEdr / PEbasic) × BIN i
Through this step, the predicted error distribution pl i of the learned driving characteristic model is
[X, y] = [BIN_MODE i , DIST i ]
Can be expressed.

(動作・作用)
従来技術のような技術では、長時間学習期間(例えば4〜5日)が必要であり、その学習期間が終了してから情報提示が可能となる。これに対し、本実施形態での運転状態の学習では、運転操作の誤差分布を使用することにより、短期間で所定精度の運転特性を収集出来る。すなわち、運転者の特性を推定、推定結果に基づき運転者の特性に近いパラメータの同定を行うことで、短時間で情報提示が可能となる。
(Operation / Action)
A technique such as the prior art requires a long learning period (for example, 4 to 5 days), and information can be presented after the learning period ends. On the other hand, in the learning of the driving state in the present embodiment, it is possible to collect driving characteristics with a predetermined accuracy in a short period of time by using the error distribution of the driving operation. That is, information can be presented in a short time by estimating the driver characteristics and identifying parameters close to the driver characteristics based on the estimation results.

ここで、各運転者の「普段(長時間)」と、「現在このとき(直近)」の運転行動指標の比較により、普段と違う運転状態を検出し、情報提示を行っている。またステアリングエントロピー法を用いて、「普段」と「現在このとき」の2つの予測誤差(PE)分布から、ふらつき運転指標を算出し、ふらつき運転状態を検出している。
一般に、「普段」の運転行動データを取得するためには、長時間のデータが必要である。しかしながら、長時間のデータが取得できるまで、効果的な情報提示をすることは困難である。長時間のデータが必要な理由は、長時間のデータはノイズの影響が小さいためである。
Here, by comparing the driving behavior indexes of “normal (long time)” and “currently (current)” of each driver, a different driving state is detected and information is presented. Further, the steering entropy method is used to calculate a wobbling driving index from two prediction error (PE) distributions of “normal” and “currently at this time”, and detect the wobbling driving state.
In general, long-term data is required to obtain “normal” driving behavior data. However, it is difficult to present information effectively until long-term data can be acquired. The reason why long time data is required is that long time data is less affected by noise.

ここで、本発明者らが実験などによって調査、検討したところ、運転操作の予測誤差分布について、ノイズの影響はある幅の外側に影響し、幅の内側ではノイズの影響を受けにくいことを見出した。
更に、ノイズの影響が少ない(長時間のデータの)PE分布における、被験者違いや走行環境違いなどの条件違いの差異は、ある特定の幅が異なるだけで分布のプロポーションは変わらないことを見出した。
以上のことから、ある代表値で正規化することで、ほぼ同じ分布プロポーションが導けることから、分布の代表値がわかれば、その分布自体を推定することができる。
Here, as a result of investigation and examination by the present inventors through experiments, etc., it was found that the influence of noise affects the outside of a certain width of the prediction error distribution of driving operation, and that the influence of noise is hardly affected by the inside of the width. It was.
Furthermore, we found that the difference in conditions such as subject differences and driving environment differences in PE distribution with little noise influence (for long-term data) does not change the distribution proportions only by a certain specific width. .
From the above, by normalizing with a certain representative value, almost the same distribution proportion can be derived. Therefore, if the representative value of the distribution is known, the distribution itself can be estimated.

上記「ノイズの影響はある幅の外側に影響し、幅の内側ではノイズの影響を受けにくい」ことについて説明する。
図9と図10にある被験者の10分間の走行データのPEと5トリップ分(およそ450分)のPE分布をそれぞれ示す。この図から分かるように、中央値から離れた値に影響を受けやすい標準偏差は両者で異なるものの、幅(90%tile幅:5〜95%tile)で見ると両者の差異は小さいことがわかる。
この特性を利用し、ある特定の幅の大きさを算出し、幅の大きさが同等なモデルを構築することで、ノイズの影響が少ない予測誤差分布を得ることができることが分かる。
The above explanation will be made that “the influence of noise affects the outside of a certain width and is less susceptible to noise inside the width”.
FIG. 9 and FIG. 10 show the PE distribution of the subject's 10-minute running data and the PE distribution for 5 trips (approximately 450 minutes), respectively. As can be seen from this figure, although the standard deviations that are easily affected by values away from the median value are different, the difference between the two is small when viewed in width (90% tile width: 5 to 95% tile). .
It can be seen that a prediction error distribution with little influence of noise can be obtained by calculating the size of a specific width using this characteristic and constructing a model having the same width.

次に、上記「被験者違いや走行環境違いなどの条件違いの差異は、ある特定の幅が異なるだけで分布のプロポーションは変わらないこと」について説明する。
図11に正規化したPEの頻度分布と普段のPEの頻度分布を示す。この図11は、被験者12名のPEの分布をそれぞれ示したものである(複数の被験者の分布を重ね合わせて表示している)。上図は、ある分布の代表値(90%tile)で正規化したPEの頻度分布であり、下図は普段の頻度分布である。この図5から分かるように、正規化した頻度分布のほうが、縦のばらつきが小さい(被験者間の差が小さい)ことが分かる。さらに各被験者のプロポーションに大きな差がないことが分かる。
Next, explanation will be made on the above-mentioned “difference in difference in conditions such as differences in subjects and driving environments is that the proportion of distribution does not change only by a certain specific width”.
FIG. 11 shows the normalized PE frequency distribution and the normal PE frequency distribution. FIG. 11 shows the distribution of PEs of 12 subjects (distributed distributions of a plurality of subjects are displayed). The upper figure shows the frequency distribution of PE normalized by a representative value (90% tile) of a certain distribution, and the lower figure shows the usual frequency distribution. As can be seen from FIG. 5, it can be seen that the normalized frequency distribution has a smaller vertical variation (difference between subjects is smaller). Furthermore, it turns out that there is no big difference in the proportion of each subject.

この特性を利用し、特定の幅がわかれば異なる分布を推定することができる。
すなわち、ある特定の幅が求められれば、ノイズの少ない分布を推定することが可能である。
そして本実施形態では、運転操作の誤差分布を使用することで、学習期間が短くなると共に、予め初期値としての分布のモデルを持つことで、さらに学習期間を短くすることが可能となる。
ここで、操舵角センサ1は運転操作情報取得手段を構成する。ステップS300は、運転特性特徴量算出手段2Abを構成する。ステップS400は、モデル取得手段2Acを構成する。ステップS600,S700は不安定運転状態判定手段2Adを構成する。
Using this characteristic, different distributions can be estimated if a specific width is known.
That is, if a specific width is obtained, a distribution with less noise can be estimated.
In this embodiment, the learning period is shortened by using the error distribution of the driving operation, and the learning period can be further shortened by having a distribution model as an initial value in advance.
Here, the steering angle sensor 1 constitutes driving operation information acquisition means. Step S300 constitutes the driving characteristic feature amount calculating means 2Ab. Step S400 constitutes model acquisition means 2Ac. Steps S600 and S700 constitute unstable operation state determination means 2Ad.

(本実施形態の効果)
(1)運転操作情報取得手段は、運転者による運転操作の情報を取得する。運転特性特徴量算出手段2Abは、運転操作情報取得手段が取得する運転操作情報に基づき、運転者の運転特性の特徴量を運転操作の予測誤差分布の情報として算出する。モデル取得手段2Acは、上記運転特性特徴量算出手段2Abが算出した運転特性の特徴量、及び、予め設定した基準となる運転操作の予測誤差分布の情報からなる運転特性モデルに基づき、学習後の運転特性モデルを取得する。不安定運転状態判定手段2Adは、モデル取得手段2Acが取得した学習後の運転特性モデルと、上記運転特性特徴量算出手段2Abが算出した現在の特徴量とに基づき、現在の運転状態が普段の運転特性とは異なる不安定運転状態か否かを判定する。
(Effect of this embodiment)
(1) The driving operation information acquisition means acquires driving operation information by the driver. The driving characteristic feature amount calculating unit 2Ab calculates the characteristic amount of the driving characteristic of the driver as information on the prediction error distribution of the driving operation based on the driving operation information acquired by the driving operation information acquiring unit. The model acquisition unit 2Ac is based on the driving characteristic model calculated from the driving characteristic feature amount calculated by the driving characteristic feature amount calculating unit 2Ab and the driving error prediction error distribution information serving as a preset reference. Obtain an operating characteristic model. The unstable driving state determination means 2Ad is based on the learned driving characteristic model acquired by the model acquisition means 2Ac and the current characteristic amount calculated by the driving characteristic feature amount calculation means 2Ab. It is determined whether or not an unstable operation state different from the operation characteristics.

運転者の運転特性の特徴量として運転操作の予測誤差分布の情報を使用する。このため、学習時間を短く出来る。また、予め基準となる運転操作の予測誤差分布の情報からなる運転特性モデルを備えることで、更に学習時間を短く出来る。
以上の事から、短時間に得られる少ないデータで、運転者の運転状態の判定が可能となる。
Information on the prediction error distribution of the driving operation is used as a feature amount of the driving characteristic of the driver. For this reason, learning time can be shortened. In addition, the learning time can be further shortened by providing a driving characteristic model including information on the prediction error distribution of the driving operation as a reference in advance.
From the above, it is possible to determine the driving state of the driver with a small amount of data obtained in a short time.

(2)上記運転特性の特徴量は、運転操作の操舵角予測誤差分布に関する情報である。
運転者の操舵操作による運転状態を判定することが可能となる。
(3)上記モデル取得手段2Acは、運転操作の予測誤差分布における、所定の累積頻度となる予測誤差を運転特性の特徴量として、学習後の運転特性モデルを取得する。
所定の累積頻度となる予測誤差を運転特性の特徴量として使用することで、より短時間で学習時間を行うことが出来る。
(2) The characteristic amount of the driving characteristic is information related to the steering angle prediction error distribution of the driving operation.
It becomes possible to determine the driving state by the steering operation of the driver.
(3) The model acquisition unit 2Ac acquires a learned driving characteristic model using a prediction error having a predetermined cumulative frequency in the prediction error distribution of the driving operation as a characteristic amount of the driving characteristic.
By using a prediction error having a predetermined cumulative frequency as a feature amount of driving characteristics, learning time can be performed in a shorter time.

(変形例)
(1)上記実施形態では、ステップS400にて、運転特性の特徴量を算出する際に、所定範囲内のPEの大きさを算出(ステップS410)している。これに代えて、運転特性特徴量を算出する際に、所定のPEにおける累積頻度を算出しても良い。この場合、ステップS410にて所定PEの累積頻度を算出し、その累積頻度によって、運転特性特徴量を算出するようにしても良い。
すなわち、上記モデル取得手段2Acは、運転操作の予測誤差分布における、所定の予測誤差の値における累積頻度を運転特性の特徴量として、学習後の運転特性モデルを取得する。
効果は、上記実施形態と同様で、短時間で運転状態の判定が可能となる。
(2)運転操作は、アクセルペダルやブレーキペダルなどの操作子の操作情報であっても良い。
(Modification)
(1) In the above embodiment, when calculating the characteristic value of the driving characteristic in step S400, the size of the PE within the predetermined range is calculated (step S410). Instead of this, the cumulative frequency in a predetermined PE may be calculated when calculating the driving characteristic feature amount. In this case, the cumulative frequency of the predetermined PE may be calculated in step S410, and the driving characteristic feature amount may be calculated based on the cumulative frequency.
That is, the model acquisition unit 2Ac acquires the learned driving characteristic model using the cumulative frequency of the predetermined prediction error value in the driving operation prediction error distribution as the characteristic amount of the driving characteristic.
The effect is the same as in the above embodiment, and the operation state can be determined in a short time.
(2) The driving operation may be operation information of an operator such as an accelerator pedal or a brake pedal.

(第2実施形態)
次に、第2実施形態について図面を参照しつつ説明する。なお、上記第1実施形態と同様な構成については同一の符号を付して説明する。
本発明の基本構成は上記実施形態と同様である。
但し、運転者推定を行う条件を限定し、より精度の高い運転特性特徴量算出を狙う。
また、車速センサ430は、車輪の回転数や変速機の出力側の回転数を計測することにより自車両の車速を検出し、検出した自車速を情報提示コントローラ2に出力する。
(Second Embodiment)
Next, a second embodiment will be described with reference to the drawings. In addition, the same code | symbol is attached | subjected and demonstrated about the structure similar to the said 1st Embodiment.
The basic configuration of the present invention is the same as that of the above embodiment.
However, the conditions for performing the driver estimation are limited and aiming for more accurate driving characteristic feature amount calculation.
The vehicle speed sensor 430 detects the vehicle speed of the host vehicle by measuring the number of rotations of the wheels and the number of rotations on the output side of the transmission, and outputs the detected host vehicle speed to the information presentation controller 2.

本実施形態の情報提示コントローラ2の処理を、図12を参照して説明する。
本実施形態では、ステップS300の処理とステップS400の処理の間にステップS1000を追加した。このステップS1000以外の他の処理は上記第1実施形態と同様であるので、説明を省略する。
ステップS1000では、特徴量算出条件を満足するか否かを判定する。特徴量算出条件を満足する場合にはステップS400に移行する。特徴量算出条件を満足しない場合には、処理を終了して復帰する。
Processing of the information presentation controller 2 of the present embodiment will be described with reference to FIG.
In the present embodiment, step S1000 is added between the process of step S300 and the process of step S400. Since other processes other than step S1000 are the same as those in the first embodiment, description thereof will be omitted.
In step S1000, it is determined whether or not a feature amount calculation condition is satisfied. If the feature amount calculation condition is satisfied, the process proceeds to step S400. If the feature amount calculation condition is not satisfied, the process ends and returns.

特徴量算出条件は、運転者特性の推定で所定以上の精度を見込める条件である。特徴量算出条件は、例えば車速40km/h以上、かつ、走行道路のカーブ割合(後述)が10%以下の場合である。カーブ割合は、例えば、操舵角が所定値(例えば12deg)以上の場合に、カーブを走行していると判定し、その時間割合を積算して求める。
なお、使用するデータは、ステップS200で記録した運転行動データ記録を参照して実施すれば良い。
The feature amount calculation condition is a condition that allows the driver characteristics to be estimated with a predetermined accuracy or higher. The feature amount calculation condition is, for example, a case where the vehicle speed is 40 km / h or more and the curve ratio (described later) of the traveling road is 10% or less. The curve ratio is obtained, for example, by determining that the vehicle is traveling on a curve when the steering angle is a predetermined value (for example, 12 degrees) or more and integrating the time ratio.
In addition, what is necessary is just to implement the data to be used with reference to the driving action data record recorded in step S200.

(動作・作用)
操舵誤差は、道路線形や運転者の負荷の状態に影響を受けるため、運転者が無負荷の状態で、学習のための操舵誤差分布を計測することが好ましい。そこで、これらの影響を受けにくくするために、直線路とみなせる道路線形走行時の操舵角データに基づいて計測した操舵角データに基づいてモデルの操舵誤差分布を更新して学習する。
ここで、車速センサ4は走行情報取得手段を構成する。
(Operation / Action)
Since the steering error is affected by the road alignment and the load state of the driver, it is preferable to measure the steering error distribution for learning while the driver is in no load state. Therefore, in order to make it less susceptible to these effects, the steering error distribution of the model is updated and learned based on the steering angle data measured based on the steering angle data during road linear traveling that can be regarded as a straight road.
Here, the vehicle speed sensor 4 constitutes a travel information acquisition unit.

(本実施形態の効果)
(1)走行情報取得手段が車両の走行情報を取得する。モデル取得手段2Acは、上記運転特性特徴量算出手段2Abが算出する運転特性の特徴量のうち、走行情報取得手段が取得する車両の走行情報に基づき所定の走行状態と判定したときの運転特性の特徴量を使用して、学習後の運転特性モデルを取得する。
モデルの学習の際に道路線形や運転者の負荷の状態に影響を受け難いデータを使用することで、モデルの精度が向上する。
(Effect of this embodiment)
(1) The travel information acquisition means acquires the travel information of the vehicle. The model acquisition means 2Ac is the driving characteristic when the driving characteristic characteristic amount calculated by the driving characteristic characteristic calculation means 2Ab is determined to be a predetermined driving state based on the vehicle driving information acquired by the driving information acquisition means. A feature model is used to obtain a driving characteristic model after learning.
When learning the model, the accuracy of the model is improved by using data that is not easily affected by the road alignment and the load state of the driver.

(第3実施形態)
次に、第3実施形態について図面を参照しつつ説明する。なお、上記第各実施形態と同様な構成については同一の符号を付して説明する。
本発明の基本構成は上記実施形態と同様である。
但し、運転特性特徴量の算出処理の一部が異なる。
すなわち、運転特性特徴量を算出した運転状態属性から、特徴量を補正し、より高い運転者特性推定を狙う。
また、車速センサ430は、車輪の回転数や変速機の出力側の回転数を計測することにより自車両の車速を検出し、検出した自車速を情報提示コントローラ2に出力する。
(Third embodiment)
Next, a third embodiment will be described with reference to the drawings. In addition, about the structure similar to said each embodiment, the same code | symbol is attached | subjected and demonstrated.
The basic configuration of the present invention is the same as that of the above embodiment.
However, a part of the calculation process of the driving characteristic feature amount is different.
That is, the feature amount is corrected from the driving state attribute for which the driving characteristic feature amount is calculated, and higher driver characteristic estimation is aimed.
The vehicle speed sensor 430 detects the vehicle speed of the host vehicle by measuring the number of rotations of the wheels and the number of rotations on the output side of the transmission, and outputs the detected host vehicle speed to the information presentation controller 2.

本実施形態の運転特性特徴量の算出処理について、図13を参照して説明する。ステップS415及びステップS425の処理を追加した点が異なる。
先ずステップS410で、所定範囲内のPEの大きさ(以下、PEdr)を求める。
次にステップS415では、走行状態属性の記録を行う。
属性として記録するのは、例えば
・カーブ走行中か直線走行か
・走行道路の種別
などを記録する。上記走行状態の属性は、運転者に普段要求される操舵量の大きさに関する属性である。
The calculation process of the driving characteristic feature amount of the present embodiment will be described with reference to FIG. The difference is that the processing of step S415 and step S425 is added.
First, in step S410, the size of PE within a predetermined range (hereinafter referred to as PEdr) is obtained.
Next, in step S415, the running state attribute is recorded.
The attributes are recorded, for example: • Whether the vehicle is traveling on a curve or traveling in a straight line • The type of traveling road is recorded. The attribute of the traveling state is an attribute related to the magnitude of the steering amount that is normally required by the driver.

なお、この走行状態の種別は、ナビ情報などから判定すれば良い。なお、ナビゲーションシステムは、GPS受信機、地図データベース、および表示モニタ等を備えており、経路探索および経路案内等を行うシステムである。ナビゲーションシステムは、GPS受信機から得られる自車両の現在位置と地図データベースに格納された道路情報に基づいて、自車両が走行する道路の種別や道路幅員等の情報を取得することができる。   The type of running state may be determined from navigation information or the like. The navigation system includes a GPS receiver, a map database, a display monitor, and the like, and performs route search, route guidance, and the like. The navigation system can acquire information such as the type of road on which the host vehicle is running and the road width based on the current position of the host vehicle obtained from the GPS receiver and the road information stored in the map database.

次にステップS420では、運転者特性の推定するためのデータ量を満たしているか判定する。データ量というのは、短時間内で運転者特性の推定できるデータ量のことを指す。データ量が満たしていない場合は、パラメータの同定する必要がないため、運転者特性推定の処理を終了して、復帰する。一方、データ量が満たしている場合は、パラメータを同定するため、ステップS430に移行する。   Next, in step S420, it is determined whether the data amount for estimating the driver characteristics is satisfied. The amount of data refers to the amount of data for which driver characteristics can be estimated within a short time. If the amount of data is not satisfied, it is not necessary to identify the parameter, so the driver characteristic estimation process is terminated and the process returns. On the other hand, when the data amount is satisfied, the process proceeds to step S430 in order to identify the parameter.

ステップS430では、ステップS410で算出したPEdrを、下記式に基づき補正する。
PEdr_modify = PEdr + ΣPE_mod
また、ΣPE_mod = Curve_mod + RoadIDX_mod
ここで、Curve_modは、カーブの割合を示し、対象とする学習時期におけるカーブの割合に基づき、図14のようなマップを使用して求める。すなわち、カーブの割合が大きいほど小さくなる。RoadIDX_modは、高速道路の割合を示し、対象とする学習時期における高速道路の割合に基づき、図15のようなマップを使用して求める。すなわち、高速道路の割合が大きいほど大きくなる
次にステップS430では、モデルの作成をする。
次にステップS440では、モデルを初期値として適合させる。モデルとして作った予測誤差分布を初期値として代入し、パラメータの同定完了とする。
In step S430, PEdr calculated in step S410 is corrected based on the following equation.
PEdr_modify = PEdr + ΣPE_mod
Also, ΣPE_mod = Curve_mod + LoadIDX_mod
Here, Curve_mod indicates the ratio of the curve, and is obtained using a map as shown in FIG. 14 based on the ratio of the curve at the target learning period. That is, the larger the curve ratio, the smaller. RoadIDX_mod indicates the ratio of the expressway, and is obtained using a map as shown in FIG. 15 based on the ratio of the expressway in the target learning period. That is, the larger the ratio of highways, the larger the model. Next, in step S430, a model is created.
Next, in step S440, the model is adapted as an initial value. The prediction error distribution created as a model is substituted as an initial value, and parameter identification is completed.

(動作・作用)
操舵誤差は、道路線形や運転者の負荷の状態に影響を受ける。この影響を受けた分だけモデルを特定するための予測誤差分布を補正する。
ここで、車速センサ4は走行情報取得手段を構成する。ステップS415,S425は運転特性特徴量補正手段を構成する。
(Operation / Action)
Steering error is affected by road alignment and driver load conditions. The prediction error distribution for specifying the model is corrected by the amount affected by this influence.
Here, the vehicle speed sensor 4 constitutes a travel information acquisition unit. Steps S415 and S425 constitute driving characteristic feature amount correcting means.

(本実施形態の効果)
(1)特性特徴量補正手段は、走行状態の属性に基づき運転特性の特徴量を補正する。 上記走行状態の属性は、普段要求される運転者の操作量の大きさに関する属性である。
これによってモデルの精度が向上する。
(2)上記走行状態の属性は、道路種別あるいは道路形状である。
これによって、操舵状態に影響のある属性を使用可能となる。
(Effect of this embodiment)
(1) The characteristic feature amount correcting unit corrects the characteristic amount of the driving characteristic based on the driving state attribute. The attribute of the driving state is an attribute related to the magnitude of the operation amount of the driver that is normally requested.
This improves the accuracy of the model.
(2) The attribute of the traveling state is a road type or a road shape.
As a result, attributes that affect the steering state can be used.

(3)上記走行状態の属性は、道路形状がカーブであるか否かである。上記運転特性特徴量補正手段は、属性が所定以上のカーブ路である時間割合に応じて補正する。
これによって、モデルへの、カーブ走行による影響を小さく出来る。
(4)上記走行状態の属性は、道路種別である。上記運転特性特徴量補正手段は、道路種別ごとの時間割合を用いて補正する。
これによって、モデルへの、道路種別による影響を小さく出来る。
(3) The attribute of the traveling state is whether or not the road shape is a curve. The driving characteristic feature amount correcting unit corrects the driving characteristic feature amount according to a time ratio in which the attribute is a curve road having a predetermined value or more.
As a result, the influence of the curve running on the model can be reduced.
(4) The attribute of the traveling state is a road type. The driving characteristic feature amount correcting unit corrects the driving characteristic feature amount using a time ratio for each road type.
As a result, the influence of the road type on the model can be reduced.

(変形例)
(1)取得する走行状態は、上記走行状態に限定しない。要求させる操舵量に関係する走行状態であれば、使用可能である。
(2)上記説明では、道路種別として高速道路を例示したが、砂利道など他の道路種別を使用しても良い。対象とする道路の走行による操舵に対する負荷量に応じて補正量を求めればよい。
(Modification)
(1) The travel state to be acquired is not limited to the travel state. Any traveling state related to the requested steering amount can be used.
(2) In the above description, a highway is exemplified as the road type, but other road types such as a gravel road may be used. What is necessary is just to obtain | require a correction amount according to the load amount with respect to the steering by driving | running | working of the object road.

1 操舵角センサ(運転操作情報取得手段)
2 情報提示コントローラ
2A 運転状態判定装置
2Aa 運転特性モデル
2Ab 運転特性特徴量算出手段
2Ac モデル取得手段
2Ad 不安定運転状態判定手段
2B 情報発信手段
3 情報提示手段
4 車速センサ(走行情報取得手段)
PE 予測誤差
分布割合
pl 分布割合
ps 分布割合
R_Hpls 相対エントロピー
θ 操舵角
θp 予測操舵角
1 Steering angle sensor (Driving operation information acquisition means)
2 Information presentation controller 2A Driving state determination device 2Aa Driving characteristic model 2Ab Driving characteristic feature amount calculation means 2Ac Model acquisition means 2Ad Unstable driving state determination means 2B Information transmission means 3 Information presentation means 4 Vehicle speed sensor (running information acquisition means)
PE prediction error p i distribution ratio pl i distribution ratio ps i distribution ratio R_Hpls relative entropy θ steering angle θp prediction steering angle

Claims (10)

運転者による運転操作の情報を取得する運転操作情報取得手段と、
運転操作情報取得手段が取得する運転操作情報に基づき、運転者の運転特性の特徴量を運転操作の予測誤差分布の情報として算出する運転特性特徴量算出手段と、
上記運転特性特徴量算出手段が算出した運転特性の特徴量、及び、予め設定した基準となる運転操作の予測誤差分布の情報からなる運転特性モデルに基づき、学習後の運転特性モデルを取得するモデル取得手段と、
モデル取得手段が取得した学習後の運転特性モデルと、上記運転特性特徴量算出手段が算出した現在の特徴量とに基づき、現在の運転状態が普段の運転特性とは異なる不安定運転状態か否かを判定する不安定運転状態判定手段と、
を備えることを特徴とする運転状態判定装置。
Driving operation information acquisition means for acquiring driving operation information by the driver;
Based on the driving operation information acquired by the driving operation information acquiring means, driving characteristic feature amount calculating means for calculating the characteristic amount of the driving characteristic of the driver as information of the prediction error distribution of the driving operation;
A model for acquiring a learned driving characteristic model based on the driving characteristic feature amount calculated by the driving characteristic feature amount calculating means and the driving characteristic model information including the prediction error distribution information of the driving operation as a reference set in advance. Acquisition means;
Whether the current driving state is an unstable driving state different from the normal driving characteristic based on the driving characteristic model after learning acquired by the model acquiring unit and the current characteristic amount calculated by the driving characteristic feature amount calculating unit. An unstable operation state determination means for determining whether or not
The driving | running state determination apparatus characterized by the above-mentioned.
上記運転特性の特徴量は、運転操作の操舵角予測誤差分布に関する情報であることを特徴とする請求項1に記載した運転状態判定装置。   The driving state determination device according to claim 1, wherein the characteristic amount of the driving characteristic is information related to a steering angle prediction error distribution of the driving operation. 上記モデル取得手段は、運転操作の予測誤差分布における、所定の累積頻度となる予測誤差を運転特性の特徴量として、学習後の運転特性モデルを取得することを特徴とする請求項1又は請求項2に記載した運転状態判定装置。   The said model acquisition means acquires the driving | running characteristic model after learning by using the prediction error used as predetermined | prescribed accumulation frequency in the prediction error distribution of driving | operation operation as the feature-value of driving | running characteristic. The operating state determination apparatus described in 2. 上記モデル取得手段は、運転操作の予測誤差分布における、所定の予測誤差の値における累積頻度を運転特性の特徴量として、学習後の運転特性モデルを取得することを特徴とする請求項1又は請求項2に記載した運転状態判定装置。   The said model acquisition means acquires the driving characteristic model after learning by using the cumulative frequency in the value of the predetermined prediction error in the prediction error distribution of the driving operation as the characteristic amount of the driving characteristic. Item 3. The operating state determination device according to Item 2. 車両の走行情報を取得する走行情報取得手段を備え、
モデル取得手段は、上記運転特性特徴量算出手段が算出する運転特性の特徴量のうち、走行情報取得手段が取得する車両の走行情報に基づき所定の走行状態と判定したときの運転特性の特徴量を使用して、学習後の運転特性モデルを取得することを特徴とする請求項1〜請求項4のいずれか1項に記載した運転状態判定装置。
Comprising travel information acquisition means for acquiring travel information of the vehicle;
The model acquisition means is a characteristic value of the driving characteristic when the driving characteristic feature quantity calculated by the driving characteristic feature quantity calculation means is determined to be a predetermined driving state based on the vehicle driving information acquired by the driving information acquisition means. The driving | running state determination apparatus described in any one of Claims 1-4 which acquires the driving | running characteristic model after learning using.
走行情報取得手段が取得する車両の走行情報に基づき、走行状態の属性を取得する走行状態属性取得手段を備え、
上記モデル取得手段は、走行状態取得手段が取得した走行状態の属性に基づき運転特性の特徴量を補正する運転特性特徴量補正手段を備え、
上記走行状態の属性は、普段要求される運転者の操作量の大きさに関する属性であることを特徴とする請求項1〜請求項5のいずれか1項に記載した運転状態判定装置。
Based on the vehicle travel information acquired by the travel information acquisition means, comprising a travel state attribute acquisition means for acquiring the attribute of the travel state,
The model acquisition unit includes a driving characteristic feature amount correction unit that corrects the characteristic amount of the driving characteristic based on the attribute of the driving state acquired by the driving state acquisition unit,
The driving state determination device according to any one of claims 1 to 5, wherein the attribute of the traveling state is an attribute relating to a magnitude of a driver's operation amount that is normally requested.
上記走行状態の属性は、道路種別あるいは道路形状であることを特徴とする請求項6に記載した運転状態判定装置。   The driving state determination device according to claim 6, wherein the attribute of the driving state is a road type or a road shape. 上記走行状態の属性は、道路形状がカーブであるか否かであり、
上記運転特性特徴量補正手段は、属性が所定以上のカーブ路である時間割合に応じて補正することを特徴とする請求項7に記載した運転状態判定装置。
The attribute of the driving state is whether or not the road shape is a curve,
8. The driving state determination apparatus according to claim 7, wherein the driving characteristic feature amount correcting unit corrects the driving characteristic feature amount according to a time ratio of an attribute having a curve road having a predetermined value or more.
上記走行状態の属性は、道路種別であり、
上記運転特性特徴量補正手段は、道路種別ごとの時間割合を用いて補正することを特徴とする請求項7又は請求項8に記載した運転状態判定装置。
The attribute of the above running state is a road type,
The driving state determination device according to claim 7 or 8, wherein the driving characteristic feature amount correcting unit corrects the driving characteristic feature amount using a time ratio for each road type.
運転者の運転操作情報に基づき、運転者の運転特性の特徴量を運転操作の予測誤差分布の情報として算出し、算出した運転特性の特徴量、及び、予め設定した基準となる運転操作の予測誤差分布の情報からなる運転特性モデルに基づき、学習後の運転特性モデルを取得して、
取得した学習後の運転特性モデルと、現在の運転特性の特徴量とに基づき、現在の運転状態が普段の運転特性とは異なる不安定運転状態か否かを判定することを特徴とする運転状態判定方法。
Based on the driving operation information of the driver, the characteristic amount of the driving characteristic of the driver is calculated as information of the prediction error distribution of the driving operation, and the calculated characteristic value of the driving characteristic and the prediction of the driving operation as a preset reference are calculated. Based on the driving characteristic model consisting of error distribution information, obtain the driving characteristic model after learning,
A driving state characterized by determining whether or not the current driving state is an unstable driving state different from the normal driving characteristic based on the acquired driving characteristic model after learning and the characteristic amount of the current driving characteristic. Judgment method.
JP2009039631A 2009-02-23 2009-02-23 Driving state determination device and driving state determination method Active JP5267203B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2009039631A JP5267203B2 (en) 2009-02-23 2009-02-23 Driving state determination device and driving state determination method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009039631A JP5267203B2 (en) 2009-02-23 2009-02-23 Driving state determination device and driving state determination method

Publications (2)

Publication Number Publication Date
JP2010198120A true JP2010198120A (en) 2010-09-09
JP5267203B2 JP5267203B2 (en) 2013-08-21

Family

ID=42822818

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009039631A Active JP5267203B2 (en) 2009-02-23 2009-02-23 Driving state determination device and driving state determination method

Country Status (1)

Country Link
JP (1) JP5267203B2 (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012144131A1 (en) * 2011-04-20 2012-10-26 日産自動車株式会社 Information provision device for use in vehicle
WO2013183117A1 (en) * 2012-06-05 2013-12-12 トヨタ自動車株式会社 Driving characteristics estimation device and driver assistance system
KR20140002373A (en) * 2012-06-29 2014-01-08 현대자동차주식회사 Apparatus and method for monitoring driver state by driving pattern learning
JP2015143936A (en) * 2014-01-31 2015-08-06 三菱電機株式会社 Vehicle driving support device, server, vehicle driving support system, and vehicle driving support program
JP2017511916A (en) * 2014-02-07 2017-04-27 ルノー エス.ア.エス. A method for diagnosing the awakening level of a vehicle driver
JP2017126253A (en) * 2016-01-15 2017-07-20 トヨタ自動車株式会社 Operating load estimation device
WO2017203694A1 (en) * 2016-05-27 2017-11-30 日産自動車株式会社 Driving control method and driving control device
JP2019159147A (en) * 2018-03-14 2019-09-19 クラリオン株式会社 On-vehicle device, traveling state estimation method, server device, information processing method, and traveling state estimation system
JP2019156355A (en) * 2018-03-16 2019-09-19 本田技研工業株式会社 Vehicle control device
JP2021051688A (en) * 2019-09-26 2021-04-01 トヨタ自動車株式会社 Drive skill evaluation device

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007272834A (en) * 2006-03-31 2007-10-18 Equos Research Co Ltd Driver model creation device, driving support device, and driving action determination device
JP2009009495A (en) * 2007-06-29 2009-01-15 Nissan Motor Co Ltd Vehicle driving support device and vehicle including vehicle driving support device

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007272834A (en) * 2006-03-31 2007-10-18 Equos Research Co Ltd Driver model creation device, driving support device, and driving action determination device
JP2009009495A (en) * 2007-06-29 2009-01-15 Nissan Motor Co Ltd Vehicle driving support device and vehicle including vehicle driving support device

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012226602A (en) * 2011-04-20 2012-11-15 Nissan Motor Co Ltd Vehicle information providing-device
RU2549578C1 (en) * 2011-04-20 2015-04-27 Ниссан Мотор Ко., Лтд. Information representation device to be used in vehicle
WO2012144131A1 (en) * 2011-04-20 2012-10-26 日産自動車株式会社 Information provision device for use in vehicle
US9256989B2 (en) 2011-04-20 2016-02-09 Nissan Motor Co., Ltd. Information provision device for use in vehicle
US9428057B2 (en) 2011-04-20 2016-08-30 Nissan Motor Co., Ltd. Information provision device for use in vehicle
CN104334431B (en) * 2012-06-05 2017-09-29 丰田自动车株式会社 Driving performance estimating device and driving assist system
WO2013183117A1 (en) * 2012-06-05 2013-12-12 トヨタ自動車株式会社 Driving characteristics estimation device and driver assistance system
CN104334431A (en) * 2012-06-05 2015-02-04 丰田自动车株式会社 Driving characteristics estimation device and driver assistance system
JPWO2013183117A1 (en) * 2012-06-05 2016-01-21 トヨタ自動車株式会社 Driving characteristic estimation device and driving support system
KR20140002373A (en) * 2012-06-29 2014-01-08 현대자동차주식회사 Apparatus and method for monitoring driver state by driving pattern learning
JP2015143936A (en) * 2014-01-31 2015-08-06 三菱電機株式会社 Vehicle driving support device, server, vehicle driving support system, and vehicle driving support program
JP2017511916A (en) * 2014-02-07 2017-04-27 ルノー エス.ア.エス. A method for diagnosing the awakening level of a vehicle driver
JP2017126253A (en) * 2016-01-15 2017-07-20 トヨタ自動車株式会社 Operating load estimation device
US10710610B2 (en) 2016-05-27 2020-07-14 Nissan Motor Co., Ltd. Driving control method and driving control apparatus
CN109153388A (en) * 2016-05-27 2019-01-04 日产自动车株式会社 Driving control method and steering control device
JPWO2017203694A1 (en) * 2016-05-27 2019-03-22 日産自動車株式会社 Operation control method and operation control apparatus
CN109153388B (en) * 2016-05-27 2019-09-03 日产自动车株式会社 Driving control method and driving control device
WO2017203694A1 (en) * 2016-05-27 2017-11-30 日産自動車株式会社 Driving control method and driving control device
JP2019159147A (en) * 2018-03-14 2019-09-19 クラリオン株式会社 On-vehicle device, traveling state estimation method, server device, information processing method, and traveling state estimation system
WO2019176356A1 (en) * 2018-03-14 2019-09-19 クラリオン株式会社 Onboard device, traveling state estimation method, server device, information processing method, and traveling state estimation system
JP7100471B2 (en) 2018-03-14 2022-07-13 フォルシアクラリオン・エレクトロニクス株式会社 In-vehicle device, driving state estimation method, server device, information processing method, and driving state estimation system
US11498576B2 (en) 2018-03-14 2022-11-15 Clarion Co., Ltd. Onboard device, traveling state estimation method, server device, information processing method, and traveling state estimation system
JP2019156355A (en) * 2018-03-16 2019-09-19 本田技研工業株式会社 Vehicle control device
JP2021051688A (en) * 2019-09-26 2021-04-01 トヨタ自動車株式会社 Drive skill evaluation device
JP7302411B2 (en) 2019-09-26 2023-07-04 トヨタ自動車株式会社 Driving skill evaluation device

Also Published As

Publication number Publication date
JP5267203B2 (en) 2013-08-21

Similar Documents

Publication Publication Date Title
JP5267203B2 (en) Driving state determination device and driving state determination method
US9256989B2 (en) Information provision device for use in vehicle
JP5621921B2 (en) Operation instability judgment device
JP4796400B2 (en) Vehicle speed control device, target speed setting method and program in the same
JP6161942B2 (en) Curve shape modeling device, vehicle information processing system, curve shape modeling method, and curve shape modeling program
JP5402335B2 (en) Reliability evaluation apparatus, reliability evaluation method, and reliability evaluation program
JP5130181B2 (en) Wheel vibration extraction device and road surface state estimation device
US9573597B2 (en) Vehicle information processing apparatus and vehicle information processing method
WO2015145819A1 (en) Map information generation system, method, and program
JP6815406B2 (en) Speed calculator, control method, program and storage medium
KR20140093961A (en) Sensor system for independently evaluating the accuracy of the data of the sensor system
JP2016156973A (en) Map data storage device, control method, program and recording medium
US20140015663A1 (en) Vehicular information-processing device and vehicular information-processing method
WO2015151594A1 (en) Driving assistance system, method, and program
KR20190038739A (en) Method for detecting the changing point of road
JP5904226B2 (en) Vehicle behavior prediction apparatus and program
US11938952B2 (en) Traveling trajectory estimation system, traveling trajectory estimation program, and traveling trajectory estimating method
JP2012059058A (en) Risk estimation device and program
WO2019188886A1 (en) Terminal device, information processing method, and storage medium
WO2015029565A1 (en) Driving assistance system, method and program
JP6448347B2 (en) Deceleration control system, method and program
JP6008049B2 (en) Vehicle information providing device
JP2008077349A (en) Vehicle state quantity estimation device and vehicle steering control device using the device
JP7456749B2 (en) Wear amount estimation system and calculation model generation system
KR20210087681A (en) Apparatus for measuring position of vehicle and methdo thereof

Legal Events

Date Code Title Description
RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20100917

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20111219

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20130311

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20130409

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20130422

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

Ref document number: 5267203

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150