JP2689362B2 - Air-fuel ratio detection method for internal combustion engine - Google Patents

Air-fuel ratio detection method for internal combustion engine

Info

Publication number
JP2689362B2
JP2689362B2 JP3359339A JP35933991A JP2689362B2 JP 2689362 B2 JP2689362 B2 JP 2689362B2 JP 3359339 A JP3359339 A JP 3359339A JP 35933991 A JP35933991 A JP 35933991A JP 2689362 B2 JP2689362 B2 JP 2689362B2
Authority
JP
Japan
Prior art keywords
air
fuel ratio
cylinder
internal combustion
combustion engine
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
JP3359339A
Other languages
Japanese (ja)
Other versions
JPH05180059A (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.)
Honda Motor Co Ltd
Original Assignee
Honda 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 Honda Motor Co Ltd filed Critical Honda Motor Co Ltd
Priority to JP3359339A priority Critical patent/JP2689362B2/en
Priority to EP92311841A priority patent/EP0553570B1/en
Priority to DE69225212T priority patent/DE69225212T2/en
Publication of JPH05180059A publication Critical patent/JPH05180059A/en
Priority to US08/282,104 priority patent/US5524598A/en
Application granted granted Critical
Publication of JP2689362B2 publication Critical patent/JP2689362B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Description

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

【0001】[0001]

【産業上の利用分野】この発明は内燃機関の空燃比検出
方法に関し、より具体的には空燃比センサの応答遅れを
推定して空燃比をリアルタイムに求めることができる様
にした内燃機関の空燃比検出方法に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a method of detecting an air-fuel ratio of an internal combustion engine, and more specifically, to an air-fuel ratio of an internal combustion engine capable of obtaining a real-time air-fuel ratio by estimating a response delay of an air-fuel ratio sensor. The present invention relates to a fuel ratio detection method.

【0002】[0002]

【従来の技術】内燃機関の排気系に酸素濃度検出素子か
らなる空燃比センサを設けて入力した燃料の空燃比を検
出することは良く行われており、その一例として特開昭
59─101562号公報記載の技術を挙げることがで
きる。この技術においては検出精度を上げるため、基準
タイミング(第1気筒TDC)から各気筒の排気ガスが
空燃比センサに到達するまでの遅れ時間を運転状態に応
じて予め求めておき、それに基づいて気筒ごとに空燃比
を検出して目標値にフィードバック制御している。
2. Description of the Related Art It is often practiced to provide an air-fuel ratio sensor consisting of an oxygen concentration detecting element in the exhaust system of an internal combustion engine to detect the air-fuel ratio of the input fuel, and one example thereof is JP-A-59-101562. The technology described in the publication can be mentioned. In this technique, in order to improve the detection accuracy, a delay time from the reference timing (first cylinder TDC) until the exhaust gas of each cylinder reaches the air-fuel ratio sensor is obtained in advance according to the operating state, and the cylinder is based on that. The air-fuel ratio is detected every time and feedback control is performed to the target value.

【0003】[0003]

【発明が解決しようとする課題】ところで、前記した酸
素濃度検出素子からなる空燃比センサは、排気ガスに接
触する素子に生ずる化学変化を起電力に変換して検出す
るため、排気ガスがセンサに到達した後も検出値として
現れるまでに遅れがある。従って、その遅れを解明しな
い限り、正確に燃焼した混合気の空燃比を求めることが
できず、よって目標値にフィードバック制御するときな
ども高精度で収束性に優れた空燃比制御を実現すること
ができない。
By the way, in the air-fuel ratio sensor comprising the oxygen concentration detecting element described above, the chemical change occurring in the element contacting the exhaust gas is converted into an electromotive force to be detected, so that the exhaust gas acts as a sensor. There is a delay in appearing as a detected value even after reaching. Therefore, unless the delay is clarified, the air-fuel ratio of the burned air-fuel mixture cannot be accurately obtained, and therefore, it is possible to realize highly accurate and highly convergent air-fuel ratio control even when performing feedback control to the target value. I can't.

【0004】従って、本発明の目的は上記した欠点を解
消し、空燃比センサの検出応答遅れを的確に推定して燃
焼した混合気の最新の値を求めることができ、よってフ
ィードバック制御を行うときなども精度と収束性に優れ
た空燃比制御を実現できる様にした内燃機関の空燃比検
出方法を提供することにある。
Therefore, the object of the present invention is to solve the above-mentioned drawbacks and to accurately estimate the detection response delay of the air-fuel ratio sensor to obtain the latest value of the burned air-fuel mixture. It is another object of the present invention to provide an air-fuel ratio detection method for an internal combustion engine, which can realize air-fuel ratio control with excellent accuracy and convergence.

【0005】[0005]

【課題を解決するための手段】上記の目的を解決するた
めに本発明に係る内燃機関の空燃比検出方法は例えば請
求項1項に示すように、内燃機関の排気系に配置される
空燃比センサの検出値から各気筒空燃比を検出する
であって、前記センサの応答遅れを1次遅れ系で擬似
的にモデル化し、その挙動を記述する状態方程式を設定
、前記状態方程式を周期ΔTで離散化して伝達関数
(l−αハット)/(Z−αハット)を求め、前記伝達
関数の係数αハットを所定機関回転数ごとに相違させつ
つ前記伝達関数の逆伝達関数(Z−αハット)/(l−
αハット)を求め、求めた逆伝達関数を前記空燃比セン
出力に乗じて内燃機関に入力される混合気の空燃比
の推定値を求め、求めた空燃比の推定値を用いて前記内
燃機関の排気系の挙動を記述するモデルを設定すると共
に、その内部状態を観測するオブザーバを設定し、およ
び前記オブザーバの出力から各気筒の空燃比を検出する
く構成した。
Means for Solving the Problems] air-fuel ratio detection method for an internal combustion engine according to the present invention in order to solve the above object, as shown in item 1 example claims, the air-fuel ratio is disposed in an exhaust system of an internal combustion engine also detects the air-fuel ratio of each cylinder from the detected values of the sensors meet the <br/>, and pseudo-model the response delay of the sensor in first-order system, setting a state equation describing the behavior
Then , the above equation of state is discretized with a period ΔT to obtain a transfer function.
(1−α hat) / (Z−α hat) is obtained, and the coefficient α hat of the transfer function is varied for each predetermined engine speed.
One inverse transfer function (Z-alpha hat) of the transfer function / (l-
α hat) , the obtained inverse transfer function is multiplied by the output of the air-fuel ratio sensor to obtain an estimated value of the air-fuel ratio of the air-fuel mixture input to the internal combustion engine, and the obtained estimated value of the air-fuel ratio is obtained . Within the above
It is common to set up a model that describes the behavior of the exhaust system of a combustion engine.
, Set an observer to observe its internal state, and
And the air-fuel ratio of each cylinder is detected from the output of the observer.
rather it was constructed.

【0006】[0006]

【作用】センサの検出応答遅れを推定して燃焼した混合
気の空燃比を気筒ごとにリアルタイムかつ正確にめる
ことができ、気筒ごとのフィードバック制御を行うとき
なども制御の応答性が向上し、よって精度と収束性に優
れた空燃比制御を実現することが可能となる。また、気
筒ごとに空燃比センサを配置する必要がなく、単一の空
燃比センサの出力に基づいて各気筒の空燃比を検出ある
いは推定することができる。
[Action] estimates the detection response delay of the sensor air-fuel ratio of the mixture burned could determined Mel be in real time and accurate to each cylinder, the response of the well control, such as when performing feedback control for each cylinder Therefore, it becomes possible to realize the air-fuel ratio control with improved accuracy and convergence. Also damn
There is no need to arrange an air-fuel ratio sensor for each cylinder,
The air-fuel ratio of each cylinder is detected based on the output of the fuel ratio sensor.
It can be estimated.

