JP6683283B1 - 内燃機関の点火時期制御装置 - Google Patents

Abstract

【課題】点火時期の過剰遅角を防止する。【解決手段】燃焼室（５）内の圧力を検出する圧力センサ（１９）の出力値から求められるノッキング強度を代表する値を、機関本体（１）の振動を検出するノッキングセンサ（１８）の出力値から推定する第１の学習済みニューラルネットワークと、点火時期が遅角されたときのノッキング強度を代表する値の推定値の予測値又は予測低下量を推定する第２の学習済みニューラルネットワークとが記憶されている。点火時期が遅角されたときのノッキング強度を代表する値の推定値の予測値とノッキング強度を代表する値の推定値との差が小さくかつノッキング強度を代表する値の推定値が閾値を越えているときに、次のサイクルにおける点火時期の遅角制御が行われる。【選択図】図２４

Description

本発明は内燃機関の点火時期制御装置に関する。

内燃機関では、燃焼室内において、混合気が点火された後、末端ガスが自着火すると圧力波が発生し、ノッキングが生ずる。ノッキングが生ずると機関本体が振動し、この場合、ノッキング強度が高くなるほど機関本体の振動が激しくなる。そこで、機関本体に機関本体の振動を検出するノッキングセンサを取り付け、ノッキングセンサにより検出された機関本体の振動の激しさからノッキング強度を検出し、ノッキング強度が閾値を越えたときにはノッキングが発生していると判断して点火時期を遅角するようにした内燃機関が公知である（例えば特許文献１を参照）。しかしながら機関本体は、例えば、吸気弁および排気弁の着座作用や、燃料噴射弁のニードルの着座作用等の機械的動作によっても振動し、従って、機関本体の振動を検出するようにしたノッキングセンサを用いると、このような機械的動作が原因で機関本体の振動強度が激しくなったときに、ノッキング強度が高くないにも拘わらず、ノッキング強度が高いと誤判定してしまうという問題がある。

特開２００６−２２６９６７号公報

一方、ノッキングが生じたときには、ノッキング強度が高くなるほど燃焼室内の圧力の変動が激しくなる。従って、点火により発生する混合気の燃焼圧を検出可能な圧力センサを用いると、この圧力センサの出力値の変動の激しさからノッキング強度を検出できることになる。この場合、例えば、ノッキング強度が高くなると、圧力センサの出力値のピーク値が高くなり、従って、この圧力センサの出力値のピーク値は、ノッキング強度を代表する値の一つであると言える。従って、このようなノッキング強度を代表する値から、ノッキング強度を検出することができることになる。この場合、機械的動作に基づく機関本体の振動は、圧力センサの出力値に影響を与えず、従って、ノッキング強度を代表する値にも影響を与えないので、圧力センサを用いると、ノッキング強度を精度よく検出することができることになる。

しかしながら、この圧力センサは極めて高価である。しかも圧力センサは長期間使用されると、圧力センサにデポジットが徐々に付着し、この付着したデポジットによって燃焼室内における混合気の燃焼形態が変化してしまう。従って、市販車に対して、このような圧力センサを用いることは難しい。そこで、本発明では、圧力センサの出力値から求められるノッキング強度を代表する値を、ニューラルネットワークを用いて、ノッキングセンサの出力値から推定するようにしている。即ち、ノッキングセンサの出力値が入力されると、ノッキング強度を代表する推定値が出力されるように、ニューラルネットワークの重みを学習し、重みの学習が完了した学習済みニューラルネットワークを用いて、ノッキングセンサの出力値から、ノッキング強度を代表する値を推定するようにしている。

ところが、市販車の制御装置にこのような学習済みニューラルネットワークを格納し、格納された学習済みニューラルネットワークを用いて、ノッキングセンサの出力値から、ノッキング強度を代表する値を推定すると問題を生ずる。即ち、機関の構成部品には公差があり、従って、市販車ごとに機関の構成部品の寸法がばらつくために、市販車毎に異なる機関の振動が発生する。しかしながら、学習済みニューラルネットワークは、このように市販車毎に発生する異なる機関の振動については、重みの学習がなされておらず、従って、このように重みの学習がなされていない機関の振動、即ち、未学習の機関の振動が発生すると、学習済みニューラルネットワークにおいて、ノッキングにより機関本体に振動が発生していると誤判断される場合がある。

このような誤判断が行われると、点火時期が過剰に遅角される危険性がある。即ち、ノッキングが発生していることが原因でノッキング強度を代表する値が閾値を越え、それにより点火時期が遅角されたときには、点火時期の遅角作用により燃焼圧が低下するためにノッキング強度を代表する値が低下する。その結果、ノッキング強度を代表する値が閾値以下となり、点火時期の遅角作用が停止される。一方、未学習の機関の振動が発生し、それによりノッキング強度を代表する値が閾値を越えた場合でも、点火時期が遅角される。しかしながら、このとき、点火時期の影響を受けない未学習の機関の振動が発生していた場合には、点火時期が遅角されても、未学習の機関の振動は抑制されない。従って、このとき、点火時期が遅角された後も点火時期の影響を受けない未学習の機関の振動が発生し続けた場合には、ノッキング強度を代表する値が閾値を越え続けることになり、点火時期の遅角作用が続行し続けることになる。その結果、点火時期が過剰に遅角し、それにより機関の出力が大幅に低下してしまうという問題を生ずることになる。

このような問題を解決するために、本発明によれば、点火により発生する混合気の燃焼圧を検出するための圧力センサの出力値から求められるノッキング強度を代表する値を、機関本体の振動を検出するノッキングセンサの出力値から推定する第１の学習済みニューラルネットワークが記憶装置内に記憶されており、点火時期が遅角されたときのノッキング強度を代表する値の推定値の予測値又は予測低下量を推定する第２の学習済みニューラルネットワークが記憶装置内に記憶されており、第１の学習済みニューラルネットワークを用いて算出されたノッキング強度を代表する値の推定値が予め定められた閾値を越えると次のサイクルにおける点火時期が遅角され、点火時期が遅角された該次のサイクルでは、第２の学習済みニューラルネットワークを用いて算出されたノッキング強度を代表する値の推定値の予測値と第１の学習済みニューラルネットワークを用いて算出されたノッキング強度を代表する値の推定値との差に基づいて更に次のサイクルにおける点火時期の遅角制御が行われ、この遅角制御では、上記差が予め定められた設定値よりも小さいときには、ノッキング強度を代表する値の推定値が予め設定された閾値よりも大きければ、更に次のサイクルにおける点火時期の遅角作用が行われ、上記差が予め定められた設定値よりも大きいときには、たとえノッキング強度を代表する値の推定値が予め設定された閾値よりも大きかったとしても、更に次のサイクルにおける点火時期の遅角作用が行われない内燃機関の点火時期制御装置が提供される。

未学習の機関の振動の発生により、ノッキング強度を代表する値の推定値が閾値を越えて次のサイクルにおける点火時期が遅角されたときには、上述の差が予め定められた設定値よりも大きくなり、従って、このとき、たとえノッキング強度を代表する値の推定値が予め設定された閾値よりも大きかったとしても、更に次のサイクルにおける点火時期の遅角作用が行われない。従って、ノッキング強度を代表する値を、圧力センサを用いることなく、高い精度で取得しつつ、点火時期が過剰に遅角されるのを阻止することが可能となる。

図１は内燃機関の全体図である。 図２は図１に示す内燃機関の側面断面図である。 図３はニューラルネットワークの説明図である。 図４は、シリンダヘッドの底面図である。 図５は内燃機関の全体図である。 図６Ａ、図６Ｂおよび図６Ｃは、ノッキングセンサの出力値を示す図である。 図７Ａ、図７Ｂおよび図７Ｃは、圧力センサの出力値を示す図である。 図８は、本発明による実施例において用いられている第１のニューラルネットワークを示す図である。 図９は、訓練データセットを示す図である。 図１０は、学習処理ルーチンを示すフローチャートである。 図１１は、電子制御ユニットにデータを読み込むためのルーチンを示すフローチャートである。 図１２は、ノッキング処理ルーチンを示すフローチャートである。 図１３は、ノッキング強度代表値の推定ルーチンを示すフローチャートである。 図１４は、ノッキング強度代表値の推定ルーチンの別の実施例を示すフローチャートである。 図１５は、ノッキング判定ルーチンを示すフローチャートである。 図１６は、閾値Ｍij を示す図である。 図１７は、点火制御ルーチンを示すフローチャートである。 図１８は、基準点火時期のマップを示す図である。 図１９は、ノッキング処理の概要を示す図である。 図２０Ａおよび図２０Ｂは、ノッキング発生時の燃焼圧を説明するための図である。 図２１は、ノッキング判定を説明するための図である。 図２２Ａおよび図２２Ｂは、点火遅角量又は点火遅角速度と燃焼圧の上昇割合等との関係を示す図である。 図２３は、ノッキング判定を説明するための図である。 図２４Ａおよび図２４Ｂは、推定値ｙｅの予測値又は予測低下量を説明するための図である。 図２５Ａおよび図２５Ｂは、推定値ｙｅ予測値又は予測低下量を説明するための図である。 図２６は、入力パラメータの一覧表を示す図である。 図２７は、本発明による実施例において用いられている第２のニューラルネットワークを示す図である。 図２８は、データの一覧表を示す図である。 図２９は、訓練データセットを示す図である。 図３０は、ノッキング処理ルーチンを示すフローチャートである。 図３１は、ノッキング判定ルーチンを示すフローチャートである。 図３２は、ノッキング判定ルーチンを示すフローチャートである。 図３３は、点火制御ルーチンを示すフローチャートである。 図３４はリカレントニューラルネットワークの説明図である。 図３５はリカレントニューラルネットワークの説明図である。 図３６は、入力値および出力値の一覧表を示す図である。 図３７は、訓練データセットを示す図である。 図３８はリカレントニューラルネットワークの説明図である。 図３９は、ノッキング処理ルーチンを示すフローチャートである。 図４０は、ノッキング判定ルーチンを示すフローチャートである。 図４１は、ノッキング判定ルーチンを示すフローチャートである。 図４２は、点火制御ルーチンを示すフローチャートである。 図４３は、ゲート区間と機関回転数との関係を示す図である。 図４４は、ゲート区間と分割区間を示す図である。 図４５Ａ，図４５Ｂおよび図４５Ｃは、運転領域を示す図である。 図４６は内燃機関の全体図である。 図４７は、本発明による別の実施例において用いられている第１のニューラルネットワークを示す図である。 図４８は、訓練データセットを示す図である。

＜内燃機関の全体構成＞
図１および図２に内燃機関の全体図を示す。図２を参照すると、１は機関本体、２はシリンダブロック、３はシリンダヘッド、４はシリンダブロック２内で往復動するピストン、５は燃焼室、６は吸気弁、７は吸気ポート、８は排気弁、９は排気ポート、１０は各燃焼室５内に燃料、例えば、ガソリンを供給するための燃料噴射弁、１１は各燃焼室５内に配置された点火栓を夫々示す。図１および図２を参照すると、吸気ポート７は夫々対応する吸気枝管１２を介してサージタンク１３に連結され、サージタンク１３は吸気ダクト１４および吸入空気量検出器１５を介してエアクリーナ１６に連結される。

一方、排気ポート９は排気マニホルド１７に連結され、排気マニホルド１７は、排気マニホルド１７内の排気ガスをサージタンク１３内に再循環させるための排気ガス再循環（以下、ＥＧＲと称す）通路２７を介してサージタンク１３に連結される。このＥＧＲ通路２７内にはＥＧＲ制御弁２８が配置される。図１に示される実施例では、機関の運転状態に応じたＥＧＲ率（＝再循環排気ガス量/（再循環排気ガス量＋吸入空気量））が予め設定されており、ＥＧＲ制御弁２８は、ＥＧＲ率が予め設定されたＥＧＲ率となるように制御される。

図１および図２に示されるように、図１に示される実施例では、シリンダブロック２に、シリンダブロック２の振動、即ち、機関本体１の振動を検出するためのノッキングセンサ１８が取り付られている。図１に示される例では、このノッキングセンサ１８では、振動の検出素子として圧電素子が用いられており、このノッキングセンサ１８は機関本体１の振動に比例した出力電圧を発生する。なお、ノッキングが発生すると、機関本体１には、周波数が５KHzから２５KHz程度の振動が発生し、このときノッキングセンサ１８の出力電圧は、周波数が５KHzから２５KHz程度で変動する。従って、ノッキングセンサ１８の出力電圧の変動、即ち、ノッキングセンサ１８の出力値の変動から、ノッキングの発生を検出することができることになる。

一方、図１において３０は、機関の運転を制御するための電子制御ユニットを示している。図１に示されるように、電子制御ユニット３０はデジタルコンピュータからなり、双方向性バス３１によって互いに接続された記憶装置３２、即ち、メモリ３２と、ＣＰＵ（マイクロプロセッサ）３３と、入力ポート３４および出力ポート３５を具備する。入力ポート３４には、ノッキングセンサ１８の出力信号が、対応するＡＤ変換器３６および５KHzから２５KHz程度の周波数の入力信号のみを通過させるデジタルバンドパスフィルタ３７を介して入力される。また、入力ポート３４には、吸入空気量検出器１５の出力信号が、対応するＡＤ変換器３６を介して入力される。また、アクセルペダル４０にはアクセルペダル４０の踏込み量に比例した出力電圧を発生する負荷センサ４１が接続され、負荷センサ４１の出力電圧は対応するＡＤ変換器３６を介して入力ポート３４に入力される。更に入力ポート３４にはクランクシャフトが例えば３０°回転する毎に出力パルスを発生するクランク角センサ４２が接続される。ＣＰＵ３３内ではクランク角センサ４２の出力信号に基づいて機関回転数が算出される。一方、出力ポート３５は対応する駆動回路３８を介して各気筒の燃料噴射弁１０、点火栓１１およびＥＧＲ制御弁２８に接続される。
＜ニューラルネットワークの概要＞

本発明による実施例では、ニューラルネットワークを用いてノッキング強度を代表する値を推定するようにしている。そこで最初にニューラルネットワークについて簡単に説明する。図３は簡単なニューラルネットワークを示している。図３における丸印は人工ニューロンを表しており、ニューラルネットワークにおいては、この人工ニューロンは、通常、ノード又はユニットと称される（本願では、ノードと称す）。図３においてL＝１は入力層、L＝２および L＝３は隠れ層、L＝４は出力層を夫々示している。また、図３において、ｘ_１およびｘ_２ は入力層 ( L＝１) の各ノードからの出力値を示しており、ｙ_１ およびｙ_２ は出力層 ( L＝４) の各ノードからの出力値を示しており、ｚ^（２） _１、ｚ^（２） _２ およびｚ^（２） _３ は隠れ層 ( L＝２) の各ノードからの出力値を示しており、ｚ^（３） _１、ｚ^（３） _２ およびｚ^（３） _３は隠れ層 ( L＝３) の各ノードからの出力値を示している。なお、隠れ層の層数は、１個又は任意の個数とすることができ、入力層のノードの数および隠れ層のノードの数も任意の個数とすることができる。また、出力層のノードの数は１個とすることもできるし、複数個とすることもできる。

入力層の各ノードでは入力がそのまま出力される。一方、隠れ層 ( L＝２) の各ノードには、入力層の各ノードの出力値ｘ_１およびｘ_２ が入力され、隠れ層 ( L＝２) の各ノードでは、夫々対応する重みｗおよびバイアスｂを用いて総入力値u が算出される。例えば、図３において隠れ層 ( L＝２) のｚ^（２） _ｋ（ｋ＝１，２，３）で示されるノードにおいて算出される総入力値u_ｋは、次式のようになる。

次いで、この総入力値u_ｋは活性化関数ｆにより変換され、隠れ層 ( L＝２) のｚ^（２） _ｋで示されるノードから、出力値ｚ^（２） _ｋ(= f (ｕ_ｋ)) として出力される。一方、隠れ層 ( L＝３) の各ノード には、隠れ層 ( L＝２) の各ノードの出力値ｚ^（２） _１、ｚ^（２） _２ およびｚ^（２） _３が入力され、隠れ層 ( L＝3 ) の各ノードでは、夫々対応する重みｗおよびバイアスｂを用いて総入力値ｕ（Σｚ・ｗ＋ｂ）が算出される。この総入力値ｕは同様に活性化関数により変換され、隠れ層 ( L＝3 ) の各ノードから、出力値ｚ^（３） _１、ｚ^（３） _２ およびｚ^（３） _３として出力される、この活性化関数としては、例えば、シグモイド関数σが用いられる。

