JP2017049067A - 機器状態判定装置、機器状態判定方法及びノッキング判定装置 - Google Patents

Abstract

【課題】検出した信号に対して判定すべきタイミングをより正確に推定することのできる機器状態判定装置、及び、機器状態判定方法、及び、内燃機関のノッキングを判定するノッキング判定装置を提供する。【解決手段】ノッキング判定装置１０は、周期性のある動作を伴う機器であるエンジンの周期性のある動作に基づく物理量を示す信号としてサイクル毎に取得された各取得信号の特定のタイミングの信号状態に基づいてエンジンの状態を判定する。ノッキング判定装置１０は、基準信号を取得する基準信号取得部６１と、判定信号を取得する判定信号取得部６２と、基準信号から平均と分散とに基づいて基準モデルを作成するモデル作成部３１と、判定信号の基準モデルに対する一致の度合いを示す指標を算出する指標算出部３５と、一致の度合いが最大であることを示す指標に基づいて同期値を決定するシフト量決定部３６とを備える。【選択図】図４

Description

本発明は、作動する機器の状態を判定する機器状態判定装置、及び、機器状態判定方法、及び、内燃機関に生じるノッキングを判定するノッキング判定装置に関する。

ガソリンエンジンなどの内燃機関における点火時期は、出力トルクの向上を目的として、ノックキングが発生しないクランク角度の範囲内において可能な限り進角されることが一般的である。そこで、内燃機関を設計したり、調整したりする過程で点火時期をクランク角度に適合させる工程では、ノックキングが発生しているか否かがノッキング判定装置によって判定される。こうしたノッキング判定装置の一例が特許文献１に記載されている。

特許文献１に記載のノッキング判定装置は、筒内圧を検出する筒内圧センサからの出力信号と、内燃機関が発する音を検出するマイクからの出力信号との相互相関値を計算し、相互相関値をノッキングの強度として取り扱っている。

特開２００３−３１４３４９号公報

特許文献１に記載のノッキング判定装置によれば、センサで検出された筒内圧と音の相互相関に基づいてノッキングの有無を判定している。
ところで、内燃機関の構造は様々であることから、ノッキングの有無を判定しようとする内燃機関に筒内圧センサが設置されているとは限らない。例えば、筒内圧を取得する筒内圧センサが設置されていない内燃機関では、検出された音と、筒内圧が上昇する点火タイミングとの正確な相関関係の取得ができない。このため、検出した音に対して、ノッキング判定を行うべきタイミングである点火タイミングが正確に得られず、ノッキング有無の判定精度が低下するおそれがある。この課題の解決には、検出した信号に対して判定すべきタイミングを正確に考慮することが望ましい。上記ノッキング有無の他に内燃機関の周期的な運転状態の判定、また、検出した信号に対して判定すべきタイミングを正確に考慮することが望ましい。また、内燃機関の他に周期的な動作を含む機器の周期的な状態の判定でも同様である。

本発明は、このような実情に鑑みてなされたものであってその目的は、検出した信号に対して判定すべきタイミングをより正確に推定することのできる機器状態判定装置、及び、機器状態判定方法、及び、内燃機関のノッキングを判定するノッキング判定装置を提供することにある。

以下、上記課題を解決するための手段及びその作用効果について説明する。
上記課題を解決する機器状態判定装置は、周期性のある動作を伴う機器の前記周期性のある動作に基づく物理量を示す信号としてサイクル毎に取得された各取得信号の特定のタイミングの信号状態に基づいて前記機器の状態を判定する機器状態判定装置であって、前記機器の基準となる条件での動作に基づく前記取得信号を前記サイクル毎に基準信号として取得する基準信号取得部と、前記機器の動作状態を判定するサイクルでの前記取得信号を判定信号として取得する判定信号取得部と、前記各基準信号の取得されたサイクル中の前記タイミングでの前記基準信号の集合に対する平均と分散とに基づいて基準モデルを作成するモデル作成部と、前記基準モデルの作成されたサイクル中の前記タイミングに対して、前記判定信号の取得されたサイクル中の前記タイミングを順次相対移動させると共に、前記相対移動毎の前記判定信号の前記基準モデルに対する一致の度合いを示す指標を算出する指標算出部と、前記一致の度合いが最大であることを示す指標に対応する相対移動に基づいて前記基準モデルに前記判定信号が同期する同期値を決定する同期値決定部と、を備え、前記同期値で移動させた前記タイミングでの前記判定信号に基づいて前記判定を行うことを要旨とする。

上記課題を解決する機器状態判定方法は、周期性のある動作を伴う機器の前記周期性のある動作に基づく物理量を示す信号としてサイクル毎に取得した各取得信号の特定のタイミングの信号状態に基づいて前記機器の状態を判定する機器状態判定装置で実行される機器状態判定方法であって、機器状態判定装置は、前記機器の基準となる条件での動作に基づく前記取得信号を前記サイクル毎の基準信号として取得する基準信号取得工程と、前記機器の動作状態を判定するサイクルでの前記取得信号を判定信号として取得する判定信号取得工程と、前記各基準信号の取得されたサイクル中の前記タイミングでの前記基準信号の集合に対する平均と分散とに基づいて基準モデルを作成するモデル作成工程と、前記基準モデルの作成されたサイクル中の前記タイミングに対して、前記判定信号の取得されたサイクル中の前記タイミングを順次相対移動させると共に、前記相対移動毎の前記判定信号の前記基準モデルに対する一致の度合いを示す指標を算出する指標算出工程と、前記一致の度合いが最大であることを示す指標に対応する相対移動に基づいて前記基準モデルに前記判定信号が同期する同期値を決定する同期値決定工程と、を備え、前記同期値で移動させた前記タイミングでの判定信号に基づいて前記判定を行うことを要旨とする。

上記課題を解決するノッキング判定装置は、内燃機関の燃焼サイクルに基づく物理量を示す信号としてサイクル毎に取得された各取得信号の特定のタイミングの信号状態に基づいて前記内燃機関のノッキングの有無を判定するノッキング判定装置であって、前記内燃機関の基準となる条件での動作に基づく前記取得信号を前記サイクル毎に基準信号として取得する基準信号取得部と、前記内燃機関の動作状態を判定するサイクルでの前記取得信号を判定信号として取得する判定信号取得部と、前記各基準信号の取得されたサイクル中の前記タイミングでの前記基準信号の集合に対する平均と分散とに基づいて基準モデルを作成するモデル作成部と、前記基準モデルの作成されたサイクル中の前記タイミングに対して、前記判定信号の取得されたサイクル中の前記タイミングを順次相対移動させると共に、前記相対移動毎の前記判定信号の前記基準モデルに対する一致の度合いを示す指標を算出する指標算出部と、前記一致の度合いが最大であることを示す指標に対応する相対移動に基づいて前記基準モデルに前記判定信号が同期する同期値を決定する同期値決定部と、を備え、前記同期値で移動させた前記タイミングでの前記判定信号に基づいて前記判定を行うことを要旨とする。

機器に周期的に生じる状態をサイクルごとに検出して取得信号とし、この取得信号に基づいて機器の状態を判定する。このとき、機器状態を判定するために得た取得信号からなる判定信号を、対比の基準として得た取得信号からなる基準信号と比較するとき、判定信号と基準信号との各サイクルのタイミングを一致（同期）させる必要がある。そこで、上記構成によれば、判定信号の基準信号に対するタイミングのずれが基準モデルに対する一致の度合いの高いときの相対移動距離が同期値として算出される。そして、この算出された同期値で判定信号のタイミングを移動させることで、比較する基準モデルとの関係においてサイクルが同期され、サイクル中のタイミングも一致するようになる。これにより、検出した信号に対して判定すべきタイミングをより正確に推定することができる。

また、信号に含まれる大きな変化などの特徴値に着目するのではなく、平均と分散に基づいて算出される基準モデルに対する一致の度合いを指標として算出する。このため基準モデルに含まれるノイズなどの「はずれ値」などの影響を相対的に低くすることができる。これによっても、検出した信号に対して判定すべきタイミングをより正確に推定することができる。

さらに、機器をノッキング判定装置とすれば、内燃機関から得られる各取得信号についてサイクル中のタイミングを同期させ、特定のタイミングにおけるノッキングの有無の状態を判定することができる。

上記機器状態判定装置において、前記指標算出部は、前記指標を前記基準モデルの分散が大きいタイミングにおける一致の度合いよりも分散が小さいタイミングにおける一致の度合いを重視して算出する。

