(57)【要約】【課題】触媒の下流側空燃比センサ出力に基づく機関の空燃比のフィードバック制御中において、制御対象のむだ時間を逐次同定でき、同フィードバック制御において最適な応答性を常に確保可能とすること。【解決手段】この排気浄化装置は、燃料噴射量Ｆｉの燃料の噴射から下流側空燃比センサ出力Ｖｏｘｓまでの系であるプラントＰ（即ち、内燃機関１０、第１触媒５３、下流側空燃比センサ６７）を表すモデルの伝達関数Ｇ（ｓ）を一次遅れ要素とパデ近似したむだ時間要素とで表し、Ｇ（ｓ）を表す式に係数として現れ、且つこのモデルの応答特性を決定する同定パラメータ（比例定数Ａ、遅れの時定数τ、むだ時間Ｌ）を逐次同定・更新する。そして、同定された同定パラメータに基づいて、フィードバック制御系において最適な応答性を常に確保できるようにフィードバック制御ゲイン（比例・積分・微分ゲインＫｐ，Ｋｉ，Ｋｄ）を逐次変更・更新する。【選択図】  図５
2004360591.tif 000002
【特許請求の範囲】
【請求項１】
内燃機関の運転状態に応じた燃料噴射量の燃料を噴射する燃料噴射手段を備えた内燃機関に適用される排気浄化装置であって、
前記内燃機関の排気通路に配設された触媒と、
前記触媒よりも下流の前記排気通路に配設された下流側空燃比センサと、
前記下流側空燃比センサの出力値が所定の目標値になるように同下流側空燃比センサの出力に基づいてフィードバック制御量を算出し、同フィードバック制御量に基づいて前記燃料噴射手段により噴射される燃料噴射量を補正することで前記機関に供給される混合気の空燃比をフィードバック制御するフィードバック制御手段と、
フィードバック制御系における前記燃料噴射手段による燃料の噴射から前記下流側空燃比センサの出力までで構成される系の一部又は全部のむだ時間を前記フィードバック制御中において逐次同定する同定手段と、
を備えた内燃機関の排気浄化装置。
【請求項２】
請求項１に記載の内燃機関の排気浄化装置において、
前記同定手段は、前記系の一部又は全部を少なくともむだ時間要素を用いて表すとともに、同むだ時間要素をパデ近似により近似した式における前記むだ時間を示すパラメータを逐次同定することにより同むだ時間を逐次同定するように構成された内燃機関の排気浄化装置。
【請求項３】
請求項１又は請求項２に記載の内燃機関の排気浄化装置であって、
前記フィードバック制御量を算出する際に使用されるフィードバック制御ゲインを前記逐次同定されたむだ時間に基づいて前記フィードバック制御中において逐次変更するフィードバック制御ゲイン変更手段を備えた内燃機関の排気浄化装置。
【請求項４】
請求項１乃至請求項３のいずれか一項に記載の内燃機関の排気浄化装置であって、
前記逐次同定されたむだ時間と所定の基準時間との比較結果に基づいて前記触媒が劣化したか否かを判定する触媒劣化判定手段を備えた内燃機関の排気浄化装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、内燃機関の排気通路に触媒を備えるとともに同触媒の下流側の排気通路に空燃比センサを備え、同空燃比センサの出力に基づいて機関に供給される混合気の空燃比をフィードバック制御する内燃機関の排気浄化装置に関する。
【０００２】
【従来の技術】
従来より、この種の排気浄化装置として、例えば、下記特許文献１に開示された装置が知られている。この内燃機関の排気浄化装置は、内燃機関の排気通路に配設された触媒よりも上流側及び下流側の排気通路に上流側空燃比センサ及び下流側空燃比センサをそれぞれ介装していて、少なくとも下流側空燃比センサの出力値と所定の目標値との偏差に基づいて（偏差を比例・積分・微分処理（ＰＩＤ処理）して）制御器によりフィードバック制御量を算出し、同フィードバック制御量に基づいて燃料噴射量を補正することで機関に供給される混合気の空燃比をフィードバック制御するようになっている。
【０００３】
また、この装置は、フィードバック制御系（制御ループ）における上流側空燃比センサの出力から下流側空燃比センサの出力までで構成される系を所定のプラントモデルで表し、前記フィードバック制御中において、上流側空燃比センサの出力データと下流側空燃比センサの出力データとに基づいて、同プラントモデルを表す式に係数として現れる各パラメータを周知の逐次型同定アルゴリズムを用いて逐次同定する。そして、上記フィードバック制御において最適な応答性が常に確保できるように、前記フィードバック制御中において、同定された各パラメータ値を有する前記プラントモデルに基づいて、前記制御器がフィードバック制御量を算出する際に使用する比例ゲイン、積分ゲイン等のフィードバック制御ゲインを逐次変更するようになっている。
【０００４】
【特許文献１】
特開平１１−３２４７６７号公報
【０００５】
【発明が解決しようとする課題】
上記プラントモデルに対応する実際のプラントのむだ時間は、内燃機関の運転条件（例えば、エンジン回転速度）のみならず、触媒の温度、触媒の劣化の程度等によっても時々刻々と大きく変化し得る。また、実際のプラントのむだ時間が変化すると、上記フィードバック制御において最適な応答性を確保するための前記フィードバック制御ゲインの最適値も直ちに変化し得る。従って、上記フィードバック制御において最適な応答性を常に確保するためには、上記プラントモデルのむだ時間を逐次正確に同定することが要求される。
【０００６】
これに対し、上記開示された装置においては、前記プラントモデルのむだ時間は、所定の一定時間に設定されるか、或いは、所定のマップ等を用いてエンジン回転速度に応じて変更され得るようになっていて、いずれにしろ、むだ時間そのものは同定されていない。これは、前記プラントモデルを一般的なむだ時間要素を用いて表したとき、むだ時間を示すパラメータが同プラントモデルを表す式に係数として現れないから同むだ時間を同定できないことに基づく。
【０００７】
従って、上述のようなエンジン回転速度にのみ基づくマップにより同プラントモデルのむだ時間を設定すると、同設定されたむだ時間と実際のプラントのむだ時間とが大きく相違する場合があり、この結果、フィードバック制御ゲインの値が前記最適値とは異なる値となって上記フィードバック制御において最適な応答性が確保できなくなる場合があるという問題があった。
【０００８】
本発明の目的は、内燃機関の排気通路に配設された触媒の下流側の排気通路に下流側空燃比センサを備え、同空燃比センサの出力値に基づいて機関に供給される混合気の空燃比をフィードバック制御する内燃機関の排気浄化装置において、フィードバック制御中において制御対象（プラント）のむだ時間を逐次同定でき、同フィードバック制御において最適な応答性を常に確保可能なものを提供することにある。
【０００９】
【本発明の概要】
本発明の特徴は、内燃機関の運転状態に応じた燃料噴射量の燃料を噴射する燃料噴射手段を備えた内燃機関に適用される排気浄化装置が、前記内燃機関の排気通路に配設された触媒と、前記触媒よりも下流の前記排気通路に配設された下流側空燃比センサと、前記下流側空燃比センサの出力値が所定の目標値になるように同下流側空燃比センサの出力に基づいてフィードバック制御量を算出し、同フィードバック制御量に基づいて前記燃料噴射手段により噴射される燃料噴射量を補正することで前記機関に供給される混合気の空燃比をフィードバック制御するフィードバック制御手段と、フィードバック制御系における前記燃料噴射手段による燃料の噴射から前記下流側空燃比センサの出力までで構成される系の一部又は全部のむだ時間を前記フィードバック制御中において逐次同定する同定手段とを備えたことにある。
【００１０】
ここにおいて、前記同定手段は、例えば、前記燃料の噴射から前記下流側空燃比センサの出力までで構成される系（即ち、前記系の全部）のむだ時間を同定する手段であってもよいし、前記触媒よりも上流の前記排気通路に配設された上流側空燃比センサの出力から同下流側空燃比センサの出力までで構成される系（即ち、前記系の一部）のむだ時間を同定する手段であってもよい。また、前記所定の目標値は理論空燃比に相当する値に設定されることが好適である。また、前記フィードバック制御手段は、前記下流側空燃比センサの出力値と前記所定の目標値との偏差を比例・積分処理（ＰＩ処理）、又は比例・積分・微分処理（ＰＩＤ処理）することによりフィードバック制御量を算出するように構成されることが好適である。
【００１１】
これによれば、フィードバック制御中において前記系の一部又は全部（制御対象（プラント））のむだ時間（及び、その他の同定パラメータ（例えば、遅れの時定数等））を逐次同定できる。従って、フィードバック制御量を算出する際に使用されるフィードバック制御ゲイン（例えば、比例ゲイン、積分ゲイン、微分ゲイン等）を前記逐次同定されたむだ時間（及び、その他の同定パラメータ）に基づいて前記フィードバック制御中において逐次変更するように構成すれば、同フィードバック制御において最適な応答性を常に確保可能となる。
【００１２】
この場合、前記同定手段は、前記系の一部又は全部を少なくともむだ時間要素を用いて表すとともに、同むだ時間要素をパデ近似により近似した式における前記むだ時間を示すパラメータを逐次同定することにより同むだ時間を逐次同定するように構成されることが好適である。
【００１３】
むだ時間要素をパデ近似により近似した式においては、後述するように、むだ時間を示すパラメータが係数として現れる。このことに着目して、上記のように、前記系の一部又は全部を表すためのむだ時間要素をパデ近似により近似した式における前記むだ時間を示すパラメータを逐次同定するように構成すれば、同むだ時間を簡易な計算で逐次同定することができる。
【００１４】
また、上記本発明に係る排気浄化装置においては、前記逐次同定されたむだ時間と所定の基準時間との比較結果に基づいて前記触媒が劣化したか否かを判定する触媒劣化判定手段を備えることが好適である。ここにおいて、前記所定の基準時間は、例えば、前記触媒が新品であると仮定した場合に前記同定手段により同定される既知のむだ時間に基づく時間であって、これに限定されない。
【００１５】
一般に、前記系（の一部又は全部）のような触媒を含む系のむだ時間は、同触媒の劣化の程度に依存し、同触媒の劣化が進行するにつれて短くなることが知られている。従って、上記のように、（或る特定の触媒劣化判定条件（例えば、触媒温度が所定温度の範囲内にあり、車速が所定車速範囲内にあって、内燃機関が定常運転状態にあること等）が成立している場合において）前記逐次同定されたむだ時間と所定の基準時間とを比較することにより、同比較結果に基づいて前記触媒が劣化したか否かを容易に判定することができる。
【００１６】
【発明の実施の形態】
以下、本発明による内燃機関の排気浄化装置の実施形態について図面を参照しつつ説明する。
【００１７】
図１は、この実施形態による排気浄化装置を火花点火式多気筒（４気筒）内燃機関１０に適用したシステムの概略構成を示している。この内燃機関１０は、シリンダブロック、シリンダブロックロワーケース、及びオイルパン等を含むシリンダブロック部２０と、シリンダブロック部２０の上に固定されるシリンダヘッド部３０と、シリンダブロック部２０にガソリン混合気を供給するための吸気系統４０と、シリンダブロック部２０からの排気ガスを外部に放出するための排気系統５０とを含んでいる。
【００１８】
シリンダブロック部２０は、シリンダ２１、ピストン２２、コンロッド２３、及びクランク軸２４を含んでいる。ピストン２２はシリンダ２１内を往復動し、ピストン２２の往復動がコンロッド２３を介してクランク軸２４に伝達され、これにより同クランク軸２４が回転するようになっている。シリンダ２１とピストン２２のヘッドは、シリンダヘッド部３０とともに燃焼室２５を形成している。
【００１９】
シリンダヘッド部３０は、燃焼室２５に連通した吸気ポート３１、吸気ポート３１を開閉する吸気弁３２、吸気弁３２を駆動するインテークカムシャフトを含むとともに同インテークカムシャフトの位相角を連続的に変更する可変吸気タイミング装置３３、可変吸気タイミング装置３３のアクチュエータ３３ａ、燃焼室２５に連通した排気ポート３４、排気ポート３４を開閉する排気弁３５、排気弁３５を駆動するエキゾーストカムシャフト３６、点火プラグ３７、点火プラグ３７に与える高電圧を発生するイグニッションコイルを含むイグナイタ３８、及び燃料を吸気ポート３１内に噴射するインジェクタ（燃料噴射手段）３９を備えている。
【００２０】
吸気系統４０は、吸気ポート３１に連通し同吸気ポート３１とともに吸気通路を形成するインテークマニホールドを含む吸気管４１、吸気管４１の端部に設けられたエアフィルタ４２、吸気管４１内にあって吸気通路の開口断面積を可変とするスロットル弁４３、及びスロットル弁駆動手段を構成するＤＣモータからなるスロットル弁アクチュエータ４３ａを備えている。
【００２１】
排気系統５０は、排気ポート３４に連通したエキゾーストマニホールド５１、エキゾーストマニホールド５１（実際には、各排気ポート３４に連通した各々のエキゾーストマニホールド５１の下流端部が集合した集合部）に接続されたエキゾーストパイプ（排気管）５２、エキゾーストパイプ５２に配設（介装）された上流側の三元触媒５３（上流側触媒コンバータ、又はスタート・キャタリティック・コンバータとも云うが、以下「第１触媒５３」と称呼する。）、及びこの第１触媒５３の下流のエキゾーストパイプ５２に配設（介装）された下流側の三元触媒５４（車両のフロア下方に配設されるため、アンダ・フロア・キャタリティック・コンバータとも云うが、以下「第２触媒５４」と称呼する。）を備えている。排気ポート３４、エキゾーストマニホールド５１、及びエキゾーストパイプ５２は、排気通路を構成している。
【００２２】
一方、このシステムは、熱線式エアフローメータ６１、スロットルポジションセンサ６２、カムポジションセンサ６３、クランクポジションセンサ６４、水温センサ６５、第１触媒５３の上流の排気通路（本例では、上記各々のエキゾーストマニホールド５１が下流端部が集合した集合部）に配設された空燃比センサ６６（以下、「上流側空燃比センサ６６」と称呼する。）、第１触媒５３の下流であって第２触媒５４の上流の排気通路に配設された空燃比センサ６７（以下、「下流側空燃比センサ６７」と称呼する。）、及びアクセル開度センサ６８を備えている。
【００２３】
熱線式エアフローメータ６１は、吸気管４１内を流れる吸入空気の単位時間あたりの質量流量に応じた電圧Ｖｇを出力するようになっている。かかるエアフローメータ６１の出力Ｖｇと、計測された吸入空気量（流量）Ｇａとの関係は、図２に示したとおりである。スロットルポジションセンサ６２は、スロットル弁４３の開度を検出し、スロットル弁開度ＴＡを表す信号を出力するようになっている。カムポジションセンサ６３は、インテークカムシャフトが９０°回転する毎に（即ち、クランク軸２４が１８０°回転する毎に）一つのパルスを有する信号（Ｇ２信号）を発生するようになっている。クランクポジションセンサ６４は、クランク軸２４が１０°回転する毎に幅狭のパルスを有するとともに同クランク軸２４が３６０°回転する毎に幅広のパルスを有する信号を出力するようになっている。この信号は、エンジン回転速度ＮＥを表す。水温センサ６５は、内燃機関１０の冷却水の温度を検出し、冷却水温ＴＨＷを表す信号を出力するようになっている。
【００２４】
上流側空燃比センサ６６は、限界電流式の酸素濃度センサであり、図３に示したように、空燃比Ａ／Ｆに応じた電流を出力し、この電流に応じた電圧である出力値ｖａｂｙｆｓを出力するようになっていて、特に、空燃比が理論空燃比であるときには出力値ｖａｂｙｆｓは上流側目標値ｖｓｔｏｉｃｈになる。図３から明らかなように、上流側空燃比センサ６６によれば、広範囲にわたる空燃比Ａ／Ｆを精度良く検出することができる。
【００２５】
下流側空燃比センサ６７は、起電力式（濃淡電池式）の酸素濃度センサであり、図４に示したように、理論空燃比近傍において急変する電圧である出力値Ｖｏｘｓを出力するようになっている。より具体的に述べると、下流側空燃比センサ６７は、空燃比が理論空燃比よりもリーンのときは略０．１（Ｖ）、空燃比が理論空燃比よりもリッチのときは略０．９（Ｖ）、及び空燃比が理論空燃比のときは所定の目標値に相当するＶｏｘｓｒｅｆ（＝０．５（Ｖ））の電圧を出力するようになっている。アクセル開度センサ６８は、運転者によって操作されるアクセルペダル８１の操作量を検出し、同アクセルペダル８１の操作量Ａｃｃｐを表す信号を出力するようになっている。
【００２６】
電気制御装置７０は、互いにバスで接続されたＣＰＵ７１、ＣＰＵ７１が実行するルーチン（プログラム）、テーブル（ルックアップテーブル、マップ）、及び定数等を予め記憶したＲＯＭ７２、ＣＰＵ７１が必要に応じてデータを一時的に格納するＲＡＭ７３、電源が投入された状態でデータを格納するとともに同格納したデータを電源が遮断されている間も保持するバックアップＲＡＭ７４、並びにＡＤコンバータを含むインターフェース７５等からなるマイクロコンピュータである。インターフェース７５は、前記センサ６１〜６８と接続され、ＣＰＵ７１にセンサ６１〜６８からの信号を供給するとともに、同ＣＰＵ７１の指示に応じて可変吸気タイミング装置３３のアクチュエータ３３ａ、イグナイタ３８、インジェクタ３９、及びスロットル弁アクチュエータ４３ａに駆動信号を送出するようになっている。また、インターフェース７５は、ＣＰＵ７１の指示に応じて警報ランプ８２に点灯を指示するための点灯信号を供給することができるようになっている。
【００２７】
（空燃比フィードバック制御の概要）
次に、上記のように構成された本発明の実施形態による排気浄化装置（以下、「本装置」と云うこともある。）が行う機関に供給される混合気の空燃比のフィードバック制御の概要について説明する。
【００２８】
第１触媒５３（第２触媒５４も同様である。）は、同第１触媒５３に流入するガスの空燃比が理論空燃比であるときに、ＨＣ，ＣＯを酸化するとともにＮＯｘを還元し、これらの有害成分を高い効率で浄化する。また、第１触媒５３は、酸素を吸蔵・放出する機能（酸素吸蔵機能、酸素吸蔵・放出機能）を有し、この酸素吸蔵・放出機能により、空燃比が理論空燃比からある程度まで偏移したとしても、ＨＣ，ＣＯ、及びＮＯｘを浄化することができる。即ち、機関に供給される混合気の空燃比がリーンとなって第１触媒５３に流入するガスにＮＯｘが多量に含まれると、第１触媒５３はＮＯｘから酸素分子を奪って同酸素分子を吸蔵するとともに同ＮＯｘを還元し、これによりＮＯｘを浄化する。また、機関に供給される混合気の空燃比がリッチになって第１触媒５３に流入するガスにＨＣ，ＣＯが多量に含まれると、第１触媒５３は吸蔵している酸素分子を与えて（放出して）酸化し、これによりＨＣ，ＣＯを浄化する。
【００２９】
従って、第１触媒５３が連続的に流入する多量のＨＣ，ＣＯを効率的に浄化するためには、同第１触媒５３が酸素を多量に貯蔵していなければならず、逆に連続的に流入する多量のＮＯｘを効率的に浄化するためには、同第１触媒５３が酸素を十分に貯蔵し得る状態になければならないことになる。以上のことから、第１触媒５３の浄化能力は、同第１触媒５３が貯蔵し得る最大の酸素量（最大酸素吸蔵量）に依存する。
【００３０】
一方、第１触媒５３のような三元触媒は燃料中に含まれる鉛や硫黄等による被毒、或いは触媒に加わる熱により劣化し、これに伴い最大酸素吸蔵量が次第に低下してくる。このように最大酸素吸蔵量が低下した場合であっても、エミッションの排出量を継続的に抑制するには、第１触媒５３から排出されるガスの空燃比が、理論空燃比に極めて近い状態となるように制御する必要がある。
【００３１】
そこで、本発明の実施形態による排気浄化装置は、下流側空燃比センサ６７の出力値が理論空燃比に対応する所定の目標値Ｖｏｘｓｒｅｆ（０．５（Ｖ））となるように、下流側空燃比センサ６７の出力値Ｖｏｘｓに応じて機関１０に供給される混合気の空燃比を制御する。
【００３２】
より具体的に述べると、本装置は、機能ブロック図である図５に示したように、前記空燃比フィードバック制御を行うためのＡ１〜Ａ７の各手段を含んで構成されている。以下、図５を参照しながら各手段について説明していく。
【００３３】
＜燃料噴射量の算出＞
先ず、燃料噴射量の算出について説明すると、筒内吸入空気量算出手段Ａ１は、エアフローメータ６１が計測している吸入空気流量Ｇａと、クランクポジションセンサ６４の出力に基づいて得られるエンジン回転速度ＮＥと、ＲＯＭ７２が記憶しているテーブルＭＡＰＭｃとに基づき今回の吸気行程を迎える気筒の吸入空気量である筒内吸入空気量Ｍｃを求める。
【００３４】
目標空燃比設定手段Ａ２は、内燃機関１０の運転状態であるエンジン回転速度ＮＥ、及びスロットル弁開度ＴＡ等に基づいて目標空燃比ａｂｙｆｒを決定する。この目標空燃比ａｂｙｆｒは、例えば、内燃機関１０の暖機終了後においては、特殊な場合を除き理論空燃比に設定されている。
【００３５】
基本燃料噴射量算出手段Ａ３は、筒内吸入空気量算出手段Ａ１により求められた筒内吸入空気量Ｍｃを目標空燃比設定手段Ａ２により設定された目標空燃比ａｂｙｆｒで除することにより、機関に供給される混合気の空燃比を同目標空燃比ａｂｙｆｒとするための基本燃料噴射量Ｆｂａｓｅを求める。
【００３６】
燃料噴射量算出手段Ａ４は、基本燃料噴射量算出手段Ａ３により求められた基本燃料噴射量Ｆｂａｓｅに後述する制御器Ａ７により求められるフィードバック制御量ＤＦｉを加えることで、下記数１に基づいて（最終）燃料噴射量Ｆｉを求める。本装置は、このように、筒内吸入空気量算出手段Ａ１、目標空燃比設定手段Ａ２、及び基本燃料噴射量算出手段Ａ３を利用して基本燃料噴射量Ｆｂａｓｅを求めるとともに、燃料噴射量算出手段Ａ４により同基本燃料噴射量Ｆｂａｓｅをフィードバック制御量ＤＦｉに基づいて補正することにより得られる燃料噴射量Ｆｉの燃料を今回の吸気行程を迎える気筒に対してインジェクタ３９により噴射する。
【００３７】
【数１】
Ｆｉ＝Ｆｂａｓｅ＋ＤＦｉ
【００３８】
＜フィードバック制御量の算出＞
次に、フィードバック制御量の算出について説明すると、目標値設定手段Ａ５は、上述した目標空燃比設定手段Ａ２と同様、内燃機関１０の運転状態であるエンジン回転速度ＮＥ、及びスロットル弁開度ＴＡ等に基づいて下流側目標空燃比に対応する所定の目標値Ｖｏｘｓｒｅｆを決定する。この目標値Ｖｏｘｓｒｅｆは、例えば、内燃機関１０の暖機終了後においては、特殊な場合を除き理論空燃比に対応する値である０．５（Ｖ）に設定されている（図４を参照。）。また、本例では、目標値Ｖｏｘｓｒｅｆは、同目標値Ｖｏｘｓｒｅｆに対応する下流側目標空燃比が上述した目標空燃比ａｂｙｆｒと常時一致するように設定される。
【００３９】
出力偏差量算出手段Ａ６は、下記数２に基づいて、目標値設定手段Ａ５により設定されている現時点での目標値Ｖｏｘｓｒｅｆから現時点での下流側空燃比センサ６７の出力値Ｖｏｘｓを減じることにより、出力偏差量ＤＶｏｘｓを求める。
【００４０】
【数２】
ＤＶｏｘｓ＝Ｖｏｘｓｒｅｆ−Ｖｏｘｓ
【００４１】
制御器Ａ７は、出力偏差量算出手段Ａ６により求められた出力偏差量ＤＶｏｘｓを比例・積分・微分処理（ＰＩＤ処理）することで、下記数３に基づいてフィードバック制御量ＤＦｉを求める。
【００４２】
【数３】
ＤＦｉ＝Ｋｐ・ＤＶｏｘｓ＋Ｋｉ・ＳＤＶｏｘｓ＋Ｋｄ・ＤＤＶｏｘｓ
【００４３】
上記数３において、Ｋｐ，Ｋｉ，Ｋｄはフィードバック制御ゲインであり、Ｋｐは比例ゲイン（比例定数）、Ｋｉは積分ゲイン（積分定数）、Ｋｄは微分ゲイン（微分定数）である。また、ＳＤＶｏｘｓは出力偏差量ＤＶｏｘｓの時間積分値であり、ＤＤＶｏｘｓは出力偏差量ＤＶｏｘｓの時間微分値である。なお、後述するように、Ｋｐ，Ｋｉ，ＫｄはＰＩＤゲイン算出器Ａ９により逐次変更・更新されている。
【００４４】
このようにして、本装置は、目標値Ｖｏｘｓｒｅｆと下流側空燃比センサ６７の出力値Ｖｏｘｓとの偏差である出力偏差量ＤＶｏｘｓに基づいて、前記基本燃料噴射量Ｆｂａｓｅを補正して最終燃料噴射量Ｆｉを求めるためのフィードバック制御量ＤＦｉを求める。
【００４５】
例えば、機関に供給される混合気の平均的な空燃比がリーンであるために下流側空燃比センサ６７の出力値Ｖｏｘｓが理論空燃比よりもリーンである空燃比に対応した値を示すと、出力偏差量算出手段Ａ６により求められる出力偏差量ＤＶｏｘｓが正の値となるので（図４を参照。）、制御器Ａ７にて求められるフィードバック制御量ＤＦｉは正の値となる。これにより、燃料噴射量算出手段Ａ４にて求められる（最終）燃料噴射量Ｆｉは基本燃料噴射量Ｆｂａｓｅよりも大きくなって、機関に供給される混合気の空燃比がリッチとなるように制御される。
【００４６】
反対に、機関に供給される混合気の平均的な空燃比がリッチであるために下流側空燃比センサ６７の出力値Ｖｏｘｓが理論空燃比よりもリッチである空燃比に対応した値を示すと、出力偏差量算出手段Ａ６により求められる出力偏差量ＤＶｏｘｓが負の値となるので、制御器Ａ７にて求められるフィードバック制御量ＤＦｉは負の値となる。これにより、燃料噴射量算出手段Ａ４にて求められる（最終）燃料噴射量Ｆｉは基本燃料噴射量Ｆｂａｓｅよりも小さくなって、機関に供給される混合気の空燃比がリーンとなるように制御される。
【００４７】
また、制御器Ａ７は積分項Ｋｉ・ＳＤＶｏｘｓを含んでいるので、定常状態では出力偏差量ＤＶｏｘｓがゼロになることが保証される。換言すれば、前記目標値Ｖｏｘｓｒｅｆと下流側空燃比センサ６７の出力値Ｖｏｘｓとの定常偏差がゼロになる。また、定常状態では、出力偏差量ＤＶｏｘｓがゼロになることで比例項Ｋｐ・ＤＶｏｘｓ、微分項Ｋｄ・ＤＤＶｏｘｓが共にゼロとなるから、フィードバック制御係数ＤＦｉは積分項Ｋｉ・ＳＤＶｏｘｓの値そのものとなる。制御器Ａ７においてかかる積分処理が実行されることにより、前記インジェクタ３９の誤差（指示燃料噴射量である燃料噴射量Ｆｉと実際の燃料噴射量の差）、エアフローメータ６１の誤差（吸入空気流量計測値Ｇａと実際の吸入空気流量の差）が補償され、且つ、定常状態において第１触媒５３の下流の空燃比（従って、機関に供給される混合気の空燃比）が前記目標値Ｖｏｘｓｒｅｆに対応する下流側目標空燃比（即ち、理論空燃比）に収束する。以上、燃料噴射量算出手段Ａ４、目標値設定手段Ａ５、出力偏差量算出手段Ａ６、及び制御器Ａ７がフィードバック制御手段に相当する。以上が空燃比フィードバック制御の概要である。
【００４８】
（同定パラメータの同定・更新の概要）
上述した空燃比フィードバック制御系において、上記最終燃料噴射量Ｆｉの燃料の噴射（具体的には、ＣＰＵ７１によるインジェクタ３９への噴射指示）から下流側空燃比センサ６７の出力Ｖｏｘｓまでで構成される系（即ち、内燃機関１０、第１触媒５３、下流側空燃比センサ６７）の全体を一つのプラントＰとみなすと、このプラントＰの実際の応答特性（具体的には、同応答特性を決定する各パラメータ（例えば、応答遅れの時定数、むだ時間等））は、内燃機関１０の運転条件（例えば、エンジン回転速度ＮＥ）のみならず、触媒の温度、触媒の劣化の程度等によっても時々刻々と大きく変化し得る。
【００４９】
また、実際のプラントＰのかかる応答特性が変化すると、上記空燃比フィードバック制御において最適な応答性を確保するための制御器Ａ７の前記フィードバック制御ゲイン（Ｋｐ，Ｋｉ，Ｋｄ）の最適値もそれぞれ直ちに変化し得る。従って、上記フィードバック制御において最適な応答性を常に確保するためには、プラントＰの応答特性を決定する各パラメータ（例えば、応答遅れの時定数、むだ時間等）を逐次正確に同定することが要求される。一方、プラントＰは、一般に、応答遅れ要素とむだ時間要素とを用いたモデルにより簡易的に表すことができる。
【００５０】
そこで、本装置は、プラントＰを表すモデルの伝達関数Ｇ（ｓ）を下記数４に従って表す。下記数４において、右辺第１項は一次遅れ要素を表す項、右辺第２項はむだ時間要素を表す項である。Ａは一次遅れ要素における比例定数、τは一次遅れ要素における遅れの時定数、Ｌはむだ時間要素におけるむだ時間である。
【００５１】
【数４】
2004360591.tif 000003
【００５２】
上記数４において、右辺第２項におけるむだ時間要素を表す式は、一般的なむだ時間要素「ｅ−ｓＬ」を一次のパデ近似により近似した式である。このようにむだ時間要素としてパデ近似により近似した式を使用するのは、むだ時間を表すパラメータＬが伝達関数Ｇ（ｓ）を表す式に係数として現れるようになり、むだ時間を同定できるようになるからである。
【００５３】
なお、上記数４は連続値系における表現であって、これを例えば所謂オイラー法により離散値系における表現に変換すると、上記数４における「ｓ」を「（ｚ−１）／Ｔ」に置き換えて、プラントＰを表すモデルの伝達関数Ｇ（ｚ）は下記数５のように表される。下記数５において、Ｔはサンプリング周期である。
【００５４】
【数５】
2004360591.tif 000004
【００５５】
そして、本装置は、上記数４（又は上記数５）に係数として現れる比例定数Ａ、遅れの時定数τ、及びむだ時間Ｌを同定すべき同定パラメータとして選択し、同定器Ａ８（図５を参照。）は、これら各同定パラメータをＣＰＵ７１の演算周期（サンプリング周期Ｔ）の経過毎に逐次同定・更新する。かかる同定は、燃料噴射量Ｆｉの時系列データと、下流側空燃比センサ６７の出力Ｖｏｘｓの時系列データとに基づいて最小二乗法を用いた逐次型同定アルゴリズムにより行われるが、この計算の詳細は後述する。このようにして、比例定数Ａ、遅れの時定数τ、及びむだ時間Ｌをフィードバック制御中において逐次同定する手段が同定手段に相当する。以上が、同定パラメータの同定・更新の概要である。
【００５６】
（フィードバック制御ゲインの変更・更新の概要）
一方、上記空燃比フィードバック制御において最適な応答性を常に確保するためには、同定器Ａ８による上記各同定パラメータの同定・更新に伴って制御器Ａ７の前記フィードバック制御ゲイン（Ｋｐ，Ｋｉ，Ｋｄ）の値もそれぞれ変更・更新する必要がある。
【００５７】
ここで、制御器Ａ７の伝達関数をＣ（ｓ）とすると、伝達関数Ｃ（ｓ）は下記数６のように表される。
【００５８】
【数６】
2004360591.tif 000005
【００５９】
なお、上記数６は連続値系における表現であって、これを例えば所謂オイラー法により離散値系における表現に変換すると、上記数６における「ｓ」を「（ｚ−１）／Ｔ」に置き換えて、制御器Ａ７の伝達関数Ｃ（ｚ）は下記数７のように表される。下記数７において、Ｔはサンプリング周期である。
【００６０】
【数７】
2004360591.tif 000006
【００６１】
上記フィードバック制御系において、燃料噴射量Ｆｉを入力、下流側空燃比センサ出力Ｖｏｘｓを出力としたときの閉ループ系の伝達関数（閉ループ伝達関数）は、「Ｇ／（１＋Ｃ・Ｇ）」となり、これに上記数４及び数６、又は上記数５及び数７を代入すれば、閉ループ伝達関数の分母は「ｓ」又は「ｚ」の３次式となる。また、上記空燃比フィードバック制御において最適な応答性を確保するために変更・更新する必要がある前記フィードバック制御ゲイン（Ｋｐ，Ｋｉ，Ｋｄ）の個数は３つである。従って、閉ループ伝達関数の分母の次数が同分母を表す式において変更可能なパラメータの個数以下となるから、この閉ループ伝達関数の極を任意に配置することができる。換言すれば、特性方程式の根（特性根）の値を自由に設定することができる。
【００６２】
また、上記特性方程式が上記特性根として負の実根のみを有する場合、フィードバック制御系が振動的にならずに安定して収束し得る系となり、このときの収束の程度（従って、系の応答性）はかかる負の実根の値そのものに依存することが一般に知られている。
【００６３】
そこで、ＰＩＤゲイン算出器Ａ９（図５を参照。）は、同定器Ａ８による上記各同定パラメータの同定・更新が実行される毎に、上記特性方程式が所定の負の一つの実根（３重根）「−α」（αは正の実定数）を有するように、即ち、上記閉ループ伝達関数の分母が「（ｓ＋α）３」なる式にて表すことができるように、所定の３つの連立方程式を解くことで前記フィードバック制御ゲイン（Ｋｐ，Ｋｉ，Ｋｄ）の値をそれぞれ変更・更新する。このようにして、フィードバック制御ゲイン（Ｋｐ，Ｋｉ，Ｋｄ）の値をフィードバック制御中において逐次変更する手段がフィードバック制御ゲイン変更手段に相当する。以上が、フィードバック制御ゲインの変更・更新の概要である。
【００６４】
このようにして、本装置は、同定器Ａ８により、プラントＰを表すモデルの比例定数Ａ、遅れの時定数τ、及びむだ時間Ｌを同定パラメータとしてそれぞれ同定するとともに、同定器Ａ８によりこれら同定パラメータの同定・更新が実行される毎に、同定された同定パラメータの各値に基づいてフィードバック制御ゲイン（Ｋｐ，Ｋｉ，Ｋｄ）の値をそれぞれ変更・更新する。
【００６５】
（実際の作動）
次に、上記排気浄化装置の実際の作動について説明する。
＜フィードバック制御量の計算＞
ＣＰＵ７１は、図６にフローチャートにより示したフィードバック制御量ＤＦｉの計算を行うルーチンを所定時間の経過毎に繰り返し実行している。従って、所定のタイミングになると、ＣＰＵ７１はステップ６００から処理を開始し、ステップ６０５に進んでフィードバック制御条件が成立しているか否かを判定する。このフィードバック制御条件は、例えば、機関１０の冷却水温ＴＨＷが所定温度以上であって、機関１０の一回転当りの吸入空気量（負荷）が所定値以下であるときに成立する。
【００６６】
いま、フィードバック制御条件が成立しているものとして説明を続けると、ＣＰＵ７１はステップ６０５にて「Ｙｅｓ」と判定してステップ６１０に進み、上記数２に従って、下流側目標値Ｖｏｘｓｒｅｆから現時点の下流側空燃比センサ６７の出力値Ｖｏｘｓを減じることにより、出力偏差量ＤＶｏｘｓを求める。
【００６７】
次に、ＣＰＵ７１はステップ６１５に進んで、下記数８に基づき出力偏差量ＤＶｏｘｓの微分値ＤＤＶｏｘｓを求める。
【００６８】
【数８】
ＤＤＶｏｘｓ＝（ＤＶｏｘｓ−ＤＶｏｘｓｏｌｄ）／Δｔ
【００６９】
上記数８において、ＤＶｏｘｓｏｌｄは前回の本ルーチン実行時において後述するステップ６３０にて設定（更新）された出力偏差量ＤＶｏｘｓの前回値である。また、Δｔは本ルーチンの計算周期（前記所定時間）である。
【００７０】
次いで、ＣＰＵ７１はステップ６２０に進み、上記数３に従って、フィードバック制御量ＤＦｉを求める。このとき、フィードバック制御ゲインＫｐ，Ｋｉ，Ｋｄの各値として、後述するルーチンにより逐次変更・更新されている最新値が使用される。次に、ＣＰＵ７１はステップ６２５に進んで、その時点における出力偏差量の積分値ＳＤＶｏｘｓに上記ステップ６１０にて求めた出力偏差量ＤＶｏｘｓを加えて、新たな出力偏差量の積分値ＳＤＶｏｘｓを求め、続くステップ６３０にて、上記ステップ６１０にて求めた出力偏差量ＤＶｏｘｓを出力偏差量ＤＶｏｘｓの前回値ＤＶｏｘｓｏｌｄとして設定した後、ステップ６９５に進んで本ルーチンを一旦終了する。
【００７１】
一方、ステップ６０５の判定時において、フィードバック制御条件が不成立であると、ＣＰＵ７１は同ステップ６０５にて「Ｎｏ」と判定してステップ６３５に進んでフィードバック制御量ＤＦｉの値を「０」に設定し、その後ステップ６９５に進んで本ルーチンを一旦終了する。このように、フィードバック制御条件が不成立であるときは、フィードバック制御量ＤＦｉを「０」として上記フィードバック制御に基づく機関に供給される混合気の空燃比の補正を行わない。
【００７２】
＜空燃比フィードバック制御＞
また、ＣＰＵ７１は、図７にフローチャートにより示した燃料噴射量Ｆｉの計算、及び燃料噴射の指示を行うルーチンを、各気筒のクランク角が各吸気上死点前の所定クランク角度（例えば、ＢＴＤＣ９０°ＣＡ）となる毎に、繰り返し実行するようになっている。従って、任意の気筒のクランク角度が前記所定クランク角度になると、ＣＰＵ７１はステップ７００から処理を開始してステップ７０５に進み、エアフローメータ６１により計測された吸入空気流量Ｇａ、エンジン回転速度ＮＥ、及び前記目標空燃比ａｂｙｆｒに基づいて、機関に供給される混合気の空燃比を同目標空燃比ａｂｙｆｒとするための基本燃料噴射量Ｆｂａｓｅを求める。
【００７３】
次いで、ＣＰＵ７１はステップ７１０に進み、前記基本燃料噴射量Ｆｂａｓｅに前記ステップ６２０にて逐次求められているフィードバック制御量ＤＦｉの最新値を加えた値を（最終）燃料噴射量Ｆｉとして設定する。
【００７４】
次に、ＣＰＵ７１はステップ７１５に進み、同ステップ７１５にて燃料噴射量Ｆｉの燃料を噴射するための指示をインジェクタ３９に対して行った後、ステップ７９５に進み、本ルーチンを一旦終了する。以上により、フィードバック制御量ＤＦｉが反映された燃料噴射量Ｆｉの燃料が吸気行程を迎える気筒に対して噴射されることで上述の空燃比フィードバック制御が実行される。
【００７５】
＜同定パラメータの算出・更新、ＰＩＤゲインの算出・更新＞
次に、同定パラメータの算出・更新、ＰＩＤゲインの算出・更新について説明すると、ＣＰＵ７１は図８にフローチャートにより示したルーチンを所定時間（前記サンプリング周期Ｔ）の経過毎に繰り返し実行している。従って、所定のタイミングになると、ＣＰＵ７１はステップ８００から処理を開始し、ステップ８０５に進んで、現時点での燃料噴射量Ｆｉの最新値を燃料噴射量の今回値Ｆｉ（ｋ）として設定するとともに、現時点での下流側空燃比センサ６７の出力Ｖｏｘｓの最新値を下流側空燃比センサ出力の今回値Ｖｏｘｓ（ｋ）として設定する。なお、添え字の（ｋ）は今回値であることを示し、添え字の（ｋ−１）は前回値であること、添え字の（ｋ−２）は前々回値であることをそれぞれ示す。
【００７６】
次に、ＣＰＵ７１はステップ８１０に進み、燃料噴射量Ｆｉの時系列データである今回値Ｆｉ（ｋ），前回値Ｆｉ（ｋ−１）及び前々回値Ｆｉ（ｋ−２）、下流側空燃比センサ出力Ｖｏｘｓの時系列データである今回値Ｖｏｘｓ（ｋ），前回値Ｖｏｘｓ（ｋ−１）及び前々回値Ｖｏｘｓ（ｋ−２）、並びに、同定パラメータの前回値である比例定数Ａの前回値Ａ（ｋ−１）、遅れの時定数τの前回値τ（ｋ−１）、むだ時間Ｌの前回値Ｌ（ｋ−１）に基づいて、比例定数Ａの今回値Ａ（ｋ）、遅れの時定数τの今回値τ（ｋ）、むだ時間Ｌの今回値Ｌ（ｋ）をそれぞれ同定する。
【００７７】
ここで、比例定数Ａの前回値Ａ（ｋ−１）、遅れの時定数τの前回値τ（ｋ−１）、むだ時間Ｌの前回値Ｌ（ｋ−１）としては、前回の本ルーチン実行時において後述するステップ８２０にて設定・更新されている最新値が使用され、燃料噴射量Ｆｉの前回値Ｆｉ（ｋ−１）及び前々回値Ｆｉ（ｋ−２）、並びに、下流側空燃比センサ出力Ｖｏｘｓの前回値Ｖｏｘｓ（ｋ−１）及び前々回値Ｖｏｘｓ（ｋ−２）としては、前回の本ルーチン実行時において後述するステップ８２５にて設定・更新されている最新値が使用される。
【００７８】
かかる同定は以下のようにして行う。先ず、同定パラメータの前回値である比例定数Ａの前回値Ａ（ｋ−１）、遅れの時定数τの前回値τ（ｋ−１）、及びむだ時間Ｌの前回値Ｌ（ｋ−１）から、同定パラメータの仮の今回値である比例定数Ａの仮の今回値、遅れの時定数τの仮の今回値、及びむだ時間Ｌの仮の今回値をそれぞれ設定する。次いで、これら同定パラメータの仮の今回値と、燃料噴射量Ｆｉの前回値Ｆｉ（ｋ−１）及び前々回値Ｆｉ（ｋ−２）、並びに下流側空燃比センサ出力Ｖｏｘｓの前回値Ｖｏｘｓ（ｋ−１）及び前々回値Ｖｏｘｓ（ｋ−２）とから燃料噴射量Ｆｉの仮の今回値及び下流側空燃比センサ出力Ｖｏｘｓの仮の今回値を算出する。そして、燃料噴射量Ｆｉの仮の今回値と前記ステップ８０５にて取得された燃料噴射量Ｆｉの今回値Ｆｉ（ｋ）との偏差、及び下流側空燃比センサ出力Ｖｏｘｓの仮の今回値と同ステップ８０５にて取得された下流側空燃比センサ出力Ｖｏｘｓの今回値Ｖｏｘｓ（ｋ）との偏差が共に所定の許容範囲内に入るような前記同定パラメータの仮の今回値を、最小二乗法を用いた周知の逐次型同定アルゴリズムを使用して検索し、検索された同定パラメータの仮の今回値である比例定数Ａの仮の今回値、遅れの時定数τの仮の今回値、及びむだ時間Ｌの仮の今回値を、それぞれ、同定パラメータの今回値である比例定数Ａの今回値Ａ（ｋ）、遅れの時定数τの今回値τ（ｋ）、及びむだ時間Ｌの今回値Ｌ（ｋ）として同定する。
【００７９】
次いで、ＣＰＵ７１はステップ８１５に進んで、同定パラメータの値として、ステップ８１０にて同定・更新された比例定数Ａの今回値Ａ（ｋ）、遅れの時定数τの今回値τ（ｋ）、及びむだ時間Ｌの今回値Ｌ（ｋ）を使用したとき、上述のごとく、閉ループ伝達関数の分母が「（ｓ＋α）３」なる式にて表すことができるようにフィードバック制御ゲイン（Ｋｐ，Ｋｉ，Ｋｄ）をそれぞれ変更・更新する。
【００８０】
続いて、ＣＰＵ７１はステップ８２０に進み、次回の本ルーチン実行の準備として、比例定数Ａの今回値Ａ（ｋ）、遅れの時定数τの今回値τ（ｋ）、及びむだ時間Ｌの今回値Ｌ（ｋ）を、それぞれ、比例定数Ａの前回値Ａ（ｋ−１）、遅れの時定数τの前回値τ（ｋ−１）、及びむだ時間Ｌの前回値Ｌ（ｋ−１）として格納する。
【００８１】
そして、ＣＰＵ７１はステップ８２５に進み、同様に、次回の本ルーチン実行の準備として、燃料噴射量Ｆｉの前回値Ｆｉ（ｋ−１）及び今回値Ｆｉ（ｋ）を、それぞれ、燃料噴射量Ｆｉの前々回値Ｆｉ（ｋ−２）及び前回値Ｆｉ（ｋ−１）として格納するとともに、下流側空燃比センサ出力Ｖｏｘｓの前回値Ｖｏｘｓ（ｋ−１）及び今回値Ｖｏｘｓ（ｋ）を、それぞれ、下流側空燃比センサ出力Ｖｏｘｓの前々回値Ｖｏｘｓ（ｋ−２）及び前回値Ｖｏｘｓ（ｋ−１）として格納した後、ステップ８９５に進んで、本ルーチンを一旦終了する。このようにして、前記サンプリング周期Ｔの経過毎に、同定パラメータが同定されていくと共に、フィードバック制御ゲインが更新されていく。
【００８２】
（触媒の劣化判定）
次に、第１触媒５３の劣化判定について説明すると、ＣＰＵ７１は図９にフローチャートにより示したルーチンを所定時間の経過毎に繰り返し実行している。従って、所定のタイミングになると、ＣＰＵ７１はステップ９００から処理を開始し、ステップ９０５に進んで、触媒劣化判定条件が成立しているか否かを判定し、「Ｎｏ」と判定するときには、ステップ９９５に直ちに進んで本ルーチンを一旦終了する。
【００８３】
触媒劣化判定条件は、例えば、暖機運転終了時点から所定時間経過後（具体的には、機関１０の冷却水温ＴＨＷが所定温度未満の状態から同所定温度以上の状態になった時点から所定時間が経過した後）であって、図示しない車速センサにより検出される車速が所定車速範囲内にあるとき（従って、機関１０が定常運転状態にあるとき）に成立する。
【００８４】
いま、触媒劣化判定条件が成立しているものとして説明を続けると、ＣＰＵ７１はステップ９０５にて「Ｙｅｓ」と判定してステップ９１０に進み、先のステップ８１０にて逐次同定されているむだ時間Ｌ（ｋ）の最新値が所定の基準時間Ｌｒｅｆよりも短い（小さい）か否かを判定する。ここで、基準時間Ｌｒｅｆは、第１触媒５３が新品である場合であって上記触媒劣化判定条件が成立している場合において同定手段Ａ８により同定されるべき既知のむだ時間Ｌ０を所定の自然数で除した時間である。このステップ９１０が触媒劣化判定手段に相当する。
【００８５】
ＣＰＵ７１はステップ９１０にて「Ｎｏ」と判定するときにはステップ９９５に直ちに進んで本ルーチンを一旦終了する一方、「Ｙｅｓ」と判定するときには、ステップ９１５に進んで点灯指示信号を警報ランプ８２に供給した後、ステップ９９５に進んで本ルーチンを一旦終了する。これにより、警報ランプ８２が点灯し、運転者は、第１触媒５３が劣化したことを知らしめられ、同第１触媒５３を交換する等の処置を行うことが可能となる。
【００８６】
以上、説明したように、本発明による内燃機関の排気浄化装置の実施形態によれば、燃料噴射量Ｆｉの燃料の噴射から下流側空燃比センサ６７の出力Ｖｏｘｓまでで構成される系であるプラントＰ（即ち、内燃機関１０、第１触媒５３、下流側空燃比センサ６７）を表すモデルの伝達関数Ｇ（ｓ）を一次遅れ要素とパデ近似したむだ時間要素とで表し、かかるモデルの応答特性を決定するパラメータであって同伝達関数Ｇ（ｓ）を表す式に係数として現れるパラメータ（同定パラメータ、比例定数Ａ、遅れの時定数τ、むだ時間Ｌ）を空燃比フィードバック制御中において逐次同定・更新する。これにより、燃料の輸送による遅れから下流側空燃比センサ６７の遅れ等のプラントＰにおいて内在する制御の応答遅れが全て考慮されて同定パラメータが同定され得る。
【００８７】
そして、同定された同定パラメータに基づいて、フィードバック制御系において最適な応答性を常に確保できるようにフィードバック制御ゲイン（比例ゲインＫｐ、積分ゲインＫｉ、微分ゲインＫｄ）をフィードバック制御中においてそれぞれ逐次変更・更新する。従って、機関の運転状態、第１触媒５３の劣化状態等に拘わらず、常に、空燃比フィードバック制御において最適な応答性を確保することが可能となった。
【００８８】
本発明は上記各実施形態に限定されることはなく、本発明の範囲内において種々の変形例を採用することができる。例えば、上記各実施形態においては、プラントＰを燃料噴射量Ｆｉの燃料の噴射から下流側空燃比センサ６７の出力Ｖｏｘｓまでで構成される系として設定しているが、プラントＰを上流側空燃比センサ６６の出力ｖａｂｙｆｓから下流側空燃比センサ６７の出力Ｖｏｘｓまでで構成される系として設定してもよい。この場合、燃料噴射量Ｆｉの時系列データの代わりに上流側空燃比センサ出力ｖａｂｙｆｓの時系列データを参照して同定パラメータがそれぞれ同定されていく。
【図面の簡単な説明】
【図１】本発明の実施形態に係る排気浄化装置を内燃機関に適用したシステムの概略図である。
【図２】図１に示したエアフローメータの出力電圧と計測された吸入空気流量との関係を示したグラフである。
【図３】図１に示した上流側空燃比センサの出力電圧と空燃比との関係を示したグラフである。
【図４】図１に示した下流側空燃比センサの出力電圧と空燃比との関係を示したグラフである。
【図５】図１に示した排気浄化装置が空燃比フィードバック制御を実行する際の機能ブロック図である。
【図６】図１に示したＣＰＵが実行するフィードバック制御量を計算するためのルーチンを示したフローチャートである。
【図７】図１に示したＣＰＵが実行する燃料噴射量計算のためのルーチンを示したフローチャートである。
【図８】図１に示したＣＰＵが実行する同定パラメータの同定・更新、及びＰＩＤゲインの変更・更新を行うためのルーチンを示したフローチャートである。
【図９】図１に示したＣＰＵが実行する触媒の劣化判定を行うためのルーチンを示したフローチャートである。
【符号の説明】
１０…内燃機関、２５…燃焼室、３９…インジェクタ、５２…エキゾーストパイプ（排気管）、５３…三元触媒（第１触媒）、６６…上流側空燃比センサ、６７…下流側空燃比センサ、７０…電気制御装置、７１…ＣＰＵ
【図１】
2004360591.tif 000007
【図２】
2004360591.tif 000008
【図３】
2004360591.tif 000009
【図４】
2004360591.tif 000010
【図５】
2004360591.tif 000011
【図６】
2004360591.tif 000012
【図７】
2004360591.tif 000013
【図８】
2004360591.tif 000014
【図９】
2004360591.tif 000015