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