上記ノッキング判定装置において、前記指標算出部は、前記指標を前記基準モデルの分散が大きいタイミングにおける一致の度合いよりも分散が小さいタイミングにおける一致の度合いを重視して算出する。

このような構成によれば、分散の小さい部分、つまり変動が小さく値が安定しているタイミングを重視して指標が算出されることから検出した信号について同期値の算出精度の向上が図られる。

上記機器状態判定装置において、前記指標算出部は、前記判定信号の取得されたサイクル中のタイミングを前記判定信号が取得されたときに定まるタイミングを含む所定の範囲内で順次相対移動させる。

上記ノッキング判定装置において、前記指標算出部は、前記判定信号の取得されたサイクル中のタイミングを前記判定信号が取得されたときに定まるタイミングを含む所定の範囲内で順次相対移動させる。

ある程度正常に作動している機器であれば、基準モデルに対して判定信号の取得されたサイクル中のタイミングが大きくずれることは少ない。よって、上記構成によれば、相対移動される範囲が、判定信号が取得されたときに定まるタイミングを含む所定の範囲内とされることから、相対移動にかかる負荷の軽減や、処理に要する時間の削減が図られる。

上記機器状態判定装置において、前記取得信号には、取得されたサイクル中のタイミングの目安となる同期信号が対応付けられており、前記指標算出部は、前記判定信号の取得されたサイクル中のタイミングを前記判定信号の前記同期信号のタイミングと前記基準モデルの前記同期信号のタイミングとを一致させたタイミングを含む所定の範囲で順次相対移動させる。

上記機器状態判定装置において、前記機器は内燃機関であり、前記サイクルは前記内燃機関の燃焼サイクルであり、前記同期信号は前記内燃機関の点火プラグへ電力が供給されることで得られる信号である。

上記ノッキング判定装置において、前記取得信号には、取得されたサイクル中のタイミングの目安となるイグニッションパルスに対応する信号が同期信号として対応付けられており、前記指標算出部は、前記判定信号の取得されたサイクル中のタイミングを前記判定信号の前記同期信号のタイミングと前記基準モデルの前記同期信号のタイミングとを一致させたタイミングを含む所定の範囲で順次相対移動させる。

取得信号に同期信号が含まれている場合、同期信号を目安に判定信号と基準モデルとを比較することが可能になるが、同期信号の精度が高くない場合もある。よって、上記構成によれば、同期信号を含む範囲に相対移動する範囲を設定することで、同期信号の精度が高くない場合でも、基準モデルに対して判定信号を判定すべきタイミングをより正確に推定することができる。

また、機器が内燃機関であれば、内燃機関から得られる各取得信号について取得されたサイクル中のタイミングを同期させ、特定のタイミングにおける運転状態を判定することができる。

上記機器状態判定装置において、前記モデル作成部は、前記基準信号のうち所定のタイミングからなる範囲に対応する基準モデルを作成し、前記指標算出部は、前記判定信号のうち前記所定のタイミングからなる範囲について前記基準モデルに対する前記指標を算出する。

上記ノッキング判定装置において、前記モデル作成部は、前記基準信号から前記燃焼サイクル中の所定のタイミング範囲に対応する基準モデルを作成し、前記指標算出部は、前記判定信号のうち前記所定のタイミング範囲について前記基準モデルに対する前記指標を算出する。

このような構成によれば、同期値を得るために必要とされる基準信号や判定信号の範囲が所定範囲に制限されるため、基準モデルの作成や指標の算出にかかる負荷の軽減や、処理に要する時間の削減が図られる。

上記機器状態判定装置において、前記指標算出部は、前記指標を前記基準モデルの平均値のとき最大値となり、前記平均値から離れることに応じて小さくなるように算出し、前記同期値決定部は、前記指標の最大値を前記一致の度合いが最大であるものとする。

上記ノッキング判定装置において、前記指標算出部は、前記指標を前記基準モデルの平均値のとき最大値となり、前記平均値から離れることに応じて小さくなるように算出し、前記同期値決定部は、前記指標の最大値を前記一致の度合いが最大であるものとする。

このような構成によれば、基準モデルが、平均値とそこからの離間、例えば分散に応じて定められる。これにより、一致の度合いによる指標の算出が好適に行えるようになる。

このような機器状態判定装置、機器状態判定方法、及びノッキング判定装置によれば、検出した信号に対して判定すべきタイミングをより正確に推定することが可能になる。

ノッキング判定装置を具体化した一実施形態について、ノッキング判定装置によりエンジンのノッキングを判定するときの各種機器との接続態様の一例を示すブロック図。 同実施形態のノッキング判定装置で取得される複数の音圧信号をクランク軸の角度に同期させたときの一例を示すグラフ。 同実施形態のノッキング判定装置で取得される複数の音圧信号を点火時期に同期させたときの一例を示すグラフ。 同実施形態のノッキング判定装置の概略構成を示すブロック図。 同実施形態のノッキング判定装置で取得される２つの音圧信号を点火時期を同期させて重ねた態様の一例を示すグラフ。 同実施形態のノッキング判定装置が行う基準モデル作成処理について各工程を実行の順に示すフローチャート。 図５に示す２つの音圧信号にフィルタ処理を施したときの一例を示すグラフ。 図７に示す音圧信号のうち基準信号である一つを絶対値化したときの一例を示すグラフ。 図８に示す基準信号の複数について平均化したときの一例を示すグラフ。 図９に示す基準信号の複数について、（ａ）はそれらの分散の一例を示すグラフ、（ｂ）は指標算出に使う確率密度関数の一例を説明するグラフ。 図８に示す平均化されたグラフと、この平均値に対して確率毎の分散とを重ね合わせて示したグラフ。 同実施形態のノッキング判定装置が行う確率最大化処理について各工程を実行の順に示すフローチャート。 同実施形態のノッキング判定装置で取得される音圧信号のうち状態が判定される信号である判定信号の一例を示すグラフ。 同実施形態のノッキング判定装置が行うスコア算出処理について各工程を実行の順に示すフローチャート。 同実施形態のノッキング判定装置が判定信号のタイミングを複数の間隔で変化させたときのスコア算出処理で算出されるスコアの一例を示すグラフ。 同実施形態のノッキング判定装置が行うシフト量決定処理について各工程を実行の順に示すフローチャート。 ノッキング判定装置を具体化したその他の実施形態について、ノッキング判定装置によりエンジンのノッキングを判定するときの各種機器との接続態様の一例を示すブロック図。

以下、図１〜図１６を参照し、機器状態判定装置、機器状態判定方法及びノッキング判定装置を具体化した一実施形態について説明する。内燃機関の一例である車両用のエンジン１の試験では、例えば点火タイミングの進角量の調整等の試験もしくは調整された点火タイミングの確認試験が行われる。これらの点火タイミングの試験においてノッキング判定装置が用いられる。ノッキングとは、エンジン１の気筒内において発生した異常燃焼により生じる圧力変動（衝撃波）が気筒の固有振動数で増幅されてエンジン１に大きな振動が発生する現象である。つまり、エンジン１には気筒内での燃料燃焼の際に発生する圧力変動に基づく振動が生じ、この振動を物理量として取得することができる。

本実施形態では、エンジン１の複数のサイクルにおける振動に基づく値を相互比較することよりノッキングが検出される。ここで、複数のサイクルにおける振動に基づく値を比較するためには、比較する各振動のサイクルの周期的なタイミングを相互に一致（同期）させる必要がある。具体的には、各サイクルにおいてクランク角度が同じタイミングの値同士を比較する。そこで、本実施形態では、ノッキングを判定する処理に先立ち、比較する複数のサイクルの振動を測定して得られた各取得信号の周期的なタイミングを同期させる処理を行う。本実施形態では、タイミングを同期させる技術の一例として、点火プラグ７に供給される電流を電流センサ３で取得して得たイグニションパルス（ＩＰ）の立ち上がりを示す信号（ＩＰ信号）を利用して、サイクル毎の取得信号のタイミングを一致（同期）させる技術について説明する。なお、以下では、エンジン１の回転速度は一定であることから各サイクルの長さが同じであり、取得されたサイクル中のタイミングのうちの一箇所を一致（同期）させることで各取得信号が同期する。

ところで、エンジン１に周期的に生じる振動、主に機械的な振動は、クランク角度に応じて変化する機器の作動状態に起因して生じるものが多い、すなわちクランク角度と相関性を有することが多い。また、ＩＰ信号は、点火プラグ７に接続される電線から容易に取得できるが、燃焼制御において点火タイミングが進角されることによって、クランク角度に対するタイミングが変化する。そのため、ＩＰ信号は取得が容易であるが、クランク角度ではなく、ＩＰ信号で取得信号を同期させるには、各サイクルでの取得信号のタイミング同期処理が必要になる。

