JP3880861B2 - 内燃機関の空燃比制御装置 - Google Patents
内燃機関の空燃比制御装置 Download PDFInfo
- Publication number
- JP3880861B2 JP3880861B2 JP2002012854A JP2002012854A JP3880861B2 JP 3880861 B2 JP3880861 B2 JP 3880861B2 JP 2002012854 A JP2002012854 A JP 2002012854A JP 2002012854 A JP2002012854 A JP 2002012854A JP 3880861 B2 JP3880861 B2 JP 3880861B2
- Authority
- JP
- Japan
- Prior art keywords
- fuel ratio
- air
- value
- target
- target air
- 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
Links
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/02—Circuit arrangements for generating control signals
- F02D41/14—Introducing closed-loop corrections
- F02D41/1438—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor
- F02D41/1439—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the position of the sensor
- F02D41/1441—Plural sensors
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01N—GAS-FLOW SILENCERS OR EXHAUST APPARATUS FOR MACHINES OR ENGINES IN GENERAL; GAS-FLOW SILENCERS OR EXHAUST APPARATUS FOR INTERNAL COMBUSTION ENGINES
- F01N13/00—Exhaust or silencing apparatus characterised by constructional features ; Exhaust or silencing apparatus, or parts thereof, having pertinent characteristics not provided for in, or of interest apart from, groups F01N1/00 - F01N5/00, F01N9/00, F01N11/00
- F01N13/009—Exhaust or silencing apparatus characterised by constructional features ; Exhaust or silencing apparatus, or parts thereof, having pertinent characteristics not provided for in, or of interest apart from, groups F01N1/00 - F01N5/00, F01N9/00, F01N11/00 having two or more separate purifying devices arranged in series
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/02—Circuit arrangements for generating control signals
- F02D41/14—Introducing closed-loop corrections
- F02D41/1401—Introducing closed-loop corrections characterised by the control or regulation method
- F02D41/1403—Sliding mode control
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/02—Circuit arrangements for generating control signals
- F02D41/14—Introducing closed-loop corrections
- F02D41/1401—Introducing closed-loop corrections characterised by the control or regulation method
- F02D41/1408—Dithering techniques
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/02—Circuit arrangements for generating control signals
- F02D41/14—Introducing closed-loop corrections
- F02D41/1438—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor
- F02D41/1493—Details
- F02D41/1495—Detection of abnormalities in the air/fuel ratio feedback system
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D2200/00—Input parameters for engine control
- F02D2200/02—Input parameters for engine control the parameters being related to the engine
- F02D2200/04—Engine intake system parameters
- F02D2200/0404—Throttle position
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D2200/00—Input parameters for engine control
- F02D2200/02—Input parameters for engine control the parameters being related to the engine
- F02D2200/04—Engine intake system parameters
- F02D2200/0406—Intake manifold pressure
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D2200/00—Input parameters for engine control
- F02D2200/02—Input parameters for engine control the parameters being related to the engine
- F02D2200/04—Engine intake system parameters
- F02D2200/0414—Air temperature
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D2200/00—Input parameters for engine control
- F02D2200/50—Input parameters for engine control said parameters being related to the vehicle or its components
- F02D2200/501—Vehicle speed
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D2200/00—Input parameters for engine control
- F02D2200/70—Input parameters for engine control said parameters being related to the vehicle exterior
- F02D2200/703—Atmospheric pressure
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/02—Circuit arrangements for generating control signals
- F02D41/14—Introducing closed-loop corrections
- F02D41/1438—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor
- F02D41/1444—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the characteristics of the combustion gases
- F02D41/1454—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the characteristics of the combustion gases the characteristics being an oxygen content or concentration or the air-fuel ratio
- F02D41/1456—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the characteristics of the combustion gases the characteristics being an oxygen content or concentration or the air-fuel ratio with sensor output signal being linear or quasi-linear with the concentration of oxygen
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Electrical Control Of Air Or Fuel Supplied To Internal-Combustion Engine (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Description
【発明の属する技術分野】
本発明は、目標空燃比を所定の振幅で周期的に変動(振動)させる制御、すなわちパータベーション制御を実行する内燃機関の空燃比制御装置に関する。
【0002】
【従来の技術】
従来、この種の空燃比制御装置として、例えば特開昭64−66441号公報に記載されたものが知られている。この空燃比制御装置は、排気管の触媒よりも下流側に設けられ、排気ガス中の酸素濃度に応じた検出信号を出力するO2センサを備えている。この空燃比制御装置は、O2センサの検出信号の平均値を算出するとともに、この平均値に応じて、パータベーション制御の基準値を算出する。そして、この基準値に所定の振幅値を加算し、または加算しないことにより、空燃比補正係数が算出され、その結果、空燃比補正係数すなわち空燃比が、所定の振幅で矩形波状に繰り返し変動するようにパータベーション制御される。また、そのパータベーション制御における空燃比の変動周波数は、1〜4Hzの範囲内の値に設定される。
【0003】
【発明が解決しようとする課題】
以上のように、この空燃比制御装置では、パータベーション制御の空燃比の変動周波数が、1〜4Hzの範囲内の値に設定されるけれども、この種のパータベーション制御を実行した場合、触媒の排気ガス浄化率は、必ずしも一定ではない。すなわち、触媒が劣化していないときには、パータベーション制御の空燃比の変動周波数にかかわらず、良好な値を示す一方、触媒が劣化したときには、3Hzよりも低い周波数域で大幅に悪化するとともに、3Hz以上、より好ましくは5Hz以上の周波数域で良好な値を示すことが確認された(図38参照)。したがって、上記従来の空燃比制御装置によれば、触媒の劣化の有無にかかわらず、空燃比の変動周波数が1〜4Hzの範囲内の値に設定されるため、触媒が劣化したときに、触媒による排気ガス浄化率が悪化し、その結果、触媒で浄化された排気ガスの特性(以下「触媒後排気ガス特性」という)が悪化してしまうおそれがある。
【0004】
本発明は、上記課題を解決するためになされたもので、パータベーション制御を行う場合において、触媒の劣化の有無にかかわらず、排気ガス浄化率を良好に保持することができ、それにより、触媒後排気ガス特性を向上させることができる内燃機関の空燃比制御装置を提供することを目的とする。
【0005】
【課題を解決するための手段】
この目的を達成するために、請求項1に係る内燃機関3の空燃比制御装置1は、内燃機関3の排気通路(排気管7)の触媒(第1触媒装置8a)よりも上流側における排気ガスの空燃比を表す検出信号KACTを出力する上流側空燃比センサ(LAFセンサ14)と、排気通路(排気管7)の触媒(第1触媒装置8a)よりも下流側における排気ガスの空燃比を表す検出信号Voutを出力する下流側空燃比センサ(酸素濃度センサ15)と、下流側空燃比センサの出力Voutを所定の目標値Vopに収束させるための目標空燃比KCMDを、所定の振幅(ゲインKDSM)で、かつ下流側空燃比センサの出力Voutが所定の目標値Vop付近にあるときにはそれ以外のときよりも高い所定の周波数で変動するように、設定する目標空燃比設定手段(ECU2、ADSMコントローラ20、スケジュール型DSMコントローラ20A)と、上流側空燃比センサの出力KACTに基づき、触媒(第1触媒装置8a)よりも上流側における排気ガスの空燃比を目標空燃比設定手段により設定された目標空燃比KCMDに一致させるように、内燃機関3に供給される混合気の空燃比を制御する空燃比制御手段(ECU2)と、を備えることを特徴とする。
【0006】
この内燃機関の空燃比制御装置によれば、下流側空燃比センサの出力を所定の目標値に収束させるための目標空燃比が、所定の振幅でかつ下流側空燃比センサの出力が所定の目標値付近にあるときには、それ以外のときよりも高い所定の周波数で変動するように設定され、さらに、触媒よりも上流側における排気ガスの空燃比をそのように設定された目標空燃比に一致させるように、内燃機関に供給される混合気の空燃比が制御される。すなわち、パータベーション制御が実行される。したがって、例えば、所定の周波数を、触媒の良好な排気ガス浄化率が得られる、前述した3Hz以上、より好ましくは5Hz以上の周波数に設定することにより、触媒の劣化の有無にかかわらず、排気ガス浄化率を良好な状態に保持することが可能になる。これに加えて、所定の目標値を、触媒による良好な排気ガス浄化率が得られるような値(例えば図2の目標値Vop)に設定することにより、下流側空燃比センサの出力が目標値付近にあるときに、目標空燃比を上記の周波数で変動させることによって、排気ガス浄化率をさらに向上させることができる。
【0007】
一方、下流側空燃比センサの出力が所定の目標値付近から離れていて、内燃機関に供給される混合気の空燃比が排気ガス浄化率の低下を招く状態にあるときには、目標空燃比を上記3Hzよりも低い周波数で変動させることにより、下流側空燃比センサの出力を迅速に所定の目標値付近に近づけることができ、排気ガス浄化率を良好な状態に迅速に復帰させることができる。以上のように、触媒の排気ガス浄化率を良好な状態に保持することができ、触媒後排気ガス特性を向上させることができる。
【0008】
請求項2に係る発明は、請求項1に記載の内燃機関3の空燃比制御装置1において、目標空燃比設定手段(ECU2、ADSMコントローラ20、スケジュール型DSMコントローラ20A)は、Δ変調アルゴリズム、ΔΣ変調アルゴリズムおよびΣΔ変調アルゴリズムのうちのいずれか1つの変調アルゴリズムに基づき、目標空燃比KCMDを設定することを特徴とする。
【0009】
この内燃機関の空燃比制御装置によれば、目標空燃比が、Δ変調アルゴリズム、ΔΣ変調アルゴリズムおよびΣΔ変調アルゴリズムのうちのいずれか1つの変調アルゴリズムに基づいて設定される。一般に、この種の各変調アルゴリズムは、出力が値1および値−1の一方に変化する特性、すなわち出力の正負が変化する特性を有している。この出力の正負は、Δ変調アルゴリズムでは、入力と出力の積分値との偏差の正負と同じになるように決定され、ΔΣ変調アルゴリズムでは、入力と出力との偏差を積分した値の正負と同じになるように決定され、ΣΔ変調アルゴリズムでは、入力の積分値と出力の積分値との偏差の正負と同じになるように決定される。したがって、このような各変調アルゴリズムの特性を利用することにより、下流側空燃比センサの出力が所定の目標値付近に近づくと、触媒の劣化状態および内燃機関の運転状態などとは無関係に、目標空燃比を、その変動周波数がより高い値に自動的に変化するように制御することができる。これにより、目標空燃比と所定の目標値との比較結果に基づいて目標空燃比の変動周波数を切り換えるようなプログラムなどを追加することなく、排気ガス浄化率を良好かつ安定した状態に保持できるように、空燃比を制御することができる。
【0010】
請求項3に係る発明は、請求項2に記載の内燃機関3の空燃比制御装置1において、目標空燃比設定手段(ECU2、ADSMコントローラ20)は、予測アルゴリズム(式(7))に基づき、下流側空燃比センサの出力を表す値(出力偏差VO2)の予測値PREVO2を算出する予測値算出手段(ECU2、状態予測器22)と、1つの変調アルゴリズムに基づき、算出された予測値PREVO2に応じて、目標空燃比KCMDを算出する目標空燃比算出手段(ECU2、DSMコントローラ21)と、を備えることを特徴とする。
【0011】
この内燃機関の空燃比制御装置によれば、空燃比センサの出力を表す値の予測値が、予測アルゴリズムに基づいて算出され、目標空燃比が、1つの変調アルゴリズムに基づき、算出された予測値に応じて算出される。目標空燃比が以上のように算出されるので、予測値を、例えば空燃比制御における制御対象の動特性、すなわち内燃機関に供給される混合気と、下流側空燃比センサの出力との間の位相遅れやむだ時間などを反映させた値として算出することにより、空燃比制御における入出力間の制御タイミングのずれを解消することが可能になる。その結果、排気ガス浄化率をより一層、向上させること、および空燃比制御の安定性を確保することが可能になる。
【0012】
請求項4に係る発明は、請求項3に記載の内燃機関3の空燃比制御装置1において、予測アルゴリズムは、下流側空燃比センサの出力を表す値(出力偏差VO2)と、目標空燃比を表す値(空燃比偏差DKCMD)および上流側空燃比センサの出力を表す値(LAF出力偏差DKACT)の一方とを変数とする制御対象モデル(式(1))を適用したアルゴリズムであることを特徴とする。
【0013】
この内燃機関の空燃比制御装置によれば、空燃比センサの出力を表す値の予測値が、制御対象モデルを適用した予測アルゴリズムに基づいて算出されるので、この制御対象モデルを、例えば制御対象の入出力間の動特性、すなわち内燃機関に供給される混合気と、下流側空燃比センサの出力との間の位相遅れやむだ時間などを反映させたものとして定義することにより、予測値を、制御対象の動特性が適切に反映された値として算出することができ、空燃比制御における入出力間の制御タイミングのずれを適切に解消することができる。
【0014】
請求項5に係る発明は、請求項2に記載の内燃機関3の空燃比制御装置1において、目標空燃比設定手段(ECU2、ADSMコントローラ20)は、目標空燃比KCMDを表す値の時系列データ(DKCMD(k−dt))および上流側空燃比センサの出力を表す値の時系列データの一方、および下流側空燃比センサの出力Voutを表す値の時系列データ(VO2(k−1),VO2(k−2))を変数とする離散時間系の制御対象モデル(式(1))と、1つの変調アルゴリズムとに基づき、目標空燃比KCMDを算出する目標空燃比算出手段(ECU2、DSMコントローラ21、状態予測器22)と、離散時間系の制御対象モデルのモデルパラメータa1,a2,b1を逐次同定する同定手段(ECU2、オンボード同定器23)と、を備えることを特徴とする。
【0015】
この内燃機関の空燃比制御装置によれば、制御対象モデルのモデルパラメータが、逐次同定される、すなわちリアルタイムで同定されるとともに、そのようにモデルパラメータが同定される制御対象モデルと、1つの変調アルゴリズムとに基づいて、目標空燃比が設定される。それにより、制御対象の動特性、すなわち内燃機関に供給される混合気と、下流側空燃比センサの出力との間の位相遅れやむだ時間などが、内燃機関の運転状態の変化などに伴って変化したり、経年変化したりしている場合でも、それらの影響を回避しながら、制御対象モデルの動特性を制御対象の実際の動特性に適合させることができる。その結果、制御対象の動特性、例えば位相遅れやむだ時間などに起因する、空燃比制御における入出力間の制御タイミングのずれを適切に補正することができる。
【0016】
請求項6に係る発明は、請求項5に記載の内燃機関3の空燃比制御装置1において、目標空燃比算出手段(ECU2、DSMコントローラ21、状態予測器22)は、制御対象モデルを適用した予測アルゴリズムに基づき、下流側空燃比センサの出力を表す値の予測値PREVO2を算出するとともに、算出された予測値PREVO2および1つの変調アルゴリズムに基づいて、目標空燃比KCMDを算出することを特徴とする。
【0017】
この内燃機関の空燃比制御装置によれば、下流側空燃比センサの出力を表す値の予測値が、制御対象モデルを適用した予測アルゴリズムに基づいて算出されるとともに、目標空燃比が、算出された予測値および1つの変調アルゴリズムに基づいて算出される。この場合、上述したように同定されたモデルパラメータを用いることによって、制御対象モデルの動特性を実際の制御対象の動特性に適合させることができるので、そのような制御対象モデルを適用した予測アルゴリズムに基づいて算出することにより、予測値を制御対象の実際の動特性を反映させた値として算出することができる。その結果、空燃比制御における入出力間の制御タイミングのずれをより精度よく補正することができる。
【0018】
請求項7に係る発明は、請求項2に記載の内燃機関3の空燃比制御装置1において、内燃機関の運転状態を表す運転状態パラメータ(エンジン回転数NE、吸気管内絶対圧PBA)を検出する運転状態パラメータ検出手段(ECU2、吸気管内絶対圧センサ11、クランク角センサ13)をさらに備え、目標空燃比設定手段(ECU2、スケジュール型DSMコントローラ20A)は、下流側空燃比センサの出力を表す値(出力偏差VO2)と、目標空燃比を表す値(空燃比偏差DKCMD)および上流側空燃比センサの出力を表す値(LAF出力偏差DKACT)の一方とを変数とする制御対象モデル(式(1))、ならびに1つの変調アルゴリズムに基づき、目標空燃比KCMDを算出する目標空燃比算出手段(ECU2、状態予測器22、DSMコントローラ24)と、検出された運転状態パラメータ(エンジン回転数NE、吸気管内絶対圧PBA)に応じて、制御対象モデルのモデルパラメータa1,a2,b1を設定するモデルパラメータ設定手段(ECU2、パラメータスケジューラ28)と、を備えることを特徴とする。
【0019】
この内燃機関の空燃比制御装置によれば、目標空燃比が、下流側空燃比センサの出力を表す値と、目標空燃比を表す値および上流側空燃比センサの出力を表す値の一方とを変数とする制御対象モデル、ならびに1つの変調アルゴリズムに基づいて算出されるとともに、制御対象モデルのモデルパラメータが、検出された運転状態パラメータに応じて設定されるので、制御対象モデルの動特性を、実際の制御対象の動特性に迅速に適合させることができる。その結果、空燃比制御における入出力間の制御タイミングのずれを、迅速かつ適切に補正することができる。
【0020】
請求項8に係る発明は、請求項7に記載の内燃機関3の空燃比制御装置1において、目標空燃比算出手段(ECU2、スケジュール型DSMコントローラ20A)は、制御対象モデル(式(1))を適用した予測アルゴリズム(式(7))に基づき、下流側空燃比センサの出力を表す値(出力偏差VO2)の予測値PREVO2を算出するとともに、算出された予測値PREVO2および1つの変調アルゴリズムに基づいて、目標空燃比KCMDを算出することを特徴とする。
【0021】
この内燃機関の空燃比制御装置によれば、下流側空燃比センサの出力を表す値の予測値が、制御対象モデルを適用した予測アルゴリズムに基づいて算出されるとともに、目標空燃比が、算出された予測値および1つの変調アルゴリズムに基づいて算出される。この場合、上述したように設定されたモデルパラメータを用いることによって、制御対象モデルの動特性を実際の制御対象の動特性に迅速に適合させることができるので、そのような制御対象モデルを適用した予測アルゴリズムに基づいて算出することにより、予測値を制御対象の実際の動特性を反映させた値として迅速に算出することができる。その結果、空燃比制御における入出力間の制御タイミングのずれをより精度よく迅速に補正することができる。
【0022】
請求項9に係る発明は、請求項1ないし8のいずれかにおいて、内燃機関の負荷を表すパラメータ(排気ガスボリュームAB_SV)を検出する負荷パラメータ検出手段(ECU2、吸気管内絶対圧センサ11、クランク角センサ13、ステップ80)をさらに備え、目標空燃比設定手段(ECU2)は、検出された負荷パラメータ(排気ガスボリュームAB_SV)に応じて、所定の振幅(ゲインKDSM)を設定する(ステップ200)ことを特徴とする。
【0023】
この内燃機関の空燃比制御装置によれば、内燃機関の負荷を表す負荷パラメータに応じて、目標空燃比の変動の振幅が設定されるので、下流側空燃比センサの出力の応答性が負荷の変化に伴って変化するのを補償しながら、目標空燃比の振幅を設定することができる。これにより、目標空燃比の振幅を、内燃機関の負荷の変化に伴うオーバーゲイン状態などを回避しながら、適切に設定することができ、その結果、良好な排気ガス浄化率を確保することができる。
【0024】
請求項10に係る内燃機関3の空燃比制御装置1は、内燃機関3の排気通路(排気管7)の触媒(第2触媒装置8b)よりも下流側における排気ガスの空燃比を表す検出信号Voutを出力する空燃比センサ(酸素濃度センサ15)と、空燃比センサの出力を所定の目標値に収束させるための目標空燃比KCMDを、所定の振幅(ゲインKDSM)で、かつ空燃比センサの出力Voutが所定の目標値Vop付近にあるときにはそれ以外のときよりも高い所定の周波数で変動するように、設定する目標空燃比設定手段(ECU2、ADSMコントローラ20、スケジュール型DSMコントローラ20A)と、設定された目標空燃比KCMDに応じて、内燃機関3に供給される混合気の空燃比を制御する空燃比制御手段(ECU2)と、を備えることを特徴とする。
【0025】
この内燃機関の空燃比制御装置によれば、空燃比センサの出力を所定の目標値に収束させるための目標空燃比が、所定の振幅でかつ空燃比センサの出力が所定の目標値付近にあるときには、それ以外のときよりも高い所定の周波数で変動するように設定され、さらに、そのように設定された目標空燃比に応じて、内燃機関に供給される混合気の空燃比が制御される。すなわち、パータベーション制御が実行される。したがって、例えば、所定の周波数を、触媒の良好な排気ガス浄化率が得られる、前述した3Hz以上、より好ましくは5Hz以上の周波数に設定することにより、触媒の劣化の有無にかかわらず、排気ガス浄化率を良好な状態に保持することが可能になる。これに加えて、所定の目標値を、触媒による良好な排気ガス浄化率が得られるような値(例えば図2の目標値Vop)に設定することにより、空燃比センサの出力が目標値付近にあるときに、目標空燃比を上記の周波数で変動させることにより、排気ガス浄化率をさらに向上させることができる。
【0026】
一方、空燃比センサの出力が所定の目標値付近から離れていて、内燃機関に供給される混合気の空燃比が排気ガス浄化率の低下を招く状態にあるときには、目標空燃比を上記3Hzよりも低い周波数で変動させることにより、言い換えれば目標空燃比をあまり変化させないことにより、空燃比センサの出力を迅速に所定の目標値付近に近づけることができ、排気ガス浄化率を良好な状態に迅速に復帰させることができる。以上のように、触媒の排気ガス浄化率を良好な状態に保持することができ、触媒後排気ガス特性を向上させることができる。また、単一の空燃比センサを用いるだけでよいので、このような空燃比制御装置を比較的、安価に実現することができる。
【0027】
以下に述べる請求項11ないし18に係る発明はそれぞれ、前述した請求項2ないし9に係る発明とほぼ同様に構成されている。
【0028】
すなわち、請求項11に係る発明は、請求項10に記載の内燃機関3の空燃比制御装置1において、目標空燃比設定手段(ECU2、ADSMコントローラ20、スケジュール型DSMコントローラ20A)は、Δ変調アルゴリズム、ΔΣ変調アルゴリズムおよびΣΔ変調アルゴリズムのうちのいずれか1つの変調アルゴリズムに基づき、目標空燃比KCMDを設定することを特徴とする。
【0029】
また、請求項12に係る発明は、請求項11に記載の内燃機関3の空燃比制御装置1において、目標空燃比設定手段(ECU2、ADSMコントローラ20)は、予測アルゴリズム(式(6))に基づき、空燃比センサの出力を表す値(出力偏差VO2)の予測値PREVO2を算出する予測値算出手段(ECU2、状態予測器22)と、1つの変調アルゴリズムに基づき、算出された予測値PREVO2に応じて、目標空燃比KCMDを算出する目標空燃比算出手段(ECU2、DSMコントローラ21)と、を備えることを特徴とする。
【0030】
さらに、請求項13に係る発明は、請求項12に記載の内燃機関3の空燃比制御装置1において、予測アルゴリズムは、空燃比センサの出力を表す値(出力偏差VO2)と目標空燃比を表す値(空燃比偏差DKCMD)とを変数とする制御対象モデル(式(1))を適用したアルゴリズムであることを特徴とする。
【0031】
また、請求項14に係る発明は、請求項11に記載の内燃機関3の空燃比制御装置1において、目標空燃比設定手段(ECU2、ADSMコントローラ20)は、目標空燃比を表す値の時系列データ(DKCMD(k−dt))および空燃比センサの出力を表す値の時系列データ(VO2(k−1),VO2(k−2))を変数とする離散時間系の制御対象モデル(式(1))と、1つの変調アルゴリズムとに基づき、目標空燃比KCMDを算出する目標空燃比算出手段(ECU2、DSMコントローラ21、状態予測器22)と、離散時間系の制御対象モデルのモデルパラメータa1,a2,b1を逐次同定する同定手段(オンボード同定器23)と、を備えることを特徴とする。
【0032】
さらに、請求項15に係る発明は、請求項14に記載の内燃機関3の空燃比制御装置1において、目標空燃比算出手段(ECU2、DSMコントローラ21、状態予測器22)は、制御対象モデルを適用した予測アルゴリズムに基づき、空燃比センサの出力を表す値の予測値PREVO2を算出するとともに、算出された予測値PREVO2および1つの変調アルゴリズムに基づいて、目標空燃比KCMDを算出することを特徴とする。
【0033】
また、請求項16に係る発明は、請求項11に記載の内燃機関3の空燃比制御装置1において、内燃機関の運転状態を表す運転状態パラメータ(エンジン回転数NE、吸気管内絶対圧PBA)を検出する運転状態パラメータ検出手段(ECU2、吸気管内絶対圧センサ11、クランク角センサ13)をさらに備え、目標空燃比設定手段(ECU2、スケジュール型DSMコントローラ20A)は、空燃比センサの出力を表す値(出力偏差VO2)および目標空燃比を表す値(空燃比偏差DKCMD)を変数とする制御対象モデル(式(1))と、1つの変調アルゴリズムとに基づき、目標空燃比KCMDを算出する目標空燃比算出手段(ECU2、状態予測器22、DSMコントローラ24)と、検出された運転状態パラメータ(エンジン回転数NE、吸気管内絶対圧PBA)に応じて、制御対象モデルのモデルパラメータa1,a2,b1を設定するモデルパラメータ設定手段(ECU2、パラメータスケジューラ28)と、を備えることを特徴とする。
【0034】
さらに、請求項17に係る発明は、請求項16に記載の内燃機関3の空燃比制御装置1において、目標空燃比算出手段(ECU2、スケジュール型DSMコントローラ20A)は、制御対象モデル(式(1))を適用した予測アルゴリズム(式(6))に基づき、空燃比センサの出力を表す値(出力偏差VO2)の予測値PREVO2を算出するとともに、算出された予測値PREVO2および1つの変調アルゴリズムに基づいて、目標空燃比KCMDを算出することを特徴とする。
【0035】
また、請求項18に係る発明は、請求項10ないし17のいずれかに記載の内燃機関3の空燃比制御装置1において、内燃機関の負荷を表す負荷パラメータ(排気ガスボリュームAB_SV)を検出する負荷パラメータ検出手段(ECU2、吸気管内絶対圧センサ11、クランク角センサ13)をさらに備え、目標空燃比設定手段(ECU2)は、検出された負荷パラメータ(排気ガスボリュームAB_SV)に応じて、所定の振幅を設定することを特徴とする。
【0036】
以上の請求項11ないし18に係る内燃機関の空燃比制御装置によれば、請求項2ないし9に係る発明と同様の作用効果をそれぞれ得ることができる。
【0037】
【発明の実施の形態】
以下、図面を参照しながら、本発明の第1実施形態に係る内燃機関の空燃比制御装置について説明する。図1は、本実施形態の空燃比制御装置1およびこれを適用した内燃機関3の概略構成を示している。同図に示すように、この空燃比制御装置1は、ECU2を備えており、このECU2は、後述するように、内燃機関(以下「エンジン」という)3の運転状態に応じて、これに供給する混合気の空燃比を制御する。
【0038】
このエンジン3は、図示しない車両に搭載された直列4気筒型ガソリンエンジンであり、第1〜第4の4つの気筒#1〜#4を備えている。このエンジン3の吸気管4のスロットル弁5の近傍には、例えばポテンショメータなどで構成されたスロットル弁開度センサ10が設けられている。このスロットル弁開度センサ10は、スロットル弁5の開度(以下「スロットル弁開度」という)θTHを検出して、その検出信号をECU2に送る。
【0039】
さらに、吸気管4のスロットル弁5よりも下流側には、吸気管内絶対圧センサ11が設けられている。この吸気管内絶対圧センサ11(運転状態パラメータ検出手段、負荷パラメータ検出手段)は、例えば半導体圧力センサなどで構成され、吸気管4内の吸気管内絶対圧PBAを検出し、その検出信号をECU2に出力する。
【0040】
また、吸気管4は、インテークマニホールド4aの4つの分岐部4bを介して4つの気筒#1〜#4にそれぞれ接続されている。各分岐部4bには、各気筒の図示しない吸気ポートの上流側に、インジェクタ6が取り付けられている。各インジェクタ6は、エンジン3の運転時に、ECU2からの駆動信号によって、その開弁時間である最終燃料噴射量TOUTおよび噴射タイミングが制御される。
【0041】
一方、エンジン3の本体には、例えばサーミスタなどで構成された水温センサ12が取り付けられている。水温センサ12は、エンジン3のシリンダブロック内を循環する冷却水の温度であるエンジン水温TWを検出し、その検出信号をECU2に出力する。
【0042】
また、エンジン3のクランクシャフト(図示せず)には、クランク角センサ13が設けられている。このクランク角センサ13(運転状態パラメータ検出手段、負荷パラメータ検出手段)は、クランクシャフトの回転に伴い、いずれもパルス信号であるCRK信号およびTDC信号をECU2に出力する。
【0043】
CRK信号は、所定のクランク角(例えば30゜)ごとに1パルスが出力される。ECU2は、このCRK信号に応じ、エンジン3の回転数(以下「エンジン回転数」という)NEを算出する。また、TDC信号は、各気筒のピストン(図示せず)が吸気行程のTDC位置よりも若干、手前の所定のクランク角位置にあることを表す信号であり、所定クランク角ごとに1パルスが出力される。
【0044】
一方、排気管7(排気通路)のエキゾーストマニホールド7aよりも下流側には、上流側から順に第1および第2の触媒装置8a,8b(触媒)が間隔を存して設けられている。各触媒装置8は、NOx触媒と3元触媒を組み合わせたものであり、このNOx触媒は、図示しないが、イリジウム触媒(イリジウムを担持した炭化ケイ素ウイスカ粉末とシリカの焼成体)をハニカム構造の基材の表面に被覆し、その上にペロブスカイト型複酸化物(LaCoO3粉末とシリカの焼成体)をさらに被覆したものである。触媒装置8は、NOx触媒による酸化還元作用により、リーンバーン運転時の排気ガス中のNOxを浄化するとともに、3元触媒の酸化還元作用により、リーンバーン運転以外の運転時の排気ガス中のCO、HCおよびNOxを浄化する。なお、触媒装置8は、NOx触媒と3元触媒を組み合わせたものに限らず、排気ガス中のCO、HCおよびNOxを浄化できるものであればよい。例えば、触媒装置8を、ペロブスカイト型触媒などの非金属触媒および/または3元触媒などの金属触媒で構成してもよい。
【0045】
これらの第1および第2触媒装置8a,8bの間には、酸素濃度センサ(以下「O2センサ」という)15が取り付けられている。このO2センサ15(下流側空燃比センサ、空燃比センサ)は、ジルコニアおよび白金電極などで構成され、第1触媒装置8aの下流側の排気ガス中の酸素濃度に基づく出力VoutをECU2に送る。このO2センサ15の出力Voutは、理論空燃比よりもリッチな混合気が燃焼したときには、ハイレベルの電圧値(例えば0.8V)となり、混合気がリーンのときには、ローレベルの電圧値(例えば0.2V)となるとともに、混合気が理論空燃比付近のときには、ハイレベルとローレベルの間の所定の目標値Vop(例えば0.6V)となる(図2参照)。
【0046】
また、第1触媒装置8aよりも上流側のエキゾーストマニホールド7aの集合部付近には、LAFセンサ14(上流側空燃比センサ)が取り付けられている。このLAFセンサ14は、O2センサ15と同様のセンサとリニアライザなどの検出回路とを組み合わせることによって構成されており、リッチ領域からリーン領域までの広範囲な空燃比の領域において排気ガス中の酸素濃度をリニアに検出し、その酸素濃度に比例する出力KACTをECU2に送る。この出力KACTは、空燃比の逆数に比例する当量比として表される。
【0047】
次に、図2を参照しながら、第1触媒装置8aの排気ガスの浄化率とO2センサ15の出力Vout(電圧値)との関係について説明する。同図は、第1触媒装置8aが、長時間の使用により浄化能力が低下した劣化状態と、浄化能力の高い未劣化状態の場合において、LAFセンサ14の出力KACTすなわちエンジン3に供給される混合気の空燃比が理論空燃比の付近で変化したときの、2つの第1触媒装置8aのHCおよびNOxの浄化率と、O2センサ15の出力Voutをそれぞれ測定した結果の一例を示している。同図において、破線で示すデータはいずれも、第1触媒装置8aが未劣化状態の場合の測定結果であり、実線で示すデータはいずれも、第1触媒装置8aが劣化状態の場合の測定結果である。また、LAFセンサ14の出力KACTが大きいほど、混合気の空燃比がよりリッチ側であることを示している。
【0048】
同図に示すように、第1触媒装置8aが劣化している場合には、未劣化状態の場合と比べて、排気ガスの浄化能力が低下していることにより、LAFセンサ14の出力KACTがよりリーン側の値KACT1のときに、O2センサ15の出力Voutが目標値Vopを横切っている。一方、第1触媒装置8aは、その劣化・未劣化状態にかかわらず、O2センサ15の出力Voutが目標値Vopにあるときに、HCおよびNOxを最も効率よく浄化する特性を有している。したがって、O2センサ15の出力Voutが目標値Vopになるように、混合気の空燃比を制御することにより、第1触媒装置8aによって排気ガスを最も効率よく浄化できることが判る。このため、後述する空燃比制御では、O2センサ15の出力Voutが目標値Vopに収束するように、目標空燃比KCMDが制御される。
【0049】
さらに、ECU2には、アクセル開度センサ16、大気圧センサ17、吸気温センサ18および車速センサ19などが接続されている。このアクセル開度センサ16は、車両の図示しないアクセルペダルの踏み込み量(以下「アクセル開度」という)APを検出し、その検出信号をECU2に出力する。また、大気圧センサ17、吸気温センサ18および車速センサ19はそれぞれ、大気圧PA、吸気温TAおよび車速VPを検出し、その検出信号をECU2に出力する。
【0050】
次に、ECU2(目標空燃比設定手段、空燃比制御手段、予測値算出手段、目標空燃比算出手段、同定手段、運転状態パラメータ検出手段、モデルパラメータ設定手段、負荷パラメータ検出手段)について説明する。
【0051】
このECU2は、I/Oインターフェース、CPU、RAMおよびROMなどからなるマイクロコンピュータから構成されており、前述した各種のセンサ10〜19の出力に応じて、エンジン3の運転状態を判別するとともに、ROMに予め記憶された制御プログラムやRAMに記憶されたデータなどに従って、後述する適応空燃比制御処理またはマップ検索処理を実行することにより、目標空燃比KCMDを算出する(なお、この目標空燃比KCMDは、空燃比の逆数に比例する当量比として算出される)。さらに、後述するように、この目標空燃比KCMDに基づいて、インジェクタ6の最終燃料噴射量TOUTを気筒ごとに算出し、この算出した最終燃料噴射量TOUTに基づいた駆動信号で、インジェクタ6を駆動することにより、LAFセンサ14の出力KACTを目標空燃比KCMDに一致させるように、混合気の空燃比をフィードバック制御する。
【0052】
図3に示すように、空燃比制御装置1は、目標空燃比KCMDを算出するADSMコントローラ20およびPRISMコントローラ21を備えており、両コントローラ20,21はいずれも、具体的には、ECU2により構成されている。
【0053】
以下、ADSMコントローラ20(目標空燃比設定手段)について説明する。このADSMコントローラ20は、以下に述べる適応予測型ΔΣ変調制御(Adaptive prediction Delta Sigma Modulation Control:以下「ADSM」という)処理の制御アルゴリズムにより、O2センサ15の出力Voutを目標値Vopに収束させるための目標空燃比KCMDを算出するものであり、状態予測器22、オンボード同定器23およびDSMコントローラ24により構成されている。なお、このADSM処理の具体的なプログラムについては、後述する。
【0054】
まず、状態予測器22(予測値算出手段、目標空燃比算出手段)について説明する。この状態予測器22は、以下に述べる予測アルゴリズムにより、出力偏差VO2の予測値PREVO2を予測(算出)するものである。本実施形態では、制御対象への制御入力を混合気の目標空燃比KCMDとし、制御対象の出力をO2センサ15の出力Voutとし、インジェクタ6を含むエンジン3の吸気系から、第1触媒装置8aを含む排気系の第1触媒装置8aの下流側のO2センサ15までの系を、制御対象と見なすとともに、この制御対象を、下式(1)に示すように、離散時間系モデルであるARXモデル(auto-regressive model with exogeneous input:外部入力を持つ自己回帰モデル)としてモデル化する。
【0055】
VO2(k)=a1・VO2(k-1)+a2・VO2(k-2)+b1・DKCMD(k-dt) ……(1)
ここで、VO2は、O2センサ15の出力Voutと前述した目標値Vopとの偏差(Vout−Vop)である出力偏差を表し、DKCMDは、目標空燃比KCMD(=φop)と基準値FLAFBASEとの偏差(KCMD−FLAFBASE)である空燃比偏差を表し、記号kは、各データのサンプリングサイクルの順番を表している。この基準値FLAFBASEは、所定の一定値に設定される。また、a1,a2,b1はモデルパラメータを表しており、オンボード同定器23により、後述するように逐次同定される。
【0056】
さらに、上記式(1)のdtは、目標空燃比KCMDの混合気がインジェクタ6により吸気系に供給されてから、O2センサ15の出力Voutに反映されるまでの予測時間を表しており、下式(2)のように定義される。
dt=d+d'+dd ……(2)
ここで、dは、LAFセンサ14からO2センサ15までの排気系のむだ時間を、d'は、インジェクタ6からLAFセンサ14までの空燃比操作系のむだ時間を、ddは、排気系と空燃比操作系との間の位相遅れ時間をそれぞれ表している(なお、後述する適応空燃比制御処理の制御プログラムでは、ADSM処理とPRISM処理とに切り換えて目標空燃比KCMDを算出する処理を行うため、位相遅れ時間dd=0に設定されている)。
【0057】
以上のように、制御対象モデルを、出力偏差VO2(下流側空燃比センサの出力を表す値、空燃比センサの出力を表す値)の時系列データ、および空燃比偏差DKCMD(目標空燃比を表す値)で構成した理由は以下による。すなわち、一般に、制御対象モデルでは、制御対象の入出力と所定値との偏差を、入出力を表す変数として定義した場合の方が、入出力の絶対値を変数として定義した場合よりも、モデルパラメータをより正確に同定または定義できることで、制御対象モデルの動特性を制御対象の実際の動特性に適合させることができるという事実が知られている。したがって、本実施形態の空燃比制御装置1のように、制御対象モデルを、出力偏差VO2の時系列データおよび空燃比偏差DKCMDで構成することにより、O2センサ15の出力Voutおよび目標空燃比KCMDの絶対値を変数とする場合と比べて、制御対象の実際の動特性に対する制御対象モデルの動特性の適合性を向上させることができ、それにより予測値PREVO2の算出精度を向上させることができる。
【0058】
また、予測値PREVO2は、目標空燃比KCMDの混合気が吸気系に供給されてから予測時間dtが経過した後の出力偏差VO2(k+dt)を予測した値であり、上記式(1)に基づき、予測値PREVO2の算出式を導出すると、下式(3)が得られる。
【0059】
この式(3)では、出力偏差VO2(k)の未来値に相当するVO2(k+dt−1),VO2(k+dt−2)の算出が必要となり、実際にプログラム化するのは困難である。そのため、マトリクスA、Bを、モデルパラメータa1,a2,b1を用いて図4に示す式(4)、(5)のように定義するとともに、上式(3)の漸化式を繰り返し用いることにより、上式(3)を変形すると、図4に示す式(6)が得られる。予測アルゴリズムすなわち予測値PREVO2の算出式として、この式(6)を用いた場合、予測値PREVO2が、出力偏差VO2および空燃比偏差DKCMDにより算出される。
【0060】
次に、LAF出力偏差DKACT(上流側空燃比センサの出力を表す値)を、LAFセンサ14の出力KACT(=φin)と基準値FLAFBASEとの偏差(KACT−FLAFBASE)として定義すると、DKACT(k)=DKCMD(k−d')の関係が成立するので、この関係を図4の式(6)に適用すると、図4に示す式(7)が得られる。
【0061】
以上の式(6)または式(7)により算出される予測値PREVO2を用い、後述するように目標空燃比KCMDを算出することによって、制御対象の入出力間の応答遅れやむだ時間を適切に補償しながら、目標空燃比KCMDを算出することができる。特に、予測アルゴリズムとして、上記式(7)を用いた場合、予測値PREVO2が、出力偏差VO2、LAF出力偏差DKACTおよび目標空燃比KCMDにより算出されるので、第1触媒装置8aに実際に供給される排気ガスの空燃比の状態が反映された値として、予測値PREVO2を算出でき、その算出精度すなわち予測精度を上記式(6)を用いた場合よりも向上させることができる。また、式(7)を用いた場合において、d'≦1と見なせるときには、空燃比偏差DKCMDを用いることなく、出力偏差VO2およびLAF出力偏差DKACTのみにより、予測値PREVO2を算出できる。本実施形態では、LAFセンサ14がエンジン3に設けられているので、予測アルゴリズムとして上記式(7)を採用する。
【0062】
なお、前述した式(1)の制御対象モデルは、DKACT(k)=DKCMD(k−d')の関係を適用することにより、出力偏差VO2およびLAF出力偏差DKACTを変数とするモデルとして定義することも可能である。
【0063】
次に、オンボード同定器23(同定手段)について説明する。このオンボード同定器23は、以下に述べる逐次型同定アルゴリズムにより、前述した式(1)のモデルパラメータa1,a2,b1を同定(算出)するものである。具体的には、図5に示す(8),(9)により、モデルパラメータのベクトルθ(k)を算出する。同図の式(8)において、KP(k)は、ゲイン係数のベクトルであり、ide_f(k)は同定誤差フィルタ値である。また、式(9)におけるθ(k)Tは、θ(k)の転置行列を表し、a1'(k)、a2'(k)およびb1'(k)は、後述するリミット処理を施す前のモデルパラメータを表している。なお、以下の説明では、「ベクトル」という表記を適宜、省略する。
【0064】
上記式(8)の同定誤差フィルタ値ide_f(k)は、図5に示す式(11)〜(13)により算出される同定誤差ide(k)に、図5の式(10)に示す移動平均フィルタリング処理を施した値である。図5の式(10)のnは、移動平均フィルタリング処理のフィルタ次数(1以上の整数)を表しており、式(12)のVO2HAT(k)は、出力偏差VO2の同定値を表している。
【0065】
この同定誤差フィルタ値ide_f(k)を用いる理由は以下による。すなわち、本実施形態の制御対象は、目標空燃比KCMDを制御入力とし、O2センサ15の出力Voutを制御対象の出力とするものであり、その周波数特性としてはローパス特性を有している。このようなローパス特性を有する制御対象では、オンボード同定器23の同定アルゴリズム、具体的には後述する重み付き最小2乗法アルゴリズムの周波数重み特性に起因して、制御対象の高周波特性が強調された状態で、モデルパラメータが同定されるため、制御対象モデルのゲイン特性が制御対象の実際のゲイン特性よりも低くなる傾向を示す。その結果、空燃比制御装置1によりADSM処理またはPRISM処理が実行された際、オーバーゲイン状態になることで、制御系が発散状態になり、不安定になる可能性がある。
【0066】
したがって、本実施形態では、重み付き最小2乗法アルゴリズムの周波数重み特性を適切に補正し、制御対象モデルのゲイン特性を、制御対象の実際のゲイン特性に一致させるために、上記同定誤差ide(k)に移動平均フィルタリング処理を施した同定誤差フィルタ値ide_f(k)を用いるとともに、後述するように、移動平均フィルタリング処理のフィルタ次数nを、排気ガスボリュームAB_SVに応じて設定している。
【0067】
さらに、前述した図5の式(8)のゲイン係数のベクトルKP(k)は、図5の式(14)により算出される。この式(14)のP(k)は、図5の式(15)で定義される3次の正方行列である。
【0068】
以上のような同定アルゴリズムでは、式(15)の重みパラメータλ1、λ2の設定により、以下の4つの同定アルゴリズムのうちの1つが選択される。
すなわち、
λ1=1,λ2=0 ;固定ゲインアルゴリズム
λ1=1,λ2=1 ;最小2乗法アルゴリズム
λ1=1,λ2=λ ;漸減ゲインアルゴリズム
λ1=λ,λ2=1 ;重み付き最小2乗法アルゴリズム
ただし、λは、0<λ<1に設定される所定値。
【0069】
本実施形態では、これらの4つの同定アルゴリズムのうちの重み付き最小2乗法アルゴリズムを採用する。これは、重みパラメータλ1の値をエンジン3の運転状態、具体的には排気ガスボリュームAB_SVに応じて設定することにより、同定精度と、モデルパラメータの最適値への収束速度とを適切に設定できることによる。例えば、低負荷運転状態のときには、それに応じて重みパラメータλ1の値を値1に近い値に設定することで、すなわち最小2乗法アルゴリズムに近いアルゴリズムに設定することで、良好な同定精度を確保できるとともに、高負荷運転状態のときには、それに応じて重みパラメータλ1の値を低負荷運転状態のときよりも小さい値に設定することにより、モデルパラメータを迅速に最適値に収束させることができる。以上のように、重みパラメータλ1の値を排気ガスボリュームAB_SVに応じて設定することにより、同定精度と、モデルパラメータの最適値への収束速度とを適切に設定することができ、それにより、触媒装置8a,8bで浄化される排気ガスの特性、すなわち触媒後排気ガス特性を向上させることができる。
【0070】
以上の式(8)〜(15)の同定アルゴリズムにおいて、前述したDKACT(k)=DKCMD(k−d')の関係を適用すると、図6に示す式(16)〜(23)の同定アルゴリズムが得られる。本実施形態では、LAFセンサ14がエンジン3に設けられているので、これらの式(16)〜(23)を用いる。これらの式(16)〜(23)を用いた場合、前述した理由により、モデルパラメータを、第1触媒装置8aに実際に供給される排気ガスの空燃比の状態がより反映された値として同定することができ、それにより、上記式(8)〜(15)の同定アルゴリズムを用いた場合よりも、モデルパラメータの同定精度を向上させることができる。
【0071】
また、このオンボード同定器23では、以上の同定アルゴリズムにより算出されたモデルパラメータa1'(k)、a2'(k)およびb1'(k)に、後述するリミット処理を施すことにより、モデルパラメータa1(k)、a2(k)およびb1(k)が算出される。さらに、前述した状態予測器22では、このようにリミット処理を施した後のモデルパラメータa1(k)、a2(k)およびb1(k)に基づき、予測値PREVO2が算出される。
【0072】
次に、DSMコントローラ24(目標空燃比算出手段)について説明する。このDSMコントローラ24は、ΔΣ変調アルゴリズムを応用した制御アルゴリズムにより、状態予測器22で算出された予測値PREVO2に基づき、制御入力φop(k)(=目標空燃比KCMD)を生成(算出)するとともに、これを制御対象に入力することにより、制御対象の出力としてのO2センサ15の出力Voutを目標値Vopに収束させるように制御するものである。
【0073】
まず、一般的なΔΣ変調アルゴリズムについて説明する。図7は、ΔΣ変調アルゴリズムを適用したコントローラ26により、制御対象27を制御する制御系の構成を示している。同図に示すように、このコントローラ26では、差分器26aにより、参照信号r(k)と遅延素子26bで遅延されたDSM信号u(k−1)との偏差として偏差信号δ(k)が生成される。次に、積分器26cにより、偏差積分値σd(k)が、偏差信号δ(k)と遅延素子26dで遅延された偏差積分値σd(k−1)との和の信号として生成される。次いで、量子化器26e(符号関数)により、DSM信号u(k)が、この偏差積分値σd(k)を符号化した信号として生成される。そして、以上のように生成されたDSM信号u(k)が制御対象27に入力されることにより、出力信号 y(k)が制御対象27から出力される。
【0074】
以上のΔΣ変調アルゴリズムは、以下の数式(24)〜(26)で表される。
δ(k)=r(k)−u(k−1) ……(24)
σd(k)=σd(k−1)+δ(k) ……(25)
u(k)=sgn(σd(k)) ……(26)
ただし、符号関数sgn(σd(k))の値は、σd(k)≧0のときにはsgn(σd(k))=1となり、σd(k)<0のときにはsgn(σd(k))=−1となる(なお、σd(k)=0のときに、sgn(σd(k))=0と設定してもよい)。
【0075】
すなわち、このΔΣ変調アルゴリズムでは、DSM信号u(k)は、偏差積分値σd(k)が値0以上のときに値1に、偏差積分値σd(k)が値0よりも小さいときに値−1にそれぞれ設定される。
【0076】
図8は、以上の制御系の制御シミュレーション結果を示している。同図に示すように、正弦波状の参照信号r(k)を制御系に入力した場合、DSM信号u(k)が矩形波状の信号として生成され、これを制御対象27に入力することにより、参照信号r(k)と異なる振幅で同じ周波数の、ノイズを有するものの全体として同様の波形の出力信号y(k)が、制御対象27から出力される。このように、ΔΣ変調アルゴリズムの特性は、参照信号r(k)から生成されたDSM信号u(k)を制御対象27に入力したときに、制御対象27の出力y(k)が、参照信号r(k)に対して、異なる振幅で同じ周波数の、全体として同様の波形の信号となるような値として、DSM信号u(k)を生成できるという点にある。言い換えれば、DSM信号u(k)を、参照信号r(k)が制御対象27の実際の出力y(k)に再現されるような値として、生成(算出)できるという点にある。
【0077】
DSMコントローラ24は、このようなΔΣ変調アルゴリズムの特性を利用し、O2センサ15の出力Voutを目標値Vopに収束させるための制御入力φop(k)を算出するものである。その原理について説明すると、例えば図9に1点鎖線で示すように、出力偏差VO2が値0に対して揺らいでいる場合(すなわち、O2センサ15の出力Voutが目標値Vopに対して揺らいでいる場合)、出力偏差VO2を値0に収束させる(すなわち出力Voutを目標値Vopに収束させる)には、図9に破線で示す、出力偏差VO2を打ち消すような逆位相波形の出力偏差VO2*が生じるように、制御入力φop(k)を生成すればよい。
【0078】
しかし、前述したように、本実施形態の制御対象では、制御入力φop(k)としての目標空燃比KCMDが制御対象に入力されてからO2センサ15の出力Voutに反映されるまでに、予測時間dt分の時間遅れが発生するため、現在の出力偏差VO2に基づいて、制御入力φop(k)を算出した場合の出力偏差VO2#は、図9に実線で示すように、出力偏差VO2*に対して遅れを生じ、それにより、制御タイミングのずれが生じてしまう。したがって、これを補償するために、本実施形態のADSMコントローラ20におけるDSMコントローラ24では、出力偏差VO2の予測値PREVO2を用いることにより、制御入力φop(k)が、制御タイミングのずれを生じることなく、現在の出力偏差VO2を打ち消すような出力偏差(逆位相波形の出力偏差VO2*と同様の出力偏差)を生じさせる信号として生成される。
【0079】
具体的には、このDSMコントローラ24では、図10に示すように、反転増幅器24aにより、参照信号r(k)が、値−1、参照信号用のゲインGdおよび予測値PREVO2(k)を互いに乗算した信号として生成される。次に、差分器24bにより、この参照信号r(k)と遅延素子24cで遅延されたDSM信号u''(k−1)との偏差として偏差信号δ(k)が生成される。
【0080】
次いで、積分器24dにより、偏差積分値σd(k)が、偏差信号δ(k)と遅延素子24eで遅延された偏差積分値σd(k−1)との和の信号として生成され、次に、量子化器24f(符号関数)により、DSM信号u''(k)が、この偏差積分値σd(k)を符号化した値として生成される。そして、増幅器24gにより、増幅DSM信号u(k)がDSM信号u''(k)を所定のゲインFdで増幅した値として生成され、次に、加算器24hにより、この増幅DSM信号u(k)を所定の基準値FLAFBASEに加算した値として、制御入力φop(k)が生成される。
【0081】
以上のDSMコントローラ24の制御アルゴリズムは、以下の式(27)〜(32)で表される。
r(k)=−1・Gd・PREVO2(k) ……(27)
δ(k)=r(k)−u''(k−1) ……(28)
σd(k)=σd(k−1)+δ(k) ……(29)
u''(k)=sgn(σd(k)) ……(30)
u(k)=Fd・u''(k) ……(31)
φop(k)=FLAFBASE+u(k) ……(32)
ここで、Gd,Fdはゲインを表す。また、符号関数sgn(σd(k))の値は、σd(k)≧0のときにはsgn(σd(k))=1となり、σd(k)<0のときにはsgn(σd(k))=−1となる(なお、σd(k)=0のときに、sgn(σd(k))=0と設定してもよい)。
【0082】
このDSMコントローラ24では、以上の式(27)〜(32)に示す制御アルゴリズムにより、前述したように、制御入力φop(k)が、制御タイミングのずれを生じることなく、出力偏差VO2を打ち消すような出力偏差VO2*を生じさせる値として算出される。すなわち、制御入力φop(k)が、O2センサ15の出力Voutを目標値Vopに収束させることができる値として算出される。また、制御入力φop(k)が、増幅DSM信号u(k)を所定の基準値FLAFBASEに加算した値として算出されるので、制御入力φop(k)を値0を中心して正負反転する値だけでなく、基準値FLAFBASEを中心として増減を繰り返す値として算出できる。これにより、通常のΔΣ変調アルゴリズムと比べて、制御の自由度を高めることができる。
【0083】
次に、前記PRISMコントローラ21について説明する。このPRISMコントローラ21は、以下に述べるオンボード同定型スライディングモード制御処理(以下「PRISM処理」という)の制御アルゴリズムにより、O2センサ15の出力Voutを目標値Vopに収束させるための目標空燃比KCMDを算出するものであり、状態予測器22、オンボード同定器23およびスライディングモードコントローラ(以下「SLDコントローラ」という)25により構成されている。なお、このPRISM処理の具体的なプログラムについては後述する。
【0084】
このPRISMコントローラ21のうちの状態予測器22およびオンボード同定器23については、既に説明したので、ここではSLDコントローラ25についてのみ説明する。このSLDコントローラ25は、スライディングモード制御アルゴリズムに基づいてスライディングモード制御を行うものであり、以下、一般的なスライディングモード制御アルゴリズムについて説明する。このスライディングモード制御アルゴリズムでは、前述した式(1)の離散時間系モデルを制御対象モデルとして用いるため、切換関数σは、下式(33)に示すように、出力偏差VO2の時系列データの線形関数として設定される。
σ(k)=S1・VO2(k)+S2・VO2(k−1) ……(33)
ここで、S1,S2は、−1<(S2/S1)<1の関係が成立するように設定される所定の係数である。
【0085】
一般にスライディングモード制御アルゴリズムでは、切換関数σが2つの状態変数(本実施形態では出力偏差VO2の時系列データ)で構成されている場合、2つの状態変数で構成される位相空間は、これらをそれぞれ縦軸および横軸とする2次元の位相平面となるため、この位相平面上において、σ=0を満たす2つの状態変数の値の組み合わせは、切換直線と呼ばれる直線上に載ることになる。したがって、制御対象への制御入力を、2つの状態変数の組み合わせが切換直線上に収束する(載る)ように適切に決定することにより、2つの状態変数をいずれも、値0になる平衡点に収束(スライディング)させることができる。さらに、スライディングモード制御アルゴリズムでは、切換関数σの設定により、状態変数の動特性、より具体的には収束挙動や収束速度を指定することができる。例えば、本実施形態のように、切換関数σが2つの状態変数で構成されている場合には、切換直線の傾きを値1に近づけると、状態変数の収束速度が遅くなる一方、値0に近づけると、収束速度が速くなる。
【0086】
本実施形態では、前記式(33)に示すように、切換関数σが出力偏差VO2の2つの時系列データ、すなわち出力偏差VO2の今回値VO2(k)および前回値VO2(k−1)により構成されているので、これらの今回値VO2(k)および前回値VO2(k−1)の組み合わせを切換直線上に収束させるように、制御対象への制御入力すなわち目標空燃比KCMDを設定すればよい。具体的には、制御量Usl(k)を、基準値FLAFBASEとの和が目標空燃比KCMDとなる値として定義すると、今回値VO2(k)および前回値VO2(k−1)の組み合わせを切換直線上に収束させるための制御量Usl(k)は、適応スライディングモード制御アルゴリズムにより、図11に示す式(34)のように、等価制御入力Ueq(k)、到達則入力Urch(k)および適応則入力Uadp(k)の総和として設定される。
【0087】
この等価制御入力Ueq(k)は、出力偏差VO2の今回値VO2(k)および前回値VO2(k−1)の組み合わせを切換直線上に拘束しておくためのものであり、具体的には、図11に示す式(35)のように定義される。また、到達則入力Urch(k)は、外乱やモデル化誤差などにより、出力偏差VO2の今回値VO2(k)および前回値VO2(k−1)の組み合わせが切換直線上から外れた際に、これらを切換直線上に収束させるためのものであり、具体的には、図11に示す式(36)のように定義される。この式(36)において、Fはゲインを表す。
【0088】
さらに、適応則入力Uadp(k)は、制御対象の定常偏差、モデル化誤差および外乱の影響を抑制しながら、出力偏差VO2の今回値VO2(k)および前回値VO2(k−1)の組み合わせを、切換超平面上に確実に収束させるためのものであり、具体的には、図11に示す式(37)のように定義される。この式(37)において、Gはゲインを、ΔTは制御周期をそれぞれ表す。
【0089】
本実施形態のPRISMコントローラ21のSLDコントローラ25では、前述したように、出力偏差VO2に代えて予測値PREVO2を用いるので、PREVO2(k)≒VO2(k+dt)の関係を適用することにより、以上の式(33)〜(37)のアルゴリズムを、図12に示す式(38)〜(42)に書き換えて用いる。この式(38)におけるσPREは、予測値PREVO2を用いたときの切換関数(以下「予測切換関数」という)の値である。すなわち、このSLDコントローラ25では、以上のアルゴリズムで算出される制御量Usl(k)を基準値FLAFBASEに加算することによって、目標空燃比KCMDが算出される。
【0090】
以下、ECU2により実行される燃料噴射量の算出処理について、図13を参照しながら説明する。なお、以下の説明では、今回値であることを示す記号(k)を適宜、省略する。図13は、この制御処理のメインルーチンを示しており、本処理は、TDC信号の入力に同期して割り込み実行される。この処理では、後述する適応空燃比制御処理、またはマップ検索処理により算出された目標空燃比KCMDを用いることによって、燃料噴射量TOUTが気筒ごとに算出される。
【0091】
まず、ステップ1(図では「S1」と略す。以下同じ)において、前述した各種のセンサ10〜19の出力を読み込むとともに、読み込んだデータをRAM内に記憶する。
【0092】
次に、ステップ2に進み、基本燃料噴射量Timを算出する。この処理では、エンジン回転数NEおよび吸気管内絶対圧PBAに応じて、図示しないマップを検索することにより、基本燃料噴射量Timを算出する。
【0093】
次いで、ステップ3に進み、総補正係数KTOTALを算出する。この総補正係数KTOTALは、各種の運転パラメータ(例えば吸気温TAや、大気圧PA、エンジン水温TW、アクセル開度APなど)に応じて、各種のテーブルやマップを検索することで各種の補正係数を算出するとともに、これらの各種の補正係数を互いに乗算することにより、算出される。
【0094】
次に、ステップ4に進み、適応制御フラグF_PRISMONの設定処理を実行する。この処理の内容は図示しないが、具体的には、以下の(a)〜(f)の条件がいずれも成立しているときには、適応空燃比制御処理で算出された目標空燃比KCMDを使用する条件が成立しているとして、それを表すために、適応制御フラグF_PRISMONが「1」にセットされる。一方、(a)〜(f)の条件のうちの少なくとも1つが成立していないときには、適応制御フラグF_PRISMONが「0」にセットされる。
(a)LAFセンサ14およびO2センサ15がいずれも活性化していること。
(b)エンジン3がリーンバーン運転中でないこと。
(c)スロットル弁5が全開状態でないこと。
(d)点火時期の遅角制御中でないこと。
(e)フューエルカット運転中でないこと。
(f)エンジン回転数NEおよび吸気管内絶対圧PBAがいずれも、所定の範囲囲内の値であること。
【0095】
次に、ステップ5に進み、ステップ4で設定された適応制御フラグF_PRISMONが「1」であるか否かを判別する。この判別結果がYESのときには、ステップ6に進み、目標空燃比KCMDを、後述する適応空燃比制御処理で算出された適応目標空燃比KCMDSLDに設定する。
【0096】
一方、ステップ5の判別結果がNOのときには、ステップ7に進み、目標空燃比KCMDをマップ値KCMDMAPに設定する。このマップ値KCMDMAPは、エンジン回転数NEおよび吸気管内絶対圧PBAに応じて、図示しないマップを検索することにより、算出される。
【0097】
以上のステップ6または7に続くステップ8では、オブザーバフィードバック補正係数#nKLAFを気筒ごとに算出する。このオブザーバフィードバック補正係数#nKLAFは、気筒ごとの実際の空燃比のばらつきを補正するためのものであり、具体的には、オブザーバによりLAFセンサ14の出力KACTから気筒ごとの実際の空燃比を推定し、これらの推定した空燃比に応じて、PID制御により算出される。なお、このオブザーバフィードバック補正係数#nKLAFの記号#nは、気筒の番号#1〜#4を表すものであり、これは、後述する要求燃料噴射量#nTCYLおよび最終燃料噴射量#nTOUTにおいても同様である。
【0098】
次いで、ステップ9に進み、フィードバック補正係数KFBを算出する。このフィードバック補正係数KFBは、LAFセンサ14の出力KACTを目標空燃比KCMDに一致させるためのものであり、具体的には、以下のように算出される。すなわち、LAFセンサ14の出力KACTと目標空燃比KCMDとの偏差に応じて、PID制御によりフィードバック係数KLAFを算出する。また、図示しないSelf Tuning Regulator 型の適応制御器によりフィードバック補正係数KSTRを算出し、これを目標空燃比KCMDで除算することにより、フィードバック補正係数kstrを算出する。そして、エンジン3の運転状態に応じて、これらの2つのフィードバック係数KLAFおよびフィードバック補正係数kstrの一方を、フィードバック補正係数KFBとして設定する。
【0099】
次いで、ステップ10に進み、補正目標空燃比KCMDMを算出する。この補正目標空燃比KCMDMは、空燃比A/Fの変化による充填効率の変化を補償するためのものであり、前述したステップ6または7で算出された目標空燃比KCMDに応じて、図示しないテーブルを検索することにより算出される。
【0100】
次に、ステップ11に進み、以上のように算出した基本燃料噴射量Tim、総補正係数KTOTAL、オブザーバフィードバック補正係数#nKLAF、フィードバック補正係数KFB、および補正目標空燃比KCMDMを用い、下式(43)により、気筒ごとの要求燃料噴射量#nTCYLを算出する。
【0101】
次に、ステップ12に進み、要求燃料噴射量#nTCYLを付着補正することにより、最終燃料噴射量#nTOUTを算出する。この最終燃料噴射量#nTOUTは、具体的には、今回の燃焼サイクルでインジェクタ6から噴射された燃料が燃焼室の内壁面に付着する割合などを、運転状態に応じて算出し、そのように算出した割合に基づいて、要求燃料噴射量#nTCYLを補正することにより、算出される。
【0102】
次いで、ステップ13に進み、以上のように算出した最終燃料噴射量#nTOUTに基づく駆動信号を、対応する気筒のインジェクタ6に出力した後、本処理を終了する。以上により、LAFセンサ14の出力KACTを目標空燃比KCMDに一致させるように、混合気の空燃比がフィードバック制御される。
【0103】
次に、図14および図15を参照しながら、ADSM処理およびPRISM処理を含む適応空燃比制御処理について説明する。この処理は、後述する理由により、100msecの周期で実行される。また、この処理では、エンジン3の運転状態に応じて、ADSM処理、PRISM処理、またはスライディングモード制御量DKCMDSLDを所定値SLDHOLDに設定する処理により、目標空燃比KCMDが算出される。
【0104】
この処理では、まず、ステップ20において、F/C後判定処理を実行する。この処理の内容は図示しないが、この処理では、フューエルカット運転中は、それを表すためにF/C後判定フラグF_AFCが「1」にセットされ、フューエルカット運転の終了後、所定時間X_TM_AFCが経過したときには、それを表すためにF/C後判定フラグF_AFCが「0」にセットされる。
【0105】
次に、ステップ21に進み、車速VPに基づいて、エンジン3を搭載した車両が発進したか否かを判定する発進判定処理を実行する。図16に示すように、この処理では、まず、ステップ49において、アイドル運転フラグF_IDLEが「1」であるか否かを判別する。このアイドル運転フラグF_IDLEは、アイドル運転中であるときに「1」に、それ以外のときに「0」にセットされる。
【0106】
この判別結果がYESで、アイドル運転中であるときには、ステップ50に進み、車速VPが所定車速VSTART(例えば1km/h)より小さいか否かを判別する。この判別結果がYESで、停車中であるときには、ステップ51に進み、ダウンカウント式の第1発進判定タイマのタイマ値TMVOTVSTを第1所定時間TVOTVST(例えば3msec)に設定する。
【0107】
次いで、ステップ52に進み、ダウンカウント式の第2発進判定タイマのタイマ値TMVSTを、上記第1所定時間TVOTVSTよりも長い第2所定時間TVST(例えば500msec)に設定する。次いで、ステップ53,54において、第1および第2発進フラグF_VOTVST,F_VSTをいずれも「0」にセットした後、本処理を終了する。
【0108】
一方、ステップ49または50の判別結果がNOのとき、すなわちアイドル運転中でないか、または車両が発進したときには、ステップ55に進み、第1発進判定タイマのタイマ値TMVOTVSTが値0より大きいか否かを判別する。この判別結果がYESで、アイドル運転の終了後または車両の発進後、第1所定時間TVOTVSTが経過していないときには、第1発進モード中であるとして、ステップ56に進み、それを表すために第1発進フラグF_VOTVSTを「1」にセットする。
【0109】
一方、ステップ55の判別結果がNOで、アイドル運転の終了後または車両の発進後、第1所定時間TVOTVSTが経過したときには、第1発進モードが終了したとして、ステップ57に進み、第1発進フラグF_VOTVSTを「0」にセットする。
【0110】
ステップ56または57に続くステップ58では、第2発進判定タイマのタイマ値TMVSTが値0より大きいか否かを判別する。この判別結果がYESで、アイドル運転の終了後または車両の発進後、第2所定時間TVSTが経過していないときには、第2発進モード中であるとして、ステップ59に進み、それを表すために第2発進フラグF_VSTを「1」にセットした後、本処理を終了する。
【0111】
一方、ステップ59の判別結果がNOで、アイドル運転の終了後または車両の発進後、第2所定時間TVSTが経過したときには、第2発進モードが終了したとして、前記ステップ54を実行した後、本処理を終了する。
【0112】
図14に戻り、ステップ21に続くステップ22では、状態変数の設定処理を実行する。図示しないが、この処理では、RAM内に記憶されている、目標空燃比KCMD、LAFセンサ14の出力KACTおよび出力偏差VO2の時系列データをいずれも、1サンプリングサイクル分ずつ過去側にシフトさせる。その後、KCMD、KACTおよびVO2の時系列データの最新の値と、基準値FLAFBASEと、後述する適応補正項FLAFADPとに基づき、KCMD、KACTおよびVO2の今回値を算出する。
【0113】
次に、ステップ23に進み、PRISM/ADSM処理の実行判定処理を行う。この処理は、PRISM処理またはADSM処理の実行条件が成立しているか否かを判定するものであり、具体的には、図17に示すフローチャートのように実行される。
【0114】
すなわち、図17のステップ60〜63において、以下の(g)〜(j)の条件がいずれも成立しているときには、PRISM処理またはADSM処理を実行すべき運転状態にあるとして、それを表すために、ステップ64で、PRISM/ADSM実行フラグF_PRISMCALを「1」にセットした後、本処理を終了する。一方、(g)〜(j)の条件の少なくとも1つが成立していないときには、PRISM処理またはADSM処理を実行すべき運転状態にないとして、それを表すために、ステップ65で、PRISM/ADSM実行フラグF_PRISMCALを「0」にセットした後、本処理を終了する。
(g)O2センサ15が活性化していること。
(h)LAFセンサ14が活性化していること。
(i)エンジン3がリーンバーン運転中でないこと。
(j)点火時期の遅角制御中でないこと。
【0115】
図14に戻り、ステップ23に続くステップ24では、同定器演算の実行判定処理を行う。この処理は、オンボード同定器23によるパラメータ同定の実行条件が成立しているか否かを判定するものであり、具体的には、図18に示すフローチャートのように実行される。
【0116】
すなわち、図18のステップ70および71の判別結果がいずれもNOのとき、言い換えれば、スロットル弁開度θTHが全開状態でなく、かつフューエルカット運転中でないときには、パラメータ同定を実行すべき運転状態であるとして、ステップ72に進み、同定実行フラグF_IDCALを「1」にセットした後、本処理を終了する。一方、ステップ70または71の判別結果がYESのときには、パラメータ同定を実行すべき運転状態にないとして、ステップ73に進み、同定実行フラグF_IDCALを「0」にセットした後、本処理を終了する。
【0117】
図14に戻り、ステップ24に続くステップ25では、各種パラメータ(排気ガスボリュームAB_SVなど)を算出する。この処理の具体的な内容は、後述する。
【0118】
次に、ステップ26に進み、前記ステップ23で設定されたPRISM/ADSM実行フラグF_PRISMCALが「1」であるか否かを判別する。この判別結果がYESで、PRISM処理またはADSM処理の実行条件が成立しているときには、ステップ27に進み、前記ステップ24で設定された同定実行フラグF_IDCALが「1」であるか否かを判別する。
【0119】
この判別結果がYESで、オンボード同定器23によるパラメータ同定を実行すべき運転状態のときには、ステップ28に進み、パラメータ初期化フラグF_IDRSETが「1」であるか否かを判別する。この判別結果がNOで、RAMに記憶されているモデルパラメータa1,a2,b1の初期化が不要であるときには、後述するステップ31に進む。
【0120】
一方、この判別結果がYESで、モデルパラメータa1,a2,b1の初期化が必要であるときには、ステップ29に進み、モデルパラメータa1,a2,b1を、それぞれの初期値に設定した後、それを表すためにステップ30に進み、パラメータ初期化フラグF_IDRSETを「0」にセットする。
【0121】
このステップ30または28に続くステップ31では、オンボード同定器23の演算を実行し、モデルパラメータa1,a2,b1を同定した後、後述する図15のステップ32に進む。このオンボード同定器23の演算の具体的な内容については、後述する。
【0122】
一方、ステップ27の判別結果がNOで、パラメータ同定を実行すべき運転状態でないときには、以上のステップ28〜31をスキップして、図15のステップ32に進む。ステップ27または31に続くステップ32では、モデルパラメータa1,a2,b1として、同定値または所定値を選択する。この処理の内容は図示しないが、具体的には、前記ステップ24で設定された同定実行フラグF_IDCALが「1」のときには、モデルパラメータa1,a2,b1をステップ31で同定された同定値に設定する。一方、同定実行フラグF_IDCALが「0」のときには、モデルパラメータa1,a2,b1を所定値に設定する。
【0123】
次に、ステップ33に進み、後述するように、状態予測器22の演算を実行し、予測値PREVO2を算出する。その後、ステップ34に進み、後述するように、制御量Uslを算出する。
【0124】
次いで、ステップ35に進み、SLDコントローラ25の安定判別を実行する。この処理の内容は図示しないが、具体的には、予測切換関数σPREの値に基づき、SLDコントローラ25によるスライディングモード制御が安定状態にあるか否かを判別する。
【0125】
次に、ステップ36および37において、後述するように、SLDコントローラ25およびDSMコントローラ24により、スライディングモード制御量DKCMDSLDおよびΔΣ変調制御量DKCMDDSMをそれぞれ算出する。
【0126】
次いで、ステップ38に進み、後述するように、SLDコントローラ25により算出されたスライディングモード制御量DKCMDSLD、またはDSMコントローラ24により算出されたΔΣ変調制御量DKCMDDSMを用いて、適応目標空燃比KCMDSLDを算出する。この後、ステップ39に進み、後述するように、適応補正項FLAFADPを算出した後、本処理を終了する。
【0127】
一方、図14に戻り、前記ステップ26の判別結果がNOで、PRISM処理およびADSM処理の実行条件がいずれも成立していないときには、ステップ40に進み、パラメータ初期化フラグF_IDRSETを「1」にセットする。次に、図15のステップ41に進み、スライディングモード制御量DKCMDSLDを所定値SLDHOLDにセットする。次いで、前述したステップ38,39を実行した後、本処理を終了する。
【0128】
次に、図19を参照しながら、前述したステップ25の各種パラメータを算出する処理について説明する。この処理では、まず、ステップ80において、下式(44)により、排気ガスボリュームAB_SV(空間速度の推定値)を算出する。
AB_SV=(NE/1500)・PBA・X_SVPRA ……(44)
ここで、X_SVPRAは、エンジン排気量に基づいて決定される所定の係数である。
【0129】
次に、ステップ81に進み、前述した空燃比操作系のむだ時間KACT_D(=d')、排気系のむだ時間CAT_DELAY(=d)および予測時間dtを算出する。具体的には、ステップ80で算出された排気ガスボリュームAB_SVに応じて、図20に示すテーブルを検索することにより、むだ時間KACT_D,CAT_DELAYをそれぞれ算出するとともに、これらの和(KACT_D+CAT_DELAY)を予測時間dtとして設定する。すなわち、この制御プログラムでは、位相遅れ時間ddが値0に設定される。
【0130】
このテーブルでは、排気ガスボリュームAB_SVが大きいほど、むだ時間KACT_D,CAT_DELAYがより小さい値に設定されている。これは、排気ガスボリュームAB_SVが大きいほど、排気ガスの流速が大きくなることで、むだ時間KACT_D,CAT_DELAYが短くなることによる。以上のように、むだ時間KACT_D,CAT_DELAYおよび予測時間dtが、排気ガスボリュームに応じて算出されるので、これらを用いて算出した出力偏差VO2の予測値PREVO2に基づき、後述する適応目標空燃比KCMDSLDを算出することにより、制御対象の入出力間の制御タイミングのずれを解消することができる。また、モデルパラメータa1,a2,b1が、上記むだ時間CAT_DELAYを用いて同定されるので、制御対象モデルの動特性を、制御対象の実際の動特性に適合させることができ、それにより、制御対象の入出力間の制御タイミングのずれをさらに解消することができる。
【0131】
次に、ステップ82に進み、同定アルゴリズムの重みパラメータλ1,λ2の値を算出する。具体的には、重みパラメータλ2を値1に設定すると同時に、重みパラメータλ1を、排気ガスボリュームAB_SVに応じて、図21に示すテーブルを検索することにより算出する。
【0132】
このテーブルでは、排気ガスボリュームAB_SVが大きいほど、重みパラメータλ1がより小さい値に設定されており、言い換えれば、排気ガスボリュームAB_SVが小さいほど、重みパラメータλ1がより大きくかつ値1により近い値に設定されている。これは、排気ガスボリュームAB_SVが大きいほど、言い換えれば高負荷運転状態であるほど、モデルパラメータの同定をより迅速に行う必要があるので、重みパラメータλ1をより小さく設定することによって、モデルパラメータの最適値への収束速度を高めるためである。これに加えて、排気ガスボリュームAB_SVが小さいほど、すなわち低負荷運転状態であるほど、空燃比が変動しやすくなり、触媒後排気ガス特性が不安定になりやすいことで、モデルパラメータの良好な同定精度を確保する必要があるので、重みパラメータλ1を値1に近づける(最小2乗法アルゴリズムに近づける)ことによって、モデルパラメータの同定精度をより高めるためである。
【0133】
次に、ステップ83に進み、モデルパラメータa1,a2の値を制限するための下限値X_IDA2Lと、モデルパラメータb1の値を制限するための下限値X_IDB1Lおよび上限値X_IDB1Hとを、排気ガスボリュームAB_SVに応じて、図22に示すテーブルを検索することにより算出する。
【0134】
このテーブルでは、下限値X_IDA2Lは、排気ガスボリュームAB_SVが大きいほど、より大きい値に設定されている。これは、排気ガスボリュームAB_SVの変化に応じたむだ時間の増減に伴い、制御系が安定状態となるモデルパラメータa1,a2の組み合わせが変化することによる。また、下限値X_IDB1Lおよび上限値X_IDB1Hも、排気ガスボリュームAB_SVが大きいほど、より大きい値に設定されている。これは、排気ガスボリュームAB_SVが大きいほど、触媒前空燃比(第1触媒装置8aよりも上流側の排気ガスの空燃比)がO2センサ15の出力Voutに及ぼす影響の度合、すなわち制御対象のゲインがより大きくなることによる。
【0135】
次いで、ステップ84に進み、移動平均フィルタリング処理のフィルタ次数nを算出した後、本処理を終了する。この処理では、フィルタ次数nを、排気ガスボリュームAB_SVに応じて、図23に示すテーブルを検索することにより、算出する。
【0136】
このテーブルでは、排気ガスボリュームAB_SVが大きいほど、フィルタ次数nがより小さい値に設定されている。これは、以下の理由による。すなわち、前述したように、排気ガスボリュームAB_SVが変化すると、制御対象の周波数特性、特にゲイン特性が変化するので、制御対象モデルのゲイン特性を、制御対象の実際のゲイン特性に一致させるためには、重み付き最小2乗法アルゴリズムの周波数重み特性を、排気ガスボリュームAB_SVに応じて適切に補正する必要がある。したがって、移動平均フィルタリング処理のフィルタ次数nを、上記テーブルのように排気ガスボリュームAB_SVに応じて設定することにより、排気ガスボリュームAB_SVの変化にかかわらず、一定の同定重みを同定アルゴリズムにおいて確保できるとともに、制御対象モデルと制御対象との間で互いのゲイン特性を一致させることができ、これにより、同定精度を向上させることができる。
【0137】
次に、図24を参照しながら、前記ステップ31のオンボード同定器23の演算処理について説明する。同図に示すように、この処理では、まず、ステップ90において、前述した式(22)より、ゲイン係数KP(k)を算出する。次に、ステップ91に進み、前述した式(20)より、出力偏差VO2の同定値VO2HAT(k)を算出する。
【0138】
次いで、ステップ92に進み、前述した式(18)(19)より、同定誤差フィルタ値ide_f(k)を算出する。次に、ステップ93に進み、前述した式(16)より、モデルパラメータのベクトルθ(k)を算出した後、ステップ94に進み、モデルパラメータのベクトルθ(k)の安定化処理を実行する。この処理については後述する。
【0139】
次いで、ステップ95に進み、前述した式(23)より、正方行列P(k)の次回値P(k+1)を算出する。この次回値P(k+1)は、次回のループでの算出において、正方行列P(k)の値として用いられる。
【0140】
以下、図25を参照しながら、上記ステップ94におけるモデルパラメータのベクトルθ(k)の安定化処理について説明する。同図に示すように、まず、ステップ100で、3つのフラグF_A1STAB,F_A2STAB,F_B1STABをいずれも「0」にセットする。
【0141】
次に、ステップ101に進み、後述するように、a1'&a2'のリミット処理を実行する。次いで、ステップ102で、後述するように、b1'のリミット処理を実行した後、本処理を終了する。
【0142】
以下、図26を参照しながら、上記ステップ101のa1'&a2'のリミット処理について説明する。同図に示すように、まず、ステップ110において、前記ステップ93で算出したモデルパラメータの同定値a2'が、前記図19のステップ83で算出された下限値X_IDA2L以上であるか否かを判別する。この判別結果がNOのときには、ステップ111に進み、制御系を安定化させるために、モデルパラメータa2を下限値X_IDA2Lに設定すると同時に、モデルパラメータa2の安定化を実行したことを表すために、フラグF_A2STABを「1」にセットする。一方、この判別結果がYESで、a2'≧X_IDA2Lのときには、ステップ112に進み、モデルパラメータa2を同定値a2'に設定する。
【0143】
これらのステップ111または112に続くステップ113では、前記ステップ93で算出したモデルパラメータの同定値a1'が、所定の下限値X_IDA1L(例えば値−2以上で値0より小さい一定値)以上であるか否かを判別する。この判別結果がNOのときには、ステップ114に進み、制御系を安定化させるために、モデルパラメータa1を下限値X_IDA1Lに設定すると同時に、モデルパラメータa1の安定化を実行したことを表すために、フラグF_A1STABを「1」にセットする。
【0144】
一方、ステップ113の判別結果がYESのときには、ステップ115に進み、同定値a1'が、所定の上限値X_IDA1H(例えば値2)以下であるか否かを判別する。この判別結果がYESで、X_IDA1L≦a1'≦X_IDA1Hのときには、ステップ116に進み、モデルパラメータa1を同定値a1'に設定する。一方、この判別結果がNOで、X_IDA1H<a1'のときには、ステップ117に進み、モデルパラメータa1を上限値X_IDA1Hに設定すると同時に、モデルパラメータa1の安定化を実行したことを表すために、フラグF_A1STABを「1」にセットする。
【0145】
これらのステップ114、116または117に続くステップ118では、以上のように算出したモデルパラメータa1の絶対値と、モデルパラメータa2との和(|a1|+a2)が、所定の判定値X_A2STAB(例えば値0.9)以下であるか否かを判別する。この判別結果がYESのときには、モデルパラメータa1,a2の組み合わせが、制御系の安定性を確保できる範囲(図27にハッチングで示す規制範囲)内にあるとして、そのまま本処理を終了する。
【0146】
一方、ステップ118の判別結果がNOのときには、ステップ119に進み、モデルパラメータa1が、判定値X_A2STABから下限値X_IDA2Lを減算した値(X_A2STAB−X_IDA2L)以下であるか否かを判別する。この判別結果がYESのときには、ステップ120に進み、モデルパラメータa2を、判定値X_A2STABからモデルパラメータa1の絶対値を減算した値(X_A2STAB−|a1|)に設定すると同時に、モデルパラメータa2の安定化を実行したことを表すために、フラグF_A2STABを「1」にセットした後、本処理を終了する。
【0147】
一方、ステップ119の判別結果がNOで、a1>(X_A2STAB−X_IDA2L)のときには、ステップ121に進み、制御系を安定化させるために、モデルパラメータa1を、判定値X_A2STABから下限値X_IDA2Lを減算した値(X_A2STAB−X_IDA2L)に設定し、モデルパラメータa2を下限値X_IDA2Lに設定する。これと同時に、モデルパラメータa1,a2の安定化を実行したことを表すために、フラグF_A1STAB,F_A2STABをいずれも「1」にセットする。その後、本処理を終了する。
【0148】
前述したように、逐次型の同定アルゴリズムでは、制御対象の入出力が定常状態になると、自己励起条件の不足化に起因して、同定されたモデルパラメータの絶対値が増大する、いわゆるドリフト現象が発生しやすくなることで、制御系が不安定になったり、振動状態になったりすることがある。また、その安定限界も、エンジン3の運転状態に応じて変化する。例えば、低負荷運転状態のときには、排気ガスボリュームAB_SVが小さくなることで、供給された混合気に対する排気ガスの応答遅れやむだ時間などが大きくなり、それにより、O2センサ15の出力Voutが振動状態になりやすい。
【0149】
これに対して、以上のa1'&a2'のリミット処理では、モデルパラメータa1,a2の組み合わせが、図27にハッチングで示す規制範囲内の値に収まるように設定されるとともに、この規制範囲を決定する下限値X_IDA2Lが、排気ガスボリュームAB_SVに応じて設定されるので、この規制範囲をエンジン3の運転状態の変化、すなわち制御対象の動特性の変化に伴う安定限界の変化が反映された適切な安定限界の範囲として設定することができ、そのような規制範囲内に収まるように規制されたモデルパラメータa1,a2を用いることにより、上記ドリフト現象の発生を回避でき、制御系の安定性を確保することができる。これに加えて、モデルパラメータa1,a2の組み合わせを、制御系の安定性を確保できる上記規制範囲内の値として設定することにより、モデルパラメータa1およびモデルパラメータa2を単独で規制した場合における、制御系の不安定な状態の発生を回避できる。以上により、制御系の安定性を向上させることができ、触媒後排気ガス特性を向上させることができる。
【0150】
次に、図28を参照しながら、前記ステップ102のb1'のリミット処理について説明する。同図に示すように、この処理では、ステップ130において、前記ステップ93で算出されたモデルパラメータの同定値b1'が、前記図19のステップ83で算出された下限値X_IDB1L以上であるか否かを判別する。
【0151】
この判別結果がYESで、b1'≧X_IDB1Lのときには、ステップ1311に進み、モデルパラメータの同定値b1'が、前記図19のステップ83で算出された上限値X_IDB1H以下であるか否かを判別する。この判別結果がYESで、X_IDB1L≦b1'≦X_IDB1Hのときには、ステップ132に進み、モデルパラメータb1を同定値b1'に設定した後、本処理を終了する。
【0152】
一方、ステップ131の判別結果がNOで、b1'>X_IDB1Hのときには、ステップ133に進み、モデルパラメータb1を上限値X_IDB1Hに設定すると同時に、それを表すためにフラグF_B1LMTを「1」にセットした後、本処理を終了する。
【0153】
一方、ステップ130の判別結果がNOで、b1'<X_IDB1Lのときには、ステップ134に進み、モデルパラメータb1を下限値X_IDB1Lに設定すると同時に、それを表すためにフラグF_B1LMTを「1」にセットした後、本処理を終了する。
【0154】
以上のb1'のリミット処理を実行することにより、モデルパラメータb1を、X_IDB1L以上かつX_IDB1H以下の規制範囲内の値に制限することができ、それにより、逐次型の同定アルゴリズムによるドリフト現象の発生を回避できる。さらに、前述したように、これらの上下限値X_IDB1H,X_IDB1Lが、排気ガスボリュームAB_SVに応じて設定されるので、規制範囲をエンジン3の運転状態の変化、すなわち制御対象の動特性の変化に伴う安定限界の変化が反映された適切な安定限界の範囲として、設定することができ、そのような規制範囲内に規制されたモデルパラメータb1を用いることにより、制御系の安定性を確保することができる。以上により、制御系の安定性を向上させることができ、触媒後排気ガス特性を向上させることができる。
【0155】
次に、図29を参照しながら、前述したステップ33の状態予測器22の演算処理について説明する。この処理では、まず、ステップ140において、前述した式(7)の行列要素α1,α2,βi,βjを算出する。次いで、ステップ141に進み、ステップ140で算出した行列要素α1,α2,βi,βjを式(7)に適用することにより、出力偏差VO2の予測値PREVO2を算出した後、本処理を終了する。
【0156】
次に、図30を参照しながら、前述したステップ34の制御量Uslを算出する処理について説明する。この処理では、まず、ステップ150において、前述した図12の式(38)により、予測切換関数σPREを算出する。
【0157】
次に、ステップ151に進み、予測切換関数σPREの積算値SUMSIGMAを算出する。この処理では、図31に示すように、まず、ステップ160において、下記の3つの条件(l)〜(n)のうちの少なくとも1つが成立しているか否かを判別する。
(l)適応制御フラグF_PRISMONが「1」であること。
(m)後述する積算値保持フラグF_SS_HOLDが「0」であること。
(n)後述するADSM実行フラグF_KOPRが「0」であること。
【0158】
このステップ160の判別結果がYESのとき、すなわち積算値SUMSIGMAの算出条件が成立しているときには、ステップ161に進み、積算値SUMSIGMAの今回値SUMSIGMA(k)を、前回値SUMSIGMA(k−1)に、制御周期ΔTと予測切換関数σPREとの積を加算した値[SUMSIGMA(k−1)+ΔT・σPRE]に設定する。
【0159】
次いで、ステップ162に進み、ステップ161で算出した今回値SUMSIGMA(k)が所定の下限値SUMSLより大きいか否かを判別する。この判別結果がYESのときには、ステップ162に進み、今回値SUMSIGMA(k)が所定の上限値SUMSHより小さいか否かを判別する。この判別結果がYESで、SUMSL<SUMSIGMA(k)<SUMSHのときには、そのまま本処理を終了する。
【0160】
一方、ステップ163の判別結果がNOで、SUMSIGMA(k)≧SUMSHのときには、ステップ164に進み、今回値SUMSIGMA(k)を上限値SUMSHに設定した後、本処理を終了する。一方、ステップ162の判別結果がNOで、SUMSIGMA(k)≦SUMSLのときには、ステップ165に進み、今回値SUMSIGMA(k)を下限値SUMSLに設定した後、本処理を終了する。
【0161】
一方、ステップ160の判別結果がNOのとき、すなわち3つの条件(l)〜(n)がいずれも不成立で、積算値SUMSIGMAの算出条件が不成立であるときには、ステップ166に進み、今回値SUMSIGMA(k)を前回値SUMSIGMA(k−1)に設定する。すなわち、積算値SUMSIGMAをホールドする。この後、本処理を終了する。
【0162】
図30に戻り、ステップ151に続くステップ152〜154において、前述した図12の式(40)〜(42)により、等価制御入力Ueq、到達則入力Urchおよび適応則入力Uadpをそれぞれ算出する。
【0163】
次に、ステップ155に進み、これらの等価制御入力Ueq、到達則入力Urchおよび適応則入力Uadpの和を、制御量Uslとして設定した後、本処理を終了する。
【0164】
次に、図32,33を参照しながら、前述した図15のステップ36のスライディングモード制御量DKCMDSLDの算出処理について説明する。この処理では、まず、ステップ170において、制御量Uslのリミット値算出処理を実行する。この処理では、その詳細は説明は省略するが、前述したステップ35のコントローラの安定判別処理の判別結果と、後述する制御量Uslの適応上下限値Usl_ah,Usl_alとに基づいて、非アイドル運転用の上下限値Usl_ahf,Usl_alfと、アイドル運転用の上下限値Usl_ahfi,Usl_alfiとをそれぞれ算出する。
【0165】
次いで、ステップ171に進み、アイドル運転フラグF_IDLEが「0」であるか否かを判別する。この判別結果がYESで、アイドル運転中でないときには、ステップ171に進み、前述した図30の処理で算出された制御量Uslが、非アイドル運転用の下限値Usl_alf以下であるか否かを判別する。
【0166】
この判別結果がNOで、Usl>Usl_alfのときには、ステップ173に進み、制御量Uslが非アイドル運転用の上限値Usl_ahf以上であるか否かを判別する。この判別結果がNOで、Usl_alf<Usl<Usl_ahfのときには、ステップ174に進み、スライディングモード制御量DKCMDSLDを制御量Uslに設定すると同時に、積算値保持フラグF_SS_HOLDを「0」にセットする。
【0167】
次いで、ステップ175に進み、適応下限値の今回値Usl_al(k)を、前回値Usl_al(k−1)に所定の減少側値X_AL_DECを加算した値[Usl_al(k−1)+X_AL_DEC]に設定すると同時に、適応上限値の今回値Usl_ah(k)を、前回値Usl_ah(k−1)から所定の減少側値X_AL_DECを減算した値[Usl_al(k−1)−X_AL_DEC]に設定した後、本処理を終了する。
【0168】
一方、ステップ173の判別結果がYESで、Usl≧Usl_ahfのときには、ステップ176に進み、スライディングモード制御量DKCMDSLDを非アイドル運転用の適応上限値Usl_ahfに設定すると同時に、積算値保持フラグF_SS_HOLDを「1」にセットする。
【0169】
次いで、ステップ177に進み、始動後タイマのタイマ値TMACRが所定時間X_TMAWASTより小さいこと、またはF/C後判定フラグF_AFCが「1」であることが成立しているか否かを判別する。この始動後タイマは、エンジン3の始動後の経過時間を計時するアップカウント式のタイマである。
【0170】
この判別結果がYESのとき、すなわち、エンジン始動後、所定時間X_TMAWASTが経過していないか、またはフューエルカット運転の終了後、所定時間X_TM_AFCが経過していないときには、そのまま本処理を終了する。
【0171】
一方、ステップ177の判別結果がNOのとき、すなわち、エンジン始動後、所定時間X_TMAWASTが経過し、かつフューエルカット運転の終了後、所定時間X_TM_AFCが経過したときには、ステップ178に進み、適応下限値の今回値Usl_al(k)を、前回値Usl_al(k−1)に減少側値X_AL_DECを加算した値[Usl_al(k−1)+X_AL_DEC]に設定すると同時に、適応上限値の今回値Usl_ah(k)を、前回値Usl_ah(k−1)に所定の増大側値X_AL_INCを加算した値[Usl_ah(k−1)+X_AL_INC]に設定した後、本処理を終了する。
【0172】
一方、ステップ172の判別結果がYESで、Usl≦Usl_alfのときには、ステップ179に進み、スライディングモード制御量DKCMDSLDを非アイドル運転用の適応下限値Usl_alfに設定すると同時に、積算値保持フラグF_SS_HOLDを「1」にセットする。
【0173】
次いで、ステップ180に進み、第2発進フラグF_VSTが「1」であるか否かを判別する。この判別結果がYESで、車両の発進後、第2所定時間TVSTが経過しておらず、第2発進モード中であるときには、そのまま本処理を終了する。
【0174】
一方、ステップ180の判別結果がNOで、車両の発進後、第2所定時間TVSTが経過し、第2発進モードが終了したときには、ステップ181に進み、適応下限値の今回値Usl_al(k)を、前回値Usl_al(k−1)から増大側値X_AL_INCを減算した値[Usl_al(k−1)−X_AL_INC]に設定すると同時に、適応上限値の今回値Usl_ah(k)を、前回値Usl_ah(k−1)から減少側値X_AL_DECを減算した値[Usl_ah(k−1)−X_AL_DEC]に設定する。その後、本処理を終了する。
【0175】
一方、ステップ171の判別結果がNOで、アイドル運転中であるときには、図33のステップ182に進み、制御量Uslが、アイドル運転用の下限値Usl_alfi以下であるか否かを判別する。この判別結果がNOで、Usl>Usl_alfiのときには、ステップ183に進み、制御量Uslがアイドル運転用の上限値Usl_ahfi以上であるか否かを判別する。
【0176】
この判別結果がNOで、Usl_alfi<Usl<Usl_ahfiのときには、ステップ184に進み、スライディングモード制御量DKCMDSLDを制御量Uslに設定すると同時に、積算値保持フラグF_SS_HOLDを「0」にセットした後、本処理を終了する。
【0177】
一方、ステップ183の判別結果がYESで、Usl≧Usl_ahfiのときには、ステップ185に進み、スライディングモード制御量DKCMDSLDをアイドル運転用の上限値Usl_ahfiに設定すると同時に、積算値保持フラグF_SS_HOLDを「1」にセットした後、本処理を終了する。
【0178】
一方、ステップ182の判別結果がYESで、Usl≦Usl_alfiのときには、ステップ186に進み、スライディングモード制御量DKCMDSLDをアイドル運転用の下限値Usl_alfiに設定すると同時に、積算値保持フラグF_SS_HOLDを「1」にセットした後、本処理を終了する。
【0179】
次に、図34を参照しながら、前述した図15のステップ37のΔΣ変調制御量DKCMDDSMを算出する処理について説明する。この処理は、後述する理由により、100msecの周期で実行される。同図に示すように、この処理では、まず、ステップ190において、RAMに記憶されている、前回のループで算出されたDSM信号値の今回値DSMSGNS(k)[=u''(k)]を、前回値DSMSGNS(k−1)[=u''(k−1)]として設定する。
【0180】
次に、ステップ191に進み、RAMに記憶されている、前回のループで算出された偏差積分値の今回値DSMSIGMA(k)[=σd(k)]を、前回値DSMSIGMA(k−1)[=σd(k−1)]として設定する。
【0181】
次いで、ステップ192に進み、出力偏差の予測値PREVO2(k)が値0以上であるか否かを判別する。この判別結果がYESのときには、エンジン3が混合気の空燃比をリーン側に変更すべき運転状態にあるとして、ステップ193に進み、参照信号値用のゲインKRDSM(=Gd)を、リーン化用の値KRDSMLに設定した後、後述するステップ195に進む。
【0182】
一方、ステップ192の判別結果がNOのときには、エンジン3が混合気の空燃比をリッチ側に変更すべき運転状態にあるとして、ステップ194に進み、参照信号値用のゲインKRDSMを、リーン化用の値KRDSMLよりも大きいリッチ化用の値KRDSMRに設定した後、ステップ195に進む。
【0183】
このように、リーン化用の値KRDSMLおよびリッチ化用の値KRDSMRが互いに異なる値に設定されている理由は、以下による。すなわち、混合気の空燃比をリーン側に変更する際には、第1触媒装置8aのNOx浄化率を確保すべく、リーンバイアスによるNOx排出量の抑制効果を得るために、リーン化用の値KRDSMLをリッチ化用の値KRDSMRよりも小さい値に設定することで、O2センサ15の出力Voutの目標値Vopへの収束速度がリッチ側への変更時よりも遅くなるように、空燃比を制御する。一方、混合気の空燃比をリッチ側に変更する際には、第1および第2触媒装置8a,8bのNOx浄化率を十分に回復させるため、リッチ化用の値KRDSMRをリーン化用の値KRDSMLよりも大きい値に設定することで、O2センサ15の出力Voutの目標値Vopへの収束速度がリーン側への変更時よりも速くなるように、空燃比を制御する。以上により、混合気の空燃比をリッチ側およびリーン側に変更する際、良好な触媒後排気ガス特性を確保することができる。
【0184】
ステップ193または194に続くステップ195では、値−1、参照信号値用のゲインKRDSMおよび予測値の今回値PREVO2(k)を互いに乗算した値から、上記ステップ190で算出したDSM信号値の前回値DSMSGNS(k−1)を減算した値[−1・KRDSM・PREVO2(k)−DSMSGNS(k−1)]を、偏差信号値DSMDELTA[=δ(k)]として設定する。この処理は、前述した式(27),(28)に相当する。
【0185】
次いで、ステップ196に進み、偏差積分値の今回値DSMSIGMA(k)を、ステップ191で算出した前回値DSMSIGMA(k−1)と、ステップ195で算出した偏差信号値DSMDELTAとの和[DSMSIGMA(k−1)+DSMDELTA]に設定する。この処理は、前述した式(29)に相当する。
【0186】
次に、ステップ197〜199において、ステップ196で算出した偏差積分値の今回値DSMSIGMA(k)が値0以上のときには、DSM信号値の今回値DSMSGNS(k)を値1に設定し、偏差積分値の今回値DSMSIGMA(k)が値0よりも小さいときには、DSM信号値の今回値DSMSGNS(k)を値−1に設定する。以上のステップ197〜199の処理は、前述した式(30)に相当する。
【0187】
次いで、ステップ200において、排気ガスボリュームAB_SVに応じて、図35に示すテーブルを検索することにより、DSM信号値用のゲインKDSM(=Fd)を算出する。同図に示すように、このゲインKDSMは、排気ガスボリュームAB_SVが小さいほど、より大きな値に設定されている。これは、排気ガスボリュームAB_SVが小さいほど、すなわちエンジン3の負荷が小さい状態であるほど、O2センサ15の出力Voutの応答性が低下するので、それを補償するためである。このようにゲインKSDMを設定することにより、ΔΣ変調制御量DKCMDDSMを、例えばオーバーゲイン状態などを回避しながら、エンジン3の運転状態に応じて適切に算出することができ、それにより、触媒後排気ガス特性を向上させることができる。
【0188】
なお、このゲインKDSMの算出に用いるテーブルは、ゲインKDSMが排気ガスボリュームAB_SVに応じて設定されている上記テーブルに限らず、エンジン3の負荷を表すパラメータ(例えば基本燃料噴射時間Tim)に応じてゲインKDSMが予め設定されているものであればよい。また、触媒装置8a,8bの劣化判別器が設けられている場合には、この劣化判別器で判別された触媒装置8a,8bの劣化度合が大きいほど、ゲインDSMをより小さい値に補正するようにしてもよい。
【0189】
次に、ステップ201に進み、ΔΣ変調制御量DKCMDDSMを、DSM信号値用のゲインKDSMと、DSM信号値の今回値DSMSGNS(k)とを互いに乗算した値[KDSM・DSMSGNS(k)]に設定した後、本処理を終了する。この処理が、前述した式(31)に相当する。この場合、DSMSGNS(k)は、上記ステップ198,199で、値1または値−1に設定されるので、ΔΣ変調制御量DKCMDDSMは、値KDSMまたは値−KDSMに切り換えて設定される。
【0190】
次に、図36を参照しながら、前述した図15のステップ38の適応目標空燃比KCMDSLDを算出する処理について説明する。同図に示すように、この処理では、まず、ステップ210において、アイドル運転フラグF_IDLEが「1」であること、およびアイドル時ADSM実行フラグF_SWOPRIが「1」であることがいずれも成立しているか否かを判別する。このアイドル時ADSM実行フラグF_SWOPRIは、エンジン3がアイドル運転中で、かつADSM処理を実行すべき運転状態のときに「1」に、それ以外のときに「0」にセットされる。
【0191】
この判別結果がYESのとき、すなわちエンジン3がアイドル運転中でADSM処理により適応目標空燃比KCMDSLDを算出すべき運転状態のときには、ステップ211に進み、適応目標空燃比KCMDSLDを、基準値FLAFBASEにΔΣ変調制御量DKCMDDSMを加算した値[FLAFBASE+DKCMDDSM]に設定する。この処理が、前述した式(32)に相当する。この場合、FLAFBASEは一定値であるので、目標空燃比KCMDは、ΔΣ変調制御量DKCMDDSM分だけ変化するとともに、このΔΣ変調制御量DKCMDDSMは、前述したように値KDSMまたは値−KDSMに切り換えて設定されるので、この切換に伴い、目標空燃比KCMDは、パータベーション制御と同様に変動するように設定される。
【0192】
次いで、ステップ212に進み、ADSM処理を実行したことを表すために、ADSM実行済みフラグF_KOPRを「1」に設定した後、本処理を終了する。
【0193】
一方、ステップ210の判別結果がNOのときには、ステップ213に進み、触媒/O2センサフラグF_FCATDSMが「1」であるか否かを判別する。この触媒/O2センサフラグF_FCATDSMは、以下の4つの条件(o)〜(r)のうちの少なくとも1つが成立しているときに「1」に、それ以外は「0」にセットされる。
(o)第1触媒装置8aの触媒容量が所定値以上であること。
(p)第1触媒装置8aの貴金属含有量が所定値以上であること。
(q)LAFセンサ14がエンジン3の排気管7に設けられていないこと。
(r)O2センサ15が第2触媒装置8bよりも下流に設けられていること。
【0194】
この判別結果がYESのときには、ステップ214に進み、第1発進フラグF_VOTVST、および発進後ADSM実行フラグF_SWOPRVSTがいずれも「1」であるか否かを判別する。この発進後ADSM実行フラグF_SWOPRVSTは、車両の発進後で、かつエンジン3がADSM処理を実行すべき運転状態のときに「1」に、それ以外のときに「0」にセットされる。
【0195】
この判別結果がYESのとき、すなわち車両の発進後、第1所定時間TVOTVSTが経過し、かつADSM処理を実行すべき運転状態のときには、前述したように、ステップ211,212を実行した後、本処理を終了する。
【0196】
一方、ステップ214の判別結果がNOのときには、ステップ215に進み、排気ガスボリュームAB_SVが所定値OPRSVH以下であること、および小排気時ADSM実行フラグF_SWOPRSVが「1」であることがいずれも成立しているか否かを判別する。小排気時ADSM実行フラグF_SWOPRSVは、エンジン3の排気ガスボリュームAB_SVが小さい状態で、かつエンジン3がADSM処理を実行すべき運転状態のときに「1」に、それ以外のときに「0」にセットされる。
【0197】
この判別結果がYESのとき、すなわち排気ガスボリュームAB_SVが小さく、かつエンジン3がADSM処理を実行すべき運転状態のときには、前述したように、ステップ211,212を実行した後、本処理を終了する。
【0198】
一方、ステップ215の判別結果がNOのときには、エンジン3がPRISM処理を実行すべき運転状態であるとして、ステップ216に進み、適応目標空燃比KCMDSLDを、基準値FLAFBASEに適応補正項FLAFADPおよびスライディングモード制御量DKCMDSLDを加算した値[FLAFBASE+FLAFADP+DKCMDSLD]に設定する。次いで、ステップ217に進み、PRISM処理を実行したことを表すために、ADSM実行済みフラグF_KOPRを「0」にセットした後、本処理を終了する。
【0199】
一方、前記ステップ213の判別結果がNOのとき、すなわち前述した4つの条件(o)〜(r)がいずれも成立していないときには、ステップ214,215をスキップし、前述したステップ216,217を実行した後、本処理を終了する。以上のように、この適応目標空燃比KCMDSLDの算出処理では、適応目標空燃比KCMDSLDが、エンジン3の運転状態に応じて、ADSM処理またはPRISM処理に切り換えて算出される。
【0200】
次に、図37を参照しながら、以上の適応目標空燃比KCMDSLDの算出処理において、適応目標空燃比KCMDSLDすなわち目標空燃比KCMDがADSM処理で算出された場合の空燃比制御の動作の一例について説明する。
【0201】
同図に示すように、目標空燃比KCMDがADSM処理で算出された場合、前述したように、パータベーション制御と同様に、目標空燃比KCMDが制御される。具体的には、O2センサ15の出力Voutが目標値Vopから離れているときには、目標空燃比KCMDは、比較的低い周波数例えば1Hz以下の周波数で、かつゲインKDSMの振幅で変動するように制御される。一方、出力Voutが目標値Vop付近に近づくと、上記よりも高い5Hzの周波数でかつゲインKDSMを振幅として変動するように自動的に制御される。これは、後述するように、出力Voutが目標値Vop付近にある場合には、出力偏差VO2が値0に近くなることで、制御サイクル毎に、前述したステップ197の判別結果がYESとNOの間で切り換わることによる。
【0202】
以上のように目標空燃比KCMDの変動周波数を制御する理由について、図38を参照しながら説明する。同図は、O2センサ15の出力Voutが目標値Vop付近にいる場合に、目標空燃比KCMDを強制的に正弦波状に加振したときの第1触媒装置8aによる排気ガス浄化率の測定結果を示しており、破線で示すデータは、第1触媒装置8aが新品で未劣化状態の場合の測定結果であり、実線で示すデータは、第1触媒装置8aが劣化状態の場合の測定結果である。
【0203】
同図に示すように、排気ガス浄化率は、第1触媒装置8aが未劣化の場合には、目標空燃比KCMDの変動周波数にかかわらず、良好な値を示す。一方、第1触媒装置8aが劣化している場合には、排気ガス浄化率は、3Hzよりも低い周波数域で大幅に悪化するとともに、3Hz以上、より好ましくは5Hz以上の周波数域で良好な値を示すことが確認できた。したがって、目標空燃比KCMDの変動周波数を5Hz以上に制御することにより、排気ガス浄化率を良好な状態に保持できることが判る。
【0204】
本実施形態では、前述したステップ197の判別結果に基づき、目標空燃比KCMDが、基準値FLAFBASEを中心としてゲインKDSM分だけ変動するように制御されるため、O2センサ15の出力Voutが目標値Vop付近にあるときには、ステップ197の判別結果が制御サイクル毎に変化し、その結果、図34のΔΣ変調制御量DKCMDDSMの算出処理が2制御サイクル分、実行されたときに、1周期分の目標空燃比KCMDの変動波形が得られる。したがって、排気ガス浄化率を良好な状態に保持するには、図34の算出処理を、5×2=10Hz以上の周波数、すなわち100msec以下の周期で実行すればよい。この理由から、本実施形態では、前述した図14,15の制御処理および図34の算出処理の実行周期が100msecに設定されている。
【0205】
なお、同図に示すように、目標空燃比KCMDがPRISM処理で算出された場合、目標空燃比KCMDは、ほとんど変動せず、その変動周波数はほぼ値0になり、その際、排気ガス浄化率は、特例的に、目標空燃比KCMDの変動周波数が5Hz以上の場合と同レベルに保たれる。また、振動波形の入力を用いる制御手法として、PWM(Pulse wave Modulation)制御が知られているけれども、この制御手法は、振動周期が一定で、入力の振幅を変化させるものに過ぎないため、この制御手法により、目標空燃比KCMDの変動周波数を本実施形態のように制御することは不可能である。
【0206】
次に、図39を参照しながら、図15のステップ39の適応補正項FLAFADPの算出処理について説明する。同図に示すように、この処理では、まず、ステップ220において、出力偏差VO2が所定の範囲(ADL<VO2<ADH)内の値であるか否かを判別する。この判別結果がYESのとき、すなわち出力偏差VO2が小さく、O2センサ15の出力Voutが目標値Vopの近傍にあるときには、ステップ221に進み、適応則入力Uadpが所定の下限値NRLより小さいか否かを判別する。
【0207】
この判別結果がNOで、Uadp≧NRLのときには、ステップ222に進み、適応則入力Uadpが所定の上限値NRHより大きいか否かを判別する。この判別結果がNOで、NRL≦Uadp≦NRHのときには、ステップ223に進み、適応補正項の今回値FLAFADP(k)を前回値FLAFADP(k−1)に設定する。すなわち、適応補正項FLAFADPの値をホールドする。この後、本処理を終了する。
【0208】
一方、ステップ222の判別結果がYESで、Uadp>NRHのときには、ステップ224に進み、適応補正項の今回値FLAFADP(k)を、前回値FLAFADP(k−1)に所定の更新値X_FLAFDLTを加算した値[FLAFADP(k−1)+X_FLAFDLT]に設定した後、本処理を終了する。
【0209】
一方、ステップ221の判別結果がYESで、Uadp<NRLのときには、ステップ225に進み、適応補正項の今回値FLAFADP(k)を、前回値FLAFADP(k−1)から所定の更新値X_FLAFDLTを減算した値[FLAFADP(k−1)−X_FLAFDLT]に設定した後、本処理を終了する。
【0210】
以上のように、第1実施形態の空燃比制御装置1によれば、目標空燃比KCMDを制御入力とし、O2センサ15の出力Voutを出力とする制御対象において、目標空燃比KCMDを周期的に変動させるように制御する場合、すなわちパータベーション制御と同様の制御を行う場合、O2センサ15の出力Voutが目標値Vopに付近にあるときには、目標空燃比KCMDの変動周波数が、触媒の良好な排気ガス浄化率が得られる5Hz以上の周波数になるように制御される。したがって、第1触媒装置8aの劣化の有無にかかわらず、排気ガス浄化率を良好に保持することができ、それにより、触媒後排気ガス特性を向上させることができる。また、パータベーション制御と同様の制御が、O2センサ15の出力Voutが目標値Vop付近にあるとき、すなわち第1触媒装置8aによる良好な排気ガス浄化率が得られるときに実行されるので、排気ガス浄化率をさらに向上させることができる。
【0211】
一方、O2センサ15の出力Voutが目標値Vopから離れていて、排気ガス浄化率の低下を招く状態にあるときには、目標空燃比KCMDが上記周波数よりも低い1Hz以下の周波数で変動するように制御されるので、出力Voutを迅速に目標値Vop付近に近づけることができ、排気ガス浄化率を良好な状態に迅速に復帰させることができる。さらに、ΔΣ変調アルゴリズムを適用したADSMコントローラ20を用いることにより、前述したように、O2センサ15の出力Voutが目標値Vop付近に近づくと、エンジン3の運転状態などとは無関係に、目標空燃比KCMDの変動周波数をより高い値に自動的に変化させることができる。これにより、目標空燃比KCMDと目標値Vopとの比較結果に基づいて、目標空燃比KCMDの変動周波数を切り換えるようなプログラムなどを追加することなく、上記のように空燃比を制御することができる。
【0212】
また、オンボード同定器23により、モデルパラメータa1,a2,b1が逐次同定され、状態予測器により、予測時間dtと、オンボード同定器23で逐次同定されたモデルパラメータa1,a2,b1とを用いた制御対象モデルに基づいて、予測値PREVO2が算出され、この予測値PREVO2を用いて、DSMコントローラ24により目標空燃比KCMDが算出されるので、制御対象の入出力間での制御タイミングのずれを適切に解消することができる。それにより、排気ガス浄化率をさらに向上させることができ、触媒後排気ガス特性をさらに向上させることができる。
【0213】
また、目標空燃比KCMDの変動の振幅すなわちDSM信号値用のゲインKDSMが、排気ガスボリュームAB_SVに応じて設定されるので、O2センサ15の出力Voutの応答性が排気ガスボリュームAB_SVの変化に伴って変化するのを補償しながら、目標空燃比KCMDの振幅を設定することができる。これにより、目標空燃比KCMDの振幅を、排気ガスボリュームAB_SVすなわちエンジン負荷の変化に伴うオーバーゲイン状態などを回避しながら、適切に設定することができ、その結果、良好な排気ガス浄化率を確保することができる。
【0214】
以下、本発明の第2〜第8実施形態に係る空燃比制御装置について説明する。なお、以下の各実施形態の説明では、上述した第1実施形態と同じまたは同等の構成要素については、同一の参照番号を付し、その説明は適宜、省略するものとする。
【0215】
まず、図40を参照しながら、第2実施形態の空燃比制御装置について説明する。同図に示すように、この第2実施形態の空燃比制御装置1は、第1実施形態の空燃比制御装置1と比べて、オンボード同定器23のみが異なっている。具体的には、第1実施形態のオンボード同定器23では、KACT、Voutおよびφop(KCMD)に基づいて、モデルパラメータa1,a2,b1が算出されるのに対して、本実施形態のオンボード同定器23では、Voutおよびφopに基づいて、モデルパラメータa1,a2,b1が算出される。
【0216】
すなわち、このオンボード同定器23では、第1実施形態の図6の式(16)〜(23)に示す同定アルゴリズムに代えて、前述した図5の式(8)〜(15)に示す同定アルゴリズムにより、モデルパラメータの同定値a1',a2',b1'が算出されるとともに、これらに前述した図26,28のリミット処理を施すことにより、モデルパラメータa1,a2,b1が算出される。このオンボード同定器23の演算処理の具体的なプログラムは、図示しないが、第1実施形態ものとほぼ同様に構成される。以上のような本実施形態の空燃比制御装置1によれば、第1実施形態の空燃比制御装置1と同様の効果を得ることができる。
【0217】
次に、図41を参照しながら、第3実施形態の空燃比制御装置について説明する。同図に示すように、この第3実施形態の空燃比制御装置1は、第1実施形態の空燃比制御装置1と比べて、状態予測器22のみが異なっている。具体的には、第1実施形態の状態予測器22では、a1、a2、b1、KACT、Voutおよびφop(KCMD)に基づいて、予測値PREVO2が算出されるのに対して、本実施形態のオンボード同定器23では、a1、a2、b1、Voutおよびφopに基づいて、予測値PREVO2が算出される。
【0218】
すなわち、この状態予測器22では、第1実施形態の図4の式(7)に示す予測アルゴリズムに代えて、同図の式(6)に示す予測アルゴリズムにより、出力偏差VO2の予測値PREVO2が算出される。この状態予測器22の演算処理の具体的なプログラムは、図示しないが、第1実施形態のものとほぼ同様に構成される。この空燃比制御装置1によれば、第1実施形態の空燃比制御装置1と同様の効果を得ることができる。
【0219】
次に、図42を参照しながら、第4実施形態の空燃比制御装置について説明する。同図に示すように、この第4実施形態の空燃比制御装置1は、第1実施形態の空燃比制御装置1と比べると、ADSMコントローラ20、PRISMコントローラ21およびオンボード同定器23に代えて、スケジュール型DSMコントローラ20A(目標空燃比設定手段)、スケジュール型状態予測スライディングモードコントローラ21Aおよびパラメータスケジューラ28(モデルパラメータ設定手段)を用いることで、モデルパラメータa1,a2,b1を算出する点のみが異なっている。
【0220】
このパラメータスケジューラ28では、まず、前述した式(44)により、エンジン回転数NEおよび吸気管内絶対圧PBAに基づいて、排気ガスボリュームAB_SVが算出される。次いで、図43に示すテーブルにより、排気ガスボリュームAB_SVに応じて、モデルパラメータa1,a2,b1が算出される。
【0221】
このテーブルでは、モデルパラメータa1は、排気ガスボリュームAB_SVが大きいほど、より小さい値に設定されており、これとは逆に、モデルパラメータa2,b1は、排気ガスボリュームAB_SVが大きいほど、より大きい値に設定されている。これは、排気ガスボリュームAB_SVの増大に伴い、制御対象の出力すなわちO2センサ15の出力Voutが安定化する一方、排気ガスボリュームAB_SVの減少に伴い、O2センサ15の出力Voutが振動的になることによる。
【0222】
スケジュール型DSMコントローラ20Aは、以上のように算出されたモデルパラメータa1,a2,b1を用い、前述した第1実施形態と同様のDSMコントローラ24により目標空燃比KCMDを算出する。また、スケジュール型状態予測スライディングモードコントローラ21Aも、以上のように算出されたモデルパラメータa1,a2,b1を用い、前述した第1実施形態と同様のSLDコントローラ25により目標空燃比KCMDを算出する。
【0223】
この空燃比制御装置1によれば、第1実施形態の空燃比制御装置1と同様の効果を得ることができる。これに加えて、パラメータスケジューラ28を用いることにより、オンボード同定器23を用いる場合と比べて、モデルパラメータa1,a2,b1を、より迅速に算出することができる。これにより、制御の応答性を向上させることができ、良好な触媒後排気ガス特性をより迅速に確保することができる。
【0224】
次に、図44を参照しながら、第5実施形態の空燃比制御装置について説明する。この第5実施形態の空燃比制御装置1は、第1実施形態の空燃比制御装置1のDSMコントローラ24に代えて、SDMコントローラ29を用いる点のみが異なっている。このSDMコントローラ29は、ΣΔ変調アルゴリズムを適用した制御アルゴリズムにより、予測値PREVO2(k)に基づいて、制御入力φop(k)を算出するものである。
【0225】
すなわち、同図に示すように、このSDMコントローラ29では、反転増幅器29aにより、参照信号r(k)が、値−1、参照信号用のゲインGdおよび予測値PREVO2(k)を互いに乗算した信号として生成される。次に、積分器29bにより、参照信号積分値σdr(k)が、遅延素子29cで遅延された参照信号積分値σdr(k−1)と参照信号r(k)との和の信号として生成される。一方、積分器29dにより、SDM信号積分値σdu(k)が、遅延素子29eで遅延されたSDM信号積分値σdu(k−1)と、遅延素子29jで遅延されたSDM信号u''(k−1)との和の信号として生成される。そして、差分器29fにより、参照信号積分値σdr(k)とSDM信号積分値σdu(k)との偏差信号δ''(k)が生成される。
【0226】
次いで、量子化器29g(符号関数)により、SDM信号u''(k)が、この偏差信号δ''(k)を符号化した値として生成される。そして、増幅器29hにより、増幅SDM信号u(k)がSDM信号u''(k)を所定のゲインFdで増幅した値として生成され、次に、加算器29iにより、この増幅SDM信号u(k)を所定の基準値FLAFBASEに加算した値として、制御入力φop(k)が生成される。
【0227】
以上のSDMコントローラ29の制御アルゴリズムは、以下の式(45)〜(51)で表される。
r(k)=−1・Gd・PREVO2(k) ……(45)
σdr(k)=σdr(k−1)+r(k) ……(46)
σdu(k)=σdu(k−1)+u''(k−1) ……(47)
δ''(k)=σdr(k)−σdu(k) ……(48)
u''(k)=sgn(δ''(k)) ……(49)
u(k)=Fd・u''(k) ……(50)
φop(k)=FLAFBASE+u(k) ……(51)
ここで、Gd,Fdはゲインを表す。また、符号関数sgn(δ''(k))の値は、δ''(k)≧0のときにはsgn(δ''(k))=1となり、δ''(k)<0のときにはsgn(δ''(k))=−1となる(なお、δ''(k)=0のときに、sgn(δ''(k))=0と設定してもよい)。
【0228】
上記のように、このSDMコントローラ29に適用したΣΔ変調アルゴリズムでは、SDM信号u''(k)は、偏差信号δ''(k)が値0以上のときに値1に、偏差積分値σd(k)が値0よりも小さいときに値−1にそれぞれ設定される。
【0229】
以上のSDMコントローラ29の制御アルゴリズムにおけるΣΔ変調アルゴリズムの特性は、ΔΣ変調アルゴリズムと同様に、SDM信号u(k)を、これを制御対象に入力した際、参照信号r(k)が制御対象の出力に再現されるような値として、生成(算出)できるという点にある。すなわち、SDMコントローラ29は、前述したDSMコントローラ24と同様の制御入力φop(k)を生成できるという特性を備えている。したがって、このSDMコントローラ29を用いる本実施形態の空燃比制御装置1によれば、第1実施形態の空燃比制御装置1と同様の効果を得ることができる。なお、SDMコントローラ29の具体的なプログラムは図示しないが、DSMコントローラ24とほぼ同様に構成される。
【0230】
次に、図45を参照しながら、第6実施形態の空燃比制御装置について説明する。この第6実施形態の空燃比制御装置1は、第1実施形態の空燃比制御装置1のDSMコントローラ24に代えて、DMコントローラ30を用いる点のみが異なっている。このDMコントローラ30は、Δ変調アルゴリズムを適用した制御アルゴリズムにより、予測値PREVO2(k)に基づいて、制御入力φop(k)を算出するものである。
【0231】
すなわち、同図に示すように、このDMコントローラ30では、反転増幅器30aにより、参照信号r(k)が、値−1、参照信号用のゲインGdおよび予測値PREVO2(k)を互いに乗算した信号として生成される。一方、積分器30bにより、DM信号積分値σdu(k)が、遅延素子30cで遅延されたDM信号積分値σdu(k−1)と、遅延素子30hで遅延されたDM信号u''(k−1)との和の信号として生成される。そして、差分器30dにより、参照信号r(k)とDM信号積分値σdu(k)との偏差信号δ''(k)が生成される。
【0232】
次いで、量子化器30e(符号関数)により、DM信号u''(k)が、この偏差信号δ''(k)を符号化した値として生成される。そして、増幅器30fにより、増幅DM信号u(k)がDM信号u''(k)を所定のゲインFdで増幅した値として生成され、次に、加算器30gにより、この増幅DM信号u(k)を所定の基準値FLAFBASEに加算した値として、制御入力φop(k)が生成される。
【0233】
以上のDMコントローラ30の制御アルゴリズムは、以下の式(52)〜(57)で表される。
r(k)=−1・Gd・PREVO2(k) ……(52)
σdu(k)=σdu(k−1)+u''(k−1) ……(53)
δ''(k)=r(k)−σdu(k) ……(54)
u''(k)=sgn(δ''(k)) ……(55)
u(k)=Fd・u''(k) ……(56)
φop(k)=FLAFBASE+u(k) ……(57)
ここで、Gd,Fdはゲインを表す。また、符号関数sgn(δ''(k))の値は、δ''(k)≧0のときにはsgn(δ''(k))=1となり、δ''(k)<0のときにはsgn(δ''(k))=−1となる(なお、δ''(k)=0のときに、sgn(δ''(k))=0と設定してもよい)。
【0234】
上記のように、このDMコントローラ29に適用したΣΔ変調アルゴリズムでは、DM信号u''(k)は、偏差信号δ''(k)が値0以上のときに値1に、偏差積分値σd(k)が値0よりも小さいときに値−1にそれぞれ設定される。
【0235】
以上のDMコントローラ30の制御アルゴリズムすなわちΔ変調アルゴリズムの特性は、ΔΣ変調アルゴリズムおよびΣΔ変調アルゴリズムと同様に、DM信号u(k)を制御対象に入力した際、参照信号r(k)が制御対象の出力に再現されるような値として、DM信号u(k)を生成(算出)できるという点にある。すなわち、DMコントローラ30は、前述したDSMコントローラ24およびSDMコントローラ29と同様の制御入力φop(k)を生成できるという特性を備えている。したがって、このDMコントローラ30を用いる本実施形態の空燃比制御装置1によれば、第1実施形態の空燃比制御装置1と同様の効果を得ることができる。なお、DMコントローラ30の具体的なプログラムは図示しないが、DSMコントローラ24とほぼ同様に構成される。
【0236】
次に、図46および図47を参照しながら、第7実施形態の空燃比制御装置について説明する。図46に示すように、この第7実施形態の空燃比制御装置1は、第1実施形態の空燃比制御装置1と比べて、LAFセンサ14がエンジン3に設けられていないとともに、O2センサ15が第2触媒装置8bよりも下流側に設けられている点のみが異なっている。
【0237】
また、LAFセンサ14を備えていないため、この空燃比制御装置1では、図47に示すように、オンボード同定器23により、O2センサ15の出力Voutおよび制御入力φop(目標空燃比KCMD)に基づいて、モデルパラメータa1,a2,b1が算出される。すなわち、このオンボード同定器23では、前述した図5の式(8)〜(15)に示す同定アルゴリズムにより、モデルパラメータの同定値a1',a2',b1'が算出されるとともに、これらに前述したリミット処理を施すことにより、モデルパラメータa1,a2,b1が算出される。
【0238】
さらに、状態予測器22により、モデルパラメータa1,a2,b1、O2センサ15の出力Voutおよび制御入力φopに基づいて、出力偏差VO2の予測値PREVO2が算出される。すなわち、この状態予測器22では、図4の式(6)に示す予測アルゴリズムにより、出力偏差VO2の予測値PREVO2が算出される。なお、これらの状態予測器22およびオンボード同定器23の演算処理の具体的なプログラムは、図示しないが、第1実施形態のものとほぼ同様に構成され、それら以外のプログラムも、第1実施形態のものと同様に構成される。
【0239】
また、LAFセンサ14を備えていないため、この空燃比制御装置1では、図13の燃料噴射制御処理において、ステップ8,9が省略され、ステップ11において、気筒ごとの要求燃料噴射量#nTCYLが、基本燃料噴射量Tim、総補正係数KTOTALおよび補正目標空燃比KCMDMの積として算出される。
【0240】
以上のような本実施形態の空燃比制御装置1によれば、第1実施形態の空燃比制御装置1と同様の効果を得ることができる。特に、前述したように、図34のステップ192〜194において、参照信号値用のゲインKRDSMを、排気ガスをリーン側に制御する場合と、リッチ側に制御する場合とで互いに異なる値に設定し、目標空燃比KCMDの目標値Vopへの収束速度を変更することにより、本実施形態のようなO2センサ15のみで空燃比を制御する場合においても、混合気の空燃比をリッチ側およびリーン側に変更する際、良好な排気ガス浄化率を確実に得ることができ、触媒後排気ガス特性を確実に得ることができる。これに加えて、LAFセンサ14を用いることなく、良好な触媒後排気ガス特性を確保できるので、その分、製造コストを削減することができる。
【0241】
次に、図48を参照しながら、第8実施形態の空燃比制御装置について説明する。同図に示すように、この第8実施形態の空燃比制御装置1は、上記第7実施形態の空燃比制御装置1において、ADSMコントローラ20、PRISMコントローラ21およびオンボード同定器23を、前記第4実施形態のスケジュール型DSMコントローラ20A、スケジュール型状態予測スライディングモードコントローラ21Aおよびパラメータスケジューラ28に置き換えたものであり、これらのコントローラ20A,21Aおよびパラメータスケジューラ28は、第4実施形態のものと同様に構成されている。この空燃比制御装置1によれば、上記第7実施形態の空燃比制御装置1と同様の効果を得ることができる。これに加えて、パラメータスケジューラ28を用いることにより、オンボード同定器23を用いる場合と比べて、モデルパラメータa1,a2,b1を、より迅速に算出することができる。これにより、制御の応答性を向上させることができ、良好な触媒後排気ガス特性をより迅速に確保することができる。
【0242】
なお、ADSMコントローラ20およびPRISMコントローラ21を、実施形態のプログラムに代えて、電気回路により構成してもよい。
【0243】
【発明の効果】
以上のように、本発明の内燃機関の空燃比制御装置によれば、パータベーション制御を行う場合において、触媒の劣化の有無にかかわらず、排気ガス浄化率を良好に保持することができ、それにより、触媒後排気ガス特性を向上させることができる。
【図面の簡単な説明】
【図1】本発明の第1実施形態に係る空燃比制御装置およびこれを適用した内燃機関の概略構成を示す図である。
【図2】劣化状態および未劣化状態の第1触媒装置を用いた場合において、LAFセンサの出力KACTに対する、両第1触媒装置のHCおよびNOxの浄化率と、O2センサ15の出力Voutとをそれぞれ測定した結果の一例を示す図である。
【図3】第1実施形態の空燃比制御装置のADSMコントローラおよびPRISMコントローラの構成を示すブロック図である。
【図4】状態予測器の予測アルゴリズムの数式の一例を示す図である。
【図5】オンボード同定器の同定アルゴリズムの数式の一例を示す図である。
【図6】オンボード同定器の同定アルゴリズムの数式の他の一例を示す図である。
【図7】ΔΣ変調を実行するコントローラおよびこれを備えた制御系の構成を示すブロック図である。
【図8】図7の制御系の制御結果の一例を示すタイミングチャートである。
【図9】第1実施形態のADSMコントローラによる適応予測型ΔΣ変調制御の原理を説明するためのタイミングチャートである。
【図10】ADSMコントローラのうちのDSMコントローラの構成を示すブロック図である。
【図11】スライディングモード制御アルゴリズムの数式を示す図である。
【図12】PRISMコントローラのスライディングモード制御アルゴリズムの数式を示す図である。
【図13】内燃機関の燃料噴射制御処理を示すフローチャートである。
【図14】適応空燃比制御処理を示すフローチャートである。
【図15】図14の続きを示すフローチャートである。
【図16】図14のステップ21における発進判定処理を示すフローチャートである。
【図17】図14のステップ23におけるPRISM/ADSM処理の実行判定処理を示すフローチャートである。
【図18】図14のステップ24における同定器演算の実行判定処理を示すフローチャートである。
【図19】図14のステップ25における各種パラメータの算出処理を示すフローチャートである。
【図20】むだ時間CAT_DELAY,KACT_Dの算出に用いるテーブルの一例を示す図である。
【図21】重みパラメータλ1の算出に用いるテーブルの一例を示す図である。
【図22】モデルパラメータa1,a2,b1の値を制限するリミット値X_IDA2L,X_IDB1L,X_IDB1Hの算出に用いるテーブルの一例を示す図である。
【図23】フィルタ次数nの算出に用いるテーブルの一例を示す図である。
【図24】図14のステップ31における同定器の演算処理を示すフローチャートである。
【図25】図24のステップ94におけるθ(k)の安定化処理を示すフローチャートである。
【図26】図25のステップ101におけるa1'&a2'のリミット処理を示すフローチャートである。
【図27】図26の処理によりa1'&a2'の組み合わせが規制される規制範囲を示す図である。
【図28】図25のステップ102におけるb1'のリミット処理を示すフローチャートである。
【図29】図15のステップ33の状態予測器の演算処理を示すフローチャートである。
【図30】図15のステップ34の制御量Uslの算出処理を示すフローチャートである。
【図31】図30のステップ151の予測切換関数σPREの積算値算出処理を示すフローチャートである。
【図32】図15のステップ36のスライディングモード制御量DKCMDSLDの算出処理を示すフローチャートである。
【図33】図32の続きを示すフローチャートである。
【図34】図15のステップ37のΔΣ変調制御量DKCMDDSMの算出処理を示すフローチャートである。
【図35】KDSMの算出に用いるテーブルの一例を示す図である。
【図36】図15のステップ38の適応目標空燃比KCMDSLDの算出処理を示すフローチャートである。
【図37】目標空燃比KCMDがADSM処理で算出された場合の空燃比制御の動作の一例を示すタイミングチャートである。
【図38】目標空燃比KCMDの変動周波数と排気ガス浄化率の関係を示す図である。
【図39】図15のステップ39の適応補正項FLAFADPの算出処理を示すフローチャートである。
【図40】第2実施形態の空燃比制御装置の概略構成を示すブロック図である。
【図41】第3実施形態の空燃比制御装置の概略構成を示すブロック図である。
【図42】第4実施形態の空燃比制御装置の概略構成を示すブロック図である。
【図43】第4実施形態の空燃比制御装置のパラメータスケジューラにおいて、モデルパラメータの算出に用いるテーブルの一例を示す図である。
【図44】第5実施形態の空燃比制御装置のSDMコントローラの概略構成を示すブロック図である。
【図45】第6実施形態の空燃比制御装置のDMコントローラの概略構成を示すブロック図である。
【図46】第7実施形態に係る空燃比制御装置およびこれを適用した内燃機関の概略構成を示す図である。
【図47】第7実施形態の空燃比制御装置の構成を示すブロック図である。
【図48】第8実施形態の空燃比制御装置の構成を示すブロック図である。
【符号の説明】
1 空燃比制御装置
2 ECU(目標空燃比設定手段、空燃比制御手段、予測値算出手段、目標空燃比算出手段、同定手段、運転状態パラメータ検出手段、モデルパラメータ設定手段、負荷パラメータ検出手段)
3 内燃機関
7 排気管(排気通路)
8a 第1触媒装置(触媒)
8b 第2触媒装置(触媒)
11 吸気管内絶対圧センサ(運転状態パラメータ検出手段、負 荷パラメータ検出手段)
13 クランク角センサ(運転状態パラメータ検出手段、負荷パラメータ検出手段)
14 LAFセンサ(上流側空燃比センサ)
15 酸素濃度センサ(下流側空燃比センサ、空燃比センサ)
20 ADSMコントローラ(目標空燃比設定手段)
20A スケジュール型DSMコントローラ(目標空燃比設定手段)
21 DSMコントローラ(目標空燃比算出手段)
22 状態予測器(予測値算出手段、目標空燃比算出手段)
23 オンボード同定器(同定手段)
28 パラメータスケジューラ(モデルパラメータ設定手段)
NE エンジン回転数(運転状態パラメータ)
PBA 吸気管内絶対圧(運転状態パラメータ)
AB_SV 排気ガスボリューム(負荷パラメータ)
KACT LAFセンサの出力(検出信号)
DKACT LAF出力偏差(上流側空燃比センサの出力を表す値)
KCMD 目標空燃比
DKCMD 空燃比偏差(目標空燃比を表す値)
DKCMD(k+dt) 空燃比偏差の時系列データ
KDSM ゲイン(所定の振幅)
Vout 酸素濃度センサの出力(検出信号)
Vop 目標値
VO2 出力偏差(下流側空燃比センサの出力を表す値、空燃比セ ンサの出力を表す値)
VO2(k-1) 出力偏差の時系列データ
VO2(k-2) 出力偏差の時系列データ
PREVO2 出力偏差の予測値
a1 モデルパラメータ
a2 モデルパラメータ
b1 モデルパラメータ
Claims (18)
- 内燃機関の排気通路の触媒よりも上流側における排気ガスの空燃比を表す検出信号を出力する上流側空燃比センサと、
前記排気通路の前記触媒よりも下流側における排気ガスの空燃比を表す検出信号を出力する下流側空燃比センサと、
前記下流側空燃比センサの出力を所定の目標値に収束させるための目標空燃比を、所定の振幅で、かつ前記下流側空燃比センサの出力が前記所定の目標値付近にあるときにはそれ以外のときよりも高い所定の周波数で変動するように、設定する目標空燃比設定手段と、
前記上流側空燃比センサの出力に基づき、前記触媒よりも上流側における排気ガスの空燃比を前記目標空燃比設定手段により設定された目標空燃比に一致させるように、前記内燃機関に供給される混合気の空燃比を制御する空燃比制御手段と、
を備えることを特徴とする内燃機関の空燃比制御装置。 - 前記目標空燃比設定手段は、Δ変調アルゴリズム、ΔΣ変調アルゴリズムおよびΣΔ変調アルゴリズムのうちのいずれか1つの変調アルゴリズムに基づき、前記目標空燃比を設定することを特徴とする請求項1に記載の内燃機関の空燃比制御装置。
- 前記目標空燃比設定手段は、
予測アルゴリズムに基づき、前記下流側空燃比センサの出力を表す値の予測値を算出する予測値算出手段と、
前記1つの変調アルゴリズムに基づき、前記算出された予測値に応じて、前記目標空燃比を算出する目標空燃比算出手段と、
を備えることを特徴とする請求項2に記載の内燃機関の空燃比制御装置。 - 前記予測アルゴリズムは、前記下流側空燃比センサの出力を表す値と、前記目標空燃比を表す値および前記上流側空燃比センサの出力を表す値の一方とを変数とする制御対象モデルを適用したアルゴリズムであることを特徴とする請求項3に記載の内燃機関の空燃比制御装置。
- 前記目標空燃比設定手段は、
前記目標空燃比を表す値の時系列データおよび前記上流側空燃比センサの出力を表す値の時系列データの一方、および前記下流側空燃比センサの出力を表す値の時系列データを変数とする離散時間系の制御対象モデルと、前記1つの変調アルゴリズムとに基づき、前記目標空燃比を算出する目標空燃比算出手段と、
前記離散時間系の制御対象モデルのモデルパラメータを逐次同定する同定手段と、
を備えることを特徴とする請求項2に記載の内燃機関の空燃比制御装置。 - 前記目標空燃比算出手段は、前記制御対象モデルを適用した予測アルゴリズムに基づき、前記下流側空燃比センサの出力を表す値の予測値を算出するとともに、当該算出された予測値および前記1つの変調アルゴリズムに基づいて、前記目標空燃比を算出することを特徴とする請求項5に記載の内燃機関の空燃比制御装置。
- 前記内燃機関の運転状態を表す運転状態パラメータを検出する運転状態パラメータ検出手段をさらに備え、
前記目標空燃比設定手段は、
前記下流側空燃比センサの出力を表す値と、前記目標空燃比を表す値および前記上流側空燃比センサの出力を表す値の一方とを変数とする制御対象モデル、ならびに前記1つの変調アルゴリズムに基づき、前記目標空燃比を算出する目標空燃比算出手段と、
前記検出された運転状態パラメータに応じて、前記制御対象モデルのモデルパラメータを設定するモデルパラメータ設定手段と、を備えることを特徴とする請求項2に記載の内燃機関の空燃比制御装置。 - 前記目標空燃比算出手段は、前記制御対象モデルを適用した予測アルゴリズムに基づき、前記下流側空燃比センサの出力を表す値の予測値を算出するとともに、当該算出された予測値および前記1つの変調アルゴリズムに基づいて、前記目標空燃比を算出することを特徴とする請求項7に記載の内燃機関の空燃比制御装置。
- 前記内燃機関の負荷を表す負荷パラメータを検出する負荷パラメータ検出手段をさらに備え、
前記目標空燃比設定手段は、当該検出された負荷パラメータに応じて、前記所定の振幅を設定することを特徴とする請求項1ないし8のいずれかに記載の内燃機関の空燃比制御装置。 - 内燃機関の排気通路の触媒よりも下流側における排気ガスの空燃比を表す検出信号を出力する空燃比センサと、
前記空燃比センサの出力を所定の目標値に収束させるための目標空燃比を、所定の振幅で、かつ前記空燃比センサの出力が所定の目標値付近にあるときにはそれ以外のときよりも高い所定の周波数域の周波数で変動するように、設定する目標空燃比設定手段と、
当該設定された目標空燃比に応じて、前記内燃機関に供給される混合気の空燃比を制御する空燃比制御手段と、
を備えることを特徴とする内燃機関の空燃比制御装置。 - 前記目標空燃比設定手段は、Δ変調アルゴリズム、ΔΣ変調アルゴリズムおよびΣΔ変調アルゴリズムのうちのいずれか1つの変調アルゴリズムに基づき、前記目標空燃比を設定することを特徴とする請求項10に記載の内燃機関の空燃比制御装置。
- 前記目標空燃比設定手段は、
予測アルゴリズムに基づき、前記空燃比センサの出力を表す値の予測値を算出する予測値算出手段と、
前記1つの変調アルゴリズムに基づき、前記算出された予測値に応じて、前記目標空燃比を算出する目標空燃比算出手段と、
を備えることを特徴とする請求項11に記載の内燃機関の空燃比制御装置。 - 前記予測アルゴリズムは、前記空燃比センサの出力を表す値と、前記目標空燃比を表す値とを変数とする制御対象モデルを適用したアルゴリズムであることを特徴とする請求項12に記載の内燃機関の空燃比制御装置。
- 前記目標空燃比設定手段は、
前記目標空燃比を表す値の時系列データおよび前記空燃比センサの出力を表す値の時系列データを変数とする離散時間系の制御対象モデルと、前記1つの変調アルゴリズムとに基づき、前記目標空燃比を算出する目標空燃比算出手段と、
前記離散時間系の制御対象モデルのモデルパラメータを逐次同定する同定手段と、
を備えることを特徴とする請求項11に記載の内燃機関の空燃比制御装置。 - 前記目標空燃比算出手段は、前記制御対象モデルを適用した予測アルゴリズムに基づき、前記空燃比センサの出力を表す値の予測値を算出するとともに、当該算出された予測値および前記1つの変調アルゴリズムに基づいて、前記目標空燃比を算出することを特徴とする請求項14に記載の内燃機関の空燃比制御装置。
- 前記内燃機関の運転状態を表す運転状態パラメータを検出する運転状態パラメータ検出手段をさらに備え、
前記目標空燃比設定手段は、
前記空燃比センサの出力を表す値および前記目標空燃比を表す値を変数とする制御対象モデルと、前記1つの変調アルゴリズムとに基づき、前記目標空燃比を算出する目標空燃比算出手段と、
前記検出された運転状態パラメータに応じて、前記制御対象モデルのモデルパラメータを設定するモデルパラメータ設定手段と、
を備えることを特徴とする請求項11に記載の内燃機関の空燃比制御装置。 - 前記目標空燃比算出手段は、前記制御対象モデルを適用した予測アルゴリズムに基づき、前記空燃比センサの出力を表す値の予測値を算出するとともに、当該算出された予測値および前記1つの変調アルゴリズムに基づいて、前記目標空燃比を算出することを特徴とする請求項16に記載の内燃機関の空燃比制御装置。
- 前記内燃機関の負荷を表す負荷パラメータを検出する負荷パラメータ検出手段をさらに備え、
前記目標空燃比設定手段は、当該検出された負荷パラメータに応じて、前記所定の振幅を設定することを特徴とする請求項10ないし17のいずれかに記載の内燃機関の空燃比制御装置。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002012854A JP3880861B2 (ja) | 2002-01-22 | 2002-01-22 | 内燃機関の空燃比制御装置 |
US10/339,601 US7059115B2 (en) | 2002-01-22 | 2003-01-10 | Air/fuel ratio control apparatus and method for internal combustion engine and engine control unit |
DE10302126A DE10302126B4 (de) | 2002-01-22 | 2003-01-21 | Luft/Kraftstoff-Verhältnisregelvorrichtung und -verfahren für eine Brennkraftmaschine |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002012854A JP3880861B2 (ja) | 2002-01-22 | 2002-01-22 | 内燃機関の空燃比制御装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2003214227A JP2003214227A (ja) | 2003-07-30 |
JP3880861B2 true JP3880861B2 (ja) | 2007-02-14 |
Family
ID=27606055
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002012854A Expired - Fee Related JP3880861B2 (ja) | 2002-01-22 | 2002-01-22 | 内燃機関の空燃比制御装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US7059115B2 (ja) |
JP (1) | JP3880861B2 (ja) |
DE (1) | DE10302126B4 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2075662A1 (en) | 2007-12-25 | 2009-07-01 | Honda Motor Co., Ltd. | Control apparatus and method |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4121914B2 (ja) * | 2003-08-08 | 2008-07-23 | 本田技研工業株式会社 | 制御装置 |
DE102004046874A1 (de) * | 2004-09-28 | 2006-04-13 | Robert Bosch Gmbh | Verfahren zum Betreiben eines Verwaltungssystems von Funktionsmodulen |
US7216638B1 (en) | 2006-07-06 | 2007-05-15 | Brunswick Corporation | Control of exhaust gas stoichiometry with inducted secondary air flow |
JP2008123549A (ja) * | 2008-02-04 | 2008-05-29 | Honda Motor Co Ltd | 制御装置 |
JP4630354B2 (ja) * | 2008-05-28 | 2011-02-09 | 本田技研工業株式会社 | Δς変調アルゴリズムを用いてプラントを制御する制御装置 |
JP5616274B2 (ja) * | 2011-03-31 | 2014-10-29 | 本田技研工業株式会社 | 空燃比制御装置 |
US8165787B2 (en) * | 2011-04-08 | 2012-04-24 | Ford Global Technologies, Llc | Method for adjusting engine air-fuel ratio |
GB201119036D0 (en) * | 2011-11-03 | 2011-12-14 | Univ Oxford Brookes | A method of controlling a dynamic physical system |
JP5883140B2 (ja) * | 2012-07-17 | 2016-03-09 | 本田技研工業株式会社 | 内燃機関の制御装置 |
JP6468255B2 (ja) * | 2016-07-04 | 2019-02-13 | トヨタ自動車株式会社 | 排気センサ |
US10947910B2 (en) * | 2019-05-07 | 2021-03-16 | Ford Global Technologies, Llc | Method and system for catalyst feedback control |
US11624333B2 (en) | 2021-04-20 | 2023-04-11 | Kohler Co. | Exhaust safety system for an engine |
CN114320626B (zh) * | 2022-01-12 | 2022-11-01 | 山东交通学院 | 一种天燃气发动机用混合器控制方法及系统 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4024706A (en) * | 1976-01-07 | 1977-05-24 | Ford Motor Company | Method of improving the operational capacity of three-way catalysts |
JPS6466441A (en) * | 1987-09-05 | 1989-03-13 | Toyota Motor Corp | Air-fuel ratio control device of internal combustion engine |
JP3306930B2 (ja) * | 1992-07-03 | 2002-07-24 | 株式会社デンソー | 内燃機関の空燃比制御装置 |
US5622047A (en) * | 1992-07-03 | 1997-04-22 | Nippondenso Co., Ltd. | Method and apparatus for detecting saturation gas amount absorbed by catalytic converter |
US5678402A (en) * | 1994-03-23 | 1997-10-21 | Honda Giken Kogyo Kabushiki Kaisha | Air-fuel ratio control system for internal combustion engines and exhaust system temperature-estimating device applicable thereto |
JP3449011B2 (ja) * | 1994-05-31 | 2003-09-22 | 株式会社デンソー | 内燃機関の空燃比制御装置 |
JP3823384B2 (ja) | 1996-07-29 | 2006-09-20 | マツダ株式会社 | エンジンの触媒劣化判定装置 |
JP3331159B2 (ja) * | 1997-09-16 | 2002-10-07 | 本田技研工業株式会社 | プラントの制御装置 |
DE69917195T2 (de) * | 1998-12-17 | 2004-09-23 | Honda Giken Kogyo K.K. | Steuersystem für das Luft/Kraftstoffverhältnis einer Brennkraftmaschine |
JP3858554B2 (ja) * | 2000-02-23 | 2006-12-13 | 株式会社日立製作所 | エンジン排気浄化装置 |
JP4275330B2 (ja) * | 2001-06-28 | 2009-06-10 | 三菱電機株式会社 | 内燃機関の空燃比制御装置 |
-
2002
- 2002-01-22 JP JP2002012854A patent/JP3880861B2/ja not_active Expired - Fee Related
-
2003
- 2003-01-10 US US10/339,601 patent/US7059115B2/en not_active Expired - Lifetime
- 2003-01-21 DE DE10302126A patent/DE10302126B4/de not_active Expired - Fee Related
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2075662A1 (en) | 2007-12-25 | 2009-07-01 | Honda Motor Co., Ltd. | Control apparatus and method |
Also Published As
Publication number | Publication date |
---|---|
DE10302126A1 (de) | 2003-08-21 |
US20030140617A1 (en) | 2003-07-31 |
DE10302126B4 (de) | 2013-01-10 |
US7059115B2 (en) | 2006-06-13 |
JP2003214227A (ja) | 2003-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3904923B2 (ja) | 制御装置 | |
JP3973922B2 (ja) | 制御装置 | |
JP3922980B2 (ja) | 制御装置 | |
JP3926703B2 (ja) | 制御装置 | |
JP3880861B2 (ja) | 内燃機関の空燃比制御装置 | |
JP3998136B2 (ja) | 内燃機関の空燃比制御装置 | |
JP4028334B2 (ja) | 制御装置 | |
JP2003003892A (ja) | 内燃機関の空燃比制御装置 | |
JP4277958B2 (ja) | 内燃機関の空燃比制御装置 | |
JP4439508B2 (ja) | 制御装置 | |
JP4277959B2 (ja) | 制御装置 | |
JP4240325B2 (ja) | 制御装置 | |
EP1279820B1 (en) | Control apparatus, control method, and engine control unit | |
JP2008169849A (ja) | 内燃機関の制御装置 | |
JP2003074399A (ja) | 内燃機関の空燃比制御装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050208 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050405 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060207 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20061010 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20061108 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 3880861 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091117 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101117 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101117 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111117 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111117 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121117 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131117 Year of fee payment: 7 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |