(57)【要約】【課題】この発明は、スロットルバルブ異常診断装置に関し、スロットルバルブが開閉過渡状態であってもスロットルバルブ開度を高精度に推定可能なスロットルバルブ異常診断装置を提供することを目的とする。【解決手段】データクラスタリング部１２では、スロットルバルブへの指示電圧及び開度センサからのセンサ電圧からなる過去のデータを、特定した２つの確率分布に基づき２つのクラスタに分類する。分離超平面推定部１４では、スロットルバルブの開状態と閉状態との境界条件を該過去のデータで特定する。システムパラメータ推定部１６では、所定時刻における該過去のデータと、該時刻よりも１ステップ先のセンサ電圧との間で成立する関係を特定する。特定した境界条件及び該関係に現在のデータを適用して将来のセンサ電圧を推定する。【選択図】図１
2010133384.tif 000002
【特許請求の範囲】
【請求項１】
開閉制御時にスロットルバルブに入力される指示電圧及びスロットル開度センサから出力されるセンサ電圧の過去のデータが従う基準開度若しくはそれよりも開いた開状態に対応した第１確率分布の未知パラメータ及び前記基準開度もしくはそれよりも閉じた閉状態に対応した第２確率分布の未知パラメータを、これらのデータに基づいて最尤推定法により推定し、推定した未知パラメータにより前記第１確率分布及び前記第２確率分布を特定する確率分布特定手段と、
特定後の前記第１確率分布及び前記第２確率分布に基づいて、前記過去のデータの各々が前記開状態のデータであるか、或いは前記閉状態のデータであるかの判定を行い、夫々開状態クラスタ及び閉状態クラスタに分類するデータ分類手段と、
前記開状態と前記閉状態との境界で成立する境界条件を、前記過去のデータに基づいて推定し、推定した未知パラメータにより前記境界条件を特定する境界条件特定手段と、
所定の時刻においてスロットルバルブに入力される指示電圧及びスロットル開度センサから出力されるセンサ電圧のデータと、前記所定の時刻よりも１ステップ先の時刻においてスロットル開度センサから出力されるセンサ電圧のデータとの間に成立する関係の未知パラメータを、前記開状態クラスタのデータに基づいて推定し、推定した未知パラメータにより前記関係をスロットルバルブの開状態に対応した第１規則として特定する第１規則特定手段と、
所定の時刻においてスロットルバルブに入力される指示電圧及びスロットル開度センサから出力されるセンサ電圧のデータと、前記所定の時刻よりも１ステップ先の時刻においてスロットル開度センサから出力されるセンサ電圧のデータとの間に成立する関係の未知パラメータを、前記閉状態クラスタのデータに基づいて推定し、推定した未知パラメータにより前記関係をスロットルバルブの閉状態に対応した第２規則として特定する第２規則特定手段と、
特定した前記境界条件にスロットルバルブに入力された開度指示電圧及びスロットル開度センサから出力された開度センサ電圧のデータの現在値を含むデータを適用して、前記現在値を含むデータが前記開状態に区分されるか、或いは前記閉状態に区分されるかを推定する区分推定手段と、
前記第１規則又は前記第２規則に前記現在値を含むデータを適用してスロットル開度センサから出力される将来のセンサ電圧を推定するセンサ電圧推定手段と、
を備えることを特徴とするスロットルバルブ異常診断装置。
【請求項２】
前記区分推定手段による推定結果及び／又は前記将来のセンサ電圧に基づいてスロットルバルブの異常を判定する異常判定手段を備えることを特徴とする請求項１に記載のスロットルバルブ異常診断装置。
【請求項３】
特定後の前記第１確率分布及び前記第２確率分布に基づいて、スロットルバルブに入力された指示電圧及びスロットル開度センサから出力されたセンサ電圧のデータの現在までの実測値を含むデータが、前記開状態のデータである確率及び前記閉状態のデータである確率を示す帰属確率を演算する帰属確率演算手段と、
前記帰属確率と予め定めた判定確率とを比較することにより前記異常判定手段による判定結果の妥当性を評価する評価手段と、
を備えることを特徴とする請求項２に記載のスロットルバルブ異常診断装置。
【請求項４】
前記開状態と前記閉状態との境界から前記現在値を含むデータがどの程度離れているかを示す距離を演算する距離演算手段と、
前記距離と予め定めた判定値とを比較することにより前記異常判定手段による判定結果の妥当性を評価する評価第２手段と、
を備えることを特徴とする請求項２又は３に記載のスロットルバルブ異常診断装置。
【請求項５】
特定後の前記第１確率分布及び前記第２確率分布に基づいて、前記現在値を含むデータが前記開状態のデータであるか、或いは前記閉状態のデータであるかの判定を行い、前記異常判定手段による判定結果の妥当性を評価する評価第３手段を備えることを特徴とする請求項２乃至４いずれか１項に記載のスロットルバルブ異常診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、スロットルバルブ異常診断装置に関し、特に、スロットルバルブの開閉状態に対応した規則等を利用して、スロットルバルブの異常判定を行うスロットルバルブ異常診断装置に関する。
【背景技術】
【０００２】
従来、例えば特許文献１には、算出したスロットルバルブの推定開度と、検出したスロットルバルブの実開度とを比較してスロットルバルブの異常を判断する異常検出装置が開示されている。この装置では、スロットルバルブの推定開度を算出するために数学モデルが用いられている。この数学モデルは、スロットルバルブを駆動するモータの温度や、このモータに電力を供給するバッテリの電圧が内燃機関の運転状態に応じて変化をするため、これらの変化を考慮して構築されている。具体的には、モータ温度が上昇すると、スロットルバルブ開閉制御の際の目標開度に対する実開度の応答性が徐々に低下する。また、バッテリ電圧が低下すると、同様に実開度の応答性が徐々に低下する。このため、目標開度に対して応答性の低下分の補正を行い、推定開度を高精度に算出している。こうすることで、運転状態を考慮してスロットルバルブの推定開度が推定でき、運転状態を考慮しない場合に比べてスロットルバルブ異常検出を高精度に行うことができる。
【０００３】
【特許文献１】特開平０６−０８１７０２号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
このように、スロットルバルブの開度が高精度に推定できれば、スロットルバルブの異常検出も高精度に行うことができる。ところで、スロットルバルブの開度推定は、スロットルバルブの運転状態をも考慮して行うことができることが好ましい。特に、ある開度を基準にスロットルバルブが短時間で開閉制御される開閉過渡の場合をも考慮して行うことができることが好ましい。しかしながら、特許文献１の装置は開閉過渡の場合を考慮していない。このため、上述した応答性の低下分を考慮したとしても、短時間でスロットルバルブの開度が変化する場合には、高精度に開度推定できるとは限らなかった。したがって、スロットルバルブの異常検出を高精度に行うことができるとは限らなかった。
【０００５】
この発明は、上述のような課題を解決するためになされたもので、スロットルバルブが開閉過渡状態であってもスロットルバルブ開度を高精度に推定可能なスロットルバルブ異常診断装置を提供することを目的とする。
【課題を解決するための手段】
【０００６】
第１の発明は、上記の目的を達成するため、スロットルバルブ異常診断装置であって、
開閉制御時にスロットルバルブに入力される指示電圧及びスロットル開度センサから出力されるセンサ電圧の過去のデータが従う基準開度若しくはそれよりも開いた開状態に対応した第１確率分布の未知パラメータ及び前記基準開度もしくはそれよりも閉じた閉状態に対応した第２確率分布の未知パラメータを、これらのデータに基づいて最尤推定法により推定し、推定した未知パラメータにより前記第１確率分布及び前記第２確率分布を特定する確率分布特定手段と、
特定後の前記第１確率分布及び前記第２確率分布に基づいて、前記過去のデータの各々が前記開状態のデータであるか、或いは前記閉状態のデータであるかの判定を行い、夫々開状態クラスタ及び閉状態クラスタに分類するデータ分類手段と、
前記開状態と前記閉状態との境界で成立する境界条件を、前記過去のデータに基づいて推定し、推定した未知パラメータにより前記境界条件を特定する境界条件特定手段と、
所定の時刻においてスロットルバルブに入力される指示電圧及びスロットル開度センサから出力されるセンサ電圧のデータと、前記所定の時刻よりも１ステップ先の時刻においてスロットル開度センサから出力されるセンサ電圧のデータとの間に成立する関係の未知パラメータを、前記開状態クラスタのデータに基づいて推定し、推定した未知パラメータにより前記関係をスロットルバルブの開状態に対応した第１規則として特定する第１規則特定手段と、
所定の時刻においてスロットルバルブに入力される指示電圧及びスロットル開度センサから出力されるセンサ電圧のデータと、前記所定の時刻よりも１ステップ先の時刻においてスロットル開度センサから出力されるセンサ電圧のデータとの間に成立する関係の未知パラメータを、前記閉状態クラスタのデータに基づいて推定し、推定した未知パラメータにより前記関係をスロットルバルブの閉状態に対応した第２規則として特定する第２規則特定手段と、
特定した前記境界条件にスロットルバルブに入力された開度指示電圧及びスロットル開度センサから出力された開度センサ電圧のデータの現在値を含むデータを適用して、前記現在値を含むデータが前記開状態に区分されるか、或いは前記閉状態に区分されるかを推定する区分推定手段と、
前記第１規則又は前記第２規則に前記現在値を含むデータを適用してスロットル開度センサから出力される将来のセンサ電圧を推定するセンサ電圧推定手段と、を備えることを特徴とする。
【０００７】
また、第２の発明は、第１の発明において、
前記区分推定手段による推定結果及び／又は前記将来のセンサ電圧に基づいてスロットルバルブの異常を判定する異常判定手段を備えることを特徴とする。
【０００８】
また、第３の発明は、第２の発明において、
特定後の前記第１確率分布及び前記第２確率分布に基づいて、スロットルバルブに入力された指示電圧及びスロットル開度センサから出力されたセンサ電圧のデータの現在までの実測値を含むデータが、前記開状態のデータである確率及び前記閉状態のデータである確率を示す帰属確率を演算する帰属確率演算手段と、
前記帰属確率と予め定めた判定確率とを比較することにより前記異常判定手段による判定結果の妥当性を評価する評価手段と、を備えることを特徴とする。
【０００９】
また、第４の発明は、第２又は第３の発明において、
前記開状態と前記閉状態との境界から前記現在値を含むデータがどの程度離れているかを示す距離を演算する距離演算手段と、
前記距離と予め定めた判定値とを比較することにより前記異常判定手段による判定結果の妥当性を評価する評価第２手段と、を備えることを特徴とする。
【００１０】
また、第５の発明は、第２乃至第４いずれか１つの発明において、
特定後の前記第１確率分布及び前記第２確率分布に基づいて、前記現在値を含むデータが前記開状態のデータであるか、或いは前記閉状態のデータであるかの判定を行い、前記異常判定手段による判定結果の妥当性を評価する評価第３手段を備えることを特徴とする。
【発明の効果】
【００１１】
第１の発明によれば、スロットルバルブの開状態に対応した第１規則、閉状態に対応した第２規則という２つの規則から、将来のセンサ電圧を推定できる。これら２つの規則は、指示電圧及びセンサ電圧の過去のデータが従う２つの確率分布、すなわち、開状態に対応した第１確率分布及び閉状態に対応した第２確率分布にこれらの過去のデータのそれぞれを分類し、分類した後のデータ群である開状態クラスタ及び閉状態クラスタに基づいて特定されたものである。そして、この過去のデータは、スロットルバルブの運転状態に関らず蓄積されたものである。したがって、第１の発明によれば、スロットルバルブが開閉過渡状態であっても、このように特定した２つの規則で将来のセンサ電圧を推定できる。
【００１２】
また、第１の発明によれば、開状態と閉状態との境界で成立する境界条件により、指示電圧及びセンサ電圧のデータの現在値を含むデータが開状態、閉状態のどちらに区分されるかを推定することができる。そして、この境界条件による区分と、上述した２つの規則とを組み合わせることで、上記の現在値を含むデータがどちらの状態にあるかで将来のセンサ電圧を推定する規則を切り替えることが可能となる。したがって、第１の発明によれば、上記の現在値を含むデータがどちらの状態であるかを推定し、そのときの将来のセンサ電圧を高精度に推定することができる。
【００１３】
第２の発明によれば、区分推定手段による推定結果及び／又は推定した将来のセンサ電圧に基づいて、スロットルバルブの異常を判定することができる。上述したように、第１の発明によって、現在値を含むデータがどちらの状態であるかが推定され、そのときのセンサ電圧が推定されている。このため、これら２つの推定結果の一方又は両方を判断指標としてスロットルバルブの異常判定をすることができる。
【００１４】
第３の発明によれば、指示電圧及びセンサ電圧のデータの現在までの実測値を含むデータが、上記開状態のデータである帰属確率及び上記閉状態のデータである帰属確率を演算することができる。そして、この実測値を含むデータの信憑性を、演算した帰属確率を判定指標として用いることにより、評価することができる。こうすることで、スロットルバルブ異常判定の結果の妥当性を評価することができる。
【００１５】
第４の発明によれば、開状態と閉状態との境界から上記の現在値を含むデータがどの程度離れているかを示す距離を演算することができる。そして、上記スロットルバルブの異常判定を行う際に用いた上記の現在値のデータの信憑性を、この距離を判定指標として用いることにより、評価することができる。こうすることで、スロットルバルブ異常判定の結果の妥当性を評価することができる。
【００１６】
第５の発明によれば、特定後の前記第１確率分布及び前記第２確率分布に基づいて、上記の現在値を含むデータが、上記開状態のデータであるか、或いは上記閉状態のデータであるかの判定を行うことができる。そして、上記スロットルバルブの異常判定を行う際に用いた上記の現在値のデータの信憑性を、上記の現在値を含むデータがいずれの状態にあるかを判定指標として用いることにより、評価することができる。したがって、スロットルバルブ異常判定の結果の妥当性を評価することができる。
【発明を実施するための最良の形態】
【００１７】
実施形態１．
実施形態１のスロットルバルブ異常診断装置は、車両等に搭載されるＥＣＵ（Electronic Control Unit）内部に構成されている。このため、先ず、図１を用いて、本実施形態に係る装置の説明をする。図１は、ＥＣＵがスロットルバルブの異常判定を行う診断装置として機能する場合の機能ブロック図である。
【００１８】
図１に示すように、ＥＣＵは、システム同定部１０、センサ電圧推定部１８、開閉状態推定部２０及び異常判定部２２から構成されている。また、システム同定部１０は、データクラスタリング部１２、分離超平面推定部１４及びシステムパラメータ推定部１６から構成されている。
【００１９】
（システム同定部１０）
システム同定部１０は、電子制御スロットルシステムにおいて、スロットルバルブの要求開度に応じてスロットルモータを駆動するための電圧（以下、「指示電圧」という。）から、スロットルバルブの開度に応じてスロットル開度センサから検出される電圧（以下、「センサ電圧」という。）までの特性を記述する同定式を、スロットルバルブの開度により規定される２つの状態で同定式が切り替わる区分的アフィン自己回帰（Piece-Wise affine Auto-Regressive eXogeneous；ＰＷＡＲＸ）モデルとして同定を行うように構成されている。ここで、２つの状態とは、所定の基準開度を設定したときに、この基準開度よりも開いた状態（開状態）及び閉じた状態（閉状態）のそれぞれを意味するものとする。所定の基準開度は、スロットルバルブの運転状態に応じて自由に設定することができる。なお、基準開度を全開に設定した場合には、ここで言う開状態は、全開に維持されている状態に対応するものとする。また、全閉に設定した場合には、ここでいう閉状態は、全閉に維持されている状態に対応するものとする。
【００２０】
システム同定部１０は、過去の入出力に関する時系列データに基づいて、後述する式（１）及び（２）のパラメータθ1、θ2、a、bを決定し、ＰＷＡＲＸモデルの同定を行うことを特徴とする。同定は、過去の入力データである指示電圧及び出力データであるセンサ電圧を、過去の時系列データとして用いることにより行う。時刻kにおける指示電圧をu(k)、センサ電圧をy(k)とする。時刻kを現在時刻とすると、これらの過去の時系列データは、u(k-1)，u(k-2)，・・・、y(k-1),y(k-2)，・・・で表される。これら過去の時系列データは、ＥＣＵにバッファリングされているものとする。
【００２１】
過去の時系列データに基づくエンジン特性がＰＷＡＲＸモデルに従うと仮定すると、指示電圧とセンサ電圧との関係は次式（１）及び（２）で表される。
【数１】
2010133384.tif 000003
上式（１）が開状態に対応し、上式（２）が閉状態に対応するものとする。上式（１）及び（２）において、e(k)は式誤差である。また、右肩の添え字「T」は転置行列を意味する。また、x(k)は回帰ベクトルであり、過去の指示電圧及びセンサ電圧を用いて、
x(k)=[y(k-1)  ・・・  y(k-ny)   u(k-1)  ・・・  u(k-nu)]T      ・・・（３）
により与えられる。上式（３）において、kは、k=1,2,・・・,Nの離散時刻の値であり、Nはデータ数である。また、自然数ny,nuは、それぞれモデルを記述するために必要な過去の時系列データの個数を表す。また、θ1及びθ2は、これらのモデルを規定する未知のパラメータである。
【００２２】
また、上式（１）及び（２）に示すように、開状態と閉状態とは、回帰ベクトル空間上における分離超平面で分割されるものとする。分離超平面は、これら２つの状態の境界であり、次式（４）で表される。
aTx(k)＋b=0      ・・・（４）
上式（４）において、a及びbは、分離超平面を規定する未知のパラメータである。
【００２３】
（データクラスタリング部１２）
システム同定部１０内のデータクラスタリング部１２は、過去の時系列データを、開状態に対応する開状態クラスタC1又は閉状態に対応する閉状態クラスタC2に分類するように構成されている。
【００２４】
分類は、まず、観測データベクトルzk=(xT(k)  yT(k)）Tを定義する。そして、この観測データベクトルzkのN個の観測データz1，・・・，zNが、開状態及び閉状態に対応する正規分布からなる混合正規分布に従うと仮定する。
【００２５】
混合正規分布は、パラメータΦ=(α1，α2，μ1，μ2，Σ1，Σ2)を用いて次式（５）のように表される。
【数２】
2010133384.tif 000004
上式（５）において、α1、α2は、混合正規分布の混合割合であり、α1，α2≧0，α1+α2=1の関係を満たす。また上式（５）〜（７）において、μ1，μ2は、開状態、閉状態のそれぞれに対応する正規分布の平均ベクトルを示す。同様に、Σ1、Σ2は、それぞれに対応する正規分布の共分散行列を示す。添え字「nz」は、観測データベクトルzkの次元数を示す。さらにdetΣは、Σの行列式を示し、行列の右肩の添え字「-1」は逆行列を意味する。
【００２６】
続いて、混合正規分布のパラメータΦを最尤推定法により求める。尤度関数
【数３】
2010133384.tif 000005
を最大化するパラメータΦは、ＥＭアルゴリズムを用いることにより推定することができる。ＥＭアルゴリズムでは、パラメータΦを反復的に計算することができる。
【００２７】
続いて、求めたパラメータΦを用いて、観測データの分類を行う。ここで、観測データが開状態クラスタC1に帰属する帰属確率をα1p1(zk;Φ)、閉状態クラスタC2に帰属する帰属確率をα2p2(zk;Φ)とする。分類は、得られた帰属確率がより高い部分分布に対応するクラスタに分類する。すなわち、α1p1(zk;Φ)≧α2p2(zk;Φ)のデータを開状態クラスタC1に、α1p1(zk;Φ)＜α2p2(zk;Φ)のデータを閉状態クラスタC2に分類する。これにより、観測データに対応する時刻k=1,2,・・・,Nをそれぞれのクラスタに分類することができる。
【００２８】
（分離超平面推定部１４）
システム同定部１０内の分離超平面推定部１４は、開状態と閉状態とを分割する回帰空間上の分離超平面のパラメータを推定するように構成されている。上式（４）により与えられる分離超平面のパラメータa,bは、ソフトマージンサポートベクターマシンを利用することにより推定できる。
【００２９】
サポートベクターマシンは、与えられたデータが２つのクラスのどちらに属するかを判断することのできるアルゴリズムであり、全てのサンプルデータに対して正しいクラスが存在する（線形分離可能）ときに成立する超平面のパラメータa,bを調整する方法である。しかし、このような線形分離を実現する超平面は、１つに決まらない。そこで、超平面から各クラスまでの距離をマージン（1/||a||）とし、このマージンを最大にするように超平面のパラメータa,bを決定する。
【００３０】
ここで、全てのサンプルデータに対して常に正しいクラスが存在するとは限らない（線形分離可能でない）。このような場合には、サポートベクターマシンに分類誤差項を導入したソフトマージンサポートベクターマシンを利用する。すなわち、データが超平面を超えて反対側のクラスに入ってしまった場合に、その超えた距離νの総和を最小限にするように超平面のパラメータa,bを調整する。まとめると、ソフトマージンサポートベクターマシンを利用すれば、パラメータa,bは、次式（８）に示す２次最適化問題を解くことにより推定できる。
【数４】
2010133384.tif 000006
このように推定されたパラメータa,bは、ＥＣＵに記憶される。
【００３１】
（システムパラメータ推定部１６）
システム同定部１０内のシステムパラメータ推定部１６は、データクラスタリング部１２で分類した離散時刻に対応するデータを用いて、パラメータθ1及びθ2を推定するように構成されている。パラメータθ1,θ2は、次式（９）及び（１０）により推定できる。
【数５】
2010133384.tif 000007
上式（９）及び（１０）において、N1,N2は、それぞれのクラスタC1,C2に含まれるデータ数を示す。このように推定されたパラメータθ1^,θ2^は、ＥＣＵに記憶される。
【００３２】
以上により、システム同定部１０では、開状態と閉状態とで同定式が切り替わるＰＷＡＲＸモデルの同定が行われる。まとめると、現在の時刻を時刻kとすると、時刻k+1におけるセンサ電圧の推定値y^(k+1)は、時刻k-1までの指示電圧及びセンサ電圧のデータ履歴に基づいて推定したパラメータθ1^,θ2^、を用い、次式（１１）及び（１２）のように記述される。尚、ここで使用した添え字「^」は、推定値を表すものとする。
【数６】
2010133384.tif 000008
尚、上式（１１）が本発明の第１規則に、上式（１２）が本発明の第２規則にそれぞれ該当する。また、上式（１１）及び（１２）において、x(k+1)は回帰ベクトルであり、時刻kまでの指示電圧及びセンサ電圧を用いて、
x(k+1)=[y(k)  ・・・  y(k-ny+1)  u(k)  ・・・  u(k-nu+1)]T      ・・・（１３）
により与えられることになる。
【００３３】
このように、時刻k+1におけるセンサ電圧の推定値y^(k+1)は、時刻k-1までの指示電圧及びセンサ電圧のデータ履歴、センサ電圧の現在のデータを逐次帰属させることでその都度行うことができる。したがって、モデル同定の自由度が高いという効果を有する。
【００３４】
（センサ電圧推定部１８）
センサ電圧推定部１８は、パラメータa,b及び上式（１３）を用い、現在のスロットルバルブが開状態であるか、或いは閉状態であるかを推定し、その推定結果と、上式（１１）又は（１２）を用いてセンサ電圧の推定値y^(k+1)を推定するように構成されている。
【００３５】
具体的には、先ず、時刻k-1までのデータ履歴により推定したパラメータa,bと、このデータ履歴に時刻kのセンサ電圧の出力データを加えて生成した回帰ベクトルx(k+1)とから、aTx(k+1)＋bを演算する。そして、aTx(k+1)＋b≧0であれば、センサ電圧推定部１８は、現在のスロットルバルブ開度が開状態であると推定し、開閉状態予測フラグi^(k+1)=1がＥＣＵにバッファリングされる。一方、aTx(k+1)＋b＜0であれば、センサ電圧推定部１８は、現在のスロットルバルブ開度が閉状態であると推定し、開閉状態予測フラグi^(k+1)=0がＥＣＵにバッファリングされる。この開閉状態予測フラグi^(k+1)と、上式（１１）又は（１２）を用いて、時刻k+1におけるセンサ電圧の推定値y^(k+1)を推定する。推定したy^(k+1)は、ＥＣＵにバッファリングされる。
【００３６】
（開閉状態推定部２０）
開閉状態推定部２０は、パラメータa,b及び上式（１３）を用い、現在のスロットルバルブ開度が開状態であるか、或いは閉状態であるかを推定するように構成されている。この推定は、センサ電圧推定部１８で上述したように、aTx(k+1)＋bの演算結果を利用することにより行われる。すなわち、aTx(k+1)＋b≧0であれば、現在のスロットルバルブ開度は開状態であると推定され、開閉状態予測フラグi^(k+1)=1がＥＣＵにバッファリングされる。aTx(k+1)＋b＜0であれば、閉状態であると推定され、開閉状態予測フラグi^(k+1)=0がＥＣＵにバッファリングされる。
【００３７】
（異常判定部２２）
異常判定部２２は、センサ電圧推定部１８で推定したy^(k+1)及び開閉状態推定部２０で推定した開閉状態予測フラグi^(k+1)に基づいてスロットルバルブの異常判定を行うように構成されている。
【００３８】
y^(k+1)に基づく異常判定は、y^(k+1)と、時刻k+1における実際のセンサ電圧y(k+1)との偏差を演算することにより行われる。具体的には、これらの値の差Δy＝|y^(k+1)-y(k+1)|を演算し、予め定めた設定値（閾値）との比較を行う。この閾値は、スロットルバルブの判定精度の要求度に応じて自由に設定することが可能である。そして、Δyが閾値以上であれば、推定したセンサ電圧と、実際のセンサ電圧との間に乖離が見られることになる。このような場合に、異常判定部２２は、異常信号を発する。
【００３９】
i^(k+1)に基づく異常判定は、時刻k+1における実際のスロットルバルブ開度をi(k+1)としたときに、i^(k+1)とi(k+1)との比較により行われる。i(k+1)は、実際のセンサ電圧を利用してもよいが、例えばスロットルバルブの下流に配置した流量センサ等の出力等、ＥＣＵが取得可能なパラメータであって、スロットルバルブの開閉状態の代用値として用いることが可能なものであれば特に限定されない。また、このi(k+1)は、開閉状態予測フラグi^(k+1)と対応させて、開状態であればi(k+1)=1、閉状態であればi(k+1)=0とすることができる。そして、i^(k+1)の値がi(k+1)の値と異なる場合には、推定したスロットルバルブの開閉状態と、実際の状態との間に乖離が見られることになる。このような場合に、異常判定部２２は、異常信号を発する。
【００４０】
そして、異常判定は、上記の異常信号が、予め定めた判定時間内において、設定回数検出されるか否かで判定される。具体的には、Δy又はi^(k+1)に基づく上記の異常信号が検出された場合、異常判定部２２に別途内蔵した異常カウンタのカウンタ値が増加される。そして、予め決めた判定時間内において、このカウンタ値が設定値以上に到達した場合に、スロットルバルブ異常として判定される。本実施形態においては、上記の異常信号が、Δy及びi^(k+1)という２つの診断要素に基づいて検出される可能性がある。したがって、いずれか一方の判定結果が誤っていた場合に、最終判定を誤るという誤判定を防止することができる。
【００４１】
[実施形態１における具体的処理]
次に、図２を参照して本実施形態においてＥＣＵが実行するスロットルバルブ異常判定ルーチンを説明する。尚、図２のルーチンは、所定時間毎に繰り返し実行されるものとする。図２に示すルーチンによれば、まず、モデルのパラメータが読み込まれる（ステップ１００）。上述の通り、システム同定部１０では、時刻k-1まで時系列データに基づいて、ＰＷＡＲＸモデルのパラメータθ1^,θ2^が推定され、これらはＥＣＵに記憶されている。同様に、分離超平面推定部１４では、分離超平面のパラメータa,bが推定され、これらはＥＣＵに記憶されている。このため、ＥＣＵには、これらのパラメータが読み込まれる。ステップ１００に続いて、時刻kまでの指示電圧及びセンサ電圧の時系列データが読み込まれ（ステップ１１０）、回帰ベクトルx(k+1)が生成される。
【００４２】
ステップ１１０に続いて、ステップ１００で読み込んだパラメータa,bと、ステップ１１０で生成された回帰ベクトルx(k+1)とから、aTx(k+1)＋bが演算される（ステップ１２０）。そして、aTx(k+1)＋b≧0であればi^(k+1)=1が、i^(k+1)＜0であればi^(k+1)=0がＥＣＵにバッファリングされる（ステップ１３０）。
【００４３】
ステップ１３０に続いて、スロットルバルブが開状態か否かが判定される（ステップ１４０）。開状態か否かは上記ステップ１３０でバッファリングされたi^(k+1)が用いられる。i^(k+1)=1であれば、ステップ１５０に進み、開状態に対応したセンサ電圧の推定値y^(k+1)が演算される。このy^(k+1)は、ステップ１６０でＥＣＵにバッファリングされる。一方、本ステップにおいて、i^(k+1)=1でない場合には、ステップ１７０に進み、閉状態に対応したセンサ電圧の推定値y^(k+1)が演算され、ステップ１６０同様、ステップ１８０でＥＣＵにバッファリングされる。ステップ１６０又は１８０に続いて、スロットルバルブの異常判定が実行される（ステップ１９０）。
【００４４】
次に、ステップ１９０で実行されるスロットルバルブの異常判定のサブルーチンを図３のフローチャートを参照して説明する。図３のルーチンによれば、まず、異常判定部２２に内蔵した異常カウンタがONか否かが判定される（ステップ２００）。ONでないと判定された場合には、異常カウンタがONにされる（ステップ２１０）。一方、異常カウンタがONと判定された場合には、ステップ２２０に進む。
【００４５】
ステップ２２０では、時刻k+1における実際のスロットルバルブ開度i(k+1)が読み込まれ、同様に、ステップ２３０では、時刻k+1における実際のセンサ電圧y(k+1)が読み込まれる。
【００４６】
ステップ２４０では、開閉状態予測フラグi^(k+1)と、実際のスロットルバルブ開度i(k+1)とが比較される。上述の通り、開閉状態予測フラグi^(k+1)は、ステップ１３０の処理を実行することによりＥＣＵにバッファリングされている。このため、このi^(k+1)とステップ２２０で読み込んだi(k+1)とが比較される。そして、i^(k+1)=i(k+1）の場合には、ステップ２６０に進む。そうでない場合には、ステップ２５０に進み、異常カウンタのカウンタ値が１増加される。
【００４７】
ステップ２６０では、推定したセンサ電圧y^(k+1)と、実際のセンサ電圧y(k+1)とが比較される。上述の通り、センサ電圧y^(k+1)は、ステップ１６０又は１８０の処理を実行することによりＥＣＵにバッファリングされている。このため、このy^(k+1)とステップ２３０で読み込んだy(k+1)とが比較される。そして、これらの間の偏差が予め決められた閾値よりも小さい場合には、ステップ２７０に進み、異常カウンタがリセットされる。そうでない場合には、ステップ２５０に進み、異常カウンタのカウンタ値が１増加される。
【００４８】
続いて、ステップ２８０では、異常カウンタのカウンタ値が予め定めた設定値（閾値）以上か否かが判定される。そして、カウンタ値が閾値以上の場合には、ステップ３００に進み、スロットルバルブは異常と判定される。そうでない場合には、ステップ２９０に進み、スロットルバルブは正常と判定される。
【００４９】
以上説明した通り、図２及び３に示すルーチンによれば、過去の時系列データに基づいて同定したＰＷＡＲＸモデルを用いて、将来のセンサ電圧を逐次推定することができる。また、図２及び３に示すルーチンによれば、ＰＷＡＲＸモデルの同定に際して推定した分離超平面のパラメータを用いて、開閉状態を予測できる。そして、これら２つの診断要素を用いることにより、スロットルバルブの異常を高精度に判定することができる。
【００５０】
尚、本実施形態では、異常判定に際して、センサ電圧推定部１８で推定したセンサ電圧の推定値y^(k+1)及び開閉状態予測フラグi^(k+1)の値を共に用いたが、これらのうちいずれか一方のみを用いて異常判定を行ってもよい。尚、本変形例は後述する実施形態３、４においても同様に適用可能である。
【００５１】
尚、上述した実施形態においては、データクラスタリング部１２において、混合正規分布のパラメータΦを推定し、上式（５）を特定することにより前記第１の発明における「確率分布特定手段」が、データクラスタリング部１２において、観測データの分類を行うことにより前記第１の発明における「データ分類手段」が、分離超平面推定部１４において、分離超平面のパラメータa,bを推定し、上式（４）を特定することにより前記第１の発明における「境界条件特定手段」が、システムパラメータ推定部１６において、パラメータθ1^,θ2^を推定し、上式（１１）及び（１２）を特定することにより前記第１の発明における「第１規則特定手段」及び「第２規則特定手段」が、それぞれ実現されている。また、ＥＣＵが上記ステップ１２０、１３０の処理を実行することにより前記第１の発明における「区分推定手段」が、ＥＣＵが上記ステップ１５０又は１７０の処理を実行することにより前記第１の発明における「センサ電圧推定手段」が、それぞれ実現されている。
【００５２】
また、上述した実施形態においては、ステップ２４０〜２７０の一連の処理を実行することにより前記第２の発明における「異常判定手段」が実現されている。
【００５３】
実施形態２．
実施形態２のスロットルバルブ異常診断装置は、センサ電圧推定部１８で推定したセンサ電圧の推定値y^(k+1)を用いるのに加えて、現在のセンサ電圧の実測値y(k)を含む観測データが上述した開状態クラスタC1又は閉状態クラスタC2に帰属する帰属確率を用いて異常判定を行うことをその特徴とする。センサ電圧の推定値y^(k+1)を用いて異常判定を行うことについては実施形態１に既述したとおりである。このため、以下では実施形態１との相違点を中心に説明し、同様の事項については省略又は簡略する。
【００５４】
図４は、実施形態２のＥＣＵがスロットルバルブの異常判定を行う診断装置として機能する場合の機能ブロック図である。図４に示すように、ＥＣＵは、システム同定部１０、センサ電圧推定部１８、帰属確率演算部２４及び異常判定部２６から構成されている。
【００５５】
（帰属確率演算部２４）
帰属確率演算部２４は、パラメータΦ、及び上式（６）又は（７）を用い、時刻kで実測されたセンサ電圧y(k)を含む観測データが、開状態クラスタC1及び閉状態クラスタC2に帰属する帰属確率を演算するように構成されている。開状態クラスタC1及び閉状態クラスタC2は、データクラスタリング部１２の説明の際に述べたように、過去の時系列データを分類するときに用いた２つのクラスタが用いられる。
【００５６】
観測データベクトルzkは、そのベクトルを構成するy(k)にセンサ電圧の実測値を適用した観測データベクトルである。すなわち、現在時刻kにおいて取得された実際のセンサ電圧y(k)のデータと、過去の指示電圧及びセンサ電圧のデータとから、観測データベクトルzk=(xT(k)  yT(k)）Tを定義する。そして、この観測データベクトルzkと、データクラスタリング部１２で推定されたパラメータΦとを用いて、α1p1(zk;Φ)又はα2p2(zk;Φ)を演算する。α1p1(zk;Φ)又はα2p2(zk;Φ)は、開閉状態予測フラグの値に対応する帰属確率である。すなわち、センサ電圧推定部１８で推定した開閉状態予測フラグi^(k+1)=1であれば、演算される帰属確率はα1p1(zk;Φ)であり、i^(k+1)=0であればα2p2(zk;Φ)である。演算後のα1p1(zk;Φ)又はα2p2(zk;Φ)は、ＥＣＵにバッファリングされる。
【００５７】
（異常判定部２６）
異常判定部２６は、センサ電圧推定部１８で推定したy^(k+1)及び帰属確率演算部２４で演算した帰属確率に基づいてスロットルバルブの異常判定を行うように構成されている。この異常判定のうち、y^(k+1)、すなわちy^(k+1)と実際のセンサ電圧であるy(k+1)の差Δyに基づくものについては、実施形態１で既述した内容と重複することから説明を省略する。
【００５８】
帰属確率に基づく異常判定は、α1p1(zk;Φ)又はα2p2(zk;Φ)と判定のために設定された判定確率との比較により行われる。この判定確率は、時刻kの実測値を含む観測データがクラスタC1又はクラスタC2に帰属する帰属確率の許容される下限値として設定されるものであり、別途ＥＣＵに記憶されている。つまり、帰属確率α1p1(zk;Φ)又はα2p2(zk;Φ)がこの下限値よりも低い場合には、観測データに異常の可能性が認められる。観測データに異常の可能性が認められれば、観測データを構成したデータに異常の可能性が認められることになる。このため、異常判定部２６は、該当する帰属確率とこの下限値とを比較し、帰属確率が下限値よりも低い場合には、異常信号を発するように構成されている。
【００５９】
異常判定は、上記の異常信号が、予め定めた判定時間内において、設定回数検出されるか否かで判定される。具体的には、y^(k+1)及び帰属確率に基づく上記の異常信号が検出された場合、異常判定部２６に別途内蔵した異常カウンタのカウンタ値が増加される。そして、予め決めた判定時間内において、このカウンタ値が設定値以上に到達した場合に、スロットルバルブの異常として判定される。本実施形態においては、上記の異常信号が、本実施形態においては、上記の異常信号がΔyに基づいて検出されると共に、帰属確率に基づいて検出される可能性があるため、データの信憑性を確認することができる。したがって、Δyに基づく判定結果の信頼性を上げることができる。
【００６０】
[実施形態２における具体的処理]
次に、図５を参照して本実施形態においてＥＣＵが実行するスロットルバルブ異常判定ルーチンを説明する。尚、図５中、ステップ１００〜ステップ１８０は、実施形態１と同一の処理が実行されるため説明を省略する。
【００６１】
ステップ１６０に続いて、観測データベクトルzk=(xT(k)  yT(k)）TのクラスタC1に帰属する帰属確率α1p1(zk;Φ)が演算され（ステップ３１０）、ＥＣＵにバッファリングされる（ステップ３２０）。同様に、ステップ１８０に続いて、観測データベクトルzk=(xT(k)  yT(k)）TのクラスタC2に帰属する帰属確率α2p2(zk;Φ)が演算され（ステップ３３０）、ＥＣＵにバッファリングされる（ステップ３４０）。ステップ３２０又は３４０に続いて、スロットルバルブの異常判定が実行される（ステップ３５０）。
【００６２】
次に、ステップ３５０で実行されるスロットルバルブの異常判定のサブルーチンを図６のフローチャートを参照して説明する。尚、図６中、ステップ２００、２１０は、実施形態１と同一の処理が実行されるため説明を省略する。
【００６３】
ステップ３６０では、時刻k+1における実際のセンサ電圧y(k+1)が読み込まれる。そして、推定したセンサ電圧y^(k+1)と、実際のセンサ電圧y(k+1)とが比較される（ステップ３７０）。上述の通り、センサ電圧y^(k+1)は、ステップ１６０又は１８０の処理を実行することによりＥＣＵにバッファリングされている。このため、このy^(k+1)とステップ３６０で読み込んだy(k+1)とが比較される。そして、これらの間の偏差が予め決められた閾値よりも小さい場合には、ステップ３９０に進む。そうでない場合には、ステップ３８０に進み、異常カウンタのカウンタ値が１増加される。
【００６４】
ステップ３９０では、α1p1(zk;Φ)又はα2p2(zk;Φ)と予め決められた閾値との比較がされる。上述の通り、ステップ３２０又は３４０による処理の結果、開閉状態予測フラグの値に対応する帰属確率がＥＣＵにバッファリングされている。また、予め決められた閾値は、帰属確率の下限値として、別途ＥＣＵに記憶されている。このため、ステップ３９０では、この帰属確率と下限値とが比較される。そして、帰属確率が下限値以上の場合には、ステップ４００に進み、異常カウンタがリセットされる。そうでない場合には、ステップ３８０に進み、異常カウンタのカウンタ値が１増加される。
【００６５】
続いて、ステップ４１０では、異常カウンタのカウンタ値が予め定めた設定値（閾値）以上か否かが判定される。そして、カウンタ値が閾値以上の場合には、ステップ４３０に進み、スロットルバルブは異常と判定される。そうでない場合には、ステップ４２０に進み、スロットルバルブは正常と判定される。
【００６６】
以上説明した通り、図５及び６に示すルーチンによれば、過去の時系列データに基づいて同定したＰＷＡＲＸモデルを用いて、将来のセンサ電圧を逐次推定することができる。また、図５及び６に示すルーチンによれば、現在のセンサ電圧の実測値を含む観測データの帰属確率を演算することができる。そして、将来のセンサ電圧に基づく異常判定を行うと共に、帰属確率に基づいて観測データの信憑性を確認することにより、異常判定の精度を上げることができる。
【００６７】
尚、本実施形態では、異常判定に際して、センサ電圧推定部１８で推定したセンサ電圧の推定値y^(k+1)を用いたが、実施形態１同様、開閉状態予測フラグi^(k+1)を診断要素に加えて異常判定を行ってもよい。また、センサ電圧の推定値y^(k+1)の代わりに、実施形態１で説明した開閉状態予測フラグi^(k+1)を診断要素として用いて異常判定を行ってもよい。
【００６８】
また、上述した実施形態においては、ステップ３１０又は３３０の処理を実行することにより、前記第３の発明における「帰属確率演算手段」が、ステップ３９０、３８０、４００の処理を実行することにより前記第３の発明における「評価手段」が、それぞれ実現されている。
【００６９】
実施形態３．
実施形態３のスロットルバルブ異常診断装置は、実施形態１の異常判定に際して、開状態と閉状態との境界である分離超平面からの距離を用いて異常判定の妥当性を評価することをその特徴とする。このため、上述した実施形態１との相違点を中心に説明し、同様の事項については、その説明を省略又は簡略する。
【００７０】
図７は、実施形態３のＥＣＵがスロットルバルブの異常判定を行う診断装置として機能する場合の機能ブロック図である。図７に示すように、ＥＣＵは、システム同定部１０、センサ電圧推定部１８、開閉状態推定部２０、距離演算部２８及び異常判定部３０から構成されている。
【００７１】
（距離演算部２８）
距離演算部２８は、パラメータa,b及び上式（１３）を用い、dk=|aTx(k+1)＋b|を演算するように構成されている。このdkは、上式（４）で示される分離超平面からの距離に相当する。具体的には、時刻k-1までのデータ履歴により推定したパラメータa,bと、このデータ履歴に時刻kのデータを加えて生成した回帰ベクトルx(k+1)とから、距離dk=|aTx(k+1)＋b|を演算する。演算したdkは、ＥＣＵにバッファリングされる。
【００７２】
（異常判定部３０）
異常判定部３０は、センサ電圧推定部１８で推定したy^(k+1)及び開閉状態推定部２０で推定したi^(k+1)に基づいてスロットルバルブの異常判定を行うと共に、距離演算部２８で演算した距離dkに基づいてスロットルバルブの異常判定を行うように構成されている。この異常判定のうち、y^(k+1)及びi^(k+1)に基づくものについては実施形態１で既述した内容と重複することから説明を省略する。
【００７３】
距離dkに基づく異常判定は、この距離dkと、予め決められた設定値（閾値）とを比較することにより行われる。この閾値は、スロットルバルブ異常の判定精度の要求度に応じて自由に設定することが可能である。そして、距離dkが閾値未満であれば、y^(k+1)、i^(k+1)を演算する際に用いたデータの信憑性が低い可能性が認められる。異常判定部３０は、距離dkが閾値未満の場合には、異常信号を発するように構成されている。
【００７４】
異常判定は、上記の異常信号が、予め定めた判定時間内において、設定回数検出されるか否かで判定される。具体的には、Δy、i^(k+1)又は距離dkに基づく上記の異常信号が検出された場合、異常判定部３０に別途内蔵した異常カウンタのカウンタ値が増加される。そして、予め決めた判定時間内において、このカウンタ値が設定値以上に到達した場合に、スロットルバルブの異常として判定される。本実施形態においては、上記の異常信号がΔy及びi^(k+1)に基づいて検出されると共に、距離dkに基づいて検出される可能性があることから、データの信憑性を確認することができる。したがって、Δy及びi^(k+1)に基づく判定結果の信頼性を上げることができる。
【００７５】
[実施形態３における具体的処理]
次に、図８を参照して本実施形態においてＥＣＵが実行するスロットルバルブ異常判定ルーチンを説明する。尚、図８中、ステップ１００〜ステップ１８０は、実施形態１と同一の処理が実行されるため説明を省略する。
【００７６】
ステップ１６０に続いて、分離超平面からの距離dkが演算される（ステップ４４０）。上述の通り、距離演算部２８で、距離dkが演算される。この距離dkは、ステップ４５０でＥＣＵにバッファリングされる。ステップ４５０に続いて、スロットルバルブの異常判定が実行される（ステップ４６０）。
【００７７】
次に、ステップ４６０で実行されるスロットルバルブの異常判定のサブルーチンを図９のフローチャートを参照して説明する。尚、図９中、ステップ２００〜２５０は、実施形態１と同一の処理が実行されるため説明を省略する。
【００７８】
ステップ２６０では、推定したセンサ電圧y^(k+1)と、実際のセンサ電圧y(k+1)とが比較される。上述の通り、センサ電圧y^(k+1)は、ステップ１６０又は１８０の処理を実行することによりＥＣＵにバッファリングされている。このため、このy^(k+1)とステップ２３０で読み込んだy(k+1)とが比較される。そして、これらの間の偏差が予め決められた閾値よりも小さい場合には、ステップ４７０に進む。そうでない場合には、ステップ２５０に進み、異常カウンタのカウンタ値が１増加される。
【００７９】
ステップ４７０では、距離dkと予め決められた閾値との比較がされる。上述の通り、ステップ４５０の処理を実行することにより、距離dkがＥＣＵにバッファリングされている。このため、ステップ４７０では、この距離dkと閾値とが比較される。そして、距離dkが閾値以上の場合には、ステップ４８０へ進み、異常カウンタがリセットされる。そうでない場合には、ステップ２５０へ進み、異常カウンタのカウンタ値が１増加される。
【００８０】
続いて、ステップ４９０では、異常カウンタのカウンタ値が予め定めた設定値（閾値）以上か否かが判定される。そして、カウンタ値が閾値以上の場合には、ステップ５１０に進み、スロットルバルブは異常と判定される。そうでない場合には、ステップ５００に進み、スロットルバルブは正常と判定される。
【００８１】
以上説明した通り、図８及び９に示すルーチンによれば、過去の時系列データに基づいて同定したＰＷＡＲＸモデルを用いて、将来のセンサ電圧を逐次推定することができる。また、図８及び９に示すルーチンによれば、ＰＷＡＲＸモデルの同定に際して推定した分離超平面のパラメータを用いて、開閉状態を予測できる。さらに、図８及び９に示すルーチンによれば、この分離超平面のパラメータを用いて、分離超平面からの距離を演算することができる。そして、将来のセンサ電圧及び開閉状態に基づいてスロットルバルブの異常判定を行うと共に、距離dkに基づいてデータの信憑性を確認することにより、異常判定の精度を上げることができる。
【００８２】
尚、本実施形態は、実施形態２との組み合わせが可能である。すなわち、本実施形態に、実施形態２で説明した現在のセンサ電圧の実測値y(k)を含む観測データの帰属確率を加えて異常判定を行ってもよい。こうすることで、異常判定の精度を更に上げることができる。尚、本変形例は、後述する実施形態４においても同様に適用可能である。
【００８３】
尚、上述した実施形態においては、ステップ４４０の処理を実行することにより、前記第４の発明における「距離演算手段」が、ステップ４７０、２５０、４８０の処理を実行することにより前記第４の発明における「評価第２手段」が、それぞれ実現されている。
【００８４】
実施形態４．
実施形態４のスロットルバルブ異常診断装置は、実施形態１の異常判定に際して、センサ電圧の推定値y^(k+1)を含む観測データが、上述した開状態クラスタC1及び閉状態クラスタC2に帰属する帰属確率を用いて異常判定の妥当性を評価することをその特徴とする。このため、上述した実施形態１との相違点を中心に説明し、同様の事項については、その説明を省略又は簡略する。
【００８５】
図１０は、実施形態４のＥＣＵがスロットルバルブの異常判定を行う診断装置として機能する場合の機能ブロック図である。図１０に示すように、ＥＣＵは、システム同定部１０、センサ電圧推定部１８、開閉状態推定部２０、帰属確率演算部３２及び異常判定部３４から構成されている。
【００８６】
（帰属確率演算部３２）
帰属確率演算部３２は、パラメータΦ、及び上式（６）又は（７）を用い、時刻k+1におけるセンサ電圧の推定値y^(k+1)を含む観測データベクトルが、開状態クラスタC1及び閉状態クラスタC2に帰属する帰属確率を演算するように構成されている。開状態クラスタC1及び閉状態クラスタC2は、データクラスタリング部１２の説明の際に述べたように、過去の時系列データを分類するときに用いた２つのクラスタが用いられる。
【００８７】
帰属確率は、先ず、データクラスタリング部１２で定義した観測データベクトルzk同様、観測データベクトルzk+1=(xT(k+1)  yT(k+1)）Tを定義する。すなわち、本実施形態では、センサ電圧推定部１８で推定したy^(k+1)と、時刻kまでの指示電圧及びセンサ電圧のデータとから、観測データベクトルzk+1を定義する。そして、この観測データベクトルzk+1と、データクラスタリング部１２で推定されたパラメータΦとを用いて、α1p1(zk+1;Φ)及びα2p2(zk+1;Φ)を演算する。演算後のα1p1(zk+1;Φ)、α2p2(zk+1;Φ)は、ＥＣＵにバッファリングされる。
【００８８】
（異常判定部３４）
異常判定部３４は、センサ電圧推定部１８で推定したy^(k+1)及び開閉状態推定部２０で推定したi^(k+1)に基づいてスロットルバルブの異常判定を行うと共に、帰属確率演算部３２で演算した帰属確率α1p1(zk+1;Φ)、α2p2(zk+1;Φ) に基づいてスロットルバルブの異常判定を行うように構成されている。この異常判定のうち、y^(k+1)及びi^(k+1)に基づくものについては実施形態１で既述した内容と重複することから説明を省略する。
【００８９】
帰属確率に基づく異常判定は、開状態クラスタC1への帰属確率α1p1(zk+1;Φ)と閉状態クラスタC2への帰属確率α2p2(zk+1;Φ)とを比較することにより行われる。この比較は、α1p1(zk+1;Φ)≧α2p2(zk+1;Φ)の場合には、開状態である確率が高いが、α1p1(zk+1;Φ)＜α2p2(zk+1;Φ)の場合には、閉状態である確率が高いことを示す。つまり、開閉状態予測フラグがi^(k+1)=1と判定されているにもかかわらず、閉状態である確率が高いとされたときには、y^(k+1)を演算する際に用いたデータの信憑性が低い可能性が認められる。異常判定部３４は、開閉状態フラグの結果が帰属確率の比較結果と異なるような場合には、異常信号を発するように構成されている。
【００９０】
異常判定は、上記の異常信号が、予め定めた判定時間内において、設定回数検出されるか否かで判定される。具体的には、Δy、i^(k+1)又は帰属確率に基づく上記の異常信号が検出された場合、異常判定部３４に別途内蔵した異常カウンタのカウンタ値が増加される。そして、予め決めた判定時間内において、このカウンタ値が設定値以上に到達した場合に、スロットルバルブの異常として判定される。本実施形態においては、上記の異常信号がΔy及びi^(k+1)に基づいて検出されると共に、帰属確率に基づいて検出される可能性がある。したがって、Δy及びi^(k+1)に基づく判定結果の信頼性を上げることができる。
【００９１】
[実施形態４における具体的処理]
次に、図１１を参照して本実施形態においてＥＣＵが実行するスロットルバルブ異常判定ルーチンを説明する。尚、図１１中、ステップ１００〜ステップ１８０は、実施形態１と同一の処理が実行されるため説明を省略する。
【００９２】
ステップ１６０に続いて、帰属確率α1p1(zk+1;Φ)及びα2p2(zk+1;Φ)が演算される（ステップ５２０）。上述の通り、帰属確率演算部３２で、これら２つの帰属確率が演算される。このこれら２つの帰属確率は、ステップ５３０でＥＣＵにバッファリングされる。ステップ５３０に続いて、スロットルバルブの異常判定が実行される（ステップ５４０）。
【００９３】
次に、ステップ５４０で実行されるスロットルバルブの異常判定のサブルーチンを図１２のフローチャートを参照して説明する。尚、図１２中、ステップ２００〜２５０は、実施形態１と同一の処理が実行されるため説明を省略する。
【００９４】
ステップ２６０では、推定したセンサ電圧y^(k+1)と、実際のセンサ電圧y(k+1)とが比較される。上述の通り、センサ電圧y^(k+1)は、ステップ１６０又は１８０の処理を実行することによりＥＣＵにバッファリングされている。このため、このy^(k+1)とステップ２３０で読み込んだy(k+1)とが比較される。そして、これらの間の偏差が予め決められた閾値よりも小さい場合には、ステップ５５０に進む。そうでない場合には、ステップ２５０に進み、異常カウンタのカウンタ値が１増加される。
【００９５】
ステップ５５０では、開閉状態予測フラグi^(k+1)=1であるか否かが判定される。上述したように、ステップ１３０の処理を実行することにより開閉状態予測フラグi^(k+1)はＥＣＵにバッファリングされている。このため、本ステップでは、この開閉状態予測フラグi^(k+1)=1であるか否かが判定され、i^(k+1)=1の場合にはステップ５６０に進む。そうでない場合は、ステップ５７０に進む。
【００９６】
ステップ５６０では、帰属確率α1p1(zk+1;Φ)≧α2p2(zk+1;Φ)であるか否かが判定される。ステップ５６０は、ステップ５５０の処理の結果、スロットルバルブが開状態であると予測されている。そして、本ステップでα1p1(zk+1;Φ)≧α2p2(zk+1;Φ)であると判定された場合には、開状態である確率が高いことから、開閉状態予測フラグの結果と一致する。したがって、α1p1(zk+1;Φ)≧α2p2(zk+1;Φ)の場合には、ステップ５８０に進み、異常カウンタがリセットされる。一方、本ステップでα1p1(zk+1;Φ)＜α2p2(zk+1;Φ)であると判定された場合には、閉状態である確率が高いことから、開閉状態予測フラグの結果と相違する。このため、α1p1(zk+1;Φ)＜α2p2(zk+1;Φ)の場合には、ステップ２５０に進み、異常カウンタのカウンタ値が１増加される。
【００９７】
ステップ５７０では、帰属確率α1p1(zk+1;Φ)≦α2p2(zk+1;Φ)であるか否かが判定される。ステップ５７０は、ステップ５５０の処理の結果、スロットルバルブが閉状態であると予測されている。そして、本ステップでα1p1(zk+1;Φ)≦α2p2(zk+1;Φ)であると判定された場合には、開状態である確率が高いことから、開閉状態予測フラグの結果と一致する。したがって、α1p1(zk+1;Φ)≦α2p2(zk+1;Φ)の場合には、ステップ５８０に進み、異常カウンタがリセットされる。一方、本ステップでα1p1(zk+1;Φ)＞α2p2(zk+1;Φ)であると判定された場合には、開状態である確率が高いことから、開閉状態予測フラグの結果と相違する。このため、α1p1(zk+1;Φ)＞α2p2(zk+1;Φ)の場合には、ステップ２５０に進み、異常カウンタのカウンタ値が１増加される。
【００９８】
続いて、ステップ５９０では、異常カウンタのカウンタ値が予め定めた設定値（閾値）以上か否かが判定される。そして、カウンタ値が閾値以上の場合には、ステップ６１０に進み、スロットルバルブは異常と判定される。そうでない場合には、ステップ６００に進み、スロットルバルブは正常と判定される。
【００９９】
以上説明した通り、図１１及び１２に示すルーチンによれば、過去の時系列データに基づいて同定したＰＷＡＲＸモデルを用いて、将来のセンサ電圧を逐次推定することができる。また、図１１及び１２に示すルーチンによれば、ＰＷＡＲＸモデルの同定に際して推定した分離超平面のパラメータを用いて、開閉状態を予測できる。さらに、図１１及び１２に示すルーチンによれば、推定した将来のセンサ電圧を用いて定義した観測データベクトルzk+1を用いて、帰属確率を演算することができる。そして、将来のセンサ電圧及び開閉状態に基づいてスロットルバルブの異常判定を行うと共に、帰属確率に基づいてデータの信憑性を確認することにより、異常判定の精度を上げることができる。
【０１００】
尚、本実施形態は、実施形態３との組み合わせが可能である。すなわち、スロットルバルブの異常判定に際して、推定したセンサ電圧の推定値y^(k+1)、開閉状態予測フラグi^(k+1)、帰属確率α1p1(zk+1;Φ)、α2p2(zk+1;Φ)に加えて、分離超平面からの距離dkを用いてもよい。こうすることで、異常判定の精度を更に上げることができる。
【０１０１】
尚、上述した実施形態においては、ステップ５５０、５６０、２５０、５７０、５８０の処理を実行することにより前記第５の発明における「評価第３手段」が実現されている。
【図面の簡単な説明】
【０１０２】
【図１】実施形態１において、ＥＣＵがスロットルバルブの異常判定を行う診断装置として機能する場合の機能ブロック図である。
【図２】実施形態１において、ＥＣＵが実行するスロットルバルブ異常判定ルーチン説明するためのフローチャートである。
【図３】実施形態１において、ＥＣＵが実行するスロットルバルブ異常判定のサブルーチンを説明するためのフローチャートである。
【図４】実施形態２において、ＥＣＵがスロットルバルブの異常判定を行う診断装置として機能する場合の機能ブロック図である。
【図５】実施形態２において、ＥＣＵが実行するスロットルバルブ異常判定ルーチン説明するためのフローチャートである。
【図６】実施形態２において、ＥＣＵが実行するスロットルバルブ異常判定のサブルーチンを説明するためのフローチャートである。
【図７】実施形態３において、ＥＣＵがスロットルバルブの異常判定を行う診断装置として機能する場合の機能ブロック図である。
【図８】実施形態３において、ＥＣＵが実行するスロットルバルブ異常判定ルーチン説明するためのフローチャートである。
【図９】実施形態３において、ＥＣＵが実行するスロットルバルブ異常判定のサブルーチンを説明するためのフローチャートである。
【図１０】実施形態４において、ＥＣＵがスロットルバルブの異常判定を行う診断装置として機能する場合の機能ブロック図である。
【図１１】実施形態４において、ＥＣＵが実行するスロットルバルブ異常判定ルーチン説明するためのフローチャートである。
【図１２】実施形態４において、ＥＣＵが実行するスロットルバルブ異常判定のサブルーチンを説明するためのフローチャートである。
【符号の説明】
【０１０３】
１０    システム同定部
１２    データクラスタリング部
１４    分離超平面推定部
１６    システムパラメータ推定部
１８    センサ電圧推定部
２０    開閉状態推定部
２２、２６、３０、３４    異常判定部
２４、３２    帰属確率演算部
２８    距離演算部
【図１】
2010133384.tif 000009
【図２】
2010133384.tif 000010
【図３】
2010133384.tif 000011
【図４】
2010133384.tif 000012
【図５】
2010133384.tif 000013
【図６】
2010133384.tif 000014
【図７】
2010133384.tif 000015
【図８】
2010133384.tif 000016
【図９】
2010133384.tif 000017
【図１０】
2010133384.tif 000018
【図１１】
2010133384.tif 000019
【図１２】
2010133384.tif 000020