JPH11324782A - Learning method in air-fuel ratio control device of internal combustion engine - Google Patents

Learning method in air-fuel ratio control device of internal combustion engine

Info

Publication number
JPH11324782A
JPH11324782A JP10132197A JP13219798A JPH11324782A JP H11324782 A JPH11324782 A JP H11324782A JP 10132197 A JP10132197 A JP 10132197A JP 13219798 A JP13219798 A JP 13219798A JP H11324782 A JPH11324782 A JP H11324782A
Authority
JP
Japan
Prior art keywords
fuel
fuel ratio
air
teacher data
operating 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.)
Pending
Application number
JP10132197A
Other languages
Japanese (ja)
Inventor
Yuichiro Sawada
雄一郎 沢田
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.)
Yamaha Motor Co Ltd
Original Assignee
Yamaha 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 Yamaha Motor Co Ltd filed Critical Yamaha Motor Co Ltd
Priority to JP10132197A priority Critical patent/JPH11324782A/en
Publication of JPH11324782A publication Critical patent/JPH11324782A/en
Pending legal-status Critical Current

Links

Landscapes

  • Electrical Control Of Air Or Fuel Supplied To Internal-Combustion Engine (AREA)
  • Combined Controls Of Internal Combustion Engines (AREA)
  • Feedback Control In General (AREA)

Abstract

PROBLEM TO BE SOLVED: To obtain a learning method cupuble of simply and accurate by initializing by acquiring the middle value of plural teaching data as teaching data of its corresponding operating condition points. SOLUTION: The whole area of operating condition parameters forming a model is divided into plural divided operation areas, and the operating condition is changed in transient between two operating conditions a1, a2 dimarcating the respective deviled operation areas. From the amount of shift of an exhaust air fuel ratio from a target air fuel ratio caused at this time, teaching data for correcting output parameters of the model in two operating conditions a1, a2 changed to be transient is acquired. In the case where two teaching data are acquired for one operating condition, the average value is taken as teaching data to perform initial learning of the model. Thus, it is possible to obtain a model such that output parameter is smoothly changed in all areas of operating condition parameters without complicated operation.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、少なくとも燃料噴
射装置から噴射された燃料の動的挙動を燃料付着率とし
てモデル化し、該モデルを用いて空燃比の制御を行う内
燃機関の空燃比制御装置における前記モデルの学習方法
に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an air-fuel ratio control apparatus for an internal combustion engine which models at least the dynamic behavior of fuel injected from a fuel injection device as a fuel adhesion rate and controls the air-fuel ratio using the model. In the model learning method.

【0002】[0002]

【従来の技術】近年、排ガス規制の強化や、市場からの
燃費向上及び応答性の向上の要望等に対応するために内
燃機関に対する高精度な燃料供給が求められており、こ
れを達成するために、燃料噴射装置を用いて内燃機関に
燃料を供給することが主流になってきている。しかし、
燃料噴射装置は従来の気化器に比べて吸気管及び吸気弁
に付着する燃料の量が多く、過渡時の空燃比制御性が悪
化してしまうという問題があった。このような問題を解
決するために、出願人は、特願平7−271188号
で、燃料噴射装置から噴射された燃料が吸気管及び吸気
弁に付着する割合を燃料付着率とし、吸気管及び吸気弁
に付着した燃料が蒸発する速度を蒸発時定数として各々
学習機能付きフィードフォワード制御ロジックを用いて
モデル化し、この燃料系モデルを用いて燃焼室に入る燃
料量を推定し、この推定燃料量に基づいて噴射すべき燃
料量を算出する方法を既に提案している。
2. Description of the Related Art In recent years, there has been a demand for high-precision fuel supply to internal combustion engines in order to respond to stricter exhaust gas regulations and to market demands for improved fuel efficiency and responsiveness. In addition, it has become mainstream to supply fuel to an internal combustion engine using a fuel injection device. But,
The fuel injection device has a problem that the amount of fuel adhering to the intake pipe and the intake valve is larger than that of the conventional carburetor, and the controllability of the air-fuel ratio during the transition is deteriorated. In order to solve such a problem, the applicant has disclosed in Japanese Patent Application No. 7-271188 a fuel attachment ratio in which the rate at which fuel injected from a fuel injector adheres to an intake pipe and an intake valve is defined as a fuel adhesion rate. The rate at which fuel attached to the intake valve evaporates is modeled using the feedforward control logic with a learning function as the evaporation time constant, and the amount of fuel entering the combustion chamber is estimated using this fuel system model. A method for calculating the amount of fuel to be injected based on the above has already been proposed.

【0003】[0003]

【発明が解決しようとする課題】上記した制御方法によ
れば、燃料系モデルを用いて燃料付着率と蒸発時定数と
を推定でき、過渡時に燃焼室内に入る燃料量を予め推定
できるので、燃料噴射量を適当な量に制御することがで
きるが、前記燃料系モデルを用いる場合には、予め様々
な運転状態に対する燃料付着率及び蒸発時定数の最適値
を設定する初期学習を行わせる必要があった。燃料の動
的挙動はエンジンの運転状態が過渡状態にある時に空燃
比に影響を及ぼすものなので、前記した動的挙動のパラ
メータの設定も運転状態を過渡的に変化させながら行わ
なければならないが、動的挙動のパラメータの一つであ
る燃料付着率は、運転状態を過渡的に変化させた時にそ
の変化の前後の運転状態に影響を及ぼすものなので、そ
の設定が難しく、手間がかかるという問題があった。本
発明は、上記した燃料系モデルを用いた空燃比制御方法
における初期学習に関する問題点に鑑みて、簡単により
正確な初期設定を行うことができる内燃機関の空燃比制
御装置における学習方法を提供することを目的としてい
る。
According to the control method described above, the fuel adhesion rate and the evaporation time constant can be estimated using the fuel system model, and the amount of fuel entering the combustion chamber during transition can be estimated in advance. Although the injection amount can be controlled to an appropriate amount, in the case of using the fuel system model, it is necessary to perform an initial learning for setting an optimum value of a fuel adhesion rate and an evaporation time constant for various operation states in advance. there were. Since the dynamic behavior of the fuel affects the air-fuel ratio when the operating state of the engine is in the transient state, the parameters of the dynamic behavior must be set while the operating state is transiently changed. The fuel deposition rate, one of the parameters of the dynamic behavior, affects the operating state before and after the change when the operating state is transiently changed, so the setting is difficult and it takes time. there were. The present invention provides a learning method in an air-fuel ratio control apparatus for an internal combustion engine, which can easily and accurately perform an initial setting in view of the above-mentioned problems relating to an initial learning in an air-fuel ratio control method using a fuel system model. It is intended to be.

【0004】[0004]

【課題を解決するための手段】上記した目的を達成する
ために、本発明に係る内燃機関の空燃比制御装置におけ
る学習方法は、燃料噴射装置から噴射された燃料の動的
挙動を少なくとも燃料付着率としてモデル化し、該モデ
ルを用いて排気空燃比を制御する内燃機関の空燃比制御
装置において、前記運転状態に関するパラメータを複数
領域に割り付けて、一つの領域を画定する運転状態点間
で運転条件を過渡的に変化させ、その時に生じる排気空
燃比の目標空燃比に対するずれに基づいて、変化させた
二つの運転状態点の燃料付着率の教師データを獲得し、
前記教師データに基づいてモデルの学習を行うと共に、
同一の運転状態点に対して異なる領域に関する複数の異
なる教師データを獲得した場合には、それら複数の教師
データの中間値を、その対応する運転状態点の教師デー
タとして獲得することを特徴とするものである。
In order to achieve the above-mentioned object, a learning method in an air-fuel ratio control apparatus for an internal combustion engine according to the present invention comprises a method for learning at least the dynamic behavior of fuel injected from a fuel injection device. In an air-fuel ratio control device for an internal combustion engine, which models the ratio as an air condition and controls the exhaust air-fuel ratio using the model, the operating conditions are assigned to a plurality of regions, and operating conditions are defined between operating state points defining one region. Is transiently changed, and, based on the deviation of the exhaust air-fuel ratio from the target air-fuel ratio that occurs at that time, the teacher data of the fuel deposition rate at the two changed operating state points is obtained,
While learning the model based on the teacher data,
When a plurality of different teacher data for different regions are obtained for the same driving state point, an intermediate value of the plurality of teacher data is obtained as the teacher data of the corresponding driving state point. Things.

【0005】[0005]

【発明の実施の形態】以下、添付図面に示した幾つかの
実施例を参照して、本発明に係る内燃機関の空燃比制御
装置における学習方法の実施の形態について説明する。
図1は、エンジン1と本発明に係る予測学習方法を採用
したエンジン制御を実行可能な制御装置10との関係を
示す概略図である。エンジン1は、吸気管2に設けられ
たエアクリーナ3及び燃料噴射装置4を介してシリンダ
5の燃焼室内に混合気を導入し、燃焼後の排気ガスを排
気管6を介して大気中に排気する4サイクルエンジンで
あり、本図では吸気バルブや排気バルブ等の他の構成部
材については省略されている。尚、図1中、符号7はク
ランクケースを、また、符号8はスロットルバルブを各
々示している。制御装置10は、燃料噴射装置4からの
燃料噴射量を操作して排気ガス中の空燃比の値を制御す
るものである。この制御装置10は、図1に示すよう
に、クランクケース7に設けられたクランク角検知手段
13から得られるクランク角に関する情報rと、スロッ
トルバルブ8に設けられたスロットル開度検知手段12
から得られるスロットル開度に関する情報αと、吸気管
2に設けられた吸気管壁温検知手段14から得られる吸
気管壁温に関する情報Teとを入力し、これらの入力情
報に基づいて、吸気管2に設けられた燃料噴射装置4の
操作量Mf(即ち、燃料噴射量)を決定して出力すると
共に、排気管6に設けられた空燃比センサ15から得ら
れる実際の空燃比に関する検知信号Eを、必要に応じて
入力して、この情報に基づく内部モデルの出力の補正及
び内部モデルの学習を行い常時最適な制御が行えるよう
に構成されている。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS An embodiment of a learning method in an air-fuel ratio control apparatus for an internal combustion engine according to the present invention will be described below with reference to some embodiments shown in the accompanying drawings.
FIG. 1 is a schematic diagram showing the relationship between an engine 1 and a control device 10 capable of executing engine control employing a prediction learning method according to the present invention. In the engine 1, an air-fuel mixture is introduced into a combustion chamber of a cylinder 5 via an air cleaner 3 and a fuel injection device 4 provided in an intake pipe 2, and exhaust gas after combustion is exhausted to the atmosphere via an exhaust pipe 6. This is a four-cycle engine, and other components such as an intake valve and an exhaust valve are omitted in the drawing. In FIG. 1, reference numeral 7 denotes a crankcase, and reference numeral 8 denotes a throttle valve. The control device 10 controls the value of the air-fuel ratio in the exhaust gas by operating the amount of fuel injection from the fuel injection device 4. As shown in FIG. 1, the control device 10 includes information r regarding a crank angle obtained from a crank angle detecting means 13 provided in a crankcase 7 and a throttle opening detecting means 12 provided in a throttle valve 8.
, And information Te about the intake pipe wall temperature obtained from the intake pipe wall temperature detecting means 14 provided in the intake pipe 2, and based on the input information, the intake pipe The operation amount Mf (that is, the fuel injection amount) of the fuel injection device 4 provided in the exhaust pipe 6 is determined and output, and the detection signal E relating to the actual air-fuel ratio obtained from the air-fuel ratio sensor 15 provided in the exhaust pipe 6 is obtained. Is input as necessary, the output of the internal model is corrected based on this information, and the learning of the internal model is performed, so that optimal control can always be performed.

