JP2001073844A - Air-fuel ratio presuming device for each cylinder of internal-combustion engine and control device - Google Patents

Air-fuel ratio presuming device for each cylinder of internal-combustion engine and control device

Info

Publication number
JP2001073844A
JP2001073844A JP2000254398A JP2000254398A JP2001073844A JP 2001073844 A JP2001073844 A JP 2001073844A JP 2000254398 A JP2000254398 A JP 2000254398A JP 2000254398 A JP2000254398 A JP 2000254398A JP 2001073844 A JP2001073844 A JP 2001073844A
Authority
JP
Japan
Prior art keywords
air
cylinder
fuel ratio
combustion engine
histories
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2000254398A
Other languages
Japanese (ja)
Other versions
JP3361316B2 (en
Inventor
Yusuke Hasegawa
祐介 長谷川
Isao Komoriya
勲 小森谷
Eisuke Kimura
英輔 木村
Naosuke Akasaki
修介 赤崎
Yoichi Nishimura
要一 西村
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 JP2000254398A priority Critical patent/JP3361316B2/en
Publication of JP2001073844A publication Critical patent/JP2001073844A/en
Application granted granted Critical
Publication of JP3361316B2 publication Critical patent/JP3361316B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • Y02T10/18

Landscapes

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

Abstract

PROBLEM TO BE SOLVED: To enhance the presumptive accuracy of the air-fuel ratio of each cylinder by an observer and implement the feedback control of the air-fuel ratio using the obtained estimate. SOLUTION: In case the number of combustion history to be taken into consideration with an internal-combustion engine is three or more, the weighting factor of the weighted mean for the prior combustion history excluding the combustion history more than two immediately before and less than the mentioned number of combustion history to be taken into consideration is made zero. Estimate is used for feedback control of the air-fuel ratio.

Description

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

【0001】[0001]

【発明の属する技術分野】この発明は内燃機関の気筒別
空燃比推定装置および制御装置に関し、より具体的には
多気筒内燃機関の排気系集合部に空燃比センサを設け、
排気系の挙動を記述するモデルを設定してセンサ出力を
入力すると共に、その内部状態を観測するオブザーバを
設け、その出力から各気筒の空燃比を推定する内燃機関
の気筒別空燃比推定装置において、オブザーバによる各
気筒の空燃比の推定精度を向上させると共に、前記推定
値に基づいて空燃比をフィードバック制御するようにし
たものに関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a cylinder-by-cylinder air-fuel ratio estimating apparatus and control apparatus for an internal combustion engine, and more specifically, an air-fuel ratio sensor is provided in an exhaust system assembly of a multi-cylinder internal combustion engine.
A model describing the behavior of the exhaust system is set, a sensor output is input, and an observer for observing the internal state is provided, and a cylinder-by-cylinder air-fuel ratio estimating apparatus for estimating the air-fuel ratio of each cylinder from the output is provided. In addition, the present invention relates to a method of improving the accuracy of estimating the air-fuel ratio of each cylinder by an observer and performing feedback control of the air-fuel ratio based on the estimated value.

【0002】[0002]

【従来の技術】内燃機関の排気系に空燃比センサを設け
て空燃比を検出することは良く行われており、その一例
として特開昭59−101562号公報記載の技術を挙
げることができる。また、本出願人も先に特願平3−3
59339号(特開平5−180059号)において、
排気系の挙動を記述するモデルを設定して排気系集合部
に設けた単一の空燃比センサの出力を入力し、オブザー
バを介して各気筒の空燃比を推定する技術を提案してい
る。尚、そこにおいて、空燃比センサは広域空燃比セン
サ、即ち、理論空燃比で出力が反転するO2 センサでは
なく、理論空燃比の前後を通じて排気ガス中の酸素濃度
に比例した出力特性を有するものを使用している。
2. Description of the Related Art It is common practice to provide an air-fuel ratio sensor in an exhaust system of an internal combustion engine to detect an air-fuel ratio, and an example thereof is disclosed in Japanese Patent Application Laid-Open No. Sho 59-101562. In addition, the present applicant has also previously filed Japanese Patent Application No. Hei.
No. 59339 (JP-A-5-180059),
A technique has been proposed in which a model describing the behavior of the exhaust system is set, the output of a single air-fuel ratio sensor provided in the exhaust system assembly is input, and the air-fuel ratio of each cylinder is estimated via an observer. Here, the air-fuel ratio sensor is not a wide-range air-fuel ratio sensor, that is, an O 2 sensor whose output is inverted at the stoichiometric air-fuel ratio, but has an output characteristic proportional to the oxygen concentration in the exhaust gas before and after the stoichiometric air-fuel ratio. You are using

【0003】[0003]

【発明が解決しようとする課題】上記した構成によって
各気筒の空燃比を精度良く推定することができたが、運
転状態は頻繁に変化し、それに応じて演算パラメータを
変更する必要があるが、実際の制御において演算時間は
必ずしも十分に得ることができない。
Although the air-fuel ratio of each cylinder can be accurately estimated by the above-described configuration, the operating condition frequently changes, and it is necessary to change the calculation parameters accordingly. In actual control, a sufficient calculation time cannot always be obtained.

【0004】従って、この発明の目的は先に提案した技
術の改良にあり、演算を簡略化してオブザーバによる各
気筒の空燃比の推定精度を一層向上させた内燃機関の気
筒別空燃比推定装置を提供することを目的とする。
SUMMARY OF THE INVENTION Accordingly, an object of the present invention is to improve the technique proposed above, and to provide a cylinder-by-cylinder air-fuel ratio estimating apparatus for an internal combustion engine, in which the calculation is simplified and the observer is further improved in estimating the air-fuel ratio of each cylinder. The purpose is to provide.

【0005】更には、前記の如くして推定された値に基
づいて空燃比フィードバック制御を行うようにした内燃
機関の気筒別空燃比フィードバック制御装置を提供する
ことを目的とする。
It is a further object of the present invention to provide a cylinder-by-cylinder air-fuel ratio feedback control device for performing air-fuel ratio feedback control based on the value estimated as described above.

【0006】[0006]

【課題を解決するための手段】上記の目的を解決するた
めに請求項1項にあっては、多気筒内燃機関の排気系集
合部に空燃比センサを配置してその出力から各気筒の入
力混合気の空燃比を推定する装置であって、前記センサ
の出力値を各気筒の燃焼履歴に重み係数Cnを乗じた加
重平均値からなるものとみなして排気系の挙動を記述す
るモデルを設定して各気筒の空燃比を状態変数とする状
態方程式を設定し、その内部状態を観測するオブザーバ
を設定してその出力を求める第1の手段、および前記オ
ブザーバ出力から各気筒の空燃比を推定する第2の手段
を有する内燃機関の気筒別空燃比推定装置において、前
記内燃機関で考慮すべき燃焼履歴の数を3以上とすると
き、直前の2つの燃焼履歴以上で、かつ前記考慮すべき
燃焼履歴の数より少ない燃焼履歴を除く、それ以前の燃
焼履歴についての加重平均の重み係数Cnを零にする如
く構成した。
According to a first aspect of the present invention, an air-fuel ratio sensor is disposed in an exhaust system collecting portion of a multi-cylinder internal combustion engine, and an input of each cylinder is obtained from an output of the sensor. A device for estimating an air-fuel ratio of an air-fuel mixture, wherein a model describing an exhaust system behavior is set by regarding the output value of the sensor as a weighted average value obtained by multiplying a combustion history of each cylinder by a weight coefficient Cn. First means for setting the state equation using the air-fuel ratio of each cylinder as a state variable, setting an observer for observing the internal state thereof, and obtaining the output thereof, and estimating the air-fuel ratio of each cylinder from the output of the observer When the number of combustion histories to be considered in the internal combustion engine is three or more, in the cylinder-by-cylinder air-fuel ratio estimation device having the second means, From the number of combustion history Except no combustion history, and as configured to zero the weighting factors Cn weighted average of the previous combustion history.

【0007】請求項2項に係る内燃機関の気筒別空燃比
フィードバック制御装置にあっては、多気筒内燃機関の
排気系集合部に配置された空燃比センサの出力値を各気
筒の燃焼履歴に重み係数Cnを乗じた加重平均値からな
るものとみなすと共に、前記内燃機関で考慮すべき燃焼
履歴の数を3以上とするとき、直前の2つの燃焼履歴以
上で、かつ前記考慮すべき燃焼履歴の数より少ない燃焼
履歴を除く、それ以前の燃焼履歴についての加重平均の
重み係数Cnを零にして排気系の挙動を記述するモデル
を設定して各気筒の空燃比を状態変数とする状態方程式
を設定し、その内部状態を観測するオブザーバを設定し
てその出力を求める第1の手段、前記オブザーバ出力か
ら各気筒の空燃比#nA/Fを推定する第2の手段、前
記推定された各気筒の空燃比#nA/Fを目標値にフィ
ードバック制御する気筒毎フィードバック補正項#nK
LAFを演算する第3の手段、および前記気筒毎フィー
ドバック補正項#nKLAFに基づいて前記内燃機関の
空燃比をフィードバック制御するフィードバック制御手
段を備える如く構成した。
According to a second aspect of the present invention, there is provided an air-fuel ratio feedback control device for each cylinder of an internal combustion engine, wherein an output value of an air-fuel ratio sensor arranged in an exhaust system collecting section of a multi-cylinder internal combustion engine is stored in a combustion history of each cylinder. When the number of combustion histories to be considered in the internal combustion engine is 3 or more and the number of combustion histories to be considered in the internal combustion engine is 3 or more, A state equation using the air-fuel ratio of each cylinder as a state variable by setting a model that describes the behavior of the exhaust system by setting the weighted average weight coefficient Cn of the previous combustion history to zero, excluding the combustion history less than A first means for setting an observer for observing the internal state thereof and obtaining its output; a second means for estimating an air-fuel ratio # nA / F of each cylinder from the observer output; Mind Air-fuel ratio # nA / F section each cylinder feedback correction performs feedback control to the target value #nK of
A third means for calculating the LAF and a feedback control means for feedback-controlling the air-fuel ratio of the internal combustion engine based on the cylinder-by-cylinder feedback correction term #nKLAF are provided.

【0008】請求項3項にあっては、さらに、前記空燃
比センサの出力から前記排気系集合部の空燃比を検出
し、検出した排気系集合部の空燃比を目標値に制御する
ための集合部フィードバック補正項KLAFを演算する
第4の手段を備え、前記フィードバック制御手段は、前
記集合部フィードバック補正項KLAFと気筒毎フィー
ドバック補正項#nKLAFに基づいて前記内燃機関の
空燃比をフィードバック制御する如く構成した。
According to a third aspect of the present invention, the air-fuel ratio of the exhaust system collecting section is detected from the output of the air-fuel ratio sensor, and the detected air-fuel ratio of the exhaust system collecting section is controlled to a target value. A fourth means for calculating a collective feedback correction term KLAF; wherein the feedback control means performs feedback control of the air-fuel ratio of the internal combustion engine based on the collective feedback correction term KLAF and a cylinder-by-cylinder feedback correction term #nKLAF; It was configured as follows.

【0009】[0009]

【作用】請求項1項にあっては、加重平均の重み係数C
nは運転状態に応じて変化するが、例えば多気筒内燃機
関において直前の2つの燃焼履歴以前の燃焼履歴のそれ
を零にすれば、重み係数の合計値は1であるので、その
2つの燃焼履歴のいずれかについて同定し、他方は1よ
り減算することで求められる。その結果、4気筒である
とき、1気筒分の係数を同定すれば足り、3気筒分のそ
れを同定するのに比較すればオブザーバ行列の推定精度
が向上し、オブザーバの係数同定精度が向上するので、
各気筒空燃比の推定精度が向上する。
According to the first aspect, the weighting coefficient C of the weighted average
Although n changes according to the operating state, for example, in a multi-cylinder internal combustion engine, if the value of the combustion history before the two immediately preceding combustion histories is set to zero, the total value of the weighting coefficients is 1, so that the two combustion histories are 1. One of the histories is identified, and the other is obtained by subtracting 1 from the history. As a result, when the number of cylinders is four, it is sufficient to identify the coefficient for one cylinder, and when compared to identifying the coefficient for three cylinders, the estimation accuracy of the observer matrix is improved, and the coefficient identification accuracy of the observer is improved. So
The accuracy of estimating the air-fuel ratio of each cylinder is improved.

【0010】請求項1項において、「直前の2つの燃焼
履歴以上で、かつ前記考慮すべき燃焼履歴の数より少な
い燃焼履歴を除く、それ以前の燃焼履歴についての加重
平均の重み係数Cnを零にする」としたのは、推定すべ
き気筒数が多い内燃機関については、後で段落0058
で述べるように、重み係数を求めるべき気筒数を増加す
ることになる可能性があるからである。尚、オブザーバ
の次数自体は、推定気筒数の整数倍とする。
The weighting coefficient Cn of the weighted average of previous combustion histories, excluding combustion histories that are equal to or greater than the two immediately preceding combustion histories and less than the number of combustion histories to be considered, is set to zero. The reason is that, for an internal combustion engine having a large number of cylinders to be estimated,
This is because there is a possibility that the number of cylinders for which the weight coefficient is to be obtained will be increased, as will be described later. Note that the order of the observer itself is an integral multiple of the estimated number of cylinders.

【0011】請求項2項にあっては、請求項1項で述べ
たと同様の手法で推定された各気筒の空燃比#nA/F
を目標値にフィードバック制御する気筒毎フィードバッ
ク補正項#nKLAFを演算し、それに基づいて前記内
燃機関の空燃比をフィードバック制御する如く構成した
ので、請求項1項で述べた作用に加えて、単一のセンサ
を用いるのみで、各気筒の空燃比フィードバック制御を
精度良く行うことができる。
In the second aspect, the air-fuel ratio # nA / F of each cylinder estimated by the same method as described in the first aspect.
Is calculated based on the cylinder-by-cylinder feedback correction term #nKLAF, which feedback-controls the air-fuel ratio of the internal combustion engine to the target value. The air-fuel ratio feedback control of each cylinder can be performed with high accuracy only by using the above sensor.

【0012】尚、「前記内燃機関で考慮すべき燃焼履歴
の数を3以上とするとき、直前の2つの燃焼履歴以上
で、かつ前記考慮すべき燃焼履歴の数より少ない燃焼履
歴を除く、それ以前の燃焼履歴についての加重平均の重
み係数Cnを零にする」としたが、実施の形態において
は、前記内燃機関で考慮すべき燃焼履歴の数を3以上と
するとき、直前の2つの燃焼履歴を除く、それ以前の燃
焼履歴についての加重平均の重み係数Cnを零にするよ
うにした。
[0012] When the number of combustion histories to be considered in the internal combustion engine is three or more, the combustion histories excluding the combustion histories that are more than the two immediately preceding combustion histories and less than the number of combustion histories to be considered are excluded. The weighted average weight coefficient Cn of the previous combustion history is set to zero. "However, in the embodiment, when the number of combustion histories to be considered in the internal combustion engine is three or more, the two immediately preceding combustion The weight coefficient Cn of the weighted average for the previous combustion history excluding the history is set to zero.

【0013】請求項3項にあっては、さらに、前記空燃
比センサの出力から前記排気系集合部の空燃比を検出
し、検出した排気系集合部の空燃比を目標値に制御する
ための集合部フィードバック補正項KLAFを演算する
第4の手段を備え、前記フィードバック制御手段は、前
記集合部フィードバック補正項KLAFと気筒毎フィー
ドバック補正項#nKLAFに基づいて前記内燃機関の
空燃比をフィードバック制御する如く構成したので、請
求項2項で述べた作用に加えて、各気筒の空燃比は集合
部空燃比に収束し、集合部空燃比は目標空燃比に収束す
ることとなって、結果的に全ての気筒の空燃比が目標空
燃比に収束し、空燃比フィードバック制御の収束性を向
上させることができる。
According to a third aspect of the present invention, the air-fuel ratio of the exhaust system collecting section is detected from the output of the air-fuel ratio sensor, and the detected air-fuel ratio of the exhaust system collecting section is controlled to a target value. A fourth means for calculating a collective feedback correction term KLAF; wherein the feedback control means performs feedback control of the air-fuel ratio of the internal combustion engine based on the collective feedback correction term KLAF and a cylinder-by-cylinder feedback correction term #nKLAF; With such a configuration, in addition to the effect described in claim 2, the air-fuel ratio of each cylinder converges to the collecting air-fuel ratio, and the collecting air-fuel ratio converges to the target air-fuel ratio. The air-fuel ratios of all the cylinders converge to the target air-fuel ratio, and the convergence of the air-fuel ratio feedback control can be improved.

【0014】[0014]

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

【0015】図1はこの発明に係る内燃機関の気筒別空
燃比推定装置を実現するための空燃比フィードバック制
御装置を全体的に示す概略図である。図において符号1
0は4気筒の内燃機関を示しており、吸気路12の先端
に配置されたエアクリーナ14から導入された吸気は、
スロットル弁16でその流量を調節されつつインテーク
マニホルド18を経て第1ないし第4気筒に流入され
る。各気筒の吸気弁(図示せず)の付近にはインジェク
タ20が設けられて燃料を噴射する。噴射されて吸気と
一体となった混合気は、各気筒内で図示しない点火プラ
グで点火されて燃焼してピストン(図示せず)を駆動す
る。燃焼後の排気ガスは排気弁(図示せず)を介してエ
キゾーストマニホルド22に排出され、エキゾーストパ
イプ24を経て三元触媒コンバータ26で浄化されつつ
機関外に排出される。また、吸気路12には、スロット
ル弁配置位置付近に、それをバイパスするバイパス路2
8が設けられる。
FIG. 1 is a schematic diagram generally showing an air-fuel ratio feedback control device for realizing a cylinder-by-cylinder air-fuel ratio estimation device for an internal combustion engine according to the present invention. In the figure, reference numeral 1
0 indicates a four-cylinder internal combustion engine, and the intake air introduced from the air cleaner 14 disposed at the tip of the intake passage 12 is:
The gas flows into the first to fourth cylinders through the intake manifold 18 while the flow rate is adjusted by the throttle valve 16. An injector 20 is provided near an intake valve (not shown) of each cylinder to inject fuel. The air-fuel mixture injected and integrated with the intake air is ignited by an ignition plug (not shown) in each cylinder and burns to drive 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, the intake passage 12 has a bypass passage 2 that bypasses the vicinity of the throttle valve arrangement position.
8 are provided.

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

【0017】図2は制御ユニット42の詳細を示すブロ
ック図である。広域空燃比センサ40の出力は検出回路
46に入力され、そこで適当な線形化処理が行われ、理
論空燃比を中心としてリーンからリッチにわたる広い範
囲において排気ガス中の酸素濃度に比例したリニアな特
性からなる空燃比(A/F)が検出される。その詳細は
先に本出願人が提案した別の出願(特開平4−3694
71号)に述べられているので、これ以上の説明は省略
する。尚、以下の説明において、このセンサを「LAF
センサ」(リニア・エーバイエフ・センサ)と称する。
検出回路46の出力はA/D変換回路48を介してCP
U50,ROM52,RAM54などからなるマイクロ
コンピュータに取り込まれ、RAM54に格納される。
FIG. 2 is a block diagram showing details of the control unit 42. The output of the wide-range air-fuel ratio sensor 40 is input to a detection circuit 46, where appropriate linearization processing is performed, and a linear characteristic proportional to the oxygen concentration in the exhaust gas in a wide range from lean to rich around the stoichiometric air-fuel ratio. The air-fuel ratio (A / F) is detected. The details are described in another application previously proposed by the present applicant (Japanese Unexamined Patent Publication No.
No. 71), further description is omitted. In the following description, this sensor is referred to as “LAF
Sensor "(linear EV sensor).
The output of the detection circuit 46 is supplied to the CP via the A / D conversion circuit 48.
U50, ROM52, RAM54 and the like are taken into a microcomputer and stored in RAM54.

【0018】同様に、スロットル開度センサ36などの
アナログ出力はレベル変換回路56、マルチプレクサ5
8および第2のA/D変換回路60を介して、またクラ
ンク角センサ34の出力は波形整形回路62で波形整形
された後、カウンタ64で出力値がカウントされ、カウ
ント値はマイクロ・コンピュータ内に入力される。マイ
クロコンピュータにおいてCPU50は、ROM52に
格納された命令に従って検出値から制御値を演算し、駆
動回路66を介して各気筒のインジェクタ20を駆動す
ると共に、第2の駆動回路68を介して電磁弁70を駆
動し、図1に示したバイパス路28を通る2次空気量を
制御する。
Similarly, an analog output from the throttle opening sensor 36 and the like is supplied to a level conversion circuit 56 and a multiplexer 5.
After the output of the crank angle sensor 34 is shaped by the waveform shaping circuit 62, the output value is counted by the counter 64, and the count value is stored in the microcomputer. Is input to In the microcomputer, the CPU 50 calculates a control value from the detected value according to a command stored in the ROM 52, drives the injector 20 of each cylinder via the drive circuit 66, and controls the solenoid valve 70 via the second drive circuit 68. To control the amount of secondary air passing through the bypass 28 shown in FIG.

【0019】ここで、この発明の説明に入る前に、理解
の便宜上、先に提案した排気系の挙動を記述するモデル
について簡単に説明する。
Before describing the present invention, a model for describing the behavior of the exhaust system proposed earlier will be briefly described for the sake of convenience.

【0020】先ず、1個のLAFセンサの出力から各気
筒の空燃比を精度良く分離抽出するためには、LAFセ
ンサの検出応答遅れを正確に解明する必要がある。そこ
で、とりあえずこの遅れを1次遅れ系と擬似的にモデル
化し、図3に示す如きモデルを作成した。ここでLA
F:LAFセンサ出力、A/F:入力A/F、とする
と、その状態方程式は下記の数1で示すことができる。
First, in order to accurately separate and extract the air-fuel ratio of each cylinder from the output of one LAF sensor, it is necessary to accurately clarify the detection response delay of the LAF sensor. Therefore, this delay is tentatively modeled as a first-order delay system, and a model as shown in FIG. 3 is created. Here LA
When F: LAF sensor output and A / F: input A / F, the state equation can be expressed by the following equation (1).

【0021】[0021]

【数1】 (Equation 1)

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

【0023】[0023]

【数2】 (Equation 2)

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

【0025】[0025]

【数3】 (Equation 3)

【0026】[0026]

【数4】 (Equation 4)

【0027】具体的には数2をZ変換を用いて伝達関数
で示せば数5の如くになるので、その逆伝達関数を今回
のセンサ出力LAFに乗じることによって前回の入力空
燃比をリアルタイムに推定することができる。図5にそ
のリアルタイムのA/F推定器のブロック線図を示す。
More specifically, if Equation 2 is expressed as a transfer function using Z-transformation, Equation 5 is obtained, and the inverse input transfer function is multiplied by the current sensor output LAF to determine the previous input air-fuel ratio in real time. Can be estimated. FIG. 5 shows a block diagram of the real-time A / F estimator.

【0028】[0028]

【数5】 (Equation 5)

【0029】続いて、上記の如く求めた真の空燃比に基
づいて各気筒の空燃比を分離抽出する手法について説明
すると、先願でも述べたように、排気系の集合部の空燃
比を各気筒の空燃比の時間的な寄与度を考慮した加重平
均であると考え、時刻kのときの値を、数6のように表
した。尚、F(燃料量)を制御量としたため、ここでは
『燃空比F/A』を用いているが、後の説明においては
理解の便宜のため、支障ない限り「空燃比」を用いる。
尚、空燃比(ないしは燃空比)は、先に数5で求めた応
答遅れを補正した真の値を意味する。
Next, a method of separating and extracting the air-fuel ratio of each cylinder based on the true air-fuel ratio obtained as described above will be described. The value at the time k was considered as a weighted average in consideration of the temporal contribution of the air-fuel ratio of the cylinder, and was represented as in Expression 6. Note that "F / A" is used here because F (fuel amount) is a control amount, but "Air / fuel ratio" will be used in the following description for convenience of understanding unless there is a problem.
Note that the air-fuel ratio (or the fuel-air ratio) means a true value obtained by correcting the response delay previously obtained by Expression 5.

【0030】[0030]

【数6】 (Equation 6)

【0031】即ち、集合部の空燃比は、気筒ごとの過去
の燃焼履歴に重み係数Cn(例えば直近に燃焼した気筒
は40%、その前が30%...など)を乗じたものの
合算で表した。このモデルをブロック線図であらわす
と、図6のようになる。
That is, the air-fuel ratio of the collecting portion is the sum of the values obtained by multiplying the past combustion history of each cylinder by a weighting coefficient Cn (for example, 40% for the most recently burned cylinder, 30% before that, etc.). expressed. This model is represented by a block diagram as shown in FIG.

【0032】また、その状態方程式は数7のようにな
る。
The state equation is as shown in the following equation (7).

【0033】[0033]

【数7】 (Equation 7)

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

【0035】[0035]

【数8】 (Equation 8)

【0036】上記において、u(k)は観測不可能のた
め、この状態方程式からオブザーバを設計してもx
(k)は観測することができない。そこで4TDC前
(即ち、同一気筒)の空燃比は急激に変化しない定常運
転状態にあると仮定して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
(K) cannot be observed. Therefore, if x (k + 1) = x (k-3) assuming that the air-fuel ratio before 4TDC (that is, the same cylinder) is in a steady operation state in which the air-fuel ratio does not suddenly change, Equation 9 is obtained.

【0037】[0037]

【数9】 (Equation 9)

【0038】ここで、上記の如く求めたモデルについて
シミュレーション結果を示す。図7は4気筒内燃機関に
ついて3気筒の空燃比を14.7にし、1気筒だけ1
2.0にして燃料を供給した場合を示す。図8はそのと
きの集合部の空燃比を上記モデルで求めたものを示す。
同図においてはステップ状の出力が得られているが、こ
こで更にLAFセンサの応答遅れを考慮すると、センサ
出力は図9に「モデル出力値」と示すようになまされた
波形となる。図中「実測値」は同じ場合のLAFセンサ
出力の実測値であるが、これと比較し、上記モデルが多
気筒内燃機関の排気系を良くモデル化していることを検
証している。
Here, simulation results are shown for the model obtained as described above. FIG. 7 shows that the air-fuel ratio of three cylinders is set to 14.7 for a four-cylinder internal combustion engine,
The case where the fuel is supplied at 2.0 is shown. FIG. 8 shows the air-fuel ratio of the collecting portion at that time obtained by the above model.
Although a step-like output is obtained in the same figure, if the response delay of the LAF sensor is further taken into consideration, the sensor output has a waveform smoothed as shown in FIG. 9 as “model output value”. In the figure, “actual measurement value” is an actual measurement value of the LAF sensor output in the same case, and it is verified by comparison with this that the above model models the exhaust system of the multi-cylinder internal combustion engine well.

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

【0040】[0040]

【数10】 (Equation 10)

【0041】[0041]

【数11】 [Equation 11]

【0042】[0042]

【数12】 (Equation 12)

【0043】これよりA−KCを求めると、数13のよ
うになる。
When A-KC is obtained from this, the following equation (13) is obtained.

【0044】[0044]

【数13】 (Equation 13)

【0045】一般的なオブザーバの構成は図10に示さ
れるようになるが、今回のモデルでは入力u(k)がな
いので、図11に示すようにy(k)のみを入力とする
構成となり、これを数式で表すと数14のようになる。
Although the general observer configuration is as shown in FIG. 10, there is no input u (k) in the present model, so that only y (k) is input as shown in FIG. When this is expressed by a mathematical formula, it becomes as shown in Expression 14.

【0046】[0046]

【数14】 [Equation 14]

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

【0048】[0048]

【数15】 (Equation 15)

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

【0050】[0050]

【数16】 (Equation 16)

【0051】図12に上記したモデルとオブザーバを組
み合わせたものを示す。シミュレーション結果は先の出
願に示されているので省略するが、これにより集合部空
燃比より各気筒の空燃比を的確に抽出することができ
る。
FIG. 12 shows a combination of the above-described model and observer. The simulation result is omitted since it is shown in the earlier application, but the air-fuel ratio of each cylinder can be accurately extracted from the air-fuel ratio of the collecting portion.

【0052】オブザーバによって集合部空燃比より各気
筒空燃比を推定することができたことから、PIDなど
の制御則を用いて空燃比を気筒別に制御することが可能
となる。具体的には図13に示すように、センサ出力
(集合部A/F)と目標空燃比とからPID制御則を用
いて集合部フィードバック補正項KLAFを求めると共
に、オブザーバ推定値#nA/Fから気筒毎のフィード
バック補正項#nKLAF(n:気筒)を求める。気筒
毎のフィードバック補正項#nKLAFはより具体的に
は、集合部A/Fを気筒毎のフィードバック補正項#n
KLAFの平均値の前回演算値で除算して求めた目標値
と、オブザーバ推定値#nA/Fとの偏差を解消するよ
うにPID則を用いて求める。
Since the air-fuel ratio of each cylinder can be estimated from the air-fuel ratio of the collecting section by the observer, the air-fuel ratio can be controlled for each cylinder using a control law such as PID. More specifically, as shown in FIG. 13, an aggregate feedback correction term KLAF is obtained from the sensor output (aggregate A / F) and the target air-fuel ratio using a PID control law, and the observer estimated value # nA / F is used. A feedback correction term #nKLAF (n: cylinder) for each cylinder is obtained. More specifically, the feedback correction term #nKLAF for each cylinder is changed to the feedback correction term #n for each cylinder.
The target value obtained by dividing the average value of KLAF by the previous calculation value and the observer estimated value # nA / F are determined using the PID rule so as to eliminate the deviation.

【0053】これにより、各気筒の空燃比(A/F)は
集合部空燃比(A/F)に収束し、集合部空燃比(A/
F)は目標空燃比(A/F)に収束することとなって、
結果的に全ての気筒の空燃比(A/F)が目標空燃比
(A/F)に収束する。ここで、各気筒の燃料噴射量#
nTout (インジェクタ開弁時間で規定される)は、 #nTout =Tim×KCMD×KTOTAL×#nK
LAF×KLAF で求められる。上記で、Tim:基本値、KCMD:目
標空燃比、KTOTAL:その他の補正項、である。更
に、バッテリ補正などの加算項もあるが省略する。尚、
かかる制御の詳細は本出願人が先に提案した特願平5−
251138号に述べられているので、これ以上の説明
は省略する。
As a result, the air-fuel ratio (A / F) of each cylinder converges to the air-fuel ratio (A / F) at the collecting portion, and the air-fuel ratio (A / F) at the collecting portion.
F) converges to the target air-fuel ratio (A / F),
As a result, the air-fuel ratios (A / F) of all cylinders converge on the target air-fuel ratio (A / F). Here, the fuel injection amount of each cylinder #
nTout (defined by the injector valve open time) is: # nTout = Tim × KCMD × KTOTAL × # nK
It is determined by LAF × KLAF. In the above, Tim: basic value, KCMD: target air-fuel ratio, and KTOTAL: other correction terms. Further, there are addition terms such as battery correction, but they are omitted. still,
The details of such control are described in Japanese Patent Application No. Hei.
No. 251138, further description is omitted.