以下では、サイクル毎の取得信号を同期させる処理の一例について詳細に説明する。なお、ノッキングに関する処理の一例についても簡単に説明する。
図１に示すように、ノッキング判定装置１０は、データ収集装置５を介してエンジン１に関する物理量の変化に基づく情報を取得する。また、ノッキング判定装置１０は、複数の取得信号に基づく値の比較により、ノッキング判定の対象とした取得信号（判定信号）がノッキングを示す信号であるか否かを判定するための演算を実行して、その演算結果をモニタ６等に出力する。

図１に示すように、試験対象のエンジン１は、車両１００に搭載されている。エンジン１は、例えば、「吸入」、「圧縮」、「燃焼」及び「排気」の４行程を１サイクルとする４ストロークエンジンである。なお、試験対象のエンジン１は、車両１００に搭載されない単独の状態でもよい。エンジン１には、点火装置８に接続された点火プラグ７が設けられ、点火装置８にはエンジンＥＣＵ２が接続されている。点火プラグ７はその先端がエンジン１の気筒内に配置され、その基端が電線を介して点火装置８に電気的に接続される。点火装置８は、エンジンＥＣＵ２から点火タイミング毎に受ける点火指示に基づき点火プラグ７に電流を供給する。エンジンＥＣＵ２は、エンジン１の駆動を制御するエンジンコントロールユニットであって、ＣＰＵ、ＲＯＭ、ＲＡＭ、その他の記憶装置を有するコンピュータを含み構成されている。エンジンＥＣＵ２は、ＲＯＭや記憶装置に記憶されているプログラムをＣＰＵで演算処理することで、エンジン１の駆動の制御に必要な各種情報をエンジンＥＣＵ２の外部から取得しながら演算した点火指示などの出力によりエンジン１の駆動を制御する。エンジン１は、各サイクルにおいて同じクランク角度では、通常、同様の機械的な振動を生じる。よって、サイクル毎の取得信号における機械的な振動のタイミングを一致させるように、比較する取得信号の各タイミングをシフトさせることによって、サイクル毎の取得信号のタイミングが一致（同期）する。

エンジン１の近くには、音圧センサ４が設置されている。音圧センサ４は、エンジン１から発生する音を検出し、この検出した音に基づく音圧信号をデータ収集装置５に出力する。詳述すると、音圧センサ４は、エンジン１に発生する圧力変動に基づく物理量の一例である音圧を検出し、検出された音圧の大きさを示す音圧信号を生成する。つまり、音圧センサ４から出力される音圧信号には、クランク角度に相関性のある動作に基づく振動である、周期性を有する動作に基づく音圧が取得される。例えば、周期性を有する動作に基づく音圧としては、ピストンの上死点や下死点に対応する音、吸気弁や排気弁の開閉弁に対応する音、燃料の燃焼に関連して生じる音などが挙げられる。つまり、音圧センサ４から出力される音圧信号には、エンジン１の周期性を有する動作に基づく音が含まれる。また、エンジン１にノッキングが発生しているとき、音圧センサ４から出力される音圧信号には、エンジン１の周期性を有する動作に基づく音に加えて、ノッキングに相関のある音も含まれる。

点火プラグ７と点火装置８とを接続する電線には、点火タイミングにおいて点火プラグ７に供給される電流（イグニッションパルス）の有無を検出する電流センサ３が設けられている。また、電流センサ３は、電流の検出結果をデータ収集装置５に出力する。よって、データ収集装置５には、電流センサ３による電流の検出結果が、通常、１回の燃焼サイクルにおいて１回、エンジン１の点火タイミングとして入力される。

データ収集装置５は、音圧センサ４からの音圧信号を入力してＡ／Ｄ変換する。また、データ収集装置５は、電流センサ３により検出された電流をＩＰ信号として入力する。そして、データ収集装置５は、連続入力される音圧信号から、ＩＰ信号に基づいてエンジン１の１サイクル分取得し、これをサイクル毎の取得信号とする。よって、データ収集装置５は、単位時間当たりの回転速度に応じた数、回転速度が３０００ｒ／ｍｉｎであれば一分間に１５００個の取得信号を作成する。また、データ収集装置５は、作成した取得信号にＩＰ信号の情報を関連付ける。そして、データ収集装置５は、ＩＰ信号が関連付けられた取得信号をノッキング判定装置１０に出力する。なお、データ収集装置５は、生成した取得信号を一時的に保持したり、一旦蓄えてからノッキング判定装置１０に出力してもよい。

ところで、データ収集装置５は、エンジン１に常設ではなく、ノッキングを測定する必要の生じたエンジン１に対して随時取り付けられる。そのため、データ収集装置５は、エンジンＥＣＵ２に保持される情報を容易に取得できるものではない。データ収集装置５は、エンジンＥＣＵ２から情報を取得する場合、エンジンＥＣＵ２に適合する接続装置などを準備し、この装置を介して情報を収集する必要がある。一方、エンジン１の外部に露出されている点火プラグ７に接続される電線に流れる電流をＩＰ信号として取得することは、電線に測定用のクランプを取り付けることなどで比較的容易に行える。

なお、以降の説明では、音圧信号は、「吸入」、「圧縮」、「燃焼」及び「排気」の４行程からなる１サイクルのうちの一部である「圧縮」行程の一部の範囲もしくは全部の範囲を切り出したものとして説明をする。つまり、１サイクルの一部からなる音圧信号に基づいて、サイクル毎の音圧信号のタイミングを一致（同期）させる処理を行う。例えば、「燃焼」行程で生じるノッキングの有無を判定するとき、「燃焼」行程の直前の「圧縮」行程に基づいて音圧信号のタイミングを一致（同期）させる。なお、タイミングの同期に用いられる音圧信号は、「吸入」、「圧縮」、「燃焼」及び「排気」の４行程からなる１サイクルのうちの一部を切り出したものであっても、全部であってもよい。また、１サイクルのうち切り出される一部は、「圧縮」行程以外の行程、「吸入」、「燃焼」及び「排気」の各行程であってもよいし、いずれか複数の行程を含んでいてもよい。

データ収集装置５に収集されるサイクル毎の音圧信号とＩＰ信号との関係について、図２及び図３を参照して説明する。
まず、図２を参照し、サイクル毎の音圧信号と、クランク角度に同期しているクランク軸同期信号Ｔｃとの関係について説明する。

図２に示すように、音圧センサ４から進角量の異なる音圧信号が取得される。例えば、点火タイミングの進角量が「０」の音圧信号Ｌａと、第１進角量の音圧信号Ｌｂと、第２進角量の音圧信号Ｌｃとが取得される。このとき、サイクル毎の各音圧信号Ｌａ，Ｌｂ，Ｌｃのタイミングを、エンジンＥＣＵ２から取得したクランク軸同期信号Ｔｃに基づいて一致（同期）させると、所定のクランク角度において生じる特定の音Ｎ１の発生タイミングＴ０が一致（同期）する。一方、例えば、電流センサ３で取得されたＩＰ信号から推定される推定同期信号Ｔｉｐは、進角量の影響により音圧信号毎に取得されたサイクル中におけるタイミングが相違している。なお、「所定のクランク角度において生じる特定の音」とは、サイクル毎に、同じクランク角度でエンジン１から生じる相対的に大きいなど特徴的な音であって、その音の発生原因は特に限定されない。よって、エンジン１の個体差、測定時期、又は、回転速度や負荷の有無などの各種運転条件によって適切な音を、同期処理用の特定の音として任意に選択することができる。

図３は、音圧センサ４から取得した、各音圧信号Ｌａ，Ｌｂ，Ｌｃを、ＩＰ信号に基づき推定された推定同期信号Ｔｉｐで一致（同期）させた例を示す。推定同期信号Ｔｉｐは進角量の影響を受けてクランク軸同期信号のタイミングがタイミングＴｃ１，Ｔｃ２，Ｔｃ３のようにずれるため、同一のクランク角度で生じている特定の音Ｎ１が相違する発生タイミングＴ０，Ｔ１，Ｔ２で得られる。すなわち、クランク軸同期信号Ｔｃの取得に代えて、ＩＰ信号に基づく推定同期信号Ｔｉｐで各音圧信号Ｌａ，Ｌｂ，Ｌｃを比較しようとするとき、サイクル毎のタイミングに進角量に起因するずれが生じる。

