"複数気筒を有する内燃機関の排気通路において空燃比を検出する空燃比検出手段と、前記複数気筒のそれぞれに供給する燃料量を制御する燃料量制御手段とを備える、内燃機関の空燃比制御装置において、前記機関の運転状態を示す少なくとも一つの運転状態パラメータを取得する運転状態パラメータ取得手段と、故障判定期間中に前記空燃比検出手段の検出信号から特定周波数成分を抽出する抽出手段と、抽出された特定周波数成分に基づいて前記機関の空燃比制御系の故障判定を行う故障判定手段と、前記故障判定期間開始後において前記運転状態パラメータの変動状態を示す変動状態パラメータであって、前記運転状態パラメータの変動履歴が反映される変動状態パラメータを算出する変動状態パラメータ算出手段と、前記変動状態パラメータが所定閾値以上である特定変動状態が検出されたときに、前記故障判定を中断または停止する判定停止手段とを備え、前記特定変動状態は、前記抽出手段により抽出される特定周波数成分の強度に影響を与える変動状態であることを特徴とする内燃機関の空燃比制御装置。前記変動状態パラメータ算出手段は、前記故障判定期間開始後において前記運転状態パラメータまたは前記運転状態パラメータの平均値が所定上下限値の範囲内にあるときに、前記変動状態パラメータの算出を実行することを特徴とする請求項1に記載の内燃機関の空燃比制御装置。前記運転状態パラメータ取得手段は、複数の運転状態パラメータを取得し、前記変動状態パラメータ算出手段は、前記変動状態パラメータの算出に適用される運転状態パラメータとは異なる他の運転状態パラメータまたは該他の運転状態パラメータの平均値が所定上下限値の範囲内にあるときに、前記変動状態パラメータの算出を実行することを特徴とする請求項1に記載の内燃機関の空燃比制御装置。前記変動状態パラメータ算出手段は、前記故障判定期間開始後において前記運転状態パラメータの変化量絶対値を積算することにより前記変動状態パラメータを算出することを特徴とする請求項1から3の何れか1項に記載の内燃機関の空燃比制御装置。前記変動状態パラメータ算出手段は、前記故障判定期間開始後において前記運転状態パラメータの正の変化量を積算することにより増加量積算値を算出するとともに、前記運転状態パラメータの負の変化量の絶対値を積算することにより減少量積算値を算出し、前記増加量積算値及び減少量積算値の少なくとも一方を、前記変動状態パラメータとして算出することを特徴とする請求項1から3の何れか1項に記載の内燃機関の空燃比制御装置。前記変動状態パラメータ算出手段は、前記故障判定期間開始後における前記運転状態パラメータの最大値及び最小値を更新するとともに前記最大値と前記最小値との差分値を前記変動状態パラメータとして算出することを特徴とする請求項1から3の何れか1項に記載の内燃機関の空燃比制御装置。前記変動状態パラメータ算出手段は、前記故障判定期間開始後における前記運転状態パラメータに含まれる所定周波数成分を抽出するバンドパスフィルタ処理を実行し、該バンドパスフィルタ処理後の運転状態パラメータを前記変動状態パラメータとして算出することを特徴とする請求項1から3の何れか1項に記載の内燃機関の空燃比制御装置。前記特定周波数成分は、前記機関の回転速度に対応する周波数の1/2の周波数成分である0.5次周波数成分であり、前記故障判定手段は、前記0.5次周波数成分に基づいて、前記複数気筒のそれぞれに対応する空燃比が許容限度を超えてばらついているインバランス故障を判定することを特徴とする請求項1から7の何れか1項に記載の内燃機関の空燃比制御装置。前記空燃比を設定周波数で変動させる空燃比変動手段をさらに備え、前記特定周波数成分は、前記設定周波数の成分であり、前記故障判定手段は、前記設定周波数の成分に基づいて、前記空燃比検出手段の劣化故障を判定することを特徴とする請求項1から7の何れか1項に記載の内燃機関の空燃比制御装置。上記従来の装置では、機関の運転状態パラメータXOP吸入空気量、機関回転数の変化量を用いた判定は、一定時間DT例えば100msec毎の変化量DXOPと、所定量DXOPTHとを比較することにより行われるため、例えば300msec程度の期間では大きく変化しているような運転状態であっても、変化量DXOPは所定量DXOPTH以下であるような場合には、データ取得が許可される図14a参照。また運転状態パラメータXOPが大きく変動していても、所定上下限値XOPLMH,XOPLMLの範囲内の変動であれば、データ取得が許可される図14b参照。したがって、上記従来の装置で採用されるデータ取得条件では、運転状態パラメータXOPが大きく変動している状態において取得したデータに基づく判定が行われ、判定精度が低下するおそれがある。特に上記1サイクル周波数に近い周波数の変動が含まれているときに、判定精度が低下する可能性が高くなる。本発明はこの点に着目してなされたものであり、空燃比制御系の故障判定の開始後における機関運転状態の変化を的確に監視し、故障判定精度を向上させることができる内燃機関の空燃比制御装置を提供することを目的とする。上記目的を達成するため請求項1に記載の発明は、複数気筒を有する内燃機関の排気通路において空燃比KACTを検出する空燃比検出手段と、前記複数気筒のそれぞれに供給する燃料量TOUTを制御する燃料量制御手段とを備える、内燃機関の空燃比制御装置において、前記機関の運転状態を示す少なくとも一つの運転状態パラメータXOPを取得する運転状態パラメータ取得手段と、故障判定期間中に前記空燃比検出手段の検出信号から特定周波数成分周波数f1成分、0.5次周波数成分を抽出する抽出手段と、抽出された特定周波数成分に基づいて前記機関の空燃比制御系の故障判定を行う故障判定手段と、前記故障判定期間開始後において前記運転状態パラメータXOPの変動状態を示す変動状態パラメータであって、前記運転状態パラメータの変動履歴が反映される変動状態パラメータIDXOP,IDXOPP,IDXOPN,XOPMAX-XOPMIN,XOPBFAを算出する変動状態パラメータ算出手段と、前記変動状態パラメータが所定閾値以上である特定変動状態が検出されたときに、前記故障判定を中断または停止する判定停止手段とを備え、前記特定変動状態は、前記抽出手段により抽出される特定周波数成分の強度に影響を与える変動状態であることを特徴とする。請求項2に記載の発明は、請求項1に記載の内燃機関の空燃比制御装置において、前記変動状態パラメータ算出手段は、前記故障判定期間開始後において前記運転状態パラメータXOPまたは前記運転状態パラメータの平均値XOPAVEが所定上下限値XOPLMH,XOPLMLの範囲内にあるときに前記変動状態パラメータの算出を実行することを特徴とする。請求項3に記載の発明は、請求項1に記載の内燃機関の空燃比制御装置において、前記運転状態パラメータ取得手段は、複数の運転状態パラメータを取得し、前記変動状態パラメータ算出手段は、前記変動状態パラメータの算出に適用される運転状態パラメータとは異なる他の運転状態パラメータGAIRまたは該他の運転状態パラメータの平均値GAIRAVEが所定上下限値の範囲内にあるときに前記変動状態パラメータの算出を実行することを特徴とする。請求項4に記載の発明は、請求項1から3の何れか1項に記載の内燃機関の空燃比制御装置において、前記変動状態パラメータ算出手段は、前記故障判定期間開始後において前記運転状態パラメータの変化量絶対値DXOPAを積算することにより前記変動状態パラメータIDXOPを算出することを特徴とする。請求項5に記載の発明は、請求項1から3の何れか1項に記載の内燃機関の空燃比制御装置において、前記変動状態パラメータ算出手段は、前記故障判定期間開始後において前記運転状態パラメータの正の変化量を積算することにより増加量積算値IDXOPPを算出するとともに、前記運転状態パラメータの負の変化量の絶対値を積算することにより減少量積算値IDXOPNを算出し、前記増加量積算値及び減少量積算値の少なくとも一方を、前記変動状態パラメータとして算出することを特徴とする。請求項6に記載の発明は、請求項1から3の何れか1項に記載の内燃機関の空燃比制御装置において、前記変動状態パラメータ算出手段は、前記故障判定期間開始後における前記運転状態パラメータの最大値XOPMAX及び最小値XOPMINを更新するとともに前記最大値と前記最小値との差分値XOPMAX-XOPMINを前記変動状態パラメータとして算出することを特徴とする。請求項7に記載の発明は、請求項1から3の何れか1項に記載の内燃機関の空燃比制御装置において、前記変動状態パラメータ算出手段は、前記故障判定期間開始後における前記運転状態パラメータに含まれる所定周波数成分を抽出するバンドパスフィルタ処理を実行し、該バンドパスフィルタ処理後の運転状態パラメータXOPBFAを前記変動状態パラメータとして算出することを特徴とする。請求項8に記載の発明は、請求項1から7の何れか1項に記載の内燃機関の空燃比制御装置において、前記特定周波数成分は、前記機関の回転速度に対応する周波数の1/2の周波数成分である0.5次周波数成分であり、前記故障判定手段は、前記0.5次周波数成分に基づいて、前記複数気筒のそれぞれに対応する空燃比が許容限度を超えてばらついているインバランス故障を判定することを特徴とする。請求項9に記載の発明は、請求項1から7の何れか1項に記載の内燃機関の空燃比制御装置において、前記空燃比を設定周波数f1で変動させる空燃比変動手段をさらに備え、前記特定周波数成分は、前記設定周波数の成分周波数f1成分であり、前記故障判定手段は、前記設定周波数の成分に基づいて、前記空燃比検出手段の劣化故障を判定することを特徴とする。請求項1に記載の発明によれば、故障判定期間中に空燃比検出手段の検出信号から特定周波数成分が抽出され、抽出された特定周波数成分に基づいて空燃比制御系の故障判定が行われる。故障判定期間開始後において機関の運転状態を示す運転状態パラメータの変動状態を示し、かつ運転状態パラメータの変動履歴が反映される変動状態パラメータが算出され、この変動状態パラメータが所定閾値以上である特定変動状態、すなわち抽出手段によって抽出される特定周波数成分の強度に影響を与える変動状態が検出されたときに、故障判定が中断または停止される。運転状態パラメータの変動履歴が反映される変動状態パラメータを用いて特定変動状態を検出することにより、運転状態パラメータの変動が特定周波数成分の強度に影響を及ぼすことに起因する誤判定を防止し、故障判定精度を向上させることができる。なお、変動状態パラメータとして、運転状態パラメータの平均値を用いた場合には、上記特定変動状態を検出することはできないので、運転状態パラメータの平均値は、変動状態パラメータには含まれない。請求項2に記載の発明によれば、故障判定期間開始後において運転状態パラメータまたはその平均値が所定上下限値の範囲内にあるときに変動状態パラメータの算出が行われるので、運転状態パラメータが故障判定に適さない値に変化したときは、変動状態パラメータの算出が停止される。これにより、変動状態パラメータによる不適切な停止条件判定が回避され、故障判定精度の低下を防止できる。請求項3に記載の発明によれば、故障判定期間開始後において、変動状態パラメータの算出に適用される運転状態パラメータとは異なる他の運転状態パラメータまたは該他の運転状態パラメータの平均値が、所定上下限値の範囲内にあるときに変動状態パラメータの算出が行われる。他の運転状態パラメータが故障判定に適さない値に変化したときは、変動状態パラメータの算出が停止されるので、変動状態パラメータによる不適切な停止条件判定が回避され、故障判定精度の低下を防止できる。請求項4に記載の発明によれば、故障判定期間開始後において運転状態パラメータの変化量絶対値を積算することにより変動状態パラメータが算出されるので、比較的簡便な演算で、運転状態パラメータの変動履歴を的確に示す変動状態パラメータが得られる。その結果、制御装置の演算負荷を増大させることなく故障判定の実行可否をより適切に判定し、判定精度を向上させることができる。請求項5に記載の発明によれば、故障判定期間開始後において運転状態パラメータの正の変化量を積算することにより増加量積算値が算出されるとともに、運転状態パラメータの負の変化量の絶対値を積算することにより減少量積算値が算出され、増加量積算値及び減少量積算値の少なくとも一方が、変動状態パラメータとして算出される。したがって、比較的簡便な演算で、運転状態パラメータの増加及び/または減少の履歴を的確に示す変動状態パラメータが得られる。その結果、制御装置の演算負荷を増大させることなく故障判定の実行可否をより適切に判定し、判定精度を向上させることができる。請求項6に記載の発明によれば、故障判定期間開始後における運転状態パラメータの最大値及び最小値が更新されるとともに、最大値と最小値との差分値が変動状態パラメータとして算出されるので、比較的簡便な演算で、空燃比検出手段の出力特性が大きく変化する変動履歴を確実に判定することができる。その結果、制御装置の演算負荷を増大させることなく故障判定の実行可否をより適切に判定し、判定精度を向上させることができる。請求項7に記載の発明によれば、故障判定期間開始後における運転状態パラメータに含まれる所定周波数成分を抽出するバンドパスフィルタ処理が実行され、該バンドパスフィルタ処理後の運転状態パラメータが変動状態パラメータとして算出されるので、故障判定に使用される特定周波数成分に大きな影響を与える変動周波数成分が運転状態パラメータに含まれている状態を、比較的簡便な演算で確実に判定することができる。その結果、制御装置の演算負荷を増大させることなく故障判定の実行可否をより適切に判定し、判定精度を向上させることができる。請求項8に記載の発明によれば、特定周波数成分は、機関の回転速度に対応する周波数の1/2の周波数成分である0.5次周波数成分であり、0.5次周波数成分に基づいて、複数気筒のそれぞれに対応する空燃比が許容限度を超えてばらついているインバランス故障が判定される。したがって、運転状態パラメータの変動によってインバランス故障の判定精度が低下することを防止できる。請求項9に記載の発明によれば、空燃比を設定周波数で変動させる空燃比変動制御が行われ、設定周波数の成分に基づいて、空燃比検出手段の劣化故障が判定される。したがって、運転状態パラメータの変動によって空燃比検出手段の劣化故障の判定精度が低下することを防止できる。以下本発明の実施の形態を図面を参照して説明する。[第1の実施形態]図1は、本発明の一実施形態にかかる内燃機関以下エンジンという及びその空燃比制御装置の全体構成図であり、例えば4気筒のエンジン1の吸気管2の途中にはスロットル弁3が配されている。スロットル弁3にはスロットル弁開度THを検出するスロットル弁開度センサ4が連結されており、その検出信号は電子制御ユニット以下ECUという5に供給される。燃料噴射弁6はエンジン1とスロットル弁3との間かつ吸気管2の図示しない吸気弁の少し上流側に各気筒毎に設けられており、各噴射弁は図示しない燃料ポンプに接続されていると共にECU5に電気的に接続されて当該ECU5からの信号により燃料噴射弁6の開弁時間が制御される。スロットル弁3の上流側には吸入空気流量GAIRを検出する吸入空気流量センサ7が設けられている。またスロットル弁3の下流側には吸気圧PBAを検出する吸気圧センサ8、及び吸気温TAを検出する吸気温センサ9が設けられている。これらのセンサの検出信号は、ECU5に供給される。エンジン1の本体には、エンジン冷却水温TWを検出する冷却水温センサ10が装着されており、その検出信号はECU5に供給される。ECU5には、エンジン1のクランク軸図示せずの回転角度を検出するクランク角度位置センサ11が接続されており、クランク軸の回転角度に応じた信号がECU5に供給される。クランク角度位置センサ11は、エンジン1の特定の気筒の所定クランク角度位置でパルス以下CYLパルスというを出力する気筒判別センサ、各気筒の吸入行程開始時の上死点TDCに関し所定クランク角度前のクランク角度位置で4気筒エンジンではクランク角180度毎にTDCパルスを出力するTDCセンサ及びTDCパルスより短い一定クランク角周期例えば6度周期で1パルス以下CRKパルスというを発生するCRKセンサから成り、CYLパルス、TDCパルス及びCRKパルスがECU5に供給される。これらのパルスは、燃料噴射時期、点火時期等の各種タイミング制御、エンジン回転数エンジン回転速度NEの検出に使用される。排気通路13には三元触媒14が設けられている。三元触媒14は、酸素蓄積能力を有し、エンジン1に供給される混合気の空燃比が理論空燃比よりリーン側に設定され、排気中の酸素濃度が比較的高い排気リーン状態では、排気中の酸素を蓄積し、逆にエンジン1に供給される混合気の空燃比が理論空燃比よりリッチ側に設定され、排気中の酸素濃度が低く、HC、CO成分が多い排気リッチ状態では、蓄積した酸素により排気中のHC,COを酸化する機能を有する。三元触媒14の上流側であって各気筒に連通する排気マニホールドの集合部より下流側には、比例型酸素濃度センサ15以下LAFセンサ15というが装着されており、このLAFセンサ15は排気中の酸素濃度空燃比にほぼ比例した検出信号を出力し、ECU5に供給する。ECU5には、エンジン1により駆動される車両のアクセルペダルの踏み込み量以下アクセルペダル操作量というAPを検出するアクセルセンサ21及び当該車両の走行速度車速VPを検出する車速センサ22が接続されており、それらセンサの検出信号がECU5に供給される。スロットル弁3は図示しないアクチュエータにより開閉駆動され、スロットル弁開度THはアクセルペダル操作量APに応じてECU5により制御される。なお、図示は省略しているが、エンジン1には周知の排気還流機構が設けられている。ECU5は、各種センサからの入力信号波形を整形し、電圧レベルを所定レベルに修正し、アナログ信号値をデジタル信号値に変換する等の機能を有する入力回路、中央演算処理ユニット以下CPUという、該CPUで実行される各種演算プログラム及び演算結果等を記憶する記憶回路、燃料噴射弁6に駆動信号を供給する出力回路を備えている。ECU5のCPUは、上述の各種センサの検出信号に基づいて、種々のエンジン運転状態を判別するとともに、該判別されたエンジン運転状態に応じて、次式1を用いて、TDCパルスに同期して開弁作動する燃料噴射弁6の燃料噴射時間TOUTを演算する。燃料噴射時間TOUTは、噴射される燃料量にほぼ比例するので、以下燃料噴射量TOUTという。TOUT=TIM×KCMD×KAF×KTOTAL1ここに、TIMは基本燃料量、具体的には燃料噴射弁6の基本燃料噴射時間であり、吸入空気流量GAIRに応じて設定されたTIMテーブルを検索して決定される。TIMテーブルは、エンジンにおいて燃焼する混合気の空燃比AFがほぼ理論空燃比になるように設定されている。KCMDはエンジン1の運転状態に応じて設定される目標空燃比係数である。目標空燃比係数KCMDは、空燃比A/Fの逆数、すなわち燃空比F/Aに比例し、理論空燃比のとき値1.0をとるので、以下目標当量比という。後述するように、LAFセンサ15の応答特性劣化故障判定を行うときは、1.0±DAFの範囲で時間経過に伴って正弦波状に変化するように設定される。KAFは、空燃比フィードバック制御の実行条件が成立するときは、LAFセンサ15の検出値から算出される検出当量比KACTが目標当量比KCMDに一致するようにPID比例積分微分制御あるいは適応制御器SelfTuningRegulatorを用いた適応制御により算出される空燃比補正係数である。KTOTALは夫々各種エンジンパラメータ信号に応じて演算される他の補正係数エンジン冷却水温TMに応じた補正係数KTW、吸気温TAに応じた補正係数KTAなどの積である。ECU5のCPUは上述のようにして求めた燃料噴射量TOUTに基づいて燃料噴射弁6を開弁させる駆動信号を出力回路を介して燃料噴射弁6に供給する。また、ECU5のCPUは、以下に説明するようにLAFセンサ15の応答特性劣化故障判定を行う。本実施形態における応答特性劣化故障判定は、基本的には特開2010-101289号公報に示される手法と同一のものであり、エンジン運転中に空燃比を周波数f1で振動させる空燃比振動制御を実行し、LAFセンサ15の出力信号から算出される検出当量比KACTに含まれる周波数f1の成分強度MPTf1と、周波数f1の2倍の周波数である周波数f2に対応する周波数f2の成分強度MPTf2とを用いて、応答特性劣化故障が判定される。図2は、故障判定処理の全体構成を示すフローチャートである。この処理は、所定クランク角度CACAL例えば30度毎にECU5のCPUで実行される。ステップS1では、実行条件フラグFMCNDが1であるか否かを判別する。実行条件フラグFMCNDは、実行条件判定処理図示せずにおいて故障判定の実行条件が成立するとき1に設定される。具体的には下記の条件1~11がすべて満たされると、実行条件フラグFMCNDが1に設定される。条件1~11の何れかが満たされないときは、実行条件フラグFMCNDは0に維持される。1エンジン回転数NEが所定上下限値の範囲内にある。2吸気圧PBAが所定圧より高い判定に必要な排気流量が確保されている。3LAFセンサ15が活性化している。4LAFセンサ15の出力に応じた空燃比フィードバック制御が実行されている。5エンジン冷却水温TWが所定温度より高い。6エンジン回転数NEの単位時間当たりの変化量DNEが所定回転数変化量より小さい。7吸気圧PBAの単位時間当たりの変化量DPBAFが所定吸気圧変化量より小さい。8燃料の加速増量急加速時に実行されるが行われていない。9排気還流率が所定値より大きい。10LAFセンサ出力が上限値または下限値に張り付いた状態ではない。11LAFセンサの応答特性が正常である応答特性の劣化故障が発生しているとの判定が行われていない。ステップS1で実行条件フラグFMCNDが0であるときは直ちに処理を終了する。実行条件フラグFMCNDが1に設定されると、ステップS1からステップS2に進み、目標当量比KCMDを下記式2によって振動させる空燃比振動制御を実行する。空燃比振動制御実行中は、空燃比補正係数KAFを1.0または1.0以外の特定の値に固定される。式2のKf1は、例えば振動周波数f1を0.4fNEfNE=NE[rpm]/60とする場合には、0.4に設定される第1周波数係数であり、kは目標当量比KCMDの算出周期CACALで離散化した離散化時刻である。KCMD=DAF×sinKf1×CACAL×k12ステップS3では、空燃比振動制御フラグFPTが1であるか否かを判別する。空燃比振動制御フラグFPTは、空燃比振動制御の開始時点から所定安定化時間TSTBLが経過すると1に設定される。ステップS3の答が否定NOであるときは直ちに処理を終了する。ステップS3の答が肯定YESであるときは、停止条件フラグFDSTPが1であるか否かを判別するステップS4。停止条件フラグFDSTPは、図3に示す停止条件判定処理において、故障判定処理を停止すべき条件が成立するとき1に設定される。ステップS4の答が否定NOであるときは、図5に示すLAFセンサ故障判定処理を実行するステップS5。ステップS4でFDSTP=1であるときは、実行条件フラグFMCNDを0に戻しステップS6、処理を終了する。図3は、停止条件判定処理のフローチャートである。この処理は、実行条件フラグFMCNDが1であるときに、ECU5のCPUで所定時間例えば100msec毎に実行される。ステップS11では、エンジン回転数NEが所定高回転数NETHHより低いか否かを判別し、その答が肯定YESであるときは、エンジン冷却水温TWが所定低水温TWTHLより高いか否かを判別するステップS12。ステップS11またはS12の答が否定NOであるときは、故障判定を停止すべきであると判定し、停止条件フラグFDSTPを1に設定する。ここで、所定高回転数NETHHは、上記故障判定実行条件の1における上限回転数と同一またはより高い値に設定され、所定低水温TWLは、上記故障判定実行条件の5における所定温度と同一またはより低い温度に設定される。所定高回転数NETHHを実行条件1における上限回転数より高い値に設定し、所定低水温TWLを実行条件5における所定温度より低い温度に設定する場合には、故障判定が開始され易く、停止中断され難くすることができる。ステップS12の答が肯定YESであるときは、ステップS13~S17により、特定運転状態パラメータXOPに基づく停止条件判定を行う。特定運転状態パラメータXOPとしては、吸入空気流量GAIR、気筒吸入空気量GAIRCYL、エンジン回転数NE、または吸気圧PBAの何れか1つが使用される。気筒吸入空気量GAIRCYLは、吸入空気量流量GAIRに基づいて公知の手法例えば特開2011-144683号公報で算出される1TDC期間TDCパルスの発生周期当たりの気筒吸入空気量である。ステップS13では、特定運転状態パラメータXOPが所定下限値XOPLMLより大きいか否かを判別し、その答が肯定YESであるときは、特定運転状態パラメータXOPが所定上限値XOPLMHより小さいか否かを判別するステップS14。ステップS13またはS14の答が否定NOであるときは、前記ステップS19に進む。ステップS14の答が肯定YESであるときは、下記式11により、特定運転状態パラメータXOPの変化量絶対値DXOPAを算出する。式11のXOPZは、特定運転状態パラメータXOPの前回値である。DXOPA=|XOP-XOPZ|11ステップS16では、変化量絶対値DXOPAを下記式12に適用し、変化量積算値IDXOPを算出する。式12のIDXOPZは、変化量積算値IDXOPの前回値である。IDXOP=IDXOPZDXOPA12ステップS17では、変化量積算値IDXOPが所定閾値IDXOPTHより小さいか否かを判別し、その答が肯定YESであるときは、停止条件フラグFDSTPを0に設定する。したがって、LAFセンサ故障判定処理は継続される。ステップS17の答が否定NOであって変化量積算値IDXOPが所定閾値IDXOPTH以上であるときは、停止条件成立と判定し、前記ステップS19に進む。変化量積算値IDXOPは、例えば図4に示す時刻t1においては、変化量絶対値DXOPA1~DXOPA6の和に相当する値となり、特定運転状態パラメータXOPの増加方向及び減少方向の変動履歴が反映される。図5は、図2のステップS5で実行されるLAFセンサ故障判定処理のフローチャートである。ステップS101では、LAFセンサ出力から算出される検出当量比KACTについて、周波数f1成分を抽出するバンドパスフィルタ処理を行い、バンドパスフィルタ処理出力の絶対値振幅を積算することにより、周波数f1成分強度MPTf1を算出する。ステップS102では、周波数f2成分を抽出するバンドパスフィルタ処理を行い、バンドパスフィルタ処理出力の絶対値振幅を積算することにより、周波数f2成分強度MPTf2を算出する。ステップS103では、周波数成分強度の算出開始時点から所定積算時間TINTが経過したか否かを判別し、その答が否定NOである間は直ちに処理を終了する。ステップS103の答が肯定YESとなると、ステップS104に進み、周波数f1成分強度MPTf1が強度判定閾値MPTf1THより小さいか否かを判別する。ステップS104の答が肯定YESであるときは、LAFセンサ出力のリッチ側の応答特性及びリーン側の応答特性がほぼ同様に劣化する第1故障パターンの故障が発生していると判定するステップS105。ステップS104の答が否定NOであるときは、周波数f1成分強度MPTf1及び周波数f2成分強度MPTf2を下記式13に適用し、判定パラメータRTLAFを算出するステップS106。RTLAF=MPTf1/MPTf213ステップS107では、判定パラメータRTLAFが判定閾値RTLAFTHより大きいか否かを判別し、その答が否定NOであるときは、LAFセンサ出力のリッチ側の応答特性及びリーン側の応答特性が非対称に劣化する第2故障パターンの故障が発生していると判定するステップS108。ステップS107の答が肯定YESであるときは、LAFセンサ15は正常である応答特性劣化故障は発生していないと判定するステップS109。以上のように本実施形態では、故障判定期間中にLAFセンサ15の出力信号から算出される検出当量比KACTに含まれる周波数f1成分及び周波数f2成分が抽出され、これらの周波数成分に基づいてLAFセンサ15の応答特性劣化故障の判定が行われる。故障判定開始後において特定運転状態パラメータXOPの変動状態を示し、かつ特定運転状態パラメータXOPの変動履歴が反映される変化量積算値IDXOPが算出され、変化量積算値IDXOPが所定閾値IDXOPTH以上であるときに、故障判定が中断停止される。運転状態パラメータの変動履歴が反映される変化量積算値IDXOPを用いることにより、図14を参照して説明したような従来手法の課題を解決し、特定運転状態パラメータXOPの変動が、周波数f1成分強度MPTf1び周波数f2成分強度MPTf2に影響を及ぼすことに起因する誤判定を防止し、故障判定精度を向上させることができる。また故障判定開始後において特定運転状態パラメータXOPが所定上下限値XOPLMH,XOPLMLの範囲内にあるときに変化量積算値IDXOPの算出が行われるので、特定運転状態パラメータXOPが故障判定に適さない値に変化したときは、変化量積算値IDXOPが停止される。これにより、変化量積算値IDXOPによる不適切な停止条件判定が回避され、故障判定精度の低下を防止できる。また特定運転状態パラメータXOPの変動履歴が反映される変化量積算値IDXOPは、比較的簡便な演算で算出可能であり、しかも特定運転状態パラメータXOPの変動履歴を的確に示すパラメータであるので、ECU5のCPUの演算負荷を増大させることなく故障判定の実行可否をより適切に判定し、判定精度を向上させることができる。本実施形態では、LAFセンサ15が空燃比検出手段に相当し、吸入空気流量センサ7,吸気圧センサ8,及びクランク角度位置センサ11が運転状態パラメータ取得手段に相当し、燃料噴射弁6が燃料量制御手段及び空燃比変動手段の一部を構成し、ECU5が燃料量制御手段及び空燃比制御手段の一部、運転状態パラメータ取得手段の一部、抽出手段、故障判定手段、変動状態パラメータ算出手段、及び判定停止手段を構成する。具体的には、図5のステップS101,S102が抽出手段に相当し、ステップS104~S109が故障判定手段に相当し、図3のステップS13~S16が変動状態パラメータ算出手段に相当し、図3のステップS17,S19及び図2のステップS4が判定停止手段に相当する。[変形例1]図3の処理は、図6に示すように変形してもよい。図6に示す処理は、図3のステップS13及びS14をそれぞれステップS13a及び14aに変更し、かつステップS12aを追加したものである。ステップS12aでは、特定運転状態パラメータXOPの今回値を含む所定数のデータの平均値である移動平均値XOPAVEを算出する。ステップS13aでは、移動平均値XOPAVEが所定下限値XOPLMLより大きいか否かを判別し、その答が肯定YESであるときは、移動平均値XOPAVEが所定上限値XOPLMHより小さいか否かを判別するステップS14a。ステップS13aまたはS14aの答が否定NOであるときは、ステップS19に進み、ステップS14aの答が肯定YESであるときは、ステップS15に進む。この変形例によれば、変化量積算値IDXOPの算出は、移動平均値XOPAVEが所定上下限値の範囲内にあるときに行われるので、特定運転状態パラメータXOPの小さな変動の影響を排除して、停止条件判定を安定化することができる。[変形例2]図3の処理は、図7に示すように変形してもよい。図7に示す処理は、図3のステップS13及びS14をそれぞれステップS13b及び14bに変更し、かつステップS12bを追加したものである。ステップS12bでは、吸入空気流量GAIRの今回値を含む所定数のデータの移動平均値である平均吸入空気流量GAIRAVEを算出する。ステップS13bでは、平均吸入空気流量GAIRAVEが所定下限値GAIRLMLより大きいか否かを判別し、その答が肯定YESであるときは、平均吸入空気流量GAIRAVEが所定上限値GAIRLMHより小さいか否かを判別するステップS14b。ステップS13bまたはS14bの答が否定NOであるときは、ステップS19に進み、ステップS14bの答が肯定YESであるときは、ステップS15に進む。この変形例によれば、変化量積算値IDXOPの算出は、平均吸入空気流量GAIRAVEが所定上下限値の範囲内にあるときに行われるので、吸入空気流量GAIRの小さな変動の影響を排除して、停止条件判定を安定化することができる。[変形例3]上述した実施形態では、特定運転状態パラメータXOPとしては、吸入空気流量GAIR、気筒吸入空気量GAIRCYL、エンジン回転数NE、または吸気圧PBAの何れか1つが使用されるが、これらの運転状態パラメータのうちの2以上のパラメータについて、図3の停止条件判定処理を実行し、何れか1つの運転状態パラメータについて停止条件が成立したとき停止条件フラグFDSTPが1に設定されたときに、故障判定を中断停止するようにしてよい。[第2の実施形態]図8は、本発明の第2の実施形態にかかる停止条件判定処理のフローチャートである。図8の処理は、図3のステップS15~S19を、ステップS21~S29に変えたものである。本実施形態は、以下に説明する点以外は、第1の実施形態と同一である。ステップS21では、下記式21により、特定運転状態パラメータXOPの変化量DXOPを算出する。式21の右辺は、式11の絶対値記号を削除したものに相当する。DXOP=XOP-XOPZ21ステップS22では、変化量DOXPが負の値であるか否かを判別し、その答が肯定YESであるときは、下記式22により、減少量積算値IDXOPNを算出するステップS25。式22のIDXOPNZは、減少量積算値IDXOPNの前回値である。IDXOPN=IDXOPNZ|DXOP|22ステップS22の答が否定NOであるときは、変化量DOXPが正の値であるか否かを判別しステップS23、その答が肯定YESであるときは、下記式23により、増加量積算値IDXOPPを算出する。式23のIDXOPPZは、増加量積算値IDXOPPの前回値である。IDXOPP=IDXOPPZDXOP23ステップS23の答が否定NOであるとき、すなわちDXOP=0であるときは、直ちに処理を終了する。ステップS24またはS25実行後はステップS26に進み、増加量積算値IDXOPPが所定閾値IDXOPTHaより小さいか否かを判別する。この答が肯定YESであるときはさらに減少量積算値IDXOPNが所定閾値IDXOPTHaより小さいか否かを判別するステップS27。ステップS26またはS27の答が否定NOであるときは、故障判定処理を停止すべきであると判定し、停止条件フラグFSTPを1に設定するステップS29。ステップS27の答が肯定YESであるときは、停止条件フラグFSTPを0に設定するステップS28。図8の処理によれば、増加量積算値IDXOPPは、例えば図4に示す時刻t1においては、変化量絶対値DXOPA1,DXOPA3,及びDXOPA56の和に相当する値となり、減少量積算値IDXOPNは、変化量絶対値DXOPA2,DXOPA4,及びDXOPA6の和に相当する値となり、それぞれ特定運転状態パラメータXOPの増加方向及び減少方向の変動履歴が反映される。したがって、比較的簡便な演算によって特定運転状態パラメータXOPの変動履歴を的確に示す増加量積算値IDXOPP及び減少量積算値IDXOPNが得られ、ECU5のCPUの演算負荷を増大させることなく故障判定の実行可否をより適切に判定し、判定精度を向上させることができる。本実施形態では、図8のステップS13,S14,及びS21~24が変動状態パラメータ算出手段に相当し、ステップS26及びS27が判定停止手段に相当する。[変形例]ステップS26及びS27の答がともに否定NOであるとき、停止条件フラグFDSTPを1に設定するようにしてもよい。また、本実施形態も、第1の実施形態における変形例1、2、または3と同様の変形が可能である。[第3の実施形態]図9は、本発明の第3の実施形態にかかる停止条件判定処理のフローチャートである。図9の処理は、図3のステップS15~S19を、ステップS31~S37に変えたものである。本実施形態は、以下に説明する点以外は、第1の実施形態と同一である。ステップS31では、特定運転状態パラメータXOPが最大値XOPMAXより大きいか否かを判別する。最大値XOPMAXは、特定運転状態パラメータXOPが通常とることはない小さな値に初期化されており、最初はステップS31の答は肯定YESとなり、ステップS34で最大値XOPMAXを特定運転状態パラメータXOPの今回値に更新する。ステップS31の答が否定NOであるときは、特定運転状態パラメータXOPが最小値XOPMINより小さいか否かを判別するステップS32。最小値XOPMINは、特定運転状態パラメータXOPが通常とることはない大きな値に初期化されており、最初はステップS32の答は肯定YESとなり、ステップS33で最小値XOPMINを特定運転状態パラメータXOPの今回値に更新する。ステップS32の答が否定NOであるときは、直ちにステップS36に進む。ステップS33またはS34実行後はステップS35に進み、最大値XOPMAXと最小値XOPMINとの差分値が所定閾値DXMMTHより小さいか否かを判別する。この答が否定NOであるときは、故障判定処理を停止すべきであると判定し、停止条件フラグFSTPを1に設定するステップS37。ステップS35の答が肯定YESであるときは、停止条件フラグFSTPを0に設定するステップS36。図9の処理によれば、図10に示すように最大値XOPMAX及び最小値XOPMINが算出され、時刻t2における両者の差分値はDMAXMINで与えられる。したがって、特定運転状態パラメータXOPの増加方向及び減少方向の変動履歴が反映され、比較的簡便な演算によってLAFセンサ15の出力特性が大きく変化する変動履歴を的確に示す変動状態パラメータとしての差分値DMAXMINが得られる。その結果、ECU5のCPUの演算負荷を増大させることなく故障判定の実行可否をより適切に判定し、判定精度を向上させることができる。本実施形態では、図9のステップS13,S14,及びS31~34が変動状態パラメータ算出手段に相当し、ステップS35が判定停止手段に相当する。[変形例]本実施形態も、第1の実施形態における変形例1、2、または3と同様の変形が可能である。[第4の実施形態]図11は、本発明の第4の実施形態にかかる停止条件判定処理のフローチャートである。図11の処理は、図3のステップS15~S19を、ステップS41~S44に変えたものである。なお、図11の処理は所定クランク角度毎に実行される。本実施形態は、以下に説明する点以外は、第1の実施形態と同一である。ステップS41では、特定運転状態パラメータXOPについて周波数f1近傍の周波数成分を抽出するバンドパスフィルタ処理を実行し、フィルタ処理後パラメータXOPBFAを算出する。バンドパスフィルタ処理は、下記式31を用いて行われるが、フィルタ処理後パラメータXOPBFAは、式31を用いて算出されるフィルタ出力XOPBFの絶対値に相当する。式31に適用される特定運転状態パラメータXOPのデータ数N1は、例えば3以上の値に設定される。式31のN及びMは、必要なフィルタ特性に応じた値に設定されるパラメータ整数であり、a及びbは、必要なフィルタ特性に応じた値に設定されるフィルタ係数である。なお、本実施形態では、ステップS41におけるバンドパスフィルタ処理の通過帯域幅は、図5のステップS101において周波数f1成分の抽出に適用されるバンドパスフィルタ処理の通過帯域幅よりも広く設定される。ステップS42では、フィルタ処理後パラメータXOPBFAが所定閾値XOPBFTHより小さいか否かを判別する。この答が否定NOであるときは、故障判定処理を停止すべきであると判定し、停止条件フラグFSTPを1に設定するステップS44。ステップS42の答が肯定YESであるときは、停止条件フラグFSTPを0に設定するステップS43。図11のバンドパスフィルタ処理には、特定運転状態パラメータXOPの今回値とN個の過去値が適用されるので、フィルタ処理後パラメータXOPBFAには、特定運転状態パラメータXOPの変動履歴が反映される。したがって、故障判定に使用される周波数f1成分及び周波数f2成分に大きな影響を与える変動周波数成分が特定運転状態パラメータXOPに含まれている状態を、比較的簡便な演算で確実に判定することができる。その結果、ECU5のCPUの演算負荷を増大させることなく故障判定の実行可否をより適切に判定し、判定精度を向上させることができる。本実施形態では、図11のステップS13,S14,及びS41が変動状態パラメータ算出手段に相当し、ステップS42が判定停止手段に相当する。[変形例]本実施形態も、第1の実施形態における変形例1、2、または3と同様の変形が可能である。なお本発明は上述した実施形態に限るものではなく、種々の変形が可能である。例えば、上述した実施形態では、空燃比制御系の故障判定処理として、LAFセンサ15の応答特性劣化故障判定処理を示したが、本発明は、気筒毎の空燃比が許容限度を超えてばらついているインバランス故障を判定する処理の停止条件判定にも適用可能である。図12は、図2のフローチャートをインバランス故障判定処理用に変形したものであり、図2のステップS2及びS3を削除し、ステップS5をステップS5aに変えたものである。図12のステップS5aでは、LAFセンサ故障判定処理に代えて図13に示すインバランス故障判定処理が実行される。すなわち、インバランス故障判定では、空燃比振動制御は行われず、検出当量比KACTに応じた空燃比フィードバック制御実行中において、検出当量比KACTに含まれる、エンジン回転数に対応する周波数の1/2の周波数成分0.5次周波数成分の強度MIMBに基づいて、故障判定が行われる。図13のステップS111では、検出当量比KACTについて0.5次周波数成分を抽出するバンドパスフィルタ処理を行い、バンドパスフィルタ処理出力の絶対値振幅を積算することにより、0.5次周波数成分強度MIMBを算出する。ステップS112では、0.5次周波数成分強度MIMBが判定閾値MINBTHより大きいか否かを判別する。ステップS112の答が肯定YESであるときは、インバランス故障が発生していると判定するステップS113。ステップS112の答が否定NOであるときは、気筒毎の空燃比ばらつきは許容限度内にあるインバランス故障は発生していないと判定するステップS114。図13に示すインバランス故障判定手法は、基本的には特開2009-270543号公報に示される手法と同一である。本変形例によれば、特定運転状態パラメータXOPの変動によってインバランス故障の判定精度が低下することを防止できる。本変形例では、図13の処理が故障判定手段に相当する。なお、上記第4の実施形態においてこの変形例を適用する場合には、図11のステップS41では、0.5次周波数の近傍の帯域を通過帯域とするバンドパスフィルタ処理を実行する。そのバンドパスフィルタ処理の通過帯域幅は、図13のステップS111において0.5次周波数成分の抽出に適用されるバンドパスフィルタ処理の通過帯域幅よりも広く設定される。またインバランス故障の判定手法は、上述したものに限らず例えば特開2011-144754号公報に示された手法を用いてもよく、またLAFセンサ応答劣化故障の判定手法は、上述したものに限らず例えば特開2010-133418号公報に示された手法を用いてもよい。また本発明は、クランク軸を鉛直方向とした船外機などのような船舶推進機用エンジンなどの空燃比制御装置にも適用が可能である。"