一方、出力層 ( L＝４) の各ノード には、隠れ層 ( L＝３) の各ノードの出力値ｚ^（３） _１、ｚ^（３） _２ およびｚ^（３） _３が入力され、出力層 の各ノードでは、夫々対応する重みｗおよびバイアスｂを用いて総入力値ｕ（Σｚ・ｗ＋ｂ）が算出されるか、又は、夫々対応する重みｗのみを用いて総入力値ｕ（Σｚ・ｗ）が算出される。本発明による実施例では、出力層のノードでは、活性化関数として恒等関数が用いられており、従って、出力層のノードからは、出力層のノードにおいて算出された総入力値ｕが、そのまま出力値ｙとして出力される。
＜ニューラルネットワークにおける学習＞

さて、ニューラルネットワークの出力値ｙの正解値を示す教師データをｙ_ｔとすると、ニューラルネットワークにおける各重みｗおよびバイアスｂは、出力値ｙと教師データをｙ_ｔとの差が小さくなるように、誤差逆伝播法を用いて学習される。この誤差逆伝播法は周知であり、従って、誤差逆伝播法についてはその概要を以下に簡単に説明する。なお、バイアスｂは重みｗの一種なので、以下、バイアスｂも含めて重みｗと称する。さて、図３に示すようなニューラルネットワークにおいて、Ｌ＝２，Ｌ＝３又は Ｌ＝４の各層のノードへの入力値ｕ^（Ｌ）における重みをｗ^（Ｌ）で表すと、誤差関数Ｅの重みｗ^（Ｌ）による微分、即ち、勾配∂Ｅ/∂ｗ^（Ｌ）は、書き換えると、次式で示されるようになる。

ここで、ｚ^{（Ｌ−１）}・∂ｗ^（Ｌ）＝ ∂ｕ^（Ｌ）であるので、（∂Ｅ/∂ｕ^（Ｌ））＝δ^（Ｌ）とすると、上記（１）式は、次式でもって表すことができる。

ここで、ｕ^（Ｌ）が変動すると、次の層の総入力値ｕ^{（Ｌ＋１）}の変化を通じて誤差関数Ｅの変動を引き起こすので、δ^（Ｌ）は、次式で表すことができる。

ここで、z^（Ｌ）＝ｆ(u^（Ｌ）) と表すと、上記（３）式の右辺に現れる入力値ｕ_k ^{（Ｌ＋１）}は、次式で表すことができる。

ここで、上記（３）式の右辺第１項（∂Ｅ/∂ｕ^{（Ｌ＋１）}）はδ^{（Ｌ＋１）}であり、上記（３）式の右辺第２項（∂u_ｋ ^{（Ｌ＋１）} /∂u^（Ｌ））は、次式で表すことができる。

従って、δ^（Ｌ）は、次式で示される。

即ち、δ^{（Ｌ＋１）}が求まると、δ^（Ｌ）を求めることができることになる。

さて、出力層 ( L＝４) のノードが一個であって、或る入力値に対して教師データｙ_ｔが求められており、この入力値に対する出力層からの出力値がｙであった場合において、誤差関数として二乗誤差が用いられている場合には、二乗誤差Ｅは、Ｅ＝１/２(ｙ−ｙ_ｔ)^２で求められる。この場合、出力層（Ｌ＝４）のノードでは、出力値ｙ＝ ｆ(u^（Ｌ）) となり、従って、この場合には、出力層（Ｌ＝４）のノードにおけるδ^（Ｌ）の値は、次式で示されるようになる。

この場合、本発明による実施例では、前述したように、ｆ(u^（Ｌ）) は恒等関数であり、ｆ’(u^（Ｌｌ）) ＝ １となる。従って、δ^（Ｌ）＝ｙ−ｙ_ｔ となり、δ^（Ｌ）が求まる。

δ^（Ｌ）が求まると、上式（６）を用いて前層のδ^{（Ｌ−１）}が求まる。このようにして順次、前層のδが求められ、これらδの値を用いて、上式（２）から、各重みｗについて誤差関数Ｅの微分、即ち、勾配∂Ｅ/∂ｗ^（Ｌ）か求められる。勾配∂Ｅ/∂ｗ^（Ｌ）か求められると、この勾配∂Ｅ/∂ｗ^（Ｌ）を用いて、誤差関数Ｅの値が減少するように、重みｗが更新される。即ち、重みｗの学習が行われることになる。なお、図３に示されるように、出力層 ( L＝４) が複数個のノードを有する場合には、各ノードからの出力値をｙ_１、ｙ_１・・・、対応する教師データｙ_ｔ１、ｙ_ｔ２・・・とすると、誤差関数Ｅとして、次の二乗和誤差Ｅが用いられる。

この場合も、出力層 ( L＝４) の各ノードにおけるδ^（Ｌ）の値は、δ^（Ｌ）＝ｙ−ｙ_tk （ｋ＝１，２・・・ｎ）となり、これらδ^（Ｌ）の値から上式（６）を用いて前層のδ^{（Ｌ−１）}が求まる。
＜本発明による実施例＞

さて、内燃機関では、ノッキング強度が高くなるほど機関本体１の振動が激しくなる。従って、図１および図２に示されるように、機関本体１にノッキングセンサ１８を取り付け、ノッキングセンサ１８により検出された機関本体の振動の激しさからノッキング強度を検出することができる。しかしながら前述したように、機関本体１は、例えば、吸気弁６および排気弁８の着座作用や、燃料噴射弁１０のニードルの着座作用等の機械的動作によっても振動し、従って、機関本体１の振動を検出するようにしたノッキングセンサ１８を用いると、このような機械的動作が原因で機関本体の振動強度が激しくなったときに、ノッキング強度が高くないにも拘わらず、ノッキング強度が高いと誤判定してしまうという問題がある。

一方、前述したように、点火により発生する混合気の燃焼圧を検出可能な圧力センサを用いると、この圧力センサの出力値の変動の激しさからノッキング強度を検出できる。この場合、例えば、ノッキング強度が高くなると、圧力センサの出力値のピーク値が高くなり、従って、この圧力センサの出力値のピーク値は、ノッキング強度を代表する値の一つであると言える。従って、このようなノッキング強度を代表する値から、ノッキング強度を検出することができることになる。この場合、機械的動作に基づく機関本体の振動は、圧力センサの出力値に影響を与えず、従って、ノッキング強度を代表する値にも影響を与えないので、圧力センサを用いると、ノッキング強度を精度よく検出することができることになる。

しかしながら、この圧力センサは高価であり、しかも圧力センサに徐々に付着するデポジットにより燃焼室５内における混合気の燃焼形態が変化してしまうので、市販車において、このような圧力センサを用いることは難しい。そこで、本発明では、圧力センサの出力値から求められるノッキング強度を代表する値を、ニューラルネットワークを用いて、ノッキングセンサの出力値から推定するようにしている。このことを説明するに当たり、最初に、図４、図５、図６Ａから図６Ｃおよび図７Ａから図７Ｃを参照しつつ、ノッキングセンサの出力値と、圧力センサの出力値並びに圧力センサの出力値から求められるノッキング強度を代表する値について説明する。

図４および図５は夫々、ニューラルネットワークの重みの学習に必要なデータを取得するためにだけ用いられるシリンダヘッド３の底面図および機関本体１の一部を示している。図４および図５を参照すると、圧力センサの出力値から求められるノッキング強度を代表する値を、ノッキングセンサ１８の出力値から推定するために、各気筒のシリンダヘッド３の内壁面上には、夫々圧力センサ１９が配置される。各圧力センサ１９からは、夫々対応する燃焼室５内の圧力に比例した出力電圧、即ち、点火により発生する混合気の燃焼圧を示す出力値が出力される。なお、この場合、燃焼室５内の圧力は、点火栓と一体型の圧力センサを用いて計測することもできる。これら圧力センサ１９は、ニューラルネットワークの重みの学習に必要なデータを取得するためだけに用いられる。

さて、ノッキングは通常、圧縮上死点から圧縮上死点後９０°までの間に発生し、ノッキングが発生したときには、周波数が５KHzから２５KHz程度の燃焼室５内の圧力変動および周波数が５KHzから２５KHz程度の機関本体１の振動が発生する。従って、ノッキングによる機関本体１の振動のみを抽出するために、通常、ノッキングセンサ１８の出力値は、ＡＤ変換された後、例えば５KHzから２５KHz程度の周波数の入力信号のみを通過させるデジタルバンドパスフィルタに送り込まれ、このデジタルバンドパスフィルタによりフィルタリングされた後のノッキングセンサ１８の出力値に基づいて、ノッキングの強度が判別される。

同様に、ノッキングによる燃焼室５内の圧力変動のみを抽出するために、通常、圧力センサ１９の出力値も、ＡＤ変換された後、例えば５KHzから２５KHz程度の周波数の入力信号のみを通過させるデジタルバンドパスフィルタに送り込まれ、このデジタルバンドパスフィルタによりフィルタリングされた後の圧力センサ１９の出力値に基づいて、ノッキングの強度が判別される。なお、このようにデジタルバンドパスフィルタによりフィルタリングされた後のノッキングセンサ１８の出力値を、以下、単にフィルタリング後のノッキングセンサ１８の出力値と称し、デジタルバンドパスフィルタによりフィルタリングされた後の圧力センサ１９の出力値を、以下、単にフィルタリング後の圧力センサ１９の出力値と称する。

図６Ａおよび図７Ａは夫々、同一のノッキングの発生に対するフィルタリング後のノッキングセンサ１８の出力値（Ｖ）の変化と、フィルタリング後の圧力センサ１９の出力値の変化（Ｖ）とを示している。なお、図６Ａの横軸および図７Ａの横軸は共に、ＡＴＤＣで表したクランク角（圧縮上死点から圧縮上死点後９０°）を示している。

図６Ａおよび図７Ａを比較するとわかるように、図７Ａに示されるフィルタリング後の圧力センサ１９の出力値の振動強度は、ノッキングが発生すると急激に増大し、その後、次第に減少していく。即ち、圧力センサ１９の出力値は、機械的動作に基づく機関本体の振動の影響を受けないために、ノッキングの発生が圧力センサ１９の出力値に明確に表れる。これに対し、図６Ａに示されるフィルタリング後のノッキングセンサ１８の出力値の振動強度も、ノッキングが発生すると増大はするが、ノッキングセンサ１８の出力値は機械的動作に基づく機関本体の振動の影響を大きく受けるため、ノッキングセンサ１８の出力値の振動強度は、ノッキングが発生した前後においても大きな値となっている。即ち、ノッキングセンサ１８の出力値には、機械的動作に基づく機関本体の振動がノイズとなって現れるため、ノッキングの発生が、ノッキングセンサ１８の出力値に明確に表れなくなる。

従って、ノッキングが発生したか否かは、圧力センサ１９の出力値から明確に判別できることになる。この場合、圧力センサ１９の出力値には、発生したノッキングの強度も明確に現れる。次に、このことについて、図７Ｂおよび図７Ｃを参照しつつ説明する。なお、図７Ｂおよび図７Ｃは、図７Ａに示されるフィルタリング後の圧力センサ１９の出力値の波形の一部を横軸方向に引き伸ばして図解的に表した図である。

さて、発生したノッキングの強度が高くなると、フィルタリング後の圧力センサ１９の出力値のピーク値が高くなる。従って、図７Ｂにおいて丸印で示されるフィルタリング後の圧力センサ１９の出力値のピーク値は、ノッキング強度を代表する値を示していることになる。一方、発生したノッキングの強度が高くなると、図７Ｃにおいて出力値の波形と出力値が０（Ｖ）の線とで囲まれたハッチングで示される面積の総和、即ち、フィルタリング後の圧力センサ１９の出力値の積分値（負側の積分値も正の値とする）が高くなる。従って、フィルタリング後の圧力センサ１９の出力値の積分値（負側の積分値も正の値とする）も、ノッキング強度を代表する値を示していることになる。このように圧力センサ１９の出力値からノッキング強度を代表する値を取得することができる。

一方、ノッキングセンサ１８の出力値からは、このようなノッキング強度を代表する値を抽出することは難しい。そこで本発明では、ニューラルネットワークを用い、ノッキングセンサ１８の出力値に基づいて、圧力センサ１９の出力値から求められるノッキング強度を代表する値を取得し得るようにニューラルネットワークの重みを学習し、それにより圧力センサ１９の出力値から求められるノッキング強度を代表する値を、学習済みのニューラルネットワークを用いてノッキングセンサ１８の出力値から推定するようにしている。そこで次に、ノッキングセンサ１８の出力値に基づいて、圧力センサ１９の出力値から求められるノッキング強度を代表する値を取得し得るようにニューラルネットワークの重みを学習する方法について説明する。

最初に、このとき用いられる第１のニューラルネットワークについて図８を参照しつつ説明する。図８を参照すると、この第１のニューラルネットワーク２０においても、図３に示されるニューラルネットワークと同様に、L＝１は入力層、L＝２および L＝３は隠れ層、L＝４は出力層を夫々示している。図８に示されるように、入力層 ( L＝１) がｎ個のノードからなり、ｎ個の入力値ｘ_１、ｘ_２・・・ｘ_ｎ−１、ｘ_ｎが、入力層 ( L＝１) の各ノードに入力されている。一方、図８には隠れ層 ( L＝２)および隠れ層 ( L＝３)が記載されているが、これら隠れ層の層数は、１個又は任意の個数とすることができ、またこれら隠れ層のノードの数も任意の個数とすることができる。なお、出力層 ( L＝４) のノードの数は１個とされており、出力層のノードからの出力値がｙで示されている。

次に、図８における入力値ｘ_１、ｘ_２・・・ｘ_ｎ−１、ｘ_ｎと、出力値ｙについて説明する。最初に、入力値ｘ_１、ｘ_２・・・ｘ_ｎ−１、ｘ_ｎについて図６Ｂおよび図６Ｃを参照しつつ説明する。なお、図６Ｂおよび図６Ｃは、図６Ａに示されるノッキングセンサ１８の出力値の波形の一部を横軸方向に引き伸ばして図解的に表した図である。本発明による実施例では、入力値ｘ_１、ｘ_２・・・ｘ_ｎ−１、ｘ_ｎとして、図６Ｂに示される値か、或いは、図６Ｃに示される値が用いられる。

即ち、図６Ｂに示される例では、フィルタリング後のノッキングセンサ１８の出力値そのものが入力値ｘ_１、ｘ_２・・・ｘ_ｎ−１、ｘ_ｎとされる。この場合、図６Ｂにおいて黒丸で示されるように、一定時間毎の、或いは、一定クランク角毎のノッキングセンサ１８の出力値が入力値ｘ_１、ｘ_２・・・ｘ_ｎ−１、ｘ_ｎとされる。一方、図６Ｃに示される例では、フィルタリング後のノッキングセンサ１８の出力値の積分値（負側の積分値も正の値とする）であって例えば一定クランク角度内におけるフィルタリング後のノッキングセンサ１８の出力値の積分値が入力値ｘ_１、ｘ_２・・・ｘ_ｎ−１、ｘ_ｎとされる。

一方、本発明の実施例では、圧力センサ１９の出力値から求められるノッキング強度を代表する値が、図８に示される出力値ｙとされる。この場合、図７Ｂにおいて丸印で示される圧力センサ１９の出力値のピーク値はノッキング強度を表しており、また、図７Ｃに示される圧力センサ１９の出力値の積分値（負側の積分値も正の値とする）もノッキング強度を表している。従って、本発明の実施例では、図７Ｂにおいて丸印で示される圧力センサ１９の出力値のピーク値が出力値ｙとされるか、或いは、図７Ｃに示される圧力センサ１９の出力値の積分値（負側の積分値も正の値とする）が出力値ｙとされる。この場合圧力センサ１９の出力値から求められるノッキング強度を代表する値の実測値が教師データｙｔとされる。

図９は、入力値ｘ_１、ｘ_２・・・ｘ_ｍ−１、ｘ_ｍと、入力値がｘ_１、ｘ_２・・・ｘ_ｍ−１、ｘ_ｍであるときに圧力センサ１９の出力値から求められるノッキング強度を代表する値の実測値、即ち、教師データｙｔとを用いて作成された訓練データセットを示している。図９に示されるように、この訓練データセットでは、入力値ｘ_１、ｘ_２・・・ｘ_ｎ−１、ｘ_ｎと教師データｙｔとの関係を表すｍ個のデータが取得されている。例えば、２番目のデータ（No. 2）には、取得された入力値ｘ_１２、ｘ_{２２・・・}ｘ_ｍ−１２、ｘ_ｍ２と教師データｙｔ_２とが列挙されており、ｍ−１番目のデータ（No. ｍ−１）には、取得された入力パラメータの入力値ｘ_１ｍ−１、ｘ_{２ｍ−１・・・}ｘ_{ｎ−１ｍ−１}、ｘ_ｎｍ−１と教師データｙｔ_ｍ−１が列挙されている。