例えば、図３に示すように、音圧信号Ｌａと音圧信号Ｌｂとの間では特定の音Ｎ１のタイミングが「Ｔ０−Ｔ１」だけずれる。また、音圧信号Ｌａと音圧信号Ｌｃとの間では特定の音Ｎ１のタイミングが「Ｔ０−Ｔ２」だけずれる。つまり、推定同期信号Ｔｉｐでタイミングを同期させるためには、音圧信号Ｌａに対する各音圧信号Ｌｂ，Ｌｃのタイミングのずれを、各音圧信号Ｌｂ，Ｌｃの同期値としての移動量（シフト量）として算出する。そして、この算出した移動量で各音圧信号Ｌｂ，Ｌｃのサイクルのタイミングを補正することで、推定同期信号Ｔｉｐに基づいてタイミングが特定されていた各音圧信号Ｌａ，Ｌｂ，Ｌｃの各タイミングを一致（同期）させることができるようになる。

次に、ノッキング判定装置１０について説明する。
図４に示すように、ノッキング判定装置１０は、データ収集装置５から入力した取得信号を同期させる同期処理を実行し、同期した各取得信号に基づく値の比較によってノッキング判定処理を行う。

ノッキング判定装置１０は、ＣＰＵ、ＲＯＭ、ＲＡＭ、その他の記憶装置を有するコンピュータを含み構成されている。ノッキング判定装置１０は、ＲＯＭや記憶装置に記憶されているプログラムをＣＰＵで演算処理することにより、所定の処理を行う。本実施形態では、ＲＯＭや記憶装置に記憶されているプログラムをＣＰＵで演算処理することにより、取得信号の同期処理やノッキング判定処理に関する各演算処理を実行する。なお、ノッキング判定装置１０は、取得信号の同期方法やノッキング検出方法を実行するプログラムを有するパーソナルコンピュータ（ＰＣ）等であってもよい。

ノッキング判定装置１０は、データ収集装置５から出力される取得信号と、その取得信号に関連付けられたＩＰ信号とを取得する。取得信号は、エンジン１に発生する機械的な振動に基づく物理量に相当する。取得信号には、１サイクル中の一部区間又は全区間にスペクトル（周波数成分）が含まれている。なお、ノッキング音は「燃焼」行程を含む区間に比較的大きな音として含まれ、またスペクトルにも特徴量として現れる。

図４に示すように、ノッキング判定装置１０は、取得信号などを保持する記憶部２０と、同期処理に必要な信号を記憶部２０に保持させる管理部６０と、記憶部２０から同期処理に必要な信号を選択する信号同期部３０とを備える。また、ノッキング判定装置１０は、記憶部２０からノッキング判定処理に必要な信号を選択する信号選択部３９と、取得信号のスペクトルを算出するスペクトル算出部４０と、取得信号のバイスペクトルを算出するバイスペクトル算出部４１と、判定信号がノッキングを示す信号であるか否かを判定するノッキング判定部５０とを備えている。

記憶部２０は、ノッキング判定装置１０を構成する記憶装置の一部または全部から構成されており、ＩＰ信号や、取得信号より選択される基準信号や判定信号などを記憶したり、削除させたりできるようになっている。

記憶部２０には、取得信号から選択される基準信号に関する記憶領域である基準信号領域２０１と、基準信号に基づいて算出される基準モデルに関する記憶領域である基準モデル領域２０２とが設けられている。また、記憶部２０には、ノッキング判定を行うものとして取得信号から選択される判定信号に関する記憶領域である判定信号領域２０３が設けられている。さらに、記憶部２０には、基準モデルに対する判定信号の一致の度合いに関する指標としての指標値の記憶領域である指標値領域２０４と、基準モデルに判定信号を同期させるための移動量としてのシフト量に関する記憶領域であるシフト量領域２０５とが設けられている。また、記憶部２０には、ノッキング判定処理に必要な各種の信号、例えば、判定信号に基づく値に対し基準となる対象信号に基づく値に関する値なども保持される。

管理部６０は、ノッキング判定装置１０が入力した取得信号にＩＰ信号を関連付けるとともに、所定の選択を行い記憶部２０に記憶させる。
管理部６０は、取得信号から基準モデルを作成するための基準信号を取得する基準信号取得部６１と、取得信号から判定信号を取得する判定信号取得部６２と、記憶部２０に記憶されている基準信号や判定信号等の整理を行うデータ整理部６３とを備える。

基準信号取得部６１は、基準信号取得条件が成立していることを条件にノッキング判定装置１０に入力された取得信号を基準信号として取得する（基準信号取得工程）。つまり、基準信号取得部６１は、エンジン１の基準となる条件での動作に基づく取得信号をサイクル毎の基準信号として取得する。基準信号取得条件は、基準信号の取得が指示されていること、かつ、基準信号の取得に適した運転条件（以下、基準運転条件）でエンジン１が運転されていることを条件に含む。例えば、基準信号取得は、外部からの指示やプログラムなどにより指示される。また例えば、基準運転条件は、点火タイミングの進角量が「０」であり、回転速度が所定の回転速度で一定であり、負荷条件も所定の負荷条件で一定の条件である。なお、基準運転条件は、基準モデルが生成可能であれば、１つの条件で定められるものであっても、複数の条件が組み合わされて定められるものであってもよい。

そして、基準信号取得部６１は、ＩＰ信号に基づいて特定されたサイクル毎の取得信号をそれぞれ基準信号として取得（選択）し、この取得した基準信号を記憶部２０の基準信号領域２０１に保持させる。

判定信号取得部６２は、判定信号取得条件が成立していることを条件にノッキング判定装置１０に入力された取得信号を判定信号として取得する（判定信号取得工程）。判定信号取得部６２は、エンジン１の動作状態を判定するサイクルでの取得信号を判定信号として取得する。判定信号取得条件は、信号同期処理が必要とされていること、かつ、エンジン１の点火タイミングが変化することを条件に含む場合を例示する。同期処理が必要とされる条件は、例えば、ノッキング判定処理が指示されていることである。この指示は、外部からの指示、プログラムによる指示どちらでもよい。また、基準運転条件に対して、エンジン１の点火タイミングは変化するが、回転速度と負荷条件は変化しない場合を例示する。

そして、判定信号取得部６２は、ＩＰ信号に基づいて特定された１サイクル分の取得信号を判定信号として取得（選択）し、この取得した判定信号を記憶部２０の判定信号領域２０３に保持させる。

データ整理部６３は、記憶部２０に記憶されている基準信号、基準モデル、判定信号、指標値及びシフト量などを定められた条件に応じて整理する。例えば、判定信号、スコア値及びシフト量はそれぞれ、新たな判定信号の設定に応じて削除する。また例えば、新たな基準信号の取得が指示されると、基準信号、基準モデル、判定信号、指標値及びシフト量を削除する。

図５を参照し、基準信号取得部６１と判定信号取得部６２とが取得する基準信号Ｌ１と判定信号Ｌ２との関係について説明する。
基準信号取得部６１は、進角量が「０」である基準信号Ｌ１を複数取得し、判定信号取得部６２は、進角量が所定の値に変更されている判定信号Ｌ２を１つ取得する。このとき、基準信号Ｌ１も判定信号Ｌ２もＩＰ信号に基づいて判定されるサイクル毎に取得信号として取得されているため、進角量の相違に応じて各サイクルの各タイミングにずれが含まれる。

図４に示す信号同期部３０は、基準信号Ｌ１に基づき基準モデルを作成する。また、信号同期部３０は、作成した基準モデルと判定信号Ｌ２を比較し、判定信号Ｌ２の基準モデルに対する移動量（シフト量）を算出する。そこで、信号同期部３０は、基準モデルを作成するモデル作成部３１と、基準モデルと判定信号Ｌ２との一致の度合いである指標（指標値）を算出する指標算出部３５と、指標に基づいて移動量（シフト量）を決定する同期値決定部としてのシフト量決定部３６とを備える。

続いて、ノッキング判定処理を行う構成の一例について簡単に説明する。
信号選択部３９は、対象信号及び判定信号を選択する所定の条件に基づいて、取得信号から対象信号や判定信号を選択する。また、信号選択部３９は、対象信号や判定信号に含まれる音圧信号の決定された角度範囲（切出角度範囲）内の音圧信号を切り出す。