【0054】上記を前提としてこの発明を説明する。The present invention will be described on the premise of the above.

【0055】図14はこの発明の第1の実施の形態を示
す。前記した加重平均の重み係数C 1 ,C2 ,C3 ,C
4 について、今回排気(燃焼)気筒のそれをCn 、前
回、前々回および3TDC前の排気(燃焼)気筒のそれ
をそれぞれCn-1 ,Cn-2 ,C n-3 とするとき、前々回
および3TDC前の排気(燃焼)気筒についての重み係
数Cn-2 ,Cn-3 を零とし、図示の如く、直近(直前)
の2つの気筒について、即ち、直近(直前)の2つの燃
焼履歴についてのみ求める如くした。例えば、いま第1
気筒が排気(燃焼)気筒であるとすると、第3気筒と第
4気筒についての重み係数を零とするようにした。
FIG. 14 shows a first embodiment of the present invention.
You. Weighting coefficient C of the above-mentioned weighted average 1 , CTwo , CThree , C
Four About this time, it is C of the exhaust (combustion) cylinder.n,Previous
Times, 2 times before and 3TDC before exhaust (combustion) cylinder
To Cn-1 , Cn-2 , C n-3 And before and after
Weights for exhaust (combustion) cylinders before and 3TDC
Number Cn-2 , Cn-3 As zero, as shown in the figure, the most recent (immediately before)
Of the two cylinders, that is, the two most recent (immediately before)
Only the burning history was requested. For example, now the first
Assuming that the cylinder is an exhaust (combustion) cylinder, the third cylinder and the
The weight coefficient for the four cylinders was set to zero.