【0007】[0007]

【実施例】以下、添付図面に即して本発明の実施例を説
明する。
Embodiments of the present invention will be described below with reference to the accompanying drawings.

【0008】図1は本発明にかかる方法を実現するため
の内燃機関の空燃比検出・制御装置を全体的に示す概略
図である。図において、符号10は4気筒の内燃機関を
示しており、吸気路12の先端に配置されたエアクリー
ナ14から導入された吸気は、スロットル弁16でその
流量を調節されつつインテークマニホルド18を経て第
1〜第4気筒に流入される。各気筒の吸気弁(図示せ
ず)の付近にはインジェクタ20が設けられて燃料を噴
射する。噴射され吸気と一体となった混合気は各気筒内
で図示しない点火プラグで点火されて燃焼してピストン
(図示せず)を駆動する。燃焼後の排気ガスは排気弁
(図示せず)を介してエキゾーストマニホルド22に排
出され、エキゾーストパイプ24を経て三元触媒コンバ
ータ26で浄化されつつ機関外に排出される。また吸気
路12にはスロットル弁16を配置した位置付近にそれ
をバイパスするバイパス路28が設けられる。
FIG. 1 is an overall schematic view of an air-fuel ratio detection / control apparatus for an internal combustion engine for implementing the method according to the present invention. In the figure, reference numeral 10 denotes a four-cylinder internal combustion engine, and intake air introduced from an air cleaner 14 disposed at the tip of an intake passage 12 passes through an intake manifold 18 while its flow rate is adjusted by a throttle valve 16. The fuel flows into the first to fourth cylinders. An injector 20 is provided near an intake valve (not shown) of each cylinder to inject fuel. The air-fuel mixture that has been injected and integrated with the intake air is ignited by a spark plug (not shown) in each cylinder, burns, and drives a piston (not shown). The exhaust gas after the combustion is discharged to an exhaust manifold 22 through an exhaust valve (not shown), and is discharged outside the engine through an exhaust pipe 24 while being purified by a three-way catalytic converter 26. In addition, a bypass passage 28 is provided in the intake passage 12 near the position where the throttle valve 16 is disposed to bypass the throttle valve 16.

【0009】また内燃機関10のディストリビュータ
(図示せず)内にはピストン(図示せず)のクランク角
度位置を検出するクランク角センサ34が設けられると
共に、スロットル弁16の開度を検出するスロットル開
度センサ36、スロットル弁16下流の吸気圧力を絶対
圧力で検出する絶対圧センサ38も設けられる。更に、
排気系においてエキゾーストマニホルド22の下流側で
三元触媒コンバータ26の上流側には酸素濃度検出素子
からなる空燃比センサ40が設けられ、排気ガスの空燃
比を検出する。これらセンサ34などの出力は、制御ユ
ニット42に送られる。
A crank angle sensor 34 for detecting the crank angle position of a piston (not shown) is provided in a distributor (not shown) of the internal combustion engine 10, and a throttle opening for detecting the opening of the throttle valve 16 is provided. A degree sensor 36 and an absolute pressure sensor 38 for detecting the intake pressure downstream of the throttle valve 16 by absolute pressure are also provided. Furthermore,
In the exhaust system, an air-fuel ratio sensor 40 composed of an oxygen concentration detecting element is provided downstream of the exhaust manifold 22 and upstream of the three-way catalytic converter 26, and detects an air-fuel ratio of exhaust gas. Outputs of these sensors 34 and the like are sent to the control unit 42.

【0010】図2は制御ユニット42の詳細を示すブロ
ック図である。空燃比センサ40の出力は検出回路46
に入力され、そこで適宜な線型化処理が行われてリーン
からリッチにわたる広い範囲において排気ガス中の酸素
濃度に比例したリニアな特性からなる空燃比(A/F)
が検出される。尚、その詳細は先に本出願人が提案した
出願(特願平3−169456号)に述べてあるので、
これ以上の説明は省略する。また以下の説明においてこ
のセンサを「LAFセンサ」(リニア・エーバイエフ・
センサ)と称する。検出回路46の出力はA/D変換回
路48を介してCPU50,ROM52,RAM54か
らなるマイクロ・コンピュータ内に取り込まれ、RAM
54に格納される。同様にスロットル開度センサ36な
どのアナログ出力はレベル変換回路56、マルチプレク
サ58及び第2のA/D変換回路60を介して、またク
ランク角センサ34の出力は波形整形回路62で波形整
形された後、カウンタ64で出力値かカウントされ、カ
ウント値はマイクロ・コンピュータ内に入力される。マ
イクロ・コンピュータにおいてCPU50はROM52
に格納された命令に従って検出値から空燃比のフィード
バック制御値を演算し、駆動回路66を介して各気筒の
インジェクタ20を駆動すると共に、第2の駆動回路6
8を介して電磁弁70を駆動し、図1に示したバイ
路28を通る2次空気量を制御する。
FIG. 2 is a block diagram showing details of the control unit 42. The output of the air-fuel ratio sensor 40 is
The air-fuel ratio (A / F) having a linear characteristic proportional to the oxygen concentration in the exhaust gas in a wide range from lean to rich, where appropriate linearization processing is performed.
Is detected. The details are described in the application proposed by the present applicant (Japanese Patent Application No. 3-169456).
Further description is omitted. In the following description, this sensor is referred to as a “LAF sensor” (Linear AV
Sensor). The output of the detection circuit 46 is taken into a microcomputer composed of a CPU 50, a ROM 52, and a RAM 54 via an A / D conversion circuit 48, and is stored in a RAM.
54. Analog output level conversion circuit 56, such as well the throttle opening sensor 36, via the multiplexer 58 and the second A / D converter circuit 60, also the output of the crank angle sensor 34 is shaped by a waveform shaper 62 After that, the counter 64 counts the output value, and the count value is input into the microcomputer. In the microcomputer, the CPU 50 is a ROM 52
Calculates the feedback control value of the air-fuel ratio from the detected value in accordance with the command stored in the second driving circuit 6 and drives the injector 20 of each cylinder via the driving circuit 66.
8 by driving the solenoid valve 70 via, for controlling the amount of secondary air passing through the by-path passage 28 shown in FIG.

【0011】続いて、この制御装置の動作を説明する。
その動作は結局のところ空燃比を検出して目標値との偏
差を求め、その偏差を解消する様に燃料供給量(噴射
量)を制御することにつきるが、この発明の要旨は、空
燃比センサの応答遅れを求めて空燃比をリアルタイムに
推定することにあることにあるので、以下その点に焦点
をおいて説明する。また併せて検出した空燃比に基づい
て空燃比を制御する例を説明する。
Next, the operation of this control device will be described.
After all, the operation is to detect the air-fuel ratio, find the deviation from the target value, and control the fuel supply amount (injection amount) so as to eliminate the deviation. This is because the air-fuel ratio is to be estimated in real time by obtaining the response delay of the sensor, so the following description focuses on that point. An example of controlling the air-fuel ratio based on the detected air-fuel ratio will also be described.