入力値ｘ_１、ｘ_２・・・ｘ_ｎ−１、ｘ_ｎとして、フィルタリング後のノッキングセンサ１８の出力値そのものを用いる場合には、例えば、圧縮上死点から圧縮上死点後９０°までの間において一定クランク角毎に取得されたフィルタリング後のノッキングセンサ１８の出力値が、入力値ｘ_１、ｘ_２・・・ｘ_ｎ−１、ｘ_ｎとされる。この場合には、入力値の個数ｎが数百以上とされ、従って、図８に示される入力層 ( L＝１) のノードの個数ｎが数百以上とされる。

一方、本発明の実施例では、入力値ｘ_１、ｘ_２・・・ｘ_ｎ−１、ｘ_ｎとして、フィルタリング後のノッキングセンサ１８の出力値の積分値（負側の積分値も正の値とする）を用いる場合には、圧縮上死点から圧縮上死点後９０°までの間がクランク角度５°毎に１８個の区間に分割され、分割された各区間内におけるフィルタリング後のノッキングセンサ１８の出力値の積分値（負側の積分値も正の値とする）が、入力値ｘ_１、ｘ_２・・・ｘ_ｎ−１、ｘ_ｎとされる。この場合には、入力値の個数ｎは１８個となり、従って、図８に示される入力層 ( L＝１) のノードの個数ｎは１８個とされる。

一方、図９に示される訓練データセットにおける教師データｙｔは、圧力センサ１９の出力値から求められるノッキング強度を代表する値の実測値であり、この教師データｙｔとしては、図７Ｂにおいて丸印で示される圧力センサ１９の出力値のピーク値が用いられるか、或いは、図７Ｃに示される圧力センサ１９の出力値の積分値（負側の積分値も正の値とする）が用いられる。ところで、本発明による実施例では、図９に示される訓練データセットを用いて、図８に示される第１のニューラルネットワーク２０の重みの学習が行われる。そこで次に、図４および図５を参照しつつ、図９に示される訓練データセットの作成方法について説明する。

図４および図５には、図９に示される訓練データセットの作成方法の一例が示されている。図４および図５を参照すると、各気筒の燃焼室５に面するシリンダヘッド３の内壁面上には、ノッキング強度を代表する値を取得するために、夫々圧力センサ１９が配置されている。これら圧力センサ１９は、前述したように、学習に必要なデータを取得するためだけに設けられている。一方、図５に示される例では、ノッキングセンサ１８は、ノッキングセンサ１８の出力値の波形を検出可能な検出器２１、例えば、オシロスコープに接続され、各圧力センサ１９は、各圧力センサ１９の出力値の波形を検出可能な検出器２２、例えば、オシロスコープに接続されている。

検出器２１内では、ノッキングによる機関本体１の振動のみを抽出するために、ノッキングセンサ１８の出力値は、ＡＤ変換された後、例えば５KHzから２５KHz程度の周波数の入力信号のみを通過させるデジタルバンドパスフィルタに送り込まれ、このデジタルバンドパスフィルタによりフィルタリングされた後のノッキングセンサ１８の出力値の波形が検出される。一方、検出器２２では、ノッキングによる燃焼室５内の圧力変動のみを抽出するために、圧力センサ１９の出力値は、ＡＤ変換された後、例えば５KHzから２５KHz程度の周波数の入力信号のみを通過させるデジタルバンドパスフィルタに送り込まれ、このデジタルバンドパスフィルタによりフィルタリングされた後の圧力センサ１９の出力値の波形が検出される。

訓練データセットを作成する際には、機関負荷、機関回転数およびＥＧＲ率の種々の組み合わせについてノッキングの発生しない運転状態とノッキングの発生する運転状態の双方の状態が生ずるように機関を運転させ、そのとき各検出器２１、２２から得られたフィルタリング後の圧力センサ１９の出力値の波形データおよびフィルタリング後の圧力センサ１９の出力値の波形データに基づいて、図９に示されるような訓練データセットが作成される。この場合、この訓練データセットは、各検出器２１、２２から得られた波形データに基づき、手作業でもって作成することもできるし、各検出器２１、２２から得られた電子データに基づき、電子的に作成することもできる。このようにして作成された訓練データセットの電子データを用いて、図８に示される第１のニューラルネットワーク２０の重みの学習が行われる。

図５に示される例では、ニューラルネットワークの重みの学習を行うための学習装置２３が設けられている。図５に示されるように、この学習装置２３は、記憶装置２４、即ち、メモリ２４と、ＣＰＵ（マイクロプロセッサ）２５を具備している。図５に示される例では、図８に示される第１のニューラルネットワーク２０のノード数、および作成された訓練データセットの電子データが学習装置２３のメモリ２４に記憶され、ＣＰＵ２５においてニューラルネットワーク２０の重みの学習が行われる。

次に、入力値ｘ_１、ｘ_２・・・ｘ_ｎ−１、ｘ_ｎとして、上死点から上死点後９０°までの間においてクランク角度５°毎に分割された各区間内におけるフィルタリング後のノッキングセンサ１８の出力値の積分値（負側の積分値も正の値とする）を用い、教師データｙｔとして、圧力センサ１９の出力値のピーク値の実測値を用いて、学習装置２３により行われる第１のニューラルネットワーク２０の重みの学習方法について説明する。

図１０は、学習装置２３において行われる第１のニューラルネットワークの重みの学習処理ルーチンを示す。
図１０を参照すると、まず初めに、ステップ１００において、学習装置２３のメモリ２４に記憶されている訓練データセットの各データが読み込まれる。次いで、ステップ１０１において、第１のニューラルネットワーク２０の入力層 ( L＝１) のノード数、隠れ層 ( L＝２)および隠れ層 ( L＝３)のノード数および出力層 ( L＝４) のノード数が読み込まれ、これらノード数に基づいて、図８に示されるような第１のニューラルネットワーク２０が作成される。この場合、この例では、入力層 ( L＝１) のノード数は１８個とされ、出力層 ( L＝４) のノード数は１個とされる。

次いで、ステップ１０２では、第１のニューラルネットワーク２０の重みの学習が行われる。最初は、図９の１番目（No. 1）の入力値ｘ_１、ｘ_２・・・ｘ_ｎ−１、ｘ_ｎ、即ち、入力値ｘ_１、ｘ_２・・・ｘ₁₇、ｘ₁₈が第１のニューラルネットワーク２０の入力層 ( L＝１) の１８個の各ノードに入力される。次いで、このときの第１のニューラルネットワーク２０の出力値ｙと１番目（No. 1）の教師データｙｔ_１との間の二乗誤差Ｅ＝１/２(ｙ−ｙ_ｔ１)^２が算出され、この二乗誤差Ｅが小さくなるように、前述した誤差逆伝播法を用いて、第１のニューラルネットワーク２０の重みの学習が行われる。

図９の１番目（No. 1）のデータに基づく第１のニューラルネットワーク２０の重みの学習が完了すると、次に、図９の２番目（No. ２）のデータに基づく第１のニューラルネットワーク２０の重みの学習が、誤差逆伝播法を用いて行われる。同様にして、図９のｍ番目（No. ｍ）まで順次、第１のニューラルネットワーク２０の重みの学習が行われる。図９の１番目（No. 1）からｍ番目（No. ｍ）までの全てについて第１のニューラルネットワーク２０の重みの学習が完了すると、ステップ１０３に進む。

ステップ１０３では、例えば、図９の１番目（No. 1）からｍ番目（No. ｍ）までの全ての第１のニューラルネットワークの出力値ｙと教師データｙｔとの間の二乗和誤差Ｅが算出され、この二乗和誤差Ｅが、予め設定された設定誤差以下になったか否かが判別される。二乗和誤差Ｅが、予め設定された設定誤差以下になっていないと判別されたときには、ステップ１０２に戻り、再度、図９に示される訓練データセットに基づいて、第１のニューラルネットワークの重み学習が行われる。次いで、二乗和誤差Ｅが、予め設定された設定誤差以下になるまで、第１のニューラルネットワーク２０の重みの学習が続行される。ステップ１０３において、二乗和誤差Ｅが、予め設定された設定誤差以下になったと判別されたときには、ステップ１０４に進んで、第１のニューラルネットワーク２０の学習済み重みが学習装置２３のメモリ２４に記憶される。このようにしてノッキング強度を代表する値の推定モデルが作成される。

本発明による実施例では、このようにして作成されたノッキング強度を代表する値の推定モデルを用いて、機関においてノッキング処理が行われ、そのためにノッキング強度を代表する値の推定モデルが図１に示される電子制御ユニット３０に格納される。図１１は、ノッキング強度を代表する値の推定モデルを電子制御ユニット３０に格納するために、電子制御ユニット３０において行われる電子制御ユニットへのデータ読み込みルーチンを示している。

即ち、図１１に示されるように、ステップ１１０において、図８に示される第１のニューラルネットワーク２０の入力層 ( L＝１) のノード数、隠れ層 ( L＝２)および隠れ層 ( L＝３)のノード数および出力層 ( L＝４) のノード数が、電子制御ユニット３０のメモリ３２に読み込まれ、これらノード数に基づき、図８に示されるような第１のニューラルネットワーク２０が作成される。次いで、ステップ１１１において、この第１のニューラルネットワーク２０の学習済みの重みが電子制御ユニット３０のメモリ３２に読み込まれる。それにより、第１の学習済みニューラルネットワーク２０、即ち、ノッキング強度を代表する値の推定モデルが、機関の電子制御ユニット３０内に格納される。

さて、電子制御ユニット３０内に第１の学習済みニューラルネットワーク２０が作成されると、クランク角度５°毎に分割された各区間内におけるフィルタリング後のノッキングセンサ１８の出力値の積分値（負側の積分値も正の値とする）を第１の学習済みニューラルネットワーク２０に入力したときに、第１の学習済みニューラルネットワーク２０からは、圧力センサ１９の出力値から求められるノッキング強度を代表する値の推定値が出力される。このノッキング強度を代表する値の推定値は、ノッキング強度を代表する値の実測値に極めて近く、従って、ノッキングセンサ１８の出力値から、ノッキング強度を正確に推定できることになる。

このように、本発明による実施例では、機関本体１の振動を検出するノッキングセンサ１８と、内燃機関の燃焼室５内の圧力を検出する圧力センサ１９と、第１のニューラルネットワーク２０を用いて内燃機関のノッキング強度を代表する値を推定するための学習装置２３とが設けられている。圧力センサ１９の出力値からノッキング強度を代表する値が取得される。学習装置２３により、ノッキングセンサ１８により検出された機関本体１の振動を表す値を第１のニューラルネットワーク２０の入力値とすると共に、取得されたノッキング強度を代表する値を教師データとして、第１のニューラルネットワーク２０の重みが学習され、学習済みの第１のニューラルネットワーク２０を用いてノッキングセンサ１８の出力値からノッキング強度を代表する値が推定される。

また、本発明による実施例では、機関本体１の振動を検出するノッキングセンサ１８の出力値と、内燃機関の燃焼室５内の圧力を検出する圧力センサ１９の出力値とに基づき、ニューラルネットワーク２０を用いて内燃機関のノッキング強度が推定される。圧力センサ１９の出力値からノッキング強度を代表する値が求められ、ノッキングセンサ１８により検出された機関本体１の振動を表す値とノッキング強度を代表する値との関係を示す訓練データセットが作成される。ノッキングセンサ１８により検出された機関本体１の振動を表す値を第１のニューラルネットワーク２０の入力値とすると共に、取得されたノッキング強度を代表する値を教師データとして、第１のニューラルネットワーク２０の重みが学習される。学習済みの第１のニューラルネットワーク２０を用いてノッキングセンサ１８の出力値からノッキング強度を代表する値が推定される。

この場合、本発明による実施例では、ノッキング強度を代表する値が、予め設定されている期間内における圧力センサ１９の出力値のピーク値であり、或いは、ノッキング強度を代表する値が、予め設定されている期間内における圧力センサ１９の出力値の積分値である。この場合、予め設定されている期間が、一定のクランク角度範囲、例えば、圧縮上死点から圧縮上死点後９０°の範囲である。

一方、本発明による実施例では、機関本体１の振動を表す値が、予め設定されている期間内におけるノッキングセンサ１８の出力値であり、機関本体１の振動を表す値が、予め設定されている期間内において等分割された区間内におけるノッキングセンサ１８の出力値の積分値である。この場合、予め設定されている期間が、一定のクランク角度範囲、例えば、圧縮上死点から圧縮上死点後９０°の範囲である。また、ノッキングセンサ１８の出力値の積分値に関し、一例を挙げると、このノッキングセンサ１８の出力値の積分値は、クランク角度５°毎に等分割された各区間内におけるフィルタリング後のノッキングセンサ１８の出力値の積分値（負側の積分値も正の値とする）である。

次に、電子制御ユニット３０のメモリ３２内に形成された第１の学習済みニューラルネットワーク２０を用いて、機関運転時に実行されるノッキング処理について説明する。図１２は、この機関運転時に実行されるノッキング処理ルーチンを示している。図１２に示されるように、このノッキング処理ルーチンが実行されると、まず初めにステップ１２０において、圧力センサ１９の出力値から求められるノッキング強度を代表する値、即ち、ノッキング強度代表値の推定処理が行われる。次いで、ステップ１２１では、ノッキングが生じているか否かのノッキング判定処理が行われ、次いで、ステップ１２２では、点火制御が行われる。

図１３は、図１２のステップ１２０において行われるノッキング強度代表値の推定処理ルーチンの一実施例を示している。このルーチンは、一定クランク角毎、或いは、一定時間毎の割り込みにより実行される。図１３を参照すると、まず初めに、ステップ１３０において、ノッキングセンサ１８の出力信号の取り込み期間中であるか否かが判別される。この実施例では、ノッキングセンサ１８の出力信号の取り込み期間が圧縮上死点から圧縮上死点後９０°までの間とされている。従って、クランク角が圧縮上死点前であるときには、処理サイクルを終了し、クランク角が圧縮上死点を過ぎるとノッキングセンサ１８の出力信号の取り込み期間中であると判別され、ステップ１３１に進む。

ステップ１３１では、ノッキングセンサ１８の出力信号が、対応するＡＤ変換器３６および例えば５KHzから２５KHz程度の周波数の入力信号のみを通過させるデジタルバンドパスフィルタ３７を介して取り込まれる。即ち、フィルタリング後のノッキングセンサ１８の出力値が取り込まれる。次いで、ステップ１３２では、取り込まれたフィルタリング後のノッキングセンサ１８の出力値が電子制御ユニット３０のメモリ３２内に記憶される。次いで、ステップ１３３では、ノッキングセンサ１８の出力信号の取り込み期間が終了したか否かが判別される。即ち、クランク角が圧縮上死点後９０°に達したか否かが判別される。ステップ１３３において、ノッキングセンサ１８の出力信号の取り込み期間が終了していないと判別されたとき、即ち、クランク角が圧縮上死点後９０°に達していないと判別されたときには処理サイクルを終了する。

これに対し、ステップ１３３において、ノッキングセンサ１８の出力信号の取り込み期間が終了したと判別されたとき、即ち、クランク角が圧縮上死点後９０°に達したと判別されたときにはステップ１３４に進む。このときメモリ３２内には図６Ｂの黒丸で示すようなフィルタリング後のノッキングセンサ１８の出力値が記憶されている。ステップ１３４では、このメモリ３２内に記憶されているフィルタリング後のノッキングセンサ１８の出力値が、学習済みの第１のニューラルネットワーク２０の入力層 ( L＝１) の各ノードに入力される。このとき、学習済みの第１のニューラルネットワーク２０からはノッキング強度を代表する値の推定値ｙe が出力され、従って、ステップ１３５に示されているように、ノッキング強度を代表する値の推定値ｙe が算出されることになる。

図１４は、図１２のステップ１２０において行われるノッキング強度代表値の推定処理ルーチンの別の実施例を示している。このルーチンも、一定クランク角毎、或いは、一定時間毎の割り込みにより実行される。なお、図１４のステップ１４０からステップ１４３の内容は、図１３のステップ１３０からステップ１３３の内容と同一である。即ち、図１４を参照すると、まず初めに、ステップ１４０において、ノッキングセンサ１８の出力信号の取り込み期間中であるか否かが判別される。この実施例でも、ノッキングセンサ１８の出力信号の取り込み期間が圧縮上死点から圧縮上死点後９０°までの間とされている。従って、クランク角が圧縮上死点前であるときには、処理サイクルを終了し、クランク角が圧縮上死点を過ぎるとノッキングセンサ１８の出力信号の取り込み期間中であると判別され、ステップ１４１に進む。