【0006】図2は、制御装置10の構成を示す概略ブ
ロック図である。制御装置10は、目標空燃比Epに沿
ってフィードフォワード制御ロジックを用いて基本燃料
噴射量Mfnを決定するモデルベース制御部20、前記
クランク角情報rからエンジン回転数Nを算出するエン
ジン回転数演算部30、エンジン回転数Nとスロットル
開度とに基づいて予め用意したマップ又は演算式に基づ
いて前記目標空燃比Epを算出する目標空燃比算出部4
0、前記モデルベース制御部20で得られる基本燃料噴
射量Mfnをパルス信号としての燃料噴射量Mfに変換
する変換部50、及び前記モデルベース制御部20のず
れ量を検出し、そのずれを補正するための教師データの
獲得を行うモデル補正部60を備えている。
FIG. 2 is a schematic block diagram showing the configuration of the control device 10. The control device 10 includes a model-based control unit 20 that determines a basic fuel injection amount Mfn using a feedforward control logic along the target air-fuel ratio Ep, and an engine speed calculation that calculates an engine speed N from the crank angle information r. Unit 30, a target air-fuel ratio calculation unit 4 that calculates the target air-fuel ratio Ep based on a map or an arithmetic expression prepared in advance based on the engine speed N and the throttle opening degree
0, a conversion unit 50 that converts the basic fuel injection amount Mfn obtained by the model base control unit 20 into a fuel injection amount Mf as a pulse signal, and detects a deviation amount of the model base control unit 20 and corrects the deviation. And a model correction unit 60 for acquiring teacher data for performing the operation.

【0007】(モデルベース制御部について)以下、図
3〜図10を参照してモデルベース制御部20の構成に
ついて説明する。図3は、図2におけるモデルベース制
御部20の構成を示す概略ブロック図である。モデルベ
ース制御部20は、吸気管2内の空気の挙動をモデル化
した空気系順モデル21、燃料噴射装置4から噴射され
る燃料の挙動をモデル化した燃料系順モデル22、及び
各順モデル21及び22から出力される推定空気量Av
及び推定燃料量Fvに基づいて推定空燃比Evを算出す
る推定空燃比演算部23を備えている。また、モデルベ
ース制御部20は、推定空燃比演算部23から出力され
る推定空燃比Evを基本操作量演算部24にフィードバ
ックするフィードバックループを備えている。前記基本
操作量演算部24は、推定空燃比演算部23から出力さ
れる推定空燃比Evと、目標空燃比算出部30から出力
される目標空燃比Epとを入力してエンジン1の燃料噴
射装置4に対する基本操作量Mfn(基本燃料噴射量)
を算出する。この基本操作量Mfnは、モデルベース制
御部20から出力されると共に燃料系順モデル22にも
入力され、燃料系順モデル22は前記基本操作量Mfn
に基づいて推定燃料量Fvを求める。上記したように、
モデルベース制御部20では、空気系順モデル21、燃
料系順モデル22、及び推定空燃比演算部23によりエ
ンジン1の順モデルを構成し、かつ、前記燃料系順モデ
ル22、推定空燃比演算部23、及び基本操作量演算部
24を含むフィードバックループを用いて前記エンジン
1の順モデルから出力される推定空燃比Evをフィード
バックして基本操作量Mfnを出力するエンジンの逆モ
デルを構成している。
(Model Base Control Unit) The configuration of the model base control unit 20 will be described below with reference to FIGS. FIG. 3 is a schematic block diagram showing the configuration of the model base control unit 20 in FIG. The model-based control unit 20 includes an air system forward model 21 that models the behavior of air in the intake pipe 2, a fuel system forward model 22 that models the behavior of fuel injected from the fuel injection device 4, and each forward model. Estimated air amount Av output from 21 and 22
And an estimated air-fuel ratio calculation unit 23 for calculating an estimated air-fuel ratio Ev based on the estimated fuel amount Fv. Further, the model base control unit 20 includes a feedback loop for feeding back the estimated air-fuel ratio Ev output from the estimated air-fuel ratio calculation unit 23 to the basic operation amount calculation unit 24. The basic manipulated variable calculator 24 receives the estimated air-fuel ratio Ev output from the estimated air-fuel ratio calculator 23 and the target air-fuel ratio Ep output from the target air-fuel ratio calculator 30 and receives the fuel injection device of the engine 1. 4 basic operation amount Mfn (basic fuel injection amount)
Is calculated. The basic operation amount Mfn is output from the model base control unit 20 and also input to the fuel system forward model 22. The fuel system forward model 22 outputs the basic operation amount Mfn.
The estimated fuel amount Fv is obtained based on As mentioned above,
In the model base control unit 20, a forward model of the engine 1 is configured by the air system forward model 21, the fuel system forward model 22, and the estimated air-fuel ratio calculation unit 23, and the fuel system forward model 22, the estimated air-fuel ratio calculation unit 23 and a feedback loop including a basic manipulated variable operation unit 24 to constitute an inverse model of the engine that outputs the basic manipulated variable Mfn by feeding back the estimated air-fuel ratio Ev output from the forward model of the engine 1. .

【0008】(空気系順モデルについて)前記空気系順
モデル21は、スロットル開度α及びエンジン回転数N
を入力し、これらに基づいて吸入空気量を推定算出す
る。具体的には、スロットル開度及び吸気負圧を用いた
流体力学的な数式で空気量をモデル化し、さらに空気量
から算出される圧力と体積効率とに基づいて流体力学的
な数式で吸気負圧をモデル化すると共に、エンジン回転
数N及びスロットル開度αを入力とするファジィニュー
ラル回路網を用いて体積効率をモデル化し、前記吸気負
圧の初期値を設定することで、スロットル開度α及びエ
ンジン回転数Nから推定空気量Avが算出できるように
構成され得る。尚、前記体積効率のモデルを構成するフ
ァジィニューラル回路網は、エンジンが定常状態にある
時に、排気空燃比をフィードバックし、この排気空燃比
と目標空燃比とのずれを補正するよう学習可能に構成さ
れ得る。また、前記体積効率のモデル化は、ファジィニ
ューラル回路網に限定されるものではなく、ニューラル
回路網、CMAC又は単なるマップ等、学習可能なフィ
ードフォワード制御ロジックであれば任意のものが使用
され得る。
The air system forward model 21 has a throttle opening α and an engine speed N.
Is input, and the intake air amount is estimated and calculated based on these. Specifically, the air amount is modeled by a hydrodynamic formula using the throttle opening and the intake negative pressure, and the intake negative pressure is modeled by a hydrodynamic formula based on the pressure and the volumetric efficiency calculated from the air amount. The pressure is modeled, the volumetric efficiency is modeled using a fuzzy neural network that takes the engine speed N and the throttle opening α as inputs, and the initial value of the intake negative pressure is set. And the estimated air amount Av can be calculated from the engine speed N. The fuzzy neural network constituting the volumetric efficiency model is configured to be able to learn to feed back the exhaust air-fuel ratio and correct the deviation between the exhaust air-fuel ratio and the target air-fuel ratio when the engine is in a steady state. Can be done. Further, the modeling of the volumetric efficiency is not limited to a fuzzy neural network, and any feedforward control logic that can be learned, such as a neural network, a CMAC, or a simple map, may be used.

