[go: up one dir, main page]

JP5665400B2 - Vehicle control device - Google Patents

Vehicle control device Download PDF

Info

Publication number
JP5665400B2
JP5665400B2 JP2010164216A JP2010164216A JP5665400B2 JP 5665400 B2 JP5665400 B2 JP 5665400B2 JP 2010164216 A JP2010164216 A JP 2010164216A JP 2010164216 A JP2010164216 A JP 2010164216A JP 5665400 B2 JP5665400 B2 JP 5665400B2
Authority
JP
Japan
Prior art keywords
vehicle
estimated
traveling place
control algorithm
place
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.)
Expired - Fee Related
Application number
JP2010164216A
Other languages
Japanese (ja)
Other versions
JP2012025232A (en
Inventor
由幸 松岡
由幸 松岡
和夫 中澤
和夫 中澤
了 古郡
了 古郡
隆雅 末冨
隆雅 末冨
黒田 康秀
康秀 黒田
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.)
Mazda Motor Corp
Keio University
Original Assignee
Mazda Motor Corp
Keio University
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 Mazda Motor Corp, Keio University filed Critical Mazda Motor Corp
Priority to JP2010164216A priority Critical patent/JP5665400B2/en
Publication of JP2012025232A publication Critical patent/JP2012025232A/en
Application granted granted Critical
Publication of JP5665400B2 publication Critical patent/JP5665400B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/60Other road transportation technologies with climate change mitigation effect
    • Y02T10/62Hybrid vehicles

Landscapes

  • Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
  • Electric Propulsion And Braking For Vehicles (AREA)
  • Hybrid Electric Vehicles (AREA)

Description

本発明は車両用制御装置に係り、より詳細には、車両状態、車外状態及び乗員状態に応じて車両の制御アルゴリズムを動的に変化させる車両用制御装置に関する。   The present invention relates to a vehicle control device, and more particularly to a vehicle control device that dynamically changes a vehicle control algorithm in accordance with a vehicle state, a vehicle exterior state, and an occupant state.

下記の特許文献1には、条件により目的関数が変化する最適化対象の決定変数を短時間で効率良く最適化する方法が記載されている。また、下記の特許文献2には、ドライバの運転操作と外部環境との関連性を随時学習して、ドライバの普段の内部状態からの逸脱を認識する方法が記載されている。   Patent Document 1 below describes a method for efficiently optimizing a decision variable to be optimized whose objective function changes depending on conditions in a short time. Patent Document 2 below describes a method of recognizing a deviation from the driver's normal internal state by learning the relationship between the driving operation of the driver and the external environment as needed.

特開2009−99051号公報JP 2009-99051 A 特開2009−73465号公報JP 2009-73465 A

ところで、上記特許文献1に記載の技術では、条件により目的変数が変化する最適化対称の決定変数を短時間で効率よく最適化する方法が提案されているが、目的関数を運転行動から推定し自律的に変化させることは困難である。   By the way, in the technique described in Patent Document 1, a method for efficiently optimizing an optimization symmetric decision variable whose objective variable changes depending on conditions is proposed in a short time, but the objective function is estimated from driving behavior. It is difficult to change autonomously.

また、上記特許文献2に記載の技術では、ドライバによる運転操作と外部環境との関連性を随時学習し、ドライバの普段の内部状態からの逸脱を認識する方法が提案されているが、プログラムやモデル構造を動的に変更することは困難である。   Further, in the technique described in Patent Document 2, a method has been proposed in which the relationship between the driving operation by the driver and the external environment is learned at any time, and the deviation from the usual internal state of the driver is recognized. It is difficult to change the model structure dynamically.

このように、従来の車両用制御装置では、制御アルゴリズムの設計段階で想定されていなかった状況に車両が遭遇したときに、その状況に対応して車両制御を行うことは困難であった。   As described above, in the conventional vehicle control device, when the vehicle encounters a situation that was not assumed at the design stage of the control algorithm, it is difficult to perform vehicle control corresponding to the situation.

そこで、本発明は、車両が、想定されていなかった走行場に遭遇したときに、その走行場に対応した新たな制御アルゴリズムを自律的に創発することができる車両用制御装置を提供することを目的としている。   Therefore, the present invention provides a vehicle control device capable of autonomously generating a new control algorithm corresponding to a travel location when the vehicle encounters an unexpected travel location. It is aimed.

なお、車両の「走行場」とは、車両の遭遇している状況、状態、事態又は場面を含み、停止状態を含む車両の走行状態、車外環境の状況、及び車両のドライバの挙動、認識(例えば危機感)及び意図をも含む。また、制御アルゴリズムの「創発」とは、新しい制御アルゴリズムを生成することを含む。創発に当たっては、部分の性質の単純な総和にとどまらない性質が、全体として現れることがある。例えば、局所的な複数の相互作用が複雑に組織化することで、個別の要素の振る舞いからは予測できないようなシステムが構成されることもある。   The “running place” of the vehicle includes the situation, state, situation or scene that the vehicle is encountering, including the running state of the vehicle including the stop state, the situation of the outside environment, and the behavior and recognition of the driver of the vehicle ( For example, a sense of crisis) and intention. The “emergence” of the control algorithm includes generating a new control algorithm. In emergence, properties that are not just a simple sum of the properties of parts may appear as a whole. For example, a system that cannot be predicted from the behavior of individual elements may be configured by a complex organization of a plurality of local interactions.

上記の目的を達成するため、本発明の車両用制御装置は、車両を運転するドライバの状態、車外環境状態、及び車両状態の少なくとも一つを検出する検出手段と、上記検出手段によって検出された検出状態に基づいて、車両が遭遇している走行場を、恒等写像モデルを利用して推定する推定手段と、既知の走行場と車両の制御アルゴリズムとを互いに対応付けて記憶した第1記憶手段と、上記推定手段によって推定された推定走行場が、上記第1記憶手段に記憶された既知の走行場であるか、上記第1記憶手段に記憶されていない未知の走行場であるかを判定する判定手段と、上記判定手段によって推定走行場が未知の走行場であると判定された場合に、遺伝的ネットワークプログラミングを利用して、推定走行場に対応する新たな制御アルゴリズムを創発する創発手段と、上記判定手段によって推定走行場が既知の走行場であると判定された場合に、上記第1記憶手段から獲得した推定走行場に対応する制御アルゴリズムを車両の制御パラメータに変換し、推定走行場が未知の走行場であると判定された場合に、上記創発手段によって創発された新たな制御アルゴリズムを車両の制御パラメータに変換する変換手段と、を有することを特徴としている。 In order to achieve the above object, a vehicle control apparatus according to the present invention is detected by a detection means for detecting at least one of a state of a driver driving a vehicle, an external vehicle environment state, and a vehicle state, and the detection means. Based on the detection state, a first memory that stores an estimation unit that estimates a traveling place that the vehicle encounters using an identity mapping model, and a known traveling place and a vehicle control algorithm that are associated with each other. And the estimated traveling place estimated by the estimating means is a known traveling place stored in the first storage means or an unknown traveling place that is not stored in the first storage means. A new control algorithm corresponding to the estimated travel location using genetic network programming when the determination means determines that the estimated travel location is an unknown travel location. A control algorithm corresponding to the estimated travel location obtained from the first storage means when the estimated travel location is determined to be a known travel location by the emergence means for generating a rhythm and the determination means is a control parameter of the vehicle And converting means for converting a new control algorithm created by the generating means into vehicle control parameters when the estimated driving place is determined to be an unknown driving place. Yes.

このように構成された本発明によれば、車両が、想定されていなかった走行場に遭遇したときに、その走行場に対応した新たな制御アルゴリズムを自律的に創発することができる車両用制御装置を提供することができる。
また、遺伝的ネットワークプログラミングでは、有向グラフ(ネットワーク)を用いて遺伝子型が表現されるため、従来の遺伝的アルゴリズムでは表現が困難であった数式やプログラムのコードなど、構造を有するデータを表現することができる。このため、遺伝的ネットワークプログラミングは、未知の走行場に対する新たな制御アルゴリズムを創発するのに用いて好適である。
また、恒等写像モデルは、入力層と出力層に同一データを提示して、恒等写像としての中間層出力を学習させる階層型ニューラルネットワークである。恒等写像モデルでは、曖昧で不完全な多種多様な情報を入力及び出力とし、問題解決に必要な縮約された情報を自動的に抽出することができる。特に、恒等写像モデルでは、教師信号を必要としないため、未知の予め想定されていなかった未知の走行場を推定するのに用いて好適である。
According to the present invention configured as described above, when the vehicle encounters an unexpected driving field, a vehicle control capable of autonomously generating a new control algorithm corresponding to the driving field. An apparatus can be provided.
In genetic network programming, a genotype is expressed using a directed graph (network), so that data with a structure such as mathematical formulas and program codes that were difficult to express with conventional genetic algorithms is expressed. Can do. For this reason, genetic network programming is suitable for generating new control algorithms for unknown driving fields.
The identity mapping model is a hierarchical neural network that presents the same data to the input layer and the output layer and learns the intermediate layer output as the identity map. In the identity mapping model, a variety of vague and incomplete information can be input and output, and the reduced information necessary for problem solving can be automatically extracted. In particular, since the identity mapping model does not require a teacher signal, it is suitable for use in estimating an unknown traveling place that has not been assumed in advance.

また、本発明において、好ましくは、上記第1記憶手段に格納された、既知の走行場に対する制御アルゴリズムは、遺伝的ネットワークプログラミングを利用して創発したものである。
遺伝的ネットワークプログラミングでは、有向グラフ(ネットワーク)を用いて遺伝子型が表現されるため、従来の遺伝的アルゴリズムでは表現が困難であった数式やプログラムのコードなど、構造を有するデータを表現することができる。このため、遺伝的ネットワークプログラミングは、既知の走行場に対する適切な制御アルゴリズムを創発するのに用いて好適である。
In the present invention, it is preferable that the control algorithm for the known traveling place stored in the first storage means is created using genetic network programming.
In genetic network programming, genotypes are expressed using directed graphs (networks), so it is possible to express structured data such as mathematical formulas and program codes that were difficult to express with conventional genetic algorithms. . For this reason, genetic network programming is suitable for use in generating appropriate control algorithms for known driving fields.

また、本発明において、好ましくは、更に、前記創発手段によって創発された新たな制御アルゴリズムを記憶する第3記憶手段を有する。
想定外の走行場に対応した新たな制御アルゴリズムを蓄積することによって、より多様な走行場に対応して車両制御を行うことが可能となる。
In the present invention, it is preferable to further include a third storage unit that stores a new control algorithm generated by the generating unit.
By accumulating new control algorithms corresponding to unexpected driving fields, it becomes possible to perform vehicle control corresponding to more diverse driving fields.

本発明の車両用制御装置によれば、車両が、想定されていなかった走行場に遭遇したときに、その走行場に対応した新たな制御アルゴリズムを自律的に創発することができる。   According to the vehicle control device of the present invention, when a vehicle encounters a driving field that has not been assumed, a new control algorithm corresponding to the driving field can be autonomously generated.