ステップ１４１では、ノッキングセンサ１８の出力信号が、対応するＡＤ変換器３６および例えば５KHzから２５KHz程度の周波数の入力信号のみを通過させるデジタルバンドパスフィルタ３７を介して取り込まれる。即ち、フィルタリング後のノッキングセンサ１８の出力値が取り込まれる。次いで、ステップ１４２では、取り込まれたフィルタリング後のノッキングセンサ１８の出力値が電子制御ユニット３０のメモリ３２内に記憶される。次いで、ステップ１４３では、ノッキングセンサ１８の出力信号の取り込み期間が終了したか否かが判別される。即ち、クランク角が圧縮上死点後９０°に達したか否かが判別される。ステップ１４３において、ノッキングセンサ１８の出力信号の取り込み期間が終了していないと判別されたとき、即ち、クランク角が圧縮上死点後９０°に達していないと判別されたときには処理サイクルを終了する。

これに対し、ステップ１４３において、ノッキングセンサ１８の出力信号の取り込み期間が終了したと判別されたとき、即ち、クランク角が圧縮上死点後９０°に達したと判別されたときにはステップ１４４に進む。ステップ１４４では、メモリ３２内に記憶されているフィルタリング後のノッキングセンサ１８の出力値に基づいて、図６Ｃを参照しつつ説明したように、クランク角度５°毎に分割された各区間内におけるフィルタリング後のノッキングセンサ１８の出力値の積分値（負側の積分値も正の値とする）が算出される。次いで、ステップ１４５では、フィルタリング後のノッキングセンサ１８の出力値の算出された各積分値が、学習済みの第１のニューラルネットワーク２０の入力層 ( L＝１) の各ノードに入力される。このとき、学習済みのニューラルネットワーク２０からはノッキング強度を代表する値の推定値ｙe が出力され、従って、ステップ１４６に示されているように、ノッキング強度を代表する値の推定値ｙe が算出されることになる。

図１５は、図１３のステップ１３５において算出されたノッキング強度を代表する値の推定値ｙe、又は図１４のステップ１４６において算出されたノッキング強度を代表する値の推定値ｙe に基づき電子制御ユニット３０内で実行されるノッキング判定ルーチンを示している。このルーチンは、一定クランク角毎の割り込みにより実行される。図１５を参照すると、まず初めに、ステップ１５０において、ノッキング強度を代表する値に対する閾値Ｍij が読み込まれる。この閾値Ｍij は、図１６に示されるように、機関負荷Ｌと機関回転数Ｎに応じて複数個に分割された各機関運転領域に対して予め設定されている。

次いで、ステップ１５１では、算出されたノッキング強度を代表する値の推定値ｙe が閾値Ｍij よりも大きいか否かが判別される。算出されたノッキング強度を代表する値の推定値ｙe が閾値Ｍij よりも大きくないと判別されたときにはステップ１５２に進んで、ノッキングが発生していないと判別され、次いで処理サイクルを終了する。なお、以下、このようにノッキング強度を代表する値の推定値ｙe が閾値Ｍij よりも大きくないときには、たとえ弱いノッキングが発生していたとしても、処理説明の便宜上、ノッキングが発生していないと言う。これに対し、算出されたノッキング強度を代表する値の推定値ｙe が閾値Ｍij よりも大きいと判別されたときにはステップ１５３に進んで、ノッキングが発生している判別され、次いで処理サイクルを終了する。

図１７は、図１５に示されるノッキング判定ルーチンにおける判定結果に基づき電子制御ユニット３０内で実行される点火制御ルーチンを示している。このルーチンは、一定クランク角毎の割り込みにより実行される。図１７を参照すると、まず初めに、ステップ１６０において、基準点火時期ＩＧ（ＢＴＤＣ）が算出される。この基準点火時期ＩＧは、機関負荷Ｌと機関回転数Ｎの関数として、図１８に示されるようなマップの形で予めメモリ３２内に記憶されている。次いで、ステップ１６１では、図１５に示されるノッキング判定ルーチンにおける判定結果に基づき、ノッキングが発生したか否かが判別される。ノッキングが発生したと判別されたときにはステップ１６２に進む。

ステップ１６２では、点火時期を遅角させるために点火時期遅角量ΔＩＧに一定量αが加算される。次いで、ステップ１６６では、基準点火時期ＩＧから点火時期遅角量ΔＩＧを減算することによって最終的な点火時期ＩＧＯが算出される。この最終的な点火時期ＩＧＯに基づいて。点火栓１１による点火作用が制御される。このとき点火時期は一定量αだけ遅角される。一方、ステップ１６１において、ノッキングが発生していないと判別されたときにはステップ１６３に進み、点火時期を進角させるために点火時期遅角量ΔＩＧから一定量βが減算される。

次いで、ステップ１６４では、点火時期遅角量ΔＩＧか負になったか否かが判別され、点火時期遅角量ΔＩＧか負になっていないときにはステップ１６６に進んで、最終的な点火時期ＩＧＯが算出される。このときには、点火時期は一定量βだけ進角される。なお、本発明による実施例では、一定量αは一定量βよりも大きな値に設定されている。言い換えると、遅角量αは進角量βよりも大きな値に設定されている。一方、ステップ１６４において、点火時期遅角量ΔＩＧか負になったと判別されたときにはステップ１６５に進んで点火時期遅角量ΔＩＧが零とされた後、ステップ１６６に進む。このときには、点火時期が基準点火時期ＩＧとされる。

図１２から図１５および図１７に示されるノッキング処理は、気筒毎に夫々独立して実行され、図１９は、或る一つの気筒におけるノッキング処理の一例を示している。図１９には、サイクル番号と、点火の行われるタイミングと、第１の学習済みニューラルネットワーク２０への入力値ｘ_１、ｘ_２・・・ｘ_ｎ−１、ｘ_ｎの取り込み作業が行なわれる期間と、第１の学習済みニューラルネットワーク２０から出力されるノッキング強度を代表する値の推定値ｙe と、ノッキングが発生したか否かの判定が行われるノッキング判定時期と、ノッキングの判定結果と、点火時期とが示されている。

本発明による実施例では、吸気上死点ＴＤＣから次の吸気上死点ＴＤＣまでの７２０クランク角度が１サイクルとされており、機関の運転が開始されたときからのサイクル番号がｋ_１、ｋ_２、・・・ｋ_ｎ、ｋ_ｎ＋１、ｋ_ｎ＋２、ｋ_ｎ＋３・・・とされている。図１９には、これらサイクル番号が、ｋ_ｎ、ｋ_ｎ＋１、ｋ_ｎ＋２、ｋ_ｎ＋３となったときが示されており、ノッキング処理を容易に理解し得るように、以下、これらサイクル番号ｋ_ｎ、ｋ_ｎ＋１、ｋ_ｎ＋２、ｋ_ｎ＋３を用いて説明する。本発明による実施例では、ピストン４の位置が圧縮上死点ＴＤＣに達する毎に、図１２に示されるノッキング処理が実行される。ノッキング処理が実行されると、図１３又は図１４に示されるノッキング強度代表値の推定ルーチンが実行され、第１の学習済みニューラルネットワーク２０への入力値ｘ_１、ｘ_２・・・ｘ_ｎ−１、ｘ_ｎの取り込み作業が行なわれる。

第１の学習済みニューラルネットワーク２０への入力値ｘ_１、ｘ_２・・・ｘ_ｎ−１、ｘ_ｎの取り込み作業が完了すると、第１の学習済みニューラルネットワーク２０からは、ノッキング強度を代表する値の推定値ｙe が出力される。第１の学習済みニューラルネットワーク２０からノッキング強度を代表する値の推定値ｙe が出力されると、図１５に示されるノッキング判定ルーチンが実行され、ノッキングが発生したか否かの判定が行われる。このとき、ノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越えると、ノッキングが発生したと判定される。ノッキングが発生したか否かの判定が行われると、図１７に示される点火制御ルーチンが実行される。

図１９からわかるように、ノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越えると、次のサイクルでは点火時期が遅角される。点火時期が遅角されると多くの場合、ノッキング強度を代表する値の推定値ｙe は閾値Ｍij 以下まで低下する。一方、点火時期が遅角されても、ノッキング強度を代表する値の推定値ｙe が閾値Ｍij以下まで低下しない場合もある。図１９は、このように点火時期が遅角されても、ノッキング強度を代表する値の推定値ｙe が閾値Ｍij以下まで低下しない場合の例を示している。

即ち、図１９に示される例では、サイクルが番号ｋ_ｎのサイクルとなったときにノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越えたと判別され、それにより次のサイクルｋ_ｎ＋１において点火時期が遅角される。次いで、次のサイクルｋ_ｎ＋１でも、ノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越えたと判別され、それにより更に次のサイクルｋ_ｎ＋２においても点火時期が遅角される。サイクルｋ_ｎ＋２において点火時期が遅角されると、ノッキング強度を代表する値の推定値ｙe が閾値Ｍij 以下に低下し、それによりノッキングが発生していないと判別される。ノッキングが発生していないと判別されると、次のサイクルｋ_ｎ＋３では点火時期が進角され、ノッキングが発生しない限り、点火時期は進角され続ける。なお、図１９からわかるように、本発明による実施例では、遅角量に比べて進角量は小さくされている。

さて、前述したように、図５に示される内燃機関は、第１のニューラルネットワーク２０の重みを学習するために、図１に示される内燃機関に圧力センサ１９を取り付けた内燃機関である。これら図１或いは図５に示される内燃機関において、第１の学習済みニューラルネットワーク２０を用い、図１９を参照しつつ説明したようなノッキング処理を行った場合には、何ら問題は生じない。ところが他の内燃機関、例えば、市販用の内燃機関において、この第１の学習済みニューラルネットワーク２０を用い、図１９を参照しつつ説明したようなノッキング処理を行うと問題を生じる。そこで次に、どのような問題を生ずるかについて図２０Ａから図２３を参照しつつ説明する。

図２０Ａは燃焼室５内の燃焼圧とクランク角（ＡＴＤＣ）との関係を示している。なお、図２０Ａにおいて実線はノッキングが発生したときの燃焼室５内の燃焼圧の変化を示しており、破線はノッキングが発生していないときの燃焼室５内の燃焼圧の変化を示している。図２０Ａからわかるように、ノッキングが発生すると、圧縮上死点ＴＤＣ後に、燃焼室５内の燃焼圧が急激に上昇する。図２０Ｂには、このときの実際の燃焼室５内の燃焼圧の上昇量とノッキング強度との関係が各点で示されている。図２０Ｂに示されるように燃焼室５内の燃焼圧の上昇量が大きくなるほど、ノッキング強度が大きくなる。即ち、燃焼室５内の燃焼圧の上昇量が大きくなるほど、ノッキング強度を代表する値の推定値ｙe が大きくなる。従って、第１の学習済みニューラルネットワーク２０の入力値と出力値は、これら燃焼室５内の燃焼圧の上昇量とノッキング強度との関係を示していることになる。

図２１は、図５に示される内燃機関において、第１の学習済みニューラルネットワーク２０を用い、図１９に示されるノッキング処理が行われたときの圧力センサ１９の出力値と、ノッキングセンサ１８の出力値と、燃焼室５内の燃焼圧と、ノッキング強度を代表する値の推定値ｙe の夫々の値の変化を示している。なお、図２１において、（ａ）は図１９のサイクルｋ_ｎにおけるこれらの値の変化を示しており、（ｂ）は図１９のサイクルｋ_ｎ＋１におけるこれらの値の変化を示しており、（ｃ）は図１９のサイクルｋ_ｎ＋２におけるこれらの値の変化を示している。図１９からわかるように、サイクルが、ｋ_ｎ＋１、ｋ_ｎ＋２へと進むにつれて、点火時期が次第に遅角され、図２１からわかるように、サイクルが、ｋ_ｎ＋１、ｋ_ｎ＋２へと進むにつれて、圧力センサ１９の出力値、ノッキングセンサ１８の出力値、燃焼室５内の燃焼圧の上昇量およびノッキング強度を代表する値の推定値ｙe が次第に低下していく。

このように図１９および図２１から、点火時期の遅角量を増大させると、燃焼室５内の燃焼圧の上昇量およびノッキング強度を代表する値の推定値ｙe が低下することがわかる。この場合、点火時期の遅角速度を増大させると、燃焼室５内の燃焼圧の上昇量およびノッキング強度を代表する値の推定値ｙe が低下することもわかっている。即ち、図２２Ａの各点で示されるように、点火時期の遅角量又は遅角速度を増大させると、燃焼室５内の燃焼圧の上昇割合が低下し、図２２Ｂの各点で示されるように、点火時期の遅角量又は遅角速度を増大させると、ノッキング強度を代表する値の推定値ｙe が低下することがわかっている。

一方、図２３は、図５に示される内燃機関以外の他の内燃機関、例えば、市販用の内燃機関に圧力センサ１９を取り付け、第１の学習済みニューラルネットワーク２０を用いて、図１９を参照しつつ説明したようなノッキング処理が行なわれた場合の圧力センサ１９の出力値と、ノッキングセンサ１８の出力値と、燃焼室５内の燃焼圧と、ノッキング強度を代表する値の推定値ｙe の夫々の値の変化を示している。なお、この図２３においても、（ａ）、（ｂ）、（ｃ）は、図１９の各サイクルｋ_ｎ、ｋ_ｎ＋１，ｋ_ｎ＋２におけるこれらの値の変化を示している。

さて、前述したように、機関の構成部品には公差があり、従って、機関毎に機関の構成部品の寸法がばらつくために、機関毎に異なる機関の振動が発生する。しかしながら、第１の学習済みニューラルネットワーク２０は、このように機関毎に発生する異なる機関の振動については、重みの学習がなされておらず、従って、このように重みの学習がなされていない機関の振動、即ち、未学習の機関の振動が発生すると、第１の学習済みニューラルネットワーク２０において、ノッキングにより機関本体１に振動が発生していると誤判断される場合がある。図２３の（a）は、一例として、サイクルｋ_ｎにおいてノッキングおよび未学習の機関の振動が発生していてノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越え、その結果、ノッキングが発生していると判別された場合を示している。

この場合、次のサイクルｋ_ｎ＋１では、点火時期が遅角され、それにより図２３の（ｂ）で示されるように燃焼室５内の燃焼圧の上昇量は低下し、圧力センサ１９の出力値も低下する。しかしながら、点火時期の影響を受けない未学習の機関の振動が発生している場合には、点火時期が遅角されたとしても未学習の機関の振動は抑制されず、従って、図２３に示される例では、未学習の機関の振動は、点火時期が遅角された後も、図２３の（ｂ）で示されるように発生し続ける。その結果、図２３に示される例では、図２３の（ｂ）で示されるように、未学習の機関の振動によってノッキング強度を代表する値の推定値ｙe が閾値Ｍij 以上となり、その結果、依然として、ノッキングが発生していると判別されることになる。

この場合、次のサイクルｋ_ｎ＋２では、再度点火時期が遅角され、従って、図２３の（ｃ）で示されるように燃焼室５内の燃焼圧の上昇量は低下し、圧力センサ１９の出力値も低下する。しかしながらこのように点火時期が遅角されても、未学習の機関の振動が抑制されず、その結果、図２３に示される例では、点火時期が遅角された後も、図２３の（ｃ）で示されるように未学習の機関の振動が発生し続ける。その結果、図２３の（ｃ）で示されるように、未学習の機関の振動によってノッキング強度を代表する値の推定値ｙe が閾値Ｍij 以上となり、その結果、依然としてノッキングが発生していると判別されることになる。このように、未学習の機関の振動が発生し続けると、ノッキング強度を代表する値の推定値ｙe が閾値Ｍijを越え続け、それにより点火時期の遅角作用が続行され続ける場合がある。このような場合には、点火時期が過剰に遅角し、その結果、機関の出力が大幅に低下してしまうという問題を生ずることになる。

一方、ノッキングが発生していないが、未学習の機関の振動が発生したために、ノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越え、ノッキングが発生していると判別された場合も同様である。この場合も、未学習の機関の振動が発生し続けると、ノッキング強度を代表する値の推定値ｙe が閾値Ｍijを越え続け、点火時期の遅角作用が続行され続ける場合が生ずる。この場合にも、点火時期が過剰に遅角し、それにより機関の出力が大幅に低下してしまうという問題を生ずることになる。

