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