本発明の実施形態による車両用制御装置の構成を示すブロック図である。It is a block diagram which shows the structure of the control apparatus for vehicles by embodiment of this invention. 本発明の実施形態による車両用制御装置の作動を示すフローチャートである。It is a flowchart which shows the action | operation of the control apparatus for vehicles by embodiment of this invention. 自車挙動の走行場を推定するアルゴリズムのフローチャートである。It is a flowchart of the algorithm which estimates the traveling place of the own vehicle behavior. 自車挙動の走行場を推定した結果を示すグラフであり、(A)は車速の時間変化を示し、(B)は加速度の時間変化を示す。It is a graph which shows the result of having estimated the run field of the own vehicle behavior, (A) shows time change of vehicle speed, and (B) shows time change of acceleration. 相対距離及び相対速度と推定された車両関係の走行場との関係を示すグラフである。It is a graph which shows the relationship between relative distance and relative speed, and the estimated vehicle-related driving field. ドライバの危機感の走行場を推定するアルゴリズムのフローチャートである。It is a flowchart of the algorithm which estimates the driving field of a driver's sense of crisis. ドライバの危機感の走行場を推定した結果を示すグラフであり、(A)はブレーキ踏込角の時間変化を示し、(B)はブレーキ微分値の時間変化を示す。It is a graph which shows the result of having estimated the driving field of a driver's crisis feeling, (A) shows time change of a brake depression angle, and (B) shows time change of a brake differential value. 遺伝的アルゴリズム(GA)の基本的なアルゴリズムの流れを示すフローチャートである。It is a flowchart which shows the flow of the basic algorithm of a genetic algorithm (GA). 遺伝的ネットワークプログラミング(GNP)における集団と、その集団を構成する一つの個体の有向グラフの模式図である。FIG. 3 is a schematic diagram of a group in genetic network programming (GNP) and a directed graph of one individual constituting the group. 遺伝的ネットワークプログラミング(GNP)の一個体の遺伝子の表現型、遺伝子型及びノードライブラリの一例を示す図表である。It is a graph which shows an example of the phenotype of a gene of an individual of genetic network programming (GNP), a genotype, and a node library. 遺伝的ネットワークプログラミング(GNP)の交叉の例を示す模式図である。It is a schematic diagram which shows the example of the crossover of a genetic network programming (GNP). 遺伝的ネットワークプログラミング(GNP)の突然変異の例を示す模式図である。It is a schematic diagram which shows the example of the mutation of a genetic network programming (GNP). シリーズ方式ハイブリッド電気自動車の車両システムの模式図である。It is a schematic diagram of the vehicle system of a series system hybrid electric vehicle. 車両諸元のパラメータの図表である。It is a chart of a parameter of vehicle specifications. エンジン出力と燃料消費率との関係を示すグラフである。It is a graph which shows the relationship between an engine output and a fuel consumption rate. バッテリ構成の模式図である。It is a schematic diagram of a battery structure. 遺伝的ネットワークプログラミング(GNP)の流れを示すフローチャートである。It is a flowchart which shows the flow of genetic network programming (GNP). 対応テーブルを用いた制御アルゴリズムの創発処理を示すフローチャートである。It is a flowchart which shows the emergence process of the control algorithm using a correspondence table. 走行場の模式図である。It is a schematic diagram of a traveling place. 推定された走行場と、対応テーブルとの関係を示す模式図である。It is a schematic diagram which shows the relationship between the estimated driving field and a corresponding | compatible table. 想定外の走行場に対する制御アルゴリズムの創発処理を示すフローチャートである。It is a flowchart which shows the emergence process of the control algorithm with respect to an unexpected driving field. 遺伝的ネットワークプログラミング(GNP)の有向グラフ構造の模式図である。It is a schematic diagram of the directed graph structure of genetic network programming (GNP). 第2実施形態における走行場の推定に用いる恒等写像モデルの模式図である。It is a schematic diagram of the identity mapping model used for estimation of a driving field in 2nd Embodiment. 第2実施形態における恒等写像モデルの中間層出力の模式図である。It is a schematic diagram of the intermediate | middle layer output of the identity mapping model in 2nd Embodiment. (A)は、第2実施形態における想定範囲の模式図であり、(B)は、想定範囲と対応テーブルとの関係を示す模式図である。(A) is a schematic diagram of the assumed range in 2nd Embodiment, (B) is a schematic diagram which shows the relationship between an assumed range and a corresponding | compatible table.

(1.第1実施形態)
以下、添付の図面を参照して、本発明の車両用制御装置の実施形態を説明する。
まず、図1を参照して、本発明の第1実施形態による車両用制御装置の構成について説明する。図1は、実施形態による車両用制御装置の構成を示すブロック図である。
(1. First embodiment)
Hereinafter, an embodiment of a vehicle control device of the present invention will be described with reference to the accompanying drawings.
First, the configuration of the vehicle control device according to the first embodiment of the present invention will be described with reference to FIG. FIG. 1 is a block diagram illustrating a configuration of a vehicle control device according to an embodiment.

図1に示すように、実施形態による車両用制御装置は、車両を運転するドライバの状態を検出するドライバ状態モニタ11、車外環境状態を検出する車外環境モニタ12及び車両状態を検出する車両状態モニタ13を含む検出手段10と、この検出手段10によって検出された検出状態に基づいて、車両が遭遇している走行場を推定する推定部20と、過去の経験に基づく既知の走行場と車両の制御アルゴリズムとを互いに対応付けて記憶した第1記憶部30と、推定部11によって推定された推定走行場が、第1記憶部30に記憶された既知の走行場から想定される想定内の走行場であるか、想定外の走行場であるかを判定する判定部40と、この判定部40によって推定走行場が想定外の走行場であると判定された場合に、推定走行場に対応する制御アルゴリズムを創発する創発部50と、判定部手段によって推定走行場が想定内の走行場であると判定された場合に、上記第1記憶手段から獲得した推定走行場に対応する制御アルゴリズムを車両の制御パラメータに変換し、推定走行場が想定外の走行場であると判定された場合に、上記創発手段によって創発された制御アルゴリズムを車両の制御パラメータに変換する変換部60とを有している。   As shown in FIG. 1, the vehicle control apparatus according to the embodiment includes a driver state monitor 11 that detects a state of a driver who drives the vehicle, an external environment monitor 12 that detects an external environment state, and a vehicle state monitor that detects the vehicle state. 13, a detection means 10 including 13, an estimation unit 20 that estimates a traveling place where the vehicle encounters based on a detection state detected by the detection means 10, a known traveling place based on past experience, and the vehicle The first storage unit 30 that stores the control algorithm in association with each other, and the estimated traveling place estimated by the estimating unit 11 is assumed to be assumed from the known traveling place stored in the first storage unit 30. When the determination unit 40 determines whether the estimated traveling place is an unexpected traveling place or not, the determination unit 40 determines whether the traveling place is an unexpected traveling place. A control algorithm corresponding to the estimated driving field acquired from the first storage means when the generating unit 50 for generating a corresponding control algorithm and the determining unit determine that the estimated driving field is an expected driving field Is converted into a vehicle control parameter, and when it is determined that the estimated driving field is an unexpected driving field, a conversion unit 60 is provided for converting the control algorithm generated by the generating means into the vehicle control parameter. doing.

さらに、本実施形態では、過去の経験に基づかない走行場と車両の制御アルゴリズムとを互いに対応付けた対応テーブルを記憶した第2記憶部70と、創発部50によって創発された新たな制御アルゴリズムを記憶する第3記憶部80とを有している。   Furthermore, in the present embodiment, the second storage unit 70 that stores a correspondence table in which a running field that is not based on past experience and a vehicle control algorithm are associated with each other, and a new control algorithm that has been created by the emergence unit 50 are provided. And a third storage unit 80 for storing.

なお、推定部20,判定部40、創発部50及び変換部60は、例えば車載ECU(electric control unit:電子制御装置)における処理機能に相当する。これらの処理機能は、コンピュータにおいて所定のプログラムを実行することによって実現してもよいし、マイクロチップによって実現してもよい。また、第1記憶部30、第2記憶部70及び第3記憶部80は、半導体メモリチップやハードディスクといった記憶装置で構成するとよい。   The estimation unit 20, the determination unit 40, the emergence unit 50, and the conversion unit 60 correspond to, for example, a processing function in an in-vehicle ECU (electric control unit). These processing functions may be realized by executing a predetermined program in a computer, or may be realized by a microchip. The first storage unit 30, the second storage unit 70, and the third storage unit 80 may be configured by a storage device such as a semiconductor memory chip or a hard disk.

ここで、図2のフローチャートを参照して、第1実施形態による車両用制御装置の動作について説明する。
(1−1.状態検出)
まず、検出手段10が、車両を運転するドライバの状態、車外環境状態、及び車両状態車両状況を検出する(S21)。
Here, with reference to the flowchart of FIG. 2, operation | movement of the control apparatus for vehicles by 1st Embodiment is demonstrated.
(1-1. State detection)
First, the detection means 10 detects the state of the driver who drives the vehicle, the environment state outside the vehicle, and the vehicle state vehicle situation (S21).

検出部10を構成するドライバ状態モニタ11の例としては、車室内に配置されてドライバの視線や、顔の動きを撮像するモニタカメラ及び画像処理システムや、ドライバのブレーキの踏込量を検出するセンサが挙げられる。また、車外環境モニタ12の例としては、先行車両を含む障害物までの距離や相対速度を検出する車両前方に向けたレーダや、車両の側方又は後方に向けたレーダや、車外に向けたモニタカメラ及び画像処理システムが挙げられる。また、車両状態モニタ13の例としては、車速センサ、加速度センサ、ブレーキ踏込量センサ、アクセル踏込量センサ、及び操舵角センサが挙げられる。   Examples of the driver status monitor 11 that constitutes the detection unit 10 include a monitor camera and an image processing system that are arranged in a vehicle interior and that captures the driver's line of sight and facial movements, and a sensor that detects the amount of depression of the driver's brake Is mentioned. Moreover, as an example of the vehicle exterior environment monitor 12, the radar directed to the front of the vehicle for detecting the distance to the obstacle including the preceding vehicle and the relative speed, the radar directed to the side or the rear of the vehicle, or directed to the outside of the vehicle Examples include monitor cameras and image processing systems. Examples of the vehicle state monitor 13 include a vehicle speed sensor, an acceleration sensor, a brake depression amount sensor, an accelerator depression amount sensor, and a steering angle sensor.

なお、検出部10は、ドライバ状態モニタ11、車外環境モニタ12及び車両状態モニタ13の三つを常に設ける必要はなく、これらのうちの一つだけ又は二つだけを設けたものであってもよい。   The detection unit 10 does not always need to be provided with the driver status monitor 11, the vehicle exterior environment monitor 12, and the vehicle status monitor 13, and may be provided with only one or two of them. Good.

(1−2.走行場の推定)
次に、推定部20が、検出手段10によって検出された検出状態に基づいて、車両が遭遇している走行場を推定する(S22)。
(1-2. Estimation of driving field)
Next, the estimation unit 20 estimates the traveling place where the vehicle encounters based on the detection state detected by the detection means 10 (S22).

以下、車両状況に基づいて自車の挙動を推定し、車外環境状況に基づいて、先行車両と自車両との車両関係を推定し、さらに、ドライバ状態に基づいてドライバの危機感を推定する例について説明する。   Hereinafter, an example in which the behavior of the host vehicle is estimated based on the vehicle status, the vehicle relationship between the preceding vehicle and the host vehicle is estimated based on the environment outside the vehicle, and the driver's sense of crisis is estimated based on the driver state. Will be described.

(1−2−1.自車挙動の走行場の推定)
一般に、車両は、停止した状態から発進し、走行する。そして、加速及び減速を繰り返して最終的に停止する。このように、車両の挙動は、「停止」、「発進」及び「走行」というサイクルを繰り返す。このため、自車挙動についての走行場として、下記の表1に示すように、「発進」、「走行」及び「停止」の3つの走行場を推定する。
(1-2-1. Estimating the driving field of own vehicle behavior)
Generally, a vehicle starts from a stopped state and travels. Then, acceleration and deceleration are repeated to finally stop. Thus, the behavior of the vehicle repeats the cycle of “stop”, “start”, and “travel”. For this reason, as shown in Table 1 below, three driving locations of “start”, “travel”, and “stop” are estimated as the travel locations for the behavior of the host vehicle.

Figure 0005665400
Figure 0005665400

自車挙動の走行場を推定するにあたり、車速V、及びその時間微分である加速度αに着目すると、自車両は、車速V=0で「停止」している状態から動き始め、加速度α>0となって「発進」へ移行する。さらに、「発進」後に車速が極大値をとり、加速度α<0となって「走行」へ移行する。そして、「走行」から減速し、車速V=0となって「停止」へ移行する。このため、車速V及び加速度αは、自車挙動の走行場を推定するためのパラメータとして使用することができる。   In estimating the traveling place of the own vehicle behavior, paying attention to the vehicle speed V and the acceleration α which is a time derivative thereof, the own vehicle starts to move from the state of “stop” at the vehicle speed V = 0, and the acceleration α> 0. And move to "start". Further, after “starting”, the vehicle speed reaches the maximum value, and the acceleration α <0, and the routine proceeds to “running”. Then, the vehicle decelerates from “travel”, and the vehicle speed V = 0 and shifts to “stop”. For this reason, the vehicle speed V and the acceleration α can be used as parameters for estimating the traveling place of the vehicle behavior.