【0012】図3は単気筒の内燃機関を例にとり、吸気
量を一定として供給燃料量をステップ状に変化させた場
合の空燃比センサの応答の実測データを示す(図中「実
測値」)。図示の如く、空燃比をステップ状に変化させ
た場合、LAFセンサ出力の実測値は入力値に対して遅
れを持つが、この遅れはセンサの化学反応に起因するの
で、正確に解析することは困難である。そこで本発明者
達はとりあえずこの遅れを1次遅れ系と擬似的にモデル
化し、図4に示す如きモデルを作成した。ここでLA
F:LAFセンサ出力、A/F:入力A/F、とする
と、その状態方程式は下記の数1で示すことができる。
FIG. 3 shows an example of a single-cylinder internal combustion engine, and shows measured data of the response of the air-fuel ratio sensor when the supplied fuel amount is changed stepwise while the intake amount is constant ("measured value" in the figure). . As shown in the figure, when the air-fuel ratio is changed stepwise, the measured value of the LAF sensor output has a delay with respect to the input value, but this delay is caused by a chemical reaction of the sensor, so that accurate analysis cannot be performed. Have difficulty. Therefore, the present inventors have tentatively modeled this delay as a first-order delay system and created a model as shown in FIG. Here LA
Assuming that F: LAF sensor output and A / F: input A / F, the state equation can be expressed by the following equation 1.

【0013】[0013]

【数1】 (Equation 1)

【0014】これを周期ΔTで離散化すると、数2で示
す様になる。図5は数2をブロック線図で表したもので
ある。
When this is discretized with the period ΔT, it becomes as shown by the equation 2. FIG. 5 is a block diagram of Equation (2).

【0015】[0015]

【数2】 (Equation 2)

【0016】従って、数2を用いることによってセンサ
出力より真の空燃比を求めることができる。即ち、数2
を変形すれば数3に示す様になるので、時刻kのときの
値から時刻k−1のときの値を数4の様に逆算すること
ができる。
Therefore, by using the equation 2, the true air-fuel ratio can be obtained from the sensor output. That is, Equation 2
Is modified as shown in equation (3), the value at time k-1 can be inversely calculated from the value at time k as in equation (4).

【0017】[0017]

【数3】 (Equation 3)

【0018】[0018]

【数4】 (Equation 4)

【0019】具体的には数2をZ変換を用いて伝達関数
で示せば数5の如くになるので、その逆伝達関数を今回
のセンサ出力LAFに乗じることによって前回の空燃比
をリアルタイムに推定することができる。図6にそのリ
アルタイムのA/F推定器のブロック線図を示す。尚、
前述したとおりLAFセンサの応答遅れは化学反応に起
因するもので正確に解析することは困難であるが、機関
回転数との間に相関関係のあることが確認できた。した
がって伝達関数の係数は、適宜設定する所定の機関回転
数ごとに相違させることとした。よって、所定機関回転
数ごとにA/F推定器、すなわち逆伝達関数の係数を持
ち換えることにより、推定するA/F値の精度をより向
上させることができる。
Specifically, if the expression 2 is expressed by a transfer function using the Z-transform, it becomes as shown in the expression 5. Therefore, the previous air-fuel ratio is estimated in real time by multiplying the inverse transfer function by the sensor output LAF of this time. can do. FIG. 6 shows a block diagram of the real-time A / F estimator. still,
As described above, the response delay of the LAF sensor is caused by a chemical reaction and is difficult to analyze accurately, but it was confirmed that there is a correlation with the engine speed. Therefore, the coefficient of the transfer function is set to be different for each predetermined engine speed set as appropriate. Therefore, by changing the A / F estimator, that is, the coefficient of the inverse transfer function for each predetermined engine speed, the accuracy of the estimated A / F value can be further improved.

【0020】[0020]

【数5】 (Equation 5)

【0021】上記についてのシミュレーション結果を図
3(図中「シミュレーション」)、図7に示す。前述し
たとおり、図3に「実測値」とあるのはステップ状の空
燃比入力を与えた場合のセンサ出力の実測値であるが、
ここで実測値とシミュレーション結果(図5モデルにス
テップ状空燃比を入力して得た出力)とがほぼ対応して
いるのが見てとれよう。以上からセンサ応答遅れを1次
遅れと擬似的にモデル化したことの正しさが検証された
と言うことができよう。図7はセンサ出力の実測値に逆
伝達関数を乗じ、真の空燃比を推定する場合を示す。同
図において、例えば時刻Taの時の真の空燃比は、1
2.5ではなく13.2であると推定することができる
のである。尚、真の空燃比の推定値に若干の上下動があ
らわれているのは、センサ出力の実測値に細かなばらつ
きがあるためである。
Simulation results for the above are shown in FIG. 3 (“simulation” in the figure) and FIG. 7. As described above, “actual measurement value” in FIG. 3 is an actual measurement value of the sensor output when a step-like air-fuel ratio input is given.
Here, it can be seen that the measured values and the simulation results (outputs obtained by inputting the step-like air-fuel ratio to the model in FIG. 5) substantially correspond to each other. From the above, it can be said that the correctness of pseudo modeling of the sensor response delay as the first-order delay has been verified. FIG. 7 shows a case where a true air-fuel ratio is estimated by multiplying an actually measured value of the sensor output by an inverse transfer function. In the figure, for example, the true air-fuel ratio at time Ta is 1
It can be estimated to be 13.2 instead of 2.5. The reason why the estimated value of the true air-fuel ratio slightly fluctuates is that there is a small variation in the actually measured value of the sensor output.

【0022】続いて、上記の如く求めた真の空燃比に基
づいて空燃比を制御する場合について説明する。
Next, the case where the air-fuel ratio is controlled based on the true air-fuel ratio obtained as described above will be described.

【0023】先に述べた如く、多気筒内燃機関に1個の
空燃比センサを配置した場合、その出力は、排気系の集
合部において各気筒の検出値が混ざり合った値を示し、
気筒ごとの真の検出値を求めることが困難である。その
ため、各気筒のA/Fを個別に目標値に制御することが
できず、ある気筒はリーンであったり、他の気筒はリッ
チであったりする場合があって、エミッション悪化の原
因となる。それを解決するためには気筒ごとにセンサを
配置すれば良いが、それではコストの上昇を招く。そこ
で本発明者達はセンサ応答遅れを1次遅れでモデル化で
きたことによって、以下の手法から排気系の集合部に配
置した1個の空燃比センサで多気筒、実施例の場合は4
気筒の内燃機関の空燃比を気筒ごとに正確に検出できる
様にした。以下、それについて説明する。
As described above, when one air-fuel ratio sensor is arranged in the multi-cylinder internal combustion engine, the output thereof shows a value obtained by mixing the detection values of the respective cylinders at the collecting portion of the exhaust system,
It is difficult to obtain a true detection value for each cylinder. Therefore, the A / F of each cylinder cannot be individually controlled to the target value, and a certain cylinder may be lean or another cylinder may be rich, which may cause deterioration of emission. In order to solve this, a sensor may be arranged for each cylinder, but this causes an increase in cost. Therefore, the inventors of the present invention have been able to model the sensor response delay as a first-order delay. Therefore, one air-fuel ratio sensor arranged in the collecting part of the exhaust system is used for multiple cylinders in the following method, and in the case of the embodiment, 4
The air-fuel ratio of the cylinder internal combustion engine can be detected accurately for each cylinder. This will be described below.