スペクトル算出部４０は、信号選択部３９により選択された対象信号や判定信号に対して離散フーリエ変換を行い、周波数成分（スペクトル）を算出する。離散フーリエ変換は、例えば、離散フーリエ変換を高速に計算する高速フーリエ変換（ＦＦＴ：Ｆａｓｔ Ｆｏｕｒｉｅｒ Ｔｒａｎｓｆｏｒｍ）により行われる。

バイスペクトル算出部４１は、スペクトル算出部４０が算出した周波数成分からスペクトルに基づく値としてのバイスペクトルを算出する。ここで、時刻ｔにおける信号選択部３９が選択した音圧信号をｘ（ｔ）とする。この音圧信号ｘ（ｔ）にＦＦＴ処理を施すと、周波数ｆを変数とするフーリエ係数Ｘ（ｆ）が算出される。このフーリエ係数Ｘ（ｆ）に基づくバイスペクトルＢ（ｆ_１，ｆ_２）は下記の式（１）で示される。バイスペクトルＢ（ｆ_１，ｆ_２）は、スペクトルの３乗積である。

ここで、ｆ_１，ｆ_２は周波数、＊は複素共役を表している。
ノッキング判定部５０は、判定信号がノッキングを示す信号であるか否かを判定する。
ところでノッキングは、気筒内の一種の共振現象であり、ノッキングが発生すると気筒内に衝撃波が発生し、その衝撃波に起因してエンジンから発せられる振動や音には、その気筒の、複数の固有振動数（共振周波数）の成分とそれらの高調波成分（特に２倍の周波数成分）が含まれる。バイスペクトルでｆ_１＝ｆ_２のバイスペクトルは、各周波数における２次の調和成分の大きさを反映している。

エンジンの固有振動数（共振周波数）ｆ_ｐ，ｑは、下記の式（２）で示される。

ただし、Ｃは音速、Ｄはボア径、ｐ，ｑは振動モード、Ρ_ｐ，ｑは振動モード（ｐ，ｑ）における定数である。例えば、振動モード（ｐ，ｑ）が（１，０）モードのときＰ_１，０＝１．８４１である。つまりノッキング判定処理では、取得信号の全周波数領域を対象にするのではなく、各振動モードの固有振動数とその周辺や高調波とその周辺を監視することでノッキング有無の判定が適切になされるようになる。

ノッキング判定部５０は、複数の対象信号の平均と分散共分散行列を算出する統計パラメータ算出部５１と、算出された統計パラメータを用いて対象信号や判定信号を正規化する正規化部５２とを備える。また、ノッキング判定部５０は、正規化された対象信号に基づいてノッキング判定用の判定値を算出する判定値算出部５３と、ノッキング判定用の判定値と正規化した判定信号を比較する比較部５４とを備える。

まず、ノッキング判定部５０は、バイスペクトル算出部４１で算出されたバイスペクトルＢ（ｆ_１，ｆ_２）を正規化する。
統計パラメータ算出部５１は、全ての対象信号のバイスペクトルの分布が多次元確率分布であると仮定して、全ての対象信号のバイスペクトルのばらつきによる平均と分散共分散とを統計パラメータとして推定（算出）する。例えば、複数の対象信号の各々で算出した振動モード別の特徴量から絶対値を算出し、複数の対象信号に対する平均値及び分散共分散行列を算出する。この平均値及び分散共分散行列を使用することで、複数の対象信号に対するマハラノビス距離を算出できる。

正規化部５２は、対象信号や判定信号のバイスペクトルについてマハラノビス距離を算出する。マハラノビス距離の算出に用いられる分散共分散行列は、バイスペクトルの各周波数成分の平均からの偏差の積の平均値である共分散を配列した行列であることから、マハラノビス距離を算出することによって多次元において正規化が行われる。

正規化部５２は、判定値算出処理時は、複数の対象信号のバイスペクトルの全てについて複数の対象信号の平均値及び分散共分散行列を用いマハラノビス距離を算出する。そして、これら算出されたマハラノビス距離に基づいて、ノッキング判定に用いる判定値が算出される。

正規化部５２は、判定処理時には、判定信号のバイスペクトルについて、複数の対象信号に対する平均値及び分散共分散行列を用いマハラノビス距離を算出する。
判定値算出部５３は、対象信号のバイスペクトルから得られた正規化後の集合と、判定信号のバイスペクトルから得られた正規化後の値との乖離の度合いを判定する判定値を算出する。ここで、複数の対象信号のバイスペクトルから得られた正規化後の集合は、各対象信号のバイスペクトルについて算出したマハラノビス距離の集合である。また、判定信号のバイスペクトルから得られた正規化後の値は、判定信号のバイスペクトルについて算出したマハラノビス距離の値である。

判定値算出部５３は、複数の対象信号に対して得られた各マハラノビス距離に対してマージンを有するように判定値を算出する。こうして算出された判定値が、判定信号がノッキングを示す信号であるか否かの判定に用いられる。

比較部５４は、判定信号に基づくマハラノビス距離を判定値算出部５３で算出された判定値と比較することで、判定信号がノッキングを示す信号で有るか否かを判定する。比較部５４は、判定信号に基づくマハラノビス距離の値が判定値よりも大きければ、判定信号がノッキングを示す信号であると判定する。逆に、比較部５４は、判定信号に基づくマハラノビス距離の値が判定値以下であれば、判定信号がノッキングを示す信号ではないと判定する。

次に、図６〜図１６を参照し、ノッキング判定装置１０における基準モデルと判定信号との取得されたサイクル中のタイミング同期処理について説明する。なお、本実施形態では、エンジン１の回転速度が一定に保たれている場合について例示する。

まず、図６〜図１１を参照し、ノッキング判定装置１０による基準モデルの作成処理（モデル作成工程）について説明する。この基準モデルの作成処理は、ノッキング判定処理の開始に先立ち行われる。

図６〜図１１に示すように、基準モデルの作成処理が開始されると、信号同期部３０のモデル作成部３１は、記憶部２０から複数サイクル分の基準信号を取得し（ステップＳ１０）、取得した基準信号から同期処理に利用するタイミングの範囲、例えば「圧縮」行程に対応する範囲Ｄ１（図５参照）を切り出す処理を行う（ステップＳ１１）。信号同期部３０のモデル作成部３１は、基準信号Ｌ１（図５参照）に対して所定の前処理を行う（ステップＳ１２）。例えば、前処理として、周波数フィルタ処理と値の絶対値化を行う。周波数フィルタ処理では、基準信号Ｌ１を、高周波数を透過させるハイパスフィルタを通過させて信号Ｌ１１（例えば、図７の範囲Ｄ１の区間）を得る。また、値の絶対値化では、通過させた信号Ｌ１１の振幅の中心を「０」とし、負側の振幅を正側に配置させて絶対値信号Ｌ１２（図８参照）を得る。そして、ここまでの処理が複数の基準信号Ｌ１に対して行われる。なお、範囲Ｄ１を切り出す処理（ステップＳ１１）は、前処理（ステップＳ１２）の後などに行ってもよい。

前処理が終了すると、モデル作成部３１は、複数の基準信号Ｌ１から得られた複数の絶対値信号Ｌ１２に基づいて平均（μ）を算出する（ステップＳ１３）とともに、分散（σ）を算出する（ステップＳ１４）。

図９に示すように、平均は、各絶対値信号Ｌ１２の各タイミングの値に基づいて、範囲Ｄ１に含まれるタイミング毎に算出される平均値Ｌ１３である。つまり、平均値Ｌ１３は、特定のタイミングにおける各絶対値信号Ｌ１２の値の集合に対して算出される平均値を、範囲Ｄ１の区間のそれぞれのタイミングに対して算出したものである。

図１０に示すように、分散は、各絶対値信号Ｌ１２の各タイミングの値に基づいて、範囲Ｄ１に含まれるタイミング毎に算出される分散Ｌ１４（図１０（ａ）参照）である。つまり、分散Ｌ１４は、特定のタイミングにおける各絶対値信号Ｌ１２の値の集合に対して算出される確率分布（図１０（ｂ）参照）におけるσの大きさを、範囲Ｄ１の区間のそれぞれのタイミングに対して算出したものである。例えば、所定の分散ｋσ（ｋ：任意の定数）は、ばらつき（標準偏差）の小さいタイミングでは小さく、ばらつき（標準偏差）の大きいタイミングでは大きくなる。

図１１に示すように、平均と分散とが算出されると、モデル作成部３１は、平均と分散により基準モデルを作成する（ステップＳ１５）。基準モデルは、平均値Ｌ１３と、標準偏差Ｌ１４ａ，Ｌ１４ｂ，Ｌ１４ｃなどから構成される。そして、この基準モデルが、判定信号の基準モデルに対するタイミングのずれを示す移動量（シフト量）の算出に用いられる。

