(57)【要約】【課題】内燃機関のノック発生点火時期推定装置に関し、膨大な工数や時間をかけることなく高い精度でノック発生点火時期を推定できるようにする。【解決手段】異なる複数の点火時期の設定において実際に内燃機関を運転し、内燃機関の自着火に係る状態量を測定する。測定された状態量に基づいて自着火の条件が成立するクランク角を予測し、そのクランク角を第１判定値Clとして算出する。また、測定された状態量に基づいて筒内の未燃燃料量を求め、ノックの発生に必要な量の未燃燃料が存在する限界クランク角を第２判定値Cbとして算出する。点火時期の設定値と第１判定値Clとの線形関係、及び、点火時期の設定値と第２判定値Cbとの線形関係をそれぞれ求め、２つの線形関係に基づいて第１判定値Clと第２判定値Cbとが一致する点火時期をノック発生点火時期として算出する。【選択図】図４
2008095553.tif 000002
【特許請求の範囲】
【請求項１】
点火時期を含む内燃機関の運転条件を設定する運転条件設定手段と、
前記運転条件設定手段において設定された運転条件下で実際に内燃機関を運転し、内燃機関の自着火に係る状態量を測定する状態量測定手段と、
前記状態量測定手段によって測定された状態量に基づいて自着火の条件が成立するクランク角を予測し、そのクランク角を第１判定値として算出する第１判定値算出手段と、
前記状態量測定手段によって測定された状態量に基づいて筒内の未燃燃料量を求め、ノックの発生に必要な量の未燃燃料が存在する限界クランク角を第２判定値として算出する第２判定値算出手段と、
異なる複数の点火時期の設定において前記の第１判定値及び第２判定値を取得し、点火時期の設定値と第１判定値との線形関係、及び、点火時期の設定値と第２判定値との線形関係をそれぞれ求め、２つの線形関係に基づいて第１判定値と第２判定値とが一致する点火時期をノック発生点火時期として算出するノック発生点火時期算出手段と、
を備えることを特徴とする内燃機関のノック発生点火時期推定装置。
【請求項２】
前記ノック発生点火時期算出手段により算出されたノック発生点火時期を所定量だけ遅角し、内燃機関の点火時期の初期値として設定する点火時期初期設定手段と、
実際に内燃機関を運転してノックが発生するか否か検証し、ノックの発生が検知されるまで点火時期の設定値を所定量ずつ進角させていくノック検証手段と、
前記ノック検証手段によってノックの発生が検知された場合、そのときの点火時期の設定値をノック発生点火時期として確定するノック発生点火時期確定手段と、
をさらに備えることを特徴とする請求項１記載の内燃機関のノック発生点火時期推定装置。
【請求項３】
点火時期を含む内燃機関の運転条件を設定する第１のステップと、
前記第１のステップにおいて設定された運転条件下で実際に内燃機関を運転し、内燃機関の自着火に係る状態量を測定する第２のステップと、
前記第２のステップで測定された状態量に基づいて自着火の条件の成否を予測し、自着火条件の成立が予測されるクランク角を第１判定値として算出する第３のステップと、
前記第２のステップで測定された状態量に基づいて筒内の未燃燃料量を求め、ノックの発生に必要な量の未燃燃料が存在する限界クランク角を第２判定値として算出する第４のステップと、
異なる複数の点火時期の設定において前記の第１判定値及び第２判定値を取得し、点火時期の設定値と第１判定値との線形関係、及び、点火時期の設定値と第２判定値との線形関係をそれぞれ求め、２つの線形関係に基づいて第１判定値と第２判定値とが一致する点火時期をノック発生点火時期として算出する第５のステップと、
を含むことを特徴とする内燃機関のノック発生点火時期推定方法。
【請求項４】
前記第５のステップで算出されたノック発生点火時期を所定量だけ遅角し、内燃機関の点火時期の初期値として設定する第６のステップと、
実際に内燃機関を運転してノックが発生するか否か検証し、ノックの発生が検知されるまで点火時期の設定値を所定量ずつ進角させていく第７のステップと、
前記第７のステップでノックの発生が検知された場合、そのときの点火時期の設定値をノック発生点火時期として確定する第８のステップと、
をさらに含むことを特徴とする請求項３記載の内燃機関のノック発生点火時期推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は内燃機関のノック発生点火時期推定装置及び方法に関し、詳しくは、ノックモデルを利用してノック発生点火時期を推定する装置及び方法に関する。
【背景技術】
【０００２】
内燃機関の点火時期制御においては、ノックが発生しないようにしながら点火時期を限界まで進角させてＭＢＴに近付けることが行われている。ノックが発生する点火時期（以下、ノック発生点火時期）は、吸入空気温度、機関温度、機関回転数、負荷等の内燃機関の運転条件に依存する。このため、従来、様々な運転条件の下で実機による適合が行われ、運転条件毎にノック発生点火時期が求められていた。実機による適合では、ノックが発生しないことが確実な点火時期を初期値として、ノックが発生するまで点火時期を徐々に進角していく方法が採られていた。
【０００３】
また、従来、ノック発生点火時期を机上での計算により推定する方法も知られている。ノック発生点火時期の推定には、いわゆるノックモデルが用いられる。ノックモデルには、自着火が発生するタイミングを予測するモデルと、自着火発生タイミングでの未燃燃料の量を計算するモデルとが含まれる。前者のモデルとしては、Livengood-Wuモデルが一般に用いられている。後者のモデルとしては、例えばWiebeモデルを用いることができる。ノックモデルを用いてノック発生点火時期を予測する場合、まず、実機を模したエンジンモデルによって筒内圧の履歴がシミュレートされる。シミュレーションで得られた筒内圧データをノックモデルに入力することで自着火発生タイミングが計算される。自着火発生タイミングにおいてノックの発生に必要な量の未燃燃料が残っていると判定されたならば、そのときの点火時期の設定値がノック発生点火時期であると推定することができる。
【特許文献１】特開２００４−２４５１７３号公報
【特許文献２】特開２００４−３３２５８４号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、従来の実機適合による方法では、適合作業に要する工数や時間は膨大なものとなっていた。これは、予めノック発生点火時期の目安をつけることなく、十分に余裕のある点火時期を初期値としてノック発生点火時期の探索を開始していたことによる。ノック発生点火時期の探索範囲が広くなるほど、適合作業に要する工数や時間は増大してしまう。
【０００５】
一方、従来のノックモデルを用いた方法では、ノック発生点火時期を高い精度で推定することは難しかった。これは、現状のエンジンモデルでは十分に高い精度で筒内圧を計算できないことによる。ノック発生点火時期は筒内圧データを用いて予測されるため、筒内圧データの精度が低いとノック発生点火時期の推定精度も悪化してしまう。なお、エンジンモデルではなく実機を用いて筒内圧を測定することも可能であるが、点火時期を含む全ての運転条件について実機測定を行うとなると、膨大な工数や時間が必要となってノックモデルを用いるメリットが減少してしまう。
【０００６】
ところで、本発明にかかる発明者らは、本発明の創案過程において、自着火条件の成立が予測されるクランク角（自着火クランク角）と点火時期との間には、一定の線形関係が存在することを発見した。また、ノックの発生に必要な量の未燃燃料が存在する限界クランク角と点火時期との間にも、一定の線形関係が存在することを発見した。これらの発見を利用すれば、グラフの横軸に点火時期、縦軸にクランク角をとったとき、自着火クランク角と点火時期との関係を示す直線と、限界クランク角と点火時期との関係を示す直線との交点を計算することで、ノック発生点火時期を求めることができる。各直線は、それぞれ少なくとも２点のデータがあれば近似することができる。
【０００７】
本発明は、上記のような発見に基づいてなされたものであり、その目的は、膨大な工数や時間をかけることなく高い精度でノック発生点火時期を推定できるようにすることにある。
【課題を解決するための手段】
【０００８】
第１の発明は、上記の目的を達成するため、内燃機関のノック発生点火時期推定装置において、
点火時期を含む内燃機関の運転条件を設定する運転条件設定手段と、
前記運転条件設定手段において設定された運転条件下で実際に内燃機関を運転し、内燃機関の自着火に係る状態量を測定する状態量測定手段と、
前記状態量測定手段によって測定された状態量に基づいて自着火の条件が成立するクランク角を予測し、そのクランク角を第１判定値として算出する第１判定値算出手段と、
前記状態量測定手段によって測定された状態量に基づいて筒内の未燃燃料量を求め、ノックの発生に必要な量の未燃燃料が存在する限界クランク角を第２判定値として算出する第２判定値算出手段と、
異なる複数の点火時期の設定において前記の第１判定値及び第２判定値を取得し、点火時期の設定値と第１判定値との線形関係、及び、点火時期の設定値と第２判定値との線形関係をそれぞれ求め、２つの線形関係に基づいて第１判定値と第２判定値とが一致する点火時期をノック発生点火時期として算出するノック発生点火時期算出手段と、
を備えることを特徴としている。
【０００９】
また、第２の発明は、第１の発明において、
前記ノック発生点火時期算出手段により算出されたノック発生点火時期を所定量だけ遅角し、内燃機関の点火時期の初期値として設定する点火時期初期設定手段と、
実際に内燃機関を運転してノックが発生するか否か検証し、ノックの発生が検知されるまで点火時期の設定値を所定量ずつ進角させていくノック検証手段と、
前記ノック検証手段によってノックの発生が検知された場合、そのときの点火時期の設定値をノック発生点火時期として確定するノック発生点火時期確定手段と、
をさらに備えることを特徴としている。
【００１０】
第３の発明は、上記の目的を達成するため、内燃機関のノック発生点火時期推定方法において、
点火時期を含む内燃機関の運転条件を設定する第１のステップと、
前記第１のステップにおいて設定された運転条件下で実際に内燃機関を運転し、内燃機関の自着火に係る状態量を測定する第２のステップと、
前記第２のステップで測定された状態量に基づいて自着火の条件の成否を予測し、自着火条件の成立が予測されるクランク角を第１判定値として算出する第３のステップと、
前記第２のステップで測定された状態量に基づいて筒内の未燃燃料量を求め、ノックの発生に必要な量の未燃燃料が存在する限界クランク角を第２判定値として算出する第４のステップと、
異なる複数の点火時期の設定において前記の第１判定値及び第２判定値を取得し、点火時期の設定値と第１判定値との線形関係、及び、点火時期の設定値と第２判定値との線形関係をそれぞれ求め、２つの線形関係に基づいて第１判定値と第２判定値とが一致する点火時期をノック発生点火時期として算出する第５のステップと、
を含むことを特徴としている。
【００１１】
また、第４の発明は、第３の発明において、
前記第５のステップで算出されたノック発生点火時期を所定量だけ遅角し、内燃機関の点火時期の初期値として設定する第６のステップと、
実際に内燃機関を運転してノックが発生するか否か検証し、ノックの発生が検知されるまで点火時期の設定値を所定量ずつ進角させていく第７のステップと、
前記第７のステップでノックの発生が検知された場合、そのときの点火時期の設定値をノック発生点火時期として確定する第８のステップと、
をさらに含むことを特徴としている。
【発明の効果】
【００１２】
第１の発明及び第３の発明によれば、実際に内燃機関を運転して得られた状態量の測定値（実機データ）に基づいてノック発生点火時期を推定することができる。しかも、点火時期の設定値と各判定値との線形関係を把握することで、ノックが発生しない運転条件での実機データからでも、ノック発生点火時期を精度よく推定することができる。また、点火時期の設定値と各判定値との線形関係は、少なくとも２つの点火時期での実機データがあれば求めることができるので、実機データの取得に多くの工数や時間をかける必要もない。
【００１３】
第２の発明及び第４の発明によれば、設定した点火時期でノックが発生するか否か実際に内燃機関を運転して検証することで、ノック発生点火時期をより精度よく推定することができる。その際、ノック発生点火時期の探索範囲は、ノックモデルを用いて推定したノック発生点火時期を基準として決定されるので、従来の実機適合による方法に比較して、ノック発生点火時期の探索に要する工数や時間を大幅に抑えることができる。
【発明を実施するための最良の形態】
【００１４】
以下、図面を参照して本発明の実施の形態について説明する。図１は本発明の実施の形態としてのノック発生点火時期推定装置の構成を示すブロック図である。ノック発生点火時期推定装置１０は、内燃機関（以下、エンジン）２を運転して種々のデータを取得し、取得したデータを用いて当該エンジン２のノック発生点火時期を推定する装置である。ノック発生点火時期推定装置１０は、エンジン２を制御するドライバ１４と、エンジン２からのデータの取得に用いる各種センサ１６と、種々の計算を行うコンピュータ１２とを含んでいる。
【００１５】
ノック発生点火時期の推定のために必要なデータは、主としてエンジン２の自着火に係る状態量であり、データには実測値（筒内圧、当量比、オクタン価、H2O濃度など）と、実測値から算出される予測値（未燃ガス温度、燃焼割合、残留ガス割合など）とが含まれる。これらのデータは、ある時点での瞬間値ではなく、複数サイクルにわたる履歴データである。データのうち実測値はセンサ１６によってエンジン２から取得され、予測値はコンピュータ１２によって算出される。また、データ取得のための実機試験は様々な運転条件下で行われ、コンピュータ１２は取得データに基づいて運転条件毎にノック発生点火時期を計算する。コンピュータ１２は、実機試験の運転条件を設定してドライバ１４に指示する機能を有している。なお、ここで言う運転条件には、エンジン２の点火時期も含まれている。この場合、コンピュータ１２は第１の発明にかかる「運転条件設定手段」に相当する。また、コンピュータ１２、ドライバ１４及びセンサ１６によって第１の発明にかかる「状態量測定手段」が実現される。
【００１６】
図２は、本実施の形態にかかるノック発生点火時期の推定方法を説明するためのフローチャートである。ノック発生点火時期推定装置１０が図２のフローチャートに示すルーチンに従って動作することで、エンジン２のノック発生点火時期が推定されるようになっている。
【００１７】
図２に示すルーチンの最初のステップＳ００１では、運転条件のうち点火時期のみを異ならせた複数点において実機試験が行われ、各点火時期のもとで必要なデータが取得される。取得したデータは、次のステップＳ００２においてノック発生を判定するための判定値Cl、Cbの算出に用いられる。
【００１８】
判定値Cl、Cbの算出には公知のノックモデルが用いられる。ノックモデルはコンピュータ１２にプログラムとして記憶されている。ノックモデルには、自着火予測モデルと、ノック判定値算出モデルとが含まれる。
【００１９】
自着火予測モデルは、自着火の条件が成立するクランク角を予測するモデルである。本実施の形態では、自着火予測モデルとしてLivengood-Wuモデルを使用する。Livengood-Wuモデルは、圧力と温度の関数として表される自着火遅れ時間の逆数の積分式である。Livengood-Wuモデルによれば、筒内圧及び未燃ガス温度を入力値として式の積分を実行し、その積分値が１になるときのクランク角を求めることで、自着火が起きるクランク角を予測することができる。筒内圧及び未燃ガス温度以外のデータは、Livengood-Wuモデルの各種パラメータの値を決定するのに用いられる。
【００２０】
ステップＳ００２では、ステップＳ００１で取得したデータを用いてLivengood-Wuモデルによる計算を実行し、自着火が起きるクランク角（自着火クランク角）を予測する。この自着火クランク角が第１判定値Clとなる。第１判定値Clは、実機試験が行われた複数の点火時期のそれぞれについて算出される。
【００２１】
ノック判定値算出モデルは、ノックの発生に必要な量の未燃燃料が存在する限界クランク角を算出するモデルである。本実施の形態では、ノック判定値算出モデルとしてWiebeモデルを使用する。Wiebeモデルによれば、クランク角と燃焼割合（供給燃料の総発熱量に対するそのクランク角までの総発熱量）との関係を算出することができる。燃焼割合が高くなるほど筒内の未燃燃料は減少していくため、燃焼割合があまりに高い状態では、もはや燃えるべき燃料が無いために急激な燃焼は起こらずノックは発生しない。つまり、燃焼割合にはノックが発生しうる上限値が存在する。Wiebeモデルを用いることで、燃焼割合が上限値に達するクランク角、すなわち、ノックの発生に必要な量の未燃燃料が存在する限界クランク角を算出することができる。
【００２２】
ステップＳ００２では、Wiebeモデルの各パラメータの値をステップＳ００１で取得したデータに基づいてマップから決定し、Wiebeモデルによる計算によって限界クランク角を算出する。この限界クランク角が第２判定値Cbとなる。第２判定値Cbは、実機試験が行われた複数の点火時期のそれぞれについて算出される。
【００２３】
ノックが発生する条件は、自着火発生タイミングにおいてノックの発生に必要な量の未燃燃料が残っていることである。点火時期がＭＢＴよりも大きく遅角されている場合、自着火クランク角は限界クランク角よりも遅角側にあるためノックは発生しない。しかし、点火時期が進角されるに従って自着火クランク角は限界クランク角に近付き、やがて自着火クランク角が限界クランク角に一致してノックが発生する。そのときの点火時期、つまり、第１判定値Clが第２判定値Cbに一致したときの点火時期の設定値が、ノック発生点火時期であると推定することができる。
【００２４】
本実施の形態では、ステップＳ００２で算出した第１判定値Cl及び第２判定値Cbを用い、次のようにしてノック発生点火時期を推定する。図３は、点火時期と第１判定値Clとの関係、及び、点火時期と第２判定値Cbとの関係を示す図である。図の横軸は点火時期であり、縦軸はクランク角である。この図に示すように、第１判定値Clと点火時期との関係は直線L1で表すことができる。また、第２判定値Cbと点火時期との関係も直線L2で表すことができる。図に示すような一定の線形関係が各判定値Cl、Cbと点火時期との間に存在することは、実験によって確認されている。この事実を利用すれば、第１判定値Clと点火時期との関係を示す直線L1と、第２判定値Cbと点火時期との関係を示す直線L2との交点を計算することで、ノック発生点火時期を求めることができる。
【００２５】
各直線L1、L2は、図４に示すように、それぞれ少なくとも２点のデータがあれば近似することができる。ステップＳ００３では、ステップＳ００２で算出した複数点の第１判定値Clから直線L1を近似計算し、また、複数点の第２判定値Cbから直線L2を近似計算する。そして、２つの直線L1、L2が交差する点火時期をノック発生点火時期として算出する。
【００２６】
このように、本実施の形態にかかるノック発生点火時期の推定方法によれば、エンジン２の実機試験で得られたデータ（エンジン２の自着火に係る状態量）に基づいてノック発生点火時期を求めることができる。しかも、各判定値Cl、Cbと点火時期との関係を示す２本の直線L1、L2を計算することで、ノックが発生しない運転条件でのデータからでも、ノック発生点火時期を求めることができる。直線L1、L2は少なくとも２つの点火時期での実機データがあれば近似計算できるので、実機データの取得に多くの工数や時間をかける必要もない。
【００２７】
ただし、直線L1、L2の近似計算に用いるデータの数が少ないと、その分、直線L1、L2の近似精度が低下してノック発生点火時期の予測誤差は大きくなる。そこで、本実施の形態にかかるノック発生点火時期の推定方法では、ノック発生点火時期の確定値は実機による適合試験によって得るものとし、直線L1、L2の交点から算出したノック発生点火時期（以下、予測点火時期という）は、実機試験によってノック発生点火時期を探索するための基準値として用いる。以下に説明するステップＳ００４乃至Ｓ０１１の処理が実機による適合試験に該当する。
【００２８】
ステップＳ００４では、実機試験における点火時期の設定値を予測点火時期に対して遅角するか否か決定する。予測点火時期の精度が高いと推定される場合には、点火時期の遅角は行わない。具体的には、次のような条件を満たしている場合である。
条件１：直線L1、L2の近似に使用したデータの点火時期方向の間隔が大きい。
条件２：直線L1、L2の近似に使用したデータの点数が多い。
条件３：筒内圧などのノックモデル入力値の測定誤差が小さい。
【００２９】
ステップＳ００４で点火時期を遅角すると決定した場合は、ステップＳ００５の処理を行う。ステップＳ００５では、ノックモデルの精度に応じて点火時期の設定値を予測点火時期に対して遅角する。ここでは、直線L1、L2の近似に使用したデータの点火時期方向の間隔が小さいほど、直線L1、L2の近似に使用したデータの点数が少ないほど、また、ノックモデル入力値の測定誤差が大きいほど、ノックモデルの精度は低いと考える。そして、ノックモデルの精度が低いほど、予測点火時期に対する遅角幅を大きく設定する。
【００３０】
ステップＳ００６では、ステップＳ００４及びＳ００５の処理によって決定したとおりに点火時期を設定し、実機試験を実施する。この実機試験は、ステップＳ００１で実施した実機試験と同一の運転条件下で行われる。実機試験ではノックセンサを用いてノックの発生を検出する。実機試験の終了後、ステップＳ００７の判定を行う。
【００３１】
ステップＳ００７では、実機試験においてノックが発生したか否か判定する。ノックが発生していなかった場合、ステップＳ００８に進む。ステップＳ００８では、点火時期の設定値が現在値よりも所定幅だけ進角される。このときの進角幅は、ステップＳ００５で設定される遅角幅よりも格段に小さい幅に設定されている。
【００３２】
ステップＳ００７の処理後は、再びステップＳ００６に進んで実機試験を実施する。そして、ステップＳ００７でノックが発生したと判定されるまで、ステップＳ００６乃至Ｓ００８の処理を繰り返し実施する。ステップＳ００７でノックが発生したと判定されたとき、そのときの点火時期の設定値をノック発生点火時期として確定する。このように、点火時期を設定し直す毎にノックが発生するか否か実機試験によって検証することで、ノック発生点火時期をより精度よく推定することができる。
【００３３】
ステップＳ００４乃至Ｓ００８の処理では、ノック発生点火時期は図５に示すように探索される。この図に示すように、本実施の形態にかかるノック発生点火時期の推定方法によれば、ノック発生点火時期の探索範囲は、ノックモデルを用いて計算した予測点火時期を基準として決定されるので、従来の実機試験のみに頼った方法に比較して、ノック発生点火時期の探索に要する工数や時間を大幅に抑えることができる。
【００３４】
ステップＳ００９では、全ての試験が完了したか否か判定する。全ての試験が完了している場合には、本ルーチンを終了する。試験項目が残っている場合、つまり、ノック発生点火時期の推定作業が必要な運転条件（試験条件）が残っている場合は、ステップＳ０１０に進み、次の試験条件に切り替えてノック発生点火時期の推定作業を続行する。
【００３５】
また、ステップＳ０１１では、ステップＳ００６の実機試験で取得したデータをノックモデルの適合点として追加する。ノックモデルで使用するパラメータは、例えば図６に示すように、エンジンの運転状態（図６ではエンジン回転数）に応じて決定される。ノックモデルの計算に使用するパラメータ値を決定する場合、適合点が存在するエンジン回転数については適合点の値がパラメータ値として用いられる。しかし、適合点が存在しないエンジン回転数については、近傍の適合点から補間計算せざるを得ない。実機試験で得られたデータを適合点として追加していくことにすれば、パラメータ値を適合点から補間計算する場合の補間精度を高めることができ、ひいては、ノックモデルを用いたノック発生点火時期の推定精度を高めることができる。ステップＳ０１１の処理後は再びステップ００１の処理から開始する。
【００３６】
本実施の形態では、ノック発生点火時期推定装置１０によってステップＳ００２の処理が実行されることにより、第１の発明の「第１判定値算出手段」及び「第２判定値算出手段」が実現される。また、ノック発生点火時期推定装置１０によってステップＳ００３の処理が実行されることにより、第１の発明の「ノック発生点火時期算出手段」が実現される。さらに、ノック発生点火時期推定装置１０によってステップＳ００４及びＳ００５の処理が実行されることにより、第２の発明の「点火時期初期設定手段」が実現される。また、ノック発生点火時期推定装置１０によってステップＳ００６乃至Ｓ００８の処理が実行されることにより、第２の発明の「ノック検証手段」が実現され、ステップＳ００７の判定でＮＯルートが選択されることにより、第２の発明の「ノック発生点火時期決定手段」が実現される。
【００３７】
以上、本発明の実施の形態について説明したが、本発明は上記実施の形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲で種々変形して実施することができる。例えば、次のように変形して実施してもよい。
【００３８】
上記実施の形態では、ノックモデルを用いたノック発生点火時期の推定後、すなわち、ステップＳ００１乃至Ｓ００３の処理の後、実機による適合試験、すなわち、ステップＳ００４乃至Ｓ００８の処理によってノック発生点火時期の検証を行っている。しかし、ノックモデルの精度が十分に高い場合には、実機による適合試験は省略し、ノックモデルを用いて計算したノック発生点火時期を確定値としてもよい。
【００３９】
また、図１ではノック発生点火時期推定装置１０を１つのユニットとして表しているが、複数の装置からなるシステムとして構成してもよい。また、図２のフローチャートに従って実施するノック発生点火時期の推定作業は、完全に自動化するのではなく、その一部の処理を人間が実施してもよい。
【図面の簡単な説明】
【００４０】
【図１】本発明の実施の形態としてのノック発生点火時期推定装置の構成を示すブロック図である。
【図２】本発明の実施の形態にかかるノック発生点火時期の推定方法を説明するためのフローチャートである。
【図３】点火時期と第１判定値Clとの関係、及び、点火時期と第２判定値Cbとの関係を示す図である。
【図４】判定値Cl、Cbからノック発生点火時期を推定する方法について説明するための図である。
【図５】実機試験によるノック発生点火時期の探索手順について示す図である。
【図６】ノックモデルの適合点について説明するための図である。
【符号の説明】
【００４１】
２  エンジン
１０  ノック発生点火時期推定装置
１２  コンピュータ
１４  ドライバ
１６  センサ
【図１】
2008095553.tif 000003
【図２】
2008095553.tif 000004
【図３】
2008095553.tif 000005
【図４】
2008095553.tif 000006
【図５】
2008095553.tif 000007
【図６】
2008095553.tif 000008