【0056】これは先の出願にも述べたが、集合部空燃
比の挙動を示す重み係数Cnは運転状態に応じて変化さ
せる必要があるが、いずれにしても重み係数Cnの値を
正確に知ることは困難である。そこで、4気筒内燃機関
において、排気(燃焼)順序において直近(直前)の2
つの燃焼履歴を除く、他の燃焼履歴の重み係数を零とし
た。尚、重み係数は全体として1であるから、直近(直
前)の燃焼履歴のいずれかについて求め、他方のそれは
求めた値を1から減算することで求められる(これは先
の出願において、重み係数を全ての気筒について求める
とした場合でも同様である)。即ち、図6の構成を図1
4のように修正した。
As described in the previous application, the weighting coefficient Cn indicating the behavior of the air-fuel ratio at the collecting section needs to be changed according to the operating state. In any case, the value of the weighting coefficient Cn is accurately determined. It is difficult to know. Therefore, in a four-cylinder internal combustion engine, the most recent (immediate) two in the exhaust (combustion) order
The weight coefficients of the other combustion histories except for one combustion history were set to zero. Since the weighting factor is 1 as a whole, it is obtained for one of the most recent (immediate) combustion histories, and the other is obtained by subtracting the obtained value from 1 (this is the weighting factor in the earlier application). Is the same for all cylinders). That is, the configuration of FIG.
Modified as in 4.