【0009】(燃料系順モデルについて)図4は前記燃
料系順モデル22の構成を示す概略ブロック図であり、
図5は燃料系順モデルのイメージ図を示している。この
燃料系順モデル22は、図5に示すように燃料噴射装置
4から噴射された燃料の挙動を、噴射された燃料が直接
シリンダ内に入らないで吸気管壁や吸気弁等に付着する
割合x(以下、燃料付着率xと称する)と、吸気管壁や
吸気弁等に付着した燃料が蒸発する速度τ(以下、蒸発
時定数τ)とに分けてモデル化し、前記基本操作量演算
部24から入力される基本操作量Mfn(基本燃料噴射
量)と前記燃料付着率とに基づいて燃料噴射装置4から
噴射され直接シリンダ内に入る燃料量(1−x)Mfを
推定算出すると共に、基本操作量Mfn及び燃料付着率
xから得られる付着燃料Mwと蒸発時定数τとに基づい
て付着燃料Mwが蒸発してシリンダ内に入る量a・Mw
を推定算出し、これら推定算出した燃料量(1−x)M
fとa・Mwとからシリンダ内に入る総燃料量を推定す
るもので、燃料付着率推定部22a、蒸発時定数推定部
22b、非付着燃料演算部22c、付着燃料演算部22
d、及び蒸発燃料量演算部22eを備える。前記燃料付
着率推定部22aは、図6に示すように、エンジンの運
転状態(エンジン回転数N及びスロットル開度α)に対
する燃料付着率xをモデル化したマップから成り、スロ
ットル開度検知手段12及びエンジン回転数演算部30
から、その時々のエンジン回転数N及びスロットル開度
を入力し、これらに基づいて燃料付着率xを決定する。
また、前記蒸発時定数推定部22bは、エンジンの運転
状態(エンジン回転数N若しくはスロットル開度α及び
吸気管壁温Te若しくはエンジン水温)に対する蒸発時
定数τをモデル化したマップから成り、スロットル開度
検知手段12若しくはエンジン回転数演算部30及び吸
気管壁温検知手段14から、その時々のエンジン回転数
N若しくはスロットル開度及び吸気管壁温Teを入力
し、これらに基づいて燃料付着率xを決定する。非付着
燃料演算部22cは、前記燃料付着率推定部22aから
得られる燃料付着率xに基づいて、基本操作量演算部2
4から入力される基本操作量Mfn(即ち、基本燃料噴
射量)における燃料噴射装置4から直接シリンダ5の燃
焼室内に入る燃料量(1−x)Mfを算出する。付着燃
料演算部22dは、前記燃料付着率推定部22aから得
られる燃料付着率xに基づいて、基本操作量演算部24
から入力される基本操作量Mfn(基本燃料噴射量)に
おいて吸気管壁面及び吸気弁に付着する燃料量x・Mf
(Mw)を算出する。前記蒸発燃料量演算部22eは、
前記蒸発時定数推定部22bで推定された蒸発時定数τ
を用いて付着燃料演算部22bで得られる付着燃料から
蒸発する燃料を積算して蒸発燃料量a・Mwを算出す
る。燃料系順モデル22は、前記した非付着燃料演算部
22cで算出された燃料量と蒸発燃料量演算部22eで
算出された蒸発燃料量とを加算して、推定燃料量Fvを
算出する。尚、上記した燃料付着率推定部22a及び蒸
発時定数推定部22bを構成する各マップには、内燃機
関の様々な運転状態に対応できるように、予め、様々な
運転状態に対して最適な燃料付着率x又は蒸発時定数τ
を基礎データとして獲得する初期学習が施され、初期学
習後の基礎データで構成されたマップを用いて実際に内
燃機関は駆動される。しかし、予め初期学習によりマッ
プ(モデル)の最適化を行っておいても、内燃機関駆動
後に、内燃機関の経時変化や内燃機関を使用する環境の
変化等により基礎データと実際の内燃機関における動的
挙動とのずれが生じる。このずれは、後述する方法で獲
得された教師データに基づいて、モデルと実際の内燃機
関における燃料の動的挙動とのずれを補正するオンライ
ン学習により補正される。これらの初期学習のための教
師データ及びオンライン学習のための教師データの獲得
は後で詳述するモデル補正部60により行われる。
FIG. 4 is a schematic block diagram showing the structure of the fuel system forward model 22.
FIG. 5 shows an image diagram of the fuel system forward model. As shown in FIG. 5, the fuel system forward model 22 determines the behavior of the fuel injected from the fuel injection device 4 by the ratio of the injected fuel adhering to the intake pipe wall or the intake valve without directly entering the cylinder. x (hereinafter referred to as a fuel adhesion rate x) and a speed τ (hereinafter, an evaporation time constant τ) at which fuel adhering to an intake pipe wall or an intake valve evaporates and is modeled. Based on the basic operation amount Mfn (basic fuel injection amount) input from 24 and the fuel adhesion rate, a fuel amount (1-x) Mf injected from the fuel injection device 4 and directly into the cylinder is estimated and calculated, Amount a · Mw in which the deposited fuel Mw evaporates and enters the cylinder based on the deposited fuel Mw obtained from the basic operation amount Mfn and the fuel deposition rate x and the evaporation time constant τ.
And the estimated fuel amount (1-x) M
It estimates the total amount of fuel entering the cylinder from f and a · Mw, and includes a fuel adhesion rate estimating unit 22a, an evaporation time constant estimating unit 22b, a non-adhering fuel calculating unit 22c, and an adhering fuel calculating unit 22.
d and an evaporative fuel amount calculation unit 22e. As shown in FIG. 6, the fuel adhesion rate estimating unit 22a comprises a map in which the fuel adhesion rate x is modeled with respect to the operating state of the engine (engine speed N and throttle opening α). And engine speed calculation unit 30
Then, the current engine speed N and throttle opening are input, and the fuel adhesion rate x is determined based on these.
The evaporating time constant estimating unit 22b is composed of a map that models the evaporating time constant τ with respect to the operating state of the engine (the engine speed N or the throttle opening α and the intake pipe wall temperature Te or the engine water temperature). The engine speed N or the throttle opening and the intake pipe wall temperature Te at that time are input from the degree detecting means 12 or the engine speed calculating section 30 and the intake pipe wall temperature detecting means 14, and based on these, the fuel adhesion rate x To determine. The non-adhered fuel calculating unit 22c is configured to calculate the basic operation amount calculating unit 2 based on the fuel adhering rate x obtained from the fuel adhering rate estimating unit 22a.
The fuel amount (1-x) Mf that directly enters the combustion chamber of the cylinder 5 from the fuel injection device 4 at the basic operation amount Mfn (that is, the basic fuel injection amount) input from Step 4 is calculated. The adhering fuel calculating unit 22d is configured to calculate the basic operation amount calculating unit 24 based on the fuel adhering rate x obtained from the fuel adhering rate estimating unit 22a.
Amount x · Mf adhering to the intake pipe wall surface and the intake valve in the basic operation amount Mfn (basic fuel injection amount) input from
(Mw) is calculated. The evaporative fuel amount calculation unit 22e includes:
The evaporation time constant τ estimated by the evaporation time constant estimation unit 22b
Is used to calculate the evaporated fuel amount a · Mw by integrating the fuel evaporated from the attached fuel obtained by the attached fuel calculation unit 22b. The fuel system forward model 22 calculates the estimated fuel amount Fv by adding the fuel amount calculated by the non-adhered fuel calculating unit 22c and the evaporated fuel amount calculated by the evaporated fuel amount calculating unit 22e. The maps constituting the above-described fuel adhesion rate estimating section 22a and the evaporation time constant estimating section 22b include, in advance, optimal fuels for various operating states so as to correspond to various operating states of the internal combustion engine. Adhesion rate x or evaporation time constant τ
Is acquired as basic data, and the internal combustion engine is actually driven using a map composed of the basic data after the initial learning. However, even if the map (model) is optimized in advance by the initial learning, the basic data and the dynamics of the actual internal combustion engine may be changed after the internal combustion engine is driven due to a change over time of the internal combustion engine or a change in the environment in which the internal combustion engine is used. Deviation from the dynamic behavior occurs. This deviation is corrected by online learning that corrects the deviation between the model and the actual dynamic behavior of the fuel in the internal combustion engine based on the teacher data acquired by a method described later. The acquisition of the teacher data for the initial learning and the teacher data for the online learning is performed by the model correction unit 60 described later in detail.

【0010】推定空燃比演算部23は、前記空気系順モ
デル21及び燃料系順モデル22から得られる推定空気
量Av及び推定燃料量Fvに基づいて推定空燃比Evを
算出する。そして、基本操作量演算部24は、前記推定
空燃比Evと目標空燃比Epとに基づいて基本操作量
(即ち、基本燃料噴射量Mfn)を算出し、変換部50
に出力する。
The estimated air-fuel ratio calculation unit 23 calculates an estimated air-fuel ratio Ev based on the estimated air amount Av and estimated fuel amount Fv obtained from the air system forward model 21 and the fuel system forward model 22. Then, the basic operation amount calculation unit 24 calculates a basic operation amount (that is, a basic fuel injection amount Mfn) based on the estimated air-fuel ratio Ev and the target air-fuel ratio Ep, and converts the basic operation amount into a conversion unit 50.
Output to

【0011】(モデル補正部について)次に、上記した
ように構成された制御装置における燃料系順モデル22
の初期学習及びオンライン学習について説明していく。
前記燃料系順モデル22における燃料付着率推定部22
a及び蒸発時定数推定部22bの各モデルの初期学習又
はオンライン学習用の教師データはモデル補正部60で
獲得され、各推定部22a及び22bは、この教師デー
タに基づいて、その時の運転状態に対応する燃料付着率
x又は蒸発時定数τの初期学習又はオンライン学習を行
う。
(Regarding the Model Correction Unit) Next, the fuel system forward model 22 in the control device configured as described above
The following describes the initial learning and online learning.
Fuel adhesion rate estimation unit 22 in the fuel system forward model 22
The teaching data for the initial learning or online learning of each model of the e and the evaporation time constant estimating unit 22b is obtained by the model correcting unit 60, and the estimating units 22a and 22b determine the operating state at that time based on the teaching data. Initial learning or online learning of the corresponding fuel adhesion rate x or evaporation time constant τ is performed.

【0012】(初期学習について)始めに燃料付着率推
定部22aを例に挙げて初期学習について説明する。図
7は、燃料付着率推定部22aを構成するマップの初期
状態を示している。燃料付着率推定部22a(又は蒸発
時定数推定部22b)を構成するマップは、エンジンの
運転状態を表す各パラメータ(即ち、図7ではエンジン
回転数及びスロットル開度)の全域(図7では、エンジ
ン回転数N:0rpm〜6000rpm;スロットル開
度α:0度〜90度)を適当な運転状態で複数に分割し
(図7ではエンジン回転数Nは600rpm単位で10
個に分割され、スロットル開度は10度単位で9個に分
割されている。)、分割された各運転領域(以下、分割
運転領域と称する。)を画定する各運転状態毎に最適な
出力パラメータ、即ち、燃料付着率xを設定するよう構
成されたデータの集合である。本明細書では、分割運転
領域を画定する各運転状態を運転状態点と称することと
し、従って、例えば、図7における分割運転領域Zは、
エンジン回転数3600rpmでスロットル開度が50
度という運転状態点z1と、エンジン回転数3600r
pmでスロットル開度が60度という運転状態点z2
と、エンジン回転数4200rpmでスロットル開度が
50度という運転状態点z3と、エンジン回転数420
0rpmでスロットル開度が60度という運転状態点z
4とで画定されていることになる。前記したマップにお
ける燃料付着率xは、初期状態では、全ての運転状態点
で適当な一定の値(図7では20%)に設定されてい
る。初期学習では、この状態から、運転状態の一方のパ
ラメータ(エンジン回転数)を固定し、他方のパラメー
タ(スロットル開度)を二つの隣接する運転状態点間で
過渡的に変化させ、その時に生じる排気空燃比の波形を
空燃比センサ15で検出し、検出された排気空燃比の波
形に基づいて、過渡的に変化させた二つの運転状態点に
おける出力パラメータの教師データを獲得し、この教師
データに基づいて出力パラメータの学習を行う。尚、図
示していないが、蒸発時定数推定部22bも、燃料付着
率推定部22aと同様に、エンジンの運転状態に対する
蒸発時定数のマップを備えており、運転状態全域が燃料
付着率22aと対応するように複数に分割されている。
(Initial Learning) Initial learning will be described with reference to the fuel adhesion rate estimating unit 22a as an example. FIG. 7 shows an initial state of a map constituting the fuel adhesion rate estimating unit 22a. The map constituting the fuel adhesion rate estimating unit 22a (or the evaporation time constant estimating unit 22b) includes an entire region (in FIG. 7, the engine speed and the throttle opening in FIG. 7) representing the operating state of the engine. The engine speed N: 0 rpm to 6000 rpm; the throttle opening α: 0 to 90 degrees) is divided into a plurality of parts in an appropriate operation state (in FIG. 7, the engine speed N is 10 in units of 600 rpm).
And the throttle opening is divided into nine in 10-degree units. ), A set of data configured to set an optimal output parameter, that is, a fuel adhesion rate x for each operation state that defines each divided operation region (hereinafter, referred to as a divided operation region). In the present specification, each operation state that defines the divided operation region is referred to as an operation state point. Therefore, for example, the divided operation region Z in FIG.
Throttle opening is 50 at 3600 rpm
Operating state point z1 in degrees and engine speed 3600r
Operating state point z2 with 60 degree throttle opening at pm
An operating state point z3 at which the engine rotation speed is 4200 rpm and the throttle opening is 50 degrees;
Operating state point z where throttle opening is 60 degrees at 0 rpm
4 will be defined. In the initial state, the fuel adhesion rate x in the above-described map is set to an appropriate constant value (20% in FIG. 7) at all operating state points. In the initial learning, from this state, one parameter (engine speed) of the operating state is fixed, and the other parameter (throttle opening) is transiently changed between two adjacent operating state points. The waveform of the exhaust air-fuel ratio is detected by the air-fuel ratio sensor 15, and based on the detected waveform of the exhaust air-fuel ratio, teacher data of output parameters at two transiently changed operating state points is obtained. Learning of output parameters is performed based on. Although not shown, the evaporating time constant estimating unit 22b also has a map of the evaporating time constant with respect to the operating state of the engine, similarly to the fuel adhering rate estimating unit 22a. It is divided into a plurality to correspond.