【0024】先ず、図8に示す様に内燃機関の排気系を
モデル化した(以下、このモデルを「エキマニモデル」
と称する)。尚、このエキマニモデルでは、離散系のサ
ンプル時間をTDC周期(機関回転数が1500rpm
のとき0.02sec)とした。またこのエキマニモデ
ルの中ではF(燃料)を制御量としたため、空燃比をF
/Aとした。
First, an exhaust system of an internal combustion engine was modeled as shown in FIG. 8 (hereinafter, this model is referred to as an "exhaust manifold model").
Called)). In this exhaust manifold model, the sampling time of the discrete system is set to the TDC cycle (when the engine speed is 1500 rpm).
0.02 sec). In this exhaust manifold model, since F (fuel) is a controlled variable, the air-fuel ratio is
/ A.

【0025】ここで発明者達は、排気系の集合部の空燃
比(A/F)は、各気筒の空燃比の時間的な寄与度を考
慮した加重平均値であると考えた。そうであれば、時刻
kのときの集合部の空燃比は、数6の様に表すことがで
きる。
The inventors have considered that the air-fuel ratio (A / F) of the collecting portion of the exhaust system is a weighted average value in consideration of the temporal contribution of the air-fuel ratio of each cylinder. If so, the air-fuel ratio of the collecting part at the time k can be expressed as in Equation 6.

【0026】[0026]

【数6】 (Equation 6)

【0027】即ち、集合部の空燃比は、気筒ごとの過去
の燃焼履歴に重みC(例えば直近に燃焼した気筒は40
%、その前が30%...など)を乗じたものの合算で
表すことができる。ここで、集合部における各気筒の排
気の混合状態は機関の運転状態により異なってくる。即
ち、例えば機関の低回転域ではTDC周期が長いので、
各気筒からの排気が混合する度合いは高回転域と比較し
て低い。また、高負荷のときは基本的に背圧も大きく、
排気の排出圧力が大きくなるので、各気筒からの排気が
混合する度合いは低負荷のときと比較して低い。この様
に各気筒の排気の混合する度合いが低い場合には、直近
に燃焼した気筒の重みを大きくとる必要がある。従っ
て、重みCは機関の運転状態によって持ち換えることと
する。具体的には重みCを機関回転数と負荷とをパラメ
ータとして適宜設定してマップに用意しておき、それを
検索することにより行う。尚、上記で#nは気筒番号を
示し、また気筒の燃焼(点火)順序は、1,3,4,2
とする。またここで空燃比[F/A]は先に数5で求め
た応答遅れを補正した真の値を意味する。
That is, the air-fuel ratio of the collecting portion is weighted by C in the past combustion history for each cylinder (for example, the most recently burned cylinder is 40
%, Before that 30%. . . , Etc.) can be expressed as the sum of the products. Here, the mixed state of the exhaust gas of each cylinder in the collecting part differs depending on the operating state of the engine. That is, for example, the TDC cycle is long in a low engine speed range,
The degree of mixing of the exhaust from each cylinder is lower than in the high rotation range. Also, when the load is high, the back pressure is basically large,
Since the exhaust pressure of the exhaust gas increases, the degree of mixing of the exhaust gas from each cylinder is lower than when the load is low. In the case where the degree of mixing of the exhaust gases of the cylinders is low, it is necessary to increase the weight of the cylinder that has burned most recently. Therefore, the weight C is changed depending on the operating state of the engine. Specifically, the weight C is set by appropriately setting the engine speed and the load as parameters, prepared in a map, and searched for. In the above, #n indicates the cylinder number, and the combustion (ignition) order of the cylinders is 1, 3, 4, 2
And Here, the air-fuel ratio [F / A] means a true value obtained by correcting the response delay previously obtained by Expression 5.

【0028】上記を前提とすると、エキマニモデルの状
態方程式は数7の様になる。
Based on the above, the equation of state of the exhaust manifold model is as shown in equation 7.

【0029】[0029]

【数7】 (Equation 7)

【0030】また集合部の空燃比をy(k)とおくと、
出力方程式は数8の様に表すことができる。
If the air-fuel ratio of the collecting portion is y (k),
The output equation can be expressed as in Equation 8.

【0031】[0031]

【数8】 ここで、 c:0.25379,c:0.10121,c
0.46111,c:0.18389とする。
(Equation 8) Here, c 1 : 0.25379, c 2 : 0.101021, c 3 :
0.46111, c 4: to 0.18389.

【0032】上記において、u(k)は観測不可能のた
め、この状態方程式からオブザーバを設計してもx
(k)は観測することができない。そこで4TDC前
(即ち、同一気筒)の空燃比F/Aは急激に変化しない
定常運転状態にあると仮定してx(k+1)=x(k−
3)とすると、数9の様になる。
In the above, since u (k) cannot be observed, even if an observer is designed from this equation of state, x (x)
(K) cannot be observed. Therefore, assuming that the air-fuel ratio F / A before 4TDC (that is, the same cylinder) is in a steady operation state in which it does not suddenly change, x (k + 1) = x (k−
If it is 3), it becomes like Equation 9.

【0033】[0033]

【数9】 (Equation 9)

【0034】ここで、上記の如く求めたエキマニモデル
についてシミュレーション結果を示す。図9は4気筒内
燃機関について3気筒の空燃比を14.7にし、1気筒
だけ12.0にして燃料を供給した場合を示す。図10
はそのときの集合部(即ち、図1のエキゾーストマニホ
ルドパイプ24に空燃比センサ40を配置した位置)の
空燃比(A/F)を上記エキマニモデルで求めたものを
示す。図10においてはステップ状の出力が得られてい
るが、ここで更にLAFセンサの応答遅れを考慮する
と、センサ出力は図11に「シミュレーション」と示す
様になまされた波形となる。図中「実測値」は同じ場合
のLAFセンサ出力の実測値であるが、これと比較し、
上記エキマニモデルが多気筒内燃機関の排気系を良くモ
デル化していることが検証できたといえよう。
Here, simulation results of the exhaust manifold model obtained as described above will be shown. FIG. 9 shows a case where the air-fuel ratio of three cylinders is set to 14.7 and the fuel is supplied to one cylinder at 12.0 for a four-cylinder internal combustion engine. FIG.
Shows the air-fuel ratio (A / F) of the collecting portion (that is, the position where the air-fuel ratio sensor 40 is disposed on the exhaust manifold pipe 24 in FIG. 1) obtained by the above-mentioned exhaust manifold model. Although a step-like output is obtained in FIG. 10, if the response delay of the LAF sensor is further considered, the sensor output has a waveform simulated as “simulation” in FIG. In the figure, “actual measurement value” is an actual measurement value of the output of the LAF sensor in the same case.
It can be verified that the exhaust manifold model described above models the exhaust system of the multi-cylinder internal combustion engine well.

【0035】よって、数10で示される状態方程式と出
力方程式にてx(k)を観察する定常のカルマンフィル
タの問題に帰着する。その荷重行列Q,Rを数11の様
においてリカッチの方程式を解くと、ゲイン行列Kは数
12の様になる。
Therefore, the problem of a stationary Kalman filter for observing x (k) in the state equation and the output equation shown in the equation 10 is reduced. When the Riccati equation is solved by using the weighting matrices Q and R as shown in the equation 11, the gain matrix K becomes as shown in the equation 12.

【0036】[0036]

【数10】 (Equation 10)

【0037】[0037]

【数11】 [Equation 11]

【0038】[0038]

【数12】 (Equation 12)

【0039】これよりA−KCを求めると、数13の様
になる。
When A-KC is calculated from this, the result is as shown in the equation 13.

【0040】[0040]