【0057】これについて、シミュレーションで検証し
たところ、良好な結果を得た。即ち、本来的に重み係数
Cnを3個、同定するのは現実問題として極めて困難で
あり、誤った同定値を用いてしまう恐れもある。上記の
如く構成することにより、同定すべき係数は1個で足り
ることから、演算間隔も短縮し、シミュレーションでの
検証において、かえって推定精度の向上を見た。尚、こ
の問題は気筒排気ポートから集合部空燃比センサ配置ま
での距離の遠近にも関係する。その距離が比較的短い場
合には、上記の如く構成することにより、推定精度を向
上させることができる。
When this was verified by simulation, good results were obtained. That is, it is extremely difficult as a practical matter to identify three weighting factors Cn from the beginning, and there is a possibility that an incorrect identification value may be used. With the above configuration, the number of coefficients to be identified is sufficient, so that the calculation interval is shortened, and the verification accuracy in simulation has been improved. This problem is also related to the distance from the cylinder exhaust port to the arrangement of the air-fuel ratio sensor in the collecting section. When the distance is relatively short, the configuration as described above can improve the estimation accuracy.

【0058】尚、実施の形態の場合は4気筒の内燃機関
を例にとったが、5気筒ないし6気筒までは重み係数は
同様に直近の2つの気筒について求めれば足ると予想さ
れる。それ以上の気筒数を持つ内燃機関については、重
み係数を求めるべき気筒数を増加することになる可能性
がある。
In the present embodiment, a four-cylinder internal combustion engine is taken as an example. However, it is expected that the weight coefficient for five to six cylinders will be sufficient if it is similarly obtained for the two nearest cylinders. For an internal combustion engine having a larger number of cylinders, there is a possibility that the number of cylinders for which the weight coefficient is to be obtained will be increased.