【0013】以下に、具体的な例を挙げて教師データの
獲得及び学習からなる初期学習について説明していく。
始めに分割運転領域aを画定する二つの運転状態点a1
及びa2における教師データの獲得を例に挙げて説明す
ると、エンジン回転数演算部30で算出されるエンジン
回転数1000rpmになるように初期学習用のクラン
ク角度信号を制御装置10に入力し、制御装置10にお
けるエンジン回転数Nを運転状態点a1及びa2を含む
1000rpmに固定する。そして、運転状態点a1に
関する初期学習用のスロットル開度信号(即ち、0度信
号)を入力し、前記スロットル開度信号を過渡的に運転
状態点a2(即ち、10度信号)まで変化させる。図7
に示す初期状態では、マップの出力パラメータ値は適当
な初期値(本実施例では20%)で一定に設定されてい
るので、燃料系順モデル22における燃料付着率推定部
22aは初期値を出力し、最終的に、この初期値に基づ
いてモデルベース制御部20で燃料噴射装置4の基本操
作量Mfnが算出され、燃料が実際に噴射されてエンジ
ンが運転されている。この状態で前述のように一方の運
転状態パラメータ(スロットル開度)が過渡的に変化さ
せられると、その運転状態点a1及びa2に出力パラメ
ータの初期値(20%)が適合していれば、目標空燃比
に対する排気空燃比のずれは許容値以上にならないが、
適合していないと図8に符号S1で示すように排気空燃
比が目標空燃比に対して許容値以上にずれる空燃比スパ
イクが生じる。モデル補正部60は、空燃比スパイクが
生じたときに、この空燃比スパイクを燃料付着率xとの
関連が強い特性と、蒸発時定数τとの関連が強い特性と
に分割し、各特性から燃料付着率推定部22a及び蒸発
時定数推定部22bにおけるモデル(即ち、マップ)の
教師データを獲得するように構成されている。ここで、
モデル化した燃料の動的挙動、即ち、燃料付着率x及び
蒸発時定数τの特性の違いについて説明すると、燃料付
着率xは燃料噴射装置4から噴射された燃料が直接シリ
ンダ内に入らずに吸気管壁面や吸気弁に付着する割合を
示すものなので、その変化は推定燃料量の変化として応
答性よく現れるが、蒸発時定数τは付着した燃料が除々
に蒸発してシリンダ内に入る速度を示すものなので、そ
の変化が推定燃料量に現れるは、燃料付着率xに比べて
時間がかかり、また、その量も小さい。モデル補正部6
0は、上記した燃料付着率x及び蒸発時定数τの特性の
違いを利用して、図9に示すようにエンジンの運転状態
が過渡的に変化して排気空燃比スパイクが発生した時
に、その排気空燃比スパイクを、ピーク値P(即ち、モ
デルのずれ量を応答性よく表す値)と目標空燃比に対す
るずれ面積S(即ち、モデルのずれ量を時間をかけて表
す値)とに分け、ピーク値Pを燃料付着率推定部22a
を構成するモデルの補正学習用教師データを獲得するパ
ラメータとし、ずれ面積Sを蒸発時定数推定部22bを
構成するモデルの補正学習用教師データを獲得するパラ
メータとして各教師データの獲得を行う。具体的には、
図10のフローチャートに示すように、モデル補正部6
0は空燃比センサ15から得られる排気空燃比Eの波形
を読み込み(ステップ1)、目標空燃比Epと排気空燃
比Eとの差をずれ量△Eとしてセットする(ステップ
2)。そして、前記ずれ量△Eが予め設定した許容値L
を越えるか否か(ステップ3)、及び空燃比スパイク発
生中か否か(ステップ4)を監視し、前記すれ量△Eが
許容値Lより大きくなるまでステップ1及び2の処理を
繰り返す。ずれ量△Eが許容値Lを越えると、空燃比ス
パイクが発生していると判断してスパイクフラグをセッ
トし(ステップ5)、ずれ量△Eをピーク値Pとして設
定する(ステップ7)。ステップ6でピーク値Pを設定
した後、その時のずれ量△Eをずれ面積Sとして設定し
(ステップ8)、再びステップ1の処理に戻る。上記し
たステップ7のピーク値Pの設定処理は、ずれ量△Eが
前回設定したピーク値Pより小さくなるまで繰り返し行
われ(ステップ6,7)、これにより最終的に真のピー
ク値Pが得られる。また、前記ステップ7のずれ面積S
の設定処理は、ずれ量△Eが予め設定した安定判別値よ
り小さくなるまで繰り返し行われ(ステップ9,8)、
この時に前回のずれ面積Sにその時のずれ量△Eを加算
した値を新たにずれ面積Sと設定することで最終的に真
のずれ面積Sが得られる。尚、前記安定判別値は、許容
値より小さく設定される。そして、ずれ量△Eが安定判
別値より小さくなると、空燃比スパイクが無くなったと
判断して(ステップ9)、スパイクフラグをクリアし
(ステップ10)、設定したピーク値P及びずれ面積S
から教師データの獲得を行う(ステップ11)。図11
は教師データの獲得処理を示すフローチャートであり、
このフローチャートに示すように、モデル補正部60
は、スパイクフラグをクリアした後、スロットル開度
α、エンジン回転数N及び吸気管壁温Teを読み込み
(ステップa)、これらの情報とピーク値P及びずれ面
積Sとから燃料付着率推定部22aのずれを補正する教
師データxtと蒸発時定数推定部22bのずれを補正す
る教師データτtを予め決められた演算式を用いて算出
する(ステップb)。具体的には、エンジンの過渡状態
がスロットル開度が増加して発生した場合には、次式
(1)及び(2)を用いて各教師データxt及びτtを
算出し、 xt = x + Gx ×(P−E) ・・・・(1) τt = τ + Gτ × S ・・・・(2) 前記過渡状態がスロットル開度が減少して発生した場合
には、次式(3)及び(4)を用いて各教師データxt
及びτtを算出する。 xt = x − Gx ×(P−E) ・・・・(3) τt = τ − Gτ × S ・・・・(4) 上記した式(1)〜(4)において、Gx,Gτはそれ
ぞれ教師データを獲得する際のゲインであり、Gx>
0,Gτ>0に設定される。上述の式により、具体的に
は、例えば、スロットル開度が増加している時にピーク
値が目標値より大きい場合には、燃料付着率推定部22
aで燃料付着率を実際の燃料付着率より小さい値に推定
していることになるので、獲得された教師データxt
は、学習前の燃料付着率推定部22aの出力値xより大
きい値になり、逆にピーク値が目標値より小さい場合に
は、推定した燃料付着率が実際の燃料付着率より大きい
ということになるので、獲得された教師データxTは、
学習前の燃料付着率推定部22aの出力値xより小さい
値になる。また、スロットル開度が減少している時にピ
ーク値が目標値より小さい場合には、燃料付着率推定部
22aで燃料付着率を実際の燃料付着率より小さい値に
推定していることになるので、獲得された教師データx
tは、学習前の燃料付着率推定部22aの出力値xより
大きい値になり、逆にピーク値が目標値より大きい場合
には、推定した燃料付着率が実際の燃料付着率より大き
いということになるので、獲得された教師データxT
は、学習前の燃料付着率推定部22aの出力値xより小
さい値になる。さらに、スロットル開度が増加している
時にスパイクが目標値より大きい方に発生した場合に
は、蒸発時定数推定部22bで蒸発時定数を実際の蒸発
時定数より小さい値に推定していることになるので、獲
得された教師データτtは、学習前の蒸発時定数推定部
22bの出力値τより大きい値になり、逆にスパイクが
目標値より小さい方に発生した場合には、推定した蒸発
時定数が実際の蒸発時定数より大きいということになる
ので、獲得された教師データτTは、学習前の蒸発時定
数推定部22bの出力値τより小さい値になる。さら
に、スロットル開度が減少している時にスパイクが目標
値より小さい方に発生した場合には、蒸発時定数推定部
22bで蒸発時定数を実際の蒸発時定数より小さい値に
推定していることになるので、獲得された教師データτ
tは、学習前の蒸発時定数推定部22bの出力値τより
大きい値になり、逆にスパイクが目標値より大きい方に
発生した場合には、推定した蒸発時定数が実際の蒸発時
定数より大きいということになるので、獲得された教師
データτTは、学習前の蒸発時定数推定部22bの出力
値τより小さい値になる。以上のようにしてモデル補正
部60は、運転状態点a1から運転状態点a2にスロッ
トル開度が変化した時の教師データxt及びτtを獲得
する。図7及び図8に戻って説明すると、前述のように
してモデル補正部60で得られた教師データxtは、運
転状態点a1(エンジン回転数1000rpm及びスロ
ットル開度0度)から運転状態点a2(エンジン回転数
1000rpm及びスロットル開度10度)に過渡的に
変化した時、即ち、スロットルが過渡的に開けられた時
の運転状態a1及び運転状態a2に対する教師データx
t(a1;a2)及びτt(a1;a2)として獲得さ
れ記憶される。図8の実施例では、運転状態点間の相互
の変化、即ち、運転状態点a1から運転状態点a2への
変化(スロットル開度が開く変化)と、運転状態点a2
から運転状態点a1への変化(スロットル開度が閉じる
変化)との両方の教師データxt(a1;a2)及びx
t(a2;a1)、τt(a1;a2)及びτt(a
2;a1)を獲得する。ここで、過渡的に変化した二つ
の運転状態点a1及びa2に対する燃料付着率x及び蒸
発時定数τの影響について簡単に説明する。例えば、燃
料付着率推定部22aの燃料付着率xが運転状態点a2
については実際の燃料付着率と同じ値を獲得している
が、運転状態点a1については実際の燃料付着率より小
さい値になっている場合、運転状態点a1での燃料付着
率xが実際の燃料付着率より小さい値になっているた
め、蒸発燃料量a・Mwも実際の蒸発燃料量よりも小さ
い値となる。この状態で、スロットルを運転状態点a1
から運転状態点a2に過渡的に変化させると、運転状態
点a2での燃料付着率xが合っていても、蒸発燃料量の
変化が一次遅れで遅いため、小さい値に推定された蒸発
燃料量a・Mwの不足分を直接シリンダ内に入る燃料量
で補うようにモデルベース制御部は燃料噴射量(基本操
作量)を増やしてしまい、実際に必要な燃料量よりも多
くシリンダ内に燃料が流入し、空燃比波形は運転状態点
2でリッチスパイクとなる。このように、運転状態点a
1での燃料付着率xの推定量のずれは、スロットルを運
転状態点a1から運転状態点a2に過渡的に変化させる
と、運転状態点a2での空燃比波形に影響を及ぼす。一
方、蒸発時定数については、例えば、スロットルを運転
状態点a1から運転状態点a2に過渡的に変化させる時
に、運転状態点a1での推定蒸発時定数τが実際の蒸発
時定数からずれていても付着分と蒸発分が等しくなって
いれば、空燃比波形には影響は与えないが、運転状態点
a2での推定蒸発時定数τがずれていると、過渡的に変
化させた時に推定蒸発燃料量の応答が変わり、それに合
わせて燃料噴射量を決定するので、空燃比波形に影響を
与えることになる。従って、燃料付着率推定部22aに
ついては、二つの運転状態点間の相互の変化で得られた
教師データxt(a1;a2)及びxt(a2;a1)
を加算した値を各運転状態点a1及びa2の教師データ
とし、蒸発時定数推定部22bについては、運転状態点
a1から運転状態点a2に変化させた時に得られた教師
データτt(a1;a2)を運転状態点a1の教師デー
タとし、反対に変化させた時に得られた教師データτt
(a2;a1)を運転状態点a2の教師データとしてそ
れぞれ獲得する。本実施例では、モデルをマップで構成
しているので、教師データに基づく各出力パラメータの
学習では、前記した各教師データの値がそのまま、対応
する運転状態点の新しい出力パラメータとして記憶され
る。上記した処理は、運転状態点a1及びa2間で運転
状態を相互に変化させた時に、目標空燃比に対する排気
空燃比のずれが許容値以下になるまで、繰り返し行わ
れ、その都度、教師データを獲得して出力パラメータを
更新していく。初期学習を行う場合、一方の運転状態パ
ラメータを固定して、他方の運転状態パラメータを運転
状態点間で過渡的に変化させるため、一方の運転状態パ
ラメータを固定したまま、他方の運転状態パラメータを
連続する運転状態間で次々に変化させていき、各運転状
態点における出力パラメータの設定を行っていった方が
好ましい。従って、運転状態点a1及びa2の出力パラ
メータの設定が終了した後は、図7における分割運転領
域bを画定する運転状態点a2及びa3間でスロットル
開度を過渡的に変化させて出力パラメータの設定を行
う。この時、運転状態点a2における出力パラメータ値
は既に設定されているが、運転状態点a3についても、
仮に運転状態点a2における出力パラメータ値に設定し
て処理を行うと空燃比の早期安定化が図れる。上記した
状態で、運転状態点a1及びa2の出力パラメータ値を
設定した処理と同じ処理を行い、モデル補正部60で最
終的な教師データxt(a2;a3)及びxt(a3;
a2)、τt(a2;a3)及びτt(a3;a2)を
獲得する。蒸発時定数推定部22bについては、獲得し
た教師データをそのまま各運転状態点a2及びa3の教
師データとして更新するが、燃料付着率xについては、
変化させる前後の運転状態点a2及びa3の両方に影響
を及ぼすため、一つの運転状態点a2に対して二つの教
師データxt(a1;a2)+xt(a2;a1)と、
xt(a2;a3)+xt(a3;a2)が獲得される
ことになる。モデル補正部60は、上記したように一つ
の運転状態点(例えば、a2)に対して異なる分割運転
領域(例えば、A及びB)に関する二つの教師データを
獲得すると、二つの教師データの平均値をその運転状態
点の教師データとして出力する。従って、この場合に
は、目標空燃比に対する排気空燃比のずれ量が許容値以
下になるまで運転状態点a2からa3の間でスロットル
開度を過渡的に変化させ、それにより得られた教師デー
タxt(a2;a3)+xt(a3:a2)は、運転状
態点a3に対してはそのまま出力されるが、運転状態点
a2に対しては、前回の運転状態点a2に対する教師デ
ータxt(a1;a2)+xt(a2;a1)と今回得
られた教師データxt(a2;a3)+xt(a3;a
2)との平均値が教師データXTとして出力される。以
上説明した処理を全ての運転状態点に対して行い、それ
により、各運転状態点に対する最適な出力パラメータ値
(燃料付着率x)が得られると、初期学習の処理は終了
し、図5に示すようなマップが得られる。
Hereinafter, the initial learning including acquisition and learning of teacher data will be described with a specific example.
First, two operation state points a1 that define the divided operation area a
In the following description, the acquisition of teacher data in steps a and a2 will be described as an example. A crank angle signal for initial learning is input to the control device 10 so that the engine speed calculated by the engine speed calculation unit 30 becomes 1000 rpm. The engine speed N at 10 is fixed to 1000 rpm including the operating state points a1 and a2. Then, a throttle opening signal (that is, a 0-degree signal) for initial learning relating to the operating state point a1 is input, and the throttle opening signal is transiently changed to the operating state point a2 (that is, a 10-degree signal). FIG.
In the initial state shown in (1), the output parameter value of the map is set to an appropriate initial value (20% in this embodiment) and is constant, so that the fuel adhesion rate estimating unit 22a in the fuel system forward model 22 outputs the initial value. Finally, the basic operation amount Mfn of the fuel injection device 4 is calculated by the model base control unit 20 based on this initial value, and the fuel is actually injected to operate the engine. In this state, if one of the operation state parameters (throttle opening) is transiently changed as described above, if the initial value (20%) of the output parameter matches the operation state points a1 and a2, Although the deviation of the exhaust air-fuel ratio from the target air-fuel ratio does not exceed the allowable value,
If not, an air-fuel ratio spike occurs in which the exhaust air-fuel ratio deviates from the target air-fuel ratio by more than an allowable value, as indicated by reference numeral S1 in FIG. When an air-fuel ratio spike occurs, the model correction unit 60 divides the air-fuel ratio spike into a characteristic strongly related to the fuel adhesion rate x and a characteristic strongly related to the evaporation time constant τ. It is configured to acquire teacher data of a model (that is, a map) in the fuel adhesion rate estimation unit 22a and the evaporation time constant estimation unit 22b. here,
Explaining the dynamic behavior of the modeled fuel, that is, the difference between the characteristics of the fuel adhesion rate x and the evaporation time constant τ, the fuel adhesion rate x is such that the fuel injected from the fuel injection device 4 does not directly enter the cylinder. Since it shows the rate of adhering to the intake pipe wall and the intake valve, the change appears with a good response as a change in the estimated fuel amount, but the evaporation time constant τ is the speed at which the adhering fuel evaporates gradually and enters the cylinder. Since the change appears in the estimated fuel amount, it takes a longer time than the fuel adhesion rate x, and the amount is small. Model correction unit 6
0 uses the above-mentioned difference in the characteristics of the fuel adhesion rate x and the evaporation time constant τ, when the operating state of the engine changes transiently and an exhaust air-fuel ratio spike occurs as shown in FIG. The exhaust air-fuel ratio spike is divided into a peak value P (that is, a value representing the model deviation with good responsiveness) and a deviation area S to the target air-fuel ratio (that is, a value representing the model deviation over time). The peak value P is used as the fuel adhesion rate estimation unit 22a.
Are obtained as parameters for acquiring the correction learning teacher data of the model constituting, and the deviation area S is acquired as the parameter for acquiring the correction learning teacher data of the model constituting the evaporation time constant estimating unit 22b. In particular,
As shown in the flowchart of FIG.
A value of 0 reads the waveform of the exhaust air-fuel ratio E obtained from the air-fuel ratio sensor 15 (step 1), and sets the difference between the target air-fuel ratio Ep and the exhaust air-fuel ratio E as a deviation ΔE (step 2). The deviation amount ΔE is set to a predetermined allowable value L.
Are monitored (step 3) and whether an air-fuel ratio spike is occurring (step 4), and the processing of steps 1 and 2 is repeated until the slip amount ΔE becomes larger than the allowable value L. If the deviation ΔE exceeds the allowable value L, it is determined that an air-fuel ratio spike has occurred, a spike flag is set (step 5), and the deviation ΔE is set as a peak value P (step 7). After setting the peak value P in step 6, the shift amount ΔE at that time is set as the shift area S (step 8), and the process returns to step 1 again. The above-described processing of setting the peak value P in step 7 is repeatedly performed until the deviation amount ΔE becomes smaller than the previously set peak value P (steps 6 and 7), whereby the true peak value P is finally obtained. Can be In addition, the displacement area S of the step 7
Is repeatedly performed until the deviation ΔE becomes smaller than a preset stability determination value (steps 9 and 8),
At this time, a value obtained by adding the shift amount ΔE at that time to the previous shift area S is newly set as the shift area S, so that the true shift area S is finally obtained. Note that the stability determination value is set smaller than an allowable value. When the deviation ΔE is smaller than the stability determination value, it is determined that the air-fuel ratio spike has disappeared (step 9), the spike flag is cleared (step 10), and the set peak value P and deviation area S
Is obtained (step 11). FIG.
Is a flowchart showing a process of acquiring teacher data;
As shown in this flowchart, the model correction unit 60
Reads the throttle opening α, the engine speed N, and the intake pipe wall temperature Te after clearing the spike flag (step a), and from the information, the peak value P, and the deviation area S, the fuel adhesion rate estimation unit 22a The teacher data xt for correcting the shift and the teacher data τt for correcting the shift of the evaporation time constant estimating unit 22b are calculated using a predetermined arithmetic expression (step b). Specifically, when the transient state of the engine occurs due to an increase in the throttle opening, the teacher data xt and τt are calculated using the following equations (1) and (2), and xt = x + Gx × (P−E) (1) τt = τ + Gτ × S (2) When the transient state occurs due to a decrease in throttle opening, the following equation (3) and Using (4), each teacher data xt
And τt are calculated. xt = x−Gx × (PE) (3) τt = τ−Gτ × S (4) In the above equations (1) to (4), Gx and Gτ are teachers, respectively. This is the gain when acquiring data, and Gx>
0, Gτ> 0. According to the above equation, specifically, for example, when the peak value is larger than the target value when the throttle opening is increasing, the fuel adhesion rate estimating unit 22
a, the fuel deposition rate is estimated to be smaller than the actual fuel deposition rate, so that the acquired teacher data xt
Is larger than the output value x of the fuel adhesion rate estimating unit 22a before learning. Conversely, if the peak value is smaller than the target value, the estimated fuel adhesion rate is larger than the actual fuel adhesion rate. Therefore, the acquired teacher data xT is
The value becomes smaller than the output value x of the fuel adhesion rate estimation unit 22a before learning. If the peak value is smaller than the target value when the throttle opening is decreasing, the fuel adhesion rate estimating unit 22a estimates the fuel adhesion rate to a value smaller than the actual fuel adhesion rate. , Acquired teacher data x
t is a value larger than the output value x of the fuel adhesion rate estimating unit 22a before learning, and conversely, if the peak value is larger than the target value, the estimated fuel adhesion rate is larger than the actual fuel adhesion rate. , The acquired teacher data xT
Is smaller than the output value x of the fuel adhesion rate estimation unit 22a before learning. Further, if a spike occurs in a direction larger than the target value while the throttle opening is increasing, the evaporation time constant estimating unit 22b estimates the evaporation time constant to a value smaller than the actual evaporation time constant. Therefore, the acquired teacher data τt becomes a value larger than the output value τ of the evaporation time constant estimating unit 22b before learning, and conversely, if a spike occurs below the target value, the estimated evaporation data Since the time constant is larger than the actual evaporation time constant, the acquired teacher data τT is smaller than the output value τ of the evaporation time constant estimator 22b before learning. Furthermore, if a spike occurs below the target value while the throttle opening is decreasing, the evaporation time constant estimating unit 22b estimates the evaporation time constant to a value smaller than the actual evaporation time constant. And the acquired teacher data τ
t is a value larger than the output value τ of the evaporation time constant estimating unit 22b before learning. Conversely, if the spike occurs in a direction larger than the target value, the estimated evaporation time constant is smaller than the actual evaporation time constant. Therefore, the acquired teacher data τT is smaller than the output value τ of the evaporation time constant estimator 22b before learning. As described above, the model correction unit 60 acquires the teacher data xt and τt when the throttle opening changes from the operating state point a1 to the operating state point a2. Referring back to FIGS. 7 and 8, the teacher data xt obtained by the model correction unit 60 as described above is calculated from the operating state point a1 (the engine speed 1000 rpm and the throttle opening 0 degree) to the operating state point a2. (Engine speed 1000 rpm and throttle opening 10 degrees), that is, teacher data x for the operating state a1 and the operating state a2 when the throttle is transiently opened.
acquired and stored as t (a1; a2) and τt (a1; a2). In the embodiment of FIG. 8, the mutual change between the operating state points, that is, the change from the operating state point a1 to the operating state point a2 (the change in the opening of the throttle opening) and the operating state point a2
From the operating state point a1 to the operating state point a1 (change of the throttle opening closing) and xt (a1; a2) and x
t (a2; a1), τt (a1; a2) and τt (a
2; a1) is obtained. Here, the influence of the fuel adhesion rate x and the evaporation time constant τ on the two transiently changed operating state points a1 and a2 will be briefly described. For example, the fuel adhesion rate x of the fuel adhesion rate estimating unit 22a is equal to the operating state point a2.
Has the same value as the actual fuel deposition rate, but when the operating state point a1 is smaller than the actual fuel deposition rate, the fuel deposition rate x at the operating state point a1 is Since the value is smaller than the fuel adhesion rate, the amount of evaporated fuel a · Mw is also smaller than the actual amount of evaporated fuel. In this state, the throttle is moved to the operating state point a1.
Is transiently changed from the operating state point a2 to the operating state point a2, even if the fuel adhering rate x at the operating state point a2 matches, the change in the amount of evaporative fuel is slow with a first-order lag. The model-based control unit increases the fuel injection amount (basic operation amount) so that the shortage of a · Mw is compensated for by the fuel amount directly entering the cylinder, and the fuel in the cylinder becomes larger than the actually required fuel amount. Then, the air-fuel ratio waveform becomes a rich spike at the operating state point 2. Thus, the operating state point a
The deviation of the estimated amount of the fuel adhesion rate x at 1 affects the air-fuel ratio waveform at the operating state point a2 when the throttle is transiently changed from the operating state point a1 to the operating state point a2. On the other hand, regarding the evaporation time constant, for example, when the throttle is transiently changed from the operation state point a1 to the operation state point a2, the estimated evaporation time constant τ at the operation state point a1 deviates from the actual evaporation time constant. If the adhesion amount and the evaporation amount are equal, the air-fuel ratio waveform is not affected, but if the estimated evaporation time constant τ at the operating state point a2 is shifted, the estimated evaporation time Since the response of the fuel amount changes and the fuel injection amount is determined in accordance therewith, the air-fuel ratio waveform is affected. Therefore, for the fuel adhesion rate estimating unit 22a, the teacher data xt (a1; a2) and xt (a2; a1) obtained by the mutual change between the two operating state points.
Is used as teacher data for each of the operating state points a1 and a2, and the evaporating time constant estimating unit 22b obtains teacher data τt (a1; a2) obtained when the operating state point a1 is changed to the operating state point a2. ) Is used as the teacher data of the operating state point a1, and the teacher data τt obtained when it is changed conversely.
(A2; a1) are acquired as teacher data of the operating state point a2. In the present embodiment, since the model is configured by a map, in learning each output parameter based on teacher data, the value of each teacher data is stored as it is as a new output parameter of the corresponding operating state point. The above-described processing is repeatedly performed until the deviation of the exhaust air-fuel ratio from the target air-fuel ratio becomes equal to or less than the allowable value when the operation state is changed between the operation state points a1 and a2. Obtain and update output parameters. When performing the initial learning, one operating state parameter is fixed, and the other operating state parameter is transiently changed between operating state points. It is preferable to change one after another between successive operating states and to set output parameters at each operating state point. Accordingly, after the setting of the output parameters of the operating state points a1 and a2 is completed, the throttle opening is transiently changed between the operating state points a2 and a3 that define the divided operating area b in FIG. Make settings. At this time, the output parameter value at the operating state point a2 has already been set.
If the processing is performed by setting the output parameter value at the operation state point a2, the air-fuel ratio can be stabilized early. In the above-described state, the same processing as the processing for setting the output parameter values of the operating state points a1 and a2 is performed, and the model correction unit 60 finally obtains the teacher data xt (a2; a3) and xt (a3;
a2), τt (a2; a3) and τt (a3; a2) are obtained. The evaporation time constant estimating unit 22b updates the acquired teaching data as it is as the teaching data of each of the operating state points a2 and a3.
In order to affect both the operating state points a2 and a3 before and after the change, two teaching data xt (a1; a2) + xt (a2; a1) for one operating state point a2,
xt (a2; a3) + xt (a3; a2) will be obtained. As described above, when the model correction unit 60 acquires two pieces of teacher data regarding different operation areas (for example, A and B) for one operation state point (for example, a2), the average value of the two pieces of teacher data is obtained. Is output as teacher data of the operating state point. Therefore, in this case, the throttle opening is transiently changed between the operating state points a2 and a3 until the deviation amount of the exhaust air-fuel ratio from the target air-fuel ratio becomes equal to or less than the allowable value. xt (a2; a3) + xt (a3: a2) is output as it is for the operating state point a3, but for the operating state point a2, the teacher data xt (a1; a2) + xt (a2; a1) and the newly obtained teacher data xt (a2; a3) + xt (a3; a
2) is output as teacher data XT. The process described above is performed for all the operating state points, and when the optimum output parameter value (fuel adhesion rate x) for each operating state point is obtained, the initial learning processing ends and FIG. A map as shown is obtained.