図3に、車速V及び加速度αから自車挙動の走行場を推定するアルゴリズムのフローチャートを示す。入力される車速V>0となるまで(S32で「no」の場合)、走行場は「停止」と推定される(S31)。次いで、車速V>0となると(S32で「yes」の場合)、入力される加速度α<0となるまで(S34で「no」の場合)、走行場は「発進」と推定される(S33)。次いで、加速度<0となると(S34で「yes」の場合)、入力される車速V=0となるまで(S36で「no」の場合)、走行場は「走行」と推定される(S35)。そして、車速V=0となると(S36で「yes」の場合)、フローチャートはスタートへ戻り、走行場は「停止」と推定される。   FIG. 3 shows a flowchart of an algorithm for estimating the traveling place of the own vehicle behavior from the vehicle speed V and the acceleration α. Until the input vehicle speed V> 0 (in the case of “no” in S32), the traveling place is estimated to be “stop” (S31). Next, when the vehicle speed V> 0 (in the case of “yes” in S32), the driving field is estimated to be “start” until the input acceleration α <0 (in the case of “no” in S34) (S33). ). Next, when acceleration <0 (in the case of “yes” in S34), the traveling place is estimated to be “traveling” until the input vehicle speed V = 0 (in the case of “no” in S36) (S35). . Then, when the vehicle speed V = 0 (in the case of “yes” in S36), the flowchart returns to the start, and the traveling place is estimated to be “stop”.

図4に、車速V及び加速度αで表された車両の走行パターンに対して、図3に示した推定アルゴリズムを適用して自車挙動の走行場を推定した結果を示す。図4の(A)のグラフ中の線Iは、車速Vの時間変化を示し、図4(B)のグラフ中の線IIは、加速度αの時間変化を示す。   FIG. 4 shows a result of estimating the traveling place of the own vehicle behavior by applying the estimation algorithm shown in FIG. 3 to the traveling pattern of the vehicle represented by the vehicle speed V and the acceleration α. A line I in the graph of FIG. 4A indicates a time change of the vehicle speed V, and a line II in the graph of FIG. 4B indicates a time change of the acceleration α.

図4に示す例では、時刻t0までは車速V=0であるため、走行場は「停止」と推定される。そして、時刻t0に車速V>0となってから時刻t1に加速度α<0となるまで、走行場は「発進」と推定される。さらに、時刻t1に加速度α<0となってから時刻t2に車速V=0となるまで、走行場は「走行」と推定される。そして、時刻t2以降、車速V=0となり、走行場は「停止」と推定される。このようにして、車速V及び加速度αに基づいて、自車挙動の走行場を推定することができる。   In the example shown in FIG. 4, since the vehicle speed V = 0 until time t0, the traveling place is estimated to be “stop”. Then, from the vehicle speed V> 0 at time t0 to the acceleration α <0 at time t1, the running field is estimated to be “start”. Further, the traveling place is estimated to be “traveling” from the time when acceleration α <0 at time t1 until the vehicle speed V = 0 at time t2. Then, after time t2, the vehicle speed V = 0, and the traveling place is estimated to be “stopped”. In this way, it is possible to estimate the traveling place of the vehicle behavior based on the vehicle speed V and the acceleration α.

(1−2−2.車両間系の走行場の推定)
自車両と先行車両との車両関係についての走行場として、下記の表2に示すように、「追従」、「自由」、「縮まる」、「拡がる」、「接近」及び「衝突」の6つの走行場を推定する。
(1-2-2. Estimation of the driving field of the inter-vehicle system)
As shown in Table 2 below, there are six areas of “following”, “free”, “shrinking”, “expanding”, “approaching”, and “collision”. Estimate the driving field.

Figure 0005665400
Figure 0005665400

ここでは、自車両と先行車両との車間距離L、及び、自車両に対する先行車両の相対速度Vrを、車両関係の走行場を推定するためのパラメータとして使用する。具体的には、図5に示すように、車間距離Lが所定範囲内(L1≦L≦L2)であり、かつ、相対速度VrがVr=0を挟んだ所定範囲内(−V1≦Vr≦+V1)であるとき、走行場を「追従」と推定する。また、車間距離Lが「追従」の上限車間距離より長いとき(L2<L)、走行場を「自由」と推定する。また、車間距離Lは「追従」の所定範囲内(L1≦L≦L2)であるが、相対速度Vrが「追従」の下限値未満であるとき(Vr<−V1)、走行場を「縮まる」と推定する。また、車間距離Lは「追従」の所定範囲内(L1≦L≦L2)であるが、相対速度Vrが「追従」の上限値を超えるとき(+V1<Vr)、走行場を「拡がる」と推定する。また、車間距離Lが、「追従」の下限値未満であり、かつ0以上であるとき(0≦L<L1)、走行場を「接近」と推定する。そして、車間距離Lが0未満であるとき(L<0)、「衝突」と推定する。   Here, the inter-vehicle distance L between the host vehicle and the preceding vehicle and the relative speed Vr of the preceding vehicle with respect to the host vehicle are used as parameters for estimating a vehicle-related traveling place. Specifically, as shown in FIG. 5, the inter-vehicle distance L is within a predetermined range (L1 ≦ L ≦ L2), and the relative speed Vr is within a predetermined range with Vr = 0 therebetween (−V1 ≦ Vr ≦ + V1), it is estimated that the driving field is “following”. Further, when the inter-vehicle distance L is longer than the upper limit inter-vehicle distance of “follow-up” (L2 <L), the traveling place is estimated as “free”. The inter-vehicle distance L is within a predetermined range of “following” (L1 ≦ L ≦ L2), but when the relative speed Vr is less than the lower limit value of “following” (Vr <−V1), the travel field is “shrinked”. ". The inter-vehicle distance L is within a predetermined range of “following” (L1 ≦ L ≦ L2), but when the relative speed Vr exceeds the upper limit value of “following” (+ V1 <Vr), the driving field is “expanded”. presume. Further, when the inter-vehicle distance L is less than the lower limit value of “following” and 0 or more (0 ≦ L <L1), the traveling place is estimated as “approaching”. When the inter-vehicle distance L is less than 0 (L <0), it is estimated as “collision”.

図5に、相対距離L及び相対速度Vrと、推定された車両関係の走行場との関係を示す。図5のグラフの横軸は先行車両と自車両との車間距離Lを表し、縦軸は自車両の速度に対する先行車両の相対速度Vrを表す。図5に示すように、相対距離L及び相対速度Vrによって、車両関係の走行場が推定される。   FIG. 5 shows the relationship between the relative distance L and the relative speed Vr, and the estimated vehicle-related traveling place. The horizontal axis of the graph in FIG. 5 represents the inter-vehicle distance L between the preceding vehicle and the host vehicle, and the vertical axis represents the relative speed Vr of the preceding vehicle with respect to the speed of the host vehicle. As shown in FIG. 5, a vehicle-related traveling place is estimated from the relative distance L and the relative speed Vr.

(1−2−3.ドライバ危機感の走行場の推定)
ドライバの危機感についての走行場として、下記の表3に示すように、「通常」、「前兆」及び「ヒヤリハット」の3つの走行場を推定する。一般に、ドライバは、飛び出しや急な割り込みといった緊急事態に遭遇した場合に、急ブレーキをかける。そこで、ブレーキ踏込角(ブレーキ踏込量)B、及びブレーキ踏込角の時間微分値(以下、「ブレーキ微分値」ともいう。)βをドライバの危機感についての走行場を推定するためのパラメータとして使用する。
(1-2-3. Estimating the driving field with a sense of driver crisis)
As shown in Table 3 below, three driving locations of “normal”, “precursor”, and “near hat” are estimated as driving locations for the driver's sense of crisis. Generally, a driver applies a brake suddenly when encountering an emergency such as jumping out or sudden interruption. Therefore, the brake depression angle (brake depression amount) B and the brake differential angle time differential value (hereinafter also referred to as “brake differential value”) β are used as parameters for estimating the driving field for the driver's sense of crisis. To do.

Figure 0005665400
Figure 0005665400

図6に、ブレーキ踏込角B及びブレーキ微分値βからドライバの危機感の走行場を推定するアルゴリズムのフローチャートを示す。入力されるブレーキ微分βが所定の閾値βTを超えてβ>βTとなるまで(S62で「no」の場合)、走行場は「通常」と推定される(S61)。次いで、ブレーキ微分値β>βTとなると(S62で「yes」の場合)、入力されるブレーキ踏込角Bが所定の閾値BTを超えてB>BTとなるまで(S654で「no」の場合)、走行場は「前兆」と推定される(S63)。   FIG. 6 shows a flowchart of an algorithm for estimating a driver's critical driving field from the brake depression angle B and the brake differential value β. Until the input brake differential β exceeds a predetermined threshold βT and β> βT is satisfied (in the case of “no” in S62), the traveling field is estimated to be “normal” (S61). Next, when the brake differential value β> βT is satisfied (in the case of “yes” in S62), the input brake depression angle B exceeds the predetermined threshold value BT and becomes B> BT (in the case of “no” in S654). The traveling place is estimated to be “a precursor” (S63).

なお、ブレーキ踏込角B=0となると(S64で「yes」の場合)、フローチャートはスタートへ戻り、走行場は「停止」と推定される。すなわち、全ての「前兆」が「ヒヤリハット」へ移行するわけではない。これは、ブレーキ踏込量Bが小さくても、ブレーキ踏込量の時間微分値βだけが大きな値となることがあるためである。   Note that when the brake depression angle B becomes 0 (in the case of “yes” in S64), the flowchart returns to the start, and the traveling place is estimated to be “stop”. That is, not all “precursors” will shift to “near-miss”. This is because even if the brake depression amount B is small, only the time differential value β of the brake depression amount may be a large value.

次いで、ブレーキ踏込角B>BTとなると(S65で「yes」の場合)、入力されるブレーキ踏込角B<BTとなるまで(S67で「yes」の場合)、走行場は「ヒヤリハット」即ちドライバが危機感を感じている状態と推定される。次いで、入力されるブレーキ踏込角B<BTとなると(S67で「no」の場合)、走行場は再び「前兆」と推定される(S63)。   Next, when the brake depression angle B> BT is satisfied (in the case of “yes” in S65), the traveling place is “near-miss”, that is, the driver until the input brake depression angle B <BT (in the case of “yes” in S67). Is estimated to be in a state of crisis. Next, when the input brake depression angle B <BT is satisfied (in the case of “no” in S67), the traveling place is estimated again as “precursor” (S63).

図7に、ブレーキ踏込角B及びブレーキ微分値βで表された走行パターンに対して、図6に示したドライバの危機感の推定アルゴリズムを適用して走行場を推定した結果を示す。図7(A)のグラフ中の線Iは、ブレーキ踏込角Bの時間変化を示し、図7(B)のグラフ中の線IIは、ブレーキ微分値βの時間変化を示す。   FIG. 7 shows the result of estimating the driving field by applying the driver's criticality estimation algorithm shown in FIG. 6 to the driving pattern represented by the brake depression angle B and the brake differential value β. A line I in the graph of FIG. 7A shows a time change of the brake depression angle B, and a line II in the graph of FIG. 7B shows a time change of the brake differential value β.

図7に示す例では、時刻t0まではブレーキ微分値β<βTであるため、走行場は「通常」と推定される。そして、時刻t0にブレーキ微分値β>βTとなってから時刻t1にブレーキ踏込角B>BTとなるまで、走行場は「前兆」と推定される。さらに、時刻t1にブレーキ踏込角B>BTとなってから時刻t2にブレーキ踏込角B<BTとなるまで、走行場は「ヒヤリハット」と推定される。さらに、時刻t2にブレーキ踏込角B<BTとなってから時刻t3にブレーキ踏込角B=0となるまで、走行場は「前兆」と推定される。そして、時刻t3にブレーキ踏込角B=0となってから時刻t4に再びブレーキ微分値β>βTとなるまで、走行場は「通常」と推定される。   In the example shown in FIG. 7, since the brake differential value β <βT until time t0, the traveling place is estimated to be “normal”. The traveling place is estimated to be a “precursor” from the time when the brake differential value β> βT at time t0 to the time when the brake depression angle B> BT at time t1. Further, the traveling place is estimated to be “near hat” until the brake depression angle B> BT at time t1 and the brake depression angle B <BT at time t2. Further, the traveling place is estimated to be “precursor” until the brake depression angle B <BT at time t2 and the brake depression angle B = 0 at time t3. The traveling field is estimated to be “normal” from when the brake depression angle B = 0 at time t3 until the brake differential value β> βT again at time t4.