【数13】 (Equation 13)

【0041】一般的なオブザーバの構成は図12に示さ
れる様になるが、今回のモデルでは入力u(k)がない
ので、図13で示す様にy(k)のみを入力とする構成
となり、これを数式で表すと数14の様になる。
The structure of a general observer is as shown in FIG. 12, but since there is no input u (k) in this model, the structure is such that only y (k) is input as shown in FIG. When this is expressed by a mathematical expression, it becomes as shown in Expression 14.

【0042】[0042]

【数14】 [Equation 14]

【0043】ここでy(k)を入力とするオブザーバ、
即ちカルマンフィルタのシステム行列は数15の様に表
される。
Here, an observer whose input is y (k),
That is, the system matrix of the Kalman filter is expressed as in Equation 15.

【0044】[0044]

【数15】 (Equation 15)

【0045】今回のモデルで、リカッチ方程式の荷重配
分Rの要素:Qの要素=1:1のとき、カルマンフィル
タのシステム行列Sは、数16で与えられる。
In the model this time, when the element of the weight distribution R of the Riccati equation: the element of Q = 1: 1, the system matrix S of the Kalman filter is given by the equation 16.

【0046】[0046]

【数16】 (Equation 16)

【0047】続いて、シミュレーション上で気筒別空燃
比の波形を正確に作成し、それを前記エキマニモデルに
入力し、集合部空燃比を得る。それをオブザーバに入力
し、気筒別空燃比が推定されていることを検証する。ま
た荷重行列と推定値の傾向を考察する。
Then, the waveform of the cylinder-by-cylinder air-fuel ratio is accurately created on the simulation, and the waveform is input to the exhaust manifold model to obtain the collective portion air-fuel ratio. This is input to the observer, and it is verified that the cylinder-by-cylinder air-fuel ratio is estimated. We also consider the tendency of the weight matrix and the estimated value.

【0048】今回のモデルにおいては、数17に示す様
であるので、荷重行列Qは、要素が全て同じ対角行列と
なる。
In the model this time, as shown in Equation 17, the weight matrix Q is a diagonal matrix in which all the elements are the same.

【0049】[0049]

【数17】 [Equation 17]