【0014】(オンライン学習について)全ての初期学
習が終了すると、制御装置10は、エンジンに実装され
て使用される。しかし、上記したように予め初期学習で
全ての出力パラメータを最適値に設定しても、内燃機関
駆動後に、内燃機関の経時変化や内燃機関を使用する環
境の変化等により基礎データと実際の内燃機関における
動的挙動とのずれが生じる。従って、モデル補正部60
は、エンジンに実装されて制御装置10が作動している
間も、各検知手段12,13,14からの情報と空燃比
センサ15に基づいて実際の運転状態パラメータ(エン
ジン回転数及びスロットル開度)と、排気空燃比とを監
視し、運転状態が過渡的に変化して空燃比スパイクが生
じた時に、過渡的に変化する前後の各運転状態に一番近
い運転状態点に対する教師データを前述の初期学習と同
じ方法で獲得して各運転状態点に対してオンラインでの
補正学習を行わせる。
(On-line learning) When all initial learning is completed, the control device 10 is mounted on an engine and used. However, even if all the output parameters are set to the optimum values in the initial learning as described above, the basic data and the actual internal combustion engine may be changed after the internal combustion engine is driven due to a change over time of the internal combustion engine or a change in the environment in which the internal combustion engine is used. Deviation from the dynamic behavior in the engine occurs. Therefore, the model correction unit 60
While the control device 10 is mounted on the engine and the control device 10 is operating, the actual operating state parameters (engine speed and throttle opening degree) are determined based on the information from the detecting means 12, 13 and 14 and the air-fuel ratio sensor 15. ) And the exhaust air-fuel ratio are monitored, and when the operating state changes transiently and an air-fuel ratio spike occurs, the teaching data for the operating state point closest to each operating state before and after the transient change is described above. In the same manner as in the initial learning, and online correction learning is performed for each operating state point.