また、時刻t4に再びブレーキ微分値β>βとなってから、時刻t5にブレーキ踏込量B=0となるまで、走行場は「前兆」と推定される。そして、時刻t5にブレーキ踏込量B=0となり、走行場は「通常」と推定される。このようにして、ブレーキ踏込角B及びブレーキ微分値βに基づいて、ドライバの危機感の走行場を推定することができる。   Further, the traveling place is estimated to be a “precursor” until the brake differential value β> β again at time t4 until the brake depression amount B = 0 at time t5. At time t5, the brake depression amount B becomes 0, and the traveling place is estimated to be “normal”. In this way, it is possible to estimate the driving field of the driver's sense of crisis based on the brake depression angle B and the brake differential value β.

なお、上記の自車挙動、車両関係及びドライバの危機感についての走行場の他に、道路形状に関する走行場を推定してもよい。一般に、道路形状は、道路の曲率及び勾配によって表される。このため、道路形状に関する走行場として、曲率について「直線」、「右カーブ」及び「左カーブ」の3つの走行場を推定し、勾配について「平坦」、「上り坂」及び「下り坂」の3つの走行場を推定する。   In addition to the above-described driving field regarding the vehicle behavior, vehicle relations, and driver's sense of crisis, a driving field related to the road shape may be estimated. In general, the road shape is represented by the curvature and slope of the road. For this reason, as the driving field related to the road shape, three driving fields of “straight line”, “right curve” and “left curve” are estimated for the curvature, and “flat”, “uphill” and “downhill” for the gradient. Estimate three driving fields.

道路形状に関する走行場は、全地球測位システム(GPS:Global Positioning System)を利用して得た自車両の現在位置と、カーナビゲーションシステムの地図データベースとの照合によって推定するとよい。また、道路形状の走行場は、例えば、所定の走行区間における勾配、勾配の組合せ、又は勾配と曲率の組合せとして推定するとよい。さらに、例えば、自車挙動、車両関係、ドライバの危機感、及び道路形状についての走行場の打ちの任意の2つ以上を組み合わせて、走行場として推定するようにしてもよい。   The driving field related to the road shape may be estimated by comparing the current position of the host vehicle obtained by using a global positioning system (GPS) with a map database of the car navigation system. Further, the road-shaped traveling place may be estimated as, for example, a gradient in a predetermined traveling section, a combination of gradients, or a combination of gradient and curvature. Further, for example, the traveling place may be estimated by combining any two or more of the own vehicle behavior, the vehicle relation, the driver's sense of crisis, and the hitting of the running place with respect to the road shape.

(1−3.走行場の判定)
次に、判定部40が、推定部11によって推定された推定走行場が、第1記憶部30に記憶された既知の走行場から想定される想定内の走行場であるか、想定外の走行場であるかを判定する(図2のS23)。
(1-3. Judgment of running place)
Next, the determination unit 40 determines whether the estimated traveling place estimated by the estimating unit 11 is an expected traveling place that is assumed from the known traveling places stored in the first storage unit 30 or an unexpected traveling. It is determined whether it is a place (S23 in FIG. 2).

(1−3−1.制御アルゴリズム)
第1記憶部30には、過去の経験に基づく既知の走行場と車両の制御アルゴリズムとが互いに対応付けて記憶されている。本実施形態では、第1記憶部30に格納された、既知の走行場に対する制御アルゴリズムは、遺伝的ネットワークプログラミング(GNP:Genetic Network Programming)を利用して創発されたものである。
(1-3-1. Control algorithm)
The first storage unit 30 stores a known driving field based on past experience and a vehicle control algorithm in association with each other. In the present embodiment, the control algorithm for the known driving field stored in the first storage unit 30 is created using genetic network programming (GNP).

以下、まず、遺伝的ネットワークプログラミングについて説明し、次いで、その遺伝的ネットワークプログラミングによって創発された制御アルゴリズムの例を説明する。   Hereinafter, first, genetic network programming will be described, and then an example of a control algorithm created by the genetic network programming will be described.

遺伝的ネットワークプログラミングは、生物の進化過程を模倣して最適解又は準最適解を求める進化論的計算手法の一つであって、遺伝的アルゴリズム(GA:Genetic Algorithm)を拡張した遺伝的プログラミング(GP:Genetic Programming )を更に拡張したものである。   Genetic network programming is an evolutionary calculation method that seeks an optimal solution or a sub-optimal solution by imitating the evolution process of a living organism, and is a genetic programming (GP) that extends a genetic algorithm (GA). : Genetic Programming).

(1−3−1−1.遺伝的アルゴリズム(GA))
遺伝的アルゴリズム(GA)は、生物の進化過程において見られる、染色体の選択、交差及び突然変異を模倣し、工学的に応用した計算手法である。遺伝的アルゴリズムでは、設計変数を遺伝子とみなし、遺伝子をビット列構造で表現する。
(1-3-1-1. Genetic Algorithm (GA))
A genetic algorithm (GA) is a computational technique that mimics the selection, crossover and mutation of chromosomes found in the evolutionary process of organisms and is applied engineeringly. In the genetic algorithm, a design variable is regarded as a gene, and the gene is expressed by a bit string structure.

ここで、図8のフローチャートを参照して、遺伝的アルゴリズムの基本的なアルゴリズムの流れを説明する。
先ず、初期集団を作成する(S81)。初期集団は、個体の集合であり、各個体は、解候補を記号化したものであって、それぞれが一つの遺伝子を有する。
Here, the basic algorithm flow of the genetic algorithm will be described with reference to the flowchart of FIG.
First, an initial group is created (S81). The initial population is a set of individuals, and each individual is a solution candidate symbolized, each having one gene.

次いで、各個体における評価値を計算する(S82)。
評価値とは、各個体の有効性を定量的に示す指標をいい、例えば、燃費の良さを評価値とすることができる。
Next, an evaluation value for each individual is calculated (S82).
The evaluation value refers to an index that quantitatively indicates the effectiveness of each individual. For example, good fuel efficiency can be used as the evaluation value.

次いで、評価値に基づいて個体を選択する(S83)。
個体の選択とは、次世代に残す個体を選出することをいう。
Next, an individual is selected based on the evaluation value (S83).
Individual selection refers to selecting individuals to be left in the next generation.

次いで、選択された個体の交叉及び突然変異を行う(S84)。
交叉とは、親世代の遺伝子を交換し、子世代として新しい遺伝子を作成することをいう。また、突然変異とは、親世代に遺伝子の一部をランダムに書き換え、子世代とすることをいう。
Next, the selected individual is crossed and mutated (S84).
Crossing means exchanging the genes of the parent generation and creating a new gene as a child generation. Mutation means that a part of a gene is randomly rewritten to a parent generation to become a child generation.

次いで、最終世代の判定を行う(S85)。
最終世代になるまで、上記のS82〜S84のステップを繰り返す。
このようにして、最適解又は準最適解を探索する。
Next, the final generation is determined (S85).
The above steps S82 to S84 are repeated until the final generation.
In this way, an optimal solution or a sub-optimal solution is searched.

(1−3−1−2.遺伝的プログラミング(GP))
遺伝的プログラミング(GP)のアルゴリズムの流れも、図8に示したものと基本的に同じであるが、遺伝的プログラミング(GP)では、遺伝子を、ノード関数を用いた木構造で表現する。そして、ノード関数を用いてif−then文による条件判断を行い、処理を実行する。これにより、遺伝的プログラミングでは(GP)では、遺伝的アルゴリズム(GA)のビット列構造では扱えなかった構造的なシステムを扱うことができる。
(1-3-1-2. Genetic programming (GP))
The flow of the algorithm of genetic programming (GP) is basically the same as that shown in FIG. 8, but in genetic programming (GP), a gene is expressed by a tree structure using a node function. Then, the node function is used to determine the condition based on an if-then statement, and the process is executed. Thereby, in genetic programming (GP), it is possible to handle a structural system that cannot be handled by the bit string structure of the genetic algorithm (GA).

(1−3−1−3.遺伝的ネットワークプログラミング(GNP))
遺伝的ネットワークプログラミング(GNP)のアルゴリズムの流れも、図8に示したものと基本的に同じであるが、遺伝的ネットワークプログラミング(GNP)では、図9に示すように、集団を構成する各個体が、ノード関数が有向枝によって相互に接続された有向グラフ(ネットワーク)構造を有する。ノード関数は、判定ノードと処理ノードの二種類のノードから構成される。判定ノードは、if−then文による条件判断及び分岐を担う。一方、処理ノードは、定められた処理(行動)を実行する役割を有する。また、判定ノードからは、分岐の数だけ有向枝が伸びる。一方、処理ノードからは、一本の有向枝のみが伸びる。
なお、遺伝的ネットワークプログラミング(GNP)では、どの個体でもノード数が一定であるため、ブロート(解の膨張)が発生しない。
(1-3-1-3. Genetic Network Programming (GNP))
The flow of genetic network programming (GNP) algorithm is basically the same as that shown in FIG. 8, but in genetic network programming (GNP), as shown in FIG. Have a directed graph (network) structure in which node functions are connected to each other by directed edges. The node function is composed of two types of nodes, a determination node and a processing node. The decision node is responsible for condition judgment and branching by an if-then sentence. On the other hand, the processing node has a role of executing a predetermined process (action). Further, the directional branch extends from the determination node by the number of branches. On the other hand, only one directional branch extends from the processing node.
In Genetic Network Programming (GNP), since the number of nodes is constant in any individual, bloat (solution expansion) does not occur.

さらに、図10に、遺伝的ネットワークプログラミング(GNP)の一個体の遺伝子の一例を示す。遺伝子は、図10に示す表現型のようにネットワーク状に表現され、プログラム上では、図10に示す遺伝子型のように記述される。遺伝子型の「NIDi」は、ノードiのIDであり、設定したノードの数だけ存在する。また、遺伝子型の「NTi」は、ノード関数を表し、「NTi=1」ならば判定ノードを表し、「NTi=2」ならば処理ノードを表す。遺伝子型の「IDi」は、ノード関数IDを表し、判定ノード及び処理ノード内容を記述したノートライブラリのラベルを示す。そして、「NTi」と「IDi」によって、ノートライブラリに記述されたノード内容のうち、実行されるノード内容が決定される。遺伝子型の「Ci1」及び「Ci2」は、ノードiから伸びる有向枝の接続先情報であり、接続するノードIDが記載される。なお、処理ノードは、接続先が一つしかないため、処理ノードにCi2は存在しない。   Further, FIG. 10 shows an example of an individual gene of genetic network programming (GNP). The gene is expressed in a network form like the phenotype shown in FIG. 10, and is described as the genotype shown in FIG. 10 on the program. The genotype “NIDi” is the ID of the node i, and there are as many as the set number of nodes. The genotype “NTi” represents a node function, and “NTi = 1” represents a determination node, and “NTi = 2” represents a processing node. The genotype “IDi” represents a node function ID, and indicates a label of a note library that describes the contents of a determination node and a processing node. Of the node contents described in the note library, the node contents to be executed are determined by “NTi” and “IDi”. The genotypes “Ci1” and “Ci2” are connection destination information of a directional branch extending from the node i, and a node ID to be connected is described. Since the processing node has only one connection destination, Ci2 does not exist in the processing node.

このような各個体の評価値は、個体の達成度や有効性を数値化したものであり、その値は、評価換算によって求められる。評価関数の計算式は、問題に依存して設計者によって設定される。そして、個体の評価値に基づいて遺伝子操作が行われるので、遺伝的ネットワークプログラミング(GNP)の最適化は、評価値の最適化に対応する。   Such an evaluation value of each individual is obtained by quantifying the achievement level and effectiveness of the individual, and the value is obtained by evaluation conversion. The calculation formula of the evaluation function is set by the designer depending on the problem. Since genetic manipulation is performed based on the evaluation value of the individual, optimization of genetic network programming (GNP) corresponds to optimization of the evaluation value.

