"内燃機関の排気通路に配設されるとともに同機関の排気中に含まれるパティキュレートを捕集するパティキュレートフィルタと、前記排気通路に配設されるとともに同配設された部位の温度に応じて変化する出力値を出力する温度センサと、前記温度センサの出力値をセンサ検出温度に変換する温度変換部と、再生開始条件成立時点以降において前記パティキュレートフィルタの温度が所定の指令温度に一致するように同パティキュレートフィルタに加えられるエネルギーの量を前記センサ検出温度に基いて制御する再生制御手段と、を備え、前記パティキュレートフィルタの温度を上昇させて同パティキュレートフィルタに捕集されたパティキュレートを燃焼させることにより同パティキュレートフィルタを再生する内燃機関の排気浄化装置において、前記再生制御手段は、前記再生開始条件成立時点以降において前記パティキュレートフィルタの温度を前記指令温度に一致させるように予め定められたフィードフォワード量のエネルギーを前記パティキュレートフィルタに与えるフィードフォワード手段と、前記再生開始条件成立時点以降において前記センサ検出温度が前記指令温度に基づいて定まる最終的な目標温度である最終目標温度と一致するように前記エネルギーの量をフィードバック制御するフィードバック手段と、を備え、更に、前記フィードバック手段は、前記再生開始条件成立時点から前記フィードフォワード量のエネルギーを前記パティキュレートフィルタに与えた場合における同再生開始条件成立時点からの経過時間に対する前記センサ検出温度を前記最終目標温度に基づいて推定するとともに、その推定したセンサ検出温度を前記フィードバック制御における目標値であるセンサ目標検出温度として発生する温度センサモデルを用いることにより、前記再生開始条件成立時点からの時間経過に応じた前記センサ目標検出温度を算出する目標検出温度算出手段と、前記センサ目標検出温度と前記センサ検出温度とが一致するように前記パティキュレートフィルタに加えられるエネルギーの量を制御することによって前記フィードバック制御を実行するフィードバック制御手段と、前記再生開始条件成立時点から所定期間が経過する特定時点までは前記センサ目標検出温度を用いた前記フィードバック制御の実行を禁止し、前記特定時点以降において前記フィードバック制御の実行を許可するフィードバック制御禁止許可手段と、を含む排気浄化装置。請求項1に記載の内燃機関の排気浄化装置において、前記温度センサモデルはパデ近似により構築されており、前記特定時点は、前記再生開始条件成立時点以降において、前記最終目標温度と前記センサ目標検出温度との差の大きさが所定の閾値温度差よりも大きい状態から同閾値温度差以下となった時点であると定められた排気浄化装置。請求項2に記載の内燃機関の排気浄化装置において、前記閾値温度差は、前記最終目標温度と、前記再生開始条件成立時点以降に最初に現れる前記センサ目標検出温度の極大値と、の差の大きさ以下の値に設定されている排気浄化装置。以下、本発明の実施形態に係る内燃機関の排気浄化装置について図面を参照しながら説明する。構成図1に示したように、この排気浄化装置10は内燃機関20に適用される。機関20はディーゼル機関である。機関20は、本体部30、インテークマニホールド及び吸気管を含む吸気通路部40、並びに、エキゾーストマニホールド及び排気管を含む排気通路部50を備える。本体部30は複数の気筒を有している。複数の気筒のそれぞれには燃料噴射弁31が配設されている。燃料噴射弁31は、後述する電気制御装置60からの指示信号に応答して高圧の燃料を燃焼室内に噴射するようになっている。吸気通路部40の吸気管には、ターボチャージャのコンプレッサ41及び吸気絞り弁42が配設されている。コンプレッサ41は後述するターボチャージャのタービンとともに回転し、吸気を圧縮するようになっている。吸気絞り弁42は吸気管に回動可能に支持されている。吸気絞り弁42は吸気絞り弁アクチュエータ43により回転され、吸気管の吸気通路断面積を変更するようになっている。排気通路部50の排気管には、排気の上流側から下流側に向けて順に、ターボチャージャのタービン51、DOC52及びDPF53が配設されている。タービン51は排気により回転させられ、コンプレッサ41を回転するようになっている。DOC52は、ディーゼル機関用酸化触媒DieselOxidationCatalystである。DOC52は、ハニカム構造体のセラミック製担持体の表面に酸化触媒を担持している。DOC52は、排気に含まれる炭化水素成分未燃成分を触媒反応により燃焼させることにより、排気温度を上昇させるようになっている。これにより、DPF53の温度床温が上昇する。DPF53は、周知の構造を有するセラミック製フィルタである。例えば、DPF53は、ハニカム構造に成形されたコーディライト等の耐熱性セラミックスにより形成される複数のガス流路を備えている。DPF53内に流入した排気は、その複数のガス流路の多孔性の隔壁を通過する。その際、排気中のPM及びアッシュ等が捕集される。排気浄化装置10の電気制御装置60は、CPU、ROM、RAM及びインターフェースを含む周知のマイクロコンピュータを主体として構成された電子制御回路である。電気制御装置60は、アクセル操作量センサ61、エアフローメータ62、機関回転速度センサ63、最上流温度センサ64、第1温度センサ65、第2温度センサ66、差圧センサ67及び車速センサ68と接続されている。アクセル操作量センサ61は、運転者によって操作されるアクセルペダルの操作量を検出し、アクセルペダル操作量AP以下、アクセル操作量APと称呼する。を出力する。エアフローメータ62は、吸気通路部40の吸気管内を流れる空気の量を検出し、その空気の量を表す信号GA以下、空気量GAと称呼する。を出力する。機関回転速度センサ63は、機関20のクランク軸が所定の回転角度だけ回転する毎にパルスを発生する。電気制御装置60は、このパルスに基づいて機関回転速度NEを取得するようになっている。なお、電気制御装置60は図示しないカムポジションセンサからのパルス及び機関回転速度センサ63からのパルスに基づいて、機関20の絶対クランク角を取得するようになっている。最上流温度センサ64は、排気通路部50の排気管であって、タービン51よりも下流位置且つDOC52よりも上流位置に配設されている。最上流温度センサ64は、その配設位置を流れる排気の温度に応じた出力値VT0を発生する。電気制御装置60は、この最上流温度センサ64の出力値VT0を温度に変換し、その温度T0act以下、最上流温度検出値T0actと称呼する。を取得するようになっている。第1温度センサ65は、排気通路部50の排気管であって、DOC52よりも下流位置且つDPF53よりも上流位置に配設されている。第1温度センサ65は、その配設位置を流れる排気の温度に応じた出力値VT1を発生する。電気制御装置60は、この第1温度センサ65の出力値VT1を温度に変換し、その温度T1act以下、第1温度センサ検出温度T1actと称呼する。を取得するようになっている。第1温度センサ検出温度T1actはDOC52から流出し且つDPF53に流入する排気の第1温度センサ65の出力値に基づく検出温度である。よって、第1温度センサ検出温度T1actはDPF53の実際の温度に応じて変化する。第2温度センサ66は、排気通路部50の排気管であって、DPF53の近傍且つDPF53よりも下流位置に配設されている。第2温度センサ66は、その配設位置を流れる排気の温度に応じた出力値VT2を発生する。電気制御装置60は、この第2温度センサ66の出力値VT2を温度に変換し、その温度T2act以下、第2温度センサ検出温度T2actと称呼する。を取得するようになっている。第2温度センサ検出温度T2actはDPF53から流出した排気の第2温度センサ66の出力値に基づく検出温度である。よって、第2温度センサ検出温度T2actはDPF53の実際の温度に応じて変化する。なお、第2温度センサ66はDPF53内に配設されていてもよい。差圧センサ67は、圧力導入管54に配設されている。圧力導入管54の一端は、排気管のタービン51とDOC52との間の位置に接続されている。圧力導入管54の他端は、排気管のDPF53の下流位置に接続されている。従って、差圧センサ67は、DOC52及びDPF53の上流位置の排気圧力と下流位置の排気圧力との差に応じた値を検出し、その検出した値Pd以下、DPF前後差圧Pdと称呼する。を出力する。車速センサ68は、機関20が搭載された車両の速度を検出し、その検出した車速SPDを出力する。電気制御装置60は、更に、吸気絞り弁アクチュエータ43及び複数の燃料噴射弁31等と接続されている。電気制御装置60は、これらに指示信号を送出するようになっている。DPF再生制御の概要次に、上記のように構成された排気浄化装置10によるDPF再生制御の概要について説明する。排気浄化装置10は、DPF53に捕集された堆積したPMの量以下、PM堆積量と称呼する。を周知の手法により推定している。例えば、排気浄化装置10は、DPF前後差圧Pdに基づいてPM堆積量を推定する。PM堆積量が大きくなるほどDPF前後差圧Pdは増大する。PM堆積量は、前回のPM再生制御終了後からの機関の空気量GAの積算値、或いは、メイン燃料噴射量Qfinの積算値に基づいて推定されることもできる。排気浄化装置10は、そのPM堆積量が制御開始許容量PMth以上となったとき即ち、DPF再生制御開始条件成立時、又は、再生開始条件成立時DPF再生制御を実行する。DPF再生制御においてDPF53の温度を上昇するためにDPF53に加えられるエネルギーは、本例においては、所謂ポスト噴射される燃料によりもたらされる。更に、排気浄化装置10は、DPF53に堆積したアッシュの量以下、アッシュ量SASHと称呼する。を後述する周知の手法により推定している。アッシュ量SASHが大きくなるとDPF53の熱容量も大きくなる。そこで、排気浄化装置10は、そのアッシュ量SASHをDPF再生制御に反映する。より具体的には、排気浄化装置10は、第1温度センサモデルを用いることによりDPF再生制御における第1温度センサ65を使用するフィードバック制御の目標値を算出する。以下、この目標値は第1温度センサ目標検出温度T1tgt、又は、第1温度センサ65のセンサ目標検出温度T1tgtと称呼される。第1温度センサモデルは、指令温度に応じて定まる所定の目標温度第1温度センサ最終目標温度T1inに対する第1温度センサ65の出力値に基づく検出温度第1温度センサ検出温度T1actを推定するモデルであり、その推定値を第1温度センサ目標検出温度T1tgtとして出力する。指令温度は、DPF再生制御におけるDPF53の温度床温の最終的な目標温度である。第1温度センサモデルは、第1温度センサ最終目標温度T1inの変化に対する第1温度センサ検出温度T1actの変化のむだ時間に対応するむだ時間パラメータ及び第1温度センサ最終目標温度T1inの変化に対する第1温度センサ検出温度T1actの変化の時定数に対応する遅れパラメータをパラメータとして含む。遅れパラメータは、一次遅れ式により第1温度センサモデルが構築されている場合には時定数Tdであり、Pade近似の二次遅れ式により第1温度センサモデルが構築されている場合には固有周波数ωnである。排気浄化装置10は、この第1温度センサモデルのむだ時間パラメータ及び遅れパラメータを空気量GAに基づいて変更する図6のA、図7のA及び図8のAを参照。。そして、排気浄化装置10は、第1温度センサ検出温度T1actが第1温度センサ目標検出温度T1tgtに一致するように、周知のPI制御又はPID制御に則ってポスト燃料噴射量のフィードバック量第1フィードバック量を算出する。その際、排気浄化装置10は、PI制御に用いるフィードバック制御定数積分ゲインKi1及び比例ゲインKp1等を空気量GAに基づいて変更する図9のA及び図10のAを参照。。同様に、排気浄化装置10は第2温度センサモデルを用いることによりDPF再生制御における第2温度センサ66を使用するフィードバック制御の目標値を算出する。以下、この目標値は第2温度センサ目標検出温度T2tgt又は、第2温度センサ66のセンサ目標検出温度T2tgtと称呼される。第2温度センサモデルは、指令温度に応じて定まる所定の目標温度第2温度センサ最終目標温度T2inに対する第2温度センサ66の出力値に基づく検出温度第2温度センサ検出温度T2actを推定するモデルであり、その推定値を第2温度センサ目標検出温度T2tgtとして出力する。第2温度センサモデルは、第1温度センサモデルと同型の伝達関数により表現される。従って、第2温度センサモデルは、第2温度センサ最終目標温度T2inの変化に対する第2温度センサ検出温度T2actの変化のむだ時間に対応するむだ時間パラメータ及び第2温度センサ最終目標温度T2inの変化に対する第2温度センサ検出温度T2actの変化の時定数に対応する遅れパラメータをパラメータとして含む。但し、これらのむだ時間パラメータ及び遅れパラメータが、第1温度センサモデルにて使用されるむだ時間パラメータ及び遅れパラメータとそれぞれ相違する。排気浄化装置10は、第2温度センサモデルのむだ時間パラメータ及び遅れパラメータを空気量GA及びアッシュ量SASHに基づいて変更する図6のB、図7のB及び図8のBを参照。。そして、排気浄化装置10は、第2温度センサ検出温度T2actが第2温度センサ目標検出温度T2tgtに一致するように、周知のPI制御又はPID制御に則ってポスト燃料噴射量のフィードバック量第2フィードバック量を算出する。その際、排気浄化装置10は、PI制御に用いるフィードバック制御定数積分ゲインKi2及び比例ゲインKp2等を空気量GA及びアッシュ量SASHに基づいて変更する図9のB及び図10のBを参照。。なお、第1温度センサ65をDPF再生制御のフィードバック制御に使用するか即ち、第1フィードバック量を最終的なフィードバック量とするか、第2温度センサ66をDPF再生制御のフィードバック制御に使用するか即ち、第2フィードバック量を最終的なフィードバック量とするかは、以下の点を考慮して決定される。DPF再生制御においてDPF53の温度が収束した場合における実際のDPF53の温度床温Tact以下、実DPF温度Tactと称呼する。と第1温度センサ検出温度T1actとの差乖離温度。DPF再生制御においてDPF53の温度が収束した場合における実DPF温度Tactと第2温度センサ検出温度T2actとの差乖離温度。第1温度センサモデルのむだ時間パラメータ及び遅れパラメータ。第2温度センサモデルのむだ時間パラメータ及び遅れパラメータ。この結果、特に、第2フィードバック量が最終的なフィードバック量として用いられる場合、アッシュ量SASHによりDPF53の熱容量が変化した場合であっても換言すると、アッシュ量SASHに関わらず、DPF53の温度を狙いとする温度に精度良く上昇させることができる。従って、DPF53の再生を効率よく且つDPF53の劣化を招くことなく行うことができる。排気浄化装置10は、上述したように、第1フィードバック量及び第2フィードバック量を算出し、その何れかを用いてポスト燃料噴射量のフィードバック制御DPF再生制御のフィードバック制御を行う。但し、排気浄化装置10は、後述する種々のフィードバック制御条件が成立しなければ、このフィードバック制御を実行しない禁止する。換言すると、フィードバック制御条件が成立したときにのみ、第1フィードバック量及び第2フィードバック量の何れかがポスト燃料噴射量の補正に用いられる。このフィードバック制御条件のうちの一つは、本例においてパデ近似の二次式により構築された第1温度センサモデルに基づいて推定算出される第1温度センサ検出温度T1act即ち、第1温度センサの目標検出温度T1tgtが精度の良い値となっているという条件である。具体的に述べると、図20に示したように、時刻t0にてDPFの再生開始条件が成立し、時刻t0以降から第1温度センサモデルによる温度推定が開始される場合、理想的な推定温度は破線により示されたT1optであるのに対し、第1温度センサモデルにより推定される推定温度は実線により示されたT1mdlのように変化する。即ち、少なくとも温度推定の開始直後である時刻t0~時刻t1の期間において推定温度T1mdlは大きく上下動し、理想的な推定温度T1optから大きく乖離する。従って、排気浄化装置10は、少なくともこのような期間即ち、DPFの再生開始条件が成立した時点から第1所定期間が経過する第1特定時点までの期間、第1フィードバック量の算出を停止するとともに上記フィードバック制御を禁止する。この第1特定時点は、DPFの再生開始条件成立時から予め定められた第1の一定時間が経過した時点と定めることができる。但し、本例においては、第1特定時点は、後述する第1センサ検出温度T1actが最終的に到達すべき温度である第1温度センサ最終目標温度T1inと第1センサモデルにより算出された第1センサ目標検出温度T1tgt=T1mdlとの差ΔT1=T1in-T1tgtの大きさが、第1の閾値温度差ΔT1thよりも大きい状態から第1の閾値温度差ΔT1th以下の状態となった時点図20の点P2及び時刻t2を参照。であると定められている。第1温度センサモデルによれば、少なくともこの第1特定時点時刻t2以降において、第1センサ目標検出温度T1tgtが理想的な第1センサ目標検出温度T1optに沿って変化するからである。同様に、前記フィードバック制御条件のうちの他の一つは、第2温度センサモデルに基づいて推定算出される第2温度センサ検出温度T2act即ち、第2温度センサ目標検出温度T2tgtが精度の良い値となっているという条件である。排気浄化装置10は、少なくともDPFの再生開始条件が成立した時点から第2所定期間が経過する第2特定時点までの期間において第2フィードバック量の算出を停止するとともに上記フィードバック制御を禁止する。この第2特定時点は、DPFの再生開始条件成立時から予め定められた第2の一定時間が経過した時点と定めることができる。但し、本例においては、第2特定時点は、後述する第2センサ検出温度T2actが最終的に到達すべき温度である第2温度センサ最終目標温度T2inと第2センサモデルにより算出された第2センサ目標検出温度T2tgtとの差ΔT2=T2in-T2tgtの大きさが、第2の閾値温度差ΔT2thよりも大きい状態から第2の閾値温度差ΔT2th以下の状態となった時点であると定められている。第2温度センサモデルによれば、第1温度センサモデルと同様、少なくともこの第2特定時点以降において、第2センサ目標検出温度T2tgtが理想的な第2センサ目標検出温度T2optに沿って変化するからである。DPF再生制御の詳細次に、排気浄化装置10によるDPF再生制御の詳細について機能ブロック図である図2を参照しながら説明する。なお、図2に示されたブロックの機能は、実際には電気制御装置60のCPUが後述する処理を行うことにより実現される。DPF指令温度設定部B1は、DPF再生制御開始条件成立時にDPF指令温度Tshireiを設定する。このDPF指令温度Tshireiは、一定値例えば、600°Cであってもよく、機関20の運転条件例えば、空気量GA及び機関回転速度NE及び/又はPM堆積量に応じて変更されてもよい。例えば、DPF指令温度Tshireiは、PM堆積量が高側閾値を超えている場合にはパティキュレートの急激な燃焼を回避するために低めの温度例えば、600°Cに設定され、PM堆積量が高側閾値よりも小さく且つ制御開始許容量PMthよりも大きい場合にはDPF再生制御に要する時間を短縮するために高めの温度例えば、650°Cに設定されてもよい。基本ポスト噴射量演算部B2は、DPF指令温度Tshirei、機関回転速度NE及びメイン燃料噴射量Qfinに基づいて基本ポスト噴射量QPbaseを決定する。基本ポスト噴射量QPbaseは、その量の燃料が図3に示した膨張下死点近傍のポスト噴射期間メイン燃料噴射による燃焼が終了した後の期間において燃料噴射弁31から噴射されたとき、実DPF温度TactがDPF指令温度Tshireiに一致するように、予め実験により定められたフィードフォワード量である。なお、メイン燃料噴射量Qfinは、例えば、アクセル操作量AP及び機関回転速度NEから推定される要求トルクに基づいて別途CPUにより算出される。メイン燃料噴射量Qfinの燃料は、図3に示したように、圧縮上死点近傍のメイン噴射期間において燃料噴射弁31から噴射される。図2に示した第1温度センサ目標検出温度演算部B3は、第1温度センサ最終目標温度T1inを入力し、その第1温度センサ最終目標温度T1inに基づき第1温度センサモデルを用いて第1温度センサ目標検出温度T1tgtを演算出力する。ここで、第1温度センサ最終目標温度T1inは、DPF指令温度Tshireiからオフセット値Offset1を減じた値T1in=Tshirei-Offset1である。このオフセット値Offset1は、図4に示したように、DPF53の実際の温度実DPF温度TactがDPF指令温度Tshireiに一致収束した場合におけるDPF指令温度Tshireiと第1温度センサ検出温度T1actとの差定常的な偏差、乖離温度である。オフセット値Offset1には一定値が設定されもよく、空気量GAに応じて変化する値が設定されてもよい図12のラインC1を参照。。ここで、第1温度センサ目標検出温度演算部B3が使用する第1温度センサモデルM1について説明する。第1温度センサモデルM1は図5に示したように、第1温度センサ最終目標温度T1inを入力とし、第1温度センサ目標検出温度T1tgtを出力するためのモデルである。ところで、図4は、時刻t0特定の時点にてDPF再生制御が開始されてDPF指令温度Tshireiが増大され、これにより、実DPF温度TactをDPF指令温度Tshireiに一致させるように予め適合されている基本ポスト噴射量QPbaseの燃料がポスト噴射された場合についての各値の挙動を示している。図4によれば、第1温度センサ検出温度T1actはむだ時間L1が経過した時刻t1にて増大し始める。むだ時間L1は、ポスト噴射により供給された燃料がDOC52にて燃焼し、それによってDPF53に流入する排気温度が上昇を開始するまでの時間に相当している。更に、図4から理解されるように、第1温度センサ検出温度T1actは、遅れを伴って増大する。即ち、第1温度センサ検出温度T1actの変化は、むだ時間と、一次遅れ又は二次遅れと、を表現する式本例において周知のPade近似二次により近似することができる。そこで、第1温度センサモデルM1は、パデ近似に基づき以下の1式の伝達関数を有するモデルとして構築される。第1温度センサ目標検出温度演算部B3は、上記1式のむだ時間Lとして第1温度センサ65のむだ時間L1を適用するとともに、固有周波数ωnとして第1温度センサ65の固有周波数ωn1を適用する。ところで、実験によれば、図6のAに示したように、むだ時間L1は空気量GAが大きくなるほど小さくなる。空気量GAが大きいほど基本ポスト噴射量QPbaseが大きくなってDOC52従って、DPF53に投入されるエネルギー量が大きくなるからである。更に、図7のAに示したように、固有周波数ωn1は空気量GAが大きくなるほど大きくなる。これは、図8のAに示したように、一次遅れの近似式で第1温度センサモデルを構築した場合の時定数Td1が空気量GAが大きいほど小さくなる即ち、第1温度センサ検出温度T1actが相対的に早く増大することに対応している。係る観点に基づき、第1温度センサ目標検出温度演算部B3は、実際の空気量GAを図6のAに示したテーブルに適用してむだ時間L1を取得し、そのむだ時間L1を上記1式のむだ時間Lに代入する。更に、第1温度センサ目標検出温度演算部B3は、実際の空気量GAを図7のAに示したテーブルに適用して固有周波数ωn1を取得し、その固有周波数ωn1を上記1式の固有周波数ωnに代入する。以上により、第1温度センサ目標検出温度T1tgtが算出される。この第1温度センサ目標検出温度T1tgtは、DPFの再生開始条件成立時点からの時間経過に応じて変化する第1温度センサ検出温度T1actの目標値である。再び図2を参照すると、温度変換部B4は、第1温度センサ65の出力値VT1を関数fを用いて第1温度センサ検出温度T1actに変換する。関数fによれば、出力値VT1が大きいほど第1温度センサ検出温度T1actは高い値となる。第1フィードバックコントローラB5は、第1温度センサ検出温度T1actが第1温度センサ目標検出温度T1tgtに一致するように、ポスト燃料噴射量QPの第1フィードバック量KFB1をPI比例積分制御に則って算出する。即ち、第1フィードバックコントローラB5は、下記の2式に従って第1フィードバック量KFB1を算出する。2式において、ΔT1は第1温度センサ目標検出温度T1tgtから第1温度センサ検出温度T1actを減じた値であるΔT1=第1温度センサ目標検出温度T1tgt-第1温度センサ検出温度T1act。上記2式における積分ゲイン積分定数Ki1は図9のAに示したように空気量GAが大きいほど大きくなるように設定される。同様に、上記2式における比例ゲイン比例定数Kp1は図10のAに示したように空気量GAが大きいほど大きくなるように設定される。積分ゲイン及び比例ゲインはフィードバック制御における制御ゲインであり、フィードバック制御の制御パラメータでもある。図2に示した第2温度センサ目標検出温度演算部B6は、第2温度センサ最終目標温度T2inを入力し、その第2温度センサ最終目標温度T2inに基づき第2温度センサモデルを用いて第2温度センサ目標検出温度T2tgtを演算出力する。ここで、第2温度センサ最終目標温度T2inは、DPF指令温度Tshireiからオフセット値Offset2を減じた値T2in=Tshirei-Offset2である。このオフセット値Offset2は、図4に示したように、DPF53の実際の温度実DPF温度TactがDPF指令温度Tshireiに一致収束した場合におけるDPF指令温度Tshireiと第2温度センサ検出温度T2actとの差定常的な偏差、乖離温度である。オフセット値Offset2には一定値が設定されもよく、空気量GAに応じて変化する値が設定されてもよい図13のラインC1を参照。。ここで、第2温度センサ目標検出温度演算部B6が使用する第2温度センサモデルM2について説明する。第2温度センサモデルM2は、図11に示したように、第2温度センサ最終目標温度T2inを入力とし、第2温度センサ目標検出温度T2tgtを出力するためのモデルである。図4に示したように、第2温度センサ検出温度T2actはむだ時間L1よりも長いむだ時間L2が経過した時刻t2にて増大し始める。むだ時間L2は、ポスト噴射により供給された燃料がDOC52にて燃焼し、それによってDPF53に流入する排気温度が上昇し、更に、その高温の排気がDPF53の温度を増大し始めるまでの時間に対応している。更に、図4から理解されるように、第2温度センサ検出温度T2actは、遅れを伴って増大する。即ち、第2温度センサ検出温度T2actの変化も、むだ時間と、一次遅れ又は二次遅れと、を表現する式本例において周知のPade近似二次により近似することができる。従って、第2温度センサモデルM2は、上述した第1温度センサモデルM1と同じ形の伝達関数即ち、上記1式の伝達関数により表される。但し、第2温度センサ目標検出温度演算部B3は、上記1式のむだ時間Lとして第2温度センサ66のむだ時間L2を適用するとともに、固有周波数ωnとして第2温度センサ66の固有周波数ωn2を適用する。ところで、実験によれば、図6のBに示したように、むだ時間L2は空気量GAが大きくなるほど小さくなる。空気量GAが大きいほど基本ポスト噴射量QPbaseが大きくなってDOC52従って、DPF53に投入されるエネルギー量が大きくなるからである。更に、むだ時間L2は、図6のBのラインC2及びC3により示したように、アッシュ量SASHが大きいほど大きくなる。アッシュ量SASHが大きいほどDPF53の熱容量が増大するとともに、DPF53内を排気が通過する時間が長くなるからである。更に、図7のBに示したように、固有周波数ωn2は空気量GAが大きくなるほど大きくなる。これは、図8のBに示したように、一次遅れの近似式で第2温度センサモデルを構築した場合の時定数Td2が空気量GAが大きいほど小さくなる即ち、第2温度センサ検出温度T2actが相対的に早く増大することに対応している。更に、固有周波数ωn2は、図7のBのラインC2及びC3により示したように、アッシュ量SASHが大きいほど小さくなる。これは、図8のBのラインC2及びC3により示したように、一次遅れの近似式で第2温度センサモデルを構築した場合の時定数Td2がアッシュ量SASHが大きいほど大きくなる即ち、第2温度センサ検出温度T2actが相対的に遅く増大することに対応している。アッシュ量SASHが大きいほどDPF53の熱容量が増大するからである。係る観点に基づき、第2温度センサ目標検出温度演算部B6は、実際の空気量GA及びアッシュ量SASH量を図6のBに示したテーブルに適用してむだ時間L2を取得し、そのむだ時間L2を上記1式のむだ時間Lに代入する。更に、第2温度センサ目標検出温度演算部B6は、実際の空気量GA及びアッシュ量SASH量を図7のBに示したテーブルに適用して固有周波数ωn2を取得し、その固有周波数ωn2を上記1式の固有周波数ωnに代入する。以上により、第2温度センサ目標検出温度T2tgtが算出される。この第1温度センサ目標検出温度T2tgtは、DPFの再生開始条件成立時点からの時間経過に応じて変化する第2温度センサ検出温度T2actの目標値である。再び図2を参照すると、温度変換部B7は、第2温度センサ66の出力値VT2を関数fを用いて第2温度センサ検出温度T2actに変換する。関数fによれば、出力値VT2が大きいほど第2温度センサ検出温度T2actは高い値となる。第2フィードバックコントローラB8は、第2温度センサ検出温度T2actが第2温度センサ目標検出温度T2tgtに一致するように、ポスト燃料噴射量QPの第2フィードバック量KFB2をPI比例積分制御に則って算出する。即ち、第2フィードバックコントローラB8は、下記の3式に従って第2フィードバック量KFB2を算出する。3式において、ΔT2は第2温度センサ目標検出温度T2tgtから第2温度センサ検出温度T2actを減じた値であるΔT2=第2温度センサ目標検出温度T2tgt-第2温度センサ検出温度T2act。上記3式における積分ゲイン積分定数Ki2は、図9のBに示したように、空気量GAが大きいほど大きくなるように、且つ、アッシュ量SASHが大きくなるほど小さくなるように設定される。同様に、上記3式における比例ゲイン比例定数Kp2は、図10のBに示したように、空気量GAが大きいほど大きくなるように、且つ、アッシュ量SASHが大きくなるほど小さくなるように設定される。積分ゲイン及び比例ゲイン等のフィードバック制御における制御ゲインフィードバック制御の制御パラメータがアッシュ量SASHが大きいほど小さい値に設定されるのは、アッシュ量SASHが大きいほどDPF53の熱容量が大きくなることに起因してフィードバック制御が発散することを回避するためである。但し、これらの制御ゲインは、アッシュ量SASHに対して適切になるように定められればよく、アッシュ量SASHが大きくなるほど大きくなるように設定されてもよい。以上のようにして算出される第1フィードバック量KFB1及び第2フィードバック量KFB2は図2のフィードバック量決定部B9に入力される。フィードバック量決定部B9は、それらの何れか一方を最終的なフィードバック量KFBとして採用して出力する。より具体的に述べると、空気量GAが閾値GAth以下である場合、フィードバック量決定部B9は第1フィードバック量KFB1を最終的なフィードバック量KFBとして採用する。これに対し、空気量GAが閾値GAthよりも大きい場合、フィードバック量決定部B9は第2フィードバック量KFB2を最終的なフィードバック量KFBとして採用する。この理由について次に述べる。図12は、第1温度センサ検出温度T1actが所定温度例えば、500°Cであるときの空気量GAに対する、実DPF温度Tactと第1温度センサ検出温度T1actとの差乖離温度を示すグラフである。実線C1は実DPF温度Tactと第1温度センサ検出温度T1actと差の設計上の値を示す。破線C2は、実DPF温度Tactと第1温度センサ検出温度T1actとの差のばらつき上限値を示し、破線C3はその差のばらつき下限値を示す。これによれば、空気量GAが閾値GAth以下である場合、ばらつき上限値及びばらつき下限値とも設計上の値に近い。これに対し、空気量GAが閾値GAthよりも大きい場合、ばらつき下限値は設計上の値に近いが、ばらつき上限値は設計上の値から大きく乖離している。この理由は、空気量GAが大きい場合にはDOC52内にて燃焼できない未燃成分が増大し、その未燃成分がDPF53内にて燃焼して発熱するからである。図13は、第2温度センサ検出温度T2actが所定温度例えば、500°Cであるときの空気量GAに対する、実DPF温度Tactと第2温度センサ検出温度T2actとの差乖離温度を示すグラフである。実線C1は、実DPF温度Tactと第2温度センサ検出温度T2actとの差の設計上の値を示す。破線C2は、実DPF温度Tactと第2温度センサ検出温度T2actとの差のばらつき上限値を示し、破線C3はその差のばらつき下限値を示す。図13から明らかなように、空気量GAに依らず、ばらつき上限値及びばらつき下限値とも設計上の値に近い。これに対し、図6のA及びBから理解できるように、第1温度センサ65のむだ時間L1は第2温度センサ66のむだ時間L2よりも小さい。但し、第2温度センサ66のむだ時間L2は、空気量GAが大きくなるほど小さくなる。更に、図7のA及びBから理解できるように、第1温度センサ65の固有周波数ωn1は第2温度センサ66の固有周波数ωn2よりも大きい。但し、第2温度センサ66の固有周波数ωn2は、空気量GAが大きくなるほど大きくなる。以上から、空気量GAが閾値空気量GAth以下であれば、第1温度センサ検出温度T1actは実DPF温度Tactと大きく乖離せず且つむだ時間及び固有周波数の点で第2温度センサ検出温度T2actよりもフィードバック制御上有利である。一方、空気量GAが閾値空気量GAthよりも大きくなると、第1温度センサ検出温度T1actは実DPF温度Tactと大きく乖離する可能性があるのに対し、第2温度センサ検出温度T2actは実DPF温度Tactと大きく乖離する可能性がない。更に、空気量GAが閾値空気量GAth以上であれば、むだ時間L2及び固有周波数ωn2がフィードバック制御上問題となるレベルとはならない。そこで、フィードバック量決定部B9は、空気量GAが閾値GAth以下である場合には第1フィードバック量KFB1を最終的なフィードバック量KFBとして採用し、空気量GAが閾値GAthよりも大きい場合には第2フィードバック量KFB2を最終的なフィードバック量KFBとして採用する。再び図2を参照すると、フィードバック量決定部B9から出力されたフィードバック量KFBは基本ポスト噴射量QPbaseに乗じられ、その積QPbaseKFBが最終的なポスト燃料噴射量QPとして採用される。即ち、そのポスト燃料噴射量QPの燃料が、図3に示した膨張下死点近傍のポスト噴射期間において燃料噴射弁31から噴射される。ところで、フィードバック制御禁止許可部B10は、フィードバック制御条件が成立しているか否かを判定している。従って、フィードバック制御禁止許可部B10は、DPFの再生開始条件成立時から前述の第1特定時点までの期間においてフィードバック制御条件は成立していないと判定する。更に、フィードバック制御禁止許可部B10は、DPFの再生開始条件成立時から前述の第2特定時点までの期間においてフィードバック制御条件が成立していないと判定する。フィードバック制御禁止許可部B10は、フィードバック制御条件が成立していないと判定しているとき、フィードバック量決定部B9から出力されるフィードバック量KFBを基本ポスト噴射量QPbaseに乗じることを禁止する。換言すると、フィードバック制御禁止許可部B10は、フィードバック制御条件が成立していないとき、フィードバック量決定部B9から出力されるフィードバック量KFBを1或いは、後述する学習値に置換する。フィードバック制御禁止許可部B10は、フィードバック制御条件が成立していると判定しているとき、フィードバック量決定部B9から出力されるフィードバック量KFBを基本ポスト噴射量QPbaseに乗じることにより最終的なポスト燃料噴射量QPを算出することを許容する。なお、フィードバック制御条件については、後に図19を参照しながら詳述される。以上が、再生制御におけるフィードバック制御の詳細である。実際の作動次に、排気浄化装置10の実際の作動について説明する。CPUは、所定時間が経過する毎に図14にフローチャートにより示したDPF再生制御ルーチンを実行するようになっている。従って、所定のタイミングになると、CPUは図14のステップ1400から処理を開始してステップ1402に進み、再生制御実行フラグXDPFの値が1であるか否かを判定する。この再生制御実行フラグXDPFの値は、図示しない車両のイグニッションキースイッチがオフ位置からオン位置へと変更されたときに実行されるイニシャルルーチンにおいて0に設定される。更に、再生制御実行フラグXDPFの値は、後述するようにDPF再生制御が実行されている期間において1に設定される図14のステップ1406及びステップ1440を参照。。いま、DPF再生制御が実行されておらず、且つ、DPF再生制御開始条件が成立していないと仮定する。この場合、CPUはステップ1402にてNoと判定してステップ1404に進み、再生制御開始条件が成立しているか否かを判定する。前述した仮定によれば、再生制御開始条件は成立していない。よって、CPUはステップ1404にてNoと判定し、ステップ1495に直接進んで本ルーチンを一旦終了する。従って、ポスト噴射は実行されない。次に、係る状態において、DPF再生制御開始条件が成立したと仮定する。この場合、CPUはステップ1402にてNoと判定した後、ステップ1404にてYesと判定し、以下に述べるステップ1406乃至ステップ1420の処理を順に行う。ステップ1406:CPUは、再生制御実行フラグXDPFの値を1に設定する。ステップ1408:CPUは、現時点がDPF再生制御開始条件が成立した直後の時点であるか否かを判定する。現時点は、DPF再生制御開始条件が成立した直後の時点であるので、CPUはステップ1408にてYesと判定してステップ1410に進み、上述したようにPM堆積量等に基づいてDPF指令温度Tshireiを決定する。ステップ1412:CPUは、図19に示したFB条件判定ルーチンにより別途算出されている第1学習値KFB1G及び第2学習値KFB2Gを第1フィードバック量KFB1及び第2フィードバック量KFB2としてそれぞれ格納する。次いで、CPUはステップ1414に進む。なお、現時点がDPF再生制御開始条件が成立した直後の時点でない場合、CPUはステップ1408にてNoと判定し、ステップ1414に直接進む。ステップ1414:CPUは、DPF指令温度Tshirei、機関回転速度NE及びメイン燃料噴射量Qfinに基づいて基本ポスト噴射量QPbaseを決定する。なお、DPF指令温度Tshireiは一定値であってもよい。ステップ1416:CPUはアッシュ量SASHを取得する。CPUは図示しないルーチンによりアッシュ量SASHを周知の方法に従って別途算出している例えば、特開2010-196498号公報を参照。。簡単に述べると、アッシュ量SASHは機関20の潤滑油エンジンオイルの消費量が大きくなるほど大きくなる。潤滑油の消費量は車両の走行距離が大きいほど大きくなる。そこで、CPUは、車速SPDを積算することにより車両の走行距離を算出し、その走行距離が大きいほどアッシュ量SASHが大きくなるようにアッシュ量SASHを算出する。なお、電気制御装置60に本体部30の図示しないオイルパン内に貯留されている潤滑油の液位オイルレベルを検出する周知のオイルレベルセンサが接続されている場合には、CPUは、機関20の始動前にオイルレベルセンサの出力値を取得し、その出力値が低くなるほどアッシュ量SASHが大きくなるようにアッシュ量SASHを算出してもよい。ステップ1418:CPUは、第1温度センサ目標検出温度T1tgtを上述したように第1温度センサモデルを用いて演算する。より具体的に述べると、CPUは図15に示したルーチンを実行することにより、第1温度センサ目標検出温度T1tgtを算出する。即ち、CPUはステップ1418に進むと、以下に述べるステップ1510乃至ステップ1540の処理を順に行い、その後、図14のステップ1420に進む。ステップ1510:CPUは、図6のAに示したテーブルMapL1GAに実際の空気量GAを適用することにより、むだ時間L1を取得する。ステップ1520:CPUは、図7のAに示したテーブルMapωn1GAに実際の空気量GAを適用することにより、固有周波数ωn1を取得する。ステップ1530:CPUは、下記の4式及び5式により計算中間値wnを算出する。これらの式において、固有周波数ωnには固有周波数ωn1が代入される。また、znは、第1温度センサ最終目標温度T1inであり、DPF指令温度Tshireiからオフセット値Offset1を減じた値T1in=Tshirei-Offset1である。なお、4及び5式において、前回値wn-1及び前々回値wn-2が存在しない場合、これらには互いに同じ適当な初期値例えば、第1温度センサ検出温度T1actが適用される。これらの式は、第1温度センサモデルの二次遅れに関する部分を離散化した式である。なお、演算周期Tは本ルーチンの実行間隔時間である。ステップ1540:CPUは、下記の6式及び7式により値ynを第1温度センサ目標検出温度T1tgtとして算出する。6式及び7式において、むだ時間Lにはむだ時間L1が代入される。なお、6及び7式において、前回値yn-1及び前々回値yn-2が存在しない場合、これらには互いに同じ適当な初期値例えば、第1温度センサ検出温度T1actが適用される。これらの式は、第1温度センサモデルのむだ時間に関する部分を離散化した式である。ステップ1420:CPUは、第2温度センサ目標検出温度T2tgtを上述した第2温度センサモデルを用いて演算する。より具体的に述べると、CPUは図16に示したルーチンを実行することにより、第2温度センサ目標検出温度T2tgtを算出する。即ち、CPUはステップ1420に進むと、以下に述べるステップ1610乃至ステップ1640の処理を順に行い、その後、図14のステップ1422に進む。ステップ1610:CPUは、図6のBに示したテーブルMapL2SASH,GAに実際のアッシュ量SASH及び空気量GAを適用することにより、むだ時間L2を取得する。ステップ1620:CPUは、図7のBに示したテーブルMapωn2SASH,GAに実際のアッシュ量SASH及び空気量GAを適用することにより、固有周波数ωn2を取得する。ステップ1630:CPUは、上記4式及び5式により計算中間値wnを算出する。これらの式において、固有周波数ωnには固有周波数ωn2が代入される。また、znは、第2温度センサ最終目標温度T2inであり、DPF指令温度Tshireiからオフセット値Offset2を減じた値T2in=Tshirei-Offset2である。この結果、これらの式は、第2温度センサモデルの二次遅れに関する部分を離散化した式となる。ステップ1640:CPUは、上記6式及び7式により値ynを第2温度センサ目標検出温度T2tgtとして算出する。これらの式において、むだ時間Lにはむだ時間L2が代入される。この結果、これらの式は、第2温度センサモデルのむだ時間に関する部分を離散化した式となる。その後、CPUは図14のステップ1422に進む。CPUはステップ1422にてフィードバック制御条件成立フラグXFB以下、FB条件成立フラグXFBと称呼する。の値が1であるか否かを判定する。FB条件成立フラグXFBは、その値が1であるときフィードバック制御条件が成立していることを示し、その値が0であるときフィードバック制御条件が成立していないことを示す。FB条件成立フラグXFBの値は、後述する図19に示したFB条件判定ルーチンにより変更される。更に、FB条件成立フラグXFBの値は上述したイニシャルルーチンにおいて0に設定される。このとき、FB条件成立フラグXFBの値が1であると、CPUはステップ1422にてYesと判定し、以下に述べるステップ1424及びステップ1426の処理を順に行い、ステップ1428に進む。ステップ1424:CPUは、第1フィードバック量KFB1を上述したようにPI制御に則って演算する。より具体的に述べると、CPUは図17に示したルーチンを実行することにより、第1フィードバック量KFB1を算出する。即ち、CPUはステップ1424に進むと、以下に述べるステップ1710乃至ステップ1760の処理を順に行い、その後、CPUは図14のステップ1426に進む。ステップ1710:CPUは、第1温度センサ65の出力値VT1を関数fを用いて第1温度センサ検出温度T1actに変換する。ステップ1720:CPUは、第1温度センサ目標検出温度T1tgtから第1温度センサ検出温度T1actを減ずることにより偏差ΔT1を算出する。ステップ1730:CPUは、前回の積算値ST1n-1その時点における積算値にステップ1720にて求めた偏差ΔT1を加えることにより、今回の積算値ST1nを算出する。なお、積算値ST1の初期値は0である。ステップ1740:CPUは、図9のAに示したテーブルMapKi1GAに実際の空気量GAを適用することにより、積分ゲインKi1を取得する。ステップ1750:CPUは、図10のAに示したテーブルMapKp1GAに実際の空気量GAを適用することにより、比例ゲインKp1を取得する。ステップ1760:CPUは、ステップ1760内に記載した式に従って第1フィードバック量KFB1を算出する。この式は、上記2式に対応する式である。この後、CPUは図14のステップ1426に進む。ステップ1426:CPUは、第2フィードバック量KFB2を上述したPI制御に則って演算する。より具体的に述べると、CPUは図18に示したルーチンを実行することにより、第2フィードバック量KFB2を算出する。即ち、CPUはステップ1426に進むと、以下に述べるステップ1810乃至ステップ1860の処理を順に行い、その後、図14のステップ1428に進む。ステップ1810:CPUは、第2温度センサ66の出力値VT2を関数fを用いて第2温度センサ検出温度T2actに変換する。ステップ1820:CPUは、第2温度センサ目標検出温度T2tgtから第2温度センサ検出温度T2actを減ずることにより偏差ΔT2を算出する。ステップ1830:CPUは、前回の積算値ST2n-1その時点における積算値にステップ1820にて求めた偏差ΔT2を加えることにより、今回の積算値ST2nを算出する。なお、積算値ST2の初期値は0である。ステップ1840:CPUは、図9のBに示したテーブルMapKi2SASH,GAに実際のアッシュ量SASH及び空気量GAを適用することにより、積分ゲインKi2を取得する。ステップ1850:CPUは、図10のBに示したテーブルMapKp2SASH,GAに実際のアッシュ量SASH及び空気量GAを適用することにより、比例ゲインKp2を取得する。ステップ1860:CPUは、ステップ1860内に記載した式に従って第2フィードバック量KFB2を算出する。この式は、上記3式に対応する式である。この後、CPUは図14のステップ1428に進む。ステップ1428:CPUは、空気量GAが閾値空気量GAth以下であるか否かを判定する。このとき、空気量GAが閾値空気量GAth以下であれば、CPUはステップ1430に進み、フィードバック量KFBとして第1フィードバック量KFB1を格納する。これに対し、空気量GAが閾値空気量GAthよりも大きければ、CPUはステップ1432に進んでフィードバック量KFBとして第2フィードバック量KFB2を格納する。ステップ1434:CPUは、基本ポスト噴射量QPbaseにフィードバック量KFBを乗じることにより、最終的なポスト燃料噴射量QPを求め、各気筒のポスト噴射期間においてポスト燃料噴射量QPの燃料を燃料噴射弁31から噴射させる。この結果、実DPF温度TactがDPF指令温度Tshireiへと精度良く且つ迅速に近づくので、DPF再生制御が円滑に行われる。一方、CPUがステップ1422の処理を行う時点において、FB条件成立フラグXFBの値が0であると、CPUはそのステップ1422にてNoと判定してステップ1436に進み、フィードバック量KFBを1に設定する。次いで、CPUはステップ1434に進む。以上から明らかなように、FB条件成立フラグXFBの値が0である場合、ステップ1424及びステップ1426の処理は実行されず、ステップ1346の処理によってフィードバック量KFBは1に設定される。よって、第1フィードバック量KFB1を用いる第1フィードバック制御、及び、第2フィードバック量を用いる第2フィードバック制御、の何れも実行されない。換言すると、FB条件成立フラグXFBの値が0であるとき、温度センサ第1温度センサ65及び第2温度センサ66の出力値を用いたポスト燃料噴射量QPのフィードバック制御は禁止される。この状態において、DPF再生制御の終了条件が成立したと仮定する。DPF再生制御の終了条件は、例えば、PM堆積量が制御開始許容量PMthよりも小さい制御終了許可量PMend以下になったときに成立する。或いは、DPF再生制御の終了条件は、DPF再生制御の開始時点から所定時間が経過した時点にて成立する。DPF再生制御の終了条件が成立すると、CPUはステップ1402に続くステップ1438にてYesと判定してステップ1440に進み、再生制御実行フラグXDPFの値を0に設定する。次いで、CPUはステップ1442に進み、ポスト燃料噴射量QPの値を0に設定し、その後、ステップ1495に進んで本ルーチンを一旦終了する。以上により、DPF再生制御が終了する。更に、前述したように、CPUは、所定時間が経過する毎に図19にフローチャートにより示したFB条件判定ルーチンを実行するようになっている。従って、所定のタイミングになると、CPUは図19のステップ1900から処理を開始してステップ1905に進み、再生制御実行フラグXDPFの値が1であるか否かを判定する。このとき、再生制御実行フラグXDPFの値が0であると即ち、DPF再生制御が開始されていないと、CPUはステップ1905にてNoと判定してステップ1910に進み、FB条件成立フラグXFBの値を0に設定する。従って、温度センサ第1温度センサ65及び第2温度センサ66の出力値を用いたポスト燃料噴射量QPのフィードバック制御は禁止される。次いで、CPUはステップ1915に進み、第1フィードバック量KFBの学習値KFB1G第1学習値KFB1Gに、その時点の第1フィードバック量KFB1を格納するとともに、第2フィードバック量KFBの学習値KFB2G第2学習値KFB2Gに、その時点の第2フィードバック量KFB2を格納する。その後、CPUはステップ1995に進み、本ルーチンを一旦終了する。一方、CPUがステップ1905の処理を行う時点において、再生制御実行フラグXDPFの値が1であると、CPUはステップ1905にてYesと判定してステップ1920に進み、DOP52の温度及び/又はDPF53の温度が許容温度以上になる可能性があるか否かを判定する。より具体的に述べると、CPUはステップ1910にて、第1温度センサ検出温度T1actが第1許容温度T1OTth以上であるか否かを判定するとともに、第2温度センサ検出温度T2actが第2許容温度T2OTth以上であるか否かを判定する。このとき、第1温度センサ検出温度T1actが第1許容温度T1OTth以上であるか、及び/又は、第2温度センサ検出温度T2actが第2許容温度T2OTth以上であると、CPUはステップ1920にてYesと判定してステップ1925に進み、FB条件成立フラグXFBの値を0に設定する。従って、温度センサ第1温度センサ65及び第2温度センサ66の出力値を用いたポスト燃料噴射量QPのフィードバック制御は禁止される。次いで、CPUはステップ1930に進み、第1学習値KFB1Gに1を格納するとともに、第2学習値KFB2Gに1を格納する。その後、CPUはステップ1995に進み、本ルーチンを一旦終了する。このように、各学習値に1を設定するのは、DOP52の温度及び/又はDPF53の温度が許容温度以上になる可能性がある場合の各フィードバック量KFB1、KFB2を学習値として保持しても無駄だからである。その後、CPUはステップ1995に進み、本ルーチンを一旦終了する。一方、CPUがステップ1920の処理を行う時点において、第1温度センサ検出温度T1actが第1許容温度T1OTth未満であり且つ第2温度センサ検出温度T2actが第2許容温度T2OTth未満であると、CPUはステップ1920にてNoと判定してステップ1935に進み、触媒失活状態にあるか否かを判定する。触媒失活状態とは、ポスト燃料噴射を行ってもDOC52の温度及び/又はPDF53の温度が上昇してこない状態、即ち、これらが活性を失っている状態を意味する。より具体的に述べると、CPUはステップ1920にて、PDF再生制御の開始時点から所定時間が経過したときの第2温度センサ検出温度T2actが許容下限値よりも低いとき、触媒失活状態であると判定する。CPUは、ステップ1935にて触媒失活状態であると判定すると、ステップ1935からステップ1925及びステップ1930に進む。従って、フィードバック制御は禁止される。また、各フィードバック量KFB1、KFB2は学習値KFB1G、KFB2Gとして保持されない。一方、CPUがステップ1935の処理を行う時点において、触媒失活状態でない場合、CPUはステップ1935にてNoと判定してステップ1940に進み、ポスト燃料噴射量を制限すべき条件が不成立であるか否かを判定する。このポスト燃料噴射量を制限すべき条件は白煙が発生する可能性が高いときに成立する条件であり、具体的には、空気量GAが所定の上限値GAHith未満であればポスト燃料噴射量を制限すべき条件は不成立であると判定され、空気量GAが上限値GAHith以上であればポスト燃料噴射量を制限すべき条件は成立していると判定される。従って、空気量GAが所定の上限値GAHith以上であるとき、CPUはステップ1940にてNoと判定してステップ1910及びステップ1915に進む。この結果、フィードバック制御は禁止される。また、各フィードバック量KFB1、KFB2は学習値KFB1G、KFB2Gとして格納される。一方、CPUがステップ1940の処理を行う時点において、空気量GAが上限値GAHith未満であれば、CPUはそのステップ1940にてYesと判定してステップ1945に進み、空気量GAが所定の下限値GALoth以上であるか否かを判定する。下限値GALothは上限値GAHithよりも当然に小さい。このとき、空気量GAが下限値GALoth未満であると、上述したむだ時間が著しく長くなるとともに、固有周波数ωnが著しく小さくなる。従って、このような場合に上述したフィードバック制御を実行することは好ましくない。そこで、CPUは、空気量GAが下限値GALoth未満であるとき、ステップ1945にてNoと判定してステップ1910及びステップ1915に進む。一方、CPUがステップ1945の処理を行う時点において、空気量GAが下限値GALoth以上であれば、CPUはそのステップ1945にてYesと判定してステップ1950に進み、現時点は機関20の始動後から所定時間以上が経過した時点であるか否かを判定する。このとき、機関20の始動後から所定時間以上が経過していなければ、CPUはステップ1950にてNoと判定してステップ1910及びステップ1915に進む。一方、CPUがステップ1950の処理を行う時点において、機関20の始動後から所定時間以上が経過していれば、CPUはステップ1950にてYesと判定してステップ1955に進み、最上流温度センサ64の出力値VT0に基づく最上流温度検出値T0actが最上流温度下限値T0Loth以上であるか否かを判定する。最上流温度検出値T0actが最上流温度下限値T0Loth未満であるとき、排気管の温度は非常に低いので、ポスト燃料噴射を行った場合に燃料が排気管に多量に付着するため、第1温度センサ検出温度T1act及び第2温度センサ検出温度T2actが上昇せず、その結果、第1及び第2フィードバック量が過大になる恐れがある。更に、排気管に付着した燃料がある時点にて大量にDOC52に流入する恐れがある。そこで、CPUは、最上流温度検出値T0actが最上流温度下限値T0Loth未満であるとき、ステップ1955にてNoと判定してステップ1910及びステップ1915に進み、フィードバック制御を禁止する。一方、CPUがステップ1955の処理を行う時点において、最上流温度検出値T0actが最上流温度下限値T0Loth以上であれば、CPUはそのステップ1955にてYesと判定してステップ1060に進み、第1温度センサ検出温度T1actが所定の温度閾値T1Hith以上であるか否かを判定する。第1温度センサ検出温度T1actが温度閾値T1Hith未満であると、フィードバック制御第1フィードバック制御開始直後における第1温度センサ検出温度T1actと第1温度センサ最終目標温度T1inとの差の大きさが非常に大きくなる。このような場合に、第1フィードバック制御を行うことは好ましくない。そこで、CPUは、第1温度センサ検出温度T1actが温度閾値T1Hith未満であるとき、ステップ1960にてNoと判定してステップ1910及びステップ1915に進み、フィードバック制御を禁止する。一方、CPUがステップ1960の処理を行う時点において、第1温度センサ検出温度T1actが温度閾値T1Hith以上であれば、CPUはそのステップ1960にてYesと判定してステップ1965に進み、下記の第1温度関連条件及び第2温度関連条件の双方が成立しているか否かを判定する。第1温度関連条件第1温度センサ目標検出温度T1tgtと第1温度センサ検出温度T1actとの差ΔT1=T1in-T1tgtの大きさが、第1の閾値温度差ΔT1thより小さい。第2温度関連条件第2温度センサ目標検出温度T2tgtと第2温度センサ検出温度T2actとの差ΔT2=T2in-T2tgtの大きさが、第2の閾値温度差ΔT2thより小さい。図20を参照しながら説明したように、第1温度関連条件が満足されない場合、第1温度センサモデルにより算出される第1温度センサ目標検出温度T1tgtの精度は良好でない。同様に、第2温度関連条件が満足されない場合、第2温度センサモデルにより算出される第2温度センサ目標検出温度T2tgtの精度は良好でない。そこで、CPUは、第1温度関連条件及び第2温度関連条件の少なくとも一方が不成立である場合、ステップ1965にてNoと判定してステップ1910及びステップ11915に進む。この結果、温度センサ第1温度センサ65及び第2温度センサ66の出力値を用いたポスト燃料噴射量QPのフィードバック制御は禁止される。一方、CPUがステップ1965の処理を行う時点において、第1温度関連条件及び第2温度関連条件の双方が成立していると、CPUはステップ1965にてYesと判定してステップ1970に進み、FB条件成立フラグXFBの値を1に設定する。この結果、図14のステップ1424乃至ステップ1432のうちの適当なステップの処理が実行されるので、温度センサ第1温度センサ65及び第2温度センサ66の出力値を用いたポスト燃料噴射量QPのフィードバック制御が許可及び実行される。なお、DPFの再生開始条件成立時点以降において、第1センサモデルが変動する振幅ΔW1図20の点P1の温度と及び点P3の温度との差は、むだ時間L1と、第1温度センサ最終目標温度T1inと再生開始条件成立時点における第1温度センサ検出温度T1actとの差SA1の大きさと、に応じて変化することが判明している。従って、この振幅ΔW1を、空気量GAと図6のAのテーブルとから取得されるむだ時間L1と、差SA1と、に基づいて求め、更に、その振幅ΔW1の半分ΔW1/2を第1温度センサ最終目標温度T1inから減ずることにより得られる値又は、その値以下の値を、第1の閾値温度差ΔT1thとして取得してもよい。同様に、DPFの再生開始条件成立時点以降において、第2センサモデルが変動する振幅ΔW2は、むだ時間L2と、第2温度センサ最終目標温度T2inと再生開始条件成立時点における第2温度センサ検出温度T2actとの差SA2の大きさと、に応じて変化することが判明している。従って、この振幅ΔW2を、空気量GA及びアッシュ量SASHと図6のBのテーブルとから取得されるむだ時間L2と、差SA2と、に基づいて求め、更に、その振幅ΔW2の半分ΔW2/2を第2温度センサ最終目標温度T2inから減ずることにより得られる値又は、その値以下の値を、第2の閾値温度差ΔT2thとして取得してもよい。更に、CPUは、DPFの再生開始条件成立時点以降において第1温度センサモデルにより算出される第1温度センサ目標検出温度T1tgtを監視し、再生開始条件成立時点以降に最初に現れる第1温度センサ目標検出温度T1tgtの極大値図20の点P1を参照。を求め、その極大値と第1温度センサ最終目標温度T1inとの差の大きさまたはその差の大きさに基づく値であって且つその差以下の値を第1の閾値温度差ΔT1thとして採用してもよい。同様に、CPUは、DPFの再生開始条件成立時点以降において第2温度センサモデルにより算出される第2温度センサ目標検出温度T2tgtを監視し、再生開始条件成立時点以降に最初に現れる第2温度センサ目標検出温度T2tgtの極大値を求め、その極大値と第2温度センサ最終目標温度T2inとの差の大きさまたはその差の大きさに基づく値であって且つその差以下の値を第2の閾値温度差ΔT2thとして採用してもよい。以上、説明したように、排気浄化装置10のDPF53を再生させる再生制御手段は、前記再生開始条件成立時点以降において前記パティキュレートフィルタの温度を前記指令温度に一致させるように予め定められたフィードフォワード量のエネルギーを前記パティキュレートフィルタに与えるフィードフォワード手段図2のB2、図14のステップ1414を参照。と、前記再生開始条件成立時点以降において前記センサ検出温度が前記指令温度に基づいて定まる最終的な目標温度である最終目標温度と一致するように前記エネルギーの量をフィードバック制御するフィードバック手段図2のB1、B3、B5、B6、B8、B9、図14のステップ1410、ステップ1412、ステップ1418乃至ステップ1434を参照。と、を備え、更に、前記フィードバック手段は、前記再生開始条件成立時点から前記フィードフォワード量のエネルギーを前記パティキュレートフィルタに与えた場合における同再生開始条件成立時点からの経過時間に対する前記センサ検出温度を前記最終目標温度に基づいて推定するとともに、その推定したセンサ検出温度を前記フィードバック制御における目標値であるセンサ目標検出温度として発生する温度センサモデル第1温度センサモデル、第2温度センサモデルを用いることにより、前記再生開始条件成立時点からの時間経過に応じた前記センサ目標検出温度を算出する目標検出温度算出手段図2のB3、B6、及び、図14のステップ1418、ステップ1420、図15、及び、図16等を参照。と、前記センサ目標検出温度と前記センサ検出温度とが一致するように前記パティキュレートフィルタに加えられるエネルギーの量を制御することによって前記フィードバック制御を実行するフィードバック制御手段図2のB5,B6,B9、及び、図14のステップ1424、ステップ1426、図17及び図18等を参照。と、前記再生開始条件成立時点から所定期間が経過する特定時点第1温度関連条件及び第2温度関連条件が成立する時点までは前記センサ目標検出温度を用いた前記フィードバック制御の実行を禁止し、前記特定時点以降において前記フィードバック制御の実行を許可するフィードバック制御禁止許可手段図2のB10、図14のステップ1422及びステップ1436、並びに、図19のステップ1965、ステップ1970及びステップ1910を参照。と、を含む。従って、センサモデルにより算出されるセンサ目標検出温度の推定精度が良好でない期間においてそのセンサ目標検出温度に基づくポスト燃料噴射量のフィードバック制御が禁止され、その後、センサ目標検出温度の推定精度が良好となった時点以降においてセンサ目標検出温度に基づくポスト燃料噴射量のフィードバック制御が実行されるので、DPF53の温度を狙いに近い値へと制御することができる。本発明は上記実施形態に限定されることはなく、本発明の範囲内において種々の変形例を採用することができる。例えば、上記第2センサモデルは、むだ時間及び固有周波数を空気量GA及びアッシュ量SASHに基づいて変更していたが、それらを空気量GAのみに基づいて変更してもよい。更に、第2フィードバック量を算出するPI制御において、フィードバック制御定数積分ゲインKi2及び比例ゲインKp2等が空気量GA及びアッシュ量SASHに基づいて変更されていたが、そのフィードバック制御定数は空気量GAのみに基づいて変更されてもよい。更に、DPF再生制御においてDPF53に加えられるエネルギー熱エネルギーは、ポスト燃料噴射量により変更されるだけではなく、メイン燃料噴射量の噴射時期を遅角すること、吸気絞り弁42の開度を変更すること、機関20がEGRシステムを備える場合にはEGR量を変更すること、及び、機関20がDOC52の直前に燃料添加弁を備える場合にはその燃料添加弁から供給される添加燃料量を変更すること、の何れか一つ又は一つ以上の組み合わせによってなされてもよい。更に、第1温度センサモデル及び第2温度センサモデルは、Padeの近似式二次により構築されていたが、むだ時間と遅れを表現する近似式例えば、ラゲールLaguerre近似式、又は、Katz近似式等に基づいて構築されてもよい。更に、DPF53は酸化触媒を担持していてもよい。この場合、DOC52を省略してもよい。また、上記実施形態においては、第1フィードバック量KFB1と、第2フィードバック量KFB2と、の何れかを最終的なフィードバック量KFBとして採用していたが、第1フィードバック量KFB1及び第2フィードバック量KFB2の何れか一方のみを最終的なフィードバック量KFBとして用いてもよい。換言すると、第1温度センサ65及び第2温度センサ66の何れか一方のみを使用してポスト燃料噴射量のフィードバック制御を行ってもよい。更に、上記実施形態におけるフィードバック制御を許可する特定時点を、DPFの再生制御開始時点DPF制御開始条件成立時点から所定時間が経過した時点として特定してもよい。この場合、所定時間は、むだ時間L1,L2に応じて定められてもよい。更に、再生制御開始時点DPF制御開始条件成立時点以降において、第1温度センサ目標検出温度T1tgtを監視し、その第1温度センサ目標検出温度T1tgtが極大値図20の点P1を通り、次に、極小値点P3を通ったことが確認された後に再生制御開始時点の第1温度センサ検出温度T1act=T0に戻った時点を、第1温度センサに係る上記特定時点として検出してもよい。同様に、再生制御開始時点DPF制御開始条件成立時点以降において、第2温度センサ目標検出温度T2tgtを監視し、その第2温度センサ目標検出温度T2tgtが極大値を通り、次に、極小値を通ったことが確認された後に再生制御開始時点の第2温度センサ検出温度T2actに戻った時点を、第2温度センサに係る上記特定時点として検出してもよい。更に、上記実施形態においては、第1温度関連条件及び第2温度関連条件の双方が成立したときに、第1フィードバック量KFB1及び第2フィードバック量KFB2の算出を開始していたが、第1温度関連条件が成立したときに第1フィードバック量KFB1の算出を開始し、その時点にて第2温度関連条件が成立していなければ、第1フィードバック量KFB1を最終的なフィードバック量KFBとして採用してポスト燃料噴射量のフィードバック制御を開始してもよい。同様に、第2温度関連条件が成立したときに第2フィードバック量KFB2の算出を開始し、その時点にて第1温度関連条件が成立していなければ、第2フィードバック量KFB2を最終的なフィードバック量KFBとして採用してポスト燃料噴射量のフィードバック制御を開始してもよい。"