【0015】上記したように、モデルを構成する運転状
態パラメータの全領域を複数の分割運転領域に分け、各
分割運転領域を画定する二つの運転状態点間で過渡的に
運転状態を変化させ、その時に生じる目標空燃比に対す
る排気空燃比のずれ量から、過渡的に変化させた二つの
運転状態点におけるモデルの出力パラメータを補正する
教師データを獲得すると共に、一つの運転状態点に対し
て二つの教師データを獲得した場合には、それらの平均
値を教師データとして獲得することでモデルの初期学習
を行うことにより、複雑な演算なしに、運転状態パラメ
ータの全領域で出力パラメータが滑らかに変化するモデ
ルを得ることができるという効果を奏する。また、各運
転状態点における教師データの獲得を、目標空燃比に対
する排気空燃比のずれ量が許容値以上になってから、安
定判別値以下になるまでの、前記ずれ量の特性をピーク
値及び積分値に分け、これらを用いて燃料付着率推定部
22a及び蒸発時定数推定部22bにおけるモデルのず
れを補正するための教師データを獲得するように構成す
ることにより、使用中に前記モデルを、適用するエンジ
ンの個体差、使用環境の変化、経時変化等に追従してそ
れらの適合させることができ、煩雑な作業なしに常に最
適なモデルで空燃比制御を行うことができるようにな
る。
As described above, the entire operating condition parameter region constituting the model is divided into a plurality of divided operating regions, and the operating state is transiently changed between two operating state points defining each divided operating region. From the deviation amount of the exhaust air-fuel ratio with respect to the target air-fuel ratio generated at that time, the teacher data for correcting the output parameters of the model at the two operating state points that have been transiently changed is obtained. When two sets of teacher data are acquired, the initial value of the model is obtained by acquiring the average value of them as teacher data, so that the output parameters change smoothly in the entire range of the operating state parameters without complicated calculations. This has the effect of obtaining a model that performs Further, the acquisition of the teacher data at each operating state point is based on the characteristic of the deviation amount between the deviation amount of the exhaust air-fuel ratio with respect to the target air-fuel ratio from the allowable value to the stability determination value and the peak value and the characteristic. By dividing the integral model and using them to acquire teacher data for correcting the model deviation in the fuel adhesion rate estimating section 22a and the evaporation time constant estimating section 22b, the model can be used during use. It is possible to adapt to individual differences of the engine to be applied, changes in the use environment, changes over time, and the like, and to always perform air-fuel ratio control with an optimal model without complicated work.