【0059】また、V型機関にあってはバンク毎に空燃
比センサが設けられることから、それぞれのバンクで排
気順序に従って重み係数を零とする気筒を選択する。例
えば、V型6気筒(右バンクの気筒を第1、第2、第
3、左バンクの気筒を第4、第5、第6とする)で、点
火順序を1(右)−4(左)−2(右)−5(左)−3
(右)−6(左)とすると、右バンクについて今第3気
筒が燃焼したとすると、第1気筒についての重み係数を
零とする。左バンクについて今第4気筒が燃焼したとす
れば、第5気筒についての重み係数を零とする。尚、こ
こで留意されるべきことは、オブザーバの次数はあくま
でも推定気筒数の整数倍、望ましくは1倍、この実施の
形態の場合は4であることである。
In the V-type engine, since an air-fuel ratio sensor is provided for each bank, a cylinder having a weight coefficient of zero in each bank is selected in accordance with the exhaust order. For example, in a V-type 6 cylinder (the right bank cylinders are first, second, and third cylinders, and the left bank cylinders are fourth, fifth, and sixth cylinders), and the ignition order is 1 (right) -4 (left ) -2 (right) -5 (left) -3
Assuming that (right) -6 (left), assuming that the third cylinder has just burned for the right bank, the weight coefficient for the first cylinder is set to zero. Assuming that the fourth cylinder has just burned for the left bank, the weight coefficient for the fifth cylinder is set to zero. It should be noted that the order of the observer is an integer multiple of the estimated number of cylinders, preferably one, and is four in this embodiment.

【0060】図15はこの発明の第2の実施の形態を示
す説明図である。
FIG. 15 is an explanatory diagram showing a second embodiment of the present invention.

【0061】オブザーバ行列は数14の式で示される。
そこで値A−KCは4×4の行列で示されるが、その中
で重み係数Cは上で述べたように機関運転状態に応じて
変化する。従って、この行列値について例えば機関回転
数と機関負荷とについて予想される変化範囲にわたって
予め解を求めておいてマップ化しておくことが考えられ
る。しかし、その場合でも、マップ値に該当する値がな
い際、補間演算で求めることは行列値として補償されて
いない故に、望ましくない。更に、数14の式は漸化式
であるので、演算の途中で係数が変化したときは演算を
やり直す必要がある。
The observer matrix is represented by the equation (14).
Therefore, the value A-KC is represented by a 4 × 4 matrix, in which the weight coefficient C changes according to the engine operating state as described above. Therefore, it is conceivable that a solution is obtained in advance over the expected change range of the engine speed and the engine load, and the matrix values are mapped. However, even in that case, when there is no value corresponding to the map value, it is not desirable to obtain the value by the interpolation operation because the matrix value is not compensated. Furthermore, since the equation of Equation 14 is a recurrence equation, it is necessary to repeat the operation when the coefficient changes during the operation.

【0062】そこで、図15に示すように、オブザーバ
行列を1からnまで複数個用意して平行的に演算させて
おき、機関運転状態に応じていずれかの行列値を選択す
るようにした。これによって演算時間を短縮することが
でき、オブザーバの各気筒空燃比の推定精度を向上させ
ることができる。尚、ここでオブザーバ行列の個数は2
個以上であればいくつでも良い。またオブザーバ行列の
演算も常に全てを同時に行う必要はなく、機関運転状態
の変化を予測して当該運転状態に対応するオブザーバ行
列についてのみ行わせても良い。
Therefore, as shown in FIG. 15, a plurality of observer matrices from 1 to n are prepared and operated in parallel, and one of the matrix values is selected according to the engine operating state. As a result, the calculation time can be shortened, and the accuracy of estimating the cylinder air-fuel ratio of the observer can be improved. Here, the number of observer matrices is 2
Any number is acceptable as long as it is equal to or more than the number. Also, it is not necessary to always perform all the operations of the observer matrix at the same time, and it is also possible to predict a change in the engine operating state and perform only the observer matrix corresponding to the operating state.

【0063】尚、先の出願ではオブザーバを設計する上
でリカッチ方程式の加重配分Q,Rを、Q/R=1/1
とした。この評価関数はオブザーバの性能上非常に重要
であるが、その後シミュレーションを通じて検証を重ね
た結果、Q/Rの値を大きくしていくと収束時間は短縮
していくが、あるレベルを超えると、ほとんど変わらな
いことが確認された。そこでその限界からQ/R≧10
0と設定したところ、良好な推定結果を得ることができ
た。即ち、評価関数は無次元数であるため、どのような
値も取ることができ、却って決定し難いものであるが、
検証を重ねた結果、Q/Rを100以上とすることで推
定精度の向上を見た。
In the prior application, when designing the observer, the weighted distributions Q and R of the Riccati equation are calculated as follows: Q / R = 1/1
And This evaluation function is very important for the performance of the observer. After verification through simulations, the convergence time decreases as the value of Q / R increases. It was confirmed that it hardly changed. Therefore, from the limit, Q / R ≧ 10
When set to 0, good estimation results could be obtained. That is, since the evaluation function is a dimensionless number, it can take any value and is rather difficult to determine.
As a result of repeated verification, it was found that the estimation accuracy was improved by setting the Q / R to 100 or more.

【0064】図16はこの発明の第3の実施の形態を示
すフロー・チャートである。第3の実施の形態は、いわ
ゆる可変バルブタイミング機構を備えた内燃機関につい
て空燃比を推定する場合に関する。
FIG. 16 is a flow chart showing a third embodiment of the present invention. The third embodiment relates to a case where the air-fuel ratio is estimated for an internal combustion engine having a so-called variable valve timing mechanism.

【0065】図17および図18を参照して可変バルブ
タイミング機構について簡単に説明する。
The variable valve timing mechanism will be briefly described with reference to FIGS.

【0066】可変バルブタイミング機構は図示例の場
合、吸気弁54と排気弁56の付近にそれぞれ設けられ
たロッカシャフト612に回転自在に配置される3個の
ロッカアーム614,616,618を備える(吸気側
と排気側とで構成が同一なため、図面で符号の後に添字
i(e)を付し、以下両者を共通に説明する)。該ロッ
カアームはそれぞれカムシャフト(図示せず)に取り付
けられた低速用の2個のカムと、低速用のカムに比して
径方向に突出した形状を備える1個の高速用のカム(共
に図示せず)に係合する。
In the case of the illustrated example, the variable valve timing mechanism includes three rocker arms 614, 616, 618 rotatably disposed on rocker shafts 612 provided near the intake valve 54 and the exhaust valve 56, respectively. Since the configuration is the same on the exhaust side and the exhaust side, a suffix i (e) is added after the reference numeral in the drawings, and both will be described in common below). Each of the rocker arms has two low-speed cams attached to a camshaft (not shown) and one high-speed cam having a shape protruding in the radial direction as compared with the low-speed cam (both figures are shown). (Not shown).

【0067】ロッカアーム614,616,618には
油路650、穴632,孔634,636およびピン6
40,642,644などからなる連結機構630を備
え、油路650と油圧源(図示せず)との間に介挿され
た油圧切り換え機構660は電磁弁680を介して油路
650に圧油を供給/停止し、それによって前記ピンが
前進/後退してロッカアームが連結/解放される。ロッ
カアームが連結されるときは高速カムの動きに基づい
て、解放されたときは低速カムの動きに基づいてバルブ
タイミングとリフト量とが決定される。尚、可変バルブ
タイミング機構の詳細は前記した従来技術に述べられて
いるので、説明はこの程度に止める。
The rocker arms 614, 616, and 618 have oil passages 650, holes 632, holes 634, 636, and pins 6.
40, 642, 644 and the like, and a hydraulic pressure switching mechanism 660 inserted between the oil passage 650 and a hydraulic source (not shown) is connected to the oil passage 650 via a solenoid valve 680. Is supplied / stopped, whereby the pin is advanced / retracted and the rocker arm is connected / released. The valve timing and the lift amount are determined based on the movement of the high-speed cam when the rocker arm is connected, and based on the movement of the low-speed cam when released. Since the details of the variable valve timing mechanism are described in the above-mentioned prior art, the description will be limited to this extent.

【0068】この連結/解放は具体的には図18に示す
如く、機関回転数Neと吸気圧力Pbとから決定され
る。以下、高速カムによる場合を『HiV/T』、低速
カムによる場合を『LoV/T』と言うが、HiV/T
においては、LoV/Tに比して、バルブタイミングが
早まって結果的にそのオーバラップ量が増加すると共
に、リフト量も増加する。尚、油圧切り換え機構660
には油圧スイッチからなるV/Tセンサ600が設けら
れ、油路650の油圧を通じてロッカアームの連結/解
放、即ち、バルブタイミングがLoV/TとHiV/T
のいずれに制御されているか検出する。図2に想像線で
示す如く、V/Tセンサ600の出力は、制御ユニット
42に入力される。
The connection / disconnection is specifically determined from the engine speed Ne and the intake pressure Pb, as shown in FIG. Hereinafter, the case using the high speed cam is referred to as “HiV / T”, and the case using the low speed cam is referred to as “LoV / T”.
In, the valve timing is advanced and the overlap amount is increased and the lift amount is increased as compared with LoV / T. The hydraulic pressure switching mechanism 660
Is provided with a V / T sensor 600 composed of a hydraulic switch, and the connection / release of the rocker arm through the hydraulic pressure of the oil passage 650, that is, the valve timing is LoV / T and HiV / T.
To which of the two is controlled. As indicated by the imaginary line in FIG. 2, the output of the V / T sensor 600 is input to the control unit 42.