本発明では、図５に示される内燃機関以外の他の内燃機関、例えば、市販用の内燃機関において、第１の学習済みニューラルネットワーク２０を用いてノッキング処理を行なった場合に、未学習の機関の振動の発生により点火時期が過剰に遅角することがないように点火時期が制御される。次に、このことについて、図２４Ａから図２５Ｂを参照しつつ説明する。なお、図２４Ａから図２５Ｂは、図１９と同様に、サイクルｋ_ｎにおいてノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越えたと判別され、それにより次のサイクルｋ_ｎ＋１において点火時期が遅角された場合におけるノッキング強度を代表する値の推定値ｙe および点火時期の変化を示している。

まず初めに、図２４Ａを参照すると、図２４Ａは、図１又は図５に示される内燃機関において、第１の学習済みニューラルネットワーク２０を用い、ノッキングセンサ１８の出力値に基づいてノッキング処理が行われたときのノッキング強度を代表する値の推定値ｙe および点火時期の変化を示している。図２４Ａに示される例は、サイクルｋ_ｎにおいてノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越えたと判別され、それにより次のサイクルｋ_ｎ＋１において点火時期がαだけ遅角され、その結果、サイクルｋ_ｎ＋１におけるノッキング強度を代表する値の推定値ｙe がΔｙｅだけ低下して、ｙｅｅになった場合を示している。なお、図２４Ａは、ノッキング強度を代表する値の推定値ｙe がｙｅｅになったときに、ノッキング強度を代表する値の推定値ｙe が閾値Ｍij 以下になった場合を示しており、この場合には、次のサイクルｋ_ｎ＋２ において、点火時期がβ１だけ進角される。

一方、図２４Ｂは、図５に示される内燃機関以外の他の内燃機関、例えば、市販用の内燃機関において、第１の学習済みニューラルネットワーク２０を用い、ノッキングセンサ１８の出力値に基づいてノッキング処理を行うようにした場合の本発明による新たなノッキング処理方法を説明するための図である。なお、図２４Ｂは、実線で示されるように、サイクルｋ_ｎにおいて、ノッキングおよび未学習の機関の振動の双方の発生によりノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越え、サイクルｋ_ｎ＋１以降は、未学習の機関の振動によりノッキング強度を代表する値の推定値ｙe が暫くの間、閾値Ｍij を越え続けている場合を示している。即ち、実際の内燃機関では、未学習の機関の振動は特定の運転状態で発生し、機関の運転状態が変化すると次第に収まるので、図２４Ｂに示されるように、ノッキング強度を代表する値の推定値ｙe は、一時的に閾値Ｍij を越え続けていたとしても、暫くすると次第に低下していく。ところが、図２４Ｂに示されるようにノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越え続けている場合、図１２から図１８に示されるルーチンによりノッキング処理が行われると、点火時期が遅角され続け、その結果、点火時期が過剰に遅角されることになる。そこで本発明では、点火時期が過剰に遅角されるのを阻止するために、点火時期が遅角されたときに、ノッキングの発生によりノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越えているのか否かを判別し、ノッキングの発生によりノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越えているのではないと判別されたときには、、次のサイクルにおいて点火時期を遅角しないようにしている。

即ち、図２４Ｂに示されるように、図５に示される内燃機関以外の他の内燃機関、例えば、市販用の内燃機関において、第１の学習済みニューラルネットワーク２０を用いてノッキング処理が行われた場合には、ノッキングのみが発生することによりノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越えたのか、未学習の機関の振動のみが発生することによりノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越えたのか、ノッキングおよび未学習の機関の振動の双方が発生することによりノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越えたのかは、わからない。

しかしながら、図２４Ａに示されるように、ノッキングの発生によりノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越え、それにより点火時期がαだけ遅角された場合には、次のサイクルにおけるノッキング強度を代表する値の推定値ｙe がΔｙｅだけ低下して、ｙｅｅになることはわかっている。

従って、図２４Ｂに示されるように、例えば、市販用の内燃機関において、第１の学習済みニューラルネットワーク２０を用いてノッキング処理を行った場合においても、ノッキングの発生のみによりノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越え、それにより点火時期がαだけ遅角された場合には、破線で示されるように、次のサイクルｋ_ｎ＋１におけるノッキング強度を代表する値の推定値ｙe はΔｙｅだけ低下して、ｙｅｅになると予測される。従って、点火時期がαだけ遅角されたときに、ノッキング強度を代表する値の推定値ｙeがこの予測値ｙeｅ 近くまで低下すれば、ノッキングが発生したと判断できることになる。また、このとき、ノッキング強度を代表する値の推定値ｙe は、ノッキング強度を正しく示していることになる。

一方、点火時期がαだけ遅角されたときに、図２４Ｂにおいて、実線で示されるように、ノッキング強度を代表する値の推定値ｙeが予測値ｙeｅ 近くまで低下しないときには、サイクルｋ_ｎにおいて、ノッキングの発生のみによるのではなく、未学習の機関の振動の発生により、或いは、ノッキングおよび未学習の機関の振動の双方の発生によりノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越えていることになる。従って、このときには、ノッキング強度を代表する値の推定値ｙe が、ノッキング強度を正しく示していないことになる。即ち、点火時期がαだけ遅角されたときのノッキング強度を代表する値の推定値ｙｅｅと第１の学習済みニューラルネットワーク２０から出力されるノッキング強度を代表する値の推定値ｙe との差が予め定められた設定値よりも小さいときには、ノッキング強度を代表する値の推定値ｙe は、ノッキング強度を正しく示しており、この差が予め定められた設定値よりも大きいときには、ノッキング強度を代表する値の推定値ｙe は、ノッキング強度を正しく示していないことになる。

そこで本発明では、市販用の内燃機関において、第１の学習済みニューラルネットワーク２０を用いてノッキング処理を行う場合には、図２４Ｂに示されるように、サイクルｋ_ｎにおいて、ノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越えたと判断された場合、次のサイクルｋ_ｎ＋１では、点火時期がαだけ遅角されたときのノッキング強度を代表する値の推定値ｙｅｅと第１の学習済みニューラルネットワーク２０から出力されるノッキング強度を代表する値の推定値ｙe との差が予め定められた設定値よりも小さいか、又は、大きいかに応じて、更に次のサイクルｋ_ｎ＋２における点火時期の遅角制御が行われる。

即ち、この差が予め定められた設定値よりも小さいときには、即ち、ノッキング強度を代表する値の推定値ｙe が、ノッキング強度を正しく示しているときには、ノッキング強度を代表する値の推定値ｙe に基づいて、点火時期の遅角制御が行われる。この場合、ノッキング強度を代表する値の推定値ｙe が、閾値Ｍij よりも大きいときには、点火時期の遅角作用が行われ、ノッキング強度を代表する値の推定値ｙe が、閾値Ｍij よりも小さいときには、点火時期の遅角作用は行われず、このとき、次のサイクルｋ_ｎ＋２ では、点火時期がβ１だけ進角される。一方、この差が予め定められた設定値よりも大きいときには、即ち、ノッキング強度を代表する値の推定値ｙe が、ノッキング強度を正しく示していないときには、点火時期の遅角制御は行われず、様子をみるために、図２４Ｂにおいて実線で示されるように点火時期がそのまま維持される。次いで、ノッキング強度を代表する値の推定値ｙe が値Ｍij よりも低下した後に、点火時期がβ１ずつ進角され始める。

なお、本発明による実施例では、図２４Ｂに示されるように、未学習の機関の振動によりノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越え続けている場合における点火時期進角量β２は、図２４Ａに示されるように、ノッキングの発生のみによりノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越えた場合における点火時期進角量β１に比べて小さくされている。即ち、未学習の機関の振動によりノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越え続けている場合には、点火時期はゆっくりと進角される。

一方、図２５Ａは、図２４Ａと同様に、図１又は図５に示される内燃機関において、第１の学習済みニューラルネットワーク２０を用い、ノッキングセンサ１８の出力値に基づいてノッキング処理が行われたときのノッキング強度を代表する値の推定値ｙe および点火時期の変化を示しており、図２５Ｂは、図２５Ａと同様に、図５に示される内燃機関以外の他の内燃機関、例えば、市販用の内燃機関において、第１の学習済みニューラルネットワーク２０を用い、ノッキングセンサ１８の出力値に基づいてノッキング処理が行われたときのノッキング強度を代表する値の推定値ｙe および点火時期の変化を示している。なお、図２５Ａおよび図２５Ｂは、ノッキング強度を代表する値の推定値ｙe が閾値Ｍij 以下となるまで点火時期がαだけ２回にわたり遅角された場合を示しており、点火時期が２回にわたり遅角されることを除けば、図２４Ａおよび図２４Ｂに示される場合と同様である。

即ち、図２５Ａに示されるように、ノッキングの発生によりノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越え、それにより点火時期がαだけ遅角された場合には、次のサイクルｋ_ｎ＋１におけるノッキング強度を代表する値の推定値ｙe がΔｙｅ１だけ低下して、ｙｅｅ１となる。このとき、まだノッキングの発生が収まっていないために、ノッキング強度を代表する値の推定値ｙe が依然として閾値Ｍij が越えており、それにより点火時期が再度αだけ遅角される。その結果、次のサイクルｋ_ｎ＋２におけるノッキング強度を代表する値の推定値ｙe がΔｙｅ２だけ低下して、ｙｅｅ２となる。

本発明では、この場合も、図２５Ｂの実線で示されるように、点火時期がαだけ遅角されたときのノッキング強度を代表する値の推定値ｙｅｅ１と第１の学習済みニューラルネットワーク２０から出力されるノッキング強度を代表する値の推定値ｙe との差が予め定められた設定値よりも大きいときには、即ち、ノッキング強度を代表する値の推定値ｙe が、ノッキング強度を正しく示していないときには、点火時期の遅角制御は行われず、様子をみるために、図２５Ｂにおいて実線で示されるように点火時期がそのまま維持される。次いで、ノッキング強度を代表する値の推定値ｙe が値Ｍij よりも低下した後に、点火時期がβ１ずつ進角され始める。

このように本発明では、市販用の内燃機関において、第１の学習済みニューラルネットワーク２０を用いてノッキング処理を行う場合、点火時期が遅角されたときのノッキング強度を代表する値の推定値ｙeの予測値ｙeｅと第１の学習済みニューラルネットワーク２０から出力されるノッキング強度を代表する値の推定値ｙe との差が予め定められた設定値よりも小さいか、又は、大きいかに応じて、次のサイクルにおける点火時期の遅角制御が行われる。この場合、本発明による実施例では、点火時期が遅角されたときのノッキング強度を代表する値の推定値ｙe の予測低下量Δｙｅ、或いはノッキング強度を代表する値の推定値ｙeの予測値ｙeｅが、第２のニューラルネットワークを用いて推定される。

図２６から図３３は、点火時期が遅角されたときのノッキング強度を代表する値の推定値ｙe の予測低下量Δｙｅを、第２のニューラルネットワークを用いて推定するようにした本発明による第１の実施例を示しており、図３４から図４２は、点火時期が遅角されたときのノッキング強度を代表する値の推定値ｙｅの予測値ｙｅｅを、第２のニューラルネットワークを用いて推定するようにした本発明による第２の実施例を示している。そこで最初に、図２６から図３４を参照しつつ、第１の実施例について説明する。

図２７は、第１の実施例において用いられている第２のニューラルネットワーク５０を示している。図２７を参照すると、この第２のニューラルネットワーク５０においても、図３に示されるニューラルネットワークと同様に、L＝１は入力層、L＝２および L＝３は隠れ層、L＝４は出力層を夫々示している。この第１の実施例では、図２７に示されるように、入力層 ( L＝１) が５個のノードからなり、５個の入力値ｘｘ_１、ｘｘ_２、ｘｘ_３、ｘｘ_４、ｘｘ_５が、入力層 ( L＝１) の各ノードに入力される。一方、図２７には隠れ層 ( L＝２)および隠れ層 ( L＝３)が記載されているが、これら隠れ層の層数は、１個又は任意の個数とすることができ、またこれら隠れ層のノードの数も任意の個数とすることができる。なお、出力層 ( L＝４) のノードの数は１個とされており、出力層のノードからの出力値がΔｙｅで示されている。

次に、図８における入力値ｘｘ_１、ｘｘ_２、ｘｘ_３、ｘｘ_４、ｘｘ_５と、出力値Δｙｅについて説明する。最初に出力値Δｙｅについて説明すると、この出力値Δｙｅは、ノッキングのみの発生により点火時期が遅角されたときのノッキング強度を代表する値の推定値ｙe の予測低下量を示している。一方、入力値ｘｘ_１、ｘｘ_２、ｘｘ_３、ｘｘ_４、ｘｘ_５は、ノッキング強度を代表する値の推定値ｙe の低下量に強い影響を与える入力パラメータを示しており、これら入力パラメータの一覧表が図２６に示されている。図２６に示されるように、入力パラメータは、機関の運転状態を表すパラメータと、点火時期の遅角量α又は進角量β（β１又はβ２）と、一つ前のサイクルにおけるノッキング強度を代表する値の推定値ｙeｏからなる。

なお、この第１の実施例では、機関の運転状態を表すパラメータが、機関回転数ＮＥと、機関負荷Ｌと、ＥＧＲ率ＥＲからなる。即ち、第１の実施例では、機関回転数ＮＥ、機関負荷Ｌ、ＥＧＲ率ＥＲ、点火時期の遅角量α又は進角量β、および一つ前のサイクルにおけるノッキング強度を代表する値の推定値ｙeｏが夫々、入力値ｘｘ_１、ｘｘ_２、ｘｘ_３、ｘｘ_４、ｘｘ_５とされ、これら入力値ｘｘ_１、ｘｘ_２、ｘｘ_３、ｘｘ_４、ｘｘ_５が、第２のニューラルネットワーク５０の入力層の各ノードに入力される。

図２８は、図９に示される訓練データセットを作成する際に用いられた図５に示される機関本体１を用いて、第２のニューラルネットワーク５０用訓練データセットを作成するために取得されたデータの一覧表を示している。これらのデータは、図５に示される機関本体１において、図１２から図１８に示されるルーチンを用いたノッキング処理を行いつつ、機関回転数ＮＥ、機関負荷Ｌ、ＥＧＲ率ＥＲおよび点火時期の遅角量αおよび進角量βの種々の組み合わせについてノッキングの発生しない運転状態とノッキングの発生する運転状態の双方の状態が生ずるように機関を運転させたときの、サイクル毎の機関回転数ＮＥ、機関負荷Ｌ、ＥＧＲ率ＥＲ、ノッキング強度を代表する値の推定値ｙe および点火時期の遅角量αおよび進角量βを示している。これらのデータは、例えば、電子制御ユニット３０のメモリ３２内に一時的に記憶される。

なお、これらデータを取得する際、本発明による実施例では、点火時期遅角量αは一定に維持されており、点火時期進角量βは、進角量β１又は進角量β２のいずれかとされる。また、図２８に示される一覧表には、点火時期が遅角量αだけ遅角されたときのノッキング強度を代表する値の推定値ｙe の低下量Δｙｅが含まれている。この低下量Δｙｅは、点火時期が遅角量αだけ遅角されたときのノッキング強度を代表する値の推定値ｙe と、一つ前のサイクルにおけるノッキング強度を代表する値の推定値ｙeｏ との差から求められる。例えば、サイクルｋ_ｎ＋１でおいて点火時期が遅角量αだけ遅角されたときには、サイクルｋ_ｎのときのノッキング強度を代表する値の推定値ｙeｎから、サイクルｋ_ｎ＋１のときのノッキング強度を代表する値の推定値ｙe n+1 を減算することによって、低下量Δｙe n+1が求められる。この低下量Δｙｅの算出は、例えば、ＣＰＵ３３において行われる。

このようにして、点火時期が遅角量αだけ遅角されたときに、図２８に示されるようにノッキング強度を代表する値の推定値ｙe の低下量Δｙｅが算出され、算出された低下量Δｙｅが、遅角量αと共に、例えば、電子制御ユニット３０のメモリ３２内に一時的に記憶される。このノッキング強度を代表する値の推定値ｙe の低下量Δｙｅは、ノッキング強度を代表する値の推定値ｙe の実際の低下量を表しており、従って、このノッキング強度を代表する値の推定値ｙe の低下量Δｙｅは、第２の学習済みニューラルネットワーク５０の重みを学習する際の正解値、即ち、教師データとされる。