そして、交叉や突然変異を実行する前に、選択により、評価値の高い個体を優先的に選出する。これにより、評価値の高い個体の遺伝子が次世代へ受け継がれる。選択方式の例として、ルーレット選択、トーナメント選択及びエリート保存選択が挙げられる。   Then, before performing crossover or mutation, an individual with a high evaluation value is preferentially selected by selection. Thereby, the gene of an individual with a high evaluation value is inherited to the next generation. Examples of selection methods include roulette selection, tournament selection, and elite preservation selection.

ルーレット選択では、交叉及び突然変異操作を加える個体を選出するため、各個体の子孫は、その評価値に比例した確率で選ばれる。i番目の個体の評価値をfiとすると、下記の(1)式により、i番目の個体が選ばれる確率piが求まる。
pi=fi/Σf ・・・(1)
In roulette selection, individuals to be subjected to crossover and mutation operations are selected, and therefore the offspring of each individual is selected with a probability proportional to the evaluation value. When the evaluation value of the i-th individual is fi, the probability pi that the i-th individual is selected is obtained by the following equation (1).
pi = fi / Σf (1)

トーナメント選択では、交叉及び突然変異操作を加える個体を選出するため、予め設計者によって決められた数だけ個体集団の中からランダムで個体を取り出し、その中で最も適応度の高い個体を選択する。これを必要な個体数が得られるまで繰り返す。   In the tournament selection, in order to select individuals to be subjected to crossover and mutation operations, individuals are randomly selected from the population of individuals determined in advance by the designer, and the individual with the highest fitness is selected. This is repeated until the required number of individuals is obtained.

エリート保存選択では、その世代における評価値の高い個体を一定数そのまま次世代へ残す操作を行う。これにより、最適個体の評価が悪化することを防ぐことができる。   In the elite preservation selection, an operation is performed to leave a certain number of individuals with high evaluation values in the generation as they are to the next generation. Thereby, it is possible to prevent the evaluation of the optimum individual from deteriorating.

交叉は、選択された2個の個体間で行われ、2個の子個体を生成する。このとき、交叉するノードどうしの遺伝子情報は全て交換される。ここで、図11を参照して、交叉の手順について説明する。   Crossover is performed between the two selected individuals to generate two child individuals. At this time, all the gene information between the intersecting nodes is exchanged. Here, the crossover procedure will be described with reference to FIG.

先ず、親個体の集団から2個の親個体を任意に選択する。図11では、「親1」及び「親2」が選択される。次いで、選択した親個体の一方「親1」において、設定した交叉確率に基づいて、交叉すべきノードを選択する。図11では、「親1」のノード「P1」が選択される。次いで、一方の親個体「親1」で選択されたノード「P1」と、他方の親個体「親2」における「P1」と同一識別番号のノード「P2」との間で、全ての遺伝子を交換する。このようにして生成された2個の個体「子1」及び「子2」が、次世代の個体となる。   First, two parent individuals are arbitrarily selected from a group of parent individuals. In FIG. 11, “parent 1” and “parent 2” are selected. Next, in one “parent 1” of the selected parent individuals, a node to be crossed is selected based on the set crossing probability. In FIG. 11, the node “P1” of “parent 1” is selected. Next, all the genes are transferred between the node “P1” selected in one parent individual “parent 1” and the node “P2” having the same identification number as “P1” in the other parent individual “parent 2”. Exchange. The two individuals “child 1” and “child 2” generated in this way become the next generation individuals.

突然変異は、1個の親個体内で行われ、1個の子個体が生成される。突然変異には、接続先ノードを変更するものと、ノード内容を変更するものとの2種類がある。ここで、図12を参照して、突然変異について説明する。
まず、接続先ノードを変更する突然変異(突然変異A)では、1個体を任意に選択し、設定した突然変異確率に基づき、突然変異させるべきノードを選択する。そして、選択されたノードの接続先をランダムに変更する。図12では、突然変異Aにおいて、ノードP1の接続先ノードが、右上の処理ノードから、左上の判断ノードへ切り替えられている。
Mutation is performed within one parent individual to generate one child individual. There are two types of mutations: one that changes the connection destination node and one that changes the node contents. Here, the mutation will be described with reference to FIG.
First, in the mutation for changing the connection destination node (mutation A), one individual is arbitrarily selected, and the node to be mutated is selected based on the set mutation probability. Then, the connection destination of the selected node is randomly changed. In FIG. 12, in the mutation A, the connection destination node of the node P1 is switched from the upper right processing node to the upper left determination node.

また、ノード内容を変更する突然変異(突然変異B)では、1個体を任意に選択し、設定した突然変異確率に基づき、突然変異させるべきノードを選択する。そして、選択されたノードのラベルをランダムに変更する。図12では、突然変異Bにおいて、選択されたノードP1がノードP2に変更されている。
このようにして、次世代の個体が生成される。
Further, in the mutation (mutation B) for changing the node contents, one individual is arbitrarily selected, and the node to be mutated is selected based on the set mutation probability. Then, the label of the selected node is randomly changed. In FIG. 12, in the mutation B, the selected node P1 is changed to the node P2.
In this way, a next generation individual is generated.

(1−3−1−4.遺伝的ネットワークプログラミング(GNP)の利用例)
以下に、遺伝子ネットワークプログラミング(GNP)を利用して、ハイブリッド電気自動車におけるエンジンとモータの分担に関する車両制御アルゴリズムを創発する例を示す。ここでは、走行場に応じて燃費が改善される車両制御アルゴリズムを創発する。
(1-3-1-4. Usage example of genetic network programming (GNP))
An example of generating a vehicle control algorithm related to the sharing of an engine and a motor in a hybrid electric vehicle by using gene network programming (GNP) will be described below. Here, a vehicle control algorithm that improves fuel efficiency in accordance with the driving field is created.

まず、制御アルゴリズムの前提となる、シリーズ方式ハイブリッド電気自動車の車両システムを図13に示す。また、車両諸元のパラメータを図14に示す。
また、重力加速度g=9.8m/s2と設定し、空気密度ρ=1.29kg/m3と設定する。
First, FIG. 13 shows a vehicle system of a series-type hybrid electric vehicle, which is a premise of the control algorithm. Moreover, the parameter of a vehicle specification is shown in FIG.
Further, the gravitational acceleration g is set to 9.8 m / s 2 and the air density ρ = 1.29 kg / m 3 is set.

モータ出力WM(W)は、下記の式(2)により求める。
M=TM・(V/rD)i ・・・(2)
ここで、TM[Nm]は、モータトルクを表し、V[m/s]は車速を表す。モータの最大出力WMmaxとモータトルクTMとの関係より、車速V=18.75m/sまでは最大トルクが400Nmであり、これ以上の速度では最大トルクは速度に反比例して低下していく。なお、モータトルクは最大トルクと正規化されたアクセル入力量との積から求める。
The motor output W M (W) is obtained by the following equation (2).
W M = T M · (V / r D ) i (2)
Here, T M [Nm] represents motor torque, and V [m / s] represents vehicle speed. From the relationship between the maximum motor output W Mmax and the motor torque T M , the maximum torque is 400 Nm up to the vehicle speed V = 18.75 m / s, and the maximum torque decreases in inverse proportion to the speed above this. . The motor torque is obtained from the product of the maximum torque and the normalized accelerator input amount.

図13に示すように、シリーズ方式のハイブリッド電気自動車において、エンジンの動力は発電機で電力に変換される。発電量WG[W]は、下記の(3)式より求める。
G[W]=ηGE ・・・(3)
ただし、WE[W]は、エンジン出力である。また、燃料消費率FE[g/kWh]は、図15のグラフにより求める。図15のグラフは、エンジン出力と燃料消費率との関係を示す。さらに、1秒当たりの燃料消費量FC[g/s]は、下記の式(4)により求める。
C=WE・FE/1000/3600 ・・・(4)
As shown in FIG. 13, in a series-type hybrid electric vehicle, engine power is converted into electric power by a generator. The power generation amount W G [W] is obtained from the following equation (3).
W G [W] = η G W E (3)
However, W E [W] is an engine output. Further, the fuel consumption rate F E [g / kWh] is obtained from the graph of FIG. The graph of FIG. 15 shows the relationship between engine output and fuel consumption rate. Further, the fuel consumption amount F C [g / s] per second is obtained by the following equation (4).
F C = W E · F E / 1000/3600 (4)

また、バッテリ構成を図16に示す。また、バッテリ内部電圧VB=400Vと設定し、バッテリ内部抵抗RB=400mVと設定する。 The battery configuration is shown in FIG. Further, the battery internal voltage V B = 400V is set, and the battery internal resistance R B = 400 mV is set.

