(57)【要約】【課題】本発明は、ノックの発生をより正確に検知可能なフィルタ、ウィンドウの適合を行うことを目的としている。【解決手段】このため、筒内圧センサとノックセンサと処理装置を備え、ノックセンサの適合値の候補選定を行うノック制御の適合方法において、筒内圧センサ波形とノックセンサ波形を同時測定するステップと、測定波形に対して時系列の周波数特性を算出するステップと、算出された周波数特性に基づく適合値の候補選定を行うステップと、適合値の候補の中から最良適合値を選定するステップとを有する。また、ノック制御の適合装置において、処理装置は、波形測定手段と周波数分析手段と適合値設定手段とノック検出性評価指標算出手段とを有する。更に、処理装置によって求められたフィルタとクランク角ウィンドウを用いて波形からノック検出するとともに点火時期制御する。【選択図】図１
2011179322.tif 000002
【特許請求の範囲】
【請求項１】
エンジンの気筒内部の圧力を検知する筒内圧センサと、前記エンジンの側壁の任意の位置に取り付けて振動を検知するノックセンサと、これら筒内圧センサとノックセンサとから出力させる検知信号を波形信号として測定し処理する処理装置とを備え、その処理に基づいてノックセンサの適合値の候補の選定を行うノック制御の適合方法において、筒内圧センサ波形とノックセンサ波形を同時測定するステップと、これら測定された波形それぞれに対して時系列の周波数特性を算出するステップと、この算出された周波数特性に基づく適合値の候補の選定を行うステップと、最終的に適合値の候補の中から最良な適合値を選定するステップとを有することを特徴とするノック制御の適合方法。
【請求項２】
同時測定されたノックセンサ波形における適合値の候補に対する振動強度を算出するステップと、筒内圧センサ波形に対し可聴域および広帯域のフィルタまたはその他の通過帯域を持つフィルタを用いてそれぞれの波形領域に応じた複数の圧力変動強度を算出するステップと、この複数の圧力変動強度と振動強度との相関係数を算出するステップとを行うことを特徴とする請求項１に記載のノック制御の適合方法。
【請求項３】
相関係数を用いて式４によりノック検出性評価指標を算出するステップと、このノック検出性評価指標を用いて適合値の候補の中から最良適合値を選定するステップとを有することを特徴とする請求項２に記載のノック制御の適合方法。
【請求項４】
エンジンの気筒内部の圧力を検知する筒内圧センサと、前記エンジンの側壁の所定の位置に取り付けて振動を検知するノックセンサと、これら筒内圧センサとノックセンサとから出力させる検知信号を波形信号として測定し処理し、その波形信号の処理に基づいてノックセンサの適合値の候補の選定が可能な処理装置とを備えるノック制御の適合装置において、前記処理装置は、筒内圧センサ波形とノックセンサ波形を同時に測定する波形測定手段と、これら測定された波形それぞれに対して時系列の周波数特性を算出する周波数分析手段と、この算出された周波数特性に基づく適合値の候補の選定を行う適合値設定手段と、この適合値設定手段からの出力に基づき最終的に適合条件の候補の中から最良な適合値を選定するノック検出性評価指標算出手段とを有することを特徴とするノック制御の適合装置。
【請求項５】
前記処理装置は、同時測定されたノックセンサ波形に対する適合値の候補に対応する振動強度を算出する振動強度算出手段と、筒内圧センサ波形に対し可聴域あるいは広帯域のフィルタまたはその他の通過帯域を持つフィルタを用いてそれぞれの波形領域に応じた複数の圧力変動強度を算出する圧力変動強度算出手段と、この複数の圧力変動強度と振動強度との相関係数を算出する相関係数算出手段とを有することを特徴とする請求項４に記載のノック制御の適合装置。
【請求項６】
前記ノック検出性評価指標算出手段は、ノック検出性評価指標の算出と、このノック検出性評価指標を用いて適合値の候補の中から最良適合値の選定とを行うことを特徴とする請求項５に記載のノック制御の適合装置。
【請求項７】
請求項４ないし６に記載された前記処理装置によって最良適合値として求められたフィルタとクランク角ウィンドウを用いて前記ノックセンサによって検出される波形からノック検出するとともに点火時期制御することを特徴とするノック制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明はノック制御の適合方法、適合装置、および、ノック制御装置に係り、特に、エンジンのシリンダブロック表面に設置されたノックセンサによりノックの発生を正確に検知するための技術であって、制御パラメータを適合させる技術に関するものである。
【背景技術】
【０００２】
エンジンにノックが発生すると、エンジン特有の複数の共鳴周波数を持った振動が発生する。
このノックを検出するために、エンジンのシリンダブロック表面にノックセンサを設置している。
このとき、精度の高いノック検出を行うためには、ノックセンサに関するバンドパスフィルタの通過周波数およびウィンドウの抽出時間帯が重要となっている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開平０４−８１５４６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ところで、従来の一般的なノック制御では、シリンダブロック表面に設置されたノックセンサで振動を検出し、これにノック特有の周波数のみを取り出すためのバンドパスフィルタ（「フィルタ」ともいう。）、およびノックの発生時間帯の振動波形のみを取り出すためのウインドウ処理を行い、得られた波形の最大値や積分値をノック強度の指標として、ノックの発生が認められた場合には点火時期をリタードさせることによりノックを抑制する。
このとき、様々な機械的振動が混在する振動波形からノックに起因する振動を精度良く抽出し、ノックの発生を正確に検知するには、バンドパスフィルタの通過周波数およびウィンドウの抽出時間帯の適合が重要となる。
【０００５】
従来、ノックセンサで検出されるシリンダブロック表面の振動特性のみを評価して、フィルタとウィンドウの適合が行われることが多く、その場合はノックに起因する振動あるいは機械的振動のどちらを検出したのかを特定する術が無いという不都合がある。
結果として、ノックを精度良く検知できないことによりノックを十分に抑制できず、エンジンの耐久性の低下や運転者にノック音が聞こえることによる製品品質の低下、または機械的振動を誤検知して点火時期をリタードさせてしまうことによる出力および燃費性能の低下や排気温度の過度の上昇を招いていた。
【０００６】
このような問題に対し、機械的振動の影響をほとんど受けないとされる筒内圧センサを用いて、ノックによるエンジンの燃焼室内の圧力変動の大きさを測定し、これとノックセンサで検出されるシリンダブロック表面の振動の大きさとの相関係数を求めて対処する方法が上記の特許文献１に開示されている。
しかし、この特許文献１に示されているものは、以下の点で不十分であった。
（１）特許文献１では、筒内圧センサで検出される圧力変動の大きさ（以下、「圧力変動強度」という。）とノックセンサで検出される振動の大きさ（以下、「振動強度」という。）の相関関係から、ノック強度を求める上でノック周波数を特定するとの記述があるのみで、ウィンドウの最適化については全く開示されていない。
（２）圧力変動強度を求める際に、筒内圧センサでは検出が不完全なノック共鳴モードの存在が考慮されていない。
（３）「耳に聞こえるノック」（「可聴域ノック」という。）と「筒内圧センサの波形そのものが大きいノック」（「広帯域ノック」という。）とのように複数のノックを定義して、運転条件等の諸条件により各々のノックの定義のうちあるノックの定義を重視してフィルタ・ウィンドウ、つまり適合値の最適化を行いたい場合があるが、上記の特許文献１には、このような場合に対処する方法が開示されていない。
【０００７】
この発明は、上述した問題点を解決し、ノックの発生をより正確に検知可能なフィルタ、ウィンドウの適合を行うことを目的としている。
【課題を解決するための手段】
【０００８】
そこで、この発明は、上述不都合を除去するために、エンジンの気筒内部の圧力を検知する筒内圧センサと、前記エンジンの側壁の任意の位置に取り付けて振動を検知するノックセンサと、これら筒内圧センサとノックセンサとから出力させる検知信号を波形信号として測定し処理する処理装置とを備え、その処理に基づいてノックセンサの適合値の候補の選定を行うノック制御の適合方法において、筒内圧センサ波形とノックセンサ波形を同時測定するステップと、これら測定された波形それぞれに対して時系列の周波数特性を算出するステップと、この算出された周波数特性に基づく適合値の候補の選定を行うステップと、最終的に適合値の候補の中から最良な適合値を選定するステップとを有することを特徴とする。
また、エンジンの気筒内部の圧力を検知する筒内圧センサと、前記エンジンの側壁の所定の位置に取り付けて振動を検知するノックセンサと、これら筒内圧センサとノックセンサとから出力させる検知信号を波形信号として測定し処理し、その波形信号の処理に基づいて複数のノックセンサの適合値の候補の選定が可能な処理装置とを備えるノック制御の適合装置において、前記処理装置は、筒内圧センサ波形とノックセンサ波形を同時に測定する波形測定手段と、これら測定された波形それぞれに対して時系列の周波数特性を算出する周波数分析手段と、この算出された周波数特性に基づく適合値の候補の選定を行う適合値設定手段と、この適合値設定手段からの出力に基づき最終的に適合条件の候補の中から最良な適合値を選定するノック検出性評価指標算出手段とを有することを特徴とする。
更に、請求項４ないし６に記載された前記処理装置によって最良適合値として求められたフィルタとクランク角ウィンドウを用いて前記ノックセンサによって検出される波形からノック検出するとともに点火時期制御することを特徴とする。
【発明の効果】
【０００９】
以上詳細に説明した如くこの発明によれば、エンジンの気筒内部の圧力を検知する筒内圧センサと、エンジンの側壁の任意の位置に取り付けて振動を検知するノックセンサと、これら筒内圧センサとノックセンサとから出力させる検知信号を波形信号として測定し処理する処理装置とを備え、その処理に基づいてノックセンサの適合値の候補の選定を行うノック制御の適合方法において、筒内圧センサ波形とノックセンサ波形を同時測定するステップと、これら測定された波形それぞれに対して時系列の周波数特性を算出するステップと、この算出された周波数特性に基づく適合値の候補の選定を行うステップと、最終的に適合値の候補の中から最良な適合値を選定するステップとを有する。
また、エンジンの気筒内部の圧力を検知する筒内圧センサと、エンジンの側壁の所定の位置に取り付けて振動を検知するノックセンサと、これら筒内圧センサとノックセンサとから出力させる検知信号を波形信号として測定し処理し、その波形信号の処理に基づいてノックセンサの適合値の候補の選定が可能な処理装置とを備えるノック制御の適合装置において、前記処理装置は、筒内圧センサ波形とノックセンサ波形を同時に測定する波形測定手段と、これら測定された波形それぞれに対して時系列の周波数特性を算出する周波数分析手段と、この算出された周波数特性に基づく適合値の候補の選定を行う適合値設定手段と、この適合値設定手段からの出力に基づき最終的に適合条件の候補の中から最良な適合値を選定するノック検出性評価指標算出手段とを有する。
従って、ノック周波数とノック発生時間の全体の傾向を把握し、詳細にノック検出性の評価を行う適合値の候補を絞り込むことで適合工数の削減が可能となる。
更に、請求項４ないし６に記載された前記処理装置によって最良適合値として求められたフィルタとクランク角ウィンドウを用いて前記ノックセンサによって検出される波形からノック検出するとともに点火時期制御する。
従って、最適なフィルタおよびウィンドウ条件を用いることにより高いノック制御性を得ることができる。
【図面の簡単な説明】
【００１０】
【図１】図１はこの発明の第１実施例を示すノック制御の適合方法の処理内容のフローチャートである。（実施例１）
【図２】図２はノック制御の適合装置の構成図である。（実施例１）
【図３】図３はエンジンに取り付けられたセンサと検出される信号を示す図である。（実施例１）
【図４】図４は適合手段の構成を示す図である。（実施例１）
【図５】図５はノック発生波形とノック非発生波形とを示し、（ａ）はノック発生波形を示す図、（ｂ）はノック非発生波形を示す図である。（実施例１）
【図６】図６は周波数分析手段による時系列の周波数分析のウィンドウの設定を示す図である。（実施例１）
【図７】図７は周波数分析手段による一波形に対するウィンドウ別の周波数特性を示す図である。（実施例１）
【図８】図８は各ウィンドウで平均した周波数特性を示す図である。（実施例１）
【図９】図９は周波数分析手段による周波数分析用フローチャートである。（実施例１）
【図１０】図１０は適合値設定手段におけるノック共鳴モードと共鳴周波数の例を示す図である。（実施例１）
【図１１】図１１は適合値設定手段におけるノック周波数の判断用フローチャートである。（実施例１）
【図１２】図１２は適合値設定手段における適合値の候補選定の判断用フローチャートである。（実施例１）
【図１３】図１３は適合値設定手段におけるノックセンサに関する適合値の候補の選定例を示す図である。（実施例１）
【図１４】図１４は振動強度算出手段における適合値の候補に対応する振動強度算出の例を示すタイムチャートである。（実施例１）
【図１５】図１５は振動強度算出手段における振動強度算出処理用のフローチャートである。（実施例１）
【図１６】図１６は圧力変動強度算出手段における可聴域と広帯域のフィルタの特性を示し、（ａ）は可聴域フィルタを示す図、（ｂ）は広帯域フィルタを示す図である。（実施例１）
【図１７】図１７は圧力変動強度算出手段における圧力変動強度算出の例を示すタイムチャートである。（実施例１）
【図１８】図１８は筒内圧センサ波形のノック周波数特性と圧力変動強度を示す図である。（実施例１）
【図１９】図１９は圧力変動強度算出手段における圧力変動強度算出処理用のフローチャートである。（実施例１）
【図２０】図２０は相関係数算出手段における圧力変動強度と振動強度との相関図を示し、（ａ）は可聴域ノックに関する相関図、（ｂ）は広帯域ノックに関する相関図である。（実施例１）
【図２１】図２１は相関係数算出手段における処理用のフローチャートである。（実施例１）
【図２２】図２２はノック検出性評価指標算出手段における処理を示す図である。（実施例１）
【図２３】図２３はノック検出性評価指標算出手段における処理用のフローチャートである。（実施例１）
【図２４】図２４は筒内圧センサ波形とノックセンサ波形の周波数特性（１６００ｒｐｍ、３番気筒）を示す図である。（実施例１）
【図２５】図２５は筒内圧センサ波形のノック周波数特性と圧力変動強度Ｓｎ（１６００ｒｐｍ、３番気筒）を示す図である。（実施例１）
【図２６】図２６はノックセンサの適合値の候補に関する相関図を示し、（ａ）は可聴域ノックにおける適合値の候補（１）に対する振動強度Ｉｋ（３，１）と圧力変動強度Ｉｃｈ（３）に関する相関図、（ｂ）は広帯域ノックにおける適合値の候補（１）に対する振動強度Ｉｋ（３，１）と圧力変動強度Ｉｃｗ（３）に関する相関図、（ｃ）は可聴域ノックにおける適合値の候補（２）に対する振動強度Ｉｋ（３，２）と圧力変動強度Ｉｃｈ（３）に関する相関図、（ｄ）は広帯域ノックにおける適合値の候補（２）に対する振動強度Ｉｋ（３，２）と圧力変動強度Ｉｃｗ（３）に関する相関図である。（実施例１）
【図２７】図２７はこの発明の第２実施例を示すノック制御の適合装置の構成図である。（実施例２）
【発明を実施するための形態】
【００１１】
以下図面に基づいてこの発明の実施例を詳細に説明する。
【実施例１】
【００１２】
図１〜図２６はこの発明の第１実施例を示すものである。
図２において、１はノック制御の適合装置である。
この適合装置１は、エンジン２の気筒内部の圧力を検知する筒内圧センサ３と、前記エンジン２の側壁の所定の位置に取り付けて振動を検知するノックセンサ４と、これら筒内圧センサ３とノックセンサ４とから出力させる検知信号を波形信号として測定し処理し、その波形信号の処理に基づいてノックセンサ４の適合値（「適合条件」とも換言できる。）の候補の選定が可能な処理装置５とを備えている。
このとき、前記筒内圧センサ３は、この説明において多気筒の４気筒＃１、＃２、＃３、＃４を想定しているため、各気筒＃１、＃２、＃３、＃４の圧力を検知する第１〜第４筒内圧センサ３ａ、３ｂ、３ｃ、３ｄからなる。
また、前記適合装置１は、図２に示す如く、前記エンジン２に連絡される一方、マッチングツールとしてのＥＣＵデータ書込装置６が接続されるノック制御装置（「ＥＣＵ」ともいう。）７を備えている。
前記処理装置５は、図２に示す如く、筒内圧センサ波形とノックセンサ波形を同時に測定する波形測定手段８と、記録装置９と、適合手段１０とを有している。
このとき、波形測定手段８は、図示しないパーソナルコンピュータ（「ＰＣ」ともいう。）に接続可能なデータレコーダ（図示せず）を使用する。
更に、前記記録装置９は、ＰＣの図示しないハードディスク（「ＨＤＤ」ともいう。）やデータレコーダのメモリ（図示せず）などからなる。
前記エンジン２に各気筒の燃焼室（図示せず）内の圧力変動を検出する前記第１〜第４筒内圧センサ３ａ、３ｂ、３ｃ、３ｄを取り付けられている。
前記エンジン２のシリンダブロック２ａには、このシリンダブロック２ａ表面の振動を検出する前記ノックセンサ４が取り付けられている。
また、前記エンジン２にはクランクシャフトおよびカムシャフトの回転角を夫々検出するクランク角センサ１１およびカム角センサ１２が取り付けられている。
つまり、図２に示す如く、前記筒内圧センサ３は各気筒＃１、＃２、＃３、＃４の燃焼室内の圧力変動を検出し、前記ノックセンサ４はシリンダブロック２ａ表面の振動を検出し、前記クランク角センサ１１はクランクシャフト（図示せず）の回転角を検出するとともに、前記カム角センサ１２はカムシャフト（図示せず）の回転角を検出する。
そして、前記エンジン２を、適合を行いたい所定の条件で（例えば回転数、負荷を調整して）運転し、これらの４種のセンサで検出した波形を前記波形測定手段８で同時測定し、前記記録装置９に記録し保存するものである。
また、前記適合手段１０は、記録波形に複数の処理を行い、ノック制御の適合値の最適化を行う。
この適合手段１０では、ノックセンサ４の適合値の候補に対し最適化を行うため、同時測定した波形の測定条件（エンジン回転数と気筒）を測定波形から判定する必要がある。エンジン回転数は前記クランク角センサ１１により判定可能であり、気筒はクランク角センサ１１と前記カム角センサ１２とにより判定可能である。
各測定条件に対し、ノック発生状態とノック非発生状態とのそれぞれについて、後述する相関係数が問題なく算出できる測定波形の数を設定して、ノックが発生するタイミングで各センサの波形を同時測定する。
ノック発生状態とノック非発生状態との切り替えは、前記ＥＣＵデータ書込装置６を用いて前記ノック制御装置７のデータを変更し、点火時期や空燃比を変化させて行う。
追記すれば、「ノック発生状態」とは、図５において（ａ）ノック発生波形と（ｂ）ノック非発生波形とが混在する状態である。つまり、ノックが発生するタイミングの複数サイクルのうち、ある確率でノックが発生する状態のことである。このとき、発生するノック強度は一定ではない。
また、「ノック非発生状態」とは、ノックが発生するタイミングの複数サイクルでノックが発生せず、図５（ｂ）のノック非発生波形が継続して出現する状態のことである。
【００１３】
そして、前記記録装置９は、図４に示す如く、ノック発生時及びノック非発生時における前記筒内圧センサ３からの検知信号を記録する第１記録部１３と、ノック発生時及びノック非発生時における前記ノックセンサ４からの検知信号を記録する第２記録部１４とを備えている。
また、前記処理装置５は、筒内圧センサ波形とノックセンサ波形を同時に測定する前記波形測定手段８と、筒内圧センサ波形とノックセンサ波形とを記録可能な前記記録装置９と、これら測定された波形それぞれに対して時系列の周波数特性を算出する周波数分析手段１５と、この算出された周波数特性に基づく適合値（「適合条件」とも換言できる。）の候補の選定を行う適合値設定手段１６と、この適合値設定手段１６からの出力に基づき最終的に適合条件の候補の中から最良な適合値を選定するノック検出性評価指標算出手段１７とを有する。
従って、ノック周波数とノック発生時間の全体の傾向を把握し、詳細にノック検出性の評価を行う適合値の候補を絞り込むことで適合工数の削減が可能となる。
詳述すれば、図４から明らかなように、前記適合手段１０は、６つの前記周波数分析手段１５と、前記適合値設定手段１６と、振動強度算出手段１８と、圧力変動強度算出手段１９と、相関係数算出手段２０と、前記ノック検出性評価指標算出手段１７とを備えている。
そして、前記周波数分析手段１５は、前記記録装置９の第１記録部１３からノック発生時及びノック非発生時における前記筒内圧センサ３からの検知信号と、第２記録部１４からノック発生時及びノック非発生時における前記ノックセンサ４からの検知信号とを入力する。
前記適合値設定手段１６は、周波数分析手段１５に接続されている。
前記振動強度算出手段１８は、前記記録装置９と適合値設定手段１６とに接続されている。
前記圧力変動強度算出手段１９は、前記記録装置９に接続されている。
前記相関係数算出手段２０は、前記振動強度算出手段１８と圧力変動強度算出手段１９とに接続されている。
前記ノック検出性評価指標算出手段１７は、前記相関係数算出手段２０に接続されている。
つまり、前記記録装置９内の筒内圧センサ波形と前記ノックセンサ４で検出されたノックセンサ波形を用いて、図４に示す如く、ノック制御の最適化を行う前記適合手段１０を構成する。
【００１４】
追記すれば、前記適合装置１は、図２に示す如く、エンジン２に筒内圧センサ３とノックセンサ４とクランク角センサ１１とカム角センサ１２との４種のセンサが取り付けられており、各センサにより図３に示す信号が検出される。
前記波形測定手段８により、４種のセンサで検出した波形を同時に測定する。
前記記録装置９は、同時測定した各センサの信号を記録する。
前記適合手段１０は、記録波形に複数の処理を行い、適合値の最適化を行う。
この適合手段１０は、各測定条件について適合値の最適化を行うため、同時測定した波形の測定条件（エンジン回転数、気筒）を測定波形から判定する必要がある。
エンジン回転数は前記クランク角センサ１１により判定することが可能であり、気筒は前記クランク角センサ１１とカム角センサ１２とにより判定することが可能である。
そして、各測定条件に対し、「ノック発生状態」と「ノック非発生状態」とのそれぞれについて、相関係数が問題なく算出できる測定波形の数を設定して、ノックが発生するタイミングで各センサの波形を同時測定する。
「ノック発生状態」  ：図５（ａ）のノック発生波形と図５（ｂ）のノック非発生波形とが混在する状態。
つまり、ノックが発生するタイミングの複数サイクルのうち、ある確率でノックが発生する状態のこと。
また、発生するノックの強度は一定でない。
「ノック非発生状態」：ノックが発生するタイミングの複数サイクルでノックが発生せず、図５（ｂ）のノック非発生波形が継続して出現する状態のこと。
なお、「ノック発生状態」と「ノック非発生状態」との切換は、点火時期や空燃比を変化させて行う。
前記記録装置９内の筒内圧センサ波形とノックセンサ波形とを用いて適合値の最適化を行う前記適合手段１０の構成を図４に示す。
この適合手段１０は、図４に示す如く、前記周波数分析手段１５、前記適合値設定手段１６、前記振動強度算出手段１８、前記圧力変動強度算出手段１９、前記相関係数算出手段２０、前記ノック検出性評価指標算出手段１７の６つの手段で構成される。
【００１５】
前記処理装置５は、同時測定されたノックセンサ波形に対する適合値の候補に対応する振動強度を算出する前記振動強度算出手段１８と、筒内圧センサ波形に対し可聴域あるいは広帯域のフィルタまたはその他の通過帯域を持つフィルタを用いてそれぞれの波形領域に応じた複数の圧力変動強度を算出する前記圧力変動強度算出手段１９と、この複数の圧力変動強度と振動強度との相関係数を算出する前記相関係数算出手段２０とを有する。
従って、圧力変動強度算出時に前記筒内圧センサ３では検出が不完全なノック共鳴モードを考慮することで、相関係数算出の精度が向上する。また、可聴域と広帯域ノック、その他のノックを定義することで、可聴ノック音による品質の問題や広帯域ノックによる耐久性の問題、その他のノックによる問題に、それぞれ個別に対処することが可能である。
【００１６】
前記処理装置５によって最良適合値として求められたフィルタとクランク角ウィンドウを用いて前記ノックセンサ４によって検出される波形からノック検出するとともに点火時期制御する前記ノック制御装置７とする。
従って、最適なフィルタおよびウィンドウ条件を用いることにより高いノック制御性を得ることができる。
【００１７】
追記すれば、この発明は、上記の特許文献１に対して記載した問題点（１）〜（３）までを解決した適合値の最適化により、高いノック検出性を実現するものである。
適合値の最適化は、エンジン回転数別・気筒別に実施する。
ここで、説明の理解を容易とするために、以下の（Ａ）及び（Ｂ）の２つの前提を設ける。
（Ａ）エンジン回転数・気筒の組合せ（以下、「測定条件」という。）の数をｋ個とし、測定条件の識別番号をｉ（ｉ＝１〜ｋ）とする。
例えば、４気筒エンジンで回転数を５段階に分けて適合値の最適化を行う場合、測定条件の数は２０個（ｋ＝２０）となる。
（Ｂ）適合値の候補数をｎとし、適合値の候補の識別番号をｒ（ｒ＝１〜ｎ）とする。
また、図１に示すノック制御の最適化方法を処理フローチャート順に説明する。
（１）図１の処理（１０１）において、
Ｋ個の各測定条件について、前記筒内圧センサ３による燃焼室内の圧力変動とノックセンサ５によるシリンダブロック２ａ表面の振動を同時に測定する。
これにより、後述する相関係数の算出が可能となる。
（２）図１の処理（１０２）において、
Ｋ個の各測定条件について、筒内圧センサ波形とノックセンサ波形について時系列の周波数分析を行う。      図１の処理（１０３）において、
ノック周波数とノック発生期間の傾向から各測定条件において複数の適合値を選定する。
選定する適合値の候補は、各測定条件で異なっても良いが、以下の本文では簡略な説明とするために、各測定条件に対して共通のｎ個の適合値の候補を選定したとして説明を進める。
（３）図１の処理（１０４）において、
ｉ番目の測定条件で検出されたノックセンサ波形とｒ番目の適合値の候補から、振動強度Ｉｋ（ｉ、ｒ）を算出する（ｉ＝１〜ｋ、ｒ＝１〜ｎ）。
振動強度Ｉｋ（ｉ、ｒ）は、振動強度の指標を表す複数の値の集合である。
（４）図１の処理（１０５）において、
耳に聞こえるノックと定義される可聴域ノックと圧力変動そのものが大きいノックと定義される広帯域ノックとを定義して、ｉ番目の測定条件の筒内圧センサ波形から２つのノックの定義に対応する圧力変動強度Ｉｃｈ（ｉ）、Ｉｃｗ（ｉ）を算出する（ｉ＝１〜ｋ）。
圧力変動強度Ｉｃｈ（ｉ）、Ｉｃｗ（ｉ）は、それぞれ圧力変動強度を表す複数の値の集合である。
前記筒内圧センサ３で検出が不完全なノック共鳴モードが存在する場合には、その影響を考慮する。
また、可聴域、広帯域ノック以外にノックの定義が複数あっても構わない。
（５）図１の処理（１０６）において、
振動強度Ｉｋ（ｉ、ｒ）と可聴域ノックに対応する圧力変動強度Ｉｃｈ（ｉ）から相関係数Ｃｈ（ｉ、ｒ）を、振動強度Ｉｋ（ｉ、ｒ）と広帯域ノックに対応する圧力変動強度Ｉｃｗ（ｉ）から相関係数Ｃｗ（ｉ、ｒ）をそれぞれ算出する（ｉ＝１〜ｋ、ｒ＝１〜ｎ）。
（６）図１の処理（１０７）において、
相関係数Ｃｈ（ｉ、ｒ）、Ｃｗ（ｉ、ｒ）から、ｉ番目の測定条件、ｒ番目の適合値の候補に関するノック検出性評価指標Ｉｎ（ｉ、ｒ）を算出する（ｉ＝１〜ｋ、ｒ＝１〜ｎ）。
（７）図１の処理（１０８）において、
（ｋ  ｘ  ｎ）個の全てのノック検出性評価指標Ｉｎ（ｉ、ｒ）を総合的に比較し、最良適合値を選定する。
【００１８】
ここで、前記周波数分析手段１５について追記する。
この周波数分析手段１５では、測定条件を選択した後、ノック発生状態の筒内圧センサ波形とノックセンサ波形、ノック非発生状態の筒内圧センサ波形とノックセンサ波形の４種にに対して、時系列の周波数分析を行う。
前記周波数分析手段１５による時系列の周波数分析用処理のフローチャートを図９に開示する。
この図９の時系列の周波数分析用処理のフローチャートにおいては、図６に示す如く、測定波形１つ１つに対して時系列にウィンドウ（この例では、クランク角１０度おきに（ａ）〜（ｆ）までの６つのウィンドウを設定、ウィンドウ幅はクランク角３０度）を設定し、ＦＦＴやＤＦＴを用いて図７のようにウィンドウ別に周波数特性を求める。
次に、一例として図８に示すように、ｉ番目の測定条件に関する前記筒内圧センサ３と前記ノックセンサ４の全サイクルの波形、つまり相関係数が問題なく算出できるように設定した測定波形数の全波形の周波数分析結果を各ウィンドウで平均し、それぞれの結果についてノック発生状態とノック非発生状態とで比較する形とする。
この処理を全ｋ個の測定条件で行うことで、ノック周波数とノック発生時間の全体の傾向を把握することができ、後述する前記適合値設定手段１６で適合値の候補を絞り込むことができる。
【００１９】
前記適合値設定手段１６について追記する。
この適合値設定手段１６では、前記周波数分析手段１５で得られた結果をもとに、基本的にノックセンサ波形においてノック発生状態とノック非発生状態との差が明確に現れている周波数とクランク角範囲から適合値（フィルタとウィンドウ）の候補を設定する。
まず、ｉ番目の測定条件から得られた図８の周波数分析結果から、適合値の候補を設定する例を示す。
図８の筒内圧センサ波形の周波数分析結果と図１０のノック共鳴モードからノック周波数は、
約８［ｋＨｚ］（ρ１０モード）
約１３［ｋＨｚ］（ρ２０モード）
約１５［ｋＨｚ］（ρ０１モード）
約１７［ｋＨｚ］（ρ３０モード）
の４つと判断できる。
この判断の一例を図１１のフローチャートにて開示する。
但し、約１３［ｋＨｚ］と約１７［ｋＨｚ］とのモードは、前記筒内圧センサ３で検出が不完全なノック共鳴モードであるため、図８のスケールでは見えない。
図８のノックセンサ波形の周波数分析結果からノック発生状態で現れる周波数ピークは、
約７〜８［ｋＨｚ］；１０〜８０［°ＣＡ  ＡＴＤＣ］
約１３［ｋＨｚ］  ；１０〜４０［°ＣＡ  ＡＴＤＣ］
約１７［ｋＨｚ］  ；１０〜４０［°ＣＡ  ＡＴＤＣ］
の３つであり、これらが前記筒内圧センサ３の周波数分析結果からノックを入力源とする応答であると判断できる。
一方、
約１３［ｋＨｚ］  ；４０［°ＣＡ  ＡＴＤＣ］以降、
約１７［ｋＨｚ］  ；４０［°ＣＡ  ＡＴＤＣ］以降
は、ノック発生状態とノック非発生状態との周波数ピークの強度の差が小さいので、ノイズの影響が大きいことが判る。
この判断のフローチャートの一例を図１２に開示する。
以上のノック周波数とノック発生期間の傾向から、ｉ番目の測定条件に関する適合値の候補は、図１３のように設定することができる。
ｉ番目の測定条件と同様に、全ｋ個の各測定条件についても適合値の候補を設定する。
前記適合値設定手段１６では、各測定条件で別個に図１３に示すような適合値の候補の設定が可能であるが、ノック制御仕様によっては、適合値の設定自由度に制限があるので、適合値の候補を設定する際はその点を考慮する。
上述の「適合値の設定自由度に制限」とは、ノック制御仕様によっては、全エンジン回転数で共通の適合値、全気筒で共通の適合値などの制限が設けられている、ことを指している。
なお、この発明においては、理解を容易とするために、各測定条件について共通のｎ個の適合値の候補を設定した、として説明を進める。
【００２０】
前記振動強度算出手段１８について追記する。
この振動強度算出手段１８では、ｉ番目の測定条件で検出されたノックセンサ波形とｒ番目の適合値の候補から、振動強度Ｉｋ（ｉ、ｒ）を算出する（ｉ＝１〜ｋ、ｒ＝１〜ｎ）。
この振動強度算出手段１８による処理のフローチャートを図１５に開示する。
この図１５のフローチャートの補足説明に示すように、振動強度Ｉｋ（ｉ、ｒ）は振動強度の指標を表す複数の値の集合である。
そして、振動強度Ｉｋ（ｉ、ｒ）の指標は、使用するノック制御に合わせて選択する。
一般的には、図１４に示すように、最大値や積分値が用いられる。
【００２１】
前記圧力変動強度算出手段１９について追記する。
この圧力変動強度算出手段１９では、ノック発生状態の筒内圧センサ波形に、図１６（ａ）に示す可聴域ノックの周波数帯域と図１６（ｂ）に示す広帯域ノックの周波数帯域をそれぞれ抽出できるフィルタを用いたフィルタ処理、および適切なウィンドウ処理を行い、例として図１７に示すような最大値や積分値を用いたｉ番目の測定条件における可聴域ノックおよび広帯域ノックに対応する圧力変動強度Ｉｃｈ（ｉ）、Ｉｃｗ（ｉ）（以下、単に「圧力変動強度Ｉｃｈ（ｉ）、Ｉｃｗ（ｉ）」という。）を求める。
この圧力変動強度Ｉｃｈ（ｉ）、Ｉｃｗ（ｉ）は、ｋ個の各測定条件について算出する（ｉ＝１〜ｋ）。
前記圧力変動強度算出手段１９による処理のフローチャートを図１９に開示する。
この図１９のフローチャートの補足説明に示すように、圧力変動強度Ｉｃｈ（ｉ）、Ｉｃｗ（ｉ）は、圧力変動強度を表す複数の値の集合である。
圧力変動強度Ｉｃｈ（ｉ）、Ｉｃｗ（ｉ）を算出する際のフィルタ（可聴域・広帯域のフィルタ）およびウィンドウは、各測定条件により変更しても良いが、圧力変動として検出されるノック周波数とノック発生タイミングが各測定条件によりほとんど変化しない場合は、フィルタとウィンドウとは各測定条件で同じ設定とすれば良い。
可聴域ノックを抽出するフィルタは、実験時のノック強度評価方法として一般的な聴感評価との整合性確保、および製品段階において窓を開けた状態での走行時などにノック音がドライバーに聞こえることによる品質の問題への対応性を評価する。
広帯域ノックを抽出するフィルタは、圧力波形そのものが大きいノックによる耐久性の問題への対応性を評価する。
前記筒内圧センサ３では、検出が不完全なノック共鳴モードが存在する場合について、ｉ番目の測定条件における圧力変動強度Ｉｃｈ（ｉ）、Ｉｃｗ（ｉ）の算出方法を図１８式１、式２に開示する。
この図１８は図８の筒内圧センサ３の時系列周波数分析のウィンドウ（ｂ）を拡大したものである。
各ノック周波数８、１３、１５、１７［ｋＨｚ］に対応する圧力変動強度をＳｎ（Ｓ１、Ｓ２、Ｓ３、Ｓ４、図１８ではｎ＝１〜４）とする。
圧力変動強度Ｓｎは、各ノック周波数を抽出できるフィルタ処理および適切なウィンドウ処理後、最大値や積分値などとして算出する。
この圧力変動強度Ｓｎに掛かる重み係数αｎ（図１８ではｎ＝１〜４）とし、可聴域ノックに関連する圧力変動強度をＳ１、Ｓ２とし、広帯域ノックに関連する圧力変動強度をＳ１〜Ｓ４とすると、圧力変動強度Ｉｃｈ（ｉ）、Ｉｃｗ（ｉ）はそれぞれ以下の式１と式２とで表すことができる。
【数１】
2011179322.tif 000003
【数２】
2011179322.tif 000004
重み係数αｎは、前記筒内圧センサ３での検出度合いに応じて適切な値とする。
例えば、ノック周波数１３、１７［ｋＨｚ］は、前記適合値設定手段１６で述べたように筒内圧センサ３で検出が不完全なノック共鳴モードであるので、それぞれの検出度合いに応じてα２Ｓ２とα４Ｓ４との重み係数α２とα４とをα１とα３より小さい値とすれば、検出が不完全なノック共鳴モードの存在を考慮して圧力変動強度Ｉｃｈ（ｉ）、Ｉｃｗ（ｉ）を算出することができる。
【００２２】
前記相関係数算出手段２０について追記する。
この相関係数算出手段２０では、図２０に示すように、ｉ番目の測定条件、ｒ番目の適合値の候補に対応する振動強度Ｉｋ（ｉ、ｒ）と、ｉ番目の測定条件に対応する圧力変動強度Ｉｃｈ（ｉ）、Ｉｃｗ（ｉ）に関する２つの相関図より、可聴域ノックおよび広帯域ノックに関する２つの相関係数Ｃｈ（ｉ、ｒ）、Ｃｗ（ｉ、ｒ）が得られる（ｉ＝１〜ｋ、ｒ＝１〜ｎ）。
なお、上述の「相関図」とは、同時に計測した筒内圧センサ波形とノックセンサ波形より算出した圧力変動強度と振動強度から作成された、相関図である。そして、相関図のプロット数は、相関係数が問題なく算出できるように設定した測定波形の数に一致する。
前記相関係数算出手段２０による処理のフローチャートを図２１に開示する。
例えば、図２０に示す如く、相関係数のＣｗ（ｉ、ｒ）よりＣｈ（ｉ、ｒ）が大きな値であれば、
適合値のある候補（フィルタ・ウィンドウ）は、可聴域ノックの方が広帯域ノックより検出性が良いことを示している。
【００２３】
前記ノック検出性評価指標算出手段１７について追記する。
このノック検出性評価指標算出手段１７では、２つの相関係数Ｃｈ（ｉ、ｒ）、Ｃｗ（ｉ、ｒ）からｉ番目の測定条件、ｒ番目の適合値の候補に対応するノック検出性評価指標Ｉｎ（ｉ、ｒ）を算出する（ｉ＝１〜ｋ、ｒ＝１〜ｎ）。
図２２の式３
【数３】
2011179322.tif 000005
のように、２つの相関係数Ｃｈ（ｉ、ｒ）、Ｃｗ（ｉ、ｒ）の和をノック検出性評価指標Ｉｎ（ｉ、ｒ）としても良いが、低回転では品質に関わる可聴域ノック、高回転では耐久性に関わる広帯域ノックが問題となるケースが多く、品質と耐久性とのどちらを重視するべきかを回転数に応じて対しすることが望ましい。
このような場合、図２２の式４
【数４】
2011179322.tif 000006
のように、２つの相関係数Ｃｈ（ｉ、ｒ）、Ｃｗ（ｉ、ｒ）に回転数依存の重み係数βｈ、βｗをそれぞれ掛けて加算したものをノック検出性評価指標Ｉｎ（ｉ、ｒ）とする。
回転数依存の重み係数βｈ、βｗは、例えば図２２の式５
【数５】
2011179322.tif 000007
、式６
【数６】
2011179322.tif 000008
のように、適合回転数により可変とすれば良い。
算出した全てのノック検出性評価指標Ｉｎ（ｉ、ｒ）を総合的に比較し、最もノック検出性評価指標の高い候補を最良適合値とする。
実際には、制御仕様により適合値の設定自由度に差があり、全測定条件において最もノック検出性評価指標の高い設定が不可能な場合がある。
そのような場合は、全測定条件において平均的にノック検出性評価指標の高い条件を選定するか、あるいは実際の運転状況で重要な条件（例えば、ＣＶＴ車であれば、エンジン回転数を低く保ち、それにより負荷は高くなるため、低回転域のノック制御性が重要）において、ノック検出性評価指標が高い条件を選定すれば良い。
前記ノック検出性評価指標算出手段１７による処理のフローチャートを図２３に開示する。
【００２４】
前記ノック検出性評価指標算出手段１７は、ノック検出性評価指標の算出と、このノック検出性評価指標を用いて適合値の候補の中から最良適合値の選定とを行う。
従って、２つまたはそれ以上の相関係数に重み係数を掛けてノック検出性評価指標を算出することで、可聴域（品質）および広帯域（耐久性）ノック、またはその他の問題の原因となるノックのいずれかを重視することが可能となる。
【００２５】
以上、説明したノック制御の前記適合装置１により、この適合装置１は、前記エンジン２の気筒内部の圧力を検知する筒内圧センサ３と、前記エンジン２の側壁の任意の位置に取り付けて振動を検知するノックセンサ４と、これら筒内圧センサ３とノックセンサ４とから出力させる検知信号を波形信号として測定し処理する処理装置５とを備え、その処理に基づいてノックセンサ４の適合値の候補の選定を行うノック制御の適合方法において、筒内圧センサ波形とノックセンサ波形を同時測定するステップと、これら測定された波形それぞれに対して時系列の周波数特性を算出するステップと、この算出された周波数特性に基づく適合値の候補の選定を行うステップと、最終的に適合値の候補の中から最良な適合値を選定するステップとを有する。
従って、ノック周波数とノック発生時間の全体の傾向を把握し、詳細にノック検出性の評価を行う適合値の候補を絞り込むことで適合工数の削減が可能となる。
詳述すれば、上述のノック制御の適合方法においては、適合値の最適化により、複雑なノック制御装置を用いることなく、高いノック検出性を実現するための方法として、図１に示す適合方法とした。
【００２６】
この適合方法においては、同時測定されたノックセンサ波形における適合値の候補に対する振動強度を算出するステップと、筒内圧センサ波形に対し可聴域や広帯域のフィルタ（あるいはその他の通過帯域を持つフィルタ）を用いてそれぞれの波形領域に応じた複数の圧力変動強度を算出するステップと、この複数の圧力変動強度と振動強度との相関係数を算出するステップとを行う。
従って、圧力変動強度算出時に、前記筒内圧センサ３では検出が不完全なノック共鳴モードを考慮することで、相関係数算出の精度を向上させることができる。
また、可聴域と広帯域ノック、その他のノックを定義することで、可聴ノック音による品質の問題や広帯域ノックによる耐久性の問題、その他のノックによる問題に、それぞれ個別に対処することが可能である。
【００２７】
次に、相関係数を用いて式４によりノック検出性評価指標を算出するステップと、このノック検出性評価指標を用いて適合値の候補の中から最良適合値を選定するステップとを有する。
従って、２つまたはそれ以上の相関係数に重み係数を掛けてノック検出性評価指標を算出することで、可聴域（品質）および広帯域（耐久性）ノック、またはその他の問題の原因となるノックのいずれかを重視することが可能となる。
【００２８】
ここで、筒内圧センサ波形とノックセンサ波形とを同時計測したデータから、ノックセンサ取付位置における適合値の最適化を行うまでの実施例を説明する。
一例として、エンジン回転数１６００ｒｐｍ、３番気筒＃３の測定条件での適合値を最適化する。
この測定条件の識別番号を「ｉ＝３」とする。
【００２９】
前記周波数分析手段１５により、筒内圧センサ波形とノックセンサ波形に関して、図２４に示すノック発生状態とノック非発生状態との周波数分析結果が得られる。
【００３０】
前記適合値設定手段１６により、ノック周波数は、図１０のノック共鳴モードと図２４の筒内圧センサ波形の周波数分析結果から、
約８［ｋＨｚ］（ρ１０モード）
約１３［ｋＨｚ］（ρ２０モード）
約１５［ｋＨｚ］（ρ０１モード）
約１７［ｋＨｚ］（ρ３０モード）
の４つと判断できる。
但し、約１３［ｋＨｚ］と約１７［ｋＨｚ］のモードは、図２４より前記筒内圧センサ３で検出が不完全なノック共鳴モードと見なせる。
また、ノックセンサ波形の周波数分析結果より、適合値の候補（１）及び（２）を以下のように設定した。
このとき、各適合値の候補（１）及び（２）の識別番号をそれぞれｒ＝１、ｒ＝２とする。
適合値の候補（１）：フィルタ中心周波数８．０［ｋＨｚ］、
ウィンドウ２０〜５０［°ＣＡ  ＡＴＤＣ］
適合値の候補（２）：フィルタ中心周波数１７．０［ｋＨｚ］、
ウィンドウ２０〜５０［°ＣＡ  ＡＴＤＣ］
【００３１】
前記振動強度算出手段１８により、適合値の候補（１）及び（２）に対応する振動強度Ｉｋ（３、１）、Ｉｋ（３、２）を算出する。
【００３２】
前記圧力変動強度算出手段１９により、可聴域ノックおよび広帯域ノックに対応する圧力変動強度Ｉｃｈ（３）、Ｉｃｗ（３）を算出する。
図２５は図２４の前記筒内圧センサ３の周波数分析結果（クランク角１０〜４０［°ＣＡ  ＡＴＤＣ］）を拡大したものである。
図２５において、約８、１３、１５、１７［ｋＨｚ］の各ノック周波数に対する圧力変動強度をＳ１、Ｓ２、Ｓ３、Ｓ４とする。
図２５中の約１３、１７［ｋＨｚ］の周波数ピークのゲインは非常に小さく、ここでは前記筒内圧センサ３では全く検出することができないノック共鳴モードと判断する。
従って、圧力変動強度Ｓ２、Ｓ４に対する重み係数α２、α４を０とし、図２５中に示すように可聴域ノックおよび広帯域ノックに対応する周波数帯域を考慮すると、圧力変動強度Ｉｃｈ（３）、Ｉｃｗ（３）を算出する式は、以下の図２５中の式７と式８となる。
【数７】
2011179322.tif 000009
【数８】
2011179322.tif 000010
【００３３】
前記相関係数算出手段２０により、適合値の候補（１）に対しては振動強度Ｉｋ（３、１）と圧力変動強度Ｉｃｈ（３）、および振動強度Ｉｋ（３、１）と圧力変動強度Ｉｃｗ（３）に関するそれぞれ２つの相関図（図２６の（ａ）及び（ｂ）参照。）より、２つの相関係数Ｃｈ（３，１）、Ｃｗ（３、１）が得られる。
適合値の候補（２）に対しても同様に、振動強度Ｉｋ（３、２）と圧力変動強度Ｉｃｈ（３）、および振動強度Ｉｋ（３、２）と圧力変動強度Ｉｃｗ（３）に関するそれぞれ２つの相関図（図２６の（ｃ）及び（ｄ）参照。）より、２つの相関係数Ｃｈ（３，２）、Ｃｗ（３、２）が得られる。
【００３４】
前記ノック検出性評価指標算出手段１７により、適合値の候補（１）及び（２）に対応するノック検出性評価指標Ｉｎ（３，１）、Ｉｎ（３，２）を算出する。
ここで、使用したエンジンの最大回転数（Ｒｍａｘ）と最小回転数（Ｒｍｉｎ）とをそれぞれ以下のようにする。
Ｒｍａｘ＝６０００（ｒｐｍ）
Ｒｍｉｎ＝６００（ｒｐｍ）
また、適合回転数（Ｒａ）は、
Ｒａ＝１６００（ｒｐｍ）
なので、回転数依存の重み係数βｈ、βｗは、図２２の式５及び式６より、以下のようになる。
βｈ＝５０００／６６００
βｗ＝１６００／６６００
このとき、上記の重み係数βｈ、βｗの値と図２６に示した相関係数Ｃｈ（３、１）、Ｃｗ（３、１）およびＣｈ（３、２）、Ｃｗ（３、２）の値をそれぞれ図２２の式４に代入すると、以下のように適合値の候補（１）、（２）に対応するノック検出性評価指標Ｉｎ（３、１）、Ｉｎ（３、２）の値が求まる。
Ｉｎ（３、１）＝０．８５
Ｉｎ（３、２）＝０．７５
そして、
Ｉｎ（３、１）＞Ｉｎ（３、２）
であるので、この実施例で使用したエンジン回転数１６００（ｒｐｍ）、３番気筒＃３の測定条件については、適合値の候補（１）を採用し、最良適合値とすれば良いことが分かる。
ここで、仮に、適合回転数（Ｒａ）がＲａ＝５０００（ｒｐｍ）のときに図２６に示した相関係数Ｃｈ（３、１）、Ｃｗ（３、１）、Ｃｈ（３、２）、Ｃｗ（３、２）が得られたとする。
Ｒａ＝１６００（ｒｐｍ）のときと同様の手順で、Ｒａ＝５０００（ｒｐｍ）のときのノック検出性評価指標Ｉｎ（３、１）、Ｉｎ（３、２）の値を求めると、以下の値となる。
Ｉｎ（３、１）＝０．７５
Ｉｎ（３、２）＝０．８１
このため、
Ｉｎ（３、１）＜Ｉｎ（３、２）
であるので、エンジン回転数５０００（ｒｐｍ）、３番気筒＃３の測定条件については、適合値の候補（２）を採用し、最良適合値とすれば良い。
【００３５】
次に、作用を説明する。
【００３６】
先ず、図１のノック制御の適合方法の処理内容のフローチャートに沿って説明する。
このノック制御の適合方法の処理内容用のプログラムがスタートすると、前記波形測定手段８で各測定条件について筒内圧センサ波形と各ノックセンサ波形を同時測定する処理（１０１）を行う。
そして、この処理（１０１）によって、ノックに起因する振動の特定が可能となる。
また、前記波形測定手段８で各測定条件について筒内圧センサ波形と各ノックセンサ波形を同時測定する処理（１０１）の後には、前記周波数分析手段１５により、ｋ個の各測定条件について、筒内圧センサ波形とノックセンサ波形について時系列の周波数分析を行う処理（１０２）に移行する。
この処理（１０２）によって、周波数分析手段１５は、ノック発生状態で測定した筒内圧センサ波形、ノック非発生状態で測定した筒内圧センサ波形、ノック発生状態で測定したノックセンサ波形、ノック非発生状態で測定したノックセンサ波形に対して時系列の周波数分析を行う。
図９に開示した前記周波数分析手段１５による時系列の周波数分析用処理のフローチャートは後述する。
なお、この図９の時系列の周波数分析用処理のフローチャートにおいては、図６に示す如く、同時測定したクランク角センサ波形を用いて測定波形１つ１つに対して時系列にウィンドウ（この例では、クランク角１０度おきに６つのウィンドウを設定、ウィンドウ幅はクランク角３０度）を設定し、ＦＦＴやＤＦＴを用いて図７のように周波数特性を求める。
このとき、ウィンドウ幅や間隔、および隣り合うウィンドウを重ね合わせるどうかについては限定しない。
最後に、全サイクルの測定波形の周波数分析結果をウィンドウ毎に平均し、図８のように前記筒内圧センサ３と前記ノックセンサ４のそれぞれについてノック発生状態と非発生状態との結果を比較する形とする。
そして、前記周波数分析手段１５によりｋ個の各測定条件について、筒内圧センサ波形とノックセンサ波形について時系列の周波数分析を行う処理（１０２）の後には、前記適合値設定手段１６で周波数分析の結果から各測定条件における適合値の選定を行う処理（１０３）に移行する。
つまり、前記周波数分析手段１５による処理（１０２）の結果からノック周波数とノック発生時間の傾向を把握することができるため、この処理（１０３）において、前記適合値設定手段１６で詳細にノック検出性の評価を行う適合条件の候補を絞り込むことができる。
この周波数分析手段１５による時系列の周波数分析の処理を全ての測定条件について行う。
結果として、図８のような周波数分析結果は測定条件毎に求める。
このとき、上述の処理（１０２）及び処理（１０３）においては、適合工数の削減が可能となる。
前記適合値設定手段１６で周波数分析の結果から各測定条件における適合値の選定を行う処理（１０３）の後には、前記振動強度算出手段１８で各測定条件のノックセンサ波形と適合値の各候補から振動強度Ｉｋ（ｉ、ｒ）を算出する処理（１０４）に移行する。
この処理（１０４）において、振動強度算出手段１８は、前の適合値設定手段１６にて選定した各測定条件における適合値について、ノック発生状態で測定したノックセンサ波形から振動強度を算出する。
図１５に開示した振動強度算出手段１８による処理のフローチャートは後述する。
図１４に、図１３のノックセンサ（１）の適合値の候補に対応する振動強度算出の例を示す。
ノック発生状態で測定したノックセンサ波形に通過帯の中心周波数が７．５［ｋＨｚ］のフィルタ処理、続いて抽出範囲が１０〜８０［°ＣＡ  ＡＴＤＣ］のウィンドウ処理を行い、その結果得られた波形から振動強度を算出する。
前記振動強度算出手段１８で各測定条件のノックセンサ波形と適合値の各候補から振動強度Ｉｋ（ｉ、ｒ）を算出する処理（１０４）の後には、前記圧力変動強度算出手段１９により各測定条件で筒内圧センサ波形から可聴域ノックおよび広帯域ノックに対応する圧力変動強度Ｉｃｈ（ｉ）、Ｉｃｗ（ｉ）を算出する処理（１０５）に移行する。
この処理（１０５）において、圧力変動強度算出手段１９は、ノック発生状態およびノック非発生状態で測定した筒内圧センサ波形に、図１６（ａ）に示す可聴域フィルタおよび図１６（ｂ）に示す広帯域フィルタを用いてそれぞれの圧力変動強度Ｉｃｈ（ｉ）、Ｉｃｗ（ｉ）を算出する。
図１９に開示した前記圧力変動強度算出手段１９による処理のフローチャートは後述する。
前記圧力変動強度算出手段１９により各測定条件で筒内圧センサ波形から可聴域ノックおよび広帯域ノックに対応する圧力変動強度Ｉｃｈ（ｉ）、Ｉｃｗ（ｉ）を算出する処理（１０５）の後には、前記相関係数算出手段２０で、処理（１０４）からの前記振動強度算出手段１８による振動強度Ｉｋ（ｉ、ｒ）と処理（１０５）からの前記圧力変動強度算出手段１９による圧力変動強度Ｉｃｈ（ｉ）、Ｉｃｗ（ｉ）とから各々の相関係数Ｃｈ（ｉ、ｒ）、Ｃｗ（ｉ、ｒ）を算出する処理（１０６）に移行する。
この処理（１０６）において、相関係数算出手段２０は、図２０に示すように、ｉ番目の測定条件、ｒ番目の適合値の候補に対応する振動強度Ｉｋ（ｉ、ｒ）と、ｉ番目の測定条件に対応する圧力変動強度Ｉｃｈ（ｉ）、Ｉｃｗ（ｉ）に関する２つの相関図より、可聴域ノックおよび広帯域ノックに関する２つの相関係数Ｃｈ（ｉ、ｒ）、Ｃｗ（ｉ、ｒ）が得られる（ｉ＝１〜ｋ、ｒ＝１〜ｎ）。
図２１に開示した前記相関係数算出手段２０による処理のフローチャートは後述する。
前記相関係数算出手段２０で、処理（１０４）からの前記振動強度算出手段１８による振動強度Ｉｋ（ｉ、ｒ）と処理（１０５）からの前記圧力変動強度算出手段１９による圧力変動強度Ｉｃｈ（ｉ）、Ｉｃｗ（ｉ）とから各々の相関係数Ｃｈ（ｉ、ｒ）、Ｃｗ（ｉ、ｒ）を算出する処理（１０６）の後には、前記ノック検出性評価指標算出手段１７で、前記相関係数算出手段２０による相関係数Ｃｈ（ｉ、ｒ）、Ｃｗ（ｉ、ｒ）から各々のノック検出性評価指標Ｉｎ（ｉ、ｒ）を算出する処理（１０７）に移行する。
この処理（１０７）において、前記ノック検出性評価指標算出手段１７は、図２２の式４により相関係数Ｃｈ（ｉ、ｒ）、Ｃｗ（ｉ、ｒ）からノック検出性評価指標Ｉｎ（ｉ、ｒ）を算出する。
図２３に開示した前記ノック検出性評価指標算出手段１７による処理のフローチャートは後述する。
前記ノック検出性評価指標算出手段１７で、前記相関係数算出手段２０による相関係数Ｃｈ（ｉ、ｒ）、Ｃｗ（ｉ、ｒ）からノック検出性評価指標Ｉｎ（ｉ、ｒ）を算出する処理（１０７）の後には、前記適合装置１で、全てのノック検出性評価指標Ｉｎ（ｉ、ｒ）を比較して最良適合値を選定する処理（１０８）に移行する。
この処理（１０８）において、前記適合装置１は、（ｋ  ｘ  ｎ）個の全てのノック検出性評価指標Ｉｎ（ｉ、ｒ）を総合的に比較し、最良適合値を選定する。
【００３７】
図９に開示した前記周波数分析手段１５による周波数分析用フローチャートを説明する。
この周波数分析手段１５による周波数分析用プログラムがスタート（「処理開始」ともいう。）（２０１）すると、測定条件（ｉ＝１〜ｋ）を選択する処理（２０２）に移行する。
そして、この処理（２０２）においては、回転数、気筒によって測定条件を選択する。
このとき、測定条件として負荷を追加しても良い。
この測定条件（ｉ＝１〜ｋ）を選択する処理（２０２）の後には、波形の種類を選択する処理（２０３）に移行する。
つまり、この処理（２０３）においては、「筒内圧、ノックセンサ、各々ノック発生／非発生」、追記すれば、ノック発生状態の筒内圧センサ波形とノックセンサ波形、およびノック非発生状態の筒内圧センサ波形と各ノックセンサ波形の４種類を選択する。
上述の波形の種類を選択する処理（２０３）の後には、一波形に対して時系列にウィンドウを設定し周波数分析を行う処理（２０４）に移行する。
この処理（２０４）においては、前記周波数分析手段１５で、ノック発生状態の筒内圧センサ波形、ノック非発生状態の筒内圧センサ波形、ノック発生状態のノックセンサ波形、ノック非発生状態のノックセンサ波形に対して時系列の周波数分析を行う。
一波形に対して時系列にウィンドウを設定し周波数分析を行う処理（２０４）の後には、全サイクルの波形の周波数分析が完了したか否かの判断（２０５）に移行する。
このとき、全サイクルの波形は、相関係数が問題なく算出できるように設定した測定波形数の全波形を指す。
そして、判断（２０５）がＮＯの場合には、上述の一波形に対して時系列にウィンドウを設定し周波数分析を行う処理（２０４）に戻る。
判断（２０５）がＹＥＳの場合には、各ウィンドウについて、全サイクルの波形の周波数分析結果を平均する処理（２０６）に移行する。
この各ウィンドウについて、全サイクルの波形の周波数分析結果を平均する処理（２０６）の後には、全種類の波形の（筒内圧、ノックセンサ、各々ノック発生／非発生）周波数分析が完了したか否かの判断（２０７）に移行する。
このとき、選択する波形は、
２＋２種類
であり、ノック発生状態の筒内圧センサとノックセンサ波形およびノック非発生状態の筒内圧センサ波形とノックセンサ波形、図８の例であれば、
２＋２＝４
となる。
そして、全種類の波形の（筒内圧、ノックセンサ、各々ノック発生／非発生）周波数分析が完了したか否かの判断（２０７）において、この判断（２０７）がＮＯの場合には、上述の波形の種類を選択する処理（２０３）に戻る。
判断（２０７）がＹＥＳの場合には、全測定条件（回転数、気筒）の周波数分析が完了したか否かの判断（２０８）に移行する。
このとき、回転数６段階、４気筒エンジンならば、測定条件の数は、
６Ｘ４＝２４
となる。
また、全測定条件の（回転数、気筒）周波数分析が完了したか否かの判断（２０８）において、この判断（２０８）がＮＯの場合には、上述の測定条件を選択する処理（２０２）に戻る。
判断（２０８）がＹＥＳの場合には、前記周波数分析手段１５による周波数分析用プログラムのエンド（「処理終了」ともいう。）（２０９）に移行する。
【００３８】
図１１の前記適合値設定手段１６におけるノック周波数の判断用フローチャートに沿って説明する。
この適合条件設定手段１６におけるノック周波数の判断用プログラムがスタート（「処理開始」ともいう。）（３０１）すると、前記筒内圧センサ３の周波数分析結果から周波数ピークを抽出する処理（３０２）に移行する。
そして、この筒内圧センサ３の周波数分析結果から周波数ピークを抽出する処理（３０２）の後には、周波数ピークの強度について、ノック発生状態とノック非発生状態との差を算出する処理（３０３）に移行する。
処理（３０３）において、周波数ピークの強度について、ノック発生状態とノック非発生状態との差を算出した後には、ノック発生状態とノック非発生状態との差が規定値より大きいか否かの判断（３０４）に移行する。
なお、上述の処理（３０２）及び処理（３０３）においては、図８の筒内センサ３の周波数分析結果を参照する。
そして、このノック発生状態とノック非発生状態との差が規定値より大きいか否かの判断（３０４）において、判断（３０４）がＮＯの場合には、後述するノックではないと判断する処理（３０８）に移行する。
判断（３０４）がＹＥＳの場合には、ピーク周波数を算出する処理（３０５）を経て、ピーク周波数がノック共振モードに相当するか否かの判断（３０６）に移行する。
このピーク周波数がノック共振モードに相当するか否かの判断（３０６）において、判断（３０６）がＮＯの場合には、このノックでないと判断する処理（３０８）に移行する。
判断（３０６）がＹＥＳの場合には、ノックどあると判断する処理（３０７）に移行する。
【００３９】
図１２の前記適合値設定手段１６における適合値の候補選定の判断用フローチャートに沿って説明する。
この適合値設定手段１６における適合値の候補選定の判断用プログラムがスタート（「処理開始」ともいう。）（４０１）すると、前記ノックセンサ４の周波数分析結果から周波数ピークを抽出する処理（４０２）に移行する。
そして、このノックセンサ４の周波数分析結果から周波数ピークを抽出する処理（４０２）の後には、周波数ピークの強度について、ノック発生状態とノック非発生状態との差を算出する処理（４０３）に移行する。
処理（４０３）において、周波数ピークの強度について、ノック発生状態とノック非発生状態との差を算出した後には、ノック発生状態とノック非発生状態との差が規定値より大きいか否かの判断（４０４）に移行する。
なお、上述の処理（４０２）及び処理（４０３）においては、図８のノックセンサ４の周波数分析結果を参照する。
そして、このノック発生状態とノック非発生状態との差が規定値より大きいか否かの判断（４０４）において、判断（４０４）がＮＯの場合には、後述する適合値の候補から除外（ノイズの影響大）する処理（４０８）に移行する。
判断（４０４）がＹＥＳの場合には、ピーク周波数を算出する処理（４０５）を経て、ピーク周波数が筒内圧センサ３のピーク周波数に対応しているか否かの判断（４０６）に移行する。
このピーク周波数が筒内圧センサ３のピーク周波数に対応しているか否かの判断（４０６）において、判断（４０６）がＮＯの場合には、適合値の候補から除外（ノイズの影響大）する処理（４０８）に移行する。
判断（４０６）がＹＥＳの場合には、適合値の候補に設定する処理（４０７）に移行する。
【００４０】
図１５の前記振動強度算出手段１８における振動強度算出処理用のフローチャートに沿って説明する。
この振動強度算出手段１８における振動強度算出処理用プログラムがスタート（「処理開始」ともいう。）（５０１）すると、測定条件（ｉ＝１〜ｋ）を選択する処理（５０２）に移行する。
そして、この処理（５０２）においては、回転数、気筒によって測定条件を選択する。
なお、この処理（５０２）において、測定条件として負荷を追加しても良い。
この測定条件（ｉ＝１〜ｋ）を選択する処理（５０２）の後には、
適合値の候補（ｒ＝１〜ｎ）を選択する処理（５０３）に移行する。
そして、この適合値の候補（ｒ＝１〜ｎ）を選択する処理（５０３）の後には、
ノック発生状態の振動１波形に対してフィルタおよびウィンドウ処理を行う処理（５０４）に移行する。
このノック発生状態の振動１波形に対してフィルタおよびウィンドウ処理を行う処理（５０４）の後には、
振動強度Ｉｋ（ｉ、ｒ）を算出（最大値、積分値など）する処理（５０５）に移行する。
上述の処理（５０４）及び処理（５０５）は、図１４に示すように行われ、振動１波形に対して１つの振動強度（積分値など）が算出される。
そして、振動強度Ｉｋ（ｉ、ｒ）を算出（最大値、積分値など）する処理（５０５）の後には、全サイクルの波形の処理が完了したか否かの判断（５０６）に移行する。
この全サイクルの波形の処理が完了したか否かの判断（５０６）おいて、全サイクルの波形の処理が完了すると、振動強度Ｉｋ（ｉ、ｒ）は全サイクル数分の振動強度の集合となる。
なお、「全サイクルの波形」とは、相関係数が問題なく算出できるように設定した測定波形数の全波形を指す。
上述の全サイクルの波形の処理が完了したか否かの判断（５０６）おいて、判断（５０６）がＮＯの場合には、上述のノック発生状態の振動１波形に対してフィルタおよびウィンドウ処理を行う処理（５０４）に戻る。
判断（５０６）がＹＥＳの場合には、適合値の全候補の処理が完了したか否かの判断（５０７）に移行する。
そして、この適合値の全候補の処理が完了したか否かの判断（５０７）において、判断（５０７）がＮＯの場合には、上述の適合値の候補（ｒ＝１〜ｎ）を選択する処理（５０３）に戻る。
判断（５０７）がＹＥＳの場合には、全測定条件の（回転数、気筒）処理が完了したか否かの判断（５０８）に移行する。
この全測定条件の（回転数、気筒）処理が完了したか否かの判断（５０８）において、判断（５０８）がＮＯの場合には、上述の測定条件（ｉ＝１〜ｋ）を選択する処理（５０２）に戻る。
判断（５０８）がＹＥＳの場合には、前記振動強度算出手段１８における振動強度算出処理用プログラムのエンド（「処理終了」ともいう。）（５０９）に移行する。
【００４１】
図１９の前記圧力変動強度算出手段１９における圧力変動強度算出処理用のフローチャートに沿って説明する。
この圧力変動強度算出手段１９における圧力変動強度算出処理用プログラムがスタート（「処理開始」ともいう。）（６０１）すると、測定条件（ｉ＝１〜ｋ）を選択する処理（６０２）に移行する。
そして、この処理（６０２）においては、回転数、気筒によって測定条件を選択する。
なお、この処理（６０２）において、測定条件として負荷を追加しても良い。
この測定条件（ｉ＝１〜ｋ）を選択する処理（６０２）の後には、ノック発生状態の筒内圧１波形に対して可聴域フィルタおよびウィンドウ処理を行う処理（６０３）に移行する。
このノック発生状態の筒内圧１波形に対して可聴域フィルタおよびウィンドウ処理を行う処理（６０３）において、可聴域フィルタを図１６（ａ）に開示している。
また、ノック発生状態の筒内圧１波形に対して可聴域フィルタおよびウィンドウ処理を行う処理（６０３）の後には、圧力変動強度Ｉｃｈ（ｉ）を算出（最大値、積分値など）する処理（６０４）に移行する。
この圧力変動強度Ｉｃｈ（ｉ）を算出（最大値、積分値など）する処理（６０４）は、最大値、積分値などに関して、図１７に示すように処理される。
更に、圧力変動強度Ｉｃｈ（ｉ）を算出（最大値、積分値など）する処理（６０４）の後には、同じ筒内圧１波形に対して広帯域フィルタおよびウィンドウ処理を行う処理（６０５）に移行する。
この同じ筒内圧１波形に対して広帯域フィルタおよびウィンドウ処理を行う処理（６０５）において、広帯域フィルタを図１６（ｂ）に開示している。
更にまた、同じ筒内圧１波形に対して広帯域フィルタおよびウィンドウ処理を行う処理（６０５）の後には、圧力変動強度Ｉｃｗ（ｉ）を算出（最大値、積分値など）する処理（６０６）に移行する。
この圧力変動強度Ｉｃｗ（ｉ）を算出（最大値、積分値など）する処理（６０６）は、最大値、積分値などに関して、図１７に示すように処理される。
そして、筒内圧１波形に対して２つの圧力変動強度（積分値など）が算出される。
上述の圧力変動強度Ｉｃｗ（ｉ）を算出（最大値、積分値など）する処理（６０６）の後には、全サイクルの波形の処理が完了したか否かの判断（６０７）に移行する。
この全サイクルの波形の処理が完了したか否かの判断（６０７）おいて、全サイクルの波形の処理が完了すると、圧力変動強度Ｉｃｈ（ｉ）およびＩｃｗ（ｉ）は全サイクル数分の圧力変動強度の集合となる。
なお、「全サイクルの波形」とは、相関係数が問題なく算出できるように設定した測定波形数の全波形を指す。
そして、全サイクルの波形の処理が完了したか否かの判断（６０７）において、判断（６０７）がＮＯの場合には、上述のノック発生状態の筒内圧１波形に対して可聴域フィルタおよびウィンドウ処理を行う処理（６０３）に戻る。
判断（６０７）がＹＥＳの場合には、全測定条件の（回転数、気筒）処理が完了したか否かの判断（６０８）に移行する。
この全測定条件の（回転数、気筒）処理が完了したか否かの判断（６０８）において、判断（６０８）がＮＯの場合には、上述の測定条件（ｉ＝１〜ｋ）を選択する処理（６０２）に戻る。
判断（６０８）がＹＥＳの場合には、前記圧力変動強度算出手段１９における圧力変動強度算出処理用プログラムのエンド（「処理終了」ともいう。）（６０９）に移行する。
【００４２】
図２１の前記相関係数算出手段２０における処理用のフローチャートに沿って説明する。
この相関係数算出手段２０における処理用プログラムがスタート（「処理開始」ともいう。）（７０１）すると、振動強度Ｉｋ（ｉ、ｒ）の（ｉ、ｒ）の組合せを選択する処理（７０２）に移行する。
そして、この処理（７０２）においては、振動強度Ｉｋ（ｉ、ｒ）について、（ｉ、ｒ）の組合せを選択する。
このとき、ｉ、ｒは以下のように設定する。
ｉ：測定条件（ｉ＝１〜ｋ）
ｒ：適合値の候補（ｒ＝１〜ｎ）
この振動強度Ｉｋ（ｉ、ｒ）の（ｉ、ｒ）の組合せを選択する処理（７０２）の後には、全サイクルの波形の振動強度Ｉｋ（ｉ、ｒ）と圧力変動強度Ｉｃｈ（ｉ）より相関係数Ｃｈ（ｉ、ｒ）を算出する処理（７０３）に移行する。
このとき、「全サイクルの波形」とは、相関係数が問題なく算出できるように設定した測定波形数の全波形を指す。
また、全サイクルの波形の振動強度Ｉｋ（ｉ、ｒ）と圧力変動強度Ｉｃｈ（ｉ）より相関係数Ｃｈ（ｉ、ｒ）を算出する処理（７０３）の後には、全サイクルの波形の振動強度Ｉｋ（ｉ、ｒ）と圧力変動強度Ｉｃｗ（ｉ）より相関係数Ｃｗ（ｉ、ｒ）を算出する処理（７０４）に移行する。
そして、全サイクルの波形の振動強度Ｉｋ（ｉ、ｒ）と圧力変動強度Ｉｃｗ（ｉ）より相関係数Ｃｗ（ｉ、ｒ）を算出する処理（７０４）の後には、（ｉ、ｒ）の全ての組合せの処理が完了したか否かの判断（７０５）に移行する。
この（ｉ、ｒ）の全ての組合せの処理が完了したか否かの判断（７０５）において、判断（７０５）がＮＯの場合には、上述の振動強度Ｉｋ（ｉ、ｒ）の（ｉ、ｒ）の組合せを選択する処理（７０２）に戻る。
判断（７０５）がＹＥＳの場合には、前記相関係数算出手段２０における処理用プログラムのエンド（「処理終了」ともいう。）（７０６）に移行する。
【００４３】
図２３の前記ノック検出性評価指標算出手段１７における処理用のフローチャートに沿って説明する。
このノック検出性評価指標算出手段１７における処理用プログラムがスタート（「処理開始」ともいう。）（８０１）すると、測定条件（ｉ＝１〜ｋ）を選択する処理（８０２）に移行する。
そして、この処理（８０２）においては、回転数、気筒によって測定条件を選択する。
なお、この処理（８０２）において、測定条件として負荷を追加しても良い。
この測定条件（ｉ＝１〜ｋ）を選択する処理（８０２）の後には、適合値の候補（ｒ＝１〜ｎ）を選択する処理（８０３）に移行する。
そして、この適合値の候補（ｒ＝１〜ｎ）を選択する処理（８０３）の後には、相関係数Ｃｈ（ｉ、ｒ）、Ｃｗ（ｉ、ｒ）よりノック検出性評価指標Ｉｎ（ｉ、ｒ）を算出する処理（８０４）に移行する。
この相関係数Ｃｈ（ｉ、ｒ）、Ｃｗ（ｉ、ｒ）よりノック検出性評価指標Ｉｎ（ｉ、ｒ）を算出する処理（８０４）の後には、適合値の全候補の処理が完了したか否かの判断（８０５）に移行する。
そして、この適合値の全候補の処理が完了したか否かの判断（８０５）において、判断（８０５）がＮＯの場合には、上述の適合値の候補（ｒ＝１〜ｎ）を選択する処理（８０３）に戻る。
判断（８０５）がＹＥＳの場合には、全測定条件の（回転数、気筒）処理が完了したか否かの判断（８０６）に移行する。
この全測定条件の（回転数、気筒）処理が完了したか否かの判断（８０６）において、判断（８０６）がＮＯの場合には、上述の測定条件（ｉ＝１〜ｋ）を選択する処理（８０２）に戻る。
判断（８０６）がＹＥＳの場合には、全てのノック検出性評価指標Ｉｎ（ｉ、ｒ）を比較する処理（８０７）に移行する。
そして、この全てのノック検出性評価指標Ｉｎ（ｉ、ｒ）を比較する処理（８０７）の後には、最良適合値を選定する処理（８０８）に移行し、その後に、前記ノック検出性評価指標算出手段１７における処理用プログラムのエンド（「処理終了」ともいう。）（８０９）に移行する。
【実施例２】
【００４４】
図２７はこの発明の第２実施例を示すものである。
この第２実施例において、上述第１実施例のものと同一機能を果たす箇所には、同一符号を付して説明する。
【００４５】
第２実施例の特徴とするところは、可聴域ノックの定義のみでノックセンサの最良適合値の選定を行う場合に、エンジン２の近くにマイク２１を設置し、筒内圧センサの波形の代わりに、マイク２１の波形を計測して処理を行う点にある。
すなわち、図２７に示す如く、前記エンジン２の近傍にマイク２１を配設し、このマイク２１を前記処理装置５の波形測定手段８に接続して設けるとともに、ノックセンサ４を前記波形測定手段８に接続して設けるものである。
また、可聴域ノックと広帯域ノックとの２種類のノックの定義以外に、他の複数のノックを定義しても構わない。
例えば、新たに「耳に聞こえないノック」を定義すると、図１８に示す各ノック周波数に対する圧力変動強度Ｓ１〜Ｓ４のうちＳ３、Ｓ４を考慮すればよく、そのノックの定義に対する圧力変動強度は、式２においてα１、α２を０とすれば良い。
また、新たに導入したノックの定義の数だけ算出される相関係数の数も増えるので、増えた分の相関係数とそれらに対する重み係数を考慮した項を、ノック検出性評価指標を求める図２２の式４の右辺に第３項以降として追加すれば良い（図２２の式４’参照。）。
さすれば、可聴ノック音による品質の問題や広帯域ノックによる耐久性の問題の他にも、ノックが生じることに因る問題点の対応性を評価することが可能となる。
【符号の説明】
【００４６】
１  ノック制御の適合装置
２  エンジン
２ａ  シリンダブロック
３  筒内圧センサ
３ａ、３ｂ、３ｃ、３ｄ  第１〜第４筒内圧センサ
４  ノックセンサ
５  処理装置
６  ＥＣＵデータ書込装置
７  ノック制御装置（「ＥＣＵ」ともいう。）
８  波形測定手段
９  記録装置
１０  適合手段
１１  クランク角センサ
１２  カム角センサ
１３  第１記録部
１４  第２記録部
１５  周波数分析手段
１６  適合値設定手段
１７  ノック検出性評価指標算出手段
１８  振動強度算出手段
１９  圧力変動強度算出手段
２０  相関係数算出手段
【図１】
2011179322.tif 000011
【図２】
2011179322.tif 000012
【図３】
2011179322.tif 000013
【図４】
2011179322.tif 000014
【図５】
2011179322.tif 000015
【図６】
2011179322.tif 000016
【図７】
2011179322.tif 000017
【図８】
2011179322.tif 000018
【図９】
2011179322.tif 000019
【図１０】
2011179322.tif 000020
【図１１】
2011179322.tif 000021
【図１２】
2011179322.tif 000022
【図１３】
2011179322.tif 000023
【図１４】
2011179322.tif 000024
【図１５】
2011179322.tif 000025
【図１６】
2011179322.tif 000026
【図１７】
2011179322.tif 000027
【図１８】
2011179322.tif 000028
【図１９】
2011179322.tif 000029
【図２０】
2011179322.tif 000030
【図２１】
2011179322.tif 000031
【図２２】
2011179322.tif 000032
【図２３】
2011179322.tif 000033
【図２４】
2011179322.tif 000034
【図２５】
2011179322.tif 000035
【図２６】
2011179322.tif 000036
【図２７】
2011179322.tif 000037