【0069】集合部空燃比(A/F)の時間的な挙動を
捉えるためには、各気筒空燃比(A/F)推定オブザー
バ演算において空燃比(A/F)のサンプルタイミング
の適正化を図る必要がある。即ち、LAFセンサ出力は
明らかに機関回転数に同期した特性を示すが、制御ユニ
ット側はデジタルサンプリングであるため、タイミング
によってその挙動を捉えることができない場合が生じ
る。適正なサンプルタイミングは機関回転数や機関負荷
によって変化するが、特にバルブタイミングによる影響
が大きいと思われる。そこで、上記のごとき可変バルブ
タイミング機構を備えた内燃機関にあっては、機関回転
数、機関負荷、およびバルブタイミングとからサンプル
タイミングを決定するようにした。
In order to capture the temporal behavior of the air-fuel ratio (A / F) of the collecting section, the sampling timing of the air-fuel ratio (A / F) must be optimized in the calculation of the air-fuel ratio (A / F) estimation observer for each cylinder. It is necessary to plan. That is, although the output of the LAF sensor clearly shows a characteristic synchronized with the engine speed, since the control unit performs digital sampling, the behavior may not be able to be detected depending on the timing. Although the appropriate sample timing changes depending on the engine speed and the engine load, it is considered that the influence of the valve timing is particularly large. Therefore, in the internal combustion engine having the variable valve timing mechanism as described above, the sample timing is determined from the engine speed, the engine load, and the valve timing.

【0070】同様に、バルブタイミングにより集合部空
燃比(A/F)の挙動も変化するため、バルブタイミン
グの切り換えに同期してオブザーバ行列を変更する必要
がある。しかし、各気筒空燃比(A/F)の推定は瞬時
に行えるものではなく、オブザーバ演算の収束に演算数
回を要するため、バルブタイミング変更前のオブザーバ
行列を用いた演算と変更後のオブザーバ行列を用いた演
算とをオーバーラップして同時に行い、オブザーバ演算
が収束したと判断された後、即ち、バルブタイミング変
更動作が終了してから、各気筒空燃比(A/F)推定結
果を持ち替える。
Similarly, since the behavior of the air-fuel ratio (A / F) changes depending on the valve timing, it is necessary to change the observer matrix in synchronization with the switching of the valve timing. However, the estimation of each cylinder air-fuel ratio (A / F) cannot be performed instantaneously, and the convergence of the observer operation requires several operations. Therefore, the operation using the observer matrix before the valve timing change and the observer matrix after the change are performed. Are performed at the same time in an overlapped manner, and after it is determined that the observer calculation has converged, that is, after the valve timing changing operation is completed, the estimation result of each cylinder air-fuel ratio (A / F) is changed.

【0071】図16フロー・チャートを参照して説明す
ると、先ずS100において機関回転数Ne、吸気圧力
Pb、バルブタイミングV/Tを読み出し、S102に
進んでバルブタイミングがHiV/TかLoV/Tか判
断し、判断結果に従ってS104,106に進んでHi
ないしLoバルブタイミング用のタイミングマップを検
索する。図19にその特性を示す。
Referring to the flow chart of FIG. 16, first, at S100, the engine speed Ne, the intake pressure Pb, and the valve timing V / T are read, and the program proceeds to S102, at which the valve timing is HiV / T or LoV / T. Judge and proceed to S104, S106 according to the judgment result to Hi
Or search a timing map for Lo valve timing. FIG. 19 shows the characteristics.

【0072】LAFセンサ出力は実際の変局点に近い位
置でサンプリングするのが望ましいが、その変局点(ピ
ーク値)はセンサの反応時間を一定とすれば、機関回転
数が低くなるほど早いクランク角度で生じる。また機関
負荷が高いほど排気ガス圧力やボリュームが増加し、排
気ガスの流速が増してセンサへの到達時間が早まると予
想される。従って、マップの特性は、機関回転数Neが
低く、ないしは吸気圧力Pbが高いほど、早いクランク
角でサンプリングされた値を選択するように設定する。
ここで、「早い」とは、前のTDCにより近い位置でサ
ンプリングされた比較的古い値を意味する。
It is desirable that the LAF sensor output be sampled at a position close to the actual inflection point. However, if the inflection point (peak value) is constant for the sensor reaction time, the crank speed increases as the engine speed decreases. Occurs at an angle. Also, it is expected that the higher the engine load, the higher the exhaust gas pressure and volume, the higher the exhaust gas flow rate and the earlier the arrival time at the sensor. Therefore, the characteristics of the map are set such that the value sampled at a faster crank angle is selected as the engine speed Ne is lower or the intake pressure Pb is higher.
Here, “early” means a relatively old value sampled at a position closer to the previous TDC.

【0073】更に、バルブタイミングとの関連について
触れると、機関回転数の任意の値Ne1をLo側につい
てNe1−Lo、Hi側についてNe1−Hiとし、吸
気圧力についてもその任意の値をLo側についてPb1
−Lo、Hi側についてPb1−Hiとすると、マップ
の特性は、Pb1−Lo>Pb1−Hi、Ne1−Lo
>Ne1−Hiとする。即ち、HiV/Tにあっては排
気弁の開き時点がLoV/Tのそれより早いため、機関
回転数ないし吸気圧力の値が同一であれば、早期のサン
プリング値を選択するように、設定する。
Further, regarding the relationship with the valve timing, an arbitrary value Ne1 of the engine speed is Ne1-Lo for the Lo side, Ne1-Hi for the Hi side, and the arbitrary value of the intake pressure is Lo for the Lo side. Pb1
Assuming that Pb1−Hi for −Lo and Hi sides, the characteristics of the map are Pb1−Lo> Pb1−Hi, Ne1−Lo.
> Ne1-Hi. That is, in the case of HiV / T, the opening point of the exhaust valve is earlier than that of LoV / T, so that if the values of the engine speed or the intake pressure are the same, an earlier sampling value is selected. .

【0074】続いてS108に進んで検索したマップ値
からセンサ出力のサンプリングを行い、S110に進ん
でサンプリング値から数14に示すオブザーバ行列の演
算をHiV/Tについて行い、続いてS112に進んで
同様の演算をLoV/Tについて行う。続いてS114
に進んで再びバルブタイミングを判断し、判断結果に応
じてS116,118に進んで演算結果を選択して終わ
る。この演算結果に基づいて図13で説明した如き気筒
毎の空燃比フィードバック制御が行われる。
Then, the process proceeds to S108, where the sensor output is sampled from the retrieved map value, and the process proceeds to S110, where the observer matrix shown in Equation 14 is calculated from the sampled value for HiV / T, and then the process proceeds to S112. Is performed for LoV / T. Then S114
Then, the valve timing is determined again, and according to the determination result, the process proceeds to S116, 118 to select the calculation result and finish. Based on the calculation result, the air-fuel ratio feedback control for each cylinder as described with reference to FIG. 13 is performed.

【0075】第3の実施の形態は上記の如く構成したの
で、機関回転数、吸気圧力およびバルブタイミングを考
慮して空燃比センサ出力のサンプリング値を選択するの
で、検出精度を向上すると共に、バルブタイミングに応
じてオブザーバ行列を持ち替えるようにしたので、各気
筒の空燃比の推定精度が向上する。
Since the third embodiment is configured as described above, the sampling value of the air-fuel ratio sensor output is selected in consideration of the engine speed, the intake pressure and the valve timing. Since the observer matrix is changed according to the timing, the estimation accuracy of the air-fuel ratio of each cylinder is improved.

【0076】図20はこの発明の第4の実施の形態を示
す、図16と同様のフロー・チャートである。
FIG. 20 is a flow chart similar to FIG. 16, but showing a fourth embodiment of the present invention.

【0077】第4の実施の形態は第3の実施の形態の変
形例であり、以下説明すると、S200において機関回
転数Neなどを読み出した後S202に進んでLoV/
T用のタイミングマップを検索し、S204に進んでL
oV/T用にセンサ出力のサンプリング値を選択し、S
206に進んでLoV/T用にオブザーバ行列演算を行
う。続いてS208から212に進んでHiV/T側に
ついて同様の処理を行う。
The fourth embodiment is a modified example of the third embodiment. To be described below, after reading the engine speed Ne and the like in S200, the process proceeds to S202 and proceeds to LoV /
A timing map for T is searched, and the process proceeds to S204, where L
Select the sampling value of the sensor output for oV / T,
Proceeding to 206, an observer matrix operation is performed for LoV / T. Subsequently, the process proceeds from S208 to 212 to perform the same processing on the HiV / T side.

【0078】続いてS214に進んでバルブタイミング
を判断し、判断結果に応じてS216,218のいずれ
かに進んで演算結果を選択する。即ち、オブザーバ行列
演算において第3の実施の形態では選択されたLoない
しHiV/Tいずれかのサンプリング値を共通に使用し
て演算を行っていたのに対し、第4の実施の形態ではL
oおよびHiV/Tで別々のサンプリング値を用いて演
算を行う点で相違する。残余の構成および効果は第3の
実施の形態と同様である。
Then, the program proceeds to S214, in which the valve timing is determined, and in accordance with the result of the determination, the program proceeds to one of S216, 218 to select a calculation result. That is, in the observer matrix operation, the operation is performed by using the selected sampling value of any of Lo or HiV / T in the third embodiment, whereas in the fourth embodiment, the operation is performed in the fourth embodiment.
The difference is that calculations are performed using different sampling values for o and HiV / T. The remaining configuration and effects are the same as in the third embodiment.