基準モデルの一例は、以下の式（３）で表される。

但し、Ｐは同時確率、ｐは正規分布の確率密度関数、ｘは判定信号、μは基準モデルの平均、σは基準モデルの分散、ｔは時刻。また、Ｘはｘ_ｔを要素とする行列を、Σはσ_ｔを要素とする行列を、Ｍはμ_ｔを要素とするベクトルをそれぞれ表す。

続いて、図１２〜図１６を参照し、ノッキング判定装置１０による確率最大化処理について説明する。この確率最大化処理は、基準モデル作成処理などに続いて開始される。
ところで本実施形態では、基準モデルの分散が小さい部分を重視してスコアリングを行う。例えば、図１１において、範囲Ｄ１のうち、分散の大きい部分に対する比較結果よりも、分散の小さい部分に対する比較結果を重視したスコアリングが式（３）により行われる。すなわち、上記式（３）を用いてスコアリングを行うことで分散が小さい部分を重視したスコアリングが行えるようになる。換言すると、基準モデルに含まれる分散の大きい部分の影響を小さくすることができる。これに対し、従来の技術、例えば、相互相関によるスコアリングは、大きな値が重視されるスコアリングであるため、大きな値が重視されて外乱に弱く、強度の変動が大きい音に対して一致判定の精度が高められない。また例えば、二乗誤差によるスコアリングも、大きな値の一致度を重視したスコアリングであるため、基準モデルに含まれるはずれ値の影響が大きくなり、やはり一致判定の精度が高められない。

確率最大化処理が開始されると、信号同期部３０の確率最大化部３４は、判定信号Ｌ２（図５参照）を取得し（ステップＳ２０）、判定信号のうち同期処理に用いられるタイミングの部分を切り出し（ステップＳ２１）、前処理を行う（ステップＳ２２）。判定信号は、記憶部２０の判定信号領域２０３から取得される。なお、上述した基準モデル作成処理に対し、比較するタイミングの切り出しは、処理の対象が基準信号Ｌ１から判定信号Ｌ２に変更されたものである。また、前処理は、処理の対象が基準信号Ｌ１から判定信号Ｌ２に変更されたものであり、例えば、前処理として、周波数フィルタ処理により信号Ｌ２１（図７参照）を得て、値の絶対値化により絶対値信号Ｌ２２（図１３参照）を得る。

次に、信号同期部３０の確率最大化部３４は、指標算出処理（ステップＳ２３）を行う。指標算出処理では、基準モデルに対する判定信号の一致の度合いを示す指標を判定信号をシフト範囲Ｈ１の範囲内で変化させながら算出する。そして、信号同期部３０の確率最大化部３４は、シフト量決定処理（ステップＳ２４）を行う。シフト量決定処理では、算出された指標が最大となる判定信号のタイミングの変化量（シフト量）を判定信号のシフト量として決定する。そして、確率最大化処理が終了される。

図１４及び図１５を参照し、指標算出処理（指標算出工程）とシフト量決定処理（同期値決定工程）とについて詳述する。まず、判定信号は、そこに関連付けられている推定同期信号Ｔｉｐが基準モデルに関連付けられている推定同期信号Ｔｉｐに一致される。

図１４に示すように、指標算出処理（ステップＳ２３）が実行されると、確率最大化部３４の指標算出部３５は、シフト量に初期値を設定する（ステップＳ３０）。シフト量の初期値は、判定信号の推定同期信号Ｔｉｐを中心にしたとき、判定信号をシフト（移動）させるシフト範囲Ｈ１「−α°」から「＋α°」までの間で特定の値が設定される。例えば、シフト量の初期値として、シフト範囲Ｈ１の最小値「−α°」が選択される。確率最大化部３４の指標算出部３５は、設定されたシフト量で判定信号のタイミングをシフト（移動）させるシフト処理を行う（ステップＳ３１）。シフト処理では、判定信号のタイミングのうち比較される全てのタイミングが、例えば判定信号の推定同期信号Ｔｉｐのタイミングが、設定されたシフト量だけ基準モデルに対してシフト（移動）される。そして、確率最大化部３４の指標算出部３５は、基準モデルとシフト量だけシフト（移動）した判定信号とを比較し（ステップＳ３２）、このシフト量における指標値を算出す（ステップＳ３３）。

ここで、指標値について説明する。指標値は、基準モデル（図１１参照）に対応するシフトされた判定信号の各タイミングにおける各スコアが累積された値として算出される。判定信号の各タイミングにおけるスコアは、そのタイミングにおける基準モデルの平均値に近いほど大きく、基準モデルの平均値から離れるほど小さく算出される。例えば、図１０（ｂ）に示すように、確率密度関数に基づく正規分布における平均と分散における相対尤度がスコアとして用いられる。そして、判定信号に算出された、対象とする範囲の全タイミングのスコアを加算して指標値を算出する。

指標値が算出されると、確率最大化部３４の指標算出部３５は、シフト量と指標値とを記憶部２０の指標値領域２０４に保持し（ステップＳ３４）、シフト量を更新する（ステップＳ３５）。シフト量は予め定められたステップ量で更新される。ステップ量は、同期精度と処理時間から適切なステップ量を決定する。シフト量が更新されると、確率最大化部３４の指標算出部３５は、更新されたシフト量がシフト範囲Ｈ１の範囲内にあるか否かを判定する（ステップＳ３６）。つまり、シフト量がシフト範囲Ｈ１の最小値「−α°」以上、かつ、最大値「＋α°」以下であるとき、シフト範囲Ｈ１の範囲内にある旨が判定される。シフト量がシフト範囲Ｈ１の範囲内にあると判定された場合（ステップＳ３６でＹＥＳ）、確率最大化部３４の指標算出部３５は、処理をステップＳ３１に戻すとともに、ステップＳ３１以下の処理を実行する。これにより、シフト範囲Ｈ１の範囲内において判定信号のタイミングのシフト量がステップ量だけ変更される毎に指標値が算出される。

一方、シフト量がシフト範囲Ｈ１の範囲内にないと判定されると（ステップＳ３６でＮＯ）、信号同期部３０は、スコア算出処理を終了し、シフト量決定処理を開始させる。
図１５のグラフＬ１６は、判定信号のタイミングをシフト範囲Ｈ１の範囲内でステップ量毎に移動し算出された指標値の一例を示す。例えば、判定信号がシフト範囲Ｈ１を移動されたとき、シフトされた量に応じて算出される指標値が異なる。グラフＬ１６では、シフト範囲Ｈ１の端である「−α°」や最大値「＋α°」に近いタイミングにシフトされたときの指標値よりも、シフト範囲Ｈ１の中央に近いタイミングにシフトされたときの指標値の方が相対的に大きく得られる例である。そして、グラフＬ１６では、シフト量が「β°」のとき、指標の最大値「Ｓｍａｘ」となっている。

スコア算出処理が終了することにより、シフト量決定処理が開始される。
図１６に示すように、シフト量決定処理（ステップＳ２４）が実行されると、確率最大化部３４のシフト量決定部３６は、指標値の最大値を選択する（ステップＳ４０）。指標値の最大値は、記憶部２０の指標値領域２０４から最大のものが選択される。例えば、図１５のグラフＬ１６に示す指標値が指標値領域２０４に保持されているとき、指標の最大値「Ｓｍａｘ」が選択される。そして、確率最大化部３４のシフト量決定部３６は、選択された指標の最大値に対応するシフト量を選択する（ステップＳ４１）。すなわち、選択された指標に対応付けられているシフト量が取得されている。例えば、図１５のグラフＬ１６に示す指標の最大値「Ｓｍａｘ」に対応付けられているシフト量「β°」が選択される。

シフト量の決定の一例は、以下の式（４）で表される。

但し、Ｐは同時確率、ｐは正規分布の確率密度関数、ｘは判定信号、μは基準モデルの平均、σは基準モデルの分散、ｔは時刻、ｓはシフト量、Ｓｈｉｆｔは選択されるシフト量。

上記式（４）は、基準モデルの範囲が判定信号から切り出された範囲以上であるときに適用されるが、基準モデルの範囲が判定信号から切り出された範囲未満であるときには下記式（５）が適用される。