【0016】尚、上記したずれ量の特性の分割は、本実
施例に限定されることなく、変化が推定燃料量の変化と
して応答性よく現れるという燃料付着率の特性と、変化
が推定燃料量に現れるに時間がかかるという蒸発時定数
の特性とが反映される特性に分割できれば任意の分割方
法でよく、例えば、図12(a)〜(c)に示すよう
に、ずれ量を早期に現れる特性(即ち、前半部の積分
値)と遅く現れる特性(即ち、後半部の積分値)とに分
割して前半部の積分値に基づいて燃料付着率の教師デー
タを、後半部の積分値に基づいて蒸発時定数の教師デー
タを獲得するようにしてもよい。図12(a)は、ずれ
量の特性の分割の別の例を示す図であり、この例では、
ずれ量を、ピーク値Pを境に前半部と後半部とに分割し
ており、各々の積分値S1及びS2に基づいて各教師デ
ータの獲得を行う。図12(b)は、ずれ量の特性の分
割のさらに別の例を示す図であり、この例では、ずれ
(スパイク)発生時間を予め決められた比率(例えば、
1対3)で前半部と後半部とに分割しており、各々の積
分値S1及びS2に基づいて各教師データの獲得を行
う。図12(c)は、ずれ量の特性の分割のさらに別の
例を示す図であり、この例では、ずれ(スパイク)発生
時から所定の時間(例えば、100ms)を境にずれ量
を前半部と後半部とに分割しており、各々の積分値S1
及びS2に基づいて各教師データの獲得を行う。
The above-described division of the deviation amount characteristic is not limited to the present embodiment. The characteristic of the fuel adhesion rate, in which the change appears as a change in the estimated fuel amount with good responsiveness, and the change in the characteristic of the estimated fuel amount, Any method can be used as long as it can be divided into characteristics that take into account the characteristics of the evaporation time constant, which takes time to appear, as shown in FIGS. 12 (a) to 12 (c). The characteristic data (that is, the integral value of the first half) and the characteristic that appears late (that is, the integral value of the second half) are divided, and the teacher data of the fuel adhesion rate is converted into the integral value of the second half based on the integral value of the first half. Based on this, teacher data of the evaporation time constant may be obtained. FIG. 12A is a diagram showing another example of division of the characteristic of the shift amount. In this example, FIG.
The shift amount is divided into a first half and a second half with a peak value P as a boundary, and each teacher data is obtained based on the respective integral values S1 and S2. FIG. 12B is a diagram showing still another example of the division of the characteristic of the shift amount. In this example, the shift (spike) occurrence time is set to a predetermined ratio (for example,
It is divided into the first half and the second half by (1: 3), and each teacher data is obtained based on the respective integral values S1 and S2. FIG. 12C is a diagram showing still another example of the division of the characteristic of the shift amount. In this example, the shift amount is set in the first half after a predetermined time (for example, 100 ms) from the occurrence of the shift (spike). Part and the latter half part, and each integral value S1
And each teacher data is acquired based on and S2.

【0017】尚、以上説明した実施例では、各教師デー
タxt及びτtをスロットル開度の増加及び減少に対応
した二種類の演算式で求めているが、これは本実施例に
限定されることなく、例えば、スロットル開度の微分値
dαを算出し、この微分値dαを用いて、各教師データ
毎の二種類の式を次式(5)及び(6)のように統一し
てもよい。 xt = x + Kx × dα ×(P−E) ・・・・(5) τt = τ + Kτ × dα × S ・・・・(6)
In the embodiment described above, the teacher data xt and τt are obtained by two types of arithmetic expressions corresponding to the increase and decrease of the throttle opening. However, this is not limited to this embodiment. Instead, for example, the differential value dα of the throttle opening may be calculated, and using the differential value dα, two types of equations for each teacher data may be unified as in the following equations (5) and (6). . xt = x + Kx × dα × (PE) (5) τt = τ + Kτ × dα × S (6)

【0018】また、以上説明した実施例では、ピーク値
が一つしかない空燃比スパイクを例に挙げて教師データ
の獲得処理について説明しているが、これは本実施例に
限定されることなく、ピーク値が複数ある空燃比スパイ
クが発生した場合には、適当な空燃比スパイク(例え
ば、一番始めのピーク値又は最大ピーク値)が選択され
得る。
In the above-described embodiment, the process of acquiring teacher data is described using an air-fuel ratio spike having only one peak value as an example. However, this is not limited to this embodiment. When an air-fuel ratio spike having a plurality of peak values occurs, an appropriate air-fuel ratio spike (for example, the first peak value or the maximum peak value) can be selected.

【0019】また、上記した実施例では、各運転状態点
における出力パラメータに対する修正値をそのまま教師
データとして獲得するよう説明しているが、これは本実
施例に限定されることなく、例えば、図13に示すよう
に出力パラメータに対する修正幅を教師データとして獲
得し、各出力パラメータを更新していってもよい。さら
に、上記した実施例では、一つの運転状態点に対して、
異なる分割運転領域に関する二つの教師データを獲得し
た場合には、それらの平均値を教師データとして獲得す
るよう説明しているが、これは本実施例に限定されるこ
となく、獲得した二つの教師データの中間値であれば任
意の値から教師データは獲得され得る。さらにまた、本
実施例では、エンジンの運転状態を等間隔に分割してい
るが、これは本実施例に限定されることなく任意に分割
され得る。また、本実施例では、マップを用いて燃料の
動的挙動のモデル化を行っているが、モデル化方法は本
実施例に限定されることなく、学習可能なフィードフォ
ワード制御ロジックを用いた方法であれば、任意のモデ
ル化方法でよく、例えば、ファジィニューラル回路網や
ニューラル回路網等が用いられ得る。尚、ファジィニュ
ーラル回路網を用いる場合には、各運転状態点に対する
最適なファジィニューラル回路網の出力値を教師データ
として獲得し、それらの教師データに基づいて、ファジ
ィニューラル回路網に用いられる各層間の適当な結合係
数の初期学習又はオンライン学習が行われ得る。
Further, in the above-described embodiment, the explanation has been made such that the correction value for the output parameter at each operation state point is acquired as it is as the teacher data. However, this is not limited to the present embodiment. As shown in FIG. 13, the correction range for the output parameter may be acquired as teacher data, and each output parameter may be updated. Further, in the above-described embodiment, for one operating state point,
It has been described that when two sets of teacher data relating to different divided driving regions are acquired, an average value thereof is acquired as teacher data. However, this is not limited to the present embodiment. The teacher data can be obtained from any value as long as it is an intermediate value of the data. Furthermore, in the present embodiment, the operating state of the engine is divided at equal intervals, but this may be arbitrarily divided without being limited to the present embodiment. Further, in the present embodiment, the dynamic behavior of the fuel is modeled using the map. However, the modeling method is not limited to the present embodiment, and the method using the feedforward control logic that can be learned is used. If so, any modeling method may be used. For example, a fuzzy neural network or a neural network may be used. When a fuzzy neural network is used, an output value of the optimal fuzzy neural network for each operating state point is acquired as teacher data, and based on the teacher data, each layer used in the fuzzy neural network is obtained. Initial learning or online learning of the appropriate coupling coefficient of.

【0020】[0020]

【発明の効果】以上説明したように、本発明に係る内燃
機関の空燃比制御装置における初期学習方法は、少なく
とも学習機能付きフィードフォワード制御ロジックを用
いて、エンジンの運転状態に関するパラメータに対する
燃料噴射装置から噴射された燃料の動的挙動を燃料付着
率としてモデル化した燃料系順モデルを有する内燃機関
の順モデルを用いてエンジンの排気空燃比を目標とする
空燃比に合わせて制御する内燃機関の空燃比制御装置に
おいて、前記運転状態に関するパラメータの全域を複数
の分割運転状態領域に分割し、一つの分割運転状態領域
を画定する運転状態点間で運転状態に関するパラメータ
を過渡的に変化させ、その時に生じる排気空燃比の目標
空燃比に対するずれに基づいて、変化させた二つの運転
状態点に対する燃料付着率の教師データを獲得し、前記
教師データに基づいてモデルの前記各運転状態点に対す
る燃料付着率及び蒸発時定数の学習を行うと共に、一つ
の運転状態点に対して異なる分割領域間に関する二つの
異なる教師データを獲得した場合には、二つの教師デー
タの中間値を、その対応する運転状態点の教師データと
して獲得するので、運転状態を過渡的に変化させた時に
その変化の前後の運転状態点に影響を及ぼす燃料付着率
の最適値を複雑な演算なしに求めることができるように
なり、かつ、運転状態全域にわたって滑らかに変化する
燃料付着率を初期学習終了後に得ることができる。従っ
て、複雑なパラメータの初期設定をより正確に、かつ迅
速に行うことができるようになるという効果を奏する。
As described above, the initial learning method in the air-fuel ratio control apparatus for an internal combustion engine according to the present invention uses at least a feedforward control logic with a learning function to control a fuel injection device for a parameter relating to an engine operating state. Of an internal combustion engine that controls the exhaust air-fuel ratio of the engine to the target air-fuel ratio using a forward model of the internal combustion engine that has a fuel system forward model that models the dynamic behavior of the fuel injected from the engine as the fuel adhesion rate In the air-fuel ratio control device, the entire area of the parameters relating to the operating state is divided into a plurality of divided operating state areas, and the parameters relating to the operating state are transiently changed between operating state points that define one divided operating state area. Based on the deviation of the exhaust air-fuel ratio from the target air-fuel ratio, Obtain teacher data of the sticking rate, learn the fuel sticking rate and the evaporation time constant for each of the operating state points of the model based on the teacher data, and perform two learnings between different divided regions for one operating state point. When two different teacher data are acquired, the intermediate value of the two teacher data is acquired as the teacher data of the corresponding operating state point, so when the operating state is transiently changed, the driving before and after the change is obtained. The optimum value of the fuel adhesion rate affecting the state point can be obtained without complicated calculation, and the fuel adhesion rate that smoothly changes over the entire operating state can be obtained after the completion of the initial learning. Therefore, there is an effect that the initial setting of complicated parameters can be performed more accurately and promptly.

【図面の簡単な説明】[Brief description of the drawings]

【図1】 エンジン1と本発明に係る予測学習方法を用
いてエンジン1の空燃比制御を行う制御装置10との関
係を示す概略図である。
FIG. 1 is a schematic diagram showing a relationship between an engine 1 and a control device 10 that controls an air-fuel ratio of the engine 1 using a prediction learning method according to the present invention.

【図2】 制御装置10の構成を示す概略ブロック図で
ある。
FIG. 2 is a schematic block diagram illustrating a configuration of a control device 10.

【図3】 図2におけるモデルベース制御部20の構成
を示す概略ブロック図である。
FIG. 3 is a schematic block diagram showing a configuration of a model base control unit 20 in FIG.

【図4】 燃料系順モデル22の構成を示す概略ブロッ
ク図である。
FIG. 4 is a schematic block diagram showing a configuration of a fuel system forward model 22.

【図5】 燃料系のモデルイメージ図である。FIG. 5 is a model image diagram of a fuel system.

【図6】 燃料付着率推定部を構成するマップを示して
いる。
FIG. 6 shows a map constituting a fuel adhesion rate estimation unit.

【図7】 燃料付着率推定部を構成するマップの初期学
習前の初期状態を示している。
FIG. 7 shows an initial state of a map constituting a fuel adhesion rate estimation unit before initial learning.

【図8】 初期学習実行中のスロットル開度、空燃比波
形及び燃料付着率(教師データ)の関係を示すグラフで
ある。
FIG. 8 is a graph showing a relationship among a throttle opening, an air-fuel ratio waveform, and a fuel adhesion rate (teacher data) during execution of initial learning.

【図9】 過渡時の排気空燃比の変化の詳細の一例を示
す図である。
FIG. 9 is a diagram showing an example of details of a change in an exhaust air-fuel ratio during a transition.

【図10】 モデル補正部の処理の流れを示すフローチ
ャートである。
FIG. 10 is a flowchart illustrating a flow of processing of a model correction unit.

【図11】 モデル補正部における教師データの獲得処
理の流れを示すフローチャートである。
FIG. 11 is a flowchart illustrating a flow of teacher data acquisition processing in a model correction unit.

【図12】 (a)〜(c)は過渡時の排気空燃比ずれ
量の成分分割の別の例を各々示している。
FIGS. 12A to 12C show other examples of component division of the amount of deviation of the exhaust air-fuel ratio during a transition.

【図13】 初期学習の別の実施例を示す図8に対応す
る図である。
FIG. 13 is a diagram corresponding to FIG. 8, showing another embodiment of the initial learning.

【符号の説明】[Explanation of symbols]

1 エンジン 2 吸気管 3 エアクリーナ 4 燃料噴射装置 5 シリンダ 6 排気管 7 クランクケース 8 スロットルバルブ 10 制御装置 12 スロットル開度検知手段 13 クランク角検知手段 14 吸気管壁温検知手段 15 空燃比センサ 20 モデルベース制御部 21 空気系順モデル 22 燃料系順モデル 22a 燃料付着率推定部 22b 蒸発時定数推定部 22c 非付着燃料演算部 22d 付着燃料演算 22e 蒸発燃料量演算部 23 推定空燃比演算部 24 基本操作量演算部 30 エンジン回転数演算部 40 目標空燃比算出部 50 変換部 60 モデル補正部 N エンジン回転数 α スロットル開度 Te 吸気管壁温 E 制御量 Ep 目標制御量 Ev 推定制御量 Mf 操作量 Av 推定空気量 Fv 推定燃料量 x 燃料付着率 τ 推定蒸発時定数 xt 教師データ(燃料付着率推定部) τt 教師データ(蒸発時定数推定部) DESCRIPTION OF SYMBOLS 1 Engine 2 Intake pipe 3 Air cleaner 4 Fuel injection device 5 Cylinder 6 Exhaust pipe 7 Crankcase 8 Throttle valve 10 Control device 12 Throttle opening detection means 13 Crank angle detection means 14 Intake pipe wall temperature detection means 15 Air-fuel ratio sensor 20 Model base Control unit 21 Air system forward model 22 Fuel system forward model 22a Fuel adhesion rate estimation unit 22b Evaporation time constant estimation unit 22c Non-adhesion fuel calculation unit 22d Adhesion fuel calculation 22e Evaporation fuel amount calculation unit 23 Estimated air-fuel ratio calculation unit 24 Basic operation amount Calculation unit 30 Engine rotation speed calculation unit 40 Target air-fuel ratio calculation unit 50 Conversion unit 60 Model correction unit N Engine rotation speed Throttle opening degree Te Intake pipe wall temperature E Control amount Ep Target control amount Ev Estimated control amount Mf Operation amount Av Estimation Air volume Fv Estimated fuel volume x Fuel adhesion rate τ Estimated evaporation time constant xt teacher data (fuel adhesion rate estimation unit) τt teacher data (evaporation time constant estimation unit)

Claims (2)

【特許請求の範囲】[Claims] 【請求項1】 燃料噴射装置から噴射された燃料の動的
挙動を少なくとも燃料付着率としてモデル化し、該モデ
ルを用いて排気空燃比を制御する内燃機関の空燃比制御
装置において、 前記運転状態に関するパラメータを複数領域に割り付け
て、 一つの領域を画定する運転状態点間で運転条件を過渡的
に変化させ、 その時に生じる排気空燃比の目標空燃比に対するずれに
基づいて、変化させた二つの運転状態点の燃料付着率の
教師データを獲得し、 前記教師データに基づいてモデルの学習を行うと共に、 同一の運転状態点に対して異なる領域に関する複数の異
なる教師データを獲得した場合には、それら複数の教師
データの中間値を、その対応する運転状態点の教師デー
タとして獲得することを特徴とする内燃機関の空燃比制
御装置における学習方法。
1. An air-fuel ratio control device for an internal combustion engine that models a dynamic behavior of fuel injected from a fuel injection device as at least a fuel adhesion rate and controls an exhaust air-fuel ratio using the model. The operating conditions are transiently changed between operating state points that define one area by allocating parameters to multiple areas, and the two operating modes are changed based on the deviation of the exhaust air-fuel ratio from the target air-fuel ratio that occurs at that time. When acquiring the teacher data of the fuel adhesion rate at the state point, learning the model based on the teacher data, and acquiring a plurality of different teacher data regarding different regions for the same driving state point, An air-fuel ratio control device for an internal combustion engine, wherein an intermediate value of a plurality of teacher data is obtained as teacher data of a corresponding operating state point.習方 method.
【請求項2】 前記中間値が、複数の教師データの平均
値であることを特徴とする請求項1に記載の学習方法。
2. The learning method according to claim 1, wherein the intermediate value is an average value of a plurality of teacher data.
JP10132197A 1998-05-14 1998-05-14 Learning method in air-fuel ratio control device of internal combustion engine Pending JPH11324782A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP10132197A JPH11324782A (en) 1998-05-14 1998-05-14 Learning method in air-fuel ratio control device of internal combustion engine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP10132197A JPH11324782A (en) 1998-05-14 1998-05-14 Learning method in air-fuel ratio control device of internal combustion engine

Publications (1)

Publication Number Publication Date
JPH11324782A true JPH11324782A (en) 1999-11-26

Family

ID=15075675

Family Applications (1)

Application Number Title Priority Date Filing Date
JP10132197A Pending JPH11324782A (en) 1998-05-14 1998-05-14 Learning method in air-fuel ratio control device of internal combustion engine

Country Status (1)

Country Link
JP (1) JPH11324782A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7124015B2 (en) 2002-07-05 2006-10-17 Honda Giken Kogyo Kabushiki Kaisha Control apparatus for plant
JP2010209835A (en) * 2009-03-11 2010-09-24 Denso Corp Exhaust emission control device for internal combustion engine
JP2014096145A (en) * 2012-10-22 2014-05-22 Boeing Co Real time control system management
JP2022072037A (en) * 2020-10-29 2022-05-17 株式会社ジャパンエンジンコーポレーション Main engine control system

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7124015B2 (en) 2002-07-05 2006-10-17 Honda Giken Kogyo Kabushiki Kaisha Control apparatus for plant
JP2010209835A (en) * 2009-03-11 2010-09-24 Denso Corp Exhaust emission control device for internal combustion engine
JP2014096145A (en) * 2012-10-22 2014-05-22 Boeing Co Real time control system management
JP2022072037A (en) * 2020-10-29 2022-05-17 株式会社ジャパンエンジンコーポレーション Main engine control system

Similar Documents

Publication Publication Date Title
JP3703117B2 (en) Model-based control method and apparatus
KR0160396B1 (en) Air-fuel ratio control system
US5954783A (en) Engine control system using combination of forward model and inverse model
JP2551038B2 (en) Air-fuel ratio control device for internal combustion engine
JP4144272B2 (en) Fuel injection amount control device for internal combustion engine
EP0967534B1 (en) Online learning method
US20040158387A1 (en) Control system and method, and control unit
JPH02104930A (en) Device for controlling fuel injection of internal combustion engine
US6840214B2 (en) Air-fuel ratio control apparatus for internal combustion engine
JPH11231907A (en) Method for learning prediction in feedforward control using learnable control logic
US6466859B1 (en) Control system
JPH10184431A (en) Engine control system
US6591822B2 (en) Air-fuel ratio controller of internal combustion engines
JPH11324782A (en) Learning method in air-fuel ratio control device of internal combustion engine
JPH06323181A (en) Method and device for controlling fuel in internal combustion engine
JPH10184429A (en) Engine control system
JP3508328B2 (en) Air-fuel ratio control device
JP2002309990A (en) Control device for internal combustion engine
JPH11353006A (en) On-line learning method
JPH0874636A (en) Air-fuel ratio control device
JPH11294230A (en) Fuel injection control device for engine
JPH10274082A (en) Engine control system
JPH08232751A (en) Suction air amount estimating device for internal combustion engine
JPH11229931A (en) Air-fuel ratio controlling method and air-fuel ratio controller for internal combustion engine
JP4740272B2 (en) Intake control device for internal combustion engine