【0079】尚、第3および第4の実施の形態におい
て、排気弁と吸気弁とが共に変更される可変バルブタイ
ミング機構を例にとったが、それに限られるものではな
く、既述の如く一方についてのみ変更される機構にも妥
当する。また吸気弁の休止が可能な機構にも妥当する。
更に、バルブタイミングを検出する手段として油圧スイ
ッチを例に上げたが、それに止まるものではなく、前記
したピン640,642,644の位置を検出しても良
く、あるいは可変バルブタイミング機構の制御ユニット
のバルブタイミング切り換え指令信号を参照して検出し
ても良い。
In the third and fourth embodiments, a variable valve timing mechanism in which both the exhaust valve and the intake valve are changed has been described as an example. However, the present invention is not limited to this. This also applies to the mechanism that is changed only for. Further, the present invention is also applicable to a mechanism capable of stopping the intake valve.
Further, although a hydraulic switch has been described as an example of a means for detecting the valve timing, the present invention is not limited to this, and the position of the pins 640, 642, 644 may be detected, or the control unit of the variable valve timing mechanism may be detected. The detection may be performed with reference to a valve timing switching command signal.

【0080】尚、上記した第1ないし第4の実施の形態
は、排気系の挙動を記述するモデルを設定し、その内部
状態を観測するオブザーバを使用して空燃比制御を行う
場合を例にとって説明してきたが、この発明に係る空燃
比検出技術はそれに限定されるものではなく、空燃比セ
ンサの実測値に基づいて空燃比を制御する技術に全て妥
当する。
In the first to fourth embodiments, a case is described in which a model describing the behavior of the exhaust system is set, and the air-fuel ratio control is performed using an observer for observing the internal state of the model. As described above, the air-fuel ratio detection technology according to the present invention is not limited to this, and is applicable to any technology for controlling the air-fuel ratio based on the actual measurement value of the air-fuel ratio sensor.

【0081】更に、運転状態を機関回転数と吸気圧力な
どから検出したが、これに限られるものでなく、その他
のパラメータを追加しても良い。また機関負荷を示すパ
ラメータは吸気圧力に限られるものではなく、吸入空気
量、スロットル開度などでも良い。また大気圧で修正し
ても良い。
Further, the operating state is detected from the engine speed and the intake pressure, but the present invention is not limited to this, and other parameters may be added. The parameter indicating the engine load is not limited to the intake pressure, but may be an intake air amount, a throttle opening, or the like. Alternatively, the correction may be made at atmospheric pressure.

【0082】更には、空燃比センサとして広域空燃比セ
ンサを使用する場合を例にとって説明したが、いわゆる
2 センサを用いて空燃比を制御する場合にも妥当す
る。
Further, the case where the wide-range air-fuel ratio sensor is used as the air-fuel ratio sensor has been described as an example, but the present invention is also applicable to the case where the air-fuel ratio is controlled using a so-called O 2 sensor.

【0083】[0083]

【発明の効果】請求項1項にあっては、オブザーバ行列
の推定精度が向上し、オブザーバの係数同定精度が向上
するので、各気筒空燃比の推定精度が向上する。
According to the first aspect of the present invention, the estimation accuracy of the observer matrix is improved, and the identification accuracy of the observer coefficient is improved. Therefore, the estimation accuracy of each cylinder air-fuel ratio is improved.

【0084】請求項2項にあっては、請求項1項で述べ
た効果に加え、単一のセンサを用いるのみで、各気筒の
空燃比フィードバック制御を精度良く行うことができ
る。
According to the second aspect, in addition to the effects described in the first aspect, the air-fuel ratio feedback control of each cylinder can be accurately performed only by using a single sensor.

【0085】請求項3項にあっては、請求項2項で述べ
た効果に加え、各気筒の空燃比は集合部空燃比に収束
し、集合部空燃比は目標空燃比に収束することとなっ
て、結果的に全ての気筒の空燃比が目標空燃比に収束
し、空燃比フィードバック制御の収束性を向上させるこ
とができる。
According to a third aspect of the present invention, in addition to the effects described in the second aspect, the air-fuel ratio of each cylinder converges to the collecting air-fuel ratio, and the collecting air-fuel ratio converges to the target air-fuel ratio. As a result, the air-fuel ratios of all the cylinders converge to the target air-fuel ratio, and the convergence of the air-fuel ratio feedback control can be improved.

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

【図1】この発明に係る内燃機関の気筒別空燃比推定装
置および制御装置を実現する、内燃機関の空燃比フィー
ドバック制御装置を全体的に示すブロック図である。
FIG. 1 is an overall block diagram showing an air-fuel ratio feedback control device for an internal combustion engine that implements a cylinder-by-cylinder air-fuel ratio estimation device and control device according to the present invention.

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

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

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

【図5】空燃比センサの検出挙動をモデル化した真の空
燃比推定器を示すブロック線図である。
FIG. 5 is a block diagram illustrating a true air-fuel ratio estimator that models the detection behavior of an air-fuel ratio sensor.

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

【図7】図6に示すモデルを用いて4気筒内燃機関につ
いて3気筒の空燃比を14.7に、1気筒の空燃比を1
2.0にして燃料を供給する場合を示すデータ図であ
る。
FIG. 7 shows a four-cylinder internal combustion engine having an air-fuel ratio of 14.7 and an air-fuel ratio of one cylinder of 1 using the model shown in FIG.
FIG. 4 is a data diagram showing a case where fuel is supplied at 2.0.

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

【図9】図7に示す入力を与えたときの図6モデルの集
合部の空燃比をLAFセンサの応答遅れを考慮して表し
たデータと、同じ場合のLAFセンサ出力の実測値を比
較するグラフ図である。
9 compares data representing the air-fuel ratio of the collective part of the model of FIG. 6 when the input shown in FIG. 7 is given in consideration of the response delay of the LAF sensor, and the measured value of the LAF sensor output in the same case. FIG.

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

【図11】先の出願で用いるオブザーバの構成を示すブ
ロック線図である。
FIG. 11 is a block diagram showing a configuration of an observer used in the earlier application.

【図12】図6に示すモデルと図11に示すオブザーバ
を組み合わせた構成を示す説明ブロック図である。
12 is an explanatory block diagram showing a configuration in which the model shown in FIG. 6 and the observer shown in FIG. 11 are combined.

【図13】この発明で予定する空燃比の気筒別フィード
バック制御を示すブロック図である。
FIG. 13 is a block diagram showing cylinder-by-cylinder feedback control of an air-fuel ratio planned in the present invention.

【図14】この発明の第1の実施の形態を示す図6に部
分的に類似するモデルのブロック図である。
FIG. 14 is a block diagram of a model partially similar to FIG. 6 showing the first embodiment of the present invention;

【図15】この発明の第2の実施の形態を示すオブザー
バ行列を複数個設けた構成を示す説明図である。
FIG. 15 is an explanatory diagram showing a configuration in which a plurality of observer matrices are provided according to the second embodiment of the present invention.

【図16】この発明の第3の実施の形態を示すフロー・
チャートである。
FIG. 16 is a flowchart showing a third embodiment of the present invention;
It is a chart.

【図17】第3の実施の形態で予定する可変バルブタイ
ミング機構を説明する、油圧回路図を含む部分断面図で
ある。
FIG. 17 is a partial cross-sectional view including a hydraulic circuit diagram for explaining a variable valve timing mechanism planned in a third embodiment.

【図18】バルブタイミングの切り換え特性を示す説明
図である。
FIG. 18 is an explanatory diagram showing switching characteristics of valve timing.

【図19】図16フロー・チャートで使用する空燃比セ
ンサ出力のサンプルタイミング用のマップの特性を示す
説明図である。
FIG. 19 is an explanatory diagram showing characteristics of a map for sampling timing of the output of the air-fuel ratio sensor used in the flow chart of FIG. 16;

【図20】この発明の第4の実施の形態を示す、図16
と同様のフロー・チャートである。
FIG. 20 shows a fourth embodiment of the present invention,
It is a flow chart similar to.

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

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

───────────────────────────────────────────────────── フロントページの続き (72)発明者 木村 英輔 埼玉県和光市中央1丁目4番1号 株式会 社本田技術研究所内 (72)発明者 赤崎 修介 埼玉県和光市中央1丁目4番1号 株式会 社本田技術研究所内 (72)発明者 西村 要一 埼玉県和光市中央1丁目4番1号 株式会 社本田技術研究所内 ──────────────────────────────────────────────────続 き Continuing from the front page (72) Inventor Eisuke Kimura 1-4-1 Chuo, Wako-shi, Saitama Prefecture Inside Honda R & D Co., Ltd. (72) Shusuke Akasaki 1-4-1 Chuo, Wako-shi, Saitama Inside Honda R & D Co., Ltd. (72) Inventor Yoichi Nishimura 1-4-1 Chuo, Wako-shi, Saitama Pref.

Claims (3)