シフト量が選択されると、確率最大化部３４のシフト量決定部３６は、選択したシフト量が適正か否かを判定する（ステップＳ４２）。
選択したシフト量が適正であると判定された場合（ステップＳ４２でＹＥＳ）、確率最大化部３４のシフト量決定部３６は、選択したシフト量を判定信号に対するシフト量として決定する（ステップＳ４３）。これにより、選択したシフト量が、判定信号に対するシフト量として記憶部２０のシフト量領域２０５に保持される。そして、シフト量決定処理が終了される。

一方、選択したシフト量が適正ではないと判定された場合（ステップＳ４２でＮＯ）、確率最大化部３４のシフト量決定部３６は、選択したシフト量を判定信号に対するシフト量として決定しない（ステップＳ４４）とともに、シフト量決定処理を終了する。

本実施形態によれば、筒内圧が取得されないときであれ、取得された音圧信号とＩＰ信号とに基づいて複数の取得信号を同期させるシフト量が算出される。そして、ノッキング判定装置は、複数の取得信号をそれぞれ算出されたシフト量でシフトさせることで、複数の取得信号のタイミングを同期させ、それら取得信号に基づきノッキングの有無を判定することができる。

以上説明したように、本実施形態によれば、以下に記載の効果を奏することができる。
（１）エンジン１に周期的に生じる状態をサイクルごとに検出して取得信号とし、この取得信号に基づいてエンジン１の状態を判定する。このとき、エンジン１の状態を判定するために得た取得信号からなる判定信号を、対比の基準として得た取得信号からなる基準信号に比較するとき、判定信号と基準信号との各サイクルのタイミングを一致（同期）させる必要がある。そこで、本実施形態では、判定信号の基準信号に対するタイミングのずれが基準モデルに対する一致の度合いの高いときの相対移動からシフト値として算出した。そして、この算出したシフト値で判定信号のタイミングを移動させることで、比較する基準モデルとの関係においてサイクルが同期され、取得されたサイクル中のタイミングも一致するようになる。これにより、検出した信号に対して判定すべきタイミングをより正確に推定することができる。

また、信号に含まれる大きな変化などの特徴値に着目するのではなく、平均と分散とに基づいて算出される基準モデルに対する一致の度合いを指標として算出する。このため基準モデルに含まれるノイズなどの「はずれ値」などの影響を相対的に低くすることができる。これによっても、検出した信号に対して判定すべきタイミングをより正確に推定することができる。

（２）分散の小さい部分、つまり変動が小さく値が安定しているタイミングを重視して指標が算出されることから検出した信号について同期値の算出精度の向上が図られる。
（３）取得信号に推定同期信号Ｔｉｐが含まれていることから、推定同期信号Ｔｉｐを目安に判定信号と基準モデルとを比較することが可能になるが、推定同期信号Ｔｉｐの精度が高くない場合もある。そこで、推定同期信号Ｔｉｐを含むシフト範囲Ｈ１に相対移動量を設定することで、同期信号の精度が高くない場合でも、基準モデルに対して判定信号を判定すべきタイミングをより正確に推定することができる。

（４）指標を、基準モデルの平均値のときに最大値、平均値から離れると小さくなるように算出することから、基準モデルが、平均値とそこからの離間、例えば分散に応じて定められる。これにより、一致の度合いによる指標の算出が好適に行える。

（その他の実施形態）
なお、上記実施形態は、これを適宜変更した以下の形態にて実施することもできる。
・上記実施形態では、判定信号に対するノッキングの有無の判定が、信号選択部３９、スペクトル算出部４０、バイスペクトル算出部４１、及び、ノッキング判定部５０による処理を通じて行われる場合について例示した。しかしこれに限らず、判定信号に対してノッキングの有無を判定することができるのであれば、その他の判定方法によってノッキングの有無が判定されてもよい。例えば、最も簡単に、判定信号の特定のタイミングに含まれる信号の大きさやスペクトルの状態からノッキングの有無を判定してもよい。また、上述した信号選択部３９、スペクトル算出部４０、バイスペクトル算出部４１、及び、ノッキング判定部５０による処理のうち、バイスペクトル算出部４１による処理を割愛したり、ノッキング判定部５０における処理を変更可能な範囲で変更し、ノッキングを判定するようにしてもよい。

・上記実施形態では、シフト量決定処理は、選択したシフト量が適正であるか否かを判定する処理（ステップＳ４２〜Ｓ４４）を含む場合について例示した。しかしこれに限らず、シフト量決定処理は選択したシフト量が適正であるか否かを判定する処理を含まず、シフト量を選択する処理（ステップＳ４１）で終了してもよい。

・上記実施形態では、エンジン１に発生する圧力変動に基づく物理量を空気の振動である音として音圧センサ４によって音圧信号を取得する場合について例示した。しかしこれに限らず、内燃機関に発生する圧力変動に基づく物理量は、空気の信号である音に限られるものではなく、ノッキングなどの判定対象となる現象に相関のある物理量であれば、内燃機関の加速度（振動）でもよい。

例えば、図１７に示すように、エンジン１に加速度センサ９を設け、この加速度センサ９からの加速度に関する信号をデータ収集装置５に入力して取得信号を生成させてもよい。つまり、エンジン１から得られる様々な物理量をノッキング判定処理に用いることができる。

・上記実施形態では、比較部切出し（ステップＳ１１）の処理により、基準モデルを基準信号のサイクルの一部である範囲Ｄ１に対して算出する場合について例示した。しかしこれに限らず、基準モデルを基準信号のサイクルの全部や範囲Ｄ１と相違する範囲に対して算出してもよい。これにより、サイクルの全ての部分や範囲Ｄ１と相違する範囲に対して判定信号をシフトさせて同期値を算出することができるようになる。

・上記実施形態では、比較部切出し（ステップＳ２１）の処理により、判定信号のサイクルからその一部の範囲Ｄ１を切り出す場合について例示した。しかしこれに限らず、判定信号のサイクルの全部や範囲Ｄ１と相違する範囲を基準モデルとの比較に用いてもよい。これにより、判定信号のサイクルのどの部分であっても基準モデルに比較させることができるようになる。

・上記実施形態では、確率最大化処理において、判定信号に対応付けられた推定同期信号Ｔｉｐと基準モデルに対応付けられた推定同期信号Ｔｉｐとを一致させた上でシフト範囲Ｈ１を設定する場合について例示した。しかしこれに限らず、基準モデルと判定信号とがどこかのタイミングで同期することが明らかであれば、推定同期信号なし、又は使用しなくても基準モデルに対する対象信号のシフト範囲を基準モデルに重なりが生じる全範囲とすることで、指標の最大値が算出される。よって判定信号に対するシフト量を取得することができる。

・上記実施形態では、シフト範囲Ｈ１の「０」位置（基準位置）に対して等しい距離にシフト範囲Ｈ１を区画する最大値「＋α°」と最小値「−α°」が設定される場合について例示した。しかし、シフト範囲を区画する最大値や最小値は、判定信号のどのタイミング（位置、基準位置）を基準にするかによっても影響を受けることから、最大値と最小値との距離は異なっていてもよい。これにより、指標の最大値が得られる可能性の高いタイミングのある方向へシフト範囲を設定することができる。

・上記実施形態では、シフト範囲Ｈ１を定める場合について例示した。しかしこれに限らず、シフト範囲を定めなくてもよい。シフト範囲を定めず、判定信号を基準モデルの全範囲にシフトさせることにより、指標の最大値が算出される。これにより、判定信号に対するシフト量を取得できる。

・上記実施形態では、エンジン１の回転速度は一定であることから各サイクルの長さが同じであり、取得されたサイクル中のタイミングのうちの一箇所が一致（同期）すれば各取得信号（基準信号及び判定信号）が同期する場合について例示した。しかしこれに限らず、各取得信号のサイクルの長さが異なっていても、不足するデータを補完したり、過剰なデータを間引く、公知のデータ処理技術を用い、各取得信号の長さを同じにしてもよい。

・上記実施形態では、機器状態判定装置がノッキング判定装置である場合について例示した。つまり、判定信号と基準モデル（対象信号）とを同期させる処理をエンジン１のノッキング判定に先立ち行う場合について例示した。しかしこれに限らず、機器状態判定装置がノッキング判定装置以外の異常を判定する装置であってもよい。

つまり、図１７に示すように、機器状態判定装置１１が、エンジンの運転中に生じる異常を含む各種状態のいずれを判定する装置であるとしても、判定信号と基準モデルとを同期させる必要があれば、取得した取得信号を同期させるようにしてもよい。これにより、機器状態判定装置は、エンジンから得られる周期的な取得信号に対して判定すべきタイミングがより正確に推定されて、取得信号に基づくエンジンの運転状態の判定がより適切に行える。