バッテリ残容量EB[Wh]は、下記の式(5)より求める。
dEB/dt=WB−(WB/VB2・RB
B=WG−WM−WA ・・・(5)
ただし、WBは、バッテリ入出力である。
The remaining battery capacity E B [Wh] is obtained from the following equation (5).
dE B / dt = W B − (W B / V B ) 2 · R B
W B = W G −W M −W A (5)
However, W B is a battery input and output.

また、等価燃料消費量FB(g)は、下記の式(6)より求める。
B=ΔSOC・EBmaxE/ηM/1000 ・・・(6)
等価燃料消費量とは、走行終了後にエンジンを稼働させて、バッテリ残量SOCを書おきの残量に戻すために必要とする燃料消費量である。その際、エンジンは最高効率点FE=250g/kWhで稼働させるものとする。
ただし、ΔSOC(%)は、走行終了時のバッテリ残量SOCと、走行開始時のバッテリ残量SOCとの差である。
The equivalent fuel consumption F B (g) is obtained from the following equation (6).
F B = ΔSOC · E Bmax F E / η M / 1000 (6)
The equivalent fuel consumption amount is a fuel consumption amount required for operating the engine after the travel is completed and returning the remaining battery charge SOC to the written remaining charge. In this case, the engine is operated at the maximum efficiency point F E = 250 g / kWh.
However, ΔSOC (%) is the difference between the remaining battery charge SOC at the end of travel and the remaining battery charge SOC at the start of travel.

また、ブレーキは、回生ブレーキ、油圧ブレーキ及びエンジンブレーキから構成される。回生ブレーキは、ドライバから要求される要求ブレーキトルクが最大モータトルクTMmax以下、すなわち、400Nm以下の場合に用いられる。回生ブレーキによって得られる回生エネルギーの大きさは、上記の式(2)より求まる。回生エネルギーとモータの発電効率ηMとの積がモータの発電量となる。このとき、上記の式(5)におけるモータ出力WMは負となる。 The brake includes a regenerative brake, a hydraulic brake, and an engine brake. The regenerative brake is used when the required brake torque requested by the driver is not more than the maximum motor torque T Mmax , that is, not more than 400 Nm. The magnitude of the regenerative energy obtained by the regenerative brake can be obtained from the above equation (2). The product of the regenerative energy and the power generation efficiency η M of the motor is the power generation amount of the motor. At this time, the motor output W M in the above equation (5) is negative.

油圧ブレーキは、要求ブレーキトルクが回生ブレーキで得られるトルクを超過するとき、すなわち、400Nm以上の場合に用いられ、最大モータトルクTMmaxの超過分を補う。なお、要求ブレーキトルクは、最大ブレーキトルクと正規化されたブレーキ入力量との積から求める。 The hydraulic brake is used when the required brake torque exceeds the torque obtained by the regenerative brake, that is, when it is 400 Nm or more, and compensates for the excess of the maximum motor torque T Mmax . The required brake torque is obtained from the product of the maximum brake torque and the normalized brake input amount.

エンジンブレーキは、ブレーキとアクセルが共に入力されない場合に用いられる。シリーズ方式ハイブリッド電気自動車では、駆動軸とエンジンとが直接連結していないため、本来エンジンブレーキは存在しない。しかし、回生ブレーキトルク35Nmを与えることでそれを再現し、充電走行を行う。   The engine brake is used when neither the brake nor the accelerator is input. In a series-type hybrid electric vehicle, the drive shaft and the engine are not directly connected, and therefore there is essentially no engine brake. However, regenerative braking torque of 35 Nm is applied to reproduce it, and charging is performed.

また、車体の挙動は、下記の式(7)より求める。
m(dV/dt)={(TMiηt−TB)/rD}−R ・・・(7)
R=Ri+Rr+Rl
i=mgsinθ
r=μrmg
l=(ρ/2)CDAV2
Further, the behavior of the vehicle body is obtained from the following equation (7).
m (dV / dt) = { (T M iη t -T B) / r D} -R ··· (7)
R = R i + R r + R l
R i = mg sin θ
R r = μ r mg
R l = (ρ / 2) C D AV 2

ただし、また、m[kg]は車体重量、TMはモータトルク[Nm]、iは減速比、ηtは伝達効率、TBはブレーキトルク[Nm]を表す。また、R(N)は走行抵抗、Ri[N]は勾配抵抗、Rr[N]は転がり抵抗、Rl[N]は空気抵抗、θ[rad]は道路勾配を表す。また、μrは転がり摩擦係数、CDは抵抗係数、Aは車両の前面投影面積を表す。 However, m [kg] represents the weight of the vehicle body, T M represents the motor torque [Nm], i represents the reduction ratio, η t represents the transmission efficiency, and T B represents the brake torque [Nm]. R (N) represents running resistance, Ri [N] represents gradient resistance, Rr [N] represents rolling resistance, Rl [N] represents air resistance, and θ [rad] represents road gradient. Further, mu r is the rolling friction coefficient, C D is the drag coefficient, A represents the frontal projected area of the vehicle.

遺伝的ネットワークプログラミング(GNP)の設定パラメータを下記の表4に示し、判定ノード及び処理ノードの設定を下記の表5に示す。判定ノード及び処理ノードは、入力値の範囲を任意の間隔で刻み、閾値の異なるノードを複数生成する。判定及び処理の債には、入力値のその閾値の範囲内か否かで「yes」又は「No」の二方向出力を行う。   Genetic network programming (GNP) setting parameters are shown in Table 4 below, and determination node and processing node settings are shown in Table 5 below. The determination node and the processing node cut a range of input values at arbitrary intervals and generate a plurality of nodes having different threshold values. The determination and processing bonds are output in two directions “yes” or “No” depending on whether or not the input value is within the range of the threshold value.

Figure 0005665400
Figure 0005665400

Figure 0005665400
Figure 0005665400

次に、図17のフローチャートを参照して、上述した前提の下での遺伝的ネットワークプログラミング(GNP)の流れを説明する。
まず、ランダム初期集団として、判定ノード及び処理ノードがランダムに組み合わされたネットワークを個体の数だけ作成する(S171)。
Next, the flow of genetic network programming (GNP) under the above-described assumption will be described with reference to the flowchart of FIG.
First, as a random initial group, a network in which determination nodes and processing nodes are randomly combined is created by the number of individuals (S171).

次に、走行場ごとのアクセル踏込角及び勾配を入力する(S172)。アクセル踏込量及び勾配は、例えば、車両関係の走行場が「自由」であり、かつ、道路勾配の走行場が「平坦」、「上り坂」又は「下り坂」それぞれの場合についての入力データを使用するとよい。なお、道路勾配が「上り坂」又は「下り坂」の場合の勾配は、4.0%とする。
次に、アクセル踏込角の入力から、要求されるモータ出力及びブレーキトルクを計算する(S173)。
Next, the accelerator depression angle and gradient for each traveling place are input (S172). The accelerator depression amount and the gradient are, for example, input data when the vehicle-related driving field is “free” and the road gradient driving field is “flat”, “uphill” or “downhill”. It is good to use. The slope when the road slope is “uphill” or “downhill” is 4.0%.
Next, the required motor output and brake torque are calculated from the input of the accelerator depression angle (S173).

次に、バッテリ残量SOCを所定の範囲内(40〜70%)に維持するための条件判断を行う(S174)。バッテリ残量SOCが70%を超過した場合(S174で「yes」の場合)に、エンジンを停止させ(S175)る。また、バッテリ残量SOCが40%を下回った場合(S176で「yes」の場合)に、エンジンを稼働させる(S177)。エンジンを稼働させる際に、要求されるエンジン出力が、エンジンの最適出力35kW以下であれば、35kWを出力し、35kW以上であれば、要求通りの値を出力する(S177)。   Next, a condition determination for maintaining the remaining battery charge SOC within a predetermined range (40 to 70%) is performed (S174). When the remaining battery SOC exceeds 70% (in the case of “yes” in S174), the engine is stopped (S175). Further, when the remaining battery charge SOC is less than 40% (in the case of “yes” in S176), the engine is operated (S177). When operating the engine, if the required engine output is 35 kW or less, the engine output is 35 kW, and if it is 35 kW or more, the requested value is output (S177).

また、バッテリ残量SOCが所定の範囲内(40〜70%)の場合(S174で「no」かつS176で「no」の場合)、遺伝的ネットワークプログラミング(GNP)による制御アルゴリズムによってエンジン出力を決定する(S178)。   When the remaining battery SOC is within a predetermined range (40 to 70%) (“no” in S174 and “no” in S176), the engine output is determined by a control algorithm based on genetic network programming (GNP). (S178).

続いて、発電量を計算し、燃料消費量を計算する(S179)。さらに、バッテリ残量SOCを計算し(S180)、運動方程式を解く(S181)。そして、上記のステップS172〜S181のループを、所定の走行場における走行モードが終了するまで(S182で「no」の間)繰り返す。   Subsequently, the power generation amount is calculated, and the fuel consumption amount is calculated (S179). Further, the remaining battery charge SOC is calculated (S180), and the equation of motion is solved (S181). And the loop of said step S172-S181 is repeated until the driving | running | working mode in a predetermined driving field is complete | finished (it is "no" in S182).

走行モードの終了後(S182で「yes」の場合)に、ネットワーク(1個体)の評価値を計算する(S183)。評価値は、燃料消費量と等価燃料消費量の和として表される。上記のステップS172〜S183のループを、個体の数繰り返す。   After the travel mode ends (in the case of “yes” in S182), an evaluation value of the network (one individual) is calculated (S183). The evaluation value is expressed as the sum of the fuel consumption and the equivalent fuel consumption. The loop of steps S172 to S183 is repeated for the number of individuals.

そして、全ての個体の評価値が求められた後(S184で「yes」の場合)、最終世代であるか否かの判断を行う(S185)。最終世代でない場合(S185で「no」の場合)、個体集団に対して、表4に示した設定パラメータの確率に従って、選択、交叉、突然変異といった遺伝子操作を加えて進化を行う(S186)。また、最終世代である場合(S185で「yes」の場合)、最適個体のノード遷移ネットワークを保存して(S187)、終了する。   Then, after the evaluation values of all the individuals are obtained (in the case of “yes” in S184), it is determined whether or not it is the last generation (S185). If it is not the final generation (in the case of “no” in S185), the individual population is evolved by applying genetic operations such as selection, crossover, and mutation according to the probability of the setting parameter shown in Table 4 (S186). If it is the final generation (in the case of “yes” in S185), the node transition network of the optimum individual is stored (S187), and the process ends.

このような最適個体のノード遷移ネットワークによる制御アルゴリズムが、第1記憶部30に走行場と対応付けて記憶されている。   A control algorithm based on such a node transition network of the optimum individual is stored in the first storage unit 30 in association with the traveling place.

(1−4.推定場が想定内の場合)
次に、判定部40によって推定走行場が、第1記憶部30に記憶された既知の走行場から想定される想定内の走行場であると判定された場合(図2のS23で「yes」の場合)に、変換部60が、上記第1記憶手段から推定走行場に対応する制御アルゴリズムを読み出し(S24)、車両の制御パラメータに変換する(S25)。
(1-4. When the estimated field is within the expected range)
Next, when it is determined by the determination unit 40 that the estimated traveling place is an expected traveling place assumed from the known traveling places stored in the first storage unit 30 ("yes" in S23 of FIG. 2). In this case, the conversion unit 60 reads out a control algorithm corresponding to the estimated traveling place from the first storage means (S24), and converts it into a vehicle control parameter (S25).

例えば、図13に示したシステムのシリーズ方式ハイブリッド電気自動車について、車両関係の走行場が「自由」であり、かつ道路勾配の走行場が「平坦」であるように、想定範囲内の走行場である場合には、第1記憶部30から、その走行場に対応する制御アルゴリズムが読み出される。そして、読み出された制御アルゴリズムに従って、車両の制御パラメータとして、エンジンとモータの分担出力が計算される。   For example, in the series-type hybrid electric vehicle of the system shown in FIG. 13, the vehicle-related driving field is “free” and the road gradient driving field is “flat”. In some cases, a control algorithm corresponding to the traveling place is read from the first storage unit 30. Then, according to the read control algorithm, the shared output between the engine and the motor is calculated as the vehicle control parameter.

また、制御パラメータにしたがって、エンジン及びモータ以外に、ハンドルやブレーキといったアクチュエータを制御してもよい。例えば、制御アルゴリズムから変換された舵角に従って、ハンドルを自動操作するようにしてもよいし、制御アルゴリズムから変換された車速に従って、ブレーキを自動的に作動させてもよい。   In addition to the engine and motor, actuators such as handles and brakes may be controlled according to the control parameters. For example, the steering wheel may be automatically operated according to the steering angle converted from the control algorithm, or the brake may be automatically operated according to the vehicle speed converted from the control algorithm.

また、制御パラメータに従って、また、スピーカやディスプレイ等のヒューマン・マシン・インターフェイス(human machine interface:HMI)を制御してもよい。例えば、制御アルゴリズムから変換された舵角や車速に従って車両が走行するように、ドライバを音声で誘導するようにしてもよい。   Further, according to the control parameter, a human machine interface (HMI) such as a speaker or a display may be controlled. For example, the driver may be guided by voice so that the vehicle travels according to the steering angle and the vehicle speed converted from the control algorithm.

(1−5.推定場が想定外の場合)
一方、判定部40によって推定走行場が、第1記憶部30に記憶された既知の走行場から想定される想定内の走行場ではない、想定外の走行場であると判定された場合(図2のS23で「no」の場合)には、創発部50が、推定走行場に対応する制御アルゴリズムを創発する(S26)。
(1-5. When the estimated field is unexpected)
On the other hand, when it is determined by the determination unit 40 that the estimated traveling place is an unexpected traveling place that is not an assumed traveling place assumed from the known traveling places stored in the first storage unit 30 (FIG. 2 is “no” in S23), the emergence unit 50 creates a control algorithm corresponding to the estimated traveling place (S26).

図18のフローチャートを参照して、推定された走行場が想定外である場合の制御アルゴリズムの創発処理を説明する。図18のフローチャートは、図2のフローチャートのステップS23、S26及びS27の部分に相当する。ここでは、走行場として、図19に示すように、自車両Cの進行方向前方に、障害物として事故車両Bが停車している場合について説明する。   With reference to the flowchart of FIG. 18, the emergence process of the control algorithm in the case where the estimated traveling place is unexpected is described. The flowchart in FIG. 18 corresponds to steps S23, S26, and S27 of the flowchart in FIG. Here, the case where the accident vehicle B is stopped as an obstacle in front of the traveling direction of the host vehicle C as shown in FIG.

先ず、推定部11によって推定された走行場が、想定範囲内の走行場か否かが判定される(図18のS23)。
ここでは、走行場は、図20に示すように、事故車両の側方に残された「幅員」と、自車両が事故車両の手前で停止するのに必要な「減速度」とに基づいて推定される。「幅員」は、車外環境モニタ12により撮像された前方画像から求められる。また、「減速度」は、車外環境モニタ12により撮像された前方画像から求めた自車両の現在位置から事故車両までの距離と、車両状態モニタ10によって測定された現在の自車両の車速とから求められる。
First, it is determined whether or not the traveling place estimated by the estimating unit 11 is a traveling place within an assumed range (S23 in FIG. 18).
Here, as shown in FIG. 20, the traveling place is based on “width” left on the side of the accident vehicle and “deceleration” necessary for the own vehicle to stop before the accident vehicle. Presumed. The “width” is obtained from the front image captured by the external environment monitor 12. The “deceleration” is calculated from the distance from the current position of the host vehicle to the accident vehicle obtained from the front image captured by the outside environment monitor 12 and the current vehicle speed measured by the vehicle state monitor 10. Desired.

一方、第1記憶部30に記憶された既知の走行場も、「幅員」と「減速度」の2つのパラメータで特定されている。図20に、既知の走行場の幾つかを白丸印で模式的に示す。図20中の白丸Aは、「幅員」が十分にあり、かつ、事故車までの距離が長いため「減速度」が小さいという走行場を表している。かかる走行場に対しては、ハンドル操作により事故車両を大きく回避するという制御アルゴリズムが記憶されている。また、図20中の白丸Bは、「幅員」が十分に有り、かつ、事故車までの距離が短いため「減速度」が大きいという走行場を表している。かかる走行場に対しても、ハンドル操作による事故車を大きく回避するという制御アルゴリズムが記憶されている。また、図20中の白丸Cは、「幅員」が十分に無く、かつ、「減速度」が白丸Aと白丸Bとの中間の走行場を表している。かかる走行場に対しては、事故車の横を走行することができないため、ブレーキをかける制御アルゴリズムが記憶されている。この場合、更に、事故車との衝突に備えてハンドルをしっかり持って踏ん張るように音声でドライバに指示する制御アルゴリズムも記憶されているとよい。   On the other hand, the known driving field stored in the first storage unit 30 is also specified by two parameters of “width” and “deceleration”. FIG. 20 schematically shows some of the known running places with white circles. A white circle A in FIG. 20 represents a traveling place where “width” is sufficient and the “deceleration” is small because the distance to the accident vehicle is long. A control algorithm for largely avoiding an accident vehicle by a steering wheel operation is stored for such a driving field. Also, a white circle B in FIG. 20 represents a traveling place where “width” is sufficient and the “deceleration” is large because the distance to the accident vehicle is short. A control algorithm for largely avoiding an accident vehicle caused by a steering wheel operation is stored even for such a driving field. A white circle C in FIG. 20 indicates a traveling place where “width” is not enough and “deceleration” is between the white circle A and the white circle B. A control algorithm for applying a brake is stored for such a traveling place because the vehicle cannot travel beside the accident vehicle. In this case, it is also preferable to store a control algorithm that instructs the driver by voice to hold the steering wheel in preparation for a collision with the accident vehicle.

また、図20の全ての白丸印を囲む曲線Iは、想定内の走行場の範囲を表している。想定内の走行場の範囲には、各白丸印の既知の走行場そのものだけではなく、既知の走行場どうしの間の走行場、すなわち、既知の走行場どうしの間で内挿できる走行場も含まれる。したがって、推定された走行場が、第1記憶部30に記憶された既知の走行場どうしの間で内挿できる場合は想定内の走行場と判定され、既知の走行場から外挿が必要な場合は想定外の走行場と判定される。   In addition, a curve I surrounding all the white circles in FIG. 20 represents an expected range of the driving field. The range of the expected driving ranges includes not only the known driving points themselves indicated by the white circles but also the driving points between the known driving points, that is, the driving points that can be interpolated between the known driving points. included. Therefore, when the estimated traveling place can be interpolated between known traveling places stored in the first storage unit 30, it is determined as an assumed traveling place, and extrapolation from the known traveling place is necessary. In the case, it is determined as an unexpected traveling place.

さらに、図20中に、推定部11によって推定された走行場を黒丸印Dで示す。黒丸Dは、「幅員」が白丸Aと白丸Cとの中間であり、かつ、「減速度」が大きい走行場を表している。図20に示すように、推定された走行場は、曲線Iで囲まれた想定範囲内に含まれていない。このため、黒丸Dで示される推定走行場は、想定外の走行場であると判定される(図18のS23)。   Further, in FIG. 20, the traveling place estimated by the estimation unit 11 is indicated by a black circle D. A black circle D represents a traveling place where the “width” is between the white circle A and the white circle C and the “deceleration” is large. As shown in FIG. 20, the estimated traveling place is not included in the assumed range surrounded by the curve I. For this reason, it is determined that the estimated traveling place indicated by the black circle D is an unexpected traveling place (S23 in FIG. 18).

推定走行場が想定外であると判定された場合、第2記憶部70から、対応テーブルを読み込む(S261)。対応テーブルには、過去の経験に基づかない走行場と車両の制御アルゴリズムとを互いに対応付けて記憶されている。
ここで、下記の表6に、対応テーブルの一例を示す。
When it is determined that the estimated traveling place is unexpected, the correspondence table is read from the second storage unit 70 (S261). The correspondence table stores a driving field that is not based on past experience and a vehicle control algorithm in association with each other.
Here, an example of the correspondence table is shown in Table 6 below.

Figure 0005665400
Figure 0005665400

次に、対応テーブルから、推定走行場に最も近似する近似走行場を選択し、この近似走行場に対応する制御アルゴリズムを読み出す(S262)。
ここでは、上記の表6に示す想定が対応テーブルから、「前方衝突防止距離」が「無し」であり、「後方衝突防止距離」が「有り」であり、かつ、「側方スペース」が「有り」の場合を架設外挿点とする(S262)。
Next, from the correspondence table, an approximate traveling place that most closely approximates the estimated traveling place is selected, and a control algorithm corresponding to the approximate traveling place is read out (S262).
Here, the assumption shown in Table 6 above is that the “front collision prevention distance” is “none”, the “rear collision prevention distance” is “present”, and the “side space” is “ The “existing” case is set as a construction extrapolation point (S262).

次に、第1記憶手段に記憶された既知の走行場の中から、黒丸Dで示した推定走行場に最も近似する近似走行場を選択し、この近似走行場に対応する制御アルゴリズムを読み出す(S263)。
ここでは、図20に白丸Eで示す走行場を近似走行場として選択する。
Next, an approximate traveling place that is most approximate to the estimated traveling place indicated by the black circle D is selected from the known traveling places stored in the first storage means, and a control algorithm corresponding to this approximate traveling place is read ( S263).
Here, the traveling place indicated by the white circle E in FIG. 20 is selected as the approximate traveling place.

そして、想定内の走行場から選択した白丸Eの最近似の走行場と、対応テーブルから選択した最近似の走行場の中間の走行場を、推定場として選択し(S264)、白丸Eの走行場の制御アルゴリズムと対応テーブルから選択した走行場の制御アルゴリズムの中間の制御アルゴリズムを創発する(S265)。   Then, an intermediate driving field between the most approximate driving field of the white circle E selected from the assumed driving field and the closest driving field selected from the correspondence table is selected as the estimation field (S264), and the driving of the white circle E is performed. A control algorithm intermediate between the field control algorithm and the driving field control algorithm selected from the correspondence table is generated (S265).

ここで、図21を参照して、制御アルゴリズムの創発について説明する。
まず、幅員に余裕が有る場合(S211で「有り」の場合)、例えば、図20で想定内の白丸Bの走行場の場合、ハンドル操作によって障害物を回避する制御アルゴリズムが読み出される(S212)。さらに、ドライバに対して、徐行運転が指示される(S213)。
Here, the emergence of the control algorithm will be described with reference to FIG.
First, when there is a margin in the width (in the case of “Yes” in S211), for example, in the case of the white circle B that is assumed in FIG. 20, a control algorithm for avoiding an obstacle by a steering wheel operation is read (S212). . Further, the driver is instructed to drive slowly (S213).

また、幅員に余裕が無い場合(S211で「無し」の場合)、例えば、図20で想定内の白丸Cの走行場の場合、自車両が障害物と衝突することによって発生する被害を軽減するために、ブレーキを作動させる制御アルゴリズムが読み出され(S214)、さらに、エンジンを停止させる(S215)。   Further, when there is no margin in the width (in the case of “None” in S211), for example, in the case of the driving circle of the white circle C that is assumed in FIG. 20, the damage caused by the collision of the host vehicle with the obstacle is reduced. Therefore, a control algorithm for operating the brake is read (S214), and the engine is further stopped (S215).

そして、幅員の余裕が、「有り」と「無し」との中間である場合(S211で「中間」の場合)、例えば、図20で想定外の黒丸Dの走行場である場合、例えば、遺伝的ネットワークプログラミングを利用して、新たな制御アルゴリズムが創発される(S216)。   Then, when the margin of the width is intermediate between “present” and “absent” (in the case of “intermediate” in S211), for example, in the case of an unexpected black circle D in FIG. 20, for example, genetic A new control algorithm is created using dynamic network programming (S216).

ここで、図22に、新たな制御アルゴリズムを創発する際の、初期集団の個体の有向グラフ(ネットワーク)構造の一例を示す。初期集団の個体は、白丸Eの走行場のものと、対応テーブルから選択された走行場のものとするのがよい。ここでは、自車両が障害物と衝突することによって発生する被害の程度を最小にすることを目的として、被害度合で制御アルゴリズムが評価される。そして、白丸Eの走行場の有向グラフ構造と、対応テーブルから選択された走行場の有向グラフ構造について、交叉や突然変異といった遺伝子操作を繰り返して、新たな制御アルゴリズムを創発する。   Here, FIG. 22 shows an example of the directed graph (network) structure of the individuals in the initial group when a new control algorithm is created. The individuals of the initial group may be those of the driving field of the white circle E and those of the driving field selected from the correspondence table. Here, the control algorithm is evaluated based on the degree of damage for the purpose of minimizing the degree of damage that occurs when the host vehicle collides with an obstacle. Then, genetic operations such as crossover and mutation are repeated for the directed graph structure of the driving field of the white circle E and the directed graph structure of the driving field selected from the correspondence table to generate a new control algorithm.

そして、例えば、障害物である事故車両の側面をこすりながら自車両を停止させるという制御アルゴリズムが創発される。   Then, for example, a control algorithm is generated that stops the host vehicle while rubbing the side of the accident vehicle that is an obstacle.

そして、創発された制御アルゴリズムの良否をドライバの反応として学習させて(S266)、創発された制御アルゴリズムは第3記憶部80に記憶される。このように、車両が遭遇している走行場が自動的に推定され、車両用制御装置の制御アルゴリズムの設計段階で想定されていなかった走行場に車両が遭遇したときに、その走行場に対応した新たな制御アルゴリズムが自動的に生成される。そして、既存の制御アルゴリズムに加えて、想定外の走行場に対応した新たな制御アルゴリズムを蓄積することによって、より多様な走行場に対応して車両制御を行うことが可能となる。   Then, it learns the quality of the generated control algorithm as a response of the driver (S266), and the generated control algorithm is stored in the third storage unit 80. In this way, the place where the vehicle encounters is automatically estimated, and when the vehicle encounters a place that was not assumed at the design stage of the control algorithm of the vehicle control device, it corresponds to that place A new control algorithm is automatically generated. In addition to the existing control algorithms, by accumulating new control algorithms corresponding to unexpected driving fields, it becomes possible to perform vehicle control corresponding to more diverse driving fields.

(2.第2実施形態)
以下、本発明の第2実施形態を説明する。
第2実施形態による車両用制御装置の構成及び基本動作は、上述の第1実施形態のものと同じである。
ただし、第2実施形態では、推定部20は、車両の遭遇している走行場を、恒等写像モデルを利用して推定する。
(2. Second Embodiment)
Hereinafter, a second embodiment of the present invention will be described.
The configuration and basic operation of the vehicle control device according to the second embodiment are the same as those of the first embodiment described above.
However, in 2nd Embodiment, the estimation part 20 estimates the driving field which the vehicle has encountered using an identity mapping model.

図23に、恒等写像モデルの模式図を示す。恒等写像モデルは、入力層と出力層に同一データを提示して、恒等写像としての中間層出力を学習させる階層型ニューラルネットワークである。図23に示す例では、入力層及び出力層に、車外環境データとして、先行車両と自車両との車間距離を入力し、また、車両状態データとして、車速、ブレーキ踏込量、アクセル踏込量及び操舵角を入力している。   FIG. 23 shows a schematic diagram of the identity mapping model. The identity mapping model is a hierarchical neural network that presents the same data to an input layer and an output layer and learns an intermediate layer output as an identity map. In the example shown in FIG. 23, the inter-vehicle distance between the preceding vehicle and the host vehicle is input to the input layer and the output layer as external environment data, and the vehicle speed, the brake depression amount, the accelerator depression amount, and the steering are input as vehicle state data. The corner is entered.

図23に示す例では、中間層として、自動車挙動(停止、発進、走行)、車両関係(追従、自由、縮まる、拡がる、接近、衝突)、ドライバ反応(通常、前兆、ヒヤリハット)、及び道路勾配(平坦、上り坂、下り坂)を走行場として出力する。   In the example shown in FIG. 23, as an intermediate layer, car behavior (stop, start, running), vehicle relations (following, free, shrink, spread, approach, collision), driver reaction (usually, precursor, near-miss), and road gradient (Flat, Uphill, Downhill) is output as a driving field.

図24に、恒等写像モデルの中間層出力の模式図を示す。図24の横軸x1及び縦軸x2は、それぞれ中間層の出力パラメータである。そして、図24に示すように、中間層出力によって、走行場(例えば、追い越し、通中、縦列駐車、駐車、右折、交差点通過、左折)が特定される。   FIG. 24 shows a schematic diagram of the intermediate layer output of the identity mapping model. The horizontal axis x1 and the vertical axis x2 in FIG. 24 are output parameters of the intermediate layer, respectively. Then, as shown in FIG. 24, a traveling place (for example, overtaking, passing, parallel parking, parking, right turn, crossing an intersection, left turn) is specified by the intermediate layer output.

また、第2実施形態では、第1記憶部30に記憶されている走行場が中間層出力によって特定されている。図25の(A)及び(B)に、想定範囲の模式図を示す。図25(A)及び(B)の線IIは、第1記憶部30に記憶されている想定場の想定範囲を模式的に示す。図25の(A)及び(B)の横軸及び縦軸は、いずれも中間層出力のパラメータである。そして、図25(A)では、破線で示す推定走行場e1が、線IIで囲まれた想定範囲内に含まれている。したがって、この推定走行場e1は、想定内の走行場であると判定される。   Further, in the second embodiment, the traveling place stored in the first storage unit 30 is specified by the intermediate layer output. 25A and 25B are schematic diagrams of the assumed range. Lines II in FIGS. 25A and 25B schematically show the assumed range of the assumed field stored in the first storage unit 30. The horizontal and vertical axes in FIGS. 25A and 25B are parameters of the intermediate layer output. And in FIG. 25 (A), the estimated traveling place e1 shown with a broken line is included in the assumption range enclosed with the line II. Therefore, it is determined that the estimated traveling place e1 is an expected traveling place.

一方、図25(B)では、破線で示す推定走行場e2は、線IIで囲まれた想定範囲内に含まれていない。したがって、この推定走行場e2は、想定外の走行場であると判定される。また、図25(B)には、想定内の走行場に対する、想定が対応テーブルに挙げられている、過去の経験に基づかない走行場の範囲Tが模式的に示されている。   On the other hand, in FIG. 25B, the estimated traveling place e2 indicated by a broken line is not included in the assumed range surrounded by the line II. Therefore, it is determined that the estimated traveling place e2 is an unexpected traveling place. In addition, FIG. 25B schematically shows a range T of the driving field that is not based on the past experience, in which the assumptions for the driving field within the assumption are listed in the correspondence table.

図25(B)に示す例では、恒等写像モデルの中間層出力によって特定された想定外の推定走行場e2に最も近似する対応テーブル中の走行場として、走行場t1が選択され、一方、推定走行場e2に最も近似する想定内の走行場として、走行場bが選択される。そして、上述した第1実施形態と同様にして、想定外の走行場に対応する新たな制御アルゴリズムが創発される。   In the example shown in FIG. 25 (B), the driving field t1 is selected as the driving field in the correspondence table that most closely approximates the estimated driving field e2 that is unexpected and specified by the intermediate layer output of the identity mapping model, The driving field b is selected as the driving field within the assumption that is closest to the estimated driving field e2. In the same manner as in the first embodiment described above, a new control algorithm corresponding to an unexpected traveling place is created.

上述の実施形態においては、本発明を特定の条件で構成した例について説明したが、本発明は種々の変更及び組み合わせを行うことができ、これに限定されるものではない。   In the above-mentioned embodiment, although the example which comprised this invention on the specific conditions was demonstrated, this invention can perform a various change and combination, and is not limited to this.

10 検出部
11 ドライバ状態モニタ
12 車外環境モニタ
13 車両状態モニタ
20 推定部
30 第1記憶部
40 判定部
50 創発部
60 変換部
100 ECU
DESCRIPTION OF SYMBOLS 10 Detection part 11 Driver state monitor 12 Outside vehicle environment monitor 13 Vehicle state monitor 20 Estimation part 30 1st memory | storage part 40 Judgment part 50 Emergence part 60 Conversion part 100 ECU

Claims (3)

車両を運転するドライバの状態、車外環境状態、及び車両状態の少なくとも一つを検出する検出手段と、
上記検出手段によって検出された検出状態に基づいて、車両が遭遇している走行場を、恒等写像モデルを利用して推定する推定手段と、
既知の走行場と車両の制御アルゴリズムとを互いに対応付けて記憶した第1記憶手段と、
上記推定手段によって推定された推定走行場が、上記第1記憶手段に記憶された既知の走行場であるか、上記第1記憶手段に記憶されていない未知の走行場であるかを判定する判定手段と、
上記判定手段によって推定走行場が未知の走行場であると判定された場合に、遺伝的ネットワークプログラミングを利用して、推定走行場に対応する新たな制御アルゴリズムを創発する創発手段と、
上記判定手段によって推定走行場が既知の走行場であると判定された場合に、上記第1記憶手段から獲得した推定走行場に対応する制御アルゴリズムを車両の制御パラメータに変換し、推定走行場が未知の走行場であると判定された場合に、上記創発手段によって創発された新たな制御アルゴリズムを車両の制御パラメータに変換する変換手段と、
を有することを特徴とする車両用制御装置。
Detecting means for detecting at least one of a state of a driver driving the vehicle, an external vehicle environment state, and a vehicle state;
Based on the detection state detected by the detection means, an estimation means for estimating a running place that the vehicle encounters using an identity mapping model ;
First storage means for storing a known driving field and a vehicle control algorithm in association with each other;
Determining whether the estimated traveling place estimated by the estimating means is a known traveling place stored in the first storage means or an unknown traveling place not stored in the first storage means Means,
Emerging means for generating a new control algorithm corresponding to the estimated travel location using genetic network programming when the estimated travel location is determined to be an unknown travel location by the determination means;
When the determining means determines that the estimated traveling place is a known traveling place, the control algorithm corresponding to the estimated traveling place obtained from the first storage means is converted into a vehicle control parameter, and the estimated traveling place is A conversion means for converting a new control algorithm created by the emergence means into a vehicle control parameter when it is determined that the vehicle is an unknown traveling place;
A vehicle control device characterized by comprising:
上記第1記憶手段に格納された、既知の走行場に対する制御アルゴリズムは、遺伝的ネットワークプログラミングを利用して創発したものである
ことを特徴とする請求項記載の車両用制御装置。
The stored in the first storage means, the control algorithm for the known travel field, the vehicle control device according to claim 1, wherein a is obtained by emergent using genetic network programming.
更に、前記創発手段によって創発された新たな制御アルゴリズムを記憶する第3記憶手段を有することを特徴とする請求項1又は2記載の車両用制御装置。 Further, the vehicle control device according to claim 1, wherein further comprising a third storage means for storing a new control algorithms emergent by the emergent means.
JP2010164216A 2010-07-21 2010-07-21 Vehicle control device Expired - Fee Related JP5665400B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2010164216A JP5665400B2 (en) 2010-07-21 2010-07-21 Vehicle control device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2010164216A JP5665400B2 (en) 2010-07-21 2010-07-21 Vehicle control device

Publications (2)

Publication Number Publication Date
JP2012025232A JP2012025232A (en) 2012-02-09
JP5665400B2 true JP5665400B2 (en) 2015-02-04

Family

ID=45778717

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2010164216A Expired - Fee Related JP5665400B2 (en) 2010-07-21 2010-07-21 Vehicle control device

Country Status (1)

Country Link
JP (1) JP5665400B2 (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6011788B2 (en) * 2012-09-03 2016-10-19 マツダ株式会社 Vehicle control device
KR102028064B1 (en) * 2013-03-29 2019-10-02 르노 에스.아.에스. Method and device for controlling an energy equivalence factor in a hybrid motor propulsion plant
CN105599772B (en) * 2014-11-24 2018-06-22 上海通用汽车有限公司 Vehicular intelligent driving assistance system and intelligent driving auxiliary control method
KR101827124B1 (en) * 2016-06-13 2018-02-07 현대자동차주식회사 System and Method for recognizing driving pattern of driver
JP7382304B2 (en) * 2020-12-17 2023-11-16 株式会社日立製作所 Risk management device, risk management method and risk management system

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000066723A (en) * 1998-08-25 2000-03-03 Kobe Steel Ltd Route constitution design device for mobile object
JP2002264687A (en) * 2001-03-07 2002-09-18 Nissan Motor Co Ltd Vehicle-to-vehicle distance control device
JP4175923B2 (en) * 2003-03-13 2008-11-05 株式会社エクォス・リサーチ Traveling speed pattern estimation device
JP4728839B2 (en) * 2006-03-02 2011-07-20 株式会社デンソーアイティーラボラトリ In-vehicle device controller
JP2008018872A (en) * 2006-07-13 2008-01-31 Nissan Motor Co Ltd Driving support device, automobile and driving support method
JP2008179195A (en) * 2007-01-23 2008-08-07 Toyota Motor Corp Vehicle control device and arithmetic storage device

Also Published As

Publication number Publication date
JP2012025232A (en) 2012-02-09

Similar Documents

Publication Publication Date Title
Huang et al. Personalized trajectory planning and control of lane-change maneuvers for autonomous driving
CN114312830B (en) Intelligent vehicle coupling decision model and method considering dangerous driving conditions
JP5665400B2 (en) Vehicle control device
US20210086798A1 (en) Model-free reinforcement learning
JP6011788B2 (en) Vehicle control device
Wu et al. Deep reinforcement learning-based energy-efficient decision-making for autonomous electric vehicle in dynamic traffic environments
JP2020059497A (en) Method and device of controlling automatic driving introducing &#34;detection acquisition&#34; model, and program constituted thereby
Selvaraj et al. An ML-aided reinforcement learning approach for challenging vehicle maneuvers
CN115257819A (en) Decision-making method for safe driving of large-scale commercial vehicle in urban low-speed environment
CN113722835A (en) Modeling method for anthropomorphic random lane change driving behavior
CN112078371A (en) Energy recovery method and energy recovery device of composite power electric vehicle
JP7350188B2 (en) Driving support device, learning device, driving support method, driving support program, learned model generation method, learned model generation program
Li et al. Energy-efficient autonomous vehicle control using reinforcement learning and interactive traffic simulations
Ziadia et al. An adaptive regenerative braking strategy design based on naturalistic regeneration performance for intelligent vehicles
Feiyan et al. A review of machine learning on energy management strategy for hybrid electric vehicles
Tang et al. Personalized decision-making and control for automated vehicles based on generative adversarial imitation learning
van der Ploeg et al. Connecting the dots: Context-driven motion planning using symbolic reasoning
Chevrant-Breton et al. Search space reduction in dynamic programming using monotonic heuristics in the context of model predictive optimization
Zhang et al. Game Theory-Based Lane Change Decision-Making Considering Vehicle’s Social Value Orientation
Smith et al. Fuzzy rule-based energy management strategy for a parallel mild-hybrid electric bus
Chen et al. Decision making for overtaking of unmanned vehicle based on deep Q-learning
CN117227752A (en) Vehicle autopilot control using planning models and cognitive learning models
Kwon et al. Driving Profile Optimization Using a Deep Q‐Network to Enhance Electric Vehicle Battery Life
Xu et al. Estimating risk levels of driving scenarios through analysis of driving styles for autonomous vehicles
Duan et al. Novel adaptive path-smoothening optimization method for mobile robots

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20130605

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130605

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20140115

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20140116

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140313

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20140402

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140530

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: 20141112

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20141209

R150 Certificate of patent or registration of utility model

Ref document number: 5665400

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees