JP3748418B2 - Air-fuel ratio control device for internal combustion engine - Google Patents
Air-fuel ratio control device for internal combustion engine Download PDFInfo
- Publication number
- JP3748418B2 JP3748418B2 JP2002125424A JP2002125424A JP3748418B2 JP 3748418 B2 JP3748418 B2 JP 3748418B2 JP 2002125424 A JP2002125424 A JP 2002125424A JP 2002125424 A JP2002125424 A JP 2002125424A JP 3748418 B2 JP3748418 B2 JP 3748418B2
- Authority
- JP
- Japan
- Prior art keywords
- fuel ratio
- air
- hat
- value
- internal combustion
- 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 - Lifetime
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A50/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather
- Y02A50/20—Air quality improvement or preservation, e.g. vehicle emission control or emission reduction by using catalytic converters
Landscapes
- Exhaust Gas After Treatment (AREA)
- Electrical Control Of Air Or Fuel Supplied To Internal-Combustion Engine (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
- Feedback Control In General (AREA)
Description
【0001】
【発明の属する技術分野】
本発明は、内燃機関の空燃比制御装置に関する。
【0002】
【従来の技術】
例えば自動車等、内燃機関の排ガスを三元触媒等の触媒装置により浄化して放出するシステムでは、触媒装置に進入する排ガスの空燃比を、触媒装置の排ガス浄化能力を良好に発揮させるように制御することが環境保護の観点から望まれている。
【0003】
このような制御を行うものとしては、従来、例えば特開平5−321721号公報に見られるように、触媒装置を通過した排ガスの酸素濃度を触媒装置の下流側に配置した排ガスセンサ(酸素濃度センサ)により検出して、その検出値が所定の適正値になるようにPID制御を用いて触媒装置の上流側の排ガスの目標空燃比を決定し、その目標空燃比に従って内燃機関を制御することで、触媒装置に進入する排ガスの空燃比を、触媒装置の良好な浄化性能を発揮できる所定のウィンドウ内に収めるようにしたものが本願出願人により提案されている。
【0004】
一方、本願発明者等のさらなる検討によって、内燃機関の運転状態や、触媒装置の経時劣化等によらずに、触媒装置の浄化能力を可能な限り最大限に発揮させるためには、触媒装置を通過した排ガス中の酸素濃度等の特定成分の濃度を高精度で安定して所定の適正値(一定値)に整定させる必要があることが判明した。そして、前述のようにPID制御を用いた従来の制御手法では、外乱や、触媒装置を含む排気系に存する無駄時間等の影響で上記のように触媒装置を通過した排ガス中の酸素濃度等を安定して高精度で所定の適正値(一定値)に整定させることが困難であることが判明した。
【0005】
このため、本願発明者等は、触媒装置の上流側から下流側にかけての排気系を連続系(詳しくは連続時間系)でモデル化し、そのモデルに基づいて触媒装置を通過した排ガス中の酸素濃度等が所定の適正値になるように触媒装置に進入する排ガスの空燃比を制御するシステムを先に考案した(例えば特願平9−67591号、特願平8−84048号参照)。
【0006】
上記のようなモデル化によって、触媒装置を含む排気系に存する無駄時間の影響を補償(排除)したり、外乱に対する制御の安定性がPID制御等に比して高い制御手法(例えばスライディングモード制御)を用いた制御システムの構築が可能となり、ひいては、排気系の空燃比制御の精度や安定性を高めることが可能となる。
【0007】
このような空燃比制御装置では、内燃機関の運転状態によらずに制御の安定性を高めることが望まれる。
【0008】
【発明が解決しようとする課題】
本発明はかかる背景に鑑み、排ガス中の特定成分の濃度が所定の目標値になるように排ガスの空燃比を制御する場合に、その制御を安定して高精度で行うことができる内燃機関の空燃比制御装置を提供することを目的とする。
【0009】
【課題を解決するための手段】
かかる目的を達成するために、本発明の内燃機関の空燃比制御装置の第1の態様は、内燃機関の排ガス中の特定成分の濃度を検出すべく該内燃機関の排気系に設けられた排ガスセンサを備え、該排ガスセンサの出力が所定の目標値になるように内燃機関の排ガスの空燃比を制御する内燃機関の空燃比制御装置において、前記内燃機関の排ガスの空燃比を規定する操作量をスライディングモード制御を用いて決定する空燃比操作量決定手段と、該空燃比操作量決定手段が制御対象とする系のモデルの設定すべきパラメータを同定する同定手段とを備え、前記空燃比操作量決定手段は、前記モデルのパラメータの値として、前記同定手段が同定した前記パラメータの同定値と所定値とを内燃機関の運転状態に応じて選択的に用い、該モデルのパラメータの値により規定される係数を含むスライディングモード制御の所定の演算式により前記操作量を決定するようにしたことを特徴とするものである。
【0010】
かかる本発明の第1の態様によれば、空燃比操作量決定手段は、スライディングモード制御により、前記モデルのパラメータの値を用いて前記操作量を決定する。この場合、このパラメータの値としては、内燃機関の運転状態に応じて、前記同定手段による同定値と、その同定値とは別の所定値とが選択的に用いられる。この結果、内燃機関の運転状態によらずに排ガスセンサの出力が所定の目標値になるようにする排ガスの空燃比の制御を安定して行うことが可能となる。
【0011】
また、本発明の内燃機関の空燃比制御装置の第2の態様は、内燃機関の排ガス中の特定成分の濃度を検出すべく該内燃機関の排気系に設けられた排ガスセンサを備え、該排ガスセンサの出力が所定の目標値になるように内燃機関の排ガスの空燃比を制御する内燃機関の空燃比制御装置において、前記内燃機関の排ガスの空燃比を規定する操作量を決定する空燃比操作量決定手段と、該空燃比操作量決定手段が制御対象とする系であって前記排ガスセンサの出力を生成する系が有する無駄時間後における該排ガスセンサの出力の推定値を示すデータを生成する推定手段と、前記制御対象の系のモデルの設定すべきパラメータを同定する同定手段とを備え、前記推定手段は、前記モデルのパラメータの値として、前記同定手段が同定した前記パラメータの同定値と所定値とを内燃機関の運転状態に応じて選択的に用い、該モデルのパラメータの値により規定される係数を含む推定用の所定の演算式により前記排ガスセンサの出力の推定値を示すデータを生成し、前記空燃比操作量決定手段は、該推定手段が生成したデータにより示される推定値が前記目標値になるように前記操作量を決定することを特徴とする。
【0012】
かかる本発明の第2の態様によれば、前記推定手段は、前記モデルのパラメータの値を用いて前記無駄時間後の排ガスセンサの出力の推定値を表すデータ生成する。この場合、このパラメータの値としては、内燃機関の運転状態に応じて、前記同定手段による同定値と、その同定値とは別の所定値とが選択的に用いられる。そして、前記空燃比操作量決定手段は、推定手段が生成したデータにより示される排ガスセンサの出力の推定値が前記目標値になるように前記操作量を決定する。この結果、内燃機関の運転状態によらずに排ガスセンサの出力が所定の目標値になるようにする排ガスの空燃比の制御を安定して行うことが可能となる。
【0013】
かかる本発明の第1及び第2の態様では、前記内燃機関の運転状態は、該内燃機関のスロットル弁が全開であるか否かの状態、内燃機関の燃料供給の停止中であるか否かの状態、内燃機関のアイドル運転中であるか否かの状態であることが好適である。
【0014】
【発明の実施の形態】
本発明の第1の実施形態を図1乃至図23を参照して説明する。
【0015】
図1は本実施形態の制御装置の全体構成をブロック図で表したものであり、同図1において、1は例えば4気筒のエンジン(内燃機関)である。このエンジン1の各気筒毎に生成される排ガスは、エンジン1の近傍で共通の排気管2に集合され、該排気管2を介して大気中に放出される。そして、排気管2には、排ガスを浄化するために、三元触媒を用いた二つの触媒装置3,4が該排気管2の上流側から順に介装されている。尚、下流側の触媒装置4はこれを省略してもよい。
【0016】
本実施形態の制御装置は、エンジン1の排気系の空燃比を制御するもので、触媒装置3の上流側(より詳しくはエンジン1の各気筒毎の排ガスの集合箇所)で排気管2に設けられた第2排ガスセンサとしての広域空燃比センサ5と、触媒装置3の下流側(触媒装置4の上流側)で排気管2に設けられた第1排ガスセンサとしてのO2センサ(酸素濃度センサ)6と、これらのセンサ5,6の出力等に基づき後述の制御処理を行う制御ユニット7とにより構成されている。尚、制御ユニット7には、前記広域空燃比センサ5やO2センサ6の出力の他に、エンジン1の動作状態を検出するための図示しない回転数センサや吸気圧センサ、冷却水温センサ等、各種のセンサの検出信号が与えられるようになっている。
【0017】
広域空燃比センサ5は、O2センサを用いて構成されたものであり、触媒装置3に進入するエンジン1の排ガスの空燃比(これは触媒装置3に進入する排ガスの酸素濃度により示され、エンジン1内で燃焼させる混合気の空燃比に相当する)に応じたレベルの出力を生成する。この場合、広域空燃比センサ5(以下、LAFセンサ5と称する)は、該センサ5を構成するO2センサの出力から図示しないリニアライザ等の検出回路によって、触媒装置3に進入する排ガスの空燃比の広範囲にわたって、それに比例したレベルの出力KACT、すなわち、該排ガスの空燃比の検出値を示す出力KACTを生成するものである。このようなLAFセンサ5は本願出願人が特開平4−369471号公報にて詳細に開示しているので、ここではさらなる説明を省略する。
【0018】
また、触媒装置3の下流側のO2センサ6は、触媒装置3を通過した排ガス中の酸素濃度に応じたレベルの出力VO2/OUT、すなわち、該排ガス中の酸素濃度の検出値を示す出力VO2/OUTを通常的なO2センサと同様に生成する。このO2センサ6の出力VO2/OUTは、図2に示すように、触媒装置3を通過した排ガスの空燃比(酸素濃度)が理論空燃比近傍の範囲Δに存するような状態で、該排ガスの酸素濃度にほぼ比例した高感度な変化を生じるものとなる。
【0019】
制御ユニット7はマイクロコンピュータを用いて構成されたものであり、その主要な機能的構成として、エンジン1への基本燃料噴射量Timを求める基本燃料噴射量算出部8と、基本燃料噴射量Timを補正するための第1補正係数KTOTAL及び第2補正係数KCMDMをそれぞれ求める第1補正係数算出部9及び第2補正算出部10とを具備する。
【0020】
前記基本燃料噴射量算出部8は、エンジン1の回転数NEと吸気圧PBとから、それらにより規定されるエンジン1の基準の燃料噴射量をあらかじめ設定されたマップを用いて求め、その基準の燃料噴射量をエンジン1の図示しないスロットル弁の有効開口面積に応じて補正することで基本燃料噴射量Timを算出するものである。
【0021】
また、第1補正係数算出部9が求める第1補正係数KTOTALは、エンジン1の排気還流率(エンジン1の吸入空気中に含まれる排気ガスの割合)や、エンジン1の図示しないキャニスタのパージ時にエンジン1に供給される燃料のパージ量、エンジン1の冷却水温、吸気温等を考慮して前記基本燃料噴射量Timを補正するためのものである。
【0022】
また、第2補正係数算出部10が求める第2補正係数KCMDMは、後述する空燃比操作量決定部13によって決定される目標空燃比KCMDに対応してエンジン1へ流入する燃料の冷却効果による吸入空気の充填効率を考慮して基本燃料噴射量Ti mを補正するためのものである。
【0023】
そして、これらの第1補正係数KTOTAL及び第2補正係数KCMDMによる基本燃料噴射量Timの補正は、第1補正係数KTOTAL及び第2補正係数KCMDMを基本燃料噴射量Timに乗算することで行われ、この補正によりエンジン1の要求燃料噴射量Tcylが得られる。
【0024】
尚、前記基本燃料噴射量Timや、第1補正係数KTOTAL、第2補正係数KCMDMのより具体的な算出手法は、特開平5−79374号公報等に本願出願人が開示しているので、ここでは詳細な説明を省略する。
【0025】
制御ユニット7は、上記の機能的構成の他、LAFセンサ5の出力KACTと所定の基準値FLAF/BASE(本実施形態ではこの基準値FLAF/BASEは空燃比換算で約「1」(一定値)とされている)との偏差kact(=KACT−FLAF/BASE)を求める減算処理部11と、O2センサ6の出力VO2/OUTとその目標値VO2/TARGET(本実施形態ではこの目標値VO2/TARGETは触媒装置3の最適な浄化性能が得られる所定の一定値とされている)との偏差VO2(=VO2/OUT−VO2/TARGET)を求める減算処理部12と、これらの偏差kact,VO2のデータをそれぞれLAFセンサ5の出力及びO2センサ6の出力を示すデータとして用い(以下、偏差kact,VO2をそれぞれLAFセンサ5の偏差出力kact及びO2センサ6の偏差出力VO2と称する)、LAFセンサ5の箇所の排ガスの目標空燃比KCMDを触媒装置3に進入するエンジン1の排ガスの空燃比を規定する操作量として決定する空燃比操作量決定部13と、この目標空燃比KCMDにLAFセンサ5の出力KACT(触媒装置3に進入する排ガスの検出空燃比)を一致(収束)させるようにエンジン1の燃料噴射量(燃料供給量)をフィードバック制御するフィードバック制御部14とを備えている。
【0026】
前記空燃比操作量決定部13は、その詳細は後述するが、排気管2のLAFセンサ5の箇所からO2センサ6の箇所にかけての触媒装置3を含む排気系(図1で参照符号Eを付した部分)を制御対象とし、その対象排気系E(プラント)に存する無駄時間や該対象排気系Eの挙動変化等を考慮しつつ、スライディングモード制御(より詳しくは適応スライディングモード制御)を用いてO2センサ6の出力VO2/OUTをその目標値VO2/TARGETに整定させるように(O2センサ6の偏差出力VO2を「0」に収束させるように)、LAFセンサ5の箇所の目標空燃比KCMDを逐次決定するものである。
【0027】
また、フィードバック制御部14は、本実施形態では、エンジン1の各気筒への全体的な燃料噴射量をフィードバック制御する大局的フィードバック制御部15と、エンジン1の各気筒毎の燃料噴射量をフィードバック制御する局所的フィードバック制御部16とにより構成されている。
【0028】
前記大局的フィードバック制御部15は、LAFセンサ5の出力KACTが前記目標空燃比KCMDに収束するように、前記要求燃料噴射量Tcylを補正する(要求燃料噴射量Tcylに乗算する)フィードバック補正係数KFBを逐次求めるものであり、LAFセンサ5の出力KACTと目標空燃比KCMDとの偏差に応じて周知のPID制御を用いて前記フィードバック補正係数KFBとしてのフィードバック操作量KLAFを生成するPID制御器17と、LAFセンサ5の出力KACTと目標空燃比KCMDとからエンジン1の運転状態の変化や特性変化等を考慮して前記フィードバック補正係数KFBを規定するフィードバック操作量KSTRを適応的に求める漸化式形式の制御器である適応制御器18(図ではSTRと称している)とをそれぞれ独立的に具備している。
【0029】
ここで、本実施形態では、前記PID制御器17が生成するフィードバック操作量KLAFは、LAFセンサ5の出力KACT(検出空燃比)が目標空燃比KCMDに一致している状態で「1」となり、該操作量KLAFをそのまま前記フィードバック補正係数KFBとして使用できるようになっている。一方、適応制御器18が生成するフィードバック操作量KSTRはLAFセンサ5の出力KACTが目標空燃比KCMDに一致する状態で「目標空燃比KCMD」となるもので、該フィードバック操作量KSTRを除算処理部19で目標空燃比KCMDにより除算してなるフィードバック操作量kstr(=KSTR/KCMD)が前記フィードバック補正係数KFBとして使用できるようになっている。
【0030】
そして、大局的フィードバック制御部15は、PID制御器17により生成されるフィードバック操作量KLAFと、適応制御器18が生成するフィードバック操作量KSTRを目標空燃比KCMDにより除算してなるフィードバック操作量kstrとを切換部20で適宜、択一的に選択して、いずれか一方のフィードバック操作量KLAF又はkstrを前記フィードバック補正係数KFBとして使用し、該補正係数KFBを前記要求燃料噴射量Tcylに乗算することにより該要求燃料噴射量Tcylを補正する。尚、かかる大局的フィードバック制御部15(特に適応制御器18)については後にさらに詳細に説明する。
【0031】
また、前記局所的フィードバック制御部16は、LAFセンサ5の出力KACTから各気筒毎の実空燃比#nA/F(n=1,2,3,4)を推定するオブザーバ21と、このオブザーバ21により推定された各気筒毎の実空燃比#nA/Fから各気筒毎の空燃比のばらつきを解消するよう、PID制御を用いて各気筒毎の燃料噴射量のフィードバック補正係数#nKLAFをそれぞれ求める複数(気筒数個)のPID制御器22とを具備する。
【0032】
ここで、オブザーバ21は、それを簡単に説明すると、各気筒毎の実空燃比#nA/Fの推定を次のように行うものである。すなわち、エンジン1からLAFセンサ5の箇所(各気筒毎の排ガスの集合部)にかけてのシステムを、各気筒毎の実空燃比#nA/FからLAFセンサ5で検出される排ガスの空燃比を生成するシステムと考え、これを、LAFセンサ5の検出応答遅れ(例えば一次遅れ)や、各気筒毎の排ガスの集合部における空燃比に対する各気筒毎の空燃比の時間的寄与度を考慮してモデル化する。そして、そのモデルの基で、LAFセンサ5の出力KACT(検出空燃比)から、逆算的に各気筒毎の実空燃比#nA/Fを推定する。
尚、このようなオブザーバ21は、本願出願人が例えば特開平7−83094号公報に詳細に開示しているので、ここでは、さらなる説明を省略する。
【0033】
また、局所的フィードバック制御部16の各PID制御器22は、LAFセンサ5の出力KACTを、前回の制御サイクルで各PID制御器22により求められたフィードバック補正係数#nKLAFの全気筒についての平均値により除算してなる値を各気筒の空燃比の目標値として、その目標値とオブザーバ21により求められた各気筒毎の実空燃比#nA/Fとの偏差が解消するように、今回の制御サイクルにおける、各気筒毎のフィードバック補正係数#nKLAFを求める。
そして、局所的フィードバック制御部16は、前記要求燃料噴射量Tcylに大局的フィードバック制御部15のフィードバック補正係数KFBを乗算してなる値に、各気筒毎のフィードバック補正係数#nKLAFを乗算することで、各気筒の出力燃料噴射量#nTout(n=1,2,3,4)を求める。
このようにして求められる各気筒の出力燃料噴射量#nToutは、制御ユニット7に備えた各気筒毎の付着補正部23により吸気管の壁面付着を考慮した補正が各気筒毎になされた後、エンジン1の図示しない燃料噴射装置に与えられ、その付着補正がなされた出力燃料噴射量#nToutで、エンジン1の各気筒への燃料噴射が行われるようになっている。
【0034】
尚、上記付着補正については、本願出願人が例えば特開平8−21273号公報に詳細に開示しているので、ここではさらなる説明を省略する。また、図1において、参照符号24を付したセンサ出力選択処理部は、前記オブザーバ21による各気筒毎の実空燃比#nA/Fの推定に適したLAFセンサ5の出力KACTをエンジン1の運転状態に応じて選択するもので、これについては、本願出願人が特開平7−259488号公報にて詳細に開示しているので、ここではさらなる説明を省略する。
【0035】
次に、前記空燃比操作量決定部13を詳細に説明する。
【0036】
前述の如く、空燃比操作量決定部13は、前記対象排気系Eに存する無駄時間や該排気系Aの挙動変化等を考慮しつつ、適応スライディングモード制御を用いてO2 センサ6の出力VO2/OUTをその目標値VO2/TARGETに整定させるようにLAFセンサ5の箇所の排ガスの目標空燃比KCMDを逐次決定するものである。そして、このような制御処理を行うために、本実施形態では、あらかじめ前記対象排気系Eを、前記LAFセンサ5の出力KACT(触媒装置3に進入する排ガスの空燃比)から無駄時間要素及び応答遅れ要素を介してO2センサ6の出力VO2/OUT(触媒装置3を通過した排ガス中の酸素濃度)を生成するプラントと見なし、それを離散系でモデル化している。
【0037】
この場合、本実施形態では、空燃比操作量決定部13による処理の簡素化を図るために、LAFセンサ5の出力KACT及びO2センサ6の出力VO2/OUTの代わりに、LAFセンサ5の前記偏差出力kact(=KACT−FLAF/BASE)とO2センサ6の前記偏差出力VO2(=VO2/OUT−VO2/TARGET)とを用いて、対象排気系Eの離散系モデルを次式(1)により表す。
【0038】
【数1】
【0039】
この式(1)は対象排気系EがLAFセンサ5の偏差出力kactから、無駄時間要素及び応答遅れ要素を介してO2センサ6の偏差出力VO2を生成するプラントであるとみなして、該対象排気系Eを離散系でモデル化してなるもの(無駄時間を有する自己回帰モデル)であり、上式(1)において、「k」は離散時間的な制御サイクルの番数を示し、「d」は対象排気系Eの無駄時間を制御サイクル数で表したものである。この場合、本実施形態では、対象排気系Eの無駄時間は、例えば制御サイクルの周期を30〜100msとして、d制御サイクル分の時間(d=3〜10)とされている。また、上式(1)の右辺第1項及び第2項はそれぞれ対象排気系Eの応答遅れ要素に対応するもので、第1項は1次目の自己回帰項、第2項は2次目の自己回帰項である。そして、「a1」、「a2」はそれぞれ1次目の自己回帰項、2次目の自己回帰項のゲイン係数である。さらに、上式(1)の右辺第3項は対象排気系Eの無駄時間要素に対応するもので、「b1」はその無駄時間要素に係わるゲイン係数である。これらのゲイン係数a1,a2,b1は離散系モデルを規定するパラメータである。
【0040】
本実施形態における前記空燃比操作量決定部13は、式(1)により表される離散系モデルに基づき、所定(一定)の制御サイクルで前述のような制御処理を行うもので、その機能的構成は、図3に示すように大別される。
【0041】
すなわち、空燃比操作量決定部13は、LAFセンサ5の偏差出力kact及びO2センサ6の偏差出力VO2のデータから、前記離散系モデルの設定すべきパラメータである前記ゲイン係数a1,a2,b1の値を制御サイクル毎に逐次同定する同定器25と、LAFセンサ5の偏差出力kact及びO2センサ6の偏差出力VO2のデータから、前記同定器25により同定された前記ゲイン係数a1,a2,b1の同定値a1ハット,a2ハット,b1ハット(以下、同定ゲイン係数a1ハット,a2ハット,b1ハットという)を用いて対象排気系Eの無駄時間d後のO2センサ6の偏差出力VO2の推定値VO2バー(以下、推定偏差出力VO2バーという)を制御サイクル毎に逐次求める推定器26と、該推定器26により求められたO2センサ6の推定偏差出力VO2バーから、前記同定ゲイン係数a1ハット,a2ハット,b1ハットを用いて適応スライディングモード制御によりLAFセンサ5の箇所の排ガス(触媒装置3に進入する排ガス)の目標空燃比KCMDを制御サイクル毎に逐次決定するスライディングモード制御器27とにより構成されている。
【0042】
これらの同定器25、推定器26及びスライディングモード制御器27による演算処理のアルゴリズムは前記離散系モデルに基づいて以下のように構築されている。
【0043】
まず、前記同定器25に関し、前記離散系モデルのゲイン係数a1,a2,b1に対応する実際の対象排気系Eのゲイン係数は一般に該対象排気系Eの挙動状態や経時的な特性変化等によって変化する。従って、前記離散系モデルの実際の対象排気系Eに対するモデル化誤差を極力少なくして該離散系モデルの精度を高めるためには、離散系モデルのゲイン係数a1,a2,b1を実際の対象排気系Eの挙動状態等に則して適宜、リアルタイムで同定することが好ましい。
【0044】
前記同定器25は、上記のように離散系モデルのモデル化誤差を極力小さくするために、前記ゲイン係数a1,a2,b1をリアルタイムで逐次同定するものであり、その同定処理は次のように行われる。
【0045】
すなわち、同定器25は、所定の制御サイクル毎に、まず、今現在設定されている離散系モデルの同定ゲイン係数a1ハット,a2ハット,b1ハット、すなわち前回の制御サイクルで決定した同定ゲイン係数a1ハット(k-1),a2ハット(k-1),b1ハット(k-1)と、LAFセンサ5の偏差出力kact及びO2センサ6の偏差出力VO2の過去に得られたデータとを用いて、次式(2)により今現在設定されている離散系モデル上でのO2センサ6の今現在の偏差出力VO2の同定値VO2ハット(以下、同定偏差出力VO2ハットという)を求める。
【0046】
【数2】
【0047】
尚、この式(2)は、式(1)を1制御サイクル分、過去側にシフトし、ゲイン係数a1,a2,b1を同定ゲイン係数a1ハット(k-1),a2ハット(k-1),b1ハット(k-1)で置き換えたものである。また、式(2)の第3項で用いる「d」は、対象排気系Eの無駄時間の設定値(より詳しくは無駄時間の設定値を制御サイクル数で表したもの)であり、その設定値は対象排気系Eの実際の無駄時間と等しいか、もしくはそれよりも若干長い時間になるように設定されている。
【0048】
ここで、次式(3),(4)で定義されるベクトルΘ及びξを導入すると(式(3),(4)中の添え字「T」は転置を意味する。以下同様。)、
【0049】
【数3】
【0050】
【数4】
【0051】
前記式(2)は、次式(5)により表される。
【0052】
【数5】
【0053】
さらに同定器25は、前記式(2)あるいは式(5)により求められるO2センサ6の同定偏差出力VO2ハットと今現在のO2センサ6の偏差出力VO2との偏差id/eを離散系モデルの実際の対象排気系Eに対するモデル化誤差を表すものとして次式(6)により求める(以下、偏差id/eを同定誤差id/eという)。
【0054】
【数6】
【0055】
そして、同定器25は、上記同定誤差id/eを最小にするように新たな同定ゲイン係数a1(k)ハット,a2(k)ハット,b1(k)ハット、換言すれば、これらの同定ゲイン係数を要素とする新たな前記ベクトルΘ(k)(以下、このベクトルを同定ゲイン係数ベクトルΘという)を求めるもので、その算出を、次式(7)により行う。すなわち、同定器25は、前回の制御サイクルで決定した同定ゲイン係数a1ハット(k-1),a2ハット(k-1),b1ハット(k-1)を、同定誤差id/eに比例させた量だけ変化させることで新たな同定ゲイン係数a1(k)ハット,a2(k)ハット,b1(k)ハットを求める。
【0056】
【数7】
【0057】
ここで、式(7)中の「Kθ」は次式(8)により決定される三次のベクトル(各同定ゲイン係数a1ハット,a2ハット,b1ハットの同定誤差id/eに応じた変化度合いを規定するゲイン係数ベクトル)である。
【0058】
【数8】
【0059】
また、上式(8)中の「P」は次式(9)の漸化式により決定される三次の正方行列である。
【0060】
【数9】
【0061】
尚、式(9)中の「λ1」、「λ2」は0<λ1≦1及び0≦λ2<2の条件を満たすように設定され、また、「P」の初期値P(0)は、その各対角成分を正の数とする対角行列である。
【0062】
この場合、式(9)中の「λ1」、「λ2」の設定の仕方によって、固定ゲイン法、漸減ゲイン法、重み付き最小二乗法、最小二乗法、固定トレース法等、各種の具体的なアルゴリズムが構成され、本実施形態では、例えば最小二乗法(この場合、λ1=λ2=1)を採用している。
【0063】
本実施形態における同定器25は基本的には前述のようなアルゴリズム(演算処理)によって、前記同定誤差id/eを最小化するように離散系モデルの前記同定ゲイン係数a1ハット,a2ハット,b1ハットを制御サイクル毎に逐次求めるもので、このような処理によって、実際の対象排気系Eに適合した同定ゲイン係数a1ハット,a2ハット,b1ハットが逐次得られる。
【0064】
尚、本実施形態における同定器25は、前記同定誤差id/eの算出に際して、O2センサ6の前記同定偏差出力VO2ハットとO2センサ6の偏差出力VO2とにフィルタリング処理を施したり、ゲイン係数a1,a2,b1の同定(同定ゲイン係数a1ハット,a2ハット,b1ハットの更新)を対象排気系Eの特定の挙動状態において行ったりするのであるが、これについては後述する。
【0065】
次に、前記推定器26は、後に詳細を説明するスライディングモード制御器27による目標空燃比KCMDの決定処理に際しての対象排気系Eの無駄時間dの影響を補償するために、該無駄時間d後のO2センサ6の偏差出力VO2の推定値である前記推定偏差出力VO2バーを制御サイクル毎に逐次求めるものであり、その推定処理は次のように行われる。
【0066】
まず、前記式(1)で表される離散系モデルにおいて、次式(10)により定義されるベクトルXを導入すると、
【0067】
【数10】
【0068】
式(1)は次式(11)に書き換えられる。
【0069】
【数11】
【0070】
ここで、式(11)の漸化式を繰り返し用いると、無駄時間d後のX(k+d)は、式(11)中で定義した行列A及びベクトルBやLAFセンサ5の偏差出力kactの時系列データkact(k-j)(j=1,2, d)を用いて、次式(12)により表される。
【0071】
【数12】
【0072】
この場合、式(12)の左辺の第1行成分が無駄時間d後のO2センサ6の偏差出力VO2(k+d)であるから、その推定値(推定偏差出力)VO2(k+d)バーは、式(12)の右辺の第1行成分を演算することで求めることができる。
【0073】
そこで、式(12)の両辺の第1行成分に着目し、右辺第1項の行列Adの第1行第1列成分及び第1行第2列成分をそれぞれα1,α2とおき、右辺第2項のベクトルAj-1・B(j=1,2, d)の第1行成分をそれぞれβj(j=1,2, d)とおくと、O2センサ6の推定偏差出力VO2(k+d)バーは、O2センサ6の偏差出力VO2の時系列データVO2(k)及びVO2(k-1)と、LAFセンサ5の偏差出力kactの時系列データkact(k-j)(j=1,2, d)とを用いて次式(13)により求めることができる。
【0074】
【数13】
【0075】
また、式(13)中の係数値α1,α2及びβj(j=1,2, d)は、行列A及びベクトルBの成分(式(11)参照)を構成するゲイン係数a1,a2,b1として、前記同定器25により求められた同定ゲイン係数a1ハット,a2ハット,b1ハットを用い、それらの行列A及びベクトルBから前記式(12)中の行列Ad及びベクトルAj-1・B(j=1,2, d)を求めることで決定することができる。
【0076】
よって、本実施形態における推定器26は、基本的には同定器25により求められる前記同定ゲイン係数a1ハット,a2ハット,b2ハット(詳しくは現在の制御サイクルで求められた同定ゲイン係数a1(k)ハット,a2(k)ハット,b2(k)ハット)を用いて、式(13)中の係数値α1,α2及びβj(j=1,2, d)を算出する。さらに、その算出した係数値α1,α2及びβjと、O2センサ6の偏差出力VO2の現在以前の時系列データVO2(k)及びVO2(k-1)と、LAFセンサ5の偏差出力kactの過去の時系列データkact(k-j)(j=1,2, d)とから式(13)の演算を行うことで、O2センサ6の推定偏差出力VO2(k+d)バーを求める。これが推定器26における基本的な演算処理(推定アルゴリズム)である。
【0077】
尚、本実施形態における推定器26は、O2センサ6の推定偏差出力VO2(k+d)バーの算出に際して、LAFセンサ5の偏差出力kactにフィルタリング処理を施すのであるがこれについては後述する。
【0078】
次に、前記スライディングモード制御器27を詳細に説明する。
【0079】
ここで、まず、一般的なスライディングモード制御について図4を参照して簡単に説明しておく。
【0080】
スライディングモード制御は、可変構造型のフィードバック制御手法であり、この制御手法においては、例えば制御対象の制御すべき状態量をx1,x2の二つとした場合、これらの状態量x1,x2を変数とする線形関数σ=s1x1+s2x2(s1,s2は係数)を用いて、σ=0により表される超平面をあらかじめ設計しておく。この超平面σ=0は位相空間が二次系の場合(状態量が二つの場合)は、しばしば切換線と呼ばれ、線形関数σは切換関数と呼ばれている。位相空間の次数がさらに大きくなると、切換線から切換面となり、さらには幾何学的に図示できなくなる超平面になる。尚、超平面はすべり面と呼ばれることもある。本明細書においては、これらを代表して線形関数及び超平面と表現した。
【0081】
そして、このスライディングモード制御は、例えば図4の点Pで示すように、状態量x1,x2がσ≠0となっている場合に、所謂、到達則に従って、状態量x1,x2をハイゲイン制御によって超平面σ=0上に高速で収束させ(モード1)、さらに所謂、等価制御入力によって状態量x1,x2を超平面σ=0上に拘束しつつ超平面σ=0上の平衡点(x1=x2=0の点)に収束させる(モード2)ものである。
【0082】
このようなスライディングモード制御においては、状態量x1,x2を超平面σ=0上に収束させさえすれば、等価制御入力によって、外乱等の影響を受けることなく、極めて安定に状態量x1,x2を超平面σ=0上に拘束して、該超平面σ=0の平衡点に収束させることができるという特性をもっている。尚、外乱や制御対象のモデル化誤差があると、状態量x1,x2は厳密には上記平衡点(x1=x2=0の点)には収束せず、該平衡点の近傍に収束する。
【0083】
かかるスライディングモード制御では、特に、上記モード1において状態量x1,x2をいかにして安定に超平面σ=0上に収束させるかが重要な課題となる。この場合、外乱等の影響があると、一般には、前記到達則だけでは、状態量x1,x2を超平面σ=0上に安定に収束させることが困難である。このため、近年では、例えばコロナ社により1994年10月20日に発刊された「スライディングモード制御 −非線形ロバスト制御の設計理論−」と題する文献の第134頁〜第135頁に見られるように、到達則に加えて、外乱の影響を排除しつつ状態量を超平面上に収束させるための適応則を用いた適応スライディングモード制御という手法が提案されている。
【0084】
本実施形態の前記スライディングモード制御器27は、このような適応スライディングモード制御を用いて、O2センサ6の出力VO2/OUTをその目標値VO2/TARGETに整定させるように(O2センサ6の偏差出力VO2を「0」に収束させるように)、制御対象である前記対象排気系Eに与えるべき入力(詳しくは、LAFセンサ5で検出される排ガスの空燃比と前記基準値FLAF/BASEとの偏差で、これはLAFセンサ5の偏差出力kactに相当する。以下、この入力をSLD操作入力uslと称する)を決定し、その決定したSLD操作入力uslから前記目標空燃比KCMDを決定するものである。そして、その処理のためのアルゴリズムは次のように構築されている。
【0085】
まず、スライディングモード制御器27の適応スライディングモード制御に必要な超平面の構築について説明する。
【0086】
スライディングモード制御器27は、O2センサ6の偏差出力VO2を「0」に収束させるように制御を行うものであるので、O2センサ6の偏差出力VO2の時系列データを「0」に収束させるように対象排気系Eに与えるべき前記SLD操作入力uslを決定すればよい。
【0087】
そこで、本実施形態におけるスライディングモード制御の基本的な考え方としては、制御すべき状態量として、例えば各制御サイクルで得られたO2センサ6の偏差出力VO2(k)と、その1制御サイクル前に得られた偏差出力VO2(k-1)とを用い、スライディングモード制御用の超平面を規定する線形関数σを次式(14)により設定する。
【0088】
【数14】
【0089】
尚、本実施形態では、線形関数の変数である状態量として、実際には前記推定器26により求められる前記推定偏差出力VO2バーの時系列データを用いるのであるがこれについては後述する。
【0090】
上記のように線形関数σを定義したとき、スライディングモード制御用の超平面はσ=0により表され(この場合、状態量は二つであるので超平面は直線となる。図4参照)、この超平面σ=0を規定する線形関数σの係数s1,s2(式(14)参照)は、本実施形態ではあらかじめ次のように設定する。
【0091】
すなわち、前記状態量VO2(k),VO2(k-1)を成分とする式(14)中のベクトルX(以下、単に状態量Xという)が超平面σ=0上に収束した状態では、線形関数σの値が「0」であるので、これと式(14)とから次式(15)が得られる。
【0092】
【数15】
【0093】
ここで、式(15)により表される系は、入力の無い一次遅れ系であるので、状態量Xが超平面σ=0の平衡点(VO2(k)=VO2(k-1)=0となる点)に安定に収束するための条件は、式(15)により表される系の極(この場合、この極は「−s2/s1」である)が単位円内に存在することとなる。
【0094】
従って、本実施形態では、線形関数σの係数s1,s2は、次式(16)の条件を満たすように設定する。
【0095】
【数16】
【0096】
尚、本実施形態では、簡略化のために係数s1=1とし(この場合、s2/s1=s2である)、−1<s2<1の条件を満たすように係数s2の値を設定する。
【0097】
一方、上記のように設定された超平面σ=0の平衡点に前記状態量Xを収束させるためにスライディングモード制御器27が適応スライディングモード制御により生成すべき前記SLD操作入力usl(LAFセンサ5で検出される排ガスの空燃比と前記基準値FLAF/BASE との偏差)は、前記状態量Xを超平面σ=0上に拘束するための制御則に従って対象排気系Eに与えるべき等価制御入力ueqと、状態量Xを超平面σ=0に収束させるための到達則に従って対象排気系Eに与えるべき入力urch(以下、到達則入力urchという)と、外乱等の影響を補償して状態量Xを超平面σ=0に収束させるための適応則に従って対象排気系Eに与えるべき入力uadp(以下、適応則入力uadpという)との総和により表される(次式(17)参照)。
【0098】
【数17】
【0099】
そして、これらの等価制御入力ueq、到達則入力urch及び適応則入力uadpは、本実施形態では、前記式(1)あるいは式(11)により表される離散系モデルに基づいて、次のように求めることができる。
【0100】
まず、等価制御入力ueqに関し、前記状態量Xが超平面σ=0上に留まる条件は、σ(k+1)=σ(k)=0であり、この条件は、前記式(11)及び式(14)を用いて、次式(18)に書き換えられる。
【0101】
【数18】
【0102】
ここで、等価制御入力ueqは、状態量Xを超平面σ=0に拘束するために対象排気系Eに与えるべき入力(LAFセンサ5で検出される排ガスの空燃比と前記基準値FLAF/BASEとの偏差)であるので、上記式(18)の条件を満たすLAFセンサ5の偏差出力kactが等価制御入力ueqである。
【0103】
従って、式(18)から次式(19)が得られ、
【0104】
【数19】
【0105】
さらに、この式(19)の両辺を無駄時間d分シフトすることで、次式(20)が得られる。
【0106】
【数20】
【0107】
この式(20)が本実施形態において、制御サイクル毎に等価制御入力ueqを求めるための基本式である。
【0108】
次に、前記到達則入力urchは、本実施形態では、基本的には次式(21)により決定するものとする。
【0109】
【数21】
【0110】
すなわち、到達則入力urchは、無駄時間dの影響を考慮し、無駄時間d後の線形関数σの値σ(k+d)に比例させるように決定する。
【0111】
この場合、式(21)中の係数F(これは到達則のゲインを規定する)は、次にように設定される。
【0112】
すなわち、前記式(11)において、kact(k)=ueq(k)+urch(k)とし、さらに式(14)、(20)、(21)を用いると、次式(22)が得られる。
【0113】
【数22】
【0114】
ここで、式(22)により表される系は、入力の無い一次遅れ系であるので、線形関数σの値が超平面σ=0に安定に収束する(状態量Xが超平面σ=0に安定に収束する)ための条件は、式(22)により表される系の極(この場合、この極は「1−F」である)が単位円内に存在することとなる。
【0115】
従って、本実施形態では、到達則入力urchを規定する係数Fは、次式(23)の条件を満たすように設定する。
【0116】
【数23】
【0117】
尚、線形関数σの値の挙動に関しては、該線形関数σの値が超平面σ=0に対して振動的な変化(所謂チャタリング)を生じる虞れがあり、このチャタリングを抑制するためには、式(22)により表される系の極「1−F」が上記の条件に加えて0<1−F<1なる条件を満たすことが好ましい。
【0118】
従って、到達則入力urchを規定する係数Fは、より好ましくは、次式(24)の条件を満たすように設定する。
【0119】
【数24】
【0120】
次に、前記適応則入力uad pは、本実施形態では、基本的には次式(25)により決定するものとする(式(25)中のΔTは制御サイクルの周期である)。
【0121】
【数25】
【0122】
すなわち、適応則入力uadpは、無駄時間dの影響を考慮し、無駄時間d後までの線形関数σの値の制御サイクル毎の積算値(これは線形関数σの値の積分値に相当する)に比例させるように決定する。
【0123】
この場合、式(25)中の係数G(これは適応則のゲインを規定する)は、次にように設定される。
【0124】
まず、LAFセンサ5で検出される空燃比と目標空燃比KCMDとの間の外乱等の影響による誤差成分をvとすると、LAFセンサ5の偏差出力kactは、前記等価制御入力ueq、到達則入力urch及び適応則入力uadp、並びに上記誤差成分vを用いて、次式(26)により表現することができる。
【0125】
【数26】
【0126】
そして、この式(26)を前記式(11)に適用し、さらに式(14)、(20)、(21)、(25)を用いると、次式(27)が得られる。
【0127】
【数27】
【0128】
ここで、式(27)の両辺をZ変換すると次式(28)となり、
【0129】
【数28】
【0130】
さらにこの式(28)を変形して整理すると、次式(29)が得られる。
【0131】
【数29】
【0132】
尚、式(28)、(29)における「Σ」及び「V」はそれぞれ線形関数σ及び前記誤差成分vをZ変換したものである。また、式(29)のM(Z)は誤差成分vに対する線形関数σのパルス伝達関数で、式(29)の上段の分数式により表されるものである。
【0133】
この場合、線形関数σが誤差成分v(外乱)に対して、安定となる条件は、前記パルス伝達関数M(Z)の極、すなわち、次式(30)により示される特性方程式の解(この解は二つある)が単位円内に存在することであり、
【0134】
【数30】
【0135】
上記特性方程式(30)の二つの解はそれをλm1、λm2とすると、次式(31)、(32)により与えられる。
【0136】
【数31】
【0137】
【数32】
【0138】
従って、線形関数σが誤差成分v(外乱)に対して安定となる条件は、上式(31)、(32)により与えられるλm1、λm2が単位円内に存在することである。
【0139】
そこで、本実施形態ではこの条件を満たすために、前記係数Gは、次式(33)により設定する。
【0140】
【数33】
【0141】
本実施形態におけるスライディングモード制御器27は、基本的には前記式(20)、(21)、(25)により決定される等価制御入力ueq、到達則入力urch及び適応則入力uadpの総和(ueq+urc h+uadp)を対象排気系Eに与えるべきSLD操作入力uslとして決定するのであるが、前記式(20)、(21)、(25)で使用するO2センサ6の偏差出力VO2(k+d),VO2(k+d-1)や、線形関数σの値σ(k+d) 等は未来値であるので実際には得られないものである。
【0142】
そこで、本実施形態では、スライディングモード制御器27は、実際には、前記式(20)により前記等価制御入力ueqを決定するためのO2センサ6の偏差出力VO2(k+d),VO2(k+d-1)の代わりに、前記推定器26で求められる推定偏差出力VO2(k+d)バー,VO2(k+d-1)バーを用い、次式(34)により制御サイクル毎の等価制御入力ueqを算出する。
【0143】
【数34】
【0144】
また、本実施形態では、実際には、推定器26により前述の如く逐次求められた推定偏差出力VO2バーの時系列データを制御すべき状態量とし、前記式(14)により設定された線形関数σに代えて、次式(35)により線形関数σバーを定義する(この線形関数σバーは、前記式(14)の偏差出力VO2の時系列データを推定偏差出力VO2バーの時系列データで置き換えたものに相当する)。
【0145】
【数35】
【0146】
そして、スライディングモード制御器27は、前記式(21)により前記到達則入力urchを決定するための線形関数σの値の代わりに、前記式(35)により表される線形関数σバーの値を用いて次式(36)により制御サイクル毎の到達則入力urchを算出する。
【0147】
【数36】
【0148】
同様に、スライディングモード制御器27は、前記式(25)により前記適応則入力uadpを決定するための線形関数σの値の代わりに、前記式(35)により表される線形関数σバーの値を用いて次式(37)により制御サイクル毎の適応則入力uadpを算出する。
【0149】
【数37】
【0150】
尚、前記式(34),(36),(37)により等価制御入力ueq、到達則入力urch及び適応則入力uadpを算出する際に必要となる前記ゲイン係数a1,a2,b1としては、本実施形態では基本的には前記同定器25により求められた最新の同定ゲイン係数a1(k)ハット,a2(k)ハット,b1(k)ハットを用いる。
【0151】
そして、スライディングモード制御器27は、前記式(34)、(36)、(37)によりそれぞれ求められる等価制御入力ueq、到達則入力urch及び適応則入力uadpの総和を対象排気系Eに与えるべき前記SLD操作入力uslとして求める(前記式(17)を参照)。尚、この場合において、前記式(34)、(36)、(37)中で用いる前記係数s1,s2,F,Gの設定条件は前述の通りである。
【0152】
これが、本実施形態において、スライディングモード制御器27により、対象排気系Eに与えるべきSLD操作入力uslを制御サイクル毎に決定するための基本的なアルゴリズムである。このようにしてSLD操作入力uslを決定することで、該SLD操作入力uslは、O2センサ6の推定偏差出力VO2バーを「0」に収束させるように(結果的にはO2センサ6の出力VO2バーを目標値VO2/TARGETに収束させるように)決定される。
【0153】
ところで、本実施形態におけるスライディングモード制御器27は最終的には前記目標空燃比KCMDを制御サイクル毎に逐次求めるものあるが、前述のように求められるSLD操作入力uslは、LAFセンサ5で検出される排ガスの空燃比と前記基準値FLAF/BASEとの偏差の目標値を意味する。このため、スライディングモード制御器27は、最終的には、次式(38)に示すように、制御サイクル毎に、前述の如く求めたSLD操作入力uslに前記基準値FL AF/BASEを加算することで、目標空燃比KCMDを決定する。
【0154】
【数38】
【0155】
以上が本実施形態でスライディングモード制御器27により目標空燃比KCMDを決定するための基本的アルゴリズムである。
【0156】
尚、本実施形態では、スライディングモード制御器27による適応スライディングモード制御の処理の安定性を判別して、前記SLD操作入力uslの値を制限したりするのであるが、これについては後述する。
【0157】
次に、前記大局的フィードバック制御部15、特に前記適応制御器18をさらに説明する。
【0158】
前記図1を参照して、大局的フィードバック制御部15は、前述のようにLAFセンサ5の出力KACT(検出空燃比)を目標空燃比KCMDに収束させるようにフィードバック制御を行うものであるが、このとき、このようなフィードバック制御を周知のPID制御だけで行うようにすると、エンジン1の運転状態の変化や経年的特性変化等、動的な挙動変化に対して、安定した制御性を確保することが困難である。
【0159】
前記適応制御器18は、上記のようなエンジン1の動的な挙動変化を補償したフィードバック制御を可能とするもので、I.D.ランダウ等により提唱されているパラメータ調整則を用いて、図5に示すように、複数の適応パラメータを設定するパラメータ調整部28と、設定された適応パラメータを用いて前記フィードバック操作量KSTRを算出する操作量算出部29とにより構成されている。
【0160】
ここで、パラメータ調整部28について説明すると、ランダウ等の調整則では、離散系の制御対象の伝達関数B(Z-1)/A(Z-1)の分母分子の多項式を一般的に下記の式(39),(40)のようにおいたとき、パラメータ調整部28が設定する適応パラメータθハット(j)(jは制御サイクルの番数を示す)は、式(41)のようにベクトル(転置ベクトル)で表される。また、パラメータ調整部28への入力ζ(j)は、式(42)のように表される。この場合、本実施形態では、大局的フィードバック制御部15の制御対象であるエンジン1が一次系で3制御サイクル分の無駄時間dp(エンジン1の燃焼サイクルの3サイクル分の時間)を持つプラントと考え、式(39)〜式(42)でm=n=1,dp=3とし、設定する適応パラメータはs0,r1,r2,r3,b0の5個とした(図5参照)。尚、式(42)の上段式及び中段式におけるus,ysは、それぞれ、制御対象への入力(操作量)及び制御対象の出力(制御量)を一般的に表したものであるが、本実施形態では、上記入力をフィードバック操作量KSTR、制御対象(エンジン1)の出力を前記LAFセンサ5の出力KACT(検出空燃比)とし、パラメータ調整部28への入力ζ(j)を、式(42)の下段式により表す(図5参照)。
【0161】
【数39】
【0162】
【数40】
【0163】
【数41】
【0164】
【数42】
【0165】
ここで、前記式(41)に示される適応パラメータθハットは、適応制御器18のゲインを決定するスカラ量要素b0ハット-1(Z-1,j)、操作量を用いて表現される制御要素BRハット(Z-1,j)、及び制御量を用いて表現される制御要素S(Z-1,j)からなり、それぞれ、次式(43)〜(45)により表現される(図5の操作量算出部29のブロック図を参照)。
【0166】
【数43】
【0167】
【数44】
【0168】
【数45】
【0169】
パラメータ調整部28は、これらのスカラ量要素や制御要素の各係数を設定して、それを式(26)に示す適応パラメータθハットとして操作量算出部29に与えるもので、現在から過去に渡るフィードバック操作量KSTRの時系列データとLAFセンサ5の出力KACTとを用いて、該出力KACTが前記目標空燃比KCMDに一致するように、適応パラメータθハットを算出する。
【0170】
この場合、具体的には、適応パラメータθハットは、次式(46)により算出する。
【0171】
【数46】
【0172】
同式(46)において、Γ(j)は、適応パラメータθハットの設定速度を決定するゲイン行列(この行列の次数はm+n+dp)、eアスタリスク(j)は、適応パラメータθハットの推定誤差を示すもので、それぞれ式(47),(48)のような漸化式で表される。
【0173】
【数47】
【0174】
【数48】
【0175】
ここで、式(48)中の「D(Z-1)」は、収束性を調整するための、漸近安定な多項式であり、本実施形態ではD(Z-1)=1としている。
【0176】
尚、式(47)のλ1(j),λ2(j)の選び方により、漸減ゲインアルゴリズム、可変ゲインアルゴリズム、固定トレースアルゴリズム、固定ゲインアルゴリズム等の種々の具体的なアルゴリズムが得られる。エンジン1の燃料噴射あるいは空燃比等の時変プラントでは、漸減ゲインアルゴリズム、可変ゲインアルゴリズム、固定ゲインアルゴリズム、および固定トレースアルゴリズムのいずれもが適している。
【0177】
前述のようにパラメータ調整部28により設定される適応パラメータθハット(s0,r1,r2,r3,b0)と、前記空燃比操作量決定部13により決定される目標空燃比KCMDとを用いて、操作量算出部29は、次式(49)の漸化式により、フィードバック操作量KSTRを求める。図5の操作量算出部29は、同式(49)の演算をブロック図で表したものである。
【0178】
【数49】
【0179】
尚、式(49)にり求められるフィードバック操作量KSTRは、LAFセンサ5の出力KACTが目標空燃比KCMDに一致する状態において、「目標空燃比KCMD」となる。このために、前述の如く、フィードバック操作量KSTRを除算処理部19によって目標空燃比KCMDで除算することで、前記フィードバック補正係数KFBとして使用できるフィードバック操作量kstrを求めるようにしている。
【0180】
このように構築された適応制御器18は、前述したことから明らかなように、制御対象であるエンジン1の動的な挙動変化を考慮した漸化式形式の制御器であり、換言すれば、エンジン1の動的な挙動変化を補償するために、漸化式形式で記述された制御器である。そして、より詳しくは、漸化式形式の適応パラメータ調整機構を備えた制御器と定義することができる。
【0181】
尚、この種の漸化式形式の制御器は、所謂、最適レギュレータを用いて構築する場合もあるが、この場合には、一般にはパラメータ調整機構は備えられておらず、エンジン1の動的な挙動変化を補償する上では、前述のように構成された適応制御器18が好適である。
【0182】
以上が、本実施形態で採用した適応制御器18の詳細である。
【0183】
尚、適応制御器18と共に、大局的フィードバック制御部15に具備したPID制御器17は、一般のPID制御と同様に、LAFセンサ5の出力KACT(検出空燃比)と、その目標空燃比KCMDとの偏差から、比例項(P項)、積分項(I項)及び微分項(D項)を算出し、それらの各項の総和をフィードバック操作量KLAFとして算出する。この場合、本実施形態では、積分項(I項)の初期値を“1”とすることで、LAFセンサ5の出力KACTが目標空燃比KCMDに一致する状態において、フィードバック操作量KLAFが“1”になるようにし、該フィードバック操作量KLAFをそのまま燃料噴射量を補正するための前記フィードバック補正係数KFBとして使用することができるようしている。また、比例項、積分項及び微分項のゲインは、エンジン1の回転数と吸気圧とから、あらかじめ定められたマップを用いて決定される。
【0184】
また、大局的フィードバック制御部15の前記切換部20は、エンジン1の冷却水温の低温時や、高速回転運転時、吸気圧の低圧時等、エンジン1の燃焼が不安定なものとなりやすい場合、あるいは、目標空燃比KCMDの変化が大きい時や、空燃比のフィードバック制御の開始直後等、これに応じたLAFセンサ6の出力KACTが、そのLAFセンサ5の応答遅れ等によって、信頼性に欠ける場合、あるいは、エンジン1のアイドル運転時のようエンジン1の運転状態が極めて安定していて、適応制御器18による高ゲイン制御を必要としない場合には、PID制御器17により求められるフィードバック操作量KLAFを燃料噴射量を補正するためのフィードバック補正量数KFBとして出力する。そして、上記のような場合以外の状態で、適応制御器18により求められるフィードバック操作量KSTRを目標空燃比KCMDで除算してなるフィードバック操作量kstrを燃料噴射量を補正するためのフィードバック補正係数KFBとして出力する。これは、適応制御器18が、高ゲイン制御で、LAFセンサ5の出力KACT(検出空燃比)を急速に目標空燃比KCMDに収束させるように機能するため、上記のようにエンジン1の燃焼が不安定となったり、LAFセンサ5の出力KACTの信頼性に欠ける等の場合に、適応制御器18のフィードバック操作量KSTRを用いると、かえって空燃比の制御が不安定なものとなる虞れがあるからである。
【0185】
このような切換部20の作動は、例えば特開平8−105345号公報に本願出願人が詳細に開示しているので、ここでは、さらなる説明を省略する。
【0186】
次に本実施形態の装置の作動の詳細を説明する。
【0187】
ここで、まず、制御ユニット7が行う処理の制御サイクルについて説明しておく。前記エンジン1の燃料供給量(燃料噴射量)の制御は、該エンジン1の回転数に同期させる必要があり、このため、本実施形態では、前記基本燃料噴射量算出部8、第1補正係数算出部9、第2補正係数算出部10、及びフィードバック制御部14の処理は、エンジン1のクランク角周期(所謂TDC)に同期した制御サイクルで行うようにしている。また、この場合、LAFセンサ5やO2センサ6等の各種センサの出力データの読込もクランク角周期(所謂TDC)に同期した制御サイクルで行うようにしている。
【0188】
一方、前記空燃比操作量決定部13による触媒装置3の上流側の排ガスの目標空燃比KCMDの決定処理は、触媒装置3に存する無駄時間や演算負荷等を考慮すると一定周期の制御サイクルで行うことが好ましい。このため、本実施形態では、空燃比操作量決定部13における前述したような処理やその処理のために必要な前記偏差出力kact,VO2を算出する前記減算処理部11,12の処理は一定周期(例えば30〜100ms)の制御サイクルで行うようにしている。
【0189】
尚、この一定周期は、制御対象である触媒装置3の種類や反応速度、容積等に応じて決定すればよい。また、本実施形態では、前記空燃比操作量決定部13による処理を行っているような運転状態(より具体的にはエンジン回転数の状態)において、上記一定周期の時間間隔が前記クランク角周期(TDC)の時間間隔よりも大きくなるように設定している。
【0190】
以上のことを前提として、まず、図6及び図7のフローチャートを参照して、エンジン1の燃料供給量の制御のためのエンジン1の各気筒毎の出力燃料噴射量#nTout(n=1,2,3,4)の算出処理について説明する。制御ユニット7は、各気筒毎の出力燃料噴射量#nToutをエンジン1のクランク角周期と同期した制御サイクルで次のような処理を行う。
【0191】
まず、図6を参照して、制御ユニット7は前記LAFセンサ5及びO2センサ6を含む各種センサの出力を読み込む(STEPa)。この場合、LAFセンサ5の出力KACT及びO2センサ6の出力VO2/OUTはそれぞれ過去に得られたものを含めて時系列的に図示しないメモリに記憶保持される。
【0192】
次いで、基本燃料噴射量算出部8によって、前述の如くエンジン1の回転数NE及び吸気圧PBに対応する燃料噴射量をスロットル弁の有効開口面積に応じて補正してなる基本燃料噴射量Timが求められ(STEPb)、さらに、第1補正係数算出部9によって、エンジン1の冷却水温やキャニスタのパージ量等に応じた第1補正係数KTOTALが算出される(STEPc)。
【0193】
次いで、制御ユニット7は、空燃比操作量決定部13で生成される目標空燃比KCMDを使用するか否か(ここでは、空燃比操作量決定部13のON/OFFという)の判別処理を行って、空燃比操作量決定部13のON/OFFを規定するフラグf/prism/onの値を設定する(STEPd)。尚、フラグf/prism/onの値は、それが「0」のとき、空燃比操作量決定部13で生成される目標空燃比KCMDを使用しないこと(OFF)を意味し、「1」のとき、空燃比操作量決定部13で生成される目標空燃比KCMDを使用すること(ON)を意味する。
【0194】
上記の判別処理では、図7に示すように、O2センサ6及びLAFセンサ5が活性化しているか否かの判別(STEPd−1,d−2)が行われ、いずれかが活性化していない場合には、空燃比操作量決定部13の処理に使用するO2センサ6及びLAFセンサ5の検出データを精度よく得ることができないため、フラグf/prism/onの値を「0」にセットする(STEPd−10)。
【0195】
また、エンジン1のリーン運転中(希薄燃焼運転)であるか否か(STEPd−3)、エンジン1の始動直後の触媒装置3の早期活性化を図るためにエンジン1の点火時期が遅角側に制御されているか否か(STEPd−4)、エンジン1のスロットル弁が全開であるか否か(STEPd−5)、及びエンジン1への燃料供給の停止中であるか否か(STEPd−6)の判別が行われ、これらのいずれかの条件が成立している場合には、空燃比操作量決定部13で生成される目標空燃比KCMDを使用してエンジン1の燃料供給を制御することは好ましくないので、フラグf/prism/onの値を「0」にセットする(STEPd−10)。
【0196】
さらに、エンジン1の回転数NE及び吸気圧PBがそれぞれ所定範囲内にあるか否かの判別が行われ(STEPd−7,d−8)、いずれかが所定範囲内にない場合には、空燃比操作量決定部13で生成される目標空燃比KCMDを使用してエンジン1の燃料供給を制御することは好ましくないので、フラグf/prism/onの値を「0」にセットする(STEPd−10)。
【0197】
そして、STEPd−1,d−2,d−7,d−8の条件が満たされ、且つ、STEPd−3,d−4,d−5,d−6の条件が成立していない場合に、空燃比操作量決定部13で生成される目標空燃比KCMDをエンジン1の燃料供給の制御に使用すべく、フラグf/prism/onの値を「1」にセットする(STEPd−9)。
【0198】
図6に戻って、上記のようにフラグf/prism/onの値を設定した後、制御ユニット7は、フラグf/prism/onの値を判断し(STEPe)、f/prism/on=1である場合には、空燃比操作量決定部13で生成された最新の目標空燃比KCMDを読み込む(STEPf)。また、f/prism/on=0である場合には、目標空燃比KCMDを所定値に設定する(STEPg)。この場合、目標空燃比KCMDとして設定する所定値は、例えばエンジン1の回転数NEや吸気圧PBからあらかじめ定めたマップ等を用いて決定する。
【0199】
次いで、制御ユニット7は、前記局所的フィードバック制御部16において、前述の如くオブザーバ21によりLAFセンサ5の出力KACTから推定した各気筒毎の実空燃比#nA/Fに基づき、PID制御器22により、各気筒毎のばらつきを解消するようにフィードバック補正係数#nKLAFを算出し(STEPh)、さらに、大局的フィードバック制御部15により、フィードバック補正係数KFBを算出する(STEPi)。
【0200】
この場合、大局的フィードバック制御部15は、前述の如く、PID制御器17により求められるフィードバック操作量KLAFと、適応制御器18により求められるフィードバック操作量KSTRを目標空燃比KCMDで除算してなるフィードバック操作量kstrとから、切換部20によってエンジン1の運転状態等に応じていずれか一方のフィードバック操作量KLAF又はkstrを選択し(通常的には適応制御器18側のフィードバック操作量kstrを選択する)、それを燃料噴射量を補正するためのフィードバック補正量数KFBとして出力する。 尚、フィードバック補正係数KFBを、PID制御器17側のフィードバック操作量KLAFから適応制御器18側のフィードバック操作量kstrに切り換える際には、該補正係数KFBの急変を回避するために、適応制御器18は、その切換えの際の制御サイクルに限り、補正係数KFBを前回の補正係数KFB(=KLAF)に保持するように、フィードバック操作量KSTRを求める。同様に、補正係数KFBを、適応制御器18側のフィードバック操作量kstrからPID制御器17側のフィードバック操作量KLAFに切り換える際には、PID制御器17は、自身が前回の制御サイクルで求めたフィードバック操作量KLAFが、前回の補正係数KFB(=kstr)であったものとして、今回の補正係数KLAFを算出する。
【0201】
上記のようにしてフィードバック補正係数KFBが算出された後、さらに、前記STEPfあるいはSTEPgで決定された目標空燃比KCMDに応じた第2補正係数KCMDMが第2補正係数算出部10により算出される(STEPj)。
【0202】
次いで、制御ユニット7は、前述のように求められた基本燃料噴射量Timに、第1補正係数K TOTAL、第2補正係数KCMDM 、フィードバック補正係数KFB、及び各気筒毎のフィードバック補正係数#nKLAFを乗算することで、各気筒毎の出力燃料噴射量#nToutを求める(STEPk)。そして、この各気筒毎の出力燃料噴射量#nToutが、付着補正部23によって、エンジン1の吸気管の壁面付着を考慮した補正を施された後(STEPm)、エンジン1の図示しない燃料噴射装置に出力される(STEPn)。
【0203】
そして、エンジン1にあっては、各気筒毎の出力燃料噴射量#nToutに従って、各気筒への燃料噴射が行われる。
【0204】
以上のような各気筒毎の出力燃料噴射量#nToutの算出及びそれに応じたエンジン1への燃料噴射がエンジン1のクランク角周期に同期したサイクルタイムで逐次行われ、これによりLAFセンサ5の出力KACT(触媒装置3に進入する排ガスの検出空燃比)が、目標空燃比KCMDに収束するように、エンジン1の燃料供給量(燃料噴射量)が制御される。この場合、特に、フィードバック補正係数KFBとして、適応制御部18側のフィードバック操作量kstrを使用している状態では、エンジン1の運転状態の変化や特性変化等の挙動変化に対して、高い安定性を有して、LAFセンサ5の出力KACTが迅速に目標空燃比KCMDに収束制御される。
【0205】
一方、前述のようなエンジン1の燃料制御と並行して、前記空燃比操作量決定部13は、一定周期の制御サイクルで図8のフローチャートに示すメインルーチン処理を行う。
【0206】
すなわち、図8のフローチャートを参照して、空燃比操作量決定部13は、まず、自身の演算処理(前記同定器25、推定器26及びスライディングモード制御器27の演算処理)を実行するか否かの判別処理を行って、その実行の可否を規定するフラグf/prism/calの値を設定する(STEP1)。尚、フラグf/prism/calの値は、それが「0」のとき、空燃比操作量決定部13における演算処理を行わないことを意味し、「1」のとき、空燃比操作量決定部13における演算処理を行うことを意味する。
【0207】
上記の判別処理は、図9のフローチャートに示すように行われる。
【0208】
すなわち、前記図6のSTEPdの場合と同様に、O2センサ6及びLAFセンサ5が活性化しているか否かの判別(STEP1−1,1−2)が行われ、いずれかが活性化していない場合には、空燃比操作量決定部13の演算処理に使用するO2センサ6及びLAFセンサ5の検出データを精度よく得ることができないため、フラグf/prism/calの値を「0」にセットする(STEP1−6)。さらにこのとき、同定器25の後述する初期化を行うために、その初期化を行うか否かを規定するフラグf/id/resetの値を「1」にセットする(STEP1−7)。ここで、フラグf/id/resetの値は、それが「1」であるとき、初期化を行うことを意味し、「0」であるとき、初期化を行わないことを意味する。
【0209】
また、エンジン1のリーン運転中(希薄燃焼運転)であるか否か(STEP1−3)、及びエンジン1の始動直後の触媒装置3の早期活性化を図るためにエンジン1の点火時期が遅角側に制御されているか否か(STEP1−4)の判別が行われ、これらのいずれかの条件が成立している場合には、O2センサ6の出力VO2/OUTを目標値VO2/TARGETに整定させるような目標空燃比KCMDを算出しても、それがエンジン1の燃料制御に使用されることはないので、フラグf/prism/calの値を「0」にセットし(STEP1−6)、さらに同定器25の初期化を行うために、フラグf/id/resetの値を「1」にセットする(STEP1−7)。
【0210】
図8に戻って、上記のような判別処理を行った後、空燃比操作量決定部13は、さらに、同定器25による前記ゲイン係数a1,a2,b1の同定(更新)処理を実行するか否かの判別処理を行って、その実行の可否を規定するフラグf/id/calの値を設定する(STEP2)。尚、フラグf/id/calの値は、それが「0」のとき、同定器25による前記ゲイン係数a1,a2,b1の同定(更新)処理を行わないことを意味し、「1」のとき、同定(更新)処理を行うことを意味する。
【0211】
このSTEP2の判別処理は、図10のフローチャートに示すように行われる。
【0212】
すなわち、エンジン1のスロットル弁が全開であるか否か(STEP2−1)、エンジン1への燃料供給の停止中であるか否か(STEP2−2)、及びエンジン1のアイドル運転中であるか否か(STEP2−3)の判別が行われ、これらのいずれかの条件が成立している場合には、前記ゲイン係数a1,a2,b1を適正に同定することが困難であるため、フラグf/id/calの値を「0」にセットする(STEP2−5)。そして、STEP2−1〜2−3のいずれの条件も成立していない場合には、同定器25による前記ゲイン係数a1,a2,b1の同定(更新)処理を実行すべくフラグf/id/calの値を「1」にセットする(STEP2−4)。
【0213】
図8に戻って、空燃比操作量決定部13は、次に、前記減算処理部11,12からそれぞれ最新の前記偏差出力kact(k)(=KACT−FLAF/BASE)及びVO2(k)(=VO2/OUT−VO2/TARGET)を取得する(STEP3)。この場合、減算処理部11,12は、前記図6のSTEPaにおいて取り込まれて図示しないメモリに記憶されたLAFセンサ5の出力KACT及びO2センサ6の出力VO2/OUTの時系列データの中から、最新のものを選択して前記偏差出力kact(k)及びVO2(k)を算出し、それを空燃比操作量決定部13に与える。そして、該空燃比操作量決定部13に与えられた偏差出力kact(k)及びVO2(k)は、該空燃比操作量決定部13内において、過去に与えられたものを含めて時系列的に図示しないメモリに記憶保持される。
【0214】
次いで、空燃比操作量決定部13は、前記STEP1で設定されたフラグf/prism/calの値を判断し(STEP4)、f/prism/cal=0である場合、すなわち、空燃比操作量決定部13の演算処理を行わない場合には、スライディングモード制御器27で求めるべき前記対象排気系EへのSLD操作入力uslを強制的に所定値に設定する(STEP12)。この場合、該所定値は、例えばあらかじめ定めた固定値(例えば「0」)あるいは前回の制御サイクルで決定したSLD操作入力uslの値とする。尚、このようにSLD操作入力uslを所定値とした場合において、空燃比操作量決定部13は、その所定値のSLD操作入力uslに前記基準値FLAF/BASEを加算することで、今回の制御サイクルにおける目標空燃比KCMDを決定し(STEP13)、今回の制御サイクルの処理を終了する。
【0215】
一方、STEP4の判断で、f/prism/cal=1である場合、すなわち、空燃比操作量決定部13の演算処理を行う場合には、空燃比操作量決定部13は、前記同定器25による演算処理を行う(STEP5)。
【0216】
この同定器25による演算処理は図11のフローチャートに示すように行われる。
【0217】
すなわち、同定器25は、まず、前記STEP2で設定されたフラグf/id/calの値を判断する(STEP5−1)。このときf/id/cal=0であれば、前述の通り同定器25によるゲイン係数a1,a2,b1の同定処理を行わないので、直ちに図8のメインルーチンに復帰する。
【0218】
一方、f/id/cal=1であれば、同定器25は、さらに該同定器25の初期化に係わる前記フラグf/id/resetの値(これは、前記STEP1等でその値が設定される)を判断し(STEP5−2)、f/id/reset=1である場合には、同定器25の初期化を行う(STEP5−3)。この初期化では、前記同定ゲイン係数a1ハット,a2ハット,b1ハットの各値があらかじめ定めた初期値に設定され(式(3)の同定ゲイン係数ベクトルΘの初期化)、また、前記式(9)の行列P(対角行列)の各成分があらかじめ定めた初期値に設定される。さらに、フラグf/id/resetの値は「0」にリセットされる。
【0219】
次いで、同定器25は、現在の同定ゲイン係数a1(k-1)ハット,a2(k-1)ハット,b1(k-1)ハットを用いて表される対象排気系Eの離散系モデル(前記式(2)参照)におけるO2センサ6の前記同定偏差出力VO2(k)ハットを、前記STEP3で制御サイクル毎に取得される偏差出力VO2及びkactの過去のデータVO2(k-1),VO2(k-2),kact(k-d-1)と、上記同定ゲイン係数a1(k-1)ハット,a2(k-1)ハット,b1(k-1) ハットとを用いて前記式(2)あるいはこれと等価の前記式(5)により算出する(STEP5−4)。
【0220】
さらに同定器25は、新たな同定ゲイン係数a1ハット,a2ハット,b1ハットを決定する際に使用する前記ベクトルKθ(k)を式(8)により算出した後(STEP5−5)、以下に説明する同定器25のマネージメント処理を行う(STEP5−6)。
【0221】
ここで、まず、O2センサ6の出力VO2/OUTあるいは偏差出力VO2の挙動と、前記対象排気系Eの離散系モデルのゲイン係数a1,a2,b1の同定器25による同定との関係について説明しておく。
【0222】
図12を参照して、O2センサ6の出力VO2/OUTあるいは偏差出力VO2は、触媒装置3を通過した排ガスの空燃比を示すものであり、この排ガスの空燃比は、一般に、図示のようにリーン側からリッチ側への変化が比較的急激に行われ(O2センサ6の出力VO2/OUTあるいは偏差出力VO2の時間的な変化度合いが比較的大きい)、リッチ側からリーン側への変化は比較的緩やかに行われる(O2センサ6の出力VO2/OUTあるいは偏差出力VO2の時間的な変化度合いが比較的小さい)。そして、本願発明者等の知見によれば、対象排気系Eの離散系モデルのゲイン係数a1,a2,b1をO2センサ6の出力VO2/OUTあるいは偏差出力VO2を用いて同定する場合、O2センサ6の出力VO2/OUTあるいは偏差出力VO2の時間的変化度合いが比較的小さい状態では、ゲイン係数a1,a2,b1の同定値が小さくなり過ぎる等して、適正な同定ゲイン係数a1ハット,a2ハット,b1ハットが得られない場合が生じやすい。
【0223】
そこで、本実施形態では、同定器25による前記ゲイン係数a1,a2,b1の同定(更新)を、O2センサ6の出力VO2/OUTあるいは偏差出力VO2により示される空燃比が、リーン側からリッチ側に変化する挙動状態において行うようにしており、前記マネージメント処理は、上記のような挙動状態を特定するための処理である。
【0224】
一方、図13を参照して、適応スライディングモード制御を用いた本実施形態の制御によれば、O2センサ6の偏差出力VO2の前記状態量X(VO2(k),VO2(k-1))は、その状態量Xの初期状態が例えば図中の点Qであるとしたとき、該状態量Xは、前記超平面σ=0に対して軌跡線Wで示すように変化する。そして、この場合、同図において、基本的には状態量Xが超平面σ=0の上側で変化している状態(このとき状態量Xにより規定される線形関数σの値は正となる)が、触媒装置3を通過した排ガスの空燃比のリーン側からリッチ側への変化状態であり、状態量Xが超平面σ=0の下側で変化している状態(このとき状態量Xにより規定される線形関数σの値は負となる)が、リッチ側からリーン側への変化状態である。
【0225】
従って、触媒装置3を通過した排ガスの空燃比がリーン側からリッチ側に変化する挙動状態であるか否かの判断は、基本的には、線形関数σの値が正であるか否かによって判断することができる。但し、このように線形関数σの値が正であるか否かによって排ガスの空燃比がリーン側からリッチ側に変化する挙動状態であるか否かを判断するようにすると、状態量Xが超平面σ=0上から僅かに変化しただけで、排ガスの空燃比がリーン側からリッチ側に変化する挙動状態であるか否かの判断結果が変わってしまい、その判断結果に応じて前記ゲイン係数a1,a2,b1の同定(更新)処理を安定して行う上では好ましくない。
【0226】
このため、本実施形態では、次式(50)により偏差出力VO2の時系列データを用いて定義されるマネージメント関数γを導入し、
【0227】
【数50】
【0228】
このマネージメント関数γの係数m1,m2,m3を、γ=0により表されるマネージメント用超平面(この場合は直線)が、前記図13に示したように、スライディングモード制御用の超平面σ=0から若干上側(σ>0の領域)に存するように設定した。尚、本実施形態では、線形関数σの係数s1を「1」に設定していることに合わせて、マネージメント関数γの係数m1は「1」に設定している。
【0229】
このようなマネージメント関数γを導入すると、γ≧0となる状態では、確実に排ガスの空燃比がリーン側からリッチ側に変化する挙動状態となり、この挙動状態であるか否かの判断は、マネージメント関数γの値が正(「0」を含む)であるか否かによって安定して行うことができる。
【0230】
前記STEP5−6のマネージメント処理は、上記のように定義されたマネージメント関数γを用いて、O2センサ6の偏差出力VO2により示される排ガスの空燃比がリーン側からリッチ側に変化する挙動状態、すなわち、同定器25による前記ゲイン係数a 1,a2,b1の同定(更新)に好適な挙動状態であるか否かの判断を行うものであり、その処理は具体的には次のように行われる。
【0231】
すなわち、図14のフローチャートを参照して、同定器25は、前記STEP3(図8参照)で取得された最新の偏差出力VO2(k)と前回の制御サイクルにおける偏差出力VO2(k-1)とを用いて、式(50)によりマネージメント関数γの値を算出する(STEP5−6−1)。
【0232】
次いで、同定器25は、γ≧0であるか否かを判断し(STEP5−6−2)、γ≧0である場合には、排ガスの空燃比がリーン側からリッチ側に変化する挙動状態であるか否かを示すフラグf/id/mngの値を「1」に設定し(STEP5−6−3)、γ<0である場合には、フラグf/id/mngの値を「0」に設定する(STEP5−6−4)。
【0233】
これにより、排ガスの空燃比がリーン側からリッチ側に変化する挙動状態であるか否か、すなわち、同定器25による前記ゲイン係数a1,a2,b1の同定(更新)に好適な挙動状態であるか否かが、f/id/mngの値により示されることとなる。
【0234】
図11の説明に戻って、同定器25は、前述のようにマネージメント処理を行った後、その処理において設定されるフラグf/id/mngの値を判断し(STEP5−7)、f/id/mng=1である場合、すなわち、触媒装置3を通過した排ガスの空燃比がリーン側からリッチ側に変化する挙動状態(ゲイン係数a1,a2,b1の同定(更新)に好適な挙動状態)である場合には、前記同定誤差id/e(離散系モデル上でのO2センサの同定偏差出力VO2ハットと、実際の偏差出力VO2との偏差。式(6)参照)を算出し(STEP5−8)、f/id/mng=0である場合には、前記同定誤差id/eの値を強制的に「0」とする(STEP5−9)。
【0235】
そして、同定器25は、STEP5−8あるいはSTEP5−9で得られた同定誤差id/eと、前記STEP5−5で算出されたKθとを用いて前記式(7)により新たな同定ゲイン係数ベクトルΘ(k)、すなわち、新たな同定ゲイン係数a1(k)ハット,a2(k)ハット,b1(k)ハットを算出する(STEP5−10)。
【0236】
ここで、前記STEP5−8における同定誤差id/eは、基本的には、前記式(6)に従って算出すればよいのであるが、本実施形態では、例えば図15(a)にブロック図で示すように前記STEP3(図8参照)で制御サイクル毎に取得する偏差出力VO2と、前記STEP5−4で制御サイクル毎に算出する同定偏差出力VO2ハットとにそれぞれ同一特性のフィルタリングを施した上で、STEP5−8における同定誤差id/eの算出を行う。
【0237】
すなわち、図16を参照して、前記触媒装置3を含む対象排気系Eの入力変化(LAFセンサ5の出力KACTあるいは偏差出力kactの変化)に対する、該対象排気系Eの出力変化(O2センサ6の出力VO2/OUTあるいは偏差出力VO2の変化)のゲインの周波数特性は、一般に図に実線で示すように低周波数帯Cの周波数通過特性を有するローパス特性となる。従って、O2センサ6の出力VO2/OUT(対象排気系Eの出力)を目標値VO2/TARGETに制御すべく前記スライディングモード制御器27により目標空燃比KCMD(対象排気系Eの入力の目標値)を決定する上では、上記低周波数帯Cを重視する必要がある。
【0238】
そして、スライディングモード制御器27は、基本的には、前述の通り同定器25で前記式(7)により同定した離散系モデルのゲイン係数a1,a2,b1、すなわち同定ゲイン係数a1ハット,a2ハット,b1ハットを用いて目標空燃比KCMDを求めるものであるため、該同定ゲイン係数a1ハット,a2ハット,b1ハットにより定まる離散系モデルの周波数特性も、実際の対象排気系Eの周波数特性と同じような傾向の周波数特性(低周波数帯Cの周波数通過特性を有するローパス特性)となることが好ましい。
【0239】
一方、同定器25による離散系モデルのゲイン係数a1,a2,b1の同定のための演算処理(式(7)〜(9)を参照)は、対象排気系Eがローパス特性を有するため、例えば図16に仮想線で示すように前記低周波数帯Cよりも高周波側に重みを有する。このため、前記同定偏差出力VO2ハット及び偏差出力VO2をそのまま用いて求めた同定誤差id/eに応じて同定ゲイン係数a1ハット,a2ハット,b1ハットを求めるようにすると、その同定ゲイン係数a1ハット,a2ハット,b1ハットにより定まる離散系モデルの周波数特性が、実際の対象排気系Eの周波数特性と適合せず、前記低周波数帯C外での対象排気系Eのゲイン特性を重視した特性となる。特に、該低周波数帯Cにおける離散系モデルのゲインが実際の対象排気系Eのゲインよりも小さなものとなりやすい。
【0240】
そこで、本実施形態では、図16に一点鎖線で示すように低周波数帯Cに重みを有する特性(ローパス特性)のフィルタリングを、偏差出力VO2と同定偏差出力VO2ハットとに施した上で、STEP5−8における同定誤差id/eの算出を行う。
【0241】
尚、このようなローパス特性のフィルタリング処理は、前記STEP3(図8)で取得する偏差出力VO2と前記STEP5−4で算出する同定偏差出力VO2ハットとをそれぞれ時系列的に記憶保持しておき、該偏差出力VO2及び同定偏差出力VO2ハットの時系列データのそれぞれについて、制御サイクル毎に、現在から過去に逆上った所定数のデータの加算平均あるいは重み付き加算平均を算出することで行われる。これは、ディジタルフィルタの一手法で、一般に移動平均処理といわれる手法である。そして、STEP5−8における同定誤差id/eの算出は、上記のような移動平均処理で得られた偏差出力VO2のフィリタリング値から同定偏差出力VO2ハットのフィルタリング値を減算することで行われる。
【0242】
このようなフィルタリング処理を行うことによって、該同定誤差id/eから前記式(7)によって求められる同定ゲイン係数a1ハット,a2ハット,b1ハットにより定まる離散系モデルの周波数特性を、例えば図17(a)に示すように、実際の対象排気系Eの周波数特性と同じような傾向の周波数特性にすることができる。
【0243】
この場合、本実施形態では、さらに、同図17(a)に示すように、離散系モデルの各周波数におけるゲインが対象排気系Eの各周波数におけるゲインよりも全体的に若干大きくなるように前記フィルタリングの重み特性を設定しておく。このようにすることで、離散系モデル及び対象排気系Eにおいて、ある出力変化(具体的にはO2センサ6の出力VO2/OUTを目標値VO2/TARGETに一致させるような出力変化)を生ぜしめる入力変化(LAFセンサ5の出力KACTあるいは偏差出力kactの変化)は、離散系モデルの方が対象排気系Eよりも小さくなる。このため、このような離散系モデルの同定ゲイン係数a1ハット,a2ハット,b1ハットを用いて前記スライディングモード制御器27により対象排気系Eに与えるべき入力として求められる前記SLD操作入力uslは各周波数において、小さめの値となり、O2センサ6の出力VO2/OUTの目標値VO2/TARGETへの収束制御の安定性を高めることができる。
【0244】
尚、本実施形態では、図17(a)のように、離散系モデルの各周波数におけるゲインが対象排気系Eの各周波数におけるゲインよりも全体的に若干大きくなるように前記フィルタリングの重み特性を設定したが、例えば図17(b)に示すように、前記低周波数帯Cにおける離散系モデルのゲインの周波数特性が対象排気系Eのゲインの周波数特性と略同一となり、且つ、低周波数帯Cよりも高周波側の周波数帯では、図17(a)の場合と同様に、離散系モデルのゲインが対象排気系Eのゲインよりも若干大きくなるように前記フィルタリングの重み特性を設定しておくようにしてもよい。
【0245】
このようにすると、低周波数帯Cよりも高周波側(対象排気系Eのゲインが比較的小さい周波数帯)では、図17(a)の場合と同様にO2センサ6の出力VO2/OUTの目標値VO2/TARGETへの収束制御の安定性が高まると同時に、その制御上重要な低周波数帯Cでは、実際の対象排気系Eの特性に適合した前記SLD操作入力uslをスライディングモード制御器27によって求めることができ、O2センサ6の出力VO2/OUTの目標値VO2/TARGETへの収束の迅速な追従性(速応性)を確保することができる。
【0246】
また、本実施形態では、前記図15(a)にブロック図で示したように、偏差出力VO2 と同定偏差出力VO2 ハットとに前記フィルタリングを施した上で、同定誤差id/eの算出を行うようにしたが、例えば図15(b)に示すように、前記STEP5−4で同定偏差出力VO2ハットを算出する前に、その算出に使用する偏差出力kact,VO2に同一特性のフィルタリングを施しておき、それらのフィルタリング値から前記式(5)により算出した同定偏差出力VO2ハットと、先にフィルタリングを施した偏差出力VO2のフィルタリング値とから同定誤差id/eを算出するようにしてもよい。あるいは、図15(c)に示すように偏差出力kact,VO2をそのまま用いて算出した同定偏差出力VO2ハットと、偏差出力VO2から前記式(6)をそのまま用いて同定誤差id/eを算出した後に、該同定誤差id/eにフィルタリングを施すようにしてもよい。つまり、同定誤差id/eの算出に際しての前述したようなフィルタリング処理は、結果的に偏差出力VO2と同定偏差出力VO2ハットとに同一特性のフィルタリングが施されていればよく、そのフィルタリング処理のタイミングは、任意に選択することが可能である。
【0247】
図11の説明に戻って、同定器25は、同定ゲイン係数a1ハット,a2ハット,b1ハット(同定ゲイン係数ベクトルΘ)を算出した後、その評価処理を行う(STEP5−11)。この評価処理では、図18のフローチャートに示すように、同定器25は、まず同定ゲイン係数a1ハット,a2ハット,b1ハットのうち、同定ゲイン係数b1ハットの値が所定範囲内にあるか否かを判断することで、同定器25の同定処理の安定性を判断し(STEP5−11−1)、同定ゲイン係数b1ハットの値が所定範囲内にある場合には、さらに、前記STEP5−8で算出された同定誤差id/eの大きさが所定値ε0以下の十分小さなものとなったか否か(id/eがほぼ「0」に収束して、同定ゲイン係数a1ハット,a2ハット,b1ハットがほぼ確定した状態になったか否か)を判断する(STEP5−11−3)。このとき、|id/e|≦ε0であれば、そのまま図11のフローチャートの処理に復帰する。
【0248】
一方、STEP5−11−1の判断で、同定ゲイン係数b1ハットの値が所定範囲内に無い場合には、同定器25によるゲイン係数a1,a2,b1の同定処理が不安定で、適正な同定ゲイン係数a1ハット,a2ハット,b1ハットを求めることが困難な状態であると考えられるので、前記STEP5−3の場合と同様に同定器25の初期化を行い(STEP5−11−2)、さらに、後述のスライディングモード制御の安定性の判断の際に使用するタイマカウンタtm(カウントダウンタイマ)の値を所定の初期値TMにセットする(タイマカウンタtmの起動。STEP5−11−4)。また、STEP5−11−3の判断で、|id/e|>ε0である場合、すなわち、同定ゲイン係数a1ハット,a2ハット,b1ハットがまだ十分に確定していない状態では、STEP5−11−4の処理を行って、前記タイマカウンタtmの値を初期値TMにセットする(タイマカウンタtmの起動)。
【0249】
尚、本実施形態では、同定ゲイン係数a1ハット,a2ハット,b1ハットのうち、同定ゲイン係数b1ハットの値が所定範囲内にあるか否かにより同定器25の同定処理の安定性を判断したが、他の同定ゲイン係数a1ハット,a2ハットについても同様の評価を行い、それらの値が不適切なものである場合にも同定器25の同定処理が不安定であるとして、STEP5−11−2及びSTEP5−11−4の処理を行うようにしてもよい。
【0250】
図11の説明に戻って、前述のように同定ゲイン係数ベクトルΘの評価処理を行った後、同定器25は、前記STEP5−6で設定されるフラグf/id/mngの値を判断し(STEP5−12)、f/id/mng=1である場合、すなわち、触媒装置3を通過した排ガスの空燃比がリーン側からリッチ側に変化する挙動状態(ゲイン係数a1,a2,b1の同定(更新)に好適な挙動状態)である場合には、次回の制御サイクルの処理のために前記行列P(k)を前記式(9)により算出し(STEP5−13)、図8のメインルーチンの処理に復帰する。また、f/id/mng=0である場合には、次回の制御サイクルの処理の際に使用する前記行列P(k)を今現在の行列P(k-1) に維持し(STEP5−14)、図8のメインルーチンの処理に復帰する。
【0251】
以上が図8のSTEP5における同定器25の演算処理である。
【0252】
図8のメインルーチン処理の説明に戻って、前述の通り同定器25の演算処理が行われた後、空燃比操作量決定部13はゲイン係数a1,a2,b1を決定する(STEP6)。この処理では、図19のフローチャートに示すように、前記STEP2で設定されたフラグf/id/calの値が判断され(STEP6−1)、f/id/cal=1である場合、すなわち、同定器25によるゲイン係数a1,a2,b1の同定処理を行った場合には、ゲイン係数a1,a2,b1の値として、それぞれ前記STEP5−10(図11参照)で前述の通り同定器25により求められた同定ゲイン係数a1ハット,a2ハット,b1ハットにそれぞれ所定のスケーリング係数g1,g2,g3によりスケーリングを施したものを設定する(STEP6−2)。尚、本実施形態ではスケーリング係数g1,g2,g3の値は、いずれも「1」としている。
【0253】
また、f/id/cal=0である場合、すなわち、同定器25によるゲイン係数a1,a2,b1の同定処理を行わなかった場合には、ゲイン係数a1,a2,b1の値をそれぞれあらかじめ定めた所定値とする(STEP6−3)。
【0254】
次いで、空燃比操作量決定部13は、図8のメインルーチンにおいて、前記推定器26による演算処理(推定偏差出力VO2バーの算出処理)を行う(STEP7)。
【0255】
この推定器26の演算処理は図20のフローチャートに示すように行われる。すなわち、推定器26は、前記STEP6で決定されたゲイン係数a1,a2,b1を用いて、前記式(13)で使用する係数α1,α2,βj(j=1〜d)を前述したように算出する(式(11)、(12)参照)(STEP7−1)。さらに、推定器26は、式(13)で使用するLAFセンサ5の偏差出力kact(前記図8のSTEP3で取得されたもの)のフィルタリング処理(ローパス特性のフィルタリング)を行った後(STEP7−2)、その偏差出力kactのフィルタリング値の時系列データと、O2センサの偏差出力VO2の時系列データ(前記図8のSTEP3で取得されたもの)とSTEP7−1で算出した係数α1,α2,βjとを用いて前記式(13)により、推定偏差出力VO2(k+d)バー(今回の制御サイクルの時点から無駄時間d後の偏差出力VO2の推定値)を算出する(STEP7−3)。
【0256】
ここで、前記STEP7−2においてLAFセンサ5の偏差出力kactのフィルタリングを行うのは次の理由による。すなわち、前述したように触媒装置3を含む対象排気系Eはローパス特性の周波数特性を有するため(図16参照)、該対象排気系Eの出力としての前記O2センサ6の出力VO2/OUTを目標値VO2/TARGETに制御する上では、前記低周波数帯C(図6参照)を重視する必要がある。従って、スライディングモード制御器27が前記SLD操作入力uslを決定するために用いる推定偏差出力VO2バーを前記式(13)により求めるに際しても、低周波数帯C(図6参照)を重視することが好ましい。この場合、推定偏差出力VO2バーを求めるために式(13)で使用するO2センサ6の偏差出力VO2及びLAFセンサ5の偏差出力kactのうち、偏差出力VO2は、対象排気系Eがローパス特性であるために、高周波成分をほとんど含まないが、偏差出力kactは一般に高周波成分も含みやすい。このために、本実施形態では、LAFセンサ5の偏差出力kactのフィルタリング、すなわち、偏差出力kactの高周波成分の除去を行った上で、前記式(13)により推定偏差出力VO2バーを求めるようにしている。
【0257】
尚、上記のようなフィルタングは、前記同定器25におけるフィルタリングの場合と同様に、移動平均処理によって行われる。
【0258】
図8の説明に戻って、空燃比操作量決定部13は、次に、スライディングモード制御器27によって、前記SLD操作入力uslを算出する(STEP8)。
【0259】
このSLD操作入力uslの算出は、図21のフローチャートに示すように行われる。
【0260】
すなわち、スライディングモード制御器27は、まず、前記STEP8で推定器2により求められた推定偏差出力VO2バーの時系列データ(詳しくは、今回の制御サイクルで求められたVO2(k+d)バーと、前回の制御サイクルで求められたVO2(k+d-1)バー)を用いて、前記式(35)により定義された線形関数σバーの今回の制御サイクルから無駄時間d後の値σ(k+d)バー(これは、式(14)で定義された線形関数σの無駄時間d後の推定値に相当する)を算出する(STEP8−1)。
【0261】
次いで、スライディングモード制御器27は、上記STEP8−1で制御サイクル毎に算出されるσ(k+d)バーを累積的に加算していく(前回の制御サイクルで求められた加算結果に今回の制御サイクルで算出されたσ(k+d)を加算する)ことで、σ(k+d)バーの積算値(これは式(37)の右端の項に相当する)を算出する(STEP8−2)。尚、この場合、本実施形態では、σ(k+d)バーの積算値があらかじめ定めた所定範囲内に収まるようにし、σ(k+d)バーの積算値が所定の上限値又は下限値を超えた場合には、それぞれσ(k+d)バーの積算値を該上限値又は下限値に制限するようにしている。これは、σ(k+d)バーの積算値の大きさが過大になると、前記式(37)により求められる適応則入力uadpが過大となって、制御性が損なわれる虞れがあるからである。
【0262】
次いで、スライディングモード制御器27は、前記図8のSTEP6で決定されたゲイン係数b1のリミット処理を次のように行う(STEP8−3)。
【0263】
すなわち、図22のフローチャートを参照して、スライディングモード制御器27はゲイン係数b1の大きさがあらかじめ定めた所定値ε1よりも小さいか否かを判断し(STEP3−1)、|b1|≧ε1である場合には、そのまま図21のフローチャートの処理に復帰する。
【0264】
一方、|b1|<ε1である場合(b1≒0の場合)には、ゲイン係数b1の符号が正(b1=0の場合を含む)であるか否かに応じて、それぞれゲイン係数b1の値をあらかじめ定めた正の所定値(≧ε1)及び負の所定値(≦−ε1)に強制的に制限する。このように、ゲイン係数b1の大きさを制限して、該ゲイン係数b1が過小なものとなるのを防止するのは、前記式(34)、(36)、(37)を参照して明らかなように、ゲイン係数b1が、等価制御入力ueq、到達則入力urch及び適応則入力uadpの算出する際に、分母項として使用するものであるため、該ゲイン係数b1の大きさが小さ過ぎると、過大な等価制御入力ueq、到達則入力urch及び適応則入力uadpが算出されてしまうからである。
【0265】
尚、本実施形態では、ゲイン係数b1の値のみを制限するようにしたが、他のゲイン係数a1,a2の値も制限するようにしてもよい。
【0266】
図21の説明に戻って、上記のようにゲイン係数b1のリミット処理を行った後、スライディングモード制御器27は前記STEP8で推定器2により求められた推定偏差出力VO2バーの時系列データVO2(k+d)バー,VO2(k+d-1)バーと、STEP8−1及び8−2でそれぞれ求められた線形関数の値σ(k+d)バー及びその積算値と、STEP6で決定したゲイン係数a1ハット,a2ハット及びSTEP8−3で決定したゲイン係数b1ハットとを用いて、前記式(34)、(36)、(37)に従って、それぞれ等価制御入力ueq、到達則入力urch及び適応則入力uadpを算出する(STEP8−4)。
【0267】
さらにスライディングモード制御器27は、STEP8−4で求めた等価制御入力ueq、到達則入力urch及び適応則入力uadpを加算することで、前記SLD操作入力usl、すなわち、O2センサ6の出力VO2/OUTを目標値VO2/TARGETに収束させるために必要な対象排気系Eへの入力(LAFセンサ5で検出される排ガスの空燃比と基準値FLAF/BASEとの偏差)を算出する(STEP8−5)。
【0268】
図8に戻って、上記のようにSLD操作入力uslを算出した後、空燃比操作量決定部13は、スライディングモード制御器27による適応スライディングモード制御の安定性の判別処理を行って、該適応スライディングモード制御が安定であるか否かを示すフラグf/sld/stbの値を設定する(STEP8−6)。
【0269】
この安定性の判別処理は図23のフローチャートに示すように行われる。
【0270】
すなわち、空燃比操作量決定部13は、まず、前記STEP8−1で算出される線形関数σバーの今回値σ(k+d)バーと前回値σ(k+d-1)バーとの偏差Δσバー(これは線形関数のσバーの変化速度に相当する)を算出する(STEP9−1)。
【0271】
次いで、空燃比操作量決定部13は、STEP9−1で算出した偏差Δσバーと線形関数σバーの今回値σ(k+d)バーとの積Δσバー・σ(k+d)バー(これはσバーに関するリアプノフ関数σバー2/2の時間微分関数に相当する)があらかじめ定めた所定値ε2(≧0)以下であるか否かを判断する(STEP9−2)。
【0272】
この場合、Δσバー・σ(k+d)バー>ε2となる状態は、σバー2が増加する側で、前記推定偏差出力VO2(k+d),VO2(k+d-1)が前記超平面σ=0から離間する方向へ変移している状態であるので、適応スライディングモード制御が不安定(前記STEP8で算出されるSLD操作入力uslが不適切)であると考えられる。このため、STEP9−2の判断で、Δσバー・σ(k+d)バー>ε2である場合には、適応スライディングモード制御が不安定であるとして、同定器25の初期化を行うべく前記フラグf/id/resetの値を「1」に設定し(STEP9−4)、さらに、前記STEP8で算出されるSLD操作入力uslを用いた目標空燃比KCMDの決定を所定時間、禁止するためにタイマカウンタtm(カウントダウンタイマ)の値を所定の初期値TMにセットする(タイマカウンタtmの起動。STEP9−5)。そして、前記フラグf/sld/stbの値を「0」(f/sld/stb=0は適応スライディングモード制御が不安定であることを示す)に設定する(STEP9−6)。
【0273】
尚、STEP9−2の判断で使用する所定値ε2は理論上は「0」でよいが、確率的外乱の影響を考慮すると、「0」よりも若干大きな値とすることが好ましい。
【0274】
一方、前記STEP9−2の判断で、Δσバー・σ(k+d)バー≦ε2である場合には、空燃比操作量決定部13は、線形関数σバーの今回値σ(k+d)バーがあらかじめ定めた所定範囲内にあるか否かを判断する(STEP9−3)。
【0275】
この場合、線形関数σバーの今回値σ(k+d)バーが、所定範囲内に無い状態は、前記推定偏差出力VO2(k+d),VO2(k+d-1)が前記超平面σ=0から大きく離間している状態であるので、適応スライディングモード制御が不安定(前記STEP8で算出されるSLD操作入力uslが不適切)であると考えられる。このため、STEP9−2の判断で、線形関数σバーの今回値σ(k+d)バーが、所定範囲内に無い場合には、適応スライディングモード制御が不安定であるとして、前述の場合と同様に、STEP9−5〜9−6の処理を行って、同定器25の初期化を行うべく前記フラグf/id/resetの値を「1」に設定すると共に、タイマカウンタtm を起動する。
【0276】
また、STEP9−3の判断で、線形関数σバーの今回値σ(k+d)バーが、所定範囲内にある場合には、空燃比操作量決定部13は、前記タイマカウンタtmを所定時間Δtm分、カウントダウンし(STEP9−7)、さらに、該タイマカウンタtmの値が「0」以下であるか否か、すなわち、タイマカウンタtmを起動してから前記初期値TM分の所定時間が経過したか否かを判断する(STEP9−8)。
【0277】
このとき、tm>0である場合、すなわち、タイマカウンタtmが計時動作中でまだタイムアップしていない場合は、STEP9−2あるいはSTEP9−3の判断で適応スライディングモード制御が不安定であると判断されてから、さほど時間を経過していない状態で、適応スライディングモード制御が不安定なものとなりやすいので、前記STEP9−6の処理を行って前記フラグf/sld/stbの値を「0」に設定する。
【0278】
そして、STEP9−8の判断でtm≦0である場合、すなわち、タイマカウンタtmがタイムアップしている場合には、適応スライディングモード制御が安定であるとして、フラグf/sld/stbの値を「1」(f/sld/stb=1は適応スライディングモード制御が安定であることを示す)に設定する(STEP9−9)。
【0279】
尚、前記タイマカウンタtmは、前記同定器25における同定ゲイン係数a1ハット,a2ハット,b1ハットの前述の評価処理(前記図18のフローチャートの処理)において、前記同定誤差id/eが未収束状態で、同定ゲイン係数a1ハット,a2ハット,b1ハットがまだ十分に確定していない場合でも起動される。このため、STEP9−2あるいはSTEP9−3の条件が満たされた場合であっても、STEP9−8の判断でtm>0となる場合があり、前記フラグf/sld/stbの値が「0」に設定される。これは、同定ゲイン係数a1ハット,a2ハット,b1ハットがまだ十分に確定していない段階では、該同定ゲイン係数a1ハット,a2ハット,b1ハットを用いた適応スライディングモード制御が不安定なものとなり易いからである。
【0280】
以上のような処理によって、スライディングモード制御器27による適応スライディングモード制御の安定性が判断され、不安定であると判断した場合には、フラグf/sld/stbの値が「0」に設定され、安定であると判断した場合には、フラグf/sld/stbの値が「1」に設定される。
【0281】
尚、本実施形態では、適応スライディングモード制御の安定性の判断は、基本的には、前記STEP9−2及び9−3の条件判断で行うようにしたが、いずれか一方の条件判断で行うようにしてもよく、あるいは、線形関数σバーの変化速度に相当する前記偏差Δσバーの大きさ(絶対値)を所定値と比較することで、適応スライディングモード制御の安定性の判断を行うようにすることも可能である。
【0282】
図8に戻って、上記のようにスライディングモード制御器27による適応スライディングモード制御の安定性を示すフラグf/sld/stbの値を設定した後、空燃比操作量決定部13は、フラグf/sld/stbの値を判断する(STEP10)。このとき、f/sld/stb=1である場合、すなわち、適応スライディングモード制御が安定であると判断された場合には、スライディングモード制御器27によって、前記STEP8で算出されたSLD操作入力uslのリミット処理が行われる(STEP11)。このリミット処理では、SLD操作入力uslの値や、その値の変化幅が所定範囲に制限され、STEP8で算出されたSLD操作入力uslの今回値usl(k)が所定の上限値又は下限値を超えている場合には、それぞれ、SLD操作入力uslの値が強制的に該上限値又は下限値に設定される。また、STEP8で算出されたSLD操作入力uslの今回値usl(k)の前回値usl(k-1)からの変化量が所定量を超えている場合には、SLD操作入力uslの値が強制的に前回値usl(k-1) に該所定量を加えた値に設定される。
【0283】
そして、空燃比操作量決定部13は、上記のようなSLD操作入力uslのリミット処理の後、スライディングモード制御器27によって、前記式(38)に従って前記目標空燃比KCMDを算出せしめ(STEP13)、今回の制御サイクルの処理終了する。
【0284】
また、前記STEP10の判断でf/sld/stb=0である場合、すなわち、適応スライディングモード制御が不安定であると判断された場合には、空燃比操作量決定部13は、今回の制御サイクルにおけるSLD操作入力uslの値を強制的に所定値(固定値あるいはSLD操作入力uslの前回値)に設定した後(STEP12)、スライディングモード制御器27によって、前記式(38)に従って前記目標空燃比KCMDを算出せしめ(STEP13)、今回の制御サイクルの処理終了する。
【0285】
尚、STEP13で最終的に決定される目標空燃比KCMDは、制御サイクル毎に図示しないメモリに時系列的に記憶保持される。そして、前記大局的フィードバック制御器17等が、空燃比操作量決定部13で決定された目標空燃比KCMDを用いるに際しては(図6のSTEPfを参照)、上記のように時系列的に記憶保持された目標空燃比KCMDの中から最新のものが選択される。
【0286】
以上説明した内容が本実施形態の装置の詳細な作動である。
【0287】
すなわち、その作動を要約すれば、基本的には空燃比操作量決定部13によって、触媒装置3の下流側のO2センサ6の出力VO2/OUT(これはプラントとしての対象排気系Eの出力に相当する)を目標値VO2/TARGETに収束(整定)させるように、触媒装置3に進入する排ガスの目標空燃比KCMD(これは、対象排気系Eの目標入力に相当する)が逐次決定され、この目標空燃比KCMDに従って、対象排気系Eへの入力(排ガスの空燃比)を生成するアクチュエータとしてのエンジン1の燃料供給量が該目標空燃比KCMD及び触媒装置3の上流側のLAFセンサ5の出力KACTに基づきフィードバック制御される。そして、上記のように触媒装置3の下流側のO2センサ6の出力VO2/OUT を目標値VO2/TARGETに整定させることで、触媒装置3の経時劣化等によらずに、触媒装置3の最適な排ガス浄化性能を確保することができる。
【0288】
この場合、空燃比操作量決定部13は、本来的に外乱等の影響を受けにくいという特性を有するスライディングモード制御を用い、特に外乱等の影響を極力排除するための適応則を付加した適応スライディングモード制御を用いて前記目標空燃比KCMD(対象排気系Eの目標入力)を決定するため、O2センサ6の出力VO2/OUT(対象排気系Eの出力)を目標値VO2/TARGETに整定させる上で的確な目標空燃比KCMDを外乱等の影響を極力抑えて安定して求めることができ、ひいては、O2センサ6の出力VO2/OUTの目標値VO2/TARGETへの制御を安定して精度よく行うことができる。
【0289】
また、空燃比操作量決定部13のスライディングモード制御器27が適応スライディングモード制御により目標空燃比KCMDを決定するに際しては、推定器26により求められた推定偏差出力VO2、すなわち対象排気系Eの無駄時間d後のO2センサ6の偏差出力VO2の推定値を用い、その推定偏差出力VO2により示される無駄時間d後のO2センサ6の出力VO2/OUT(対象排気系Eの出力)の推定値を目標値VO2/TARGETに収束させるように目標空燃比KCMD(対象排気系Eの目標入力)が決定される。このため、対象排気系Eに存する無駄時間dの影響が補償(排除)され、これによっても、O2センサ6の出力VO2/OUTの目標値VO2/TARGETへの収束制御の安定性を高めることができる。
【0290】
さらに、本実施形態では、スライディングモード制御器27により制御すべき状態量Xとして、O2センサ6の偏差出力VO2 の現在以前の時系列データVO2(k),VO2(k-1)(より詳しくは、O2センサ6の推定偏差出力VO2バーの最新値以前の時系列データVO2(k+d)バー,VO2(k+d-1)バー)を用いることで、スライディングモード制御器27の演算処理を対象排気系Eの離散系モデル上で構築することができ、スライディングモード制御器27の演算処理を離散時間的なコンピュータ処理に適した簡素なものとすることができる。
【0291】
さらに、上記のようにスライディングモード制御器27の演算処理のために対象排気系Eの離散系モデルを用いることで、該離散系モデルの設定すべきパラメータとしての前記ゲイン係数a1,a2,b1を、本実施形態のような同定器25を用いてリアルタイムで同定して、離散系モデルの実際の対象排気系Eに対するモデル化誤差を実際の対象排気系Eの挙動状態に則して最小限に留めることができる。そして、該同定器25で同定したゲイン係数a1,a2,b1を用いてスライディングモード制御器27の演算処理を行って目標空燃比KCMD(対象排気系Eの目標入力)を決定することで、O2センサ6の出力VO2/OUTを目標値VO 2/TARGETに収束制御する上で、的確な目標空燃比KCMDを対象排気系Eの実際の挙動状態に則して決定することができ、ひいては、O2センサ6の出力VO2/OUTの目標値VO2/TARGETへの収束制御の精度を高めることができる。特に、スライディングモード制御器27により制御すべき状態量Xが前記超平面σ=0に収束していない段階でのモデル化誤差の影響が極力抑えられるため、O2センサ6の出力VO2/OUTの目標値VO2/TARGETへの収束制御の安定性を高めることができる。
【0292】
また、前記推定器26にあっても、対象排気系Eの離散系モデルを用いて演算処理を行うことで、その演算処理をコンピュータ処理に適した簡素なものとすることができる。そして、前記同定器25によりリアルタイムで同定された離散系モデルのゲイン係数a1,a2,b1を用いて推定器25の演算処理を行うことで、対象排気系Eの無駄時間d後のO2 センサ6の出力VO2/OUTの推定値を表す前記推定偏差出力VO2バーの精度を高めることができ、このような推定偏差出力VO2バーを用いてスライディングモード制御器27により目標空燃比KCMDを決定することで、無駄時間dの影響を確実に排除して、O2 センサ6の出力VO2/OUTの目標値VO2/TARGETへの収束制御を行うことができる。
【0293】
また、本実施形態では、同定器25による離散系モデルのゲイン係数a1,a2,b1の同定処理に際して、その処理に用いる前記同定偏差出力VO2ハット(これは対象排気系Eの離散系モデル上での出力に相当する)と、前記偏差出力VO2(これは対象排気系Eの実際の出力に相当する)とに、対象排気系Eの入力変化に対する出力変化のゲインが比較的大きなものとなる周波数帯(図16の低周波数帯C)に重みを有するフィリタリングを施すことによって、離散系モデルの周波数特性が実際の対象排気系Eの周波数特性に適合するように前記同定ゲイン係数a1ハット、a2ハット、b1ハットを算出することができる。そして、このような同定ゲイン係数a1ハット、a2ハット、b1ハットを用いてスライディングモード制御器27によって目標空燃比KCMDを決定すると共に推定器26による推定偏差出力VO2バーの算出処理を行うことで、O2センサ6の出力VO2/OUTを目標値VO2/TARGETに収束制御する上で重要な周波数帯における制御性を高めることができる。また、上記のフィルタリングの重み特性を適切に設定することで、O2センサ6の出力VO2/OUTの目標値VO2/TARGETへの収束制御の安定性や、速応性を高めることができる。
【0294】
また、本実施形態では、O2センサ6の偏差出力VO2により把握される対象排気系Eの出力の挙動状態が、触媒装置3を通過した排ガスの空燃比(これはO2センサ6の出力VO2/OUTに相当する)がリーン側からリッチ側に変化する挙動状態、すなわち、ゲイン係数a1,a2,b1の同定(更新)に好適な特定の挙動状態である場合において、離散系モデルのゲイン係数a1,a2,b1の同定処理を行うようにしているため、スライディングモード制御器27による目標空燃比KCMDの決定処理や推定器26による推定偏差出力VO2バーの算出処理を的確に行う上で好適な同定ゲイン係数a1ハット、a2ハット、b1ハットを算出することができ、ひいては、O2センサ6の出力VO2/OUTの目標値VO2/TARGETへの収束制御を確実に行うことができる。
【0295】
また、本実施形態では、同定器25による同定処理の安定性やスライディングモード制御器27による適応スライディングモード制御の安定性を判断し、それらが不安定と判断される場合に、同定器25の初期化を行うため、不適正な同定ゲイン係数a1ハット、a2ハット、b1ハットを用いて、スライディングモード制御器27により不適正な目標空燃比KCMDが決定されたり、推定器26により不適正な推定偏差出力VO2バーが算出されたりするような事態を回避することができる。
【0296】
また、本実施形態では、スライディングモード制御器27による適応スライディングモード制御が不安定であると判断された場合や、その判断後、前記タイマカウンタtmの初期値TM分の所定時間が経過するまでは、前記SLD操作入力uslを所定値として目標空燃比KCMDを決定するため、O2センサ6の出力VO2/OUTが異常な状態に制御されるような事態を確実に排除することができる。
【0297】
次に、本発明の第2の実施形態を図24乃至図26を参照して説明する。尚、本実施形態は、前述の第1の実施形態において前記同定器25による前記STEP5−11の処理、すなわち同定ゲイン係数ベクトルΘの評価処理(図11及び図18参照)の内容のみが、第1の実施形態のものと相違するものである。従って、本実施形態の説明では、第1の実施形態と同一構成部分については第1の実施形態と同一の参照符号を用いて説明を省略する。
【0298】
本実施形態は、前記同定器25が同定するパラメータであるゲイン係数a1,a2,b1の値、すなわち、前記同定ゲイン係数a1ハット、a2ハット、b1ハット(同定ゲイン係数ベクトルΘの要素)を所定の条件を満たすように制限し、その制限した同定ゲイン係数a1ハット、a2ハット、b1ハットの値を用いて、前記第1の実施形態で説明した如く推定器26及びスライディングモード制御器27による演算処理を行うものである。
【0299】
この場合、同定ゲイン係数a1ハット、a2ハット、b1ハットの値を制限するための前記所定の条件は、前記式(1)で表した離散系モデルの応答遅れ要素(より詳しくは1次目の自己回帰項及び2次目の自己回帰項)に係わる同定ゲイン係数a1ハット、a2ハットの値の組み合わせを所定の組み合わせに制限するための条件(以下、第1制限条件という)と、上記離散系モデルの無駄時間要素に係わる同定ゲイン係数b1ハットの値を制限するための条件(以下、第2制限条件という)とがある。
【0300】
ここで、これらの第1及び第2制限条件を説明する前に、同定ゲイン係数a1ハット、a2ハット、b1ハットの値を制限する理由を説明しておく。
【0301】
前述の第1の実施形態では、同定器25が同定ゲイン係数a1ハット、a2ハット、b1ハットを求めるに際して、それらの値に特別な制限を加えず、同定器25が不安定であると判断した場合に同定器25の初期化を行うようにしている。ところが、本願発明者等の知見によれば、O2センサ6の出力VO2/OUT(触媒装置3を通過した排ガスの酸素濃度)がその目標値VO2/TARGETに安定して制御されている状態で、スライディングモード制御器27により前述の如く求められる目標空燃比KCMD(空燃比の操作量)が平滑的な時間変化を呈する状況と、高周波振動的な時間変化を呈する状況とが生じることが判明した。この場合、いずれの状況においても、O2センサ6の出力VO2/OUTをその目標値VO2/TARGETに制御する上では支障がないものの、目標空燃比KCMDが高周波振動的な時間変化を呈する状況は、該目標空燃比KCMDに基づいて制御されるエンジン1の円滑な運転を行う上では、あまり好ましくない。
【0302】
そして、上記の現象について本願発明者等が検討したところ、スライディングモード制御器27が求める目標空燃比KCMDが平滑的なものとなるか高周波振動的なものとなるかは、同定器25により同定するゲイン係数a1,a2の値の組み合わせや、ゲイン係数b1の値の影響を受けることが判明した。
【0303】
そこで、本実施形態は、前記第1制限条件と第2制限条件とを適切に設定し、これらの条件により、同定ゲイン係数a1ハット、a2ハットの値の組み合わせや、同定ゲイン係数b1ハットの値を適切に制限することで、O2センサ6の出力VO2/OUTの目標値VO2/TARGETへの安定した制御性を確保しつつ、目標空燃比KCMDが高周波振動的なものとなるような状況を排除するものである。
【0304】
このような背景に鑑み、本実施形態では前記第1制限条件及び第2制限条件は次のように設定する。
【0305】
まず、同定ゲイン係数a1ハット、a2ハットの値の組み合わせを制限するための第1制限条件に関し、本願発明者等の検討によれば、目標空燃比KCMDが平滑的なものとなるか高周波振動的なものとなるかは、ゲイン係数a1,a2の値により定まる前記式(13)の係数値α1,α2、すなわち、前記推定器26が前記推定偏差出力VO2(k+d)バーを前記式(13)に従って求めるために使用する前記係数値α1,α2(これらの係数値α1,α2は前記式(11)で定義した行列Aの巾乗Adの第1行第1列成分及び第1行第2列成分である)の組み合わせが密接に関連している。
【0306】
具体的には、図24に示すように係数値α1,α2をそれぞれ成分とする座標平面を設定したとき、係数値α1,α2の組により定まる該座標平面上の点が図24の斜線を付した領域(三角形Q1Q2Q3で囲まれた領域(境界を含む)。以下この領域を推定係数安定領域という)に存するとき、目標空燃比KCMDの時間的変化が平滑的なものとなりやすい。逆に、係数値α1,α2の組により定まる点が上記の推定係数安定領域を逸脱しているような場合には、目標空燃比KCMDの時間的変化が高周波振動的なものとなったり、あるいは、O2センサ6の出力VO2/OUTの目標値VO2/TARGETへの制御性が悪化しやすい。
【0307】
従って、同定器25により同定するゲイン係数a1,a2の値、すなわち同定ゲイン係数a1ハット、a2ハットの値の組み合わせは、これらの値により定まる係数値α1,α2の組に対応する図24の座標平面上の点が上記推定係数安定領域内に存するように制限することが好ましい。
【0308】
尚、図24において、上記推定係数安定領域を含んで座標平面上に表した三角形領域Q1Q4Q3は、次式(51)により定義される系、すなわち、前記式(13)の右辺のVO2(k)及びVO2(k-1)をそれぞれVO2(k)バー及びVO2(k-1)バー(これらのVO2(k)バー及びVO2(k-1)バーは、それぞれ、推定器26により無駄時間d前に求められる推定偏差出力及びその1制御サイクル前に求められる推定偏差出力を意味する)により置き換えてなる式により定義される系が、理論上、安定となるような係数値α1,α2の組み合わせを規定する領域である。
【0309】
【数51】
【0310】
すなわち、式(51)により表される系が安定となる条件は、その系の極(これは、次式(52)により与えられる)が複素平面上の単位円内に存在することである。
【0311】
【数52】
【0312】
そして、図24の三角形領域Q1Q4Q3は、上記の条件を満たす係数値α1,α2の組み合わせを規定する領域である。従って、前記推定係数安定領域は、前記式(51)により表される系が安定となるような係数値α1,α2の組み合わせのうち、α1≧0となる組み合わせを表す領域である。
【0313】
一方、係数値α1,α2は、ゲイン係数a1,a2の値の組み合わせにより定まるので、逆算的に、係数値α1,α2の組み合わせからゲイン係数a1,a2の値の組み合わせも定まる。従って、係数値α1,α2の好ましい組み合わせを規定する図24の推定係数安定領域は、ゲイン係数a1,a2を座標成分とする図25の座標平面上に変換することができ、この変換を行うと、該推定係数安定領域は、図25の座標平面上では、例えば図25の仮想線で囲まれた領域(下部に凹凸を有する大略三角形状の領域。以下、同定係数安定領域という)に変換される。すなわち、ゲイン係数a1,a2の値の組により定まる図25の座標平面上の点が、同図の仮想線で囲まれた同定係数安定領域に存するとき、それらのゲイン係数a1,a2の値により定まる係数値α1,α2の組に対応する図24の座標平面上の点が前記推定係数安定領域内に存することとなる。
【0314】
従って、同定器25により求める同定ゲイン係数a1ハット、a2ハットの値を制限するための前記第1制限条件は、基本的には、それらの値により定まる図25の座標平面上の点が前記同定係数安定領域に存することとして設定することが好ましい。
【0315】
但し、図25に仮想線で示した同定係数安定領域の境界の一部(図の下部)は凹凸を有する複雑な形状を呈しているため、実用上、同定ゲイン係数a1ハット、a2ハットの値により定まる図25の座標平面上の点を同定係数安定領域内に制限するための処理が煩雑なものとなりやすい。
【0316】
そこで、本実施形態では、同定係数安定領域を、例えば図25の実線で囲まれた四角形Q5Q6Q7Q8の領域(境界を直線状に形成した領域。以下、同定係数制限領域という)により大略近似する。この場合、この同定係数制限領域は、図示の如く、|a1|+a2=1なる関数式により表される折れ線(線分Q5Q6及び線分Q5Q8を含む線)と、a1=A1L(A1L:定数)なる定値関数式により表される直線(線分Q6Q7を含む直線)と、a2=A2L(A2L:定数)なる定値関数式により表される直線(線分Q7Q8を含む直線)とにより囲まれた領域である。そして、同定ゲイン係数a1ハット、a2ハットの値を制限するための前記第1制限条件を、それらの値により定まる図25の座標平面上の点が上記同定係数制限領域に存することとして設定し、同定ゲイン係数a1ハット、a2ハットの値により定まる点が同定係数制限領域に存するようにそれらの値を制限する。この場合、同定係数制限領域の下辺部の一部は、前記同定係数安定領域を逸脱しているものの、現実には同定器25が求める同定ゲイン係数a1ハット、a2ハットの値により定まる点は上記の逸脱領域には入らないことを実験的に確認している。従って、上記の逸脱領域があっても、実用上は支障がない。
【0317】
尚、このような同定係数制限領域の設定の仕方は例示的なもので、該同定係数制限領域は、基本的には、前記同定係数安定領域に等しいか、もしくは該同定係数安定領域を大略近似し、あるいは、同定係数制限領域の大部分もしくは全部が同定係数安定領域に属するように設定すれば、どのような形状のものに設定してもよい。つまり、同定係数制限領域は、同定ゲイン係数a1ハット、a2ハットの値の制限処理の容易さ、実際上の制御性等を考慮して種々の設定が可能である。例えば本実施形態では、同定係数制限領域の上半部の境界を|a1|+a2=1なる関数式により規定しているが、この関数式を満たすゲイン係数a1,a2の値の組み合わせは、前記式(52)により与えられる系の極が複素平面上の単位円周上に存するような理論上の安定限界の組み合わせである。従って、同定係数制限領域の上半部の境界を例えば|a1|+a2=r(但し、rは上記の安定限界に対応する「1」よりも若干小さい値で、例えば0.99)なる関数式により規定し、制御の安定性をより高めるようにしてもよい。
【0318】
また、前記同定係数制限領域の基礎となる図25の同定係数安定領域も例示的なものであり、図24の推定係数安定領域に対応する同定係数安定領域は、係数値α1,α2の定義から明らかなように(式(12)、(13)を参照)、無駄時間d(より正確にはその設定値)の影響も受け、該無駄時間dの値によって、同定係数安定領域の形状が変化する。この場合、同定係数安定領域がどのような形状のものであっても、前記同定係数制限領域は、同定係数安定領域の形状に合わせて前述の如く設定すればよい。
【0319】
次に、同定器25が同定する前記ゲイン係数b1の値、すなわち同定ゲイン係数b1ハットの値を制限するための前記第2制限条件は本実施形態では次のように設定する。
【0320】
すなわち、本願発明者等の知見によれば、前記目標空燃比KCMDの時間的変化が高周波振動的なものとなる状況は、同定ゲイン係数b1ハットの値が過大もしくは過小となるような場合にも生じ易い。そこで、本実施形態では、同定ゲイン係数b1ハットの値の上限値B1H及び下限値B1L(B1H>B1L>0)をあらかじめ実験やシミュレーションを通じて定めておく。そして、前記第2制限条件を、同定ゲイン係数b1ハットの値が上限値B1H以下で且つ下限値B1L以上の値になること(B1L≦b1ハット≦B1Hの不等式を満たすこと)として設定する。
【0321】
以上説明した如く設定した第1制限条件及び第2制限条件により同定ゲイン係数a1ハット、a2ハット、b1ハットの値を制限するための処理は、本実施形態では、前記STEP5−11における同定ゲイン係数ベクトルΘの評価処理において、次のように行われる。
【0322】
すなわち、図26のフローチャートを参照して、同定器25は、前記図11のSTEP5−10で前述の如く求めた同定ゲイン係数a1(k)ハット、a2(k)ハット、b1(k)ハットについて、まず、同定ゲイン係数a1(k)ハット、a2(k)ハットの値の組み合わせを前記第1制限条件により制限するための処理をSTEP5−11−5〜5−11−1で行う。
【0323】
具体的には、同定器25は、まず、STEP5−10で求めた同定ゲイン係数a2(k)ハットの値が、前記同定係数制限領域におけるゲイン係数a2の下限値A2L(図25参照)以上の値であるか否かを判断する(STEP5−11−5)。
【0324】
このとき、a2(k)ハット<A2Lであれば、同定ゲイン係数a1(k)ハット、a2(k)ハットの値の組により定まる図25の座標平面上の点(以下、この点を(a1(k)ハット,a2(k)ハット)で表す)が同定係数制限領域から逸脱しているので、a2(k)ハットの値を強制的に上記下限値A2Lに変更する(STEP5−11−6)。この処理により、図25の座標平面上の点(a1(k)ハット,a2(k)ハット)は、少なくともa2=A2Lにより表される直線(線分Q7Q8を含む直線)の上側(該直線上を含む)の点に制限される。
【0325】
次いで、同定器25は、STEP5−10で求めた同定ゲイン係数a1(k)ハットの値が、前記同定係数制限領域におけるゲイン係数a1の下限値A1L(図25参照)以上の値であるか否か、並びに、同定係数制限領域におけるゲイン係数a1の上限値A1H(図25参照)以下の値であるか否かを順次判断する(STEP5−11−7、5−11−9)。尚、同定係数制限領域におけるゲイン係数a1の上限値A1Hは、図25から明らかなように折れ線|a1|+a2=1(但しa1>0)と、直線a2=A2Lとの交点Q8のa1座標成分であるので、A1H=1−A2Lである。
【0326】
このとき、a1(k)ハット<A1Lである場合、あるいは、a1(k)ハット>A1Hである場合には、図25の座標平面上の点(a1(k)ハット,a2(k)ハット)が同定係数制限領域から逸脱しているので、a1(k)ハットの値をそれぞれの場合に応じて、強制的に上記下限値A1Lあるいは上限値A1Hに変更する(STEP5−11−8、5−11−10)。
【0327】
この処理により、図25の座標平面上の点(a1(k)ハット,a2(k)ハット)は、a1=A1Lにより表される直線(線分Q6Q7を含む直線)と、a1=A1Hにより表される直線(点Q8を通ってa1軸に直行する直線)との間の領域(両直線上を含む)に制限される。
【0328】
尚、STEP5−11−7及び5−11−8の処理と、STEP5−11−9及び5−11−10の処理とは順番を入れ換えてもよい。また、前記STEP5−11−5及び5−11−6の処理は、STEP5−11−7〜5−11−10の処理の後に行うようにしてもよい。
【0329】
次いで、同定器25は、前記STEP5−11−5〜5−11−10を経た今現在のa1(k)ハット,a2(k)ハットの値が|a1|+a2≦1なる不等式を満たすか否か、すなわち、点(a1(k)ハット,a2(k)ハット)が|a1|+a2=1なる関数式により表される折れ線(線分Q5Q6及び線分Q5Q8を含む線)の下側(折れ線上を含む)にあるか上側にあるかを判断する(STEP5−11−11)。
【0330】
このとき、|a1|+a2≦1なる不等式が成立しておれば、前記STEP5−11−5〜5−11−10を経たa1(k)ハット,a2(k)ハットの値により定まる点(a1(k)ハット,a2(k)ハット)は、同定係数制限領域(その境界を含む)に存している。
【0331】
一方、|a1|+a2>1である場合は、点(a1(k)ハット,a2(k)ハット)が、同定係数制限領域からその上方側に逸脱している場合であり、この場合には、a2(k)ハットの値を強制的に、a1(k)ハットの値に応じた値(1−|a1(k)ハット|)に変更する(STEP5−11−12)。換言すれば、a1(k)ハットの値を現状に保持したまま、点(a1(k)ハット,a2(k)ハット)を|a1|+a2=1なる関数式により表される折れ線上(同定係数制限領域の境界である線分Q5Q6上、もしくは線分Q5Q8上)に移動させる。
【0332】
以上のようなSTEP5−11−5〜5−11−12の処理によって、同定ゲイン係数a1(k)ハット,a2(k)ハットの値は、それらの値により定まる点(a1(k)ハット,a2(k)ハット)が同定係数制限領域内に存するように制限される。尚、前記STEP5−10で求められた同定ゲイン係数a1(k)ハット,a2(k)ハットの値に対応する点(a1(k)ハット,a2(k)ハット)が同定係数制限領域内に存する場合は、それらの値は保持される。
【0333】
この場合、前述の処理によって、前記離散系モデルの1次目の自己回帰項に係わる同定ゲイン係数a1(k)ハットに関しては、その値が、同定係数制限領域における下限値A1L及び上限値A1Hの間の値となっている限り、その値が強制的に変更されることはない。また、a1(k)ハット<A1Lである場合、あるいは、a1(k)ハット>A1Hである場合には、それぞれ、同定ゲイン係数a1(k)ハットの値は、同定係数制限領域においてゲイン係数a1が採りうる最小値である下限値A1Lと、同定係数制限領域においてゲイン係数a1が採りうる最大値である下限値A1Hとに強制的に変更されるので、これらの場合における同定ゲイン係数a1(k)ハットの値の変更量は最小なものとなる。つまり、STEP5−10で求められた同定ゲイン係数a1(k)ハット,a2(k)ハットの値に対応する点(a1(k)ハット,a2(k)ハット)が同定係数制限領域から逸脱している場合には、同定ゲイン係数a1(k)ハットの値の強制的な変更は最小限に留められる。
【0334】
このようにして、同定ゲイン係数a1(k)ハット,a2(k)ハットの値を制限したのち、同定器25は、同定ゲイン係数b1(k)ハットの値を前記第2制限条件に従って制限する処理をSTEP5−11−13〜5−11−16で行う。
【0335】
すなわち、同定器25は、前記STEP5−10で求めた同定ゲイン係数b1(k)ハットの値が、前記下限値B1L以上であるか否かを判断し(STEP5−11−13)、B1L>b1(k)ハットである場合には、b1(k)ハットの値を強制的に上記下限値B1Lに変更する(STEP5−11−14)。
【0336】
さらに、同定器25は、同定ゲイン係数b1(k)ハットの値が、前記上限値B1H以上であるか否かを判断し(STEP5−11−15)、B1H<b1(k)ハットである場合には、b1(k)ハットの値を強制的に上記上限値B1Hに変更する(STEP5−11−16)。
【0337】
このようなSTEP5−11−13〜5−11−16の処理によって、同定ゲイン係数b1(k)ハットの値は、下限値B1L及び上限値B1Hの間の範囲の値に制限される。
【0338】
このようにして、同定ゲイン係数a1(k)ハット,a2(k)ハットの値の組み合わせと同定ゲイン係数b1(k)ハットの値とを制限した後には、同定器25は、前記図18のSTEP5−11−3及び5−11−4と同じ処理を行う。すなわち、前記STEP5−8(図11参照)で算出された同定誤差id/eの大きさが所定値ε0以下の十分小さなものとなったか否か(id/eがほぼ「0」に収束して、同定ゲイン係数a1ハット,a2ハット,b1ハットがほぼ確定した状態になったか否か)を判断する(STEP5−11−17)。このとき、|id/e|≦ε0であれば、そのまま図11のフローチャートの処理に復帰する。
【0339】
一方、STEP5−11−17の判断で、|id/e|>ε0である場合には、前述の第1の実施形態で説明した如くスライディングモード制御の安定性の判断の際に使用するタイマカウンタtm(カウントダウンタイマ)の値を所定の初期値TMにセットした後(タイマカウンタtmの起動。STEP5−11−18)、図11のフローチャートの処理に復帰する。
【0340】
以上説明した以外の他の作動は、前述の第1の実施形態と全く同一である。この場合において、図11のSTEP5−10で同定ゲイン係数a1(k)ハット,a2(k)ハット,b1(k)ハットを求めるために使用する同定ゲイン係数の前回値a1(k-1)ハット,a2(k-1)ハット,b1(k-1)ハットは、前回の制御サイクルにおけるSTEP5−11の処理で前述の如く第1及び第2制限条件により制限を行った同定ゲイン係数の値である。このため、各制御サイクルにおいて前記STEP5−10で求められる同定ゲイン係数a1(k)ハット,a2(k)ハット,b1(k)ハットの値は前記第1及び第2の制限条件を満たす値に収まりやすくなる。
【0341】
尚、本実施形態では、同定ゲイン係数b1の値は前記STEP5−11−13〜5−11−16の処理により、その上限及び下限が正の値に制限されるので、前記図21のSTEP8−3の処理を省略してもよい。
【0342】
以上説明した本実施形態によれば、第1の実施形態と同様の作用効果を奏することはもちろんであるが、同定器25により求める同定ゲイン係数a1ハット,a2ハット,b1ハットの値を前述の如く設定した第1及び第2制限条件に従って制限することで、特に、スイディングモード制御器27が求める目標空燃比KCMDが高周波振動的な変化を生じるのを確実に排除することができ、エンジン1の円滑な運転を行いつつ、触媒装置3の下流側のO2センサ6の出力VO2/OUT を目標値VO2/TARGETに高精度で制御することができる。すなわち、エンジン1の円滑な運転を行いつつ、触媒装置3の最適な浄化性能を確保することができる。
【0343】
この場合、特に対象排気系Eの離散系モデルの応答遅れ要素に係わる同定ゲイン係数a1ハット,a2ハットについては、それらの値を個別に制限するのではなく、それらの値を、両者の値の相関性をもった組み合わせにより制限することで、O2センサ6の出力VO2/OUTを目標値VO2/TARGETに制御し、また、目標空燃比KCMDが高周波振動的な変化を生じるのを確実に排除する上で最適な同定ゲイン係数a1ハット,a2ハットの値を得ることができる。
【0344】
また、同定ゲイン係数a1ハット,a2ハットの値の組み合わせの制限に際しては、式(1)の右辺の自己回帰項のうちの低次側の自己回帰項(1次目の自己回帰項)に係わる同定ゲイン係数a1ハット、別の言い方をすれば式(1)により表現した離散系モデルにおいて、O2センサ6のより新しい出力VO2/OUTもしくは偏差出力VO2 に係わる同定ゲイン係数a1ハットの値の変更量が最小となるようにa1ハット,a2ハットの値の組み合わせの制限を行うことで、より信頼性の高い目標空燃比KCMDを算出することができ、安定した制御を行うことができる。
【0345】
さらに、同定ゲイン係数a1ハット,a2ハットの値の組み合わせを制限するための前記同定係数制限領域(図25参照)は、その境界を直線状に設定したため、a1ハット,a2ハットの値を制限するための処理を容易に行うことができる。
【0346】
尚、本発明、特に本発明の内燃機関の排気系の空燃比制御装置は、前述した第1及び第2の実施形態に限定されるものではなく、例えば次のような変形態様も可能である。
【0347】
すなわち、前記第1及び第2の実施形態では、第2排ガスセンサとして、LAFセンサ(広域空燃比センサ)5を用いたが、第2排気ガスセンサは排ガスの空燃比を検出できるものであれば、通常のO2センサ等、他の形式のセンサを用いてもよい。
【0348】
また、前記第1及び第2の実施形態では、第1排ガスセンサとしてO2センサ6を用いたが、第1排ガスセンサは、制御すべき触媒装置下流の排ガスの特定成分の濃度を検出できるセンサであれば、他のセンサを用いてもよい。すなわち、例えば触媒装置下流の排ガス中の一酸化炭素(CO)を制御する場合はCOセンサ、窒素酸化物(NOX)を制御する場合にはNOXセンサ、炭化水素(HC)を制御する場合にはHCセンサを用いる。三元触媒装置を使用した場合には、上記のいずれのガス成分の濃度を検出するようにしても、触媒装置の浄化性能を最大限に発揮させるように制御することができる。また、還元触媒装置や酸化触媒装置を用いた場合には、浄化したいガス成分を直接検出することで、浄化性能の向上を図ることができる。
【0349】
また、第1及び第2の実施形態では、対象排気系Eの離散系モデルや、同定器25、推定器26、スライディングモード制御器27の演算処理において、LAFセンサ5の偏差出力kactやO2センサ6の偏差出力VO2を用いたが、LAFセンサ5の出力KACTやO2センサ6の出力VO2/OUTをそのまま用いて、対象排気系Eの離散系モデルを構築したり、同定器25、推定器26、スライディングモード制御器27の演算処理を行うようにしてもよい。但し、離散系モデルの簡素化や同定器25、推定器26、スライディングモード制御器27の演算処理の簡素化を図る上では、本実施形態のように偏差出力kact,VO2を用いることが好ましい。また、この場合において、偏差出力kact(=KACT−FLAF/BASE)に係わる前記基準値FLAF/BASEは必ずしも一定値とする必要はなく、該基準値FLAF/BASEをエンジン1の回転数NEや吸気圧PB等に応じて設定するようにしてもよい。
【0350】
また、第1及び第2の実施形態では、空燃比操作量決定部13により決定する操作量を触媒装置3に進入する排ガスの目標空燃比KCMD(対象排気系Eの目標入力)とし、その目標空燃比KCMDに従ってエンジン1の燃料供給量をフィードバック制御するようにしたが、例えばエンジン1の燃料供給量の補正量を空燃比操作量決定部13により決定するようにすることも可能であり、また、目標空燃比KCMDからフィードフォワード的にエンジン1の燃料供給量を制御するようにすることも可能である。
【0351】
また、第1及び第2の実施形態では、スライディングモード制御器27は、外乱の影響を考慮した適応則を有する適応スライディングモード制御を用いたが、該適応則を省略した一般のスライディングモード制御を用いるようにしてもよい。
【0352】
また、第1及び第2の実施形態では、スライディングモード制御器27は、制御すべき状態量を二つの偏差出力VO2(k),VO2(k-1)としたが、さらに多くの偏差出力(例えばVO2(k),VO2(k-1),VO2(k-2)等)を制御すべき状態量として用いるようにしてもよい。
【0353】
また、第1及び第2の実施形態では、適応スライディングモード制御が不安定であると判断した場合に、前記SLD操作入力uslを強制的に所定値とし、従って、目標空燃比KCMDも所定値となるようにしたが、適応スライディングモード制御が不安定であると判断した場合に、例えばPID制御器等、空燃比操作量決定部13とは別に備えた制御器を用いて、暫定的にO2センサ6の出力VO2/OUTが目標値VO2/TARGETに収束するように目標空燃比KCMDを決定するようにしてもよい。
【0354】
また、第1及び第2の実施形態では、推定器26による演算処理とスライディングモード制御器27の演算処理とを、前記式(1)により表される対象排気系Eの同一の離散系モデルに基づいて行うようにしたが、推定器26とスライディングモード制御器27とで各別の離散系モデルに基づいて演算処理を行うようにしてもよい。さらにこの場合、推定器26用の離散系モデルのパラメータは、あらかじめ定めた所定値に保持したり、エンジン1の運転状態や触媒装置3の劣化状態に応じてマップ等を用いて適宜設定するようにしてもよく、さらには、推定器26による演算処理は、対象排気系Eの連続系モデルに基づいて行うようにしてもよい。
【0355】
さらに、排気系の無駄時間が十分に小さいような場合にあっては、推定器26を省略するようにしてもよく、この場合には、例えば前記第1の実施形態において、無駄時間d=0として、推定器26の処理を省略すればよい。この場合、スライディングモード制御器は、前記式(20)、(21)、(25)において、d=0とした式によって、目標空燃比KCMDを決定するための等価制御入力ueq、到達則入力urch及び適応則入力uadpを求めるようにすればよい。また、この場合において、第2の実施形態のように同定器25により同定するパラメータの値を制限する場合には、その制限条件は、推定器26の処理と無関係に、制御の安定性等を考慮し、各種実験やシミュレーションを通じて設定すればよい。
【0356】
また、第1及び第2の実施形態では、スライディングモード制御器27により同定器25で同定された離散系モデルのゲイン係数a1,a2,b1を用いて目標空燃比KCMDを決定するようにしたが、適応制御器等の他の漸化式形式の制御器により同定されたゲイン係数a1,a2,b1を用いて目標空燃比KCMDを決定するようにしてもよく、さらには、同定されたゲイン係数a1,a2,b1を用いて目標空燃比KCMDを決定し得るものであれば、ファジー制御器やニューラルネットワーク型の制御器を用いてよい。
【0357】
また、第1及び第2の実施形態では、対象排気系Eの無駄時間dをあらかじめ定めた値に設定したが、ゲイン係数a1,a2,b1と共に該無駄時間dを同定するようにすることも可能である。そして、この場合において、同定する無駄時間dの値を前記第2の実施形態と同様に適当な条件によって制限するようにしてもよい。
【0358】
また、第1及び第2の実施形態では、同定器25によるゲイン係数a1,a2,b1の同定を触媒装置3を通過した排ガスの空燃比がリーン側からリッチ側に変化する挙動状態において行うようにしたが、該空燃比がリッチ側からリーン側に変化する挙動状態において同定を行うようにしてもよく、さらには、該挙動状態を区別することなく任意の挙動状態において、逐次ゲイン係数a1,a2,b1の同定を行ったり、あるいは排ガスの空燃比がリーン側からリッチ側に変化する挙動状態と排ガスの空燃比がリッチ側からリーン側に変化する挙動状態とで各別にゲイン係数a1,a2,b1の同定を行うようにしてもよい。このような場合には、推定器26とスライディングモード制御器27とで用いるゲイン係数a1,a2,b1は、同定器25で実際に今回の制御サイクルで同定したゲイン係数を用いることとなる。
【図面の簡単な説明】
【図1】本発明の内燃機関の空燃比制御装置の第1の実施形態の全体的システム構成図。
【図2】図1の装置で使用するO2センサの出力特性図。
【図3】図1の装置の空燃比操作量決定部の基本構成を示すブロック図。
【図4】図1の装置で用いるスライディングモード制御を説明するための説明図。
【図5】図1の装置の適応制御器の基本構成を示すブロック図。
【図6】図1の装置のエンジンの燃料制御に係わる処理を説明するためのフローチャート。
【図7】図6のフローチャートにおけるサブルーチン処理を説明するためのフローチャート。
【図8】図1の装置の空燃比操作量決定部の全体的処理を説明するためのフローチャート。
【図9】図8のフローチャートのサブルーチン処理を説明するためのフローチャート。
【図10】図8のフローチャートのサブルーチン処理を説明するためのフローチャート。
【図11】図8のフローチャートのサブルーチン処理を説明するためのフローチャート。
【図12】図11のフローチャートのサブルーチン処理を説明するための説明図。
【図13】図11のフローチャートのサブルーチン処理を説明するための説明図。
【図14】図11のフローチャートのサブルーチン処理を説明するためのフローチャート。
【図15】図11のフローチャートのサブルーチン処理を説明するための説明図。
【図16】図11のフローチャートのサブルーチン処理を説明するための説明図。
【図17】図11のフローチャートのサブルーチン処理を説明するための説明図。
【図18】図11のフローチャートのサブルーチン処理を説明するためのフローチャート。
【図19】図8のフローチャートのサブルーチン処理を説明するためのフローチャート。
【図20】図8のフローチャートのサブルーチン処理を説明するためのフローチャート。
【図21】図8のフローチャートのサブルーチン処理を説明するためのフローチャート。
【図22】図21のフローチャートのサブルーチン処理を説明するためのフローチャート。
【図23】図8のフローチャートのサブルーチン処理を説明するためのフローチャート。
【図24】本発明の内燃機関の空燃比制御装置の第2の実施形態を説明するための説明図。
【図25】
本発明の内燃機関の空燃比制御装置の第2の実施形態を説明するための説明図
。
【図26】
本発明の内燃機関の空燃比制御装置の第2の実施形態の作動を説明するための
フローチャート。
【符号の説明】
1…エンジン(内燃機関)、E…対象排気系(制御対象の系)、6…O2センサ(排ガスセンサ)、13…空燃比操作量決定部、25…同定器、26…推定器、27…スライディングモード制御器。[0001]
BACKGROUND OF THE INVENTION
The present invention relates to an air-fuel ratio control apparatus for an internal combustion engine.
[0002]
[Prior art]
For example, in a system in which exhaust gas from an internal combustion engine such as an automobile is purified and released by a catalyst device such as a three-way catalyst, the air-fuel ratio of the exhaust gas entering the catalyst device is controlled so that the exhaust gas purification ability of the catalyst device is exhibited well. This is desired from the viewpoint of environmental protection.
[0003]
In order to perform such control, conventionally, as seen in, for example, Japanese Patent Laid-Open No. 5-321721, an exhaust gas sensor (oxygen concentration sensor) in which the oxygen concentration of exhaust gas that has passed through the catalyst device is disposed downstream of the catalyst device. ), The target air-fuel ratio of the exhaust gas upstream of the catalyst device is determined using PID control so that the detected value becomes a predetermined appropriate value, and the internal combustion engine is controlled according to the target air-fuel ratio. The applicant of the present invention has proposed that the air-fuel ratio of exhaust gas entering the catalyst device is set within a predetermined window that can exhibit good purification performance of the catalyst device.
[0004]
On the other hand, by further study by the inventors of the present application, in order to maximize the purification capacity of the catalytic device as much as possible without depending on the operating state of the internal combustion engine, deterioration with time of the catalytic device, etc., the catalytic device is It has been found that it is necessary to stabilize the concentration of a specific component such as oxygen concentration in the exhaust gas that has passed through to a predetermined appropriate value (a constant value) with high accuracy and stability. As described above, in the conventional control method using PID control, the oxygen concentration or the like in the exhaust gas that has passed through the catalyst device as described above is affected by disturbance, dead time existing in the exhaust system including the catalyst device, or the like. It has been found that it is difficult to stably set the predetermined value (a constant value) with high accuracy.
[0005]
For this reason, the inventors of the present application modeled the exhaust system from the upstream side to the downstream side of the catalyst device in a continuous system (specifically, a continuous time system), and based on the model, the oxygen concentration in the exhaust gas that passed through the catalyst device Have previously devised a system for controlling the air-fuel ratio of the exhaust gas entering the catalyst device so that the value of the above becomes a predetermined appropriate value (see, for example, Japanese Patent Application Nos. 9-67591 and 8-84048).
[0006]
The above-described modeling compensates (eliminates) the effects of dead time existing in the exhaust system including the catalyst device, and provides a control method (for example, sliding mode control) that has higher control stability against disturbance than PID control. ) Can be constructed, and as a result, the accuracy and stability of the air-fuel ratio control of the exhaust system can be improved.
[0007]
In such an air-fuel ratio control device, it is desired to increase the stability of control regardless of the operating state of the internal combustion engine.
[0008]
[Problems to be solved by the invention]
In view of such a background, the present invention provides an internal combustion engine capable of stably and highly accurately controlling the air-fuel ratio of exhaust gas so that the concentration of a specific component in the exhaust gas becomes a predetermined target value. An object is to provide an air-fuel ratio control device.
[0009]
[Means for Solving the Problems]
In order to achieve such an object, the first aspect of the air-fuel ratio control apparatus for an internal combustion engine of the present invention is an exhaust gas provided in an exhaust system of the internal combustion engine to detect the concentration of a specific component in the exhaust gas of the internal combustion engine. In an air-fuel ratio control apparatus for an internal combustion engine that includes a sensor and controls the air-fuel ratio of the exhaust gas of the internal combustion engine so that the output of the exhaust gas sensor becomes a predetermined target value, an operation amount that defines the air-fuel ratio of the exhaust gas of the internal combustion engine An air-fuel ratio manipulated variable determining means that uses a sliding mode control, and an identifying means for identifying a parameter to be set in a model of a system that is controlled by the air-fuel ratio manipulated variable determiner, Quantity determination means , As the parameter value of the model, The identification value of the parameter identified by the identification means and the predetermined value are selectively used according to the operating state of the internal combustion engine. , By a predetermined arithmetic expression of the sliding mode control including the coefficient defined by the value of the parameter of the model The operation amount is determined.
[0010]
According to the first aspect of the present invention, the air-fuel ratio manipulated variable determining means determines the manipulated variable by using a parameter value of the model by sliding mode control. In this case, as the value of this parameter, an identification value obtained by the identification means and a predetermined value different from the identification value are selectively used according to the operating state of the internal combustion engine. As a result, it is possible to stably control the air-fuel ratio of the exhaust gas so that the output of the exhaust gas sensor becomes a predetermined target value regardless of the operating state of the internal combustion engine.
[0011]
A second aspect of the air-fuel ratio control apparatus for an internal combustion engine of the present invention includes an exhaust gas sensor provided in an exhaust system of the internal combustion engine for detecting the concentration of a specific component in the exhaust gas of the internal combustion engine, and the exhaust gas In an air-fuel ratio control apparatus for an internal combustion engine that controls the air-fuel ratio of the exhaust gas of the internal combustion engine so that the output of the sensor becomes a predetermined target value, the air-fuel ratio operation for determining an operation amount that defines the air-fuel ratio of the exhaust gas of the internal combustion engine The amount determining means and the air-fuel ratio manipulated variable determining means generate data indicating the estimated value of the exhaust gas sensor output after the dead time of the system that is the control target and that generates the exhaust gas sensor output. An estimation unit; and an identification unit that identifies a parameter to be set in a model of the system to be controlled; , As the parameter value of the model, The identification value of the parameter identified by the identification means and the predetermined value are selectively used according to the operating state of the internal combustion engine. , By a predetermined calculation formula for estimation including a coefficient defined by the parameter value of the model Data indicating an estimated value of the output of the exhaust gas sensor is generated, and the air-fuel ratio manipulated variable determining means determines the manipulated variable so that the estimated value indicated by the data generated by the estimating means becomes the target value. It is characterized by that.
[0012]
According to the second aspect of the present invention, the estimation means generates data representing an estimated value of the output of the exhaust gas sensor after the dead time using the parameter value of the model. In this case, as the value of this parameter, an identification value obtained by the identification means and a predetermined value different from the identification value are selectively used according to the operating state of the internal combustion engine. The air-fuel ratio manipulated variable determining means determines the manipulated variable so that the estimated value of the output of the exhaust gas sensor indicated by the data generated by the estimating means becomes the target value. As a result, it is possible to stably control the air-fuel ratio of the exhaust gas so that the output of the exhaust gas sensor becomes a predetermined target value regardless of the operating state of the internal combustion engine.
[0013]
In the first and second aspects of the present invention, the operating state of the internal combustion engine includes whether the throttle valve of the internal combustion engine is fully open, whether fuel supply to the internal combustion engine is stopped, or not. In this state, it is preferable that the engine is idling.
[0014]
DETAILED DESCRIPTION OF THE INVENTION
A first embodiment of the present invention will be described with reference to FIGS.
[0015]
FIG. 1 is a block diagram showing the overall configuration of the control device of this embodiment. In FIG. 1,
[0016]
The control device of the present embodiment controls the air-fuel ratio of the exhaust system of the
[0017]
The wide area air-
[0018]
Further, O downstream of the
[0019]
The
[0020]
The basic fuel injection
[0021]
Further, the first correction coefficient KTOTAL obtained by the first correction
[0022]
The second correction coefficient KCMDM obtained by the second correction
[0023]
The basic fuel injection amount Tim is corrected by the first correction coefficient KTOTAL and the second correction coefficient KCMDM by multiplying the basic fuel injection amount Tim by the first correction coefficient KTOTAL and the second correction coefficient KCMDM. By this correction, the required fuel injection amount Tcyl of the
[0024]
Since the applicant of the present invention has disclosed a more specific calculation method of the basic fuel injection amount Tim, the first correction coefficient KTOTAL, and the second correction coefficient KCMDM in Japanese Patent Laid-Open No. 5-79374, etc. Then, detailed description is abbreviate | omitted.
[0025]
In addition to the functional configuration described above, the
[0026]
Although the details of the air-fuel ratio manipulated
[0027]
Further, in this embodiment, the
[0028]
The global feedback control unit 15 corrects the required fuel injection amount Tcyl (multiplies the required fuel injection amount Tcyl) so that the output KACT of the
[0029]
In this embodiment, the feedback manipulated variable KLAF generated by the
[0030]
Then, the global feedback control unit 15 generates a feedback operation amount KLAF generated by the
[0031]
The local feedback control unit 16 estimates an actual air-fuel ratio # nA / F (n = 1, 2, 3, 4) for each cylinder from the output KACT of the
[0032]
Here, in brief, the
Such an
[0033]
Further, each
The local feedback control unit 16 multiplies the value obtained by multiplying the required fuel injection amount Tcyl by the feedback correction coefficient KFB of the global feedback control unit 15 by the feedback correction coefficient #nKLAF for each cylinder. Then, output fuel injection amount #nTout (n = 1, 2, 3, 4) of each cylinder is obtained.
The output fuel injection amount #nTout of each cylinder obtained in this way is corrected for each cylinder by the
[0034]
The adhesion correction is disclosed in detail in, for example, JP-A-8-21273 by the applicant of the present application, and further description thereof is omitted here. In FIG. 1, the sensor output selection processing unit denoted by
[0035]
Next, the air-fuel ratio manipulated
[0036]
As described above, the air-fuel ratio manipulated
[0037]
In this case, in this embodiment, in order to simplify the processing by the air-fuel ratio manipulated
[0038]
[Expression 1]
[0039]
This equation (1) indicates that the target exhaust system E is detected from the deviation output kact of the
[0040]
The air-fuel ratio manipulated
[0041]
That is, the air-fuel ratio manipulated
[0042]
The algorithm of the arithmetic processing by the
[0043]
First, with respect to the
[0044]
The
[0045]
That is, for each predetermined control cycle, the
[0046]
[Expression 2]
[0047]
This equation (2) is obtained by shifting equation (1) to the past by one control cycle, and changing the gain coefficients a1, a2, and b1 to the identification gain coefficients a1 hat (k-1), a2 hat (k-1). ), B1 hat (k-1). Further, “d” used in the third term of the expression (2) is a set value of the dead time of the target exhaust system E (more specifically, the set value of the dead time is expressed by the number of control cycles), and the setting is made. The value is set to be equal to or slightly longer than the actual dead time of the target exhaust system E.
[0048]
Here, when the vectors Θ and ξ defined by the following equations (3) and (4) are introduced (subscript “T” in equations (3) and (4) means transposition, the same applies hereinafter).
[0049]
[Equation 3]
[0050]
[Expression 4]
[0051]
The formula (2) is represented by the following formula (5).
[0052]
[Equation 5]
[0053]
Further, the
[0054]
[Formula 6]
[0055]
Then, the
[0056]
[Expression 7]
[0057]
Here, “Kθ” in equation (7) is a third-order vector determined by the following equation (8) (the degree of change corresponding to the identification error id / e of each identification gain coefficient a1 hat, a2 hat, b1 hat). (Gain coefficient vector to be defined).
[0058]
[Equation 8]
[0059]
Further, “P” in the above equation (8) is a cubic square matrix determined by the recurrence equation of the following equation (9).
[0060]
[Equation 9]
[0061]
In the equation (9), “λ 1 ”,“ Λ 2 Is 0 <λ 1 ≦ 1 and 0 ≦ λ 2 The initial value P (0) of “P” is a diagonal matrix in which each diagonal component is a positive number.
[0062]
In this case, “λ in equation (9) 1 ”,“ Λ 2 Depending on the setting method, various specific algorithms such as a fixed gain method, a gradually decreasing gain method, a weighted least square method, a least square method, a fixed trace method, and the like are configured. In the present embodiment, for example, the least square method ( In this case, λ 1 = Λ 2 = 1) is adopted.
[0063]
In the present embodiment, the
[0064]
Note that the
[0065]
Next, in order to compensate for the influence of the dead time d of the target exhaust system E when the target air-fuel ratio KCMD is determined by the sliding
[0066]
First, when a vector X defined by the following equation (10) is introduced into the discrete system model represented by the equation (1),
[0067]
[Expression 10]
[0068]
Expression (1) can be rewritten as the following expression (11).
[0069]
## EQU11 ##
[0070]
Here, when the recurrence formula of the equation (11) is repeatedly used, X (k + d) after the dead time d is the deviation output kact of the matrix A and the vector B defined in the equation (11) and the
[0071]
[Expression 12]
[0072]
In this case, the first row component on the left side of Equation (12) is O after the dead time d. 2 Since it is the deviation output VO2 (k + d) of the
[0073]
Therefore, focusing on the first row components on both sides of Equation (12), the matrix A of the first term on the right side d The first row, first column component and the first row, second column component of .alpha. j-1 If the first row component of B (j = 1, 2, d) is set as βj (j = 1, 2, d), O 2 The estimated deviation output VO2 (k + d) bar of the
[0074]
[Formula 13]
[0075]
Also, the coefficient values α1, α2 and βj (j = 1, 2, d) in the equation (13) are gain coefficients a1, a2, b1 constituting the components of the matrix A and the vector B (see equation (11)). The identification gain coefficients a1 hat, a2 hat, b1 hat obtained by the
[0076]
Therefore, the
[0077]
Note that the
[0078]
Next, the sliding
[0079]
First, general sliding mode control will be briefly described with reference to FIG.
[0080]
The sliding mode control is a variable structure type feedback control technique. In this control technique, for example, the state quantity of the controlled object to be controlled is represented by x. 1 , X 2 These state quantities x 1 , X 2 Is a linear function σ = s 1 x 1 + S 2 x 2 (S 1 , S 2 Is used in advance to design a hyperplane represented by σ = 0. This hyperplane σ = 0 is often called a switching line when the phase space is a quadratic system (when there are two state quantities), and the linear function σ is called a switching function. When the order of the phase space is further increased, the switching line becomes a switching surface, and further, a hyperplane that cannot be geometrically illustrated. The hyperplane is sometimes called a slip plane. In the present specification, these are represented as a linear function and a hyperplane.
[0081]
Then, this sliding mode control is performed, for example, as shown by a point P in FIG. 1 , X 2 When σ ≠ 0, the state quantity x follows the so-called reaching law 1 , X 2 Is converged at high speed on the hyperplane σ = 0 by high gain control (mode 1), and the state quantity x is further controlled by so-called equivalent control input. 1 , X 2 On the hyperplane σ = 0 and the equilibrium point (x 1 = X 2 (Point of 0)) (mode 2).
[0082]
In such sliding mode control, the state quantity x 1 , X 2 As long as it converges on the hyperplane σ = 0, the state quantity x is extremely stable without being affected by disturbance or the like by the equivalent control input. 1 , X 2 Is constrained on the hyperplane σ = 0 and converged to the equilibrium point of the hyperplane σ = 0. If there is a disturbance or a modeling error of the controlled object, the state quantity x 1 , X 2 Strictly speaking, the above equilibrium point (x 1 = X 2 = 0), it converges in the vicinity of the equilibrium point.
[0083]
In such sliding mode control, in particular, in the
[0084]
The sliding
[0085]
First, the construction of the hyperplane necessary for the adaptive sliding mode control of the sliding
[0086]
The sliding
[0087]
Therefore, as a basic idea of the sliding mode control in the present embodiment, as the state quantity to be controlled, for example, the O obtained in each control cycle. 2 Using the deviation output VO2 (k) of the
[0088]
[Expression 14]
[0089]
In the present embodiment, the time series data of the estimated deviation output VO2 bar obtained by the
[0090]
When the linear function σ is defined as described above, the hyperplane for sliding mode control is represented by σ = 0 (in this case, since there are two state quantities, the hyperplane is a straight line, see FIG. 4). The coefficient s of the linear function σ defining this hyperplane σ = 0 1 , S 2 In the present embodiment, (see formula (14)) is set in advance as follows.
[0091]
That is, in a state where the vector X (hereinafter simply referred to as the state quantity X) in the equation (14) having the state quantities VO2 (k) and VO2 (k-1) as components converges on the hyperplane σ = 0, Since the value of the linear function σ is “0”, the following equation (15) is obtained from this and equation (14).
[0092]
[Expression 15]
[0093]
Here, since the system represented by Expression (15) is a first-order lag system with no input, the equilibrium point (VO2 (k) = VO2 (k-1) = 0 where the state quantity X is hyperplane σ = 0. The condition for stable convergence to the point is that the pole of the system represented by the equation (15) (in this case, this pole is “−s 2 / S 1 Is present in the unit circle.
[0094]
Therefore, in this embodiment, the coefficient s of the linear function σ 1 , S 2 Is set so as to satisfy the condition of the following equation (16).
[0095]
[Expression 16]
[0096]
In the present embodiment, the coefficient s is used for simplification. 1 = 1 (in this case s 2 / S 1 = S 2 ) <-1 <s 2 <Coefficient s to satisfy condition of 1 2 Set the value of.
[0097]
On the other hand, the SLD operation input usl (LAF sensor 5) to be generated by the sliding
[0098]
[Expression 17]
[0099]
In the present embodiment, the equivalent control input ueq, the reaching law input urch, and the adaptive law input uadp are based on the discrete system model represented by the above formula (1) or (11) as follows. Can be sought.
[0100]
First, regarding the equivalent control input ueq, the condition for the state quantity X to remain on the hyperplane σ = 0 is σ (k + 1) = σ (k) = 0. Using equation (14), the following equation (18) can be rewritten.
[0101]
[Formula 18]
[0102]
Here, the equivalent control input ueq is an input to be given to the target exhaust system E in order to constrain the state quantity X to the hyperplane σ = 0 (the exhaust gas air-fuel ratio detected by the
[0103]
Therefore, the following equation (19) is obtained from the equation (18):
[0104]
[Equation 19]
[0105]
Furthermore, the following equation (20) is obtained by shifting both sides of the equation (19) by the dead time d.
[0106]
[Expression 20]
[0107]
This equation (20) is a basic equation for obtaining the equivalent control input ueq for each control cycle in this embodiment.
[0108]
Next, in the present embodiment, the reaching law input urch is basically determined by the following equation (21).
[0109]
[Expression 21]
[0110]
That is, the reaching law input urch is determined so as to be proportional to the value σ (k + d) of the linear function σ after the dead time d in consideration of the influence of the dead time d.
[0111]
In this case, the coefficient F in equation (21) (which defines the reaching law gain) is set as follows.
[0112]
That is, in the above equation (11), when kact (k) = ueq (k) + urch (k) and further using equations (14), (20), and (21), the following equation (22) is obtained.
[0113]
[Expression 22]
[0114]
Here, since the system represented by Expression (22) is a first-order lag system with no input, the value of the linear function σ converges stably on the hyperplane σ = 0 (the state quantity X is hyperplane σ = 0). The condition is that the pole of the system represented by the equation (22) (in this case, this pole is “1-F”) is present in the unit circle.
[0115]
Therefore, in this embodiment, the coefficient F that defines the reaching law input urch is set so as to satisfy the condition of the following equation (23).
[0116]
[Expression 23]
[0117]
Regarding the behavior of the value of the linear function σ, there is a possibility that the value of the linear function σ may cause a vibrational change (so-called chattering) with respect to the hyperplane σ = 0. The pole “1-F” of the system represented by the formula (22) preferably satisfies the condition of 0 <1-F <1 in addition to the above condition.
[0118]
Therefore, the coefficient F that defines the reaching law input urch is more preferably set so as to satisfy the condition of the following equation (24).
[0119]
[Expression 24]
[0120]
Next, in this embodiment, the adaptive law input uad p is basically determined by the following equation (25) (ΔT in equation (25) is the cycle of the control cycle).
[0121]
[Expression 25]
[0122]
That is, the adaptive law input uadp considers the effect of the dead time d, and the integrated value for each control cycle of the value of the linear function σ until after the dead time d (this corresponds to the integrated value of the value of the linear function σ). To be proportional to
[0123]
In this case, the coefficient G in equation (25) (which defines the gain of the adaptive law) is set as follows.
[0124]
First, assuming that an error component due to the influence of disturbance or the like between the air-fuel ratio detected by the
[0125]
[Equation 26]
[0126]
Then, when this equation (26) is applied to the above equation (11) and further equations (14), (20), (21), and (25) are used, the following equation (27) is obtained.
[0127]
[Expression 27]
[0128]
Here, when both sides of the equation (27) are Z-transformed, the following equation (28) is obtained.
[0129]
[Expression 28]
[0130]
Further, when this equation (28) is modified and arranged, the following equation (29) is obtained.
[0131]
[Expression 29]
[0132]
In the equations (28) and (29), “Σ” and “V” are obtained by Z-transforming the linear function σ and the error component v, respectively. Further, M (Z) in the equation (29) is a pulse transfer function of the linear function σ with respect to the error component v, and is represented by the upper partial expression of the equation (29).
[0133]
In this case, the condition for the linear function σ to be stable with respect to the error component v (disturbance) is the pole of the pulse transfer function M (Z), that is, the solution of the characteristic equation represented by the following equation (30) (this There are two solutions) in the unit circle,
[0134]
[30]
[0135]
The two solutions of the characteristic equation (30) are given by the following equations (31) and (32), where λm1 and λm2 are the solutions.
[0136]
[31]
[0137]
[Expression 32]
[0138]
Accordingly, the condition that the linear function σ is stable with respect to the error component v (disturbance) is that λm1 and λm2 given by the above equations (31) and (32) exist in the unit circle.
[0139]
Therefore, in this embodiment, in order to satisfy this condition, the coefficient G is set by the following equation (33).
[0140]
[Expression 33]
[0141]
The sliding
[0142]
Therefore, in the present embodiment, the sliding
[0143]
[Expression 34]
[0144]
In the present embodiment, actually, the time series data of the estimated deviation output VO2 bar sequentially obtained by the
[0145]
[Expression 35]
[0146]
Then, the sliding
[0147]
[Expression 36]
[0148]
Similarly, the sliding
[0149]
[Expression 37]
[0150]
The gain coefficients a1, a2, and b1 required for calculating the equivalent control input ueq, the reaching law input urch, and the adaptive law input uadp by the equations (34), (36), and (37) are as follows: In the embodiment, the latest identification gain coefficients a1 (k) hat, a2 (k) hat, and b1 (k) hat obtained by the
[0151]
The sliding
[0152]
This is a basic algorithm for determining the SLD operation input usl to be given to the target exhaust system E for each control cycle by the sliding
[0153]
By the way, although the sliding
[0154]
[Formula 38]
[0155]
The above is the basic algorithm for determining the target air-fuel ratio KCMD by the sliding
[0156]
In the present embodiment, the stability of the adaptive sliding mode control process by the sliding
[0157]
Next, the global feedback control unit 15, particularly the
[0158]
Referring to FIG. 1, the global feedback control unit 15 performs feedback control so as to converge the output KACT (detected air-fuel ratio) of the
[0159]
The
[0160]
Here, the
[0161]
[39]
[0162]
[Formula 40]
[0163]
[Expression 41]
[0164]
[Expression 42]
[0165]
Here, the adaptive parameter θ hat shown in the equation (41) is a scalar quantity element b that determines the gain of the
[0166]
[Equation 43]
[0167]
(44)
[0168]
[Equation 45]
[0169]
The
[0170]
In this case, specifically, the adaptive parameter θ hat is calculated by the following equation (46).
[0171]
[Equation 46]
[0172]
In the equation (46), Γ (j) is a gain matrix that determines the setting speed of the adaptive parameter θ hat (the order of this matrix is m + n + d p ), E asterisk (j) indicates an estimation error of the adaptive parameter θ hat, and is represented by a recurrence formula as shown in equations (47) and (48), respectively.
[0173]
[Equation 47]
[0174]
[Formula 48]
[0175]
Here, “D (Z -1 ) "Is an asymptotically stable polynomial for adjusting the convergence. In this embodiment, D (Z -1 ) = 1.
[0176]
Note that λ in equation (47) 1 (j), λ 2 Depending on how (j) is selected, various specific algorithms such as a gradual decrease gain algorithm, a variable gain algorithm, a fixed trace algorithm, and a fixed gain algorithm can be obtained. In a time-varying plant such as fuel injection or air-fuel ratio of the
[0177]
As described above, the adaptive parameter θ hat (s 0 , R 1 , R 2 , R Three , B 0 ) And the target air-fuel ratio KCMD determined by the air-fuel ratio manipulated
[0178]
[Formula 49]
[0179]
The feedback manipulated variable KSTR obtained by the equation (49) becomes “target air-fuel ratio KCMD” when the output KACT of the
[0180]
As is apparent from the above, the
[0181]
Note that this type of recurrence type controller may be constructed using a so-called optimum regulator. In this case, however, the parameter adjusting mechanism is generally not provided, and the dynamics of the
[0182]
The above is the details of the
[0183]
The
[0184]
In addition, the switching
[0185]
Such an operation of the switching
[0186]
Next, details of the operation of the apparatus of the present embodiment will be described.
[0187]
Here, first, a control cycle of processing performed by the
[0188]
On the other hand, the determination processing of the target air-fuel ratio KCMD of the exhaust gas upstream of the
[0189]
In addition, what is necessary is just to determine this fixed period according to the kind, reaction rate, volume, etc. of the
[0190]
Assuming the above, first, referring to the flowcharts of FIGS. 6 and 7, the output fuel injection amount #nTout (n = 1, n) for each cylinder of the
[0191]
First, referring to FIG. 6, the
[0192]
Next, the basic fuel injection amount Tim obtained by correcting the fuel injection amount corresponding to the engine speed NE and the intake pressure PB according to the effective opening area of the throttle valve by the basic fuel injection
[0193]
Next, the
[0194]
In the above discrimination processing, as shown in FIG. 2 It is determined whether or not the
[0195]
Whether or not the
[0196]
Further, it is determined whether or not the engine speed NE and the intake pressure PB are within a predetermined range (STEP d-7, d-8). Since it is not preferable to control the fuel supply of the
[0197]
When the conditions of STEPd-1, d-2, d-7, and d-8 are satisfied and the conditions of STEPd-3, d-4, d-5, and d-6 are not satisfied, The value of the flag f / prism / on is set to “1” in order to use the target air-fuel ratio KCMD generated by the air-fuel ratio manipulated
[0198]
Returning to FIG. 6, after setting the value of the flag f / prism / on as described above, the
[0199]
Next, the
[0200]
In this case, as described above, the global feedback control unit 15 divides the feedback manipulated variable KLAF obtained by the
[0201]
After the feedback correction coefficient KFB is calculated as described above, the second correction coefficient KCMDM corresponding to the target air-fuel ratio KCMD determined in STEPf or STEPg is further calculated by the second correction coefficient calculator 10 ( STEPj).
[0202]
Next, the
[0203]
In the
[0204]
The calculation of the output fuel injection amount #nTout for each cylinder as described above and the fuel injection to the
[0205]
On the other hand, in parallel with the fuel control of the
[0206]
That is, referring to the flowchart of FIG. 8, the air-fuel ratio manipulated
[0207]
The above determination processing is performed as shown in the flowchart of FIG.
[0208]
That is, as in the case of STEPd in FIG. 2 It is determined whether or not the
[0209]
Further, whether or not the
[0210]
Returning to FIG. 8, after performing the discrimination processing as described above, the air-fuel ratio manipulated
[0211]
The determination process of STEP2 is performed as shown in the flowchart of FIG.
[0212]
That is, whether the throttle valve of the
[0213]
Returning to FIG. 8, the air-fuel ratio manipulated
[0214]
Next, the air / fuel ratio manipulated
[0215]
On the other hand, if it is determined in
[0216]
The arithmetic processing by the
[0217]
That is, the
[0218]
On the other hand, if f / id / cal = 1, the
[0219]
Next, the
[0220]
Further, the
[0221]
Here, first, O 2 The relationship between the behavior of the output VO2 / OUT or the deviation output VO2 of the
[0222]
Referring to FIG. 2 The output VO2 / OUT or the deviation output VO2 of the
[0223]
Therefore, in this embodiment, the identification (update) of the gain coefficients a1, a2, and b1 by the
[0224]
On the other hand, referring to FIG. 13, according to the control of the present embodiment using the adaptive sliding mode control, O 2 The state quantity X (VO2 (k), VO2 (k-1)) of the deviation output VO2 of the
[0225]
Therefore, the determination as to whether or not the air-fuel ratio of the exhaust gas that has passed through the
[0226]
For this reason, in the present embodiment, a management function γ defined using time series data of the deviation output VO2 is introduced by the following equation (50),
[0227]
[Equation 50]
[0228]
The management hyperplane (in this case, a straight line) represented by the coefficients m1, m2, and m3 of the management function γ is represented by γ = 0, as shown in FIG. 13, the sliding plane control hyperplane σ = It was set to be slightly above 0 (region where σ> 0). In this embodiment, the coefficient s of the linear function σ 1 Is set to “1”, the coefficient m1 of the management function γ is set to “1”.
[0229]
When such a management function γ is introduced, in a state where γ ≧ 0, the air-fuel ratio of the exhaust gas surely changes from the lean side to the rich side. This can be performed stably depending on whether the value of the function γ is positive (including “0”).
[0230]
The management processing in STEP 5-6 is performed using the management function γ defined as described above. 2 A behavior state in which the air-fuel ratio of the exhaust gas indicated by the deviation output VO2 of the
[0231]
That is, referring to the flowchart of FIG. 14, the
[0232]
Next, the
[0233]
Thereby, it is a behavior state suitable for the identification (update) of the gain coefficients a1, a2, and b1 by the
[0234]
Returning to the description of FIG. 11, after performing the management process as described above, the
[0235]
Then, the
[0236]
Here, the identification error id / e in STEP 5-8 may be basically calculated according to the equation (6). In the present embodiment, for example, FIG. 15A shows a block diagram. As described above, after filtering the same characteristics to the deviation output VO2 acquired for each control cycle in STEP3 (see FIG. 8) and the identification deviation output VO2 hat calculated for each control cycle in STEP5-4, The identification error id / e is calculated in STEP5-8.
[0237]
That is, referring to FIG. 16, the output change (O of the target exhaust system E with respect to the input change of the target exhaust system E including the catalyst device 3 (change in the output KACT or deviation output kact of the LAF sensor 5). 2 The gain frequency characteristic of the output VO2 / OUT or the deviation output VO2 of the
[0238]
The sliding
[0239]
On the other hand, the arithmetic processing (see formulas (7) to (9)) for identifying the gain coefficients a1, a2, and b1 of the discrete system model by the
[0240]
Therefore, in the present embodiment, as shown by the alternate long and short dash line in FIG. 16, after filtering the characteristic having a weight in the low frequency band C (low-pass characteristic) on the deviation output VO2 and the identification deviation output VO2 hat, STEP5 The identification error id / e at -8 is calculated.
[0241]
In this low pass characteristic filtering process, the deviation output VO2 acquired in STEP 3 (FIG. 8) and the identification deviation output VO2 hat calculated in STEP 5-4 are respectively stored in time series. For each of the time-series data of the deviation output VO2 and the identified deviation output VO2 hat, the calculation is performed by calculating the addition average or the weighted addition average of a predetermined number of data that has been reversed from the present to the past for each control cycle. . This is a technique of a digital filter and is a technique generally called moving average processing. Then, the calculation of the identification error id / e in STEP5-8 is performed by subtracting the filtering value of the identification deviation output VO2 hat from the filtering value of the deviation output VO2 obtained by the above moving average process.
[0242]
By performing such filtering processing, the frequency characteristics of the discrete system model determined by the identification gain coefficients a1 hat, a2 hat, and b1 hat obtained from the identification error id / e by the above equation (7) are shown in FIG. As shown to a), it can be set as the frequency characteristic of the tendency similar to the frequency characteristic of the actual target exhaust system E. FIG.
[0243]
In this case, in this embodiment, as shown in FIG. 17A, the gain at each frequency of the discrete system model is further slightly larger than the gain at each frequency of the target exhaust system E as a whole. Set weighting characteristics for filtering. By doing so, in the discrete system model and the target exhaust system E, a certain output change (specifically, O 2 An input change (change in the output KACT or deviation output kact of the LAF sensor 5) that causes the output VO2 / OUT of the
[0244]
In this embodiment, as shown in FIG. 17A, the weighting characteristic of the filtering is set so that the gain at each frequency of the discrete system model is slightly larger overall than the gain at each frequency of the target exhaust system E. For example, as shown in FIG. 17B, the frequency characteristic of the gain of the discrete system model in the low frequency band C is substantially the same as the frequency characteristic of the gain of the target exhaust system E, and the low frequency band C In the frequency band on the higher frequency side, the filtering weight characteristic is set so that the gain of the discrete system model is slightly larger than the gain of the target exhaust system E, as in the case of FIG. It may be.
[0245]
In this way, on the higher frequency side than the low frequency band C (frequency band in which the gain of the target exhaust system E is relatively small), as in the case of FIG. 2 The stability of the convergence control of the output VO2 / OUT of the
[0246]
Further, in the present embodiment, as shown in the block diagram of FIG. 15A, after the filtering is performed on the deviation output VO2 and the identification deviation output VO2 hat, the identification error id / e is calculated. However, as shown in FIG. 15B, for example, before calculating the identification deviation output VO2 hat in STEP 5-4, the deviation outputs kact and VO2 used for the calculation are filtered with the same characteristics. Alternatively, the identification error id / e may be calculated from the identification deviation output VO2 hat calculated from the filtering values according to the equation (5) and the filtering value of the deviation output VO2 previously filtered. Alternatively, as shown in FIG. 15 (c), the identification error id / e is calculated from the identification deviation output VO2 hat calculated using the deviation outputs kact and VO2 as they are and the deviation output VO2 using the equation (6) as it is. Later, the identification error id / e may be filtered. In other words, the filtering process as described above in calculating the identification error id / e only requires that the deviation output VO2 and the identification deviation output VO2 hat are filtered with the same characteristics, and the timing of the filtering process. Can be arbitrarily selected.
[0247]
Returning to the description of FIG. 11, the
[0248]
On the other hand, if the value of the identification gain coefficient b1 hat is not within the predetermined range as determined in STEP5-11-1, the identification process of the gain coefficients a1, a2, b1 by the
[0249]
In the present embodiment, the stability of the identification process of the
[0250]
Returning to the description of FIG. 11, after the identification gain coefficient vector Θ is evaluated as described above, the
[0251]
The above is the arithmetic processing of the
[0252]
Returning to the description of the main routine process of FIG. 8, after the calculation process of the
[0253]
In addition, when f / id / cal = 0, that is, when the identification processing of the gain coefficients a1, a2, and b1 by the
[0254]
Next, in the main routine of FIG. 8, the air-fuel ratio manipulated
[0255]
The calculation process of the
[0256]
Here, the reason why the deviation output kact of the
[0257]
The filtering as described above is performed by a moving average process as in the case of filtering in the
[0258]
Returning to the description of FIG. 8, the air-fuel ratio manipulated
[0259]
The calculation of the SLD operation input usl is performed as shown in the flowchart of FIG.
[0260]
That is, the sliding
[0261]
Next, the sliding
[0262]
Next, the sliding
[0263]
That is, referring to the flowchart of FIG. 22, the sliding
[0264]
On the other hand, when | b1 | <ε1 (when b1≈0), each of the gain coefficients b1 depends on whether the sign of the gain coefficient b1 is positive (including the case where b1 = 0). The value is forcibly limited to a predetermined positive value (≧ ε1) and a negative predetermined value (≦ −ε1). As described above, it is apparent with reference to the equations (34), (36), and (37) that the magnitude of the gain coefficient b1 is limited to prevent the gain coefficient b1 from becoming excessively small. Since the gain coefficient b1 is used as a denominator term when calculating the equivalent control input ueq, the reaching law input urch, and the adaptive law input uadp, if the magnitude of the gain coefficient b1 is too small, This is because an excessive equivalent control input ueq, reaching law input urch, and adaptive law input uadp are calculated.
[0265]
In the present embodiment, only the value of the gain coefficient b1 is limited. However, the values of the other gain coefficients a1 and a2 may be limited.
[0266]
Returning to the description of FIG. 21, after the limit processing of the gain coefficient b1 as described above, the sliding
[0267]
Further, the sliding
[0268]
Returning to FIG. 8, after calculating the SLD operation input usl as described above, the air-fuel ratio operation
[0269]
This stability determination process is performed as shown in the flowchart of FIG.
[0270]
That is, the air-fuel ratio manipulated
[0271]
Next, the air-fuel ratio manipulated
[0272]
In this case, the state where Δσ bar · σ (k + d) bar>
[0273]
It should be noted that the predetermined value ε2 used in the determination of STEP 9-2 may theoretically be “0”, but is preferably a value slightly larger than “0” in consideration of the influence of stochastic disturbance.
[0274]
On the other hand, if it is determined in STEP 9-2 that Δσ bar · σ (k + d) bar ≦
[0275]
In this case, when the current value σ (k + d) bar of the linear function σ bar is not within the predetermined range, the estimated deviation outputs VO2 (k + d) and VO2 (k + d-1) are the hyperplanes. Since it is in a state far away from σ = 0, it is considered that the adaptive sliding mode control is unstable (the SLD operation input usl calculated in
[0276]
If the current value σ (k + d) bar of the linear function σ bar is within the predetermined range as determined in STEP 9-3, the air-fuel ratio manipulated
[0277]
At this time, if tm> 0, that is, if the timer counter tm has been timed and has not yet timed out, it is determined that the adaptive sliding mode control is unstable according to STEP 9-2 or STEP 9-3. Since the adaptive sliding mode control is likely to become unstable in a state where not much time has elapsed since then, the processing of STEP 9-6 is performed to set the value of the flag f / sld / stb to “0”. Set.
[0278]
If tm ≦ 0 in the judgment of STEP 9-8, that is, if the timer counter tm has timed up, the value of the flag f / sld / stb is set to “ 1 ”(f / sld / stb = 1 indicates that the adaptive sliding mode control is stable) (STEP 9-9).
[0279]
The timer counter tm indicates that the identification error id / e has not converged in the evaluation process of the identification gain coefficients a1 hat, a2 hat, b1 hat in the identifier 25 (the process of the flowchart of FIG. 18). Thus, the identification gain coefficients a1 hat, a2 hat, and b1 hat are activated even when they are not yet fully determined. For this reason, even if the conditions of STEP 9-2 or STEP 9-3 are satisfied, there is a case where tm> 0 in the determination of STEP 9-8, and the value of the flag f / sld / stb is “0”. Set to This is because the adaptive sliding mode control using the identified gain coefficients a1 hat, a2 hat, and b1 hat becomes unstable at the stage where the identified gain coefficients a1 hat, a2 hat, and b1 hat are not yet fully determined. It is easy.
[0280]
Through the processing described above, the stability of the adaptive sliding mode control by the sliding
[0281]
In the present embodiment, the determination of the stability of the adaptive sliding mode control is basically performed based on the condition determination in STEPs 9-2 and 9-3. Alternatively, the stability of the adaptive sliding mode control may be determined by comparing the magnitude (absolute value) of the deviation Δσ bar corresponding to the changing speed of the linear function σ bar with a predetermined value. It is also possible to do.
[0282]
Returning to FIG. 8, after setting the value of the flag f / sld / stb indicating the stability of the adaptive sliding mode control by the sliding
[0283]
Then, after the limit processing of the SLD operation input usl as described above, the air-fuel ratio manipulated
[0284]
Further, when f / sld / stb = 0 is determined in
[0285]
Note that the target air-fuel ratio KCMD finally determined in
[0286]
What has been described above is the detailed operation of the apparatus of the present embodiment.
[0287]
That is, if the operation is summarized, basically, the air-fuel ratio manipulated
[0288]
In this case, the air-fuel ratio manipulated
[0289]
Further, when the sliding
[0290]
Furthermore, in this embodiment, the state quantity X to be controlled by the sliding
[0291]
Further, by using the discrete system model of the target exhaust system E for the arithmetic processing of the sliding
[0292]
Further, even in the
[0293]
Further, in the present embodiment, when the
[0294]
In this embodiment, O 2 The behavior state of the output of the target exhaust system E grasped by the deviation output VO2 of the
[0295]
In the present embodiment, the stability of the identification process by the
[0296]
In this embodiment, when it is determined that the adaptive sliding mode control by the sliding
[0297]
Next, a second embodiment of the present invention will be described with reference to FIGS. In this embodiment, only the contents of the processing of STEP 5-11 by the
[0298]
In this embodiment, values of gain coefficients a1, a2, and b1, which are parameters identified by the
[0299]
In this case, the predetermined condition for limiting the values of the identification gain coefficients a1 hat, a2 hat, and b1 hat is the response delay element (more specifically, the first order) of the discrete system model expressed by the above equation (1). A condition for limiting the combination of the values of the identification gain coefficients a1 hat and a2 hat related to the autoregressive term and the second order autoregressive term to a predetermined combination (hereinafter referred to as a first limiting condition), and the discrete system There is a condition for limiting the value of the identification gain coefficient b1 hat related to the dead time element of the model (hereinafter referred to as a second limiting condition).
[0300]
Here, before explaining these first and second limiting conditions, the reason for limiting the values of the identification gain coefficients a1 hat, a2 hat, and b1 hat will be described.
[0301]
In the first embodiment described above, when the
[0302]
The inventors of the present invention have examined the above phenomenon, and it is identified by the
[0303]
Therefore, in the present embodiment, the first limiting condition and the second limiting condition are appropriately set, and the combination of the identification gain coefficient a1 hat and a2 hat values and the identification gain coefficient b1 hat value are determined according to these conditions. By properly limiting the 2 While ensuring stable controllability of the output VO2 / OUT of the
[0304]
In view of such a background, in the present embodiment, the first restriction condition and the second restriction condition are set as follows.
[0305]
First, regarding the first limiting condition for limiting the combination of the values of the identification gain coefficients a1 hat and a2 hat, according to the study by the inventors of the present application, the target air-fuel ratio KCMD becomes smooth or high-frequency oscillating. Is determined by the values of the gain coefficients a1 and a2, that is, the coefficient values α1 and α2 of the equation (13), that is, the
[0306]
Specifically, as shown in FIG. 24, when coordinate planes having coefficient values α1 and α2 as components are set, points on the coordinate plane determined by the set of coefficient values α1 and α2 are hatched in FIG. Area (triangle Q 1 Q 2 Q Three The area enclosed by (including the border). When this region is hereinafter referred to as an estimation coefficient stable region), the temporal change in the target air-fuel ratio KCMD tends to be smooth. On the contrary, when the point determined by the set of coefficient values α1, α2 deviates from the estimated coefficient stable region, the temporal change in the target air-fuel ratio KCMD becomes high-frequency oscillation, or , O 2 Controllability to the target value VO2 / TARGET of the output VO2 / OUT of the
[0307]
Therefore, the values of the gain coefficients a1 and a2 identified by the
[0308]
In FIG. 24, the triangular area Q represented on the coordinate plane including the estimated coefficient stable area. 1 Q Four Q Three Is a system defined by the following equation (51), that is, VO2 (k) and VO2 (k-1) on the right side of the equation (13) are represented by VO2 (k) bar and VO2 (k-1) bar ( These VO2 (k) bar and VO2 (k-1) bar mean the estimated deviation output obtained by the
[0309]
[Formula 51]
[0310]
That is, the condition for the system represented by the equation (51) to be stable is that the pole of the system (which is given by the following equation (52)) exists in the unit circle on the complex plane.
[0311]
[Formula 52]
[0312]
And the triangular area Q in FIG. 1 Q Four Q Three Is an area that defines a combination of coefficient values α1 and α2 that satisfy the above conditions. Therefore, the estimated coefficient stable region is a region representing a combination of α1 ≧ 0 among combinations of coefficient values α1, α2 that stabilize the system represented by the equation (51).
[0313]
On the other hand, since the coefficient values α1 and α2 are determined by the combination of the gain coefficients a1 and a2, the combination of the gain coefficients a1 and a2 is determined from the combination of the coefficient values α1 and α2 in a reverse calculation. Therefore, the estimated coefficient stable region in FIG. 24 that defines a preferable combination of the coefficient values α1 and α2 can be converted to the coordinate plane in FIG. 25 having the gain coefficients a1 and a2 as coordinate components. On the coordinate plane of FIG. 25, the estimated coefficient stable region is converted into, for example, a region surrounded by an imaginary line in FIG. 25 (a roughly triangular region having irregularities at the bottom, hereinafter referred to as an identification coefficient stable region). The That is, when the point on the coordinate plane of FIG. 25 determined by the pair of gain coefficients a1 and a2 exists in the identification coefficient stable region surrounded by the phantom line in FIG. 25, the values of the gain coefficients a1 and a2 A point on the coordinate plane of FIG. 24 corresponding to a set of the determined coefficient values α1 and α2 exists in the estimated coefficient stable region.
[0314]
Therefore, the first limiting condition for limiting the values of the identification gain coefficients a1 hat and a2 hat obtained by the
[0315]
However, since a part of the boundary of the identification coefficient stable region shown in the phantom line in FIG. 25 (lower part of the figure) has a complex shape with irregularities, the values of the identification gain coefficients a1 hat and a2 hat are practically used. Therefore, the process for limiting the points on the coordinate plane of FIG.
[0316]
Therefore, in the present embodiment, the identification coefficient stable region is defined by, for example, the quadrangle Q surrounded by the solid line in FIG. Five Q 6 Q 7 Q 8 (A region in which the boundary is formed in a straight line; hereinafter referred to as an identification coefficient limited region). In this case, the identification coefficient limited region is a polygonal line (line segment Q) represented by a function expression of | a1 | + a2 = 1 as shown in the figure. Five Q 6 And line segment Q Five Q 8 And a straight line (line segment Q) expressed by a constant function equation a1 = A1L (A1L: constant) 6 Q 7 And a straight line (line segment Q) expressed by a constant function expression a2 = A2L (A2L: constant) 7 Q 8 A straight line including Then, the first restriction condition for restricting the values of the identification gain coefficients a1 hat and a2 hat is set as a point on the coordinate plane of FIG. 25 determined by these values exists in the identification coefficient restriction region, The values are limited so that the points determined by the values of the identification gain coefficients a1 hat and a2 hat exist in the identification coefficient restriction region. In this case, although a part of the lower side of the identification coefficient limiting area deviates from the identification coefficient stable area, the point determined by the values of the identification gain coefficients a1 hat and a2 hat obtained by the
[0317]
Note that the method of setting the identification coefficient restriction area is exemplary, and the identification coefficient restriction area is basically equal to the identification coefficient stability area or approximately approximate the identification coefficient stability area. Alternatively, any shape may be set as long as most or all of the identification coefficient limited region belongs to the identification coefficient stable region. That is, the identification coefficient restriction region can be set in various ways in consideration of the ease of restriction processing of the values of the identification gain coefficients a1 and a2 and the actual controllability. For example, in this embodiment, the boundary of the upper half part of the identification coefficient restriction region is defined by a function expression of | a1 | + a2 = 1. The combination of the values of the gain coefficients a1 and a2 satisfying this function expression is as described above. This is a combination of theoretical stability limits such that the poles of the system given by equation (52) lie on the unit circumference on the complex plane. Therefore, the boundary of the upper half of the identification coefficient limiting region is, for example, a function expression of | a1 | + a2 = r (where r is a value slightly smaller than “1” corresponding to the stability limit, for example, 0.99). The stability of control may be further increased.
[0318]
Also, the identification coefficient stability region of FIG. 25 that is the basis of the identification coefficient restriction region is also exemplary, and the identification coefficient stability region corresponding to the estimated coefficient stability region of FIG. 24 is derived from the definitions of the coefficient values α1 and α2. As is clear (see equations (12) and (13)), the shape of the identification coefficient stable region changes depending on the value of the dead time d (more precisely, the set value). To do. In this case, regardless of the shape of the identification coefficient stable region, the identification coefficient restriction region may be set as described above according to the shape of the identification coefficient stable region.
[0319]
Next, the second limiting condition for limiting the value of the gain coefficient b1 identified by the
[0320]
That is, according to the knowledge of the inventors of the present application, the situation in which the temporal change of the target air-fuel ratio KCMD is high-frequency oscillation is also the case where the value of the identification gain coefficient b1 hat is too large or too small. It is likely to occur. Therefore, in the present embodiment, the upper limit value B1H and the lower limit value B1L (B1H>B1L> 0) of the value of the identification gain coefficient b1 hat are determined in advance through experiments and simulations. Then, the second limiting condition is set such that the value of the identification gain coefficient b1 hat is not more than the upper limit value B1H and not less than the lower limit value B1L (satisfying the inequality of B1L ≦ b1 hat ≦ B1H).
[0321]
The processing for limiting the values of the identification gain coefficients a1 hat, a2 hat, and b1 hat by the first restriction condition and the second restriction condition set as described above is the identification gain coefficient in STEP 5-11 in this embodiment. The evaluation process of the vector Θ is performed as follows.
[0322]
That is, referring to the flowchart of FIG. 26, the
[0323]
Specifically, the
[0324]
At this time, if a2 (k) hat <A2L, a point on the coordinate plane of FIG. 25 determined by a set of identification gain coefficient a1 (k) hat and a2 (k) hat values (hereinafter, this point (a1 (k) hat and a2 (k) hat) deviate from the identification coefficient restriction region, and the value of a2 (k) hat is forcibly changed to the lower limit value A2L (STEP 5-11-6). ). By this processing, the points (a1 (k) hat, a2 (k) hat) on the coordinate plane in FIG. 25 are at least a straight line (line segment Q) represented by a2 = A2L. 7 Q 8 To a point on the upper side (including the straight line).
[0325]
Next, the
[0326]
At this time, if a1 (k) hat <A1L, or if a1 (k) hat> A1H, the points on the coordinate plane of FIG. 25 (a1 (k) hat, a2 (k) hat) Deviates from the identification coefficient limit region, and the value of the a1 (k) hat is forcibly changed to the lower limit value A1L or the upper limit value A1H according to each case (STEPs 5-11-8, 5- 11-10).
[0327]
With this processing, the points (a1 (k) hat, a2 (k) hat) on the coordinate plane in FIG. 25 are represented by a straight line (line segment Q) represented by a1 = A1L. 6 Q 7 And a straight line represented by a1 = A1H (point Q) 8 And a straight line passing through the a1 axis) (including both straight lines).
[0328]
Note that the order of STEPs 5-11-7 and 5-11-8 and the processes of STEPs 5-11-9 and 5-11-10 may be interchanged. The processing of STEPs 5-11-5 and 5-11-6 may be performed after the processing of STEPs 5-11-7 to 5-11-10.
[0329]
Next, the
[0330]
At this time, if the inequality | a1 | + a2 ≦ 1 holds, a point (a1) determined by the values of a1 (k) hat and a2 (k) hat that have passed through the STEPs 5-11-5 to 5-11-10. (k) hat and a2 (k) hat) exist in the identification coefficient limited region (including its boundary).
[0331]
On the other hand, when | a1 | + a2> 1, the point (a1 (k) hat, a2 (k) hat) deviates from the identification coefficient restriction region to the upper side. In this case, , A2 (k) hat value is forcibly changed to a value (1- | a1 (k) hat |) corresponding to the value of a1 (k) hat (STEP 5-11-12). In other words, the point (a1 (k) hat, a2 (k) hat) is kept on the polygonal line represented by the function expression of | a1 | + a2 = 1 (identification) while keeping the value of the a1 (k) hat as it is. Line segment Q that is the boundary of the coefficient limited area Five Q 6 Top or line segment Q Five Q 8 Move to top).
[0332]
By the processing of STEPs 5-11-5 to 5-11-12 as described above, the values of the identification gain coefficients a1 (k) hat and a2 (k) hat are determined by those values (a1 (k) hat, a2 (k) hat) is limited to exist in the identification coefficient limiting region. Note that points (a1 (k) hat, a2 (k) hat) corresponding to the values of the identification gain coefficients a1 (k) hat and a2 (k) hat obtained in STEP 5-10 are within the identification coefficient restriction region. If present, those values are retained.
[0333]
In this case, the identification gain coefficient a1 (k) hat related to the first-order autoregressive term of the discrete system model is determined by the above-described processing so that the values are the lower limit value A1L and the upper limit value A1H in the identification coefficient limit region. As long as the value is between, the value is not forcibly changed. When a1 (k) hat <A1L, or when a1 (k) hat> A1H, the value of the identification gain coefficient a1 (k) hat is the gain coefficient a1 in the identification coefficient restriction region, respectively. Is forcibly changed to the lower limit value A1L, which is the minimum value that can be taken, and the lower limit value A1H, which is the maximum value that can be taken by the gain coefficient a1 in the identification coefficient restriction region, so that the identification gain coefficient a1 (k ) The amount of change in the hat value is minimal. That is, the points (a1 (k) hat, a2 (k) hat) corresponding to the values of the identification gain coefficients a1 (k) hat and a2 (k) hat obtained in STEP 5-10 deviate from the identification coefficient restriction region. The forced change of the identification gain coefficient a1 (k) hat value is kept to a minimum.
[0334]
In this way, after limiting the values of the identification gain coefficient a1 (k) hat and a2 (k) hat, the
[0335]
That is, the
[0336]
Further, the
[0337]
By such processing of STEPs 5-11-13 to 5-11-16, the value of the identification gain coefficient b1 (k) hat is limited to a value in a range between the lower limit value B1L and the upper limit value B1H.
[0338]
Thus, after limiting the combination of the identification gain coefficient a1 (k) hat and the value of the a2 (k) hat and the identification gain coefficient b1 (k) hat value, the
[0339]
On the other hand, if it is determined in STEP 5-11-17 that | id / e |> ε0, the timer counter used in the determination of the stability of the sliding mode control as described in the first embodiment. After the value of tm (countdown timer) is set to a predetermined initial value TM (timer counter tm activation. STEP 5-11-18), the process returns to the flowchart of FIG.
[0340]
Other operations than those described above are exactly the same as those of the first embodiment. In this case, the previous value a1 (k-1) hat of the identification gain coefficient used for obtaining the identification gain coefficients a1 (k) hat, a2 (k) hat, b1 (k) hat in STEP5-10 of FIG. , A2 (k-1) hat, b1 (k-1) hat are the values of the identification gain coefficients that are limited by the first and second limiting conditions as described above in the processing of STEP 5-11 in the previous control cycle. is there. Therefore, the values of the identification gain coefficients a1 (k) hat, a2 (k) hat, and b1 (k) hat obtained in STEP 5-10 in each control cycle are values satisfying the first and second limiting conditions. Easier to fit.
[0341]
In this embodiment, the value of the identification gain coefficient b1 is limited to a positive value by the processing of STEPs 5-11-13 to 5-11-16, so that STEP8- in FIG. The
[0342]
According to the present embodiment described above, the same effects as those of the first embodiment can be obtained. However, the values of the identification gain coefficients a1 hat, a2 hat, and b1 hat obtained by the
[0343]
In this case, in particular, for the identification gain coefficients a1 hat and a2 hat related to the response delay element of the discrete system model of the target exhaust system E, those values are not restricted individually, but those values are set to the values of both. By limiting by combinations with correlation, O 2 In order to control the output VO2 / OUT of the
[0344]
Further, when limiting the combination of the identification gain coefficient a1 hat and a2 hat values, it relates to the lower-order autoregressive term (first-order autoregressive term) of the autoregressive terms on the right side of Equation (1). The identification gain coefficient a1 hat, in other words, in the discrete system model expressed by equation (1), O 2 More reliable by limiting the combination of a1 hat and a2 hat values so that the amount of change in the value of the identified gain coefficient a1 hat related to the newer output VO2 / OUT or deviation output VO2 of the
[0345]
Further, since the boundary of the identification coefficient limiting region (see FIG. 25) for limiting the combination of the identification gain coefficient a1 hat and a2 hat values is set to be linear, the values of the a1 hat and a2 hat are limited. Can be easily performed.
[0346]
The air-fuel ratio control apparatus for an exhaust system of an internal combustion engine of the present invention is not limited to the first and second embodiments described above, and for example, the following modifications are possible. .
[0347]
That is, in the first and second embodiments, the LAF sensor (wide area air-fuel ratio sensor) 5 is used as the second exhaust gas sensor. However, if the second exhaust gas sensor can detect the air-fuel ratio of the exhaust gas, Normal O 2 Other types of sensors such as sensors may be used.
[0348]
In the first and second embodiments, the first exhaust gas sensor is O. 2 Although the
[0349]
In the first and second embodiments, in the discrete system model of the target exhaust system E and the arithmetic processing of the
[0350]
In the first and second embodiments, the operation amount determined by the air-fuel ratio operation
[0351]
In the first and second embodiments, the sliding
[0352]
In the first and second embodiments, the sliding
[0353]
In the first and second embodiments, when it is determined that the adaptive sliding mode control is unstable, the SLD operation input usl is forcibly set to a predetermined value. Therefore, the target air-fuel ratio KCMD is also set to a predetermined value. However, when it is determined that the adaptive sliding mode control is unstable, for example, a controller provided separately from the air-fuel ratio manipulated
[0354]
In the first and second embodiments, the arithmetic processing by the
[0355]
Furthermore, when the dead time of the exhaust system is sufficiently small, the
[0356]
In the first and second embodiments, the target air-fuel ratio KCMD is determined using the gain coefficients a1, a2, and b1 of the discrete system model identified by the
[0357]
In the first and second embodiments, the dead time d of the target exhaust system E is set to a predetermined value. However, the dead time d may be identified together with the gain coefficients a1, a2, and b1. Is possible. In this case, the value of the dead time d to be identified may be limited by an appropriate condition as in the second embodiment.
[0358]
In the first and second embodiments, the identification of the gain coefficients a1, a2, and b1 by the
[Brief description of the drawings]
FIG. 1 is an overall system configuration diagram of a first embodiment of an air-fuel ratio control apparatus for an internal combustion engine according to the present invention.
FIG. 2 shows the O used in the apparatus of FIG. 2 The output characteristic figure of a sensor.
FIG. 3 is a block diagram showing a basic configuration of an air-fuel ratio manipulated variable determining unit of the apparatus of FIG.
FIG. 4 is an explanatory diagram for explaining a sliding mode control used in the apparatus of FIG. 1;
5 is a block diagram showing a basic configuration of an adaptive controller of the apparatus of FIG.
6 is a flowchart for explaining processing related to fuel control of the engine of the apparatus of FIG. 1;
7 is a flowchart for explaining a subroutine process in the flowchart of FIG. 6;
8 is a flowchart for explaining overall processing of an air-fuel ratio manipulated variable determining unit of the apparatus of FIG.
FIG. 9 is a flowchart for explaining subroutine processing of the flowchart of FIG. 8;
10 is a flowchart for explaining subroutine processing of the flowchart of FIG. 8;
11 is a flowchart for explaining subroutine processing of the flowchart of FIG. 8;
12 is an explanatory diagram for explaining a subroutine process of the flowchart of FIG. 11;
13 is an explanatory diagram for explaining a subroutine process of the flowchart of FIG. 11. FIG.
14 is a flowchart for explaining subroutine processing of the flowchart of FIG. 11. FIG.
15 is an explanatory diagram for explaining a subroutine process of the flowchart of FIG. 11;
16 is an explanatory diagram for explaining a subroutine process of the flowchart of FIG. 11;
FIG. 17 is an explanatory diagram for explaining subroutine processing of the flowchart of FIG. 11;
18 is a flowchart for explaining subroutine processing of the flowchart of FIG. 11;
FIG. 19 is a flowchart for explaining subroutine processing of the flowchart of FIG. 8;
20 is a flowchart for explaining subroutine processing of the flowchart of FIG. 8;
FIG. 21 is a flowchart for explaining a subroutine process of the flowchart of FIG. 8;
22 is a flowchart for explaining subroutine processing of the flowchart of FIG. 21;
FIG. 23 is a flowchart for explaining subroutine processing of the flowchart of FIG. 8;
FIG. 24 is an explanatory diagram for explaining a second embodiment of the air-fuel ratio control apparatus for an internal combustion engine of the present invention.
FIG. 25
Explanatory drawing for demonstrating 2nd Embodiment of the air fuel ratio control apparatus of the internal combustion engine of this invention
.
FIG. 26
For explaining the operation of the second embodiment of the air-fuel ratio control apparatus of the internal combustion engine of the present invention
flowchart.
[Explanation of symbols]
DESCRIPTION OF
Claims (5)
前記内燃機関の排ガスの空燃比を規定する操作量をスライディングモード制御を用いて決定する空燃比操作量決定手段と、該空燃比操作量決定手段が制御対象とする系のモデルの設定すべきパラメータを同定する同定手段とを備え、前記空燃比操作量決定手段は、前記モデルのパラメータの値として、前記同定手段が同定した前記パラメータの同定値と所定値とを内燃機関の運転状態に応じて選択的に用い、該モデルのパラメータの値により規定される係数を含むスライディングモード制御の所定の演算式により前記操作量を決定するようにしたことを特徴とする内燃機関の空燃比制御装置。An exhaust gas sensor provided in the exhaust system of the internal combustion engine for detecting the concentration of a specific component in the exhaust gas of the internal combustion engine, and the air-fuel ratio of the exhaust gas of the internal combustion engine so that the output of the exhaust gas sensor becomes a predetermined target value In an air-fuel ratio control apparatus for an internal combustion engine that controls
Air-fuel ratio manipulated variable determining means for determining the manipulated variable that defines the air-fuel ratio of the exhaust gas of the internal combustion engine using sliding mode control, and parameters to be set in the model of the system to be controlled by the air-fuel ratio manipulated variable determiner a identification means for identifying, the air-fuel ratio manipulated variable determining means, as the value of the parameter of the model, according to identified values of the parameters before Symbol identifying means has identified with a predetermined value to the operating state of the internal combustion engine An air-fuel ratio control apparatus for an internal combustion engine, wherein the operation amount is determined by a predetermined arithmetic expression of sliding mode control including a coefficient defined by a parameter value of the model .
前記内燃機関の排ガスの空燃比を規定する操作量を決定する空燃比操作量決定手段と、該空燃比操作量決定手段が制御対象とする系であって前記排ガスセンサの出力を生成する系が有する無駄時間後における該排ガスセンサの出力の推定値を示すデータを生成する推定手段と、前記制御対象の系のモデルの設定すべきパラメータを同定する同定手段とを備え、
前記推定手段は、前記モデルのパラメータの値として、前記同定手段が同定した前記パラメータの同定値と所定値とを内燃機関の運転状態に応じて選択的に用い、該モデルのパラメータの値により規定される係数を含む推定用の所定の演算式により前記排ガスセンサの出力の推定値を示すデータを生成し、前記空燃比操作量決定手段は、該推定手段が生成したデータにより示される推定値が前記目標値になるように前記操作量を決定することを特徴とする内燃機関の空燃比制御装置。An exhaust gas sensor provided in the exhaust system of the internal combustion engine for detecting the concentration of a specific component in the exhaust gas of the internal combustion engine, and the air-fuel ratio of the exhaust gas of the internal combustion engine so that the output of the exhaust gas sensor becomes a predetermined target value In an air-fuel ratio control apparatus for an internal combustion engine that controls
An air-fuel ratio manipulated variable determining means for determining an operation amount that defines an air-fuel ratio of the exhaust gas of the internal combustion engine, and a system that is controlled by the air-fuel ratio manipulated variable determiner and that generates an output of the exhaust gas sensor. An estimation unit for generating data indicating an estimated value of the output of the exhaust gas sensor after the dead time, and an identification unit for identifying a parameter to be set in the model of the system to be controlled,
The estimating means, as the value of the parameter of the model, selectively used according to identified values of the parameters before Symbol identifying means has identified with a predetermined value to the operating state of the internal combustion engine, the values of the parameters of the model Data indicating an estimated value of the output of the exhaust gas sensor is generated by a predetermined calculation formula for estimation including a specified coefficient, and the air-fuel ratio manipulated variable determining means is an estimated value indicated by the data generated by the estimating means An air-fuel ratio control apparatus for an internal combustion engine, wherein the manipulated variable is determined so that becomes the target value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002125424A JP3748418B2 (en) | 1997-09-16 | 2002-04-26 | Air-fuel ratio control device for internal combustion engine |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP9-251141 | 1997-09-16 | ||
| JP25114197 | 1997-09-16 | ||
| JP2002125424A JP3748418B2 (en) | 1997-09-16 | 2002-04-26 | Air-fuel ratio control device for internal combustion engine |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP10673898A Division JP3592519B2 (en) | 1997-09-16 | 1998-04-16 | Air-fuel ratio control device for exhaust system of internal combustion engine and control device for plant |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2003013778A JP2003013778A (en) | 2003-01-15 |
| JP3748418B2 true JP3748418B2 (en) | 2006-02-22 |
Family
ID=26540070
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002125424A Expired - Lifetime JP3748418B2 (en) | 1997-09-16 | 2002-04-26 | Air-fuel ratio control device for internal combustion engine |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3748418B2 (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4674583B2 (en) * | 2006-12-13 | 2011-04-20 | トヨタ自動車株式会社 | Control device for internal combustion engine |
-
2002
- 2002-04-26 JP JP2002125424A patent/JP3748418B2/en not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JP2003013778A (en) | 2003-01-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3354088B2 (en) | Air-fuel ratio control system for exhaust system of internal combustion engine | |
| JP3592519B2 (en) | Air-fuel ratio control device for exhaust system of internal combustion engine and control device for plant | |
| JP3621839B2 (en) | Plant control device | |
| JP3992925B2 (en) | Degradation determination method of exhaust gas purification catalyst device | |
| JP3331159B2 (en) | Plant control equipment | |
| JP3967524B2 (en) | Air-fuel ratio control device for internal combustion engine | |
| JP3484074B2 (en) | Plant control equipment | |
| JP4265704B2 (en) | Air-fuel ratio control apparatus for internal combustion engine and control apparatus for plant | |
| JP3773684B2 (en) | Air-fuel ratio control device for internal combustion engine | |
| JP3373724B2 (en) | Air-fuel ratio control device for internal combustion engine | |
| JP3300598B2 (en) | Air-fuel ratio control device for internal combustion engine | |
| JP3261038B2 (en) | Air-fuel ratio control device for internal combustion engine | |
| JP4490000B2 (en) | Air-fuel ratio control device for internal combustion engine | |
| JP2000179385A (en) | Air-fuel ratio control device for internal combustion engine and control device for plant | |
| JP4312325B2 (en) | Degradation state evaluation method of exhaust gas purification catalyst device | |
| JPWO2000004282A1 (en) | Method for determining deterioration of an exhaust gas purification catalyst device | |
| JP3655145B2 (en) | Air-fuel ratio control device for multi-cylinder internal combustion engine | |
| JP3913940B2 (en) | Air-fuel ratio control device for internal combustion engine | |
| JP4437626B2 (en) | Air-fuel ratio control device for internal combustion engine | |
| JP3261059B2 (en) | Air-fuel ratio control device for internal combustion engine | |
| JP3748418B2 (en) | Air-fuel ratio control device for internal combustion engine | |
| JP3696570B2 (en) | Plant control device | |
| JP3773859B2 (en) | Plant control device and air-fuel ratio control device for internal combustion engine | |
| JP3816387B2 (en) | Control device | |
| JP4074142B2 (en) | Air-fuel ratio control device for internal combustion engine |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050107 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050322 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050520 |
|
| 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: 20051122 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20051125 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081209 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091209 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091209 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101209 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101209 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111209 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111209 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121209 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131209 Year of fee payment: 8 |
|
| EXPY | Cancellation because of completion of term |