また、機器状態判定装置は、エンジン以外の機器の状態を測定する装置であってもよい。すなわち、機器状態判定装置は、エンジンではなくとも、周期性のある動作を伴う機器について、その動作中に生じる異常を含む各種状態の判定するとき、取得した各取得信号を同期させてもよい。これにより、機器状態判定装置は、周期性のある動作を伴う機器の周期性のある動作に基づく物理量を示す取得信号に対して、判定すべきタイミングがより正確に推定されて、取得信号に基づく機器の運転状態の判定がより適切に行える。

１…エンジン、２…エンジンＥＣＵ、３…電流センサ、４…音圧センサ、５…データ収集装置、６…モニタ、７…点火プラグ、８…点火装置、９…加速度センサ、１０…ノッキング判定装置、１１…機器状態判定装置、２０…記憶部、３０…信号同期部、３１…モデル作成部、３４…確率最大化部、３５…指標算出部、３６…シフト量決定部、３９…信号選択部、４０…スペクトル算出部、４１…バイスペクトル算出部、５０…ノッキング判定部、５１…統計パラメータ算出部、５２…正規化部、５３…判定値算出部、５４…比較部、６０…管理部、６１…基準信号取得部、６２…判定信号取得部、６３…データ整理部、１００…車両、２０１…基準信号領域、２０２…基準モデル領域、２０３…判定信号領域、２０４…指標値領域、２０５…シフト量領域。

Claims (14)

1. 周期性のある動作を伴う機器の前記周期性のある動作に基づく物理量を示す信号としてサイクル毎に取得された各取得信号の特定のタイミングの信号状態に基づいて前記機器の状態を判定する機器状態判定装置であって、
   前記機器の基準となる条件での動作に基づく前記取得信号を前記サイクル毎に基準信号として取得する基準信号取得部と、
   前記機器の動作状態を判定するサイクルでの前記取得信号を判定信号として取得する判定信号取得部と、
   前記各基準信号の取得されたサイクル中の前記タイミングでの前記基準信号の集合に対する平均と分散とに基づいて基準モデルを作成するモデル作成部と、
   前記基準モデルの作成されたサイクル中の前記タイミングに対して、前記判定信号の取得されたサイクル中の前記タイミングを順次相対移動させると共に、前記相対移動毎の前記判定信号の前記基準モデルに対する一致の度合いを示す指標を算出する指標算出部と、
   前記一致の度合いが最大であることを示す指標に対応する相対移動に基づいて前記基準モデルに前記判定信号が同期する同期値を決定する同期値決定部と、を備え、
   前記同期値で移動させた前記タイミングでの前記判定信号に基づいて前記判定を行う
   機器状態判定装置。
2. 前記指標算出部は、前記指標を前記基準モデルの分散が大きいタイミングにおける一致の度合いよりも分散が小さいタイミングにおける一致の度合いを重視して算出する
   請求項１に記載の機器状態判定装置。
3. 前記指標算出部は、前記判定信号の取得されたサイクル中のタイミングを前記判定信号が取得されたときに定まるタイミングを含む所定の範囲内で順次相対移動させる
   請求項１又は２に記載の機器状態判定装置。
4. 前記取得信号には、取得されたサイクル中のタイミングの目安となる同期信号が対応付けられており、
   前記指標算出部は、前記判定信号の取得されたサイクル中のタイミングを前記判定信号の前記同期信号のタイミングと前記基準モデルの前記同期信号のタイミングとを一致させたタイミングを含む所定の範囲で順次相対移動させる
   請求項１又は２に記載の機器状態判定装置。
5. 前記モデル作成部は、前記基準信号のうち所定のタイミングからなる範囲に対応する基準モデルを作成し、
   前記指標算出部は、前記判定信号のうち前記所定のタイミングからなる範囲について前記基準モデルに対する前記指標を算出する
   請求項１〜４のいずれか一項に記載の機器状態判定装置。
6. 前記指標算出部は、前記指標を前記基準モデルの平均値のとき最大値となり、前記平均値から離れることに応じて小さくなるように算出し、
   前記同期値決定部は、前記指標の最大値を前記一致の度合いが最大であるものとする
   請求項１〜５のいずれか一項に記載の機器状態判定装置。
7. 前記機器は内燃機関であり、
   前記サイクルは前記内燃機関の燃焼サイクルであり、
   前記同期信号は前記内燃機関の点火プラグへ電力が供給されることで得られる信号である
   請求項４に記載の機器状態判定装置。
8. 周期性のある動作を伴う機器の前記周期性のある動作に基づく物理量を示す信号としてサイクル毎に取得した各取得信号の特定のタイミングの信号状態に基づいて前記機器の状態を判定する機器状態判定装置で実行される機器状態判定方法であって、
   機器状態判定装置は、
   前記機器の基準となる条件での動作に基づく前記取得信号を前記サイクル毎の基準信号として取得する基準信号取得工程と、
   前記機器の動作状態を判定するサイクルでの前記取得信号を判定信号として取得する判定信号取得工程と、
   前記各基準信号の取得されたサイクル中の前記タイミングでの前記基準信号の集合に対する平均と分散とに基づいて基準モデルを作成するモデル作成工程と、
   前記基準モデルの作成されたサイクル中の前記タイミングに対して、前記判定信号の取得されたサイクル中の前記タイミングを順次相対移動させると共に、前記相対移動毎の前記判定信号の前記基準モデルに対する一致の度合いを示す指標を算出する指標算出工程と、
   前記一致の度合いが最大であることを示す指標に対応する相対移動に基づいて前記基準モデルに前記判定信号が同期する同期値を決定する同期値決定工程と、を備え、
   前記同期値で移動させた前記タイミングでの判定信号に基づいて前記判定を行う
   機器状態判定方法。
9. 内燃機関の燃焼サイクルに基づく物理量を示す信号としてサイクル毎に取得された各取得信号の特定のタイミングの信号状態に基づいて前記内燃機関のノッキングの有無を判定するノッキング判定装置であって、
   前記内燃機関の基準となる条件での動作に基づく前記取得信号を前記サイクル毎に基準信号として取得する基準信号取得部と、
   前記内燃機関の動作状態を判定するサイクルでの前記取得信号を判定信号として取得する判定信号取得部と、
   前記各基準信号の取得されたサイクル中の前記タイミングでの前記基準信号の集合に対する平均と分散とに基づいて基準モデルを作成するモデル作成部と、
   前記基準モデルの作成されたサイクル中の前記タイミングに対して、前記判定信号の取得されたサイクル中の前記タイミングを順次相対移動させると共に、前記相対移動毎の前記判定信号の前記基準モデルに対する一致の度合いを示す指標を算出する指標算出部と、
   前記一致の度合いが最大であることを示す指標に対応する相対移動に基づいて前記基準モデルに前記判定信号が同期する同期値を決定する同期値決定部と、を備え、
   前記同期値で移動させた前記タイミングでの前記判定信号に基づいて前記判定を行う
   ノッキング判定装置。
10. 前記指標算出部は、前記指標を前記基準モデルの分散が大きいタイミングにおける一致の度合いよりも分散が小さいタイミングにおける一致の度合いを重視して算出する
    請求項９に記載のノッキング判定装置。
11. 前記指標算出部は、前記判定信号の取得されたサイクル中のタイミングを前記判定信号が取得されたときに定まるタイミングを含む所定の範囲内で順次相対移動させる
    請求項９又は１０に記載のノッキング判定装置。
12. 前記取得信号には、取得されたサイクル中のタイミングの目安となるイグニッションパルスに対応する信号が同期信号として対応付けられており、
    前記指標算出部は、前記判定信号の取得されたサイクル中のタイミングを前記判定信号の前記同期信号のタイミングと前記基準モデルの前記同期信号のタイミングとを一致させたタイミングを含む所定の範囲で順次相対移動させる
    請求項９又は１０に記載のノッキング判定装置。
13. 前記モデル作成部は、前記基準信号から前記燃焼サイクル中の所定のタイミング範囲に対応する基準モデルを作成し、
    前記指標算出部は、前記判定信号のうち前記所定のタイミング範囲について前記基準モデルに対する前記指標を算出する
    請求項９〜１２のいずれか一項に記載のノッキング判定装置。
14. 前記指標算出部は、前記指標を前記基準モデルの平均値のとき最大値となり、前記平均値から離れることに応じて小さくなるように算出し、
    前記同期値決定部は、前記指標の最大値を前記一致の度合いが最大であるものとする
    請求項９〜１３のいずれか一項に記載のノッキング判定装置。