【0050】従って、考察すべきはQとRの要素の比で
ある。QとRの要素の比を変えて求めたゲインを表1に
示す。また、それを用いて構成したオバザーバとエキマ
ニモデルを組み合わせたシミュレーションモデルを図1
4に示す。更に、このモデルを用いて気筒別空燃比を1
2.0,14.7,14.7,14.7の理想入力とし
て計算したものを図15に、またそのときのオブザーバ
の推定誤差を表2に示す。更に、空燃比を12.0±
0.2,14.7±0.2,14.7±0.2,14.
7±0.2としてそれぞれ独立に変動させて(仮想ノイ
ズ)計算したものを図16に、そのときのオブザーバの
推定誤差を表3に示す。尚、図15,16において
(a)から(e)は、共に、(a)各気筒A/F(エキ
マニモデル入力)、(b)集合部A/F(エキマニモデ
ル出力)、(c)Qの要素:Rの要素=1:10のとき
のオブザーバ出力(入力は(b)に示す)、(d)Qの
要素:Rの要素=1:1のときのオブザーバ出力(入力
は(b)に示す)、(e)Qの要素:Rの要素=10:
1のときのオブザーバ出力(入力は(b)に示す)であ
る。
Therefore, what is to be considered is the ratio of the elements of Q and R. Table 1 shows the gains obtained by changing the ratio of the Q and R elements. Fig. 1 shows a simulation model combining an observer and an exhaust manifold model constructed using the model.
It is shown in FIG. Furthermore, using this model, the air-fuel ratio for each cylinder is set to 1
FIG. 15 shows the results calculated as ideal inputs of 2.0, 14.7, 14.7, and 14.7, and Table 2 shows observer estimation errors at that time. Further, the air-fuel ratio is set to 12.0 ±
0.2, 14.7 ± 0.2, 14.7 ± 0.2, 14.
FIG. 16 shows the results of calculation (virtual noise) independently varied as 7 ± 0.2, and Table 3 shows the estimated error of the observer at that time. In FIGS. 15 and 16, (a) to (e) show (a) each cylinder A / F (exhaust model input), (b) collecting section A / F (exit model output), and (c) Q Element: observer output when R element = 1: 10 (input is shown in (b)), (d) Observer output when Q element: R element = 1: 1 (input is (b) (E) Elements of Q: Elements of R = 10:
1 is the observer output (input is shown in (b)).

【0051】[0051]

【表1】 [Table 1]

【0052】[0052]

【表2】 [Table 2]

【0053】[0053]

【表3】 [Table 3]

【0054】図15の様に各気筒空燃比を一定としたと
きはQの重みが大きいほど収束が速いことが分かる。但
し、Q/Rを10以上にしても殆ど収束性は変わらなか
った。図16において時系列に推定偏差(各気筒空燃比
−推定空燃比)を図示すると、図17の様になり、オブ
ザーバ収束後はQの要素:Rの要素=10:1と1:1
とでそれほど差がないことから、耐外乱性を考えると、
Qの要素:Rの要素=1:1の方が良いと言える。この
様に、集合部空燃比の入力に対し、カルマンフィルタの
理論を用いたオブザーバは、集合部における気筒ごとの
空燃比を精度良く推定することを可能にする。尚、荷重
行列はQ/R=1〜10が最良であったが、実データを
用いた応答状況から決める必要があると思われる。
As shown in FIG. 15, when the air-fuel ratio of each cylinder is made constant, the larger the weight of Q, the faster the convergence. However, even when Q / R was 10 or more, the convergence was hardly changed. In FIG. 16, the estimated deviation (each cylinder air-fuel ratio−estimated air-fuel ratio) is illustrated in time series as shown in FIG. 17, and after the convergence of the observer, the Q element: the R element = 10: 1 and 1: 1.
Since there is not much difference between and, considering the disturbance resistance,
It can be said that the factor of Q: the factor of R = 1: 1 is better. As described above, the observer using the Kalman filter theory with respect to the input of the air-fuel ratio of the collecting portion enables the air-fuel ratio of each cylinder in the collecting portion to be accurately estimated. Although the weight matrix was best when Q / R = 1 to 10, it seems necessary to determine it from the response situation using actual data.

【0055】続いて、実測データを先に示したLAFセ
ンサの逆伝達関数に入力して得られる実集合部空燃比デ
ータを前記オブザーバに入力し、気筒別空燃比を推定し
た結果を図18に示す。同図において(a)LAFセン
サ出力、(b)LAFセンサ逆伝達関数出力(入力は
(a)に示す)、(c)Qの要素:Rの要素=1:10
のときのオブザーバ出力(入力は(b)に示す)、
(d)Qの要素:Rの要素=1:1のときのオブザーバ
出力(入力は(b)に示す)、(e)Qの要素:Rの要
素=10:1のときのオブザーバ出力(入力は(b)に
示す)である。ここでLAFセンサ出力の測定条件は、
機関回転数=1500rpm、吸気圧力=−281.9
mmHg、A/F=12.0(#2),14.7(#
1,#3,#4)とした。また、実際の入力A/Fの真
値はわからないので、シミュレーションではおよその値
として、〔12.0/14.7/14.7/14.7〕
を用いた。同図から明らかな様に、オブザーバ出力は4
TDC周期で変化しており、入力A/Fをほぼ推定して
いる。またカルマンフィルタを用いたことにより、荷重
行列のセッティングによって2〜8周期で収束可能であ
ることが確認された。
Subsequently, the actual assembly air-fuel ratio data obtained by inputting the actually measured data to the inverse transfer function of the LAF sensor shown above is input to the observer, and the result of estimating the cylinder-by-cylinder air-fuel ratio is shown in FIG. Show. In the figure, (a) LAF sensor output, (b) LAF sensor inverse transfer function output (input is shown in (a)), (c) element of Q: element of R = 1: 10
Observer output at the time of (input is shown in (b)),
(D) Observer output when Q element: R element = 1: 1: 1 (input is shown in (b)), (e) Observer output when Q element: R element = 10: 1 (input (Shown in (b)). Here, the measurement conditions of the LAF sensor output are as follows:
Engine speed = 1500 rpm, intake pressure = -281.9
mmHg, A / F = 12.0 (# 2), 14.7 (#
1, # 3, # 4). Also, since the true value of the actual input A / F is not known, the approximate value is [12.0 / 14.7 / 14.7 / 14.7] in the simulation.
Was used. As is clear from the figure, the observer output is 4
It changes in the TDC cycle, and the input A / F is almost estimated. It was also confirmed that the use of the Kalman filter enabled convergence in 2 to 8 periods by setting the weight matrix.

【0056】続いて、上記の如く推定して得られた気筒
別空燃比を用いて、空燃比を目標値に制御する。図19
はPID手法を用いた公知の制御例を示すブロック線図
である。フィードバックが乗算項を通じてなされる点が
通常のPID制御と異なるが、この手法自体は公知なも
のであって、図示の如く、入力Ti(噴射時間)によっ
て生じる実際の空燃比の目標値に対する偏差(1−1/
λ) を気筒ごとに求め、それに応じたゲインKLAFを
乗じて目標値にフィードバック制御すれば足る。但し、
制御手法自体は公知であっても、前記の如く各気筒の空
燃比を正確に分離抽出することができたので、それぞれ
の気筒の空燃比を目標値に精度良く制御することができ
る。
Next, the air-fuel ratio is controlled to a target value by using the cylinder-by-cylinder air-fuel ratio estimated as described above. FIG.
FIG. 6 is a block diagram showing a known control example using the PID method. Although the feedback is performed through the multiplication term, which is different from the normal PID control, this method is well known in the art, and as shown in the figure, the deviation of the actual air-fuel ratio caused by the input Ti (injection time) from the target value ( 1-1 /
It is sufficient to obtain λ) for each cylinder, multiply the gain KLAF accordingly, and perform feedback control to the target value. However,
Although the control method itself is publicly known, the air-fuel ratio of each cylinder can be accurately separated and extracted as described above, so that the air-fuel ratio of each cylinder can be accurately controlled to the target value.

【0057】上記した実施例において、空燃比センサの
検出応答遅れを1次遅れで擬似的にモデル化することが
できたので、空燃比センサの入力空燃比を正確に求める
ことができた。また、それによって排気系に1個のみ設
けた空燃比センサの集合部出力から気筒毎の空燃比を抽
出することを可能としたので、それによって気筒ごとに
空燃比を目標値に精度良く制御することができた。
In the above-described embodiment, the detection response delay of the air-fuel ratio sensor can be modeled as a primary delay, so that the input air-fuel ratio of the air-fuel ratio sensor can be accurately obtained. Further, it is possible to extract the air-fuel ratio of each cylinder from the output of the collecting portion of the air-fuel ratio sensor provided only in the exhaust system, so that the air-fuel ratio of each cylinder is accurately controlled to the target value. I was able to.

【0058】尚、上記実施例において、空燃比センサの
検出応答遅れを解析して真の空燃比を求め、それに基づ
いて1個のセンサの集合部出力から空燃比を制御する例
を示したが、それに限られるものではなく、排気系に気
筒の個数だけセンサを配置し、各センサ出力からその検
出応答遅れを勘案して各センサの真の出力を求め、それ
に基づいて各気筒の空燃比を目標値に制御しても良い。
In the above embodiment, an example in which the detection response delay of the air-fuel ratio sensor is analyzed to find the true air-fuel ratio, and the air-fuel ratio is controlled from the output of one sensor based on the true air-fuel ratio has been shown. However, the number of cylinders in the exhaust system is not limited to that, and the true output of each sensor is calculated from the output of each sensor in consideration of the detection response delay, and the air-fuel ratio of each cylinder is calculated based on that. You may control to a target value.

【0059】[0059]

【発明の効果】請求項1項記載の内燃機関の空燃比検出
方法にあっては、内燃機関の排気系に配置される空燃比
センサの検出値から各気筒空燃比を検出するものであ
て、前記センサの応答遅れを1次遅れ系で擬似的にモ
デル化し、その挙動を記述する状態方程式を設定し、前
記状態方程式を周期ΔTで離散化して伝達関数(l−α
ハット)/(Z−αハット)を求め、前記伝達関数の
数αハットを所定機関回転数ごとに相違させつつ前記伝
達関数の逆伝達関数(Z−αハット)/(l−αハッ
ト)を求め、求めた逆伝達関数を前記空燃比センサ
力に乗じて内燃機関に入力される混合気の空燃比の推定
値を求め、求めた空燃比の推定値を用いて前記内燃機関
の排気系の挙動を記述するモデルを設定すると共に、そ
の内部状態を観測するオブザーバを設定し、および前記
オブザーバの出力から各気筒の空燃比を検出する如く構
成したので、燃焼した混合気の真の空燃比を気筒ごとに
リアルタイムかつ正確に求めることができ、よって気筒
ごとのフィードバック制御を行うときなども制御の応答
性が向上すると共に、精度と収束性に優れた空燃比制御
を実現することが可能となる。また、気筒ごとに空燃比
センサを配置する必要がなく、単一の空燃比センサの出
力に基づいて各気筒の空燃比を検出あるいは推定するこ
とができる。
[Effect of the Invention] In the air-fuel ratio detection method for an internal combustion engine according to claim 1 wherein, der detects the air-fuel ratio in each cylinder from the detected value of the air-fuel ratio sensor disposed in an exhaust system of an internal combustion engine
And Tsu, artificially modeled response delay of the sensor in first-order system, set the state equation describing the behavior, the transfer function by discretizing the state equation in the period [Delta] T (l-alpha
Hat) / seeking (Z-alpha hat), engagement of the transfer function
The number α hat is changed for each predetermined engine speed and
Inverse transfer function of the reaching function (Z-α hat) / (l-α hat
Seeking g) to obtain the estimated value of the air-fuel ratio of the mixture supplied to the internal combustion engine by multiplying the inverse transfer function obtained in output <br/> force of the air-fuel ratio sensor, the estimated value of the air-fuel ratio determined Using the internal combustion engine
A model describing the behavior of the exhaust system of
Set an observer to observe the internal state of
Having rather configured to detect the air-fuel ratio of each cylinder from the output of the observer, the true air-fuel ratio of the mixture burned in each cylinder
Real-time and accurate determination is possible, so cylinder
Control response even when performing feedback control for each
As a result, it is possible to realize air-fuel ratio control with improved accuracy and excellent accuracy and convergence . Also, the air-fuel ratio for each cylinder
There is no need to place a sensor, and a single air-fuel ratio sensor
The air-fuel ratio of each cylinder can be detected or estimated based on the force.
Can be.

【0060】請求項2項記載の内燃機関の空燃比検出方
法にあっては、前記周期ΔTを所定機関回転数ごとに相
違させる如く構成したので、フィードバック制御を行う
ときなどにその制御周期を機関回転数に同期させる場合
にも精度良く空燃比を目標値に収束させることができ
る。
In the method for detecting the air-fuel ratio of the internal combustion engine according to the second aspect of the present invention, since the cycle ΔT is made different for each predetermined engine speed, the control cycle is set to the engine when feedback control is performed. The air-fuel ratio can be accurately converged to the target value even when synchronized with the rotation speed.

【0061】[0061]

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

【図1】この発明に係る方法を実現する内燃機関の空燃
比検出・制御装置を全体的に示す概略図である。
FIG. 1 is a schematic diagram generally showing an air-fuel ratio detection / control device for an internal combustion engine that implements a method according to the present invention.

【図2】図1中の制御ユニットの構成を示すブロック図
である。
FIG. 2 is a block diagram showing a configuration of a control unit in FIG.

【図3】単気筒の内燃機関において吸気量を一定として
供給燃料量をステップ状に変化させた場合の空燃比セン
サの応答遅れを示すシミュレーション結果、同じ場合の
LAFセンサ出力の実測値を表すデータである。
FIG. 3 is a simulation result showing a response delay of an air-fuel ratio sensor when a supply amount of fuel is changed stepwise with a constant intake air amount in a single-cylinder internal combustion engine, and data representing an actual measurement value of an LAF sensor output in the same case; It is.

【図4】空燃比センサの検出動作をモデル化した例を示
すブロック線図である。
FIG. 4 is a block diagram showing an example in which a detection operation of an air-fuel ratio sensor is modeled.

【図5】図4に示すモデルを周期ΔTで離散化したモデ
ルである。
FIG. 5 is a model obtained by discretizing the model shown in FIG. 4 with a period ΔT.

【図6】この発明に係る空燃比センサの検出挙動をモデ
ル化した真の空燃比推定器を示すブロック線図である。
FIG. 6 is a block diagram showing a true air-fuel ratio estimator that models the detection behavior of the air-fuel ratio sensor according to the present invention.

【図7】LAFセンサ出力の実測値と、それに逆伝達関
数を乗じて求めた真の空燃比の推定値を表すグラフであ
る。
FIG. 7 is a graph showing the measured value of the LAF sensor output and the estimated value of the true air-fuel ratio obtained by multiplying the measured value by the inverse transfer function.

【図8】この発明で用いる内燃機関の排気系の挙動を示
すモデルを表すブロック線図である。
FIG. 8 is a block diagram showing a model showing a behavior of an exhaust system of an internal combustion engine used in the present invention.

【図9】図8に示すモデルを用いて4気筒内燃機関につ
いて3気筒の空燃比を14.7に、1気筒の空燃比を1
2.0にして燃料を供給する場合を示すデータである。
FIG. 9 shows an example in which the air-fuel ratio of three cylinders is set to 14.7 and the air-fuel ratio of one cylinder is set to 1 using a model shown in FIG.
It is data showing the case where the fuel is supplied at 2.0.

【図10】図9に示す入力を与えたときの図8モデルの
集合部の空燃比を表すデータである。
FIG. 10 shows data representing the air-fuel ratio of the aggregate of the model of FIG. 8 when the input shown in FIG. 9 is given.

【図11】図9に示す入力を与えたときの図8モデルの
集合部の空燃比をLAFセンサの応答遅れを補正せずに
表したデータと、同じときのLAFセンサ出力の実測値
を比較するグラフである。
FIG. 11 is a comparison of the data showing the air-fuel ratio of the collecting portion of the model of FIG. 8 when the input shown in FIG. 9 is given without correcting the response delay of the LAF sensor and the actual measurement value of the LAF sensor output at the same time. It is a graph to do.

【図12】一般的なオブザーバの構成を示すブロック線
図である。
FIG. 12 is a block diagram showing a configuration of a general observer.

【図13】この発明で用いるオブザーバの構成を示すブ
ロック線図である。
FIG. 13 is a block diagram showing a configuration of an observer used in the present invention.

【図14】図8に示すモデルと図13に示すオブザーバ
を組み合わせたシミュレーションモデルを表すブロック
線図である。
FIG. 14 is a block diagram showing a simulation model in which the model shown in FIG. 8 and the observer shown in FIG. 13 are combined.

【図15】図14のモデルを用いて4気筒内燃機関につ
いて3気筒の空燃比を14.7に、1気筒の空燃比を1
2.0にして燃料を供給した場合のシミュレーション結
果を示すデータである。
15 shows a four-cylinder internal combustion engine with an air-fuel ratio of 14.7 and an air-fuel ratio of one cylinder of 1 using the model of FIG.
It is data which shows the simulation result at the time of supplying fuel to 2.0.

【図16】図15の例に仮想ノイズを含めて燃料を供給
した場合のシミュレーション結果を示すデータである。
FIG. 16 is data showing a simulation result when fuel is supplied including virtual noise in the example of FIG. 15;

【図17】図16に示す場合について時系列に推定偏差
を示す説明図である。
FIG. 17 is an explanatory diagram showing estimated deviations in a time series in the case shown in FIG. 16;

【図18】空燃比センサ出力の実測値を図6に示す空燃
比推定器に入力して得られた実集合部空燃比を入力した
ときのオブザーバ出力を示す説明図である。
FIG. 18 is an explanatory view showing an observer output when an actual assembly air-fuel ratio obtained by inputting an actual measurement value of an air-fuel ratio sensor output to the air-fuel ratio estimator shown in FIG. 6 is input;

【図19】図13に示すオブザーバ出力から得られる値
に基づいて空燃比を目標値に気筒ごとにPID制御する
例を示すブロック線図である。
FIG. 19 is a block diagram showing an example of performing PID control for each cylinder to set an air-fuel ratio to a target value based on a value obtained from the observer output shown in FIG.

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

10 内燃機関 18 インテークマニホルド 20 インジェクタ 22 エキゾーストパイプ 40 空燃比センサ 42 制御ユニット Reference Signs List 10 internal combustion engine 18 intake manifold 20 injector 22 exhaust pipe 40 air-fuel ratio sensor 42 control unit

───────────────────────────────────────────────────── フロントページの続き (72)発明者 小森谷 勲 埼玉県和光市中央1丁目4番1号 株式 会社本田技術研究所内 (56)参考文献 特開 平2−238150(JP,A) 特開 平1−211648(JP,A) 特開 平2−301805(JP,A) 特開 平1−280647(JP,A) 特開 昭59−120726(JP,A) 特開 昭59−101562(JP,A) 特開 平2−173334(JP,A) 実開 平3−54260(JP,U) ─────────────────────────────────────────────────── ─── Continuation of the front page (72) Inventor Isao Komoritani 4-1-1 Chuo, Wako-shi, Saitama Inside Honda R & D Co., Ltd. (56) Reference JP-A-2-238150 (JP, A) JP-A 1-211648 (JP, A) JP-A 2-301805 (JP, A) JP-A 1-280647 (JP, A) JP-A 59-120726 (JP, A) JP-A 59-101562 (JP, A) A) Japanese Unexamined Patent Publication No. 2-173334 (JP, A) Actual Development No. 3-54260 (JP, U)

Claims (2)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 内燃機関の排気系に配置される空燃比セ
ンサの検出値から各気筒空燃比を検出するものであっ
て、 a.前記センサの応答遅れを1次遅れ系で擬似的にモデ
ル化し、その挙動を記述する状態方程式を設定し、 b.前記状態方程式を周期ΔTで離散化して伝達関数
(l−αハット)/(Z− αハット)を求め、 c.前記伝達関数の係数αハットを所定機関回転数ごと
に相違させつつ前記伝達 関数の逆伝達関数(Z−αハッ
ト)/(l−αハット)を求め、求めた逆伝達関数を前記空燃比センサ出力に乗じ
て内燃機関に入力される混合気の空燃比の推定値を求
求めた空燃比の推定値を用いて前記内燃機関の排気
系の挙動を記述するモデ ルを設定すると共に、その内部
状態を観測するオブザーバを設定し、 および 前記オブザーバの出力から各気筒の空燃比を検出す
る、ことからなることを特徴とする内燃機関の空燃比検
出方法。
1. A meet those for detecting the air-fuel ratio of each cylinder from the detected value of the air-fuel ratio sensor disposed in an exhaust system of an internal combustion engine Te <br/>, a. A response delay of the sensor is simulated by a first-order delay system and a state equation describing the behavior thereof is set ; b. The above equation of state is discretized with a period ΔT to transfer function
Calculate (1−α hat) / (Z− α hat) , c. The coefficient α hat of the transfer function is set for each predetermined engine speed.
Inverse transfer function of the transfer function while differences (Z-alpha Hatton
G)) / (l-α hat) , and d . The inverse transfer function obtained by multiplying the output of the air-fuel ratio sensor obtains an estimate of the air-fuel ratio of the mixture supplied to the internal combustion engine, e. Exhaust gas of the internal combustion engine using the obtained estimated value of the air-fuel ratio
It sets a model describing the behavior of the system, its interior
Set observers to observe states, and f . The air-fuel ratio of each cylinder is detected from the output of the observer.
That, the air-fuel ratio detection method for an internal combustion engine, characterized in that it consists.
【請求項2】 前記周期ΔTを所定機関回転数ごとに相
違させることを特徴とする請求項1項記載の内燃機関の
空燃比検出方法。
2. The method of detecting an air-fuel ratio of an internal combustion engine according to claim 1, wherein the cycle ΔT is made different for each predetermined engine speed.
JP3359339A 1991-12-27 1991-12-27 Air-fuel ratio detection method for internal combustion engine Expired - Fee Related JP2689362B2 (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP3359339A JP2689362B2 (en) 1991-12-27 1991-12-27 Air-fuel ratio detection method for internal combustion engine
EP92311841A EP0553570B1 (en) 1991-12-27 1992-12-29 Method for detecting and controlling air-fuel ratio in internal combustion engines
DE69225212T DE69225212T2 (en) 1991-12-27 1992-12-29 Method for determining and controlling the air / fuel ratio in an internal combustion engine
US08/282,104 US5524598A (en) 1991-12-27 1994-07-28 Method for detecting and controlling air-fuel ratio in internal combustion engine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP3359339A JP2689362B2 (en) 1991-12-27 1991-12-27 Air-fuel ratio detection method for internal combustion engine

Publications (2)

Publication Number Publication Date
JPH05180059A JPH05180059A (en) 1993-07-20
JP2689362B2 true JP2689362B2 (en) 1997-12-10

Family

ID=18464008

Family Applications (1)

Application Number Title Priority Date Filing Date
JP3359339A Expired - Fee Related JP2689362B2 (en) 1991-12-27 1991-12-27 Air-fuel ratio detection method for internal combustion engine

Country Status (1)

Country Link
JP (1) JP2689362B2 (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0802316B1 (en) * 1994-02-04 2000-04-12 Honda Giken Kogyo Kabushiki Kaisha Air/fuel ratio estimation system for internal combustion engine
DE69507060T2 (en) * 1994-02-04 1999-05-20 Honda Motor Co Ltd Air / fuel ratio estimation system for an internal combustion engine
US5600056A (en) * 1994-06-20 1997-02-04 Honda Giken Kogyo Kabushiki Kaisha Air/fuel ratio detection system for multicylinder internal combustion engine
US5806506A (en) * 1996-08-01 1998-09-15 Honda Giken Kogyo Kabushiki Kaisha Cylinder-by-cylinder air-fuel ratio-estimating system for internal combustion engines
JP3913986B2 (en) 2001-01-09 2007-05-09 三菱電機株式会社 Air-fuel ratio control device for internal combustion engine
JP4697201B2 (en) 2007-07-19 2011-06-08 トヨタ自動車株式会社 Abnormality detection device for internal combustion engine
JP4985350B2 (en) * 2007-11-27 2012-07-25 トヨタ自動車株式会社 Air-fuel ratio sensor abnormality diagnosis device

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01211648A (en) * 1988-02-17 1989-08-24 Nissan Motor Co Ltd Fuel injection controller of internal combustion engine
JP2754676B2 (en) * 1989-03-10 1998-05-20 トヨタ自動車株式会社 Fuel injection amount control device for internal combustion engine
JPH02301805A (en) * 1989-05-17 1990-12-13 Nippondenso Co Ltd Digital controller

Also Published As

Publication number Publication date
JPH05180059A (en) 1993-07-20

Similar Documents

Publication Publication Date Title
JP2717744B2 (en) Air-fuel ratio detection and control method for internal combustion engine
EP0553570B1 (en) Method for detecting and controlling air-fuel ratio in internal combustion engines
JP3162521B2 (en) Air-fuel ratio estimator for each cylinder of internal combustion engine
EP1705359B1 (en) Method of feedforward controlling a multi-cylinder internal combustion engine and relative feedforward fuel injection control system
EP0688945B1 (en) Air/fuel ratio detection system for multicylinder internal combustion engine
EP0582085A2 (en) Fuel metering control system and cylinder air flow estimation method in internalcombustion engine
US7287525B2 (en) Method of feedforward controlling a multi-cylinder internal combustion engine and associated feedforward fuel injection control system
JP2689364B2 (en) Fuel injection amount control device for internal combustion engine
JP2832944B2 (en) Measurement data delay compensation method
JP2684011B2 (en) Internal combustion engine abnormality determination device
EP0643211B1 (en) Air-fuel ratio estimator for internal combustion engine
JP2689362B2 (en) Air-fuel ratio detection method for internal combustion engine
JP2683985B2 (en) Fuel injection amount control device for internal combustion engine
JPH0674076A (en) Intake air amount calculating method for internal combustion engine
JP2683974B2 (en) Air-fuel ratio control method for internal combustion engine
EP0026642B1 (en) A method for controlling the supply of fuel to an internal combustion engine
JP3162567B2 (en) Air-fuel ratio estimator for each cylinder of internal combustion engine
JP2689368B2 (en) Air-fuel ratio estimating device for each cylinder of internal combustion engine
JP3361316B2 (en) Cylinder air-fuel ratio estimating device and control device for internal combustion engine
JP2683986B2 (en) Fuel injection amount control device for internal combustion engine
JP2576184B2 (en) Fuel injection amount control device for internal combustion engine
JPH109038A (en) Air-fuel ratio detecting device of engine
JP3683355B2 (en) Cylinder air-fuel ratio estimation device for internal combustion engine
JP2589078B2 (en) Air-fuel ratio controller for multi-cylinder engine
JP2745800B2 (en) Idling speed control system for internal combustion engine

Legal Events

Date Code Title Description
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 19970715

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080829

Year of fee payment: 11

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090829

Year of fee payment: 12

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100829

Year of fee payment: 13

LAPS Cancellation because of no payment of annual fees