【特許請求の範囲】[Claims] 【請求項1】 多気筒内燃機関の排気系集合部に空燃比
センサを配置してその出力から各気筒の入力混合気の空
燃比を推定する装置であって、 a.前記センサの出力値を各気筒の燃焼履歴に重み係数
Cnを乗じた加重平均値からなるものとみなして排気系
の挙動を記述するモデルを設定して各気筒の空燃比を状
態変数とする状態方程式を設定し、その内部状態を観測
するオブザーバを設定してその出力を求める第1の手
段、および b.前記オブザーバ出力から各気筒の空燃比を推定する
第2の手段、を有する内燃機関の気筒別空燃比推定装置
において、前記内燃機関で考慮すべき燃焼履歴の数を3
以上とするとき、直前の2つの燃焼履歴以上で、かつ前
記考慮すべき燃焼履歴の数より少ない燃焼履歴を除く、
それ以前の燃焼履歴についての加重平均の重み係数Cn
を零にすることを特徴とする内燃機関の気筒別空燃比推
定装置。
1. An apparatus for arranging an air-fuel ratio sensor in an exhaust system collecting portion of a multi-cylinder internal combustion engine and estimating an air-fuel ratio of an input air-fuel mixture of each cylinder from an output of the sensor. Assuming that the output value of the sensor is constituted by a weighted average value obtained by multiplying the combustion history of each cylinder by a weight coefficient Cn, a model is set which describes the behavior of the exhaust system, and the air-fuel ratio of each cylinder is set as a state variable. First means for setting an equation, setting an observer for observing its internal state and determining its output, and b. A second means for estimating the air-fuel ratio of each cylinder from the observer output, wherein the number of combustion histories to be considered in the internal combustion engine is 3
When the above, excluding combustion histories that are equal to or greater than the two immediately preceding combustion histories and less than the number of combustion histories to be considered,
Weight coefficient Cn of the weighted average of the combustion history before that
A cylinder-by-cylinder air-fuel ratio estimating device for an internal combustion engine.
【請求項2】a.多気筒内燃機関の排気系集合部に配置
された空燃比センサの出力値を各気筒の燃焼履歴に重み
係数Cnを乗じた加重平均値からなるものとみなすと共
に、前記内燃機関で考慮すべき燃焼履歴の数を3以上と
するとき、直前の2つの燃焼履歴以上で、かつ前記考慮
すべき燃焼履歴の数より少ない燃焼履歴を除く、それ以
前の燃焼履歴についての加重平均の重み係数Cnを零に
して排気系の挙動を記述するモデルを設定して各気筒の
空燃比を状態変数とする状態方程式を設定し、その内部
状態を観測するオブザーバを設定してその出力を求める
第1の手段、 b.前記オブザーバ出力から各気筒の空燃比#nA/F
を推定する第2の手段、 c.前記推定された各気筒の空燃比#nA/Fを目標値
にフィードバック制御する気筒毎フィードバック補正項
#nKLAFを演算する第3の手段、および d.前記気筒毎フィードバック補正項#nKLAFに基
づいて前記内燃機関の空燃比をフィードバック制御する
フィードバック制御手段、を備えたことを特徴とする内
燃機関の気筒別空燃比フィードバック制御装置。
2. a. The output value of the air-fuel ratio sensor arranged in the exhaust system collecting part of the multi-cylinder internal combustion engine is regarded as consisting of a weighted average value obtained by multiplying the combustion history of each cylinder by a weight coefficient Cn, and the combustion to be considered in the internal combustion engine is considered. When the number of histories is three or more, the weighting coefficient Cn of the weighted average of the previous combustion histories excluding the combustion histories that are equal to or more than the two immediately preceding combustion histories and less than the number of the combustion histories to be considered is set to zero. A first means for setting a model describing the behavior of the exhaust system, setting a state equation using the air-fuel ratio of each cylinder as a state variable, setting an observer for observing the internal state thereof, and obtaining the output thereof, b. From the observer output, the air-fuel ratio of each cylinder # nA / F
A second means for estimating c. Third means for calculating a cylinder-by-cylinder feedback correction term #nKLAF for feedback-controlling the estimated air-fuel ratio # nA / F of each cylinder to a target value; and d. Feedback control means for performing feedback control of the air-fuel ratio of the internal combustion engine based on the cylinder-by-cylinder feedback correction term #nKLAF.
【請求項3】 さらに、 e.前記空燃比センサの出力から前記排気系集合部の空
燃比を検出し、検出した排気系集合部の空燃比を目標値
に制御するための集合部フィードバック補正項KLAF
を演算する第4の手段、を備え、前記フィードバック制
御手段は、前記集合部フィードバック補正項KLAFと
気筒毎フィードバック補正項#nKLAFに基づいて前
記内燃機関の空燃比をフィードバック制御する如く構成
したことを特徴とする請求項2項記載の内燃機関の気筒
別空燃比フィードバック制御装置。
3. Further, e. An assembly feedback correction term KLAF for detecting the air-fuel ratio of the exhaust system assembly from the output of the air-fuel ratio sensor and controlling the detected air-fuel ratio of the exhaust system assembly to a target value.
The feedback control means is configured to feedback-control the air-fuel ratio of the internal combustion engine based on the collective feedback correction term KLAF and the cylinder-by-cylinder feedback correction term #nKLAF. The cylinder-by-cylinder air-fuel ratio feedback control device for an internal combustion engine according to claim 2, characterized in that:
JP2000254398A 2000-08-24 2000-08-24 Cylinder air-fuel ratio estimating device and control device for internal combustion engine Expired - Fee Related JP3361316B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2000254398A JP3361316B2 (en) 2000-08-24 2000-08-24 Cylinder air-fuel ratio estimating device and control device for internal combustion engine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2000254398A JP3361316B2 (en) 2000-08-24 2000-08-24 Cylinder air-fuel ratio estimating device and control device for internal combustion engine

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP12785897A Division JPH1047134A (en) 1997-05-02 1997-05-02 Individual cylinder air-fuel ratio estimating device of internal combustion engine and control device

Publications (2)

Publication Number Publication Date
JP2001073844A true JP2001073844A (en) 2001-03-21
JP3361316B2 JP3361316B2 (en) 2003-01-07

Family

ID=18743355

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2000254398A Expired - Fee Related JP3361316B2 (en) 2000-08-24 2000-08-24 Cylinder air-fuel ratio estimating device and control device for internal combustion engine

Country Status (1)

Country Link
JP (1) JP3361316B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006336645A (en) * 2005-05-30 2006-12-14 Inst Fr Petrole Method for estimating air-fuel ratio in cylinder of internal combustion engine using extended kalman filter
CN111677593A (en) * 2020-02-24 2020-09-18 山东交通学院 Air-fuel ratio control method for electric control gas engine

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006336645A (en) * 2005-05-30 2006-12-14 Inst Fr Petrole Method for estimating air-fuel ratio in cylinder of internal combustion engine using extended kalman filter
JP4703488B2 (en) * 2005-05-30 2011-06-15 イエフペ エネルジ ヌヴェル Method for estimating air-fuel ratio in cylinder of internal combustion engine using extended Kalman filter
CN111677593A (en) * 2020-02-24 2020-09-18 山东交通学院 Air-fuel ratio control method for electric control gas engine

Also Published As

Publication number Publication date
JP3361316B2 (en) 2003-01-07

Similar Documents

Publication Publication Date Title
EP0688945B1 (en) Air/fuel ratio detection system for multicylinder internal combustion engine
US5548514A (en) Air/fuel ratio estimation system for internal combustion engine
EP0553570B1 (en) Method for detecting and controlling air-fuel ratio in internal combustion engines
JP3162585B2 (en) Air-fuel ratio detection device for internal combustion engine
JP2717744B2 (en) Air-fuel ratio detection and control method for internal combustion engine
US9885306B2 (en) Control device for internal combustion engine
CA2572738A1 (en) Device, method, and program for estimating intake air amount
JPH06173755A (en) Air-fuel ratio estimating unit for each engine cylinder of internal combustion engine
JP2684011B2 (en) Internal combustion engine abnormality determination device
EP0643211B1 (en) Air-fuel ratio estimator for internal combustion engine
JP2689368B2 (en) Air-fuel ratio estimating device for each cylinder of internal combustion engine
JP2001073844A (en) Air-fuel ratio presuming device for each cylinder of internal-combustion engine and control device
JP5002987B2 (en) Engine cylinder intake gas amount measuring device
JP2689362B2 (en) Air-fuel ratio detection method for internal combustion engine
JPH1047134A (en) Individual cylinder air-fuel ratio estimating device of internal combustion engine and control device
JP2750673B2 (en) Air-fuel ratio detector for multi-cylinder internal combustion engine
JP3162567B2 (en) Air-fuel ratio estimator for each cylinder of internal combustion engine
JP2016211504A (en) Control device of internal combustion engine
JP2683974B2 (en) Air-fuel ratio control method for internal combustion engine
JP3683355B2 (en) Cylinder air-fuel ratio estimation device for internal combustion engine
JP3602575B2 (en) Air-fuel ratio estimator for each cylinder of internal combustion engine
JPH109038A (en) Air-fuel ratio detecting device of engine
JPH04365944A (en) Fuel injection amount control device for internal combustion engine
JPH07151584A (en) Suction wall face temperature inferring device and fuel injection amount control device of internal combustion engine
JPS6098143A (en) Air-fuel ratio control method of 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: 20020910

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

Free format text: PAYMENT UNTIL: 20071018

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20081018

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20081018

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20091018

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20091018

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20101018

Year of fee payment: 8

LAPS Cancellation because of no payment of annual fees