図２９は、第２のニューラルネットワーク５０の重みを学習するための訓練データセットを示している。この訓練データセットは、第１の実施例では、図２８に示される一覧表から、点火時期が遅角量αだけ遅角されたときのサイクルにおける機関回転数ＮＥ、機関負荷Ｌ、ＥＧＲ率ＥＲ、点火時期の遅角量αに関するデータ、および点火時期が遅角量αだけ遅角されたときのサイクルの一つ前のサイクルにおけるノッキング強度を代表する値の推定値ｙeｏに関するデータを抜き出すことによって作成されており、この訓練データセットは、図２８に示される一覧表から抜き出された入力値ｘｘ_１、ｘｘ_２、ｘｘ_３、ｘｘ_４、ｘｘ_５と、教師データｙｔとの関係を表すｍ個のデータからなる。なお、この場合、図２８に示される機関回転数ＮＥ、機関負荷Ｌ、ＥＧＲ率ＥＲ、点火時期の遅角量α、および一つ前のサイクルにおけるノッキング強度を代表する値の推定値ｙeｏが夫々、入力値ｘｘ_１、ｘｘ_２、ｘｘ_３、ｘｘ_４、ｘｘ_５とされ、図２８に示されるノッキング強度を代表する値の推定値ｙe の低下量Δｙｅが、教師データｙｔとされる。

この第１の実施例では、図２７に示される第２のニューラルネットワーク５０の入力層 ( L＝１) のノード数、隠れ層 ( L＝２)および隠れ層 ( L＝３)のノード数および出力層 ( L＝４) のノード数、および、図２９に示される訓練データセットが、学習装置２３のメモリ２４に記憶され、学習装置２３において、図１０を参照しつつ既に説明した学習処理方法と同様な学習処理方法により、図１０に示される学習処理ルーチンを用いて、第２のニューラルネットワーク５０の重みの学習が行われる。それにより、第２の学習済みニューラルネットワーク５０、即ち、ノッキング強度を代表する値の推定値ｙe の低下量Δｙｅの推定モデルが作成される。

本発明による第１の実施例では、第１の学習済みニューラルネットワーク２０により生成された、ノッキング強度を代表する値の推定モデルと、第２の学習済みニューラルネットワーク５０により生成された、ノッキング強度を代表する値の推定値ｙe の低下量Δｙｅの推定モデルを用いて、市販車両の機関におけるノッキング処理が行われる。そのためにこれらノッキング強度を代表する値の推定モデルおよびノッキング強度を代表する値の推定値ｙe の低下量Δｙｅの推定モデル、即ち、第１の学習済みニューラルネットワーク２０および第２の学習済みニューラルネットワーク５０が、市販車両の電子制御ユニット３０内に格納される。これら推定モデルの市販車両の電子制御ユニット３０への格納作用は、図１１を参照しつつ既に説明した方法と同様な方法でもって、図１１に示される電子制御ユニットへのデータ読み込みルーチンを用いて行われる。

このように第１の学習済みニューラルネットワーク２０および第２の学習済みニューラルネットワーク５０が、市販車両の電子制御ユニット３０に格納されると、電子制御ユニット３０のメモリ３２内には、第１の学習済みニューラルネットワーク２０および第２の学習済みニューラルネットワーク５０が形成される。図３０は、市販車両の電子制御ユニット３０のメモリ３２内に形成された第１の学習済みニューラルネットワーク２０および第２の学習済みニューラルネットワーク５０を用いて、市販車両の機関運転時に実行されるノッキング処理を示している。このノッキング処理は、気筒毎に、個別に、一サイクル毎に実行される。本発明による第１の実施例では、このノッキング処理は、例えば、クランク角が圧縮上死点に達したときに開始される。

図３０を参照すると、まず初めにステップ２２０において、圧力センサ１９の出力値から求められるノッキング強度を代表する値、即ち、ノッキング強度代表値の推定処理が行われる。次いで、ステップ２２１では、ノッキングが生じているか否かのノッキング判定処理が行われ、次いで、ステップ２２２では、点火制御が行われる。ステップ２２０において行われるノッキング強度代表値の推定処理は、図１３又は図１４のいずれか一方に示されるノッキング強度代表値の推定ルーチンにより行われる。これらのノッキング強度代表値の推定ルーチンについては既に説明したので、これらノッキング強度代表値の推定ルーチンについては説明を省略する。なお、これらのノッキング強度代表値の推定ルーチンが実行されると、第１の学習済みニューラルネットワーク２０からはノッキング強度を代表する値の推定値ｙe が出力される。

図３１および図３２は、図３０のステップ２２０において第１の学習済みニューラルネットワーク２０からノッキング強度を代表する値の推定値ｙe が出力されたときに、この第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙe に基づき、図３０のステップ２２１で行われるノッキング判定ルーチンを示している。図３１を参照すると、まず初めに、ステップ２３０において、ノッキング強度を代表する値に対する閾値Ｍij が読み込まれる。この閾値Ｍij は、前述したように、図１６に示される如く、機関負荷Ｌと機関回転数Ｎに応じて複数個に分割された各機関運転領域に対して予め設定されている。

次いで、ステップ２３１では、ノッキングが発生したか否かの判定を中止すべきときにセットされる判定中止フラグがセットされているか否かが判別される。この判定中止フラグは、通常、リセットされており、従って、ステップ２３２に進む。ステップ２３２では、ノッキングが発生したと判別されたときにセットされるノッキング発生フラグがセットされているか否かが判別される。このノッキング発生フラグは、通常、リセットされており、従って、ステップ２３７にジャンプする。ステップ２３７では、第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙe が閾値Ｍij よりも大きいか否かが判別される。ノッキング強度を代表する値の推定値ｙe が閾値Ｍij よりも大きくないと判別されたときにはステップ２３９に進んで、ノッキング発生フラグがリセットされる。次いで、処理サイクルを終了する。

一方、ステップ２３７において、第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙe が閾値Ｍij よりも大きいと判別されたときにはステップ２３８に進んで、ノッキング発生フラグがセットされる。ノッキング発生フラグがセットされると、後述するように、次のサイクルにおける点火時期がαだけ遅角される。また、ノッキング発生フラグがセットされると、次のサイクルでは、ステップ２３２からステップ２３３に進む。ステップ２３３では、機関回転数ＮＥ、機関負荷Ｌ、ＥＧＲ率ＥＲ、点火時期の遅角量α、および一つ前のサイクルにおいて第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙeｏ が第２の学習済みニューラルネットワーク５０の入力層の各ノードに入力され、それにより第２の学習済みニューラルネットワーク５０からは、点火時期が遅角されたときのノッキング強度を代表する値の推定値ｙe の予測低下量Δｙｅが出力される。

次いで、ステップ２３４では、一つ前のサイクルにおいて第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙeｏから、ノッキング強度を代表する値の推定値ｙe の予測低下量Δｙｅを減算することによって、ノッキング強度を代表する値の推定値ｙe の予測値ｙｅe が算出される。次いで、ステップ２３５では、未学習の機関の振動が発生しているか否かを判別するために、ノッキング強度を代表する値の推定値ｙe の予測値ｙｅeと第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙe との差（ｙｅe ―ｙｅ）の絶対値が設定値Ｓよりも大きいか否かが判別される。

即ち、ノッキングのみが発生しているときには、ノッキング強度を代表する値の推定値ｙe の予測値ｙｅeと第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙe とは、ほぼ等しくなる筈である。従って、ノッキングのみが発生しているときには、差（ｙｅe ―ｙｅ）の絶対値が設定値Ｓよりも小さくなる。従って、このときには、ステップ２３６に進んで、未学習の機関の振動、即ち、ノイズが発生していることを示すノイズフラグがリセットされる。ノイズフラグがリセットされると、後述するように、点火時期の進角量がβ１とされる。

一方、このときには、第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙe は、正しいノッキング強度を示している。従って、このときには、ステップ２３７に進んで、ノッキング強度を代表する値の推定値ｙe が閾値Ｍij よりも大きいか否かが判別される。ノッキング強度を代表する値の推定値ｙe が閾値Ｍij よりも大きいと判別されたときには、ステップ２３８に進んで、ノッキング発生フラグがセットされる。このときには、次のサイクルにおける点火時期がαだけ遅角される。

一方、ノッキングおよび未学習の機関の振動が発生しているとき、或いは、未学習の機関の振動のみが発生しているときには、第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙe は、ノッキング強度を代表する値の推定値ｙe の予測値ｙｅeに比べて大きい筈である。従って、このときには、差（ｙｅe ―ｙｅ）の絶対値が設定値Ｓよりも大きくなる。このときには、ステップ２３５からステップ２４０に進んで、未学習の機関の振動、即ち、ノイズが発生していることを示すノイズフラグがセットされる。ノイズフラグがセットされると、後述するように、点火時期の進角量がβ２とされる。次いで、ステップ２４１に進む。

ステップ２４１では、ノッキング発生フラグがリセットされる。一方、このときには、第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙe は、ノッキング強度を正しく示していない。従って、このときには、点火時期の進角作用および遅角作用を一次的に中止するために、ステップ２４２において、判定中止フラグがセットされる。次いで、処理サイクルを終了する。

判定中止フラグがセットされると、次のサイクルでは、ステップ２３１からステップ２４３に進む。ステップ２４３およびステップ２４４では、ノッキング強度を代表する値の推定値ｙe が閾値Ｍij を越えた後、未学習の機関の振動の発生により、第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙe が閾値Ｍij よりも大きい状態が続いている場合に、点火時期の進角作用および遅角作用をに中止し続けると共に、ノッキングが発生したか否かの判定を中止し続ける処理が行われる。

即ち、ステップ２４３では、第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙe が閾値Ｍij よりも大きいか否かが判別される。ノッキング強度を代表する値の推定値ｙe が閾値Ｍij よりも大きいと判別されたときには、処理サイクルを終了する。これに対し、第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙe が閾値Ｍij よりも小さくなったと判別されたときには、ステップ２４４に進んで、判定中止フラグをリセットさせる処理が行われる。次いで、ステップ２３２に進んで、ノッキング判定処理が再開される。

図３３は、図３１および図３２に示されるノッキング判定ルーチンにおける判定結果に基づき電子制御ユニット３０内で実行される点火制御ルーチンを示している。図３３を参照すると、まず初めに、ステップ２５０において、ノッキング判定ルーチンにおいて用いられている判定中止フラグがセットされているか否かが判別される。判定中止フラグがセットされていると判別されたときには処理サイクルを終了する。これに対し、判定中止フラグがセットされていないと判別されたときには、ステップ２５１に進んで、基準点火時期ＩＧ（ＢＴＤＣ）が算出される。この基準点火時期ＩＧは、前述したように、機関負荷Ｌと機関回転数Ｎの関数として、図１８に示されるようなマップの形で予めメモリ３２内に記憶されている。次いで、ステップ２５２では、図３１および図３２に示されるノッキング判定ルーチンにおいて、ノッキング発生フラグがセットされているか否かが判別される。ノッキング発生フラグがセットされていると判別されたときにはステップ２６０に進む。

ステップ２６０では、点火時期を遅角させるために点火時期遅角量ΔＩＧに一定量αが加算される。次いで、ステップ２５９では、基準点火時期ＩＧから点火時期遅角量ΔＩＧを減算することによって、次のサイクルにおける最終的な点火時期ＩＧＯが算出される。次のサイクルでは、この最終的な点火時期ＩＧＯに、点火栓１１による点火作用が行われる。このとき点火時期は一定量αだけ遅角される。一方、ステップ２５２において、ノッキング発生フラグがリセットされていると判別されたときにはステップ２５３に進み、図３１および図３２に示されるノッキング判定ルーチンにおいて、ノイズフラグがセットされているか否かが判別される。

ノイズフラグがセットされていないと判別されたときにはステップ２５４に進んで、進角量β１がβとされ、次いでステップ２５６に進む。一方、ノイズフラグがセットされていると判別されたときにはステップ２５５に進んで、進角量β２がβとされ、次いでステップ２５６に進む。ステップ２５６では、点火時期を進角させるために点火時期遅角量ΔＩＧからβが減算される。次いで、ステップ２５７では、点火時期遅角量ΔＩＧか負になったか否かが判別され、点火時期遅角量ΔＩＧか負になっていないときにはステップ２５９に進んで、次のサイクルにおける最終的な点火時期ＩＧＯが算出される。このときには、点火時期はβだけ進角される。一方、ステップ２５７において、点火時期遅角量ΔＩＧか負になったと判別されたときにはステップ２５８に進んで点火時期遅角量ΔＩＧが零とされた後、ステップ２５９に進む。このときには、点火時期が基準点火時期ＩＧとされる。

次に、図３４から図４２を参照しつつ、点火時期が遅角されたときのノッキング強度を代表する値の推定値ｙｅの予測値ｙｅｅを、第２のニューラルネットワークを用いて推定するようにした本発明による第２の実施例について説明する。この第２の実施例では、第２のニューラルネットワークとして、リカレントニューラルネットワークが用いられている。図３４は、第２の実施例で用いられているリカレントニューラルネットワーク（ＲＮＮ）を展開して表した図を示しており、図３５は、リカレントニューラルネットワークを示している。このリカレントニューラルネットワークは周知であり、従って、以下、簡単にリカレントニューラルネットワークについて説明する。

図３４において、ｘ_１ ^ｔ−１、ｘ_２ ^ｔ−１、ｘ_１ ^ｔ、ｘ_２ ^ｔ、ｘ_１ ^ｔ＋１、ｘ_２ ^ｔ＋１、ｘ_１ ^ｔ＋２、ｘ_２ ^ｔ＋２ は、夫々時刻ｔ−１、ｔ、ｔ＋１、ｔ＋２におけるリカレントニューラルネットワークの入力層への時系列入力値を示しており、ｙ^ｔ、ｙ^ｔ＋１、ｙ^ｔ＋２、ｙ^ｔ＋３ は、夫々時刻ｔ−１、ｔ、ｔ＋１、ｔ＋２におけるリカレントニューラルネットワークの出力層からの出力値を示している。また、ｈ^ｔ−１、ｈ^ｔ、ｈ^ｔ＋１、ｈ^ｔ＋２ （ｈはベクトル）は、夫々時刻ｔ−１、ｔ、ｔ＋１、ｔ＋２におけるリカレントニューラルネットワークの隠れ層からの出力値を表しており、これら、ｈ^ｔ−１、ｈ^ｔ、ｈ^ｔ＋１、ｈ^ｔ＋２ は、隠れ状態ベクトルと称される。

一方、図３５を参照すると、このリカレントニューラルネットワークにおいて、L＝１は入力層、L＝２は隠れ層、L＝３は出力層を夫々示している。図３５において鎖線枠内のノードは実際には存在せず、説明のために記載したノードであり、従って、図３５に示される例では、入力層 ( L＝１) が２個のノードからなる。また、図３５に示される例では、隠れ層 ( L＝１) がｍ個のノードからなる（図３５には３個のみが記載されている）。一方、図３５には一つの隠れ層 ( L＝２)しか記載されていないが、この隠れ層の層数は、任意の個数とすることができる。なお、出力層 ( L＝３) のノードの数は１個とされている。

図３５は、図３４の時刻ｔにおける入力値ｘ_１ ^ｔ、ｘ_２ ^ｔと、出力値ｙ^ｔ＋１とを示している。また、図３６において、ｈ_１ ^ｔ、ｈ_２ ^ｔ・・・ｈ_ｍ ^ｔ は、図３４の時刻ｔにおける隠れ層の各ノードからの出力値、即ち、隠れ状態ベクトルを示している。リカレントニューラルネットワークでは、図３５に示されるように、一つ前の時刻ｔ−１における隠れ層の各ノードからの出力値、即ち、隠れ状態ベクトルｈ_１ ^ｔ−１、ｈ_２ ^ｔ−１・・・ｈ_ｍ ^ｔ−１ が、夫々対応する重みｗを掛けて、隠れ層の各ノードに入力される。従って、隠れ層の各ノードには、各入力値ｘ_１ ^ｔ、ｘ_２ ^ｔに夫々対応する重みｗを掛けた値と、一つ前の時刻ｔ−１における各隠れ状態ベクトルｈ_１ ^ｔ−１、ｈ_２ ^ｔ−１・・・ｈ_ｋ ^ｔ−１に夫々対応する重みｗを掛けた値が入力され、図３５の隠れ層 ( L＝２) の各ノード（ｋ＝１，２・・・m）において算出される総入力値ｕ_ｋは、次式のようになる。

隠れ層の各ノードにおいて算出された総入力値ｕ_ｋは、活性化関数により変換されて、隠れ層の各ノードから隠れ状態ベクトルｈ_ｋ ^ｔ （ｋ＝１，２・・・m）として出力される。この場合、活性化関数として、例えば、tanh 関数(双曲線正接関数)が用いられた場合には、隠れ層の各ノードから出力される隠れ状態ベクトルｈ_ｋ ^ｔ は,ｈ_ｋ ^ｔ ＝tanh (ｕ_k) となる。これらの隠れ状態ベクトルｈ_ｋ ^ｔ は, 出力層 ( L＝３) のノード に入力され、出力層 のノードでは、夫々対応する重みｗを用いて、次式で示される総入力値ｕが算出される。

本発明による実施例では、出力層のノードでは、活性化関数として恒等関数が用いられており、従って、出力層のノードからは、出力層のノードにおいて算出された総入力値ｕが、そのまま出力値ｙとして出力される。

次に、図３４に示される入力値ｘ_１ ^ｔ−１、ｘ_２ ^ｔ−１、ｘ_１ ^ｔ、ｘ_２ ^ｔ、ｘ_１ ^ｔ＋１、ｘ_２ ^ｔ＋１、ｘ_１ ^ｔ＋２、ｘ_２ ^ｔ＋２ および出力値ｙ^ｔ、ｙ^ｔ＋１、ｙ^ｔ＋２、ｙ^ｔ＋３ について説明する。図３４において時刻ｔ−１、ｔ、ｔ＋１、ｔ＋２は、同一気筒の連続する各サイクルに対応しており、本発明による実施例では、時刻ｔ−１、ｔ、ｔ＋１、ｔ＋２において、即ち、同一気筒の連続する各サイクルにおいて第１の学習済みニューラルネットワーク２０から順次出力されたノッキング強度を代表する値の推定値ｙeが夫々、入力値ｘ_１ ^ｔ−１、ｘ_１ ^ｔ、ｘ_１ ^ｔ＋１、ｘ_１ ^ｔ＋２ とされる。また、時刻ｔ−１、ｔ、ｔ＋１、ｔ＋２における、即ち、同一気筒の連続する各サイクルにおける点火時期の遅角量α又は進角量βが夫々、入力値ｘ_２ ^ｔ−１、ｘ_２ ^ｔ、ｘ_２ ^ｔ＋１ｘ_２ ^ｔ＋２ とされる。

また、この場合、入力値として、ＥＧＲ率ＥＲを追加することもできる。一方、出力値ｙ^ｔ、ｙ^ｔ＋１、ｙ^ｔ＋２、ｙ^ｔ＋３ は夫々、時刻ｔ、ｔ＋１、ｔ＋２、ｔ＋３における、即ち、次のサイクルにおけるノッキング強度を代表する値の推定値ｙe の予測値ｙｅｅである。図３６には、代表的な例として、時刻ｔにおける入力値および出力値の一覧表が示されている。

図３７は、図３４に示されるリカレントニューラルネットワークの重みを学習するための訓練データセットを示している。この訓練データセットは、第２の実施例では、図２８に示される一覧表から、全サイクルに亘る一部のデータ、即ち、第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙe と点火時期の遅角量α又は進角量βのみを抜き出すことによって作成される。この場合、この訓練データセットでは、次のサイクルにおける第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙe が、ノッキング強度を代表する値の推定値ｙe の予測値ｙｅｅとされ、このノッキング強度を代表する値の推定値ｙe の予測値ｙｅｅが、教師データｙｔとされる。

例えば、図３７において、サイクルがサイクルｋ_ｎのときには、サイクルｋ_ｎのときに第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙ_eｎ と, サイクルｋ_ｎのときの点火時期の遅角量αが、サイクルｋ_ｎのときの入力値ｘ_１ ^ｔ、ｘ_２ ^ｔとされ、サイクルｋ_ｎ＋１のときに第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙ_eｎ＋１ が、サイクルｋ_ｎのときの教師データｙｔとされる。

この第２実施例でも、図３5に示されるリカレントニューラルネットワークの入力層 ( L＝１) のノード数、隠れ層 ( L＝２) のノード数および出力層 ( L＝3) のノード数、および、図３７に示される訓練データセットが、図５に示される学習装置２３のメモリ２４に記憶され、学習装置２３において、リカレントニューラルネットワークの重みの学習が行われる。このリカレントニューラルネットワークの重みの学習は、通常、展開表示されたリカレントニューラルネットワークの一部を切り取って誤差逆伝播を行なわせるTrancated BBPT(Backpropagation Through Time) 手法を用いて行われる。

例えば、図３４が、切り取られたリカレントニューラルネットワークの一部であるとしると、図３４において、時刻ｔ−１、ｔ、ｔ＋１、ｔ＋２における入力値ｘ_１ ^ｔ−１、ｘ_２ ^ｔ−１、ｘ_１ ^ｔ、ｘ_２ ^ｔ、ｘ_１ ^ｔ＋１、ｘ_２ ^ｔ＋１、ｘ_１ ^ｔ＋２、ｘ_２ ^ｔ＋２ を、リカレントニューラルネットワークに順次入力させ、時刻ｔ＋２においてリカレントニューラルネットワークから出力される出力値ｙ^ｔ＋３ と、対応する教師データｙｔとの二乗誤差Ｅ（＝1/2 (ｙ^ｔ＋３ -ｙｔ)²）が小さくなるように、誤差逆伝播法により、リカレントニューラルネットワークの重みの学習が行われる。このときの誤差逆伝播は時間を遡るように行われるが、詳細については説明を省略する。

第２実施例では、図３７において、例えば、連続する１０サイクルのデータに対し、誤差逆伝播法により、リカレントニューラルネットワークの重みの学習が行われ、連続する１０サイクルのデータに対するリカレントニューラルネットワークの重みの学習が完了すると、次の連続する１０サイクルのデータに対し、リカレントニューラルネットワークの重みの学習が行われる。このようにして全ての連続する１０サイクルのデータに対し、重みの学習が完了するまで、リカレントニューラルネットワークの重みの学習が行われる。

図３８は、重みの学習が完了した学習済みリカレントニューラルネットワークを示している。図３８において、時刻ｔを現在の時刻とした場合、この学習済みリカレントニューラルネットワークでは、過去の時刻ｔ−３、ｔ−２、ｔ−１における入力値ｘ_１ ^ｔ−３、ｘ_２ ^ｔ−３、ｘ_１ ^ｔ−２、ｘ_２ ^ｔ−２、ｘ_１ ^ｔ―１、ｘ_２ ^ｔ―１ と、現在の時刻ｔにおける入力値ｘ_１ ^ｔ、ｘ_２ ^ｔを、リカレントニューラルネットワークに順次入力すると、リカレントニューラルネットワークからは、現在の時刻ｔにおける。ノッキング強度を代表する値の推定値ｙｅの予測値ｙｅｅが出力される。即ち、この第２実施例では、学習済みリカレントニューラルネットワークを用いて、点火時期が遅角されたときのノッキング強度を代表する値の推定値ｙｅの予測値ｙｅｅの推定モデルが生成される。

本発明による第２の実施例では、第１の学習済みニューラルネットワーク２０により生成されたノッキング強度を代表する値の推定モデルと、学習済みのリカレントニューラルネットワークにより生成されたノッキング強度を代表する値の推定値ｙｅの予測値ｙｅｅの推定モデルを用いて、市販車両の機関におけるノッキング処理が行われる。そのためにこれらノッキング強度を代表する値の推定モデルおよびノッキング強度を代表する値の推定値ｙe の予測値ｙｅｅの推定モデル、即ち、第１の学習済みニューラルネットワーク２０および学習済みのリカレントニューラルネットワークが、市販車両の電子制御ユニット３０内に格納される。これら推定モデルの市販車両の電子制御ユニット３０への格納作用は、図１１を参照しつつ既に説明したやり方と同様な方法でもって、図１１に示される電子制御ユニットへのデータ読み込みルーチンを用いて行われる。

このように第１の学習済みニューラルネットワーク２０および学習済みのリカレントニューラルネットワークが、市販車両の電子制御ユニット３０に格納されると、電子制御ユニット３０のメモリ３２内には、第１の学習済みニューラルネットワーク２０および学習済みのリカレントニューラルネットワークが形成される。図３９は、市販車両の電子制御ユニット３０のメモリ３２内に形成された第１の学習済みニューラルネットワーク２０および学習済みのリカレントニューラルネットワークを用いて、市販車両の機関運転時に実行されるノッキング処理を示している。このノッキング処理は、気筒毎に、個別に、一サイクル毎に実行される。本発明による第２の実施例でも、このノッキング処理は、例えば、クランク角が圧縮上死点に達したときに開始される。

図３９に示されるノッキング処理は、図３０に示されるノッキング処理と同様である。即ち、図３９を参照すると、まず初めにステップ２２０において、圧力センサ１９の出力値から求められるノッキング強度を代表する値、即ち、ノッキング強度代表値の推定処理が行われる。次いで、ステップ２２１では、ノッキングが生じているか否かのノッキング判定処理が行われ、次いで、ステップ２２２では、点火制御が行われる。ステップ２２０において行われるノッキング強度代表値の推定処理は、図１３又は図１４のいずれか一方に示されるノッキング強度代表値の推定ルーチンにより行われる。これらのノッキング強度代表値の推定ルーチンについては既に説明したので、これらノッキング強度代表値の推定ルーチンについては説明を省略する。なお、これらのノッキング強度代表値の推定ルーチンが実行されると、第１の学習済みニューラルネットワーク２０からはノッキング強度を代表する値の推定値ｙe が出力される。

図４０および図４１は、図３９のステップ２２０において第１の学習済みニューラルネットワーク２０からノッキング強度を代表する値の推定値ｙe が出力されたときに、この第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙe に基づき、図３９のステップ２２１で行われるノッキング判定ルーチンを示している。図４０および図４１に示されるノッキング判定ルーチンのステップ２３０から２３２およびステップ２３５から２４４は、図３１および図３２に示されるノッキング判定ルーチンのステップ２３０から２３２およびステップ２３５から２４４と同一であり、図４０および図４１に示されるノッキング判定ルーチンにおいて図３１および図３２に示されるノッキング判定ルーチンと異なるところは、図４０および図４１に示されるノッキング判定ルーチンでは、ステップ２３０の前にステップ２３０ａが追加されており、図３１および図３２に示されるノッキング判定ルーチンのステップ２３３および２３４に代えて、ステップ２３３ａが設けられたことだけである。従って、図４０および図４１については、ステップ２３０ａおよびステップ２３３ａのみについて説明し、その他のステップについての説明は省略する。

図４０および図４１を参照すると、ステップ２３０ａでは、第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙe が、同一気筒に対し、予め定められた回数の点火が行われるまで、例えば、５回点火が行われるまで記憶される。一方、ステップ２３３ａでは、予め定められた回数前の点火が行われたサイクルから、例えば、５回前に点火が行われたサイクルから現在のサイクルに至るまでの各サイクルにおいて第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙe 、およびこれら各サイクルにおける点火時期の遅角量α或いは進角量βが、順次、リカレントニューラルネットワークの入力層の各ノードに入力され、それによりリカレントニューラルネットワークからは、現在のサイクルにおけるノッキング強度を代表する値の推定値ｙeの予測値ｙｅｅが出力される。

この第２実施例でも、ステップ２３５では、未学習の機関の振動が発生しているか否かを判別するために、ノッキング強度を代表する値の推定値ｙe の予測値ｙｅeと第１の学習済みニューラルネットワーク２０から出力されたノッキング強度を代表する値の推定値ｙe との差（ｙｅe ―ｙｅ）の絶対値が設定値Ｓよりも大きいか否かが判別され、差（ｙｅe ―ｙｅ）の絶対値が設定値Ｓよりも小さいと判別されたときには、ステップ２３６に進んで、ノイズフラグがリセットされる。次いで、ステップ２３９に進んで、ノッキング強度を代表する値の推定値ｙe が閾値Ｍij よりも大きいか否かが判別される。一方、差（ｙｅe ―ｙｅ）の絶対値が設定値Ｓよりも大きいと判別されたときには、ステップ２３５からステップ２４０に進んで、ノイズフラグがセットされる。次いで、ステップ２４１において、ノッキング発生フラグがリセットされ、次いで、ステップ２４２において、判定中止フラグがセットされる。

図４２は、図４０および図４１に示されるノッキング判定ルーチンにおける判定結果に基づき電子制御ユニット３０内で実行される点火制御ルーチンを示している。図４２に示される点火制御ルーチンのステップ２５０から２５９は、図３３に示される点火制御ルーチンのステップ２５０から２５９と同一であり、図４２に示される点火制御ルーチンにおいて図３３に示される点火制御ルーチンと異なるところは、図４２に示されるノッキング判定ルーチンでは、ステップ２５９の後にステップ２５９ａが追加されていることだけである。従って、図４２については、ステップ２５９ａのみについて説明し、その他のステップについての説明は省略する。

図４２を参照すると、ステップ２５９ａでは、点火時期の遅角量α或いは進角量βが、同一気筒に対し、予め定められた回数の点火が行われるまで、例えば、５回点火が行われるまで記憶される。この第２実施例でも、ノッキング発生フラグがセットされたときには、点火時期は一定量αだけ遅角され、ノッキング発生フラグがリセットされたときには、点火時期は一定量β１又はβ２αだけ進角される。

このように、本発明によれば、点火により発生する混合気の燃焼圧を検出するための圧力センサ１９の出力値から求められるノッキング強度を代表する値を、機関本体の振動を検出するノッキングセンサ１８の出力値から推定する第１の学習済みニューラルネットワークが記憶装置３２内に記憶されている。点火時期が遅角されたときのノッキング強度を代表する値の推定値の予測値又は予測低下量を推定する第２の学習済みニューラルネットワークが記憶装置３２内に記憶されている。第１の学習済みニューラルネットワークを用いて算出されたノッキング強度を代表する値の推定値が予め定められた閾値を越えると次のサイクルにおける点火時期が遅角され、点火時期が遅角されたこの次のサイクルでは、第２の学習済みニューラルネットワークを用いて算出されたノッキング強度を代表する値の推定値の予測値と第１の学習済みニューラルネットワークを用いて算出されたノッキング強度を代表する値の推定値との差に基づいて更に次のサイクルにおける点火時期の遅角制御が行われる。この遅角制御では、この差が予め定められた設定値よりも小さいときには、ノッキング強度を代表する値の推定値が予め設定された閾値よりも大きければ、更に次のサイクルにおける点火時期の遅角作用が行われ、この差が予め定められた設定値よりも大きいときには、たとえノッキング強度を代表する値の推定値が予め設定された閾値よりも大きかったとしても、更に次のサイクルにおける点火時期の遅角作用が行われない。

一方、本発明による第１の実施例では、第２の学習済みニューラルネットワークは、機関の運転状態、点火時期の遅角量、および一つ前のサイクルにおけるノッキング強度を代表する値の推定値を入力すると、点火時期が遅角されときのノッキング強度を代表する値の推定値の予測低下量を出力する。この予測低下量から点火時期が遅角されときのノッキング強度を代表する値の推定値の予測値が求められる。この場合、本発明による第１の実施例では、機関の運転状態が、機関回転数、機関負荷およびＥＧＲ率からなる。

また、本発明による第２の実施例では、第２の学習済みニューラルネットワークは、予め定められた回数前の点火が行われたサイクルから現在のサイクルに至るまでの各サイクルにおけるノッキング強度を代表する値の推定値および点火時期の遅角量又は進角量を入力すると、現在のサイクルにおけるノッキング強度を代表する値の推定値の予測値を出力するリカレントニューラルネットワークからなる。なお、この第２の実施例では、リカレントニューラルネットワークに代えて、ゲート付きリカレントニューラルネットワーク、例えば、ＬＳＴＭ（Long short-term memory ）を用いることもできる。

また、本発明による第１および第２の実施例では、点火時期の遅角作用が行われた後、ノッキング強度を代表する値の推定値が予め定められた閾値以下となったときに、点火時期の進角作用が開始され、点火時期の遅角作用が行われたときのノッキング強度を代表する推定値とノッキング強度を代表する推定値の予測値との差が設定値よりも大きかったときには、この差が設定値以下だったときに比べて点火時期の進角量が小さくされる。

更に、本発明による実施例では、第１の学習済みニューラルネットワークは、予め設定されている期間内においてノッキングセンサ１８により検出された機関本体１の振動を表す値を入力すると、ノッキング強度を代表する値の推定値を出力する。この場合、本発明による実施例ではノッキング強度を代表する値が、圧力センサの出力値のピーク値とされるか、或いは、ノッキング強度を代表する値が、圧力センサの出力値の積分値とされる。また、本発明による実施例では、機関本体の振動を表す値が、予め設定されている期間内におけるノッキングセンサの出力値とされるか、或いは、機関本体の振動を表す値が、予め設定されている期間内において等分割された区間内におけるノッキングセンサの出力値の積分値とされる。この場合、この予め設定されている期間が、一定のクランク角度範囲とされる。

さて、前述したように、ノッキングは通常、圧縮上死点から圧縮上死点後９０°までの間で発生する。従って、これまで述べた実施例では、圧縮上死点から圧縮上死点後９０°までの間に取り込まれたノッキングセンサ１８の出力値および圧力センサ１９の出力値に基づいて、ノッキングの発生を検出している。しかしながらノッキングによる機関本体１の振動は、機関回転数によらずに、圧縮上死点後、一定時間発生する傾向がある。このような傾向を考慮すると、圧縮上死点後、一定時間の間に取り込まれたノッキングセンサ１８の出力値および圧力センサ１９の出力値に基づいて、ノッキングの発生を検出することが考えられる。

次に、図４３および図４４を参照しつつ、圧縮上死点後、一定時間の間に取り込まれたノッキングセンサ１８の出力値および圧力センサ１９の出力値に基づいて、第１のニューラルネットワーク２０の重みの学習を行うと共に、圧縮上死点後、一定時間の間に取り込まれたノッキングセンサ１８の出力値に基づいてノッキングの発生を検出するようにした実施例について説明する。機関の作動の一サイクルのうちで、ノッキングセンサ１８および圧力センサ１９の双方の出力値、或いは、ノッキングセンサ１８の出力値が取り込まれるクランク角範囲はゲート区間と称され、図４３には、この実施例におけるゲート区間（縦軸）と、機関回転数（横軸）との関係が実線で示されている。なお、横軸におけるＮ１，Ｎ２，Ｎ３は代表的な機関回転数を示している。

一方、図４４には、各機関回転数Ｎ１，Ｎ２，Ｎ３に対するゲート区間Ｇが示されている。また、図４４には、等間隔で分割された各分割区間内におけるノッキングセンサ１８の出力値の積分値（負側の積分値も正の値とする）を算出するようにした場合における各分割区間Ｓが示されている。図４３および図４４から、この実施例では、機関回転数が高くなるほど、ゲート区間Ｇが長くなることがわかる。なお、図４３および図４４に示される例では、ノッキングセンサ１８の出力値および圧力センサ１９の出力値が取り込まれる期間が予め設定されており、この予め設定されている期間が、一定時間とされている。しかしながら、この予め設定されている期間は、必ずしも一定時間とする必要はなく、機関回転数が増大するにつれてゲート区間も増大するようにゲート区間を設定することもできる。

一方、これまで述べた実施例では、機関の全運転領域に対し、一つの第１の学習済みニューラルネットワークを用いて、ノッキング強度を代表する値を推定するようにしている。しかしながら、このように機関の全運転領域に対して一つの第１の学習済みニューラルネットワークしか用いないとノッキング強度を代表する値の推定精度が低下する危険性がある。このような危険性を回避するためには、機関の運転領域を複数個の運転領域に分割し、分割された運転領域毎に夫々別個のニューラルネットワークを用いることが好ましい。図４５Ａから図４５Ｃは、このように分割された運転領域毎に夫々別個の第１のニューラルネットワークを用いるようにした実施例を示している。

即ち、この実施例では、図４５Ａから図４５Ｃに示されるように、機関の運転領域が、機関負荷Ｌおよび機関回転数Ｎに応じて、複数個の運転領域に分割されており、図４５Ａに示されるように、分割された運転領域毎に夫々別個の第１のニューラルネットワークＮＮij が形成されている。この場合には、図４５Ｂに示されるように、分割された運転領域毎に夫々、図９に示されるような訓練データセットＤＳij が作成される。分割された運転領域毎に夫々形成された第１のニューラルネットワークＮＮij の重みの学習は、対応する訓練データセットＤＳij を用い、図１０に示される学習処理ルーチンを用いて行われる。また、この実施例では、図４５Ｃに示されるように、分割された運転領域毎に、ノッキング強度代表値に対する閾値Ｍij が設定されている。この実施例では、分割された運転領域毎に、対応する第１の学習済みニューラルネットワークＮＮij を用い、図１３又は図１４に示されるノッキング強度代表値の推定ルーチンを用いて、ノッキング強度代表値の推定値ｙe が算出される。

一方、この実施例では、第２のニューラルネットワークも、図４５Ａに示される分割された運転領域毎に夫々形成されている。また、第２のニューラルネットワークに対する訓練データセットも、分割された運転領域毎に夫々作成され、分割された運転領域毎に夫々形成された第２のニューラルネットワークの重みの学習は、対応する訓練データセットを用いて行われる。或る運転領域の第１の学習済みニューラルネットワークＮＮij により、ノッキング強度代表値の推定値ｙe が算出されると、同じ運転領域に対して作成されている第２の学習済みニューラルネットワークを用い、図４５Ｃに示される対応する閾値Ｍij を用いて、ノッキングが発生したか否かが判別される。

このようにこの実施例では、機関の運転領域が複数個の運転領域に分割されると共に、分割された各運転領域に対して夫々第１の学習済みニューラルネットワークと第２の学習済みニューラルネットワークが記憶装置内に記憶されている。分割された運転領域毎に、第２の学習済みニューラルネットワークを用いて算出されたノッキング強度を代表する値の推定値の予測値と第１の学習済みニューラルネットワークを用いて算出されたノッキング強度を代表する値の推定値との差に基づいて次のサイクルにおける点火時期の遅角制御が行われる。

図４６および図４８に、更に別の実施例を示す。この実施例では、図４６に示されるように、機関本体１にノッキングセンサ１８ａおよびノッキングセンサ１８ｂが取付けられており、これらノッキングセンサ１８ａおよびノッキングセンサ１８ｂは、ノッキングセンサ１８ａの出力値の波形およびノッキングセンサ１８ｂの出力値の波形を検出可能な検出器２１、例えば、オシロスコープに接続されている。図４６において、左から順に１番気筒＃１、２番気筒＃２、３番気筒＃３、４番気筒＃４と称すると、ノッキングセンサ１８ａは１番気筒＃１および２番気筒＃２から等距離の位置に配置されており、ノッキングセンサ１８ｂは３番気筒＃３および４番気筒＃４から等距離の位置に配置されている。図４６に示されるその他の構成は、図５に示される構成と同じであり、従って、この図４６に示されるその他の構成については、説明を省略する。

図４７は、この実施例において用いられる第１のニューラルネットワーク５５を示している。図４７を参照すると、この第１のニューラルネットワーク５５においても、図３および図８に示されるニューラルネットワークと同様に、L＝１は入力層、L＝２および L＝３は隠れ層、L＝４は出力層を夫々示している。この実施例では、図４７からわかるように、入力層 ( L＝１) が２ｎ個のノードからなる。一方、図４７には隠れ層 ( L＝２)および隠れ層 ( L＝３)が記載されているが、これら隠れ層の層数は、１個又は任意の個数とすることができ、またこれら隠れ層のノードの数も任意の個数とすることができる。なお、出力層 ( L＝４) のノードの数は１個とされており、出力層のノードからの出力値がｙで示されている。

図４７において、入力値ｘ_１、ｘ_２・・・ｘ_ｎ−１、ｘ_ｎはフィルタリング後のノッキングセンサ１８ａの出力値であり、入力値ｘｘ_１、ｘｘ_２・・・ｘｘ_ｎ−１、ｘｘ_ｎはフィルタリング後のノッキングセンサ１８ｂの出力値である。即ち、この実施例では、両方のノッキングセンサ１８ａおよびノッキングセンサ１８ｂのフィルタリング後の出力値が入力層 ( L＝１) の各ノードに入力される。なお、この実施例でも、入力値ｘ_１、ｘ_２・・・ｘｘ_ｎ−１、ｘｘ_ｎとして、図６Ｂに示される値か、或いは、図６Ｃに示される値が用いられる。即ち、図６Ｂに示されるように、フィルタリングされた後のノッキングセンサ１８ａおよびノッキングセンサ１８ｂの出力値そのものが入力値ｘ_１、ｘ_２・・・ｘｘ_ｎ−１、ｘｘ_ｎとされるか、或いは、図６Ｃに示されるように、フィルタリングされた後のノッキングセンサ１８ａおよびノッキングセンサ１８ｂの出力値の積分値（負側の積分値も正の値とする）であって例えば一定クランク角度内におけるノッキングセンサ１８ａおよびノッキングセンサ１８ｂの出力値の積分値が入力値ｘ_１、ｘ_２・・・ｘｘ_ｎ−１、ｘｘ_ｎとされる。

一方、図４７に示される出力値ｙとしては、図８に示される場合と同様に、圧力センサ１９の出力値から求められるノッキング強度を代表する値が採用されている。この場合、図８に示される場合と同様に、図７Ｂにおいて丸印で示されるフィルタリング後の圧力センサ１９の出力値のピーク値が出力値ｙとして用いられるか、或いは、図７Ｃに示されるフィルタリング後の圧力センサ１９の出力値の積分値（負側の積分値も正の値とする）が出力値ｙとして用いられる。なお、この実施例においても、圧力センサ１９の出力値から求められるノッキング強度を代表する値の実測値が教師データｙｔとされる。

図４８は、入力値ｘ_１、ｘ_２・・・ｘｘ_ｍ−１、ｘｘ_ｍと、圧力センサ１９の出力値から求められるノッキング強度を代表する値の実測値、即ち、教師データｙｔとを用いて作成された訓練データセットを示している。図４８に示されるように、この訓練データセットでは、入力値ｘ_１、ｘ_２・・・ｘｘ_ｎ−１、ｘｘ_ｎと教師データｙｔとの関係を表すｍ個のデータが取得されている。この訓練データセットの作成方法は、既に説明した図９に示される訓練データセットの作成方法と同じである。即ち、訓練データセットを作成する際には、図４６において、機関負荷および機関回転数の種々の組み合わせについてノッキングの発生しない運転状態とノッキングの発生する運転状態の双方の状態が生ずるように機関を運転させ、そのとき各検出器２１、２２から得られたデータに基づいて、図４８に示されるような訓練データセットが作成される。

訓練データセットが作成されると、作成された訓練データセットの電子データを用いて、図４７に示される第１のニューラルネットワーク５５の重みの学習が行われる。この第１のニューラルネットワークの重みの学習方法は、図５を参照しつつ前述した第１のニューラルネットワークの重みの学習方法と同じである。即ち、図４７に示される第１のニューラルネットワーク５５の重みの学習を行うときには、この第１のニューラルネットワーク５５のノード数、および作成された訓練データセットの電子データが、図４６に示される学習装置２３のメモリ２４に記憶され、ＣＰＵ２５において、図１０に示される学習処理ルーチンを用いて、第１のニューラルネットワーク５５の重みの学習が行われる。

このように、この実施例では、機関本体１の振動を検出する複数個のノッキングセンサ１８ａおよび１８ｂを具備している。第１の学習済みニューラルネットワーク５５は、複数個のノッキングセンサ１８ａおよび１８ｂにより検出された機関本体１の振動を表す値を入力すると、ノッキング強度を代表する値を出力する。このように複数個のノッキングセンサ１８ａおよび１８ｂを用い、両方のノッキングセンサ１８ａおよび１８ｂの出力値を用いて、第１のニューラルネットワーク５５の重みの学習が行われると、一つのノッキングセンサを用いたときよりも、取得できる情報量が増大するために、ノッキング強度を代表する値を精度よく推定することができる。

１ 機関本体
５ 燃焼室
１１ 点火栓
１８、１８ａ、１８ｂ ノッキングセンサ
１９ 圧力センサ
２０，５５ 第１のニューラルネットワーク
２３ 学習装置
３０ 電子制御ユニット
５０ 第２のニューラルネットワーク

Claims (14)

1. 点火により発生する混合気の燃焼圧を検出するための圧力センサの出力値から求められるノッキング強度を代表する値を、機関本体の振動を検出するノッキングセンサの出力値から推定する第１の学習済みニューラルネットワークが記憶装置内に記憶されており、点火時期が遅角されたときの該ノッキング強度を代表する値の推定値の予測値又は予測低下量を推定する第２の学習済みニューラルネットワークが記憶装置内に記憶されており、該第１の学習済みニューラルネットワークを用いて算出されたノッキング強度を代表する値の推定値が予め定められた閾値を越えると次のサイクルにおける点火時期が遅角され、点火時期が遅角された該次のサイクルでは、第２の学習済みニューラルネットワークを用いて算出された該ノッキング強度を代表する値の推定値の予測値と第１の学習済みニューラルネットワークを用いて算出された該ノッキング強度を代表する値の推定値との差に基づいて更に次のサイクルにおける点火時期の遅角制御が行われ、この遅角制御では、該差が予め定められた設定値よりも小さいときには、該ノッキング強度を代表する値の推定値が予め設定された閾値よりも大きければ、該更に次のサイクルにおける点火時期の遅角作用が行われ、該差が予め定められた設定値よりも大きいときには、たとえ該ノッキング強度を代表する値の推定値が予め設定された閾値よりも大きかったとしても、該更に次のサイクルにおける点火時期の遅角作用が行われない内燃機関の点火時期制御装置。
2. 該第２の学習済みニューラルネットワークは、機関の運転状態、点火時期の遅角量、および一つ前のサイクルにおけるノッキング強度を代表する値の推定値を入力すると、点火時期が遅角されときのノッキング強度を代表する値の推定値の予測低下量を出力し、該予測低下量から点火時期が遅角されたときのノッキング強度を代表する値の推定値の予測値が求められる請求項１に記載の内燃機関の点火時期制御装置。
3. 該機関の運転状態が、機関回転数、機関負荷およびＥＧＲ率からなる請求項２に記載の内燃機関の点火時期制御装置。
4. 該第２の学習済みニューラルネットワークは、予め定められた回数前の点火が行われたサイクルから現在のサイクルに至るまでの各サイクルにおけるノッキング強度を代表する値の推定値および点火時期の遅角量又は進角量を入力すると、現在のサイクルにおけるノッキング強度を代表する値の推定値の予測値を出力するリカレントニューラルネットワークからなる請求項１に記載の内燃機関の点火時期制御装置。
5. 点火時期の遅角作用が行われた後、該ノッキング強度を代表する値の推定値が予め定められた閾値以下となったときに、点火時期の進角作用が開始され、点火時期の遅角作用が行われたときの上記差が設定値よりも大きかったときには該差が設定値以下だったときに比べて点火時期の進角量が小さくされる請求項１に記載の内燃機関の点火時期制御装置。
6. 該第１の学習済みニューラルネットワークは、予め設定されている期間内においてノッキングセンサにより検出された機関本体の振動を表す値を入力すると、該ノッキング強度を代表する値の推定値を出力する請求項１に記載の内燃機関の点火時期制御装置。
7. 該ノッキング強度を代表する値が、圧力センサの出力値のピーク値である請求項６に記載の内燃機関の点火時期制御装置。
8. 該ノッキング強度を代表する値が、圧力センサの出力値の積分値である請求項６に記載の内燃機関の点火時期制御装置。
9. 機関本体の振動を表す値が、該予め設定されている期間内におけるノッキングセンサの出力値である請求項６に記載の内燃機関の点火時期制御装置。
10. 機関本体の振動を表す値が、該予め設定されている期間内において等分割された区間内におけるノッキングセンサの出力値の積分値である請求項６に記載の内燃機関のノッキング検出装置。
11. 該予め設定されている期間が、一定のクランク角度範囲である請求項６に記載の内燃機関の点火時期制御装置。
12. 該予め設定されている期間が、一定時間である請求項６に記載の内燃機関の点火時期制御装置。
13. 機関の運転領域が複数個の運転領域に分割されると共に、分割された各運転領域に対して夫々第１の学習済みニューラルネットワークと第２の学習済みニューラルネットワークが記憶装置内に記憶されており、分割された運転領域毎に、第２の学習済みニューラルネットワークを用いて算出された該ノッキング強度を代表する値の推定値の予測値と第１の学習済みニューラルネットワークを用いて算出された該ノッキング強度を代表する値の推定値との差に基づいて上記更に次のサイクルにおける点火時期の遅角制御が行われる請求項１に記載の内燃機関の点火時期制御装置。
14. 機関本体の振動を検出する複数個のノッキングセンサを具備しており、第１の学習済みニューラルネットワークは、該複数個のノッキングセンサにより検出された機関本体の振動を表す値を入力すると、ノッキング強度を代表する値を出力する請求項１に記載の内燃機関の点火時期制御装置。

Images