(57)【要約】【課題】スタータ３２によってクランク軸２６に初期回転が付与される（クランキングが行われる）状況下、ピストン２４が圧縮上死点に到達する前に混合気の燃焼が開始される場合、クランク軸２６が逆回転する現象であるいわゆるケッチンが発生するおそれがあること。【解決手段】エンジン回転速度、スロットル開度、油温及び圧縮比を説明変数としたロジスティック回帰方程式に基づき、クランキング時においてケッチンが発生する確率であるケッチン発生率を予測する。そして、予測されたケッチン発生率が第１の閾値以上であって且つ第１の閾値よりも高い値に設定された第２の閾値未満であると判断された場合、点火プラグ２２の点火タイミングを圧縮上死点以降に遅角する処理を行う。一方、予測されたケッチン発生率が第２の閾値以上であると判断された場合、点火プラグ２２の点火を禁止する処理を行う。【選択図】      図１
2012002204.tif 000002
【特許請求の範囲】
【請求項１】
内燃機関の動作状態を制御するためのアクチュエータを備える内燃機関に適用され、
前記内燃機関の動作状態に影響を及ぼす複数のパラメータの値を検出する検出手段と、
該検出手段によって検出された複数のパラメータの値を説明変数とした重回帰方程式に基づき、前記内燃機関の動作状態を目的変数として推定する処理を行う推定手段と、
該推定手段によって推定された前記動作状態に基づき、前記アクチュエータを操作する処理を行う操作手段とを備えることを特徴とする内燃機関の制御装置。
【請求項２】
前記内燃機関の動作状態とは、該内燃機関の動作に関する異常発生率であり、
前記推定手段は、前記推定する処理として、前記異常発生率を目的変数として予測する処理を行うものであり、
前記操作手段は、前記予測された前記異常発生率が高いと判断された場合、前記操作する処理として、前記内燃機関の動作に関する異常を回避すべく前記アクチュエータの操作状態を変更する処理を行うことを特徴とする請求項１記載の内燃機関の制御装置。
【請求項３】
前記アクチュエータには、前記内燃機関の燃焼室内に突出する点火プラグが含まれ、
前記複数のパラメータの値には、機関回転速度と、前記燃焼室に供給される吸気量又はこれと相関を有するパラメータの値とが含まれ、
前記内燃機関の動作に関する異常発生率とは、ケッチン発生率であり、
前記操作手段は、前記予測されたケッチン発生率が高いと判断された場合、前記操作状態を変更する処理として、前記点火プラグの点火タイミングを遅角する処理及び該点火プラグの点火を禁止する処理のうち少なくとも１つを行うことを特徴とする請求項２記載の内燃機関の制御装置。
【請求項４】
前記操作手段は、前記予測されたケッチン発生率が高いと判断された場合、前記操作状態を変更する処理として、前記点火プラグの点火タイミングを圧縮上死点以降まで遅角する処理を行うことを特徴とする請求項３記載の内燃機関の制御装置。
【請求項５】
前記検出手段は、前記内燃機関の圧縮上死点となるタイミングを基準とした１燃焼サイクルの間に、前記複数のパラメータのそれぞれを少なくとも１回検出することを特徴とする請求項３又は４記載の内燃機関の制御装置。
【請求項６】
前記検出手段は、圧縮上死点直前における前記機関回転速度を検出することを特徴とする請求項５記載の内燃機関の制御装置。
【請求項７】
前記重回帰方程式とは、ロジスティック回帰方程式であることを特徴とする請求項１〜６のいずれか１項に記載の内燃機関の制御装置。
【請求項８】
前記説明変数としての前記複数のパラメータは、前記重回帰方程式に基づく前記内燃機関の動作状態の推定精度を表す自由度二重調整寄与率が最大となるように選択されていることを特徴とする請求項１〜７のいずれか１項に記載の内燃機関の制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、内燃機関の動作状態を制御するためのアクチュエータを備える内燃機関に適用される内燃機関の制御装置に関する。
【背景技術】
【０００２】
内燃機関の圧縮上死点に到達する前に内燃機関の燃焼室で吸気と燃料との混合気の燃焼が開始されると、例えば内燃機関の出力軸（クランク軸）に付与されるスタータ等からの初期回転力が不足する等、クランク軸を正回転させる方向のエネルギよりもクランク軸を逆回転させる方向のエネルギが大きくなることに起因して、クランク軸が逆回転する現象であるいわゆるケッチンが発生することがある。この場合、クランク軸と機械的に連結されたスタータ等の信頼性が低下するなどのおそれがある。
【０００３】
そこで従来、下記特許文献１に見られるように、内燃機関の駆動中において機関回転速度がアイドリング回転速度よりも低い所定回転速度未満になると判断された場合に、内燃機関の点火装置による点火を禁止する技術も提案されている。これにより、ケッチンが発生することによる不都合の抑制を図っている。なお、ケッチンの発生による不都合を抑制する技術としては、例えば下記特許文献２に記載されているものもある。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００８−２１５０８１号公報
【特許文献２】特許第３９２８８３９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、上記特許文献１に記載の技術によれば、ケッチンが発生したと把握された後に点火装置の点火を禁止するため、ケッチンが発生してからその後混合気の燃焼が開始されるまでの期間が長くなり、内燃機関の始動性が低下するおそれがある。そしてこの場合、その後の内燃機関の燃焼制御を速やかに開始することができなくなる等、内燃機関の動作状態を適切に制御することができなくなるおそれがある。
【０００６】
なお、上記に限らず、内燃機関の動作状態を制御するためのアクチュエータが備えられる内燃機関については、内燃機関の動作状態を適切に把握することができず、アクチュエータの操作によって内燃機関の動作状態を適切に制御することが困難となるおそれのあるこうした事情も概ね共通したものとなっている。
【０００７】
本発明は、上記課題を解決するためになされたものであり、その目的は、内燃機関の動作状態を適切に制御することのできる内燃機関の制御装置を提供することにある。
【課題を解決するための手段】
【０００８】
以下、上記課題を解決するための手段、及びその作用効果について記載する。
【０００９】
請求項１記載の発明は、内燃機関の動作状態を制御するためのアクチュエータを備える内燃機関に適用され、前記内燃機関の動作状態に影響を及ぼす複数のパラメータの値を検出する検出手段と、該検出手段によって検出された複数のパラメータの値を説明変数とした重回帰方程式に基づき、前記内燃機関の動作状態を目的変数として推定する処理を行う推定手段と、該推定手段によって推定された前記動作状態に基づき、前記アクチュエータを操作する処理を行う操作手段とを備えることを特徴とする。
【００１０】
上記発明では、内燃機関の動作状態に影響を及ぼす複数のパラメータの値を説明変数とした重回帰方程式に基づき内燃機関の動作状態を目的変数として推定することができる。そして推定された内燃機関の動作状態に基づき上記アクチュエータを操作する処理を行うため、内燃機関の動作状態を適切に把握しつつ内燃機関の動作状態を制御することができる。これにより、内燃機関の動作状態を適切に制御することができる。
【００１１】
請求項２記載の発明は、請求項１記載の発明において、前記内燃機関の動作状態とは、該内燃機関の動作に関する異常発生率であり、前記推定手段は、前記推定する処理として、前記異常発生率を目的変数として予測する処理を行うものであり、前記操作手段は、前記予測された前記異常発生率が高いと判断された場合、前記操作する処理として、前記内燃機関の動作に関する異常を回避すべく前記アクチュエータの操作状態を変更する処理を行うことを特徴とする。
【００１２】
上記発明では、内燃機関の動作に関する異常発生率に影響を及ぼす複数のパラメータの値を説明変数とした重回帰方程式に基づき、上記異常発生率を予測するため、上記異常が発生するおそれのある状況を予め把握することができる。そして予測された異常発生率が高いと判断された場合、上記態様にてアクチュエータの操作状態を変更する処理を行う。これにより、内燃機関の動作に関する異常が発生することを適切に抑制することができる。しかも、内燃機関の動作に関する異常を検出するためのセンサを備えることなく上記異常が発生するおそれのある状況を把握することができるため、上記異常を回避するための制御に要するセンサ等の部品数の増大を抑制することなどもできる。
【００１３】
請求項３記載の発明は、請求項２記載の発明において、前記アクチュエータには、前記内燃機関の燃焼室内に突出する点火プラグが含まれ、前記複数のパラメータの値には、機関回転速度と、前記燃焼室に供給される吸気量又はこれと相関を有するパラメータの値とが含まれ、前記内燃機関の動作に関する異常発生率とは、ケッチン発生率であり、前記操作手段は、前記予測されたケッチン発生率が高いと判断された場合、前記操作状態を変更する処理として、前記点火プラグの点火タイミングを遅角する処理及び該点火プラグの点火を禁止する処理のうち少なくとも１つを行うことを特徴とする。
【００１４】
機関回転速度が低いほど、内燃機関の出力軸（クランク軸）の回転エネルギが小さくなる。ここでクランク軸に初期回転が付与される（クランキングが行われる）状況下において、内燃機関の圧縮上死点に到達する前に燃焼室で混合気の燃焼が開始される場合、機関回転速度が低いほど、クランク軸を正回転させる方向の上記回転エネルギが、燃焼に伴い発生するエネルギのうちクランク軸を逆回転させる方向のエネルギよりも小さくなりやすくなり、ケッチンが発生しやすくなる。また、燃焼室に供給される吸気量が多いほど、燃焼室で混合気が圧縮されることに伴いクランク軸に作用するクランク軸の正回転を妨げる方向の回転力が大きくなる。ここでクランキングが行われる状況下において、上記圧縮上死点に到達する前に燃焼室で混合気の燃焼が開始される場合、吸気量が多いほど、クランク軸に作用するクランク軸の正回転を妨げる方向の回転力が大きくなりやすくなり、ケッチンが発生しやすくなる。
【００１５】
この点に鑑み、上記発明では、機関回転速度と、吸気量又はこれと相関を有するパラメータの値とを少なくとも説明変数とした重回帰方程式に基づき、クランキング時においてケッチンが発生する確率であるケッチン発生率を予測する。そして予測されたケッチン発生率が高いと判断された場合、点火プラグの操作状態を変更する処理として、点火プラグの点火タイミングを遅角する処理及び点火プラグの点火を禁止する処理のうち少なくとも１つを行う。点火タイミングを遅角する処理によれば、混合気の燃焼に伴いクランク軸に作用する回転力が大きくなる時期を圧縮上死点近傍に近づけたり、圧縮上死点以降としたりすることができ、混合気の燃焼に伴い発生するエネルギのうちクランク軸を逆回転させる方向のエネルギを極力小さくすることができる。また、点火プラグの点火を禁止する処理によれば、混合気の燃焼が行われず、混合気の燃焼に伴いクランク軸に作用するクランク軸を逆回転させる方向の回転力を大きく低下させることができる。このように、上記発明では、ケッチンが発生する蓋然性の高い状況下において上記態様にて点火プラグの操作状態を変更するため、内燃機関の始動性の低下と、ケッチンの発生との双方を好適に抑制することができる。
【００１６】
請求項４記載の発明は、請求項３記載の発明において、前記操作手段は、前記予測されたケッチン発生率が高いと判断された場合、前記操作状態を変更する処理として、前記点火プラグの点火タイミングを圧縮上死点以降まで遅角する処理を行うことを特徴とする。
【００１７】
上記発明では、上記態様にて点火プラグの操作状態を変更する処理を行うことで、混合気の燃焼に伴い発生するエネルギをクランク軸を正回転させる方向のエネルギとして適切に用いることができる。これにより、ケッチンの発生を好適に回避することができる。
【００１８】
請求項５記載の発明は、請求項３又は４記載の発明において、前記検出手段は、前記内燃機関の圧縮上死点となるタイミングを基準とした１燃焼サイクルの間に、前記複数のパラメータのそれぞれを少なくとも１回検出することを特徴とする。
【００１９】
上記発明では、ケッチンが内燃機関の圧縮上死点直前に燃焼が開始されることに起因して発生することに鑑み、ケッチンが発生すると想定されるタイミングまでにケッチン発生率の予測に用いる上記複数のパラメータとして極力新しい値を検出することができる。このため、ケッチン発生率を適切に予測することができる。
【００２０】
請求項６記載の発明は、請求項５記載の発明において、前記検出手段は、圧縮上死点直前における前記機関回転速度を検出することを特徴とする。
【００２１】
圧縮上死点直前における機関回転速度は、ケッチン発生率と高い相関を有する。この点に鑑み、上記発明では、圧縮上死点直前における機関回転速度の検出値を用いてケッチン発生率を予測するため、ケッチンの発生率の予測精度を向上させることができる。
【００２２】
請求項７記載の発明は、請求項１〜６のいずれか１項に記載の発明において、前記重回帰方程式とは、ロジスティック回帰方程式であることを特徴とする。
【００２３】
上記発明では、重回帰方程式としてロジスティック回帰方程式を用いることで、内燃機関の動作状態を適切に推定することができる。
【００２４】
請求項８記載の発明は、請求項１〜７のいずれか１項に記載の発明において、前記説明変数としての前記複数のパラメータは、前記重回帰方程式に基づく前記内燃機関の動作状態の推定精度を表す自由度二重調整寄与率が最大となるように選択されていることを特徴とする。
【００２５】
上記発明では、上記態様にて説明変数としての複数のパラメータが選択されているため、内燃機関の動作状態の推定精度を更に向上させることができる。
【図面の簡単な説明】
【００２６】
【図１】一実施形態にかかるシステム構成図。
【図２】一実施形態にかかるケッチンの発生に影響を及ぼす因子を列記した特性要因図。
【図３】一実施形態にかかる説明変数として用いる因子を示す図。
【図４】一実施形態にかかる自由度二重調整寄与率を極大とする因子を示す図。
【図５】一実施形態にかかるケッチン発生回避処理の手順を示すフローチャート。
【図６】一実施形態にかかるケッチン発生回避処理の一例を示すタイムチャート。
【発明を実施するための形態】
【００２７】
以下、本発明にかかる制御装置を空冷式単気筒エンジンが搭載された車両（四輪バギー：All Terrain Vehicle）に適用した一実施形態について、図面を参照しつつ説明する。
【００２８】
図１に本実施形態にかかるエンジンシステムの全体構成を示す。
【００２９】
図示されるエンジン１０は、４ストロークエンジンである。すなわち、このエンジン１０は、吸気・圧縮・膨張・排気の４行程を１燃焼サイクルとして運転される。
【００３０】
エンジン１０の吸気通路１２には、上流側から順に、スロットルバルブ１４と、吸気圧を検出する吸気圧センサ１６とが設けられている。スロットルバルブ１４は、その開度（スロットル開度）が調節されることで、エンジン１０の燃焼室１８へと供給される空気量（吸気量）を調節するためのものである。詳しくは、スロットル開度は、ユーザが車両の要求トルクを指示する図示しないアクセルグリップの操作に応じて調節される。なお、スロットルバルブ１４付近には、スロットル開度を検出するスロットルセンサ１４ａが内蔵されている。
【００３１】
吸気通路１２のうち、吸気圧センサ１６の下流側の吸気ポート近傍には、燃料を噴射供給する電磁駆動式の燃料噴射弁１９が設けられている。燃料噴射弁１９から噴射供給された燃料と吸気との混合気は、吸気バルブ２０の開動作によって燃焼室１８に供給される。
【００３２】
燃焼室１８に供給された混合気は、燃焼室１８内に突出する点火プラグ２２の放電火花によって着火され、燃焼に供される。混合気の燃焼によって発生するエネルギは、ピストン２４を介してエンジン１０の出力軸（クランク軸２６）の回転エネルギとして取り出される。なお、燃焼に供された混合気は、排気バルブ２８の開動作によって、排気として排気通路３０に排出される。
【００３３】
上記クランク軸２６には、スタータ３２が接続されている。スタータ３２は、図示しないバッテリからの給電を受けて図示しないイグニッションスイッチのオンによって始動し、エンジン１０を始動させるべくクランク軸２６に初期回転を付与する（クランキングを行う）。
【００３４】
クランク軸２６の回転角度（クランク回転角度）は、クランク軸２６付近に設けられたクランク角度センサ３４によって検出される。詳しくは、クランク角度センサ３４は、クランク軸２６と一体で回転するロータ２６ａの外周部に等間隔（例えば３０°ＣＡ毎）で複数設けられた突起がこのセンサを横切るときに、矩形状のクランク信号を出力する。ここでロータ２６ａは、突起が配置されない欠け歯部を有している。このため、クランク角度センサ３４は、この欠け歯部に対応するクランク回転角度においてクランク信号を出力しない。なお本実施形態では、ロータ２６ａは、ピストン２４が上死点直前に位置する時のクランク回転角度に対応する位置に欠け歯部を有している。このため、クランク角度センサ３４からクランク信号が出力されない時期は、吸気行程又は膨張行程のいずれかとなる。
【００３５】
上記クランク角度センサ３４や、スロットルセンサ１４ａ、吸気圧センサ１６、更にはエンジンオイルの温度（油温）を検出する油温センサ３６等の出力信号は、電子制御装置（以下、ＥＣＵ３８）に入力される。ＥＣＵ３８は、周知のＣＰＵ、ＲＯＭ、ＲＡＭ等よりなるマイクロコンピュータを主体として構成されている。ＥＣＵ３８は、上記各センサからの入力信号に基づき、ＲＯＭに記憶された各種の制御プログラムを実行することで、スタータ３２等による始動制御や、燃料噴射弁１９による燃料噴射制御、更には点火プラグ２２による点火制御等を行う。
【００３６】
上記始動制御は、クランキングの開始に伴い行程判別を行うとともに、燃料噴射制御及び点火制御によって燃焼室１８内の混合気の燃焼を開始させるための制御である。ここで行程判別とは、ピストン２４が規定位置（例えば、吸気行程における下死点）となる時のクランク回転角度を基準（０°ＣＡ）とした場合に、エンジン１０の１燃焼サイクル（７２０°ＣＡ）に対する現在のクランク回転角度を把握したり、クランク回転角度に対応する各行程（吸気・圧縮・膨張・排気）を把握したりすることである。詳しくはまず、クランク信号に基づき、クランキングが開始されてから欠け歯部が最初に検出されるタイミング付近の吸気圧と、欠け歯部が次に検出されるタイミング付近の吸気圧とを吸気圧センサ１６によって検出する。そして、これら検出された吸気圧のうち、低い方が検出されたタイミングを吸気行程として判別し、高い方が検出されたタイミングを膨張行程として判別する。これにより、行程判別を行うことが可能となる。
【００３７】
また、上記燃料噴射制御は、スロットルセンサ１４ａの出力値から算出されるスロットル開度に基づき、エンジン回転速度及び吸気圧と関係付けられた基本噴射量が規定されるマップ（ＤＪマップ）と、エンジン回転速度及びスロットル開度と関係付けられた基本噴射量が規定されるマップ（αＮマップ）とのうちいずれかを選択して用いて基本噴射量を算出し、算出された基本噴射量に基づき燃料噴射弁１９から燃料噴射させるための制御である。具体的には例えば、スロットル開度が規定開度以下となる場合、ＤＪマップを選択して用い、スロットル開度が上記規定開度を上回る場合には、αＮマップを選択して用いる。ここでＤＪマップとαＮマップとを使い分けるのは、燃焼室１８に供給される吸気量に対する基本噴射量の算出精度を向上させるためである。つまり、本実施形態では、吸気量を直接検出するセンサ（例えばエアフローメータ）が備えられていないため、基本噴射量を、吸気圧又はスロットル開度と、エンジン回転速度とに基づき定められる吸気量に応じて算出することとしている。ここでスロットル開度が小さい領域においては、スロットル開度の変化量に対する吸気量の変化量よりも、吸気圧の変化量に対する吸気量の変化量が大きくなるため、吸気圧及びエンジン回転速度に基づく吸気量の算出精度が、スロットル開度及びエンジン回転速度に基づく吸気量の算出精度よりも高くなる。一方、スロットル開度が大きい領域においては、吸気圧の変化量に対する吸気量の変化量よりも、スロットル開度の変化量に対する吸気量の変化量の方が大きいため、スロットル開度及びエンジン回転速度に基づく吸気量の算出精度が、吸気圧及びエンジン回転速度に基づく吸気量の算出精度よりも高くなる。このため、上記態様にて規定されたＤＪマップとαＮマップとをスロットル開度に応じて使い分けることで、吸気量に対する基本噴射量の算出精度を向上させることが可能となる。
【００３８】
ところで、クランキング時において、クランク軸２６を正回転させる方向にスタータ３２からクランク軸２６に付与される初期回転力が不足したり、クランク軸２６に作用するクランク軸２６の正回転を妨げる方向の回転力が大きくなったりする状況下、上記始動制御によってピストン２４が圧縮上死点に到達する前に混合気の燃焼が開始される場合、クランク軸２６が逆回転する現象であるいわゆるケッチンが発生することがある。この場合、クランク軸２６とスタータ３２との間の動力伝達経路に大きな力が作用することに起因して、動力伝達経路上に備えられるスタータ３２等の車載機器の信頼性が低下するなどの不都合が生じるおそれがある。ここでケッチンの発生による不都合を回避する手法としては例えば、クランク軸２６の回転方向を検出可能なセンサ（例えばロータリエンコーダ）を備え、このセンサの出力値に基づきクランク軸２６が逆回転したと判断される場合、点火プラグ２２による点火を禁止することも考えられる。しかしながら、上記手法を採用する場合には、部品数の増大によって車両の重量やコストが増大するおそれがある。特に、部品数の低減によって車両の重量やコストの増大を抑制するとの設計コンセプトを採用する車両においては、クランク軸２６の回転方向を検出可能なセンサを設けることは困難となる。
【００３９】
こうした問題を解決すべく、本実施形態では、クランキング時においてケッチンが発生する確率であるケッチン発生率に影響を及ぼす複数のパラメータの値を説明変数としたロジスティック回帰方程式に基づきケッチン発生率を目的変数として予測し、予測されたケッチン発生率に基づき点火制御処理を変更する処理（ケッチン発生回避処理）を行う。これにより、ケッチンの発生の回避を図る。以下、図２〜図４を用いて、ケッチン発生回避処理のうちケッチン発生率の予測手法について詳述する。
【００４０】
図２に、ケッチンの発生に影響を及ぼすと考えられる上記パラメータ（因子）を列記した特性要因図を示す。
【００４１】
図示されるように、本発明者はまず、ロジスティック回帰方程式で説明変数として用いる上記パラメータを決定すべく、ケッチンの発生に影響を及ぼすと考えられるパラメータを列記した特性要因図を作成した。そして、ケッチン発生率を予測するためのロジスティック回帰方程式で使用する複数のパラメータとして、特性要因図で列記されたパラメータの中から以下（Ａ）〜（Ｄ）のパラメータを説明変数として選択することとした。
【００４２】
（Ａ）エンジン回転速度ＮＥ：エンジン回転速度ＮＥが低いほど、クランク軸２６の回転エネルギが小さくなる。ここでクランキングが行われる状況下において、ピストン２４が圧縮上死点に到達する前に燃焼室１８で混合気の燃焼が開始される場合、エンジン回転速度ＮＥが低いほど、クランク軸２６を正回転させる方向の上記回転エネルギが、燃焼に伴い発生するエネルギのうちクランク軸２６を逆回転させる方向のエネルギよりも小さくなりやすくなり、ケッチン発生率が高くなる。このため、エンジン回転速度ＮＥは、ケッチン発生率に大きく影響を及ぼすパラメータとなる。
【００４３】
（Ｂ）スロットル開度ＶＴＡ：スロットル開度ＶＴＡが大きくなるほど、燃焼室１８に供給される吸気量が多くなる。そして吸気量が多くなるほど、燃焼室１８で混合気が圧縮されることに伴いクランク軸２６に作用するクランク軸２６の正回転を妨げる方向の回転力が大きくなる。ここでクランキングが行われる状況下において、ピストン２４が圧縮上死点に到達する前に燃焼室１８で混合気の燃焼が開始される場合、スロットル開度が大きいほど、吸気量が多くなり、クランキングによってクランク軸２６に作用するクランク軸２６を正回転させる方向の回転力よりも、クランク軸２６に作用するクランク軸２６の正回転を妨げる方向の回転力が大きくなりやすくなり、ケッチン発生率が高くなる。このため、スロットル開度ＶＴＡは、ケッチン発生率に大きく影響を及ぼすパラメータとなる。
【００４４】
なお本実施形態において、吸気量と相関を有するパラメータとして吸気圧ではなくスロットル開度ＶＴＡを用いているのは、本発明者の実施した実験等から、ケッチン発生率とスロットル開度との相関がケッチン発生率と吸気圧との相関よりも高くなるとの知見が得られたことによるものである。つまり、ケッチンが発生しやすい状況である吸気量が多くなる状況は、上述したαＮマップが選択される状況であり、スロットル開度と吸気量との相関が吸気圧と吸気量との相関よりも高くなる状況となる。
【００４５】
（Ｃ）油温ＴＨＯ：油温ＴＨＯが低いほど、エンジンオイルの粘度が高くなり、クランク軸２６に作用するクランク軸２６の正回転を妨げる方向の回転力が大きくなる。ここでクランキングが行われる状況下において、ピストン２４が圧縮上死点に到達する前に燃焼室１８で混合気の燃焼が開始される場合、油温ＴＨＯが低いほど、クランク軸２６に作用するクランク軸２６の正回転を妨げる方向の回転力が大きくなることで、ケッチン発生率が高くなる。このため、油温ＴＨＯは、ケッチン発生率に大きな影響を及ぼすパラメータとなる。
【００４６】
（Ｄ）圧縮比Ｒｃ：圧縮比Ｒｃが高いほど、燃焼室１８での混合気の圧縮に伴いクランク軸２６に作用するクランク軸２６の正回転を妨げる方向の回転力が大きくなる。ここでクランキングが行われる状況下において、ピストン２４が圧縮上死点に到達する前に燃焼室１８で混合気の燃焼が開始される場合、圧縮比Ｒｃが高いほど、クランク軸２６に作用するクランク軸２６の正回転を妨げる方向の回転力が大きくなることで、ケッチン発生率が高くなる。このため、圧縮比Ｒｃは、ケッチン発生率に大きな影響を及ぼすパラメータとなる。
【００４７】
図２に示した特性要因図に列記されたパラメータの中から上記（Ａ）〜（Ｄ）のパラメータを選択したのは、ロジスティック回帰方程式に基づくケッチン発生率の予測精度を表すパラメータである自由度二重調整寄与率が規定値（例えば０．７）以上であること及び自由度二重調整寄与率が最大となるように説明変数として用いる複数のパラメータを選択することを条件としたことによるものである。詳しくは、上記（Ａ）〜（Ｄ）のパラメータを説明変数としたケッチン発生率Ｐを予測するためのロジスティック回帰方程式は、下式（１）によって表される。
Ｉｎ｛Ｐ／（１−Ｐ）｝＝β０＋β１×Ｘ１＋β２×Ｘ２＋β３×Ｘ３＋β４×Ｘ４…（１）
上式（１）において、Ｘ１はエンジン回転速度ＮＥを示し、Ｘ２はスロットル開度ＶＴＡを示し、Ｘ３は油温ＴＨＯを示し、Ｘ４は圧縮比Ｒｃを示す。また、β０は、上記回帰方程式の切片を示し、β１〜β４はそれぞれ、エンジン回転速度ＮＥ、スロットル開度ＶＴＡ、油温ＴＨＯ、圧縮比Ｒｃについての偏回帰係数を示す。ここで上式（１）を変形すると、ケッチン発生率Ｐを下式（２）によって予測することが可能となる。
Ｐ＝１／［１＋ｅｘｐ｛―（β０＋β１×Ｘ１＋β２×Ｘ２＋β３×Ｘ３＋β４×Ｘ４）｝］…（２）
上式（２）における切片β０と、上記偏回帰係数β１〜β４とは、実機試験によって適合される。詳しくはまず、上記（Ａ）〜（Ｄ）のパラメータの水準を車両の実際の使用状態で想定される範囲内で割り振ることで、ケッチン発生率Ｐを測定するための試験条件を設定する。具体的には、図３に示すように、エンジン回転速度ＮＥとして、クランキング時に想定されるエンジン回転速度ＮＥの上下限値である１００ｒｐｍ，７００ｒｐｍの２水準を設定し、スロットル開度ＶＴＡとして、スロットル開度の最大値及び最小値の２水準を設定し、油温ＴＨＯとして、低温環境下において想定される油温（−３０℃）及びエンジン１０の暖機状態において想定される油温（８０℃）の２水準を設定する。また、圧縮比Ｒｃとして、エンジン１０のシリンダヘッドとシリンダブロックとの間に設けられるガスケットの厚さの変更によって燃焼室１８の容積を変化させることで、大小２水準を設定する。そして、測定回数の増大を抑制すべく、実験計画法に従い上記（Ａ）〜（Ｄ）のパラメータの各水準を組み合わせることで試験条件を決定する。詳しくは、直交表を用いてこれらパラメータの各水準を組み合わせることで試験条件を決定する。
【００４８】
次に、決定された試験条件に従った上記パラメータの各水準の組み合わせのそれぞれについて複数回（２０回）の実験を行い、実際にケッチンが発生した回数を上記複数回で除算した値をケッチン発生率Ｐとして算出する。そして、上記試験条件に従った上記パラメータの各水準の組み合わせに対するケッチン発生率Ｐの算出値に基づき、上記切片β０及び偏回帰係数β１〜β４を算出する。本実施形態では、切片β０は−７．８４５となり、エンジン回転速度ＮＥについての偏回帰係数β１は−０．００６となり、スロットル開度ＶＴＡについての偏回帰係数β２は０．０１７となり、油温ＴＨＯについての偏回帰係数β３は−０．０２８となり、圧縮比Ｒｃについての偏回帰係数β４は０．８１９となった。
【００４９】
ここで、これら偏回帰係数β１〜β４がケッチン発生率Ｐの予測精度に及ぼす影響について検証する。詳しくは、エンジン回転速度ＮＥについての偏回帰係数β１の符号が負となることから、上式（２）においてエンジン回転速度ＮＥが低いほどケッチン発生率Ｐの予測値が高くなり、スロットル開度ＶＴＡについての偏回帰係数β２の符号が正となることから、上式（２）においてスロットル開度ＶＴＡが高いほどケッチン発生率Ｐの予測値が高くなる。また、油温ＴＨＯについての偏回帰係数β３の符号が負となることから、上式（２）において油温ＴＨＯが低いほどケッチン発生率Ｐの予測値が高くなり、圧縮比Ｒｃについての偏回帰係数β４の符号が正となることから、上式（２）において圧縮比Ｒｃが高いほどケッチン発生率Ｐの予測値が高くなる。すなわち、エンジン回転速度ＮＥ、スロットル開度ＶＴＡ、油温ＴＨＯ及び圧縮比Ｒｃのそれぞれの変化によってケッチン発生率Ｐが変化する上述した傾向が反映されている。そして、算出された切片β０及び偏回帰係数β１〜β４を上式（２）に代入することにより決定されたロジスティック回帰方程式について、この回帰方程式に基づくケッチン発生率Ｐの予測精度を表す自由度二重調整寄与率Ｒａｄ＾２が０．７以上となって且つ最大となる結果が得られた（図４参照）。したがって、上式（２）で表されるロジスティック回帰方程式を用いることで、ケッチン発生率Ｐを適切に予測することが可能となる。
【００５０】
なお、ケッチン発生率Ｐの予測精度を表すパラメータとして自由度二重調整寄与率Ｒａｄ＾２を用いるのは、上記予測精度を向上させるべく説明変数として用いるパラメータを適切に選択するためである。つまり、自由度二重調整寄与率Ｒａｄ＾２に代えて寄与率Ｒ＾２を用いる場合、ケッチン発生率Ｐと相関が低いパラメータを説明変数として選択するときであっても、説明変数の数が増大するほど寄与率Ｒ＾２が増大することとなり（図４中、△で表記）、ケッチン発生率Ｐの予測精度を適切に検証することができなくなるおそれがある。これに対し、自由度二重調整寄与率Ｒａｄ＾２を用いる場合には、ケッチン発生率Ｐと相関の低いパラメータを説明変数として選択すると、自由度二重調整寄与率Ｒａｄ＾２が減少することとなり、ケッチン発生率Ｐの予測精度を向上させるためのパラメータを適切に選択することが可能となる。また、実機試験によって適合される上記切片β０及び偏回帰係数β１〜β４は、車両の仕様毎に異なる値になると考えられる。このため、上記切片β０及び偏回帰係数β１〜β４は、車両の仕様毎に適合されることが望ましい。
【００５１】
図５に、本実施形態にかかるケッチン発生回避処理の手順を示す。この処理は、ＥＣＵ３８によって、例えば所定周期で繰り返し実行される。なお本実施形態では、ＥＣＵ３８の不揮発性メモリに、上式（２）で表されるロジスティック回帰方程式の切片β０、偏回帰係数β１〜β４及び説明変数Ｘ４としてのエンジン１０の仕様から定まる所定の圧縮比（固定値）が予め記憶されている。また、ＥＣＵ３８の起動直後、説明変数Ｘ４に上記所定の圧縮比が格納される。
【００５２】
この一連の処理では、ステップＳ１０において、圧縮上死点となるタイミングを基準とした１燃焼サイクル（７２０°ＣＡ）の間に、上式（２）で表されるロジスティック回帰方程式の各説明変数Ｘ１〜Ｘ３として用いるエンジン回転速度ＮＥ、スロットル開度ＶＴＡ及び油温ＴＨＯの検出値を取得する処理を行う。詳しくは、説明変数Ｘ１の値として、クランク角度センサ３４の出力値から算出されるエンジン回転速度ＮＥの検出値を取得する。本実施形態では、エンジン回転速度ＮＥの検出値として、行程判別が完了していないと判断された場合には、上死点直前となるタイミング毎（３６０°ＣＡ毎）のエンジン回転速度ＮＥの検出値を取得し、行程判別が完了していると判断された場合には、圧縮上死点直前となるタイミング毎（７２０°ＣＡ）のエンジン回転速度ＮＥの検出値を取得する。これは、ケッチン発生率Ｐの予測精度を向上させるためである。つまり、ケッチンは、圧縮上死点直前に燃焼が開始されることに起因して発生するため、圧縮上死点直前におけるエンジン回転速度ＮＥは、ケッチン発生率を的確に予測するためのパラメータとなる。なお、行程判別が完了していない場合において、上死点直前となるタイミングは、クランク信号に基づき、欠け歯部が最初に検出されてからクランク軸２６が所定角度回転したタイミングとして把握すればよい。
【００５３】
一方、説明変数Ｘ２の値として、スロットル開度ＶＴＡの検出値を取得する。本実施形態では、スロットル開度ＶＴＡの検出値として、行程判別が完了していないと判断された場合には、吸気行程及び膨張行程であると想定される期間内の所定のタイミング毎（３６０°ＣＡ毎）のスロットル開度ＶＴＡの検出値を取得し、行程判別が完了していると判断された場合には、吸気バルブ２０の開弁期間内の所定のタイミング毎（７２０°ＣＡ毎）のスロットル開度ＶＴＡの検出値を取得する。これは、上記ステップＳ１０と同様に、ケッチン発生率Ｐの予測精度を向上させるためである。つまり、燃焼室１８に実際に吸気が供給される期間である吸気バルブ２０の開弁期間におけるスロットル開度ＶＴＡは、吸気量と高い相関を有する。なお、行程判別が完了していないと判断された場合において、吸気行程及び膨張行程であると想定される期間内の所定のタイミングは、クランク信号に基づき、欠け歯部が最初に検出されてからクランク軸２６が所定角度回転したタイミングとして把握すればよい。
【００５４】
他方、説明変数Ｘ３として、油温センサ３６の出力値から算出される油温ＴＨＯの検出値を取得する。本実施形態では、油温ＴＨＯのサンプリングが所定時間（例えば１００ｍｓｅｃ）毎に実施されるため、直近に検出された油温ＴＨＯを説明変数Ｘ３として用いる。なお、エンジン回転速度ＮＥやスロットル開度ＶＴＡと異なり、油温ＴＨＯは短時間に変化しない。このため、上記直近に検出された油温ＴＨＯを説明変数Ｘ３として用いる場合であっても、ケッチン発生率の予測精度に及ぼす影響は小さいと考えられる。
【００５５】
続くステップＳ１２では、複数の説明変数Ｘ１〜Ｘ４を入力として上記ロジスティック回帰方程式に基づきケッチン発生率Ｐを予測する。本実施形態では、行程判別が完了していないと判断された場合には、上死点となるタイミング毎（３６０°ＣＡ毎）にケッチン発生率Ｐを予測し、行程判別が完了したと判断された場合には、圧縮上死点直前となるタイミング毎（７２０°ＣＡ毎）にケッチン発生率Ｐを予測する。これにより、例えばケッチン発生率Ｐを都度予測することに伴うＥＣＵ３８の演算負荷の増大を抑制しつつ、ケッチン発生率Ｐと高い相関を有する圧縮上死点直前のエンジン回転速度ＮＥを用いてケッチン発生率Ｐを予測することが可能となる。
【００５６】
ステップＳ１２の処理の完了後、ステップＳ１４において、予測されたケッチン発生率Ｐが第１の閾値ＰＫ１（＞０）以上であるか否かを判断する。この処理は、現在の状況が、ケッチンが発生するおそれのある状況であるか否かを判断するための処理である。ここで上記第１の閾値ＰＫ１は、ケッチンが発生するおそれのないことを判別可能な値として予め実験等に基づき適合すればよい。
【００５７】
ステップＳ１４においてケッチン発生率Ｐが第１の閾値ＰＫ１未満であると判断された場合には、ケッチンが発生するおそれがないと判断し、ステップＳ１６において通常の点火制御処理を行う。
【００５８】
一方、上記ステップＳ１４においてケッチン発生率Ｐが第１の閾値ＰＫ１以上であると判断された場合には、ケッチンが発生するおそれがあると判断し、ステップＳ１８においてケッチン発生率Ｐが第１の閾値ＰＫ１よりも高い値として設定される第２の閾値ＰＫ２以上であるか否かを判断する。この処理は、現在の状況が、ケッチンが発生するおそれの度合いが大きい状況であるか否かを判断するための処理である。ここで上記第２の閾値ＰＫ２は、ケッチンが発生するおそれの度合いが大きい状況を判別可能な値として予め実験等に基づき適合すればよい。
【００５９】
ステップＳ１８においてケッチン発生率Ｐが第２の閾値ＰＫ２未満であると判断された場合には、ステップＳ２０に進み、点火タイミングを遅角する処理を行う。本実施形態では、上記遅角する処理として、行程判別が完了していないと判断された場合には、点火タイミングを上死点以降に設定する処理を行い、行程判別が完了したと判断された場合には、点火タイミングを圧縮上死点以降に設定する処理であるＡＴＤＣ点火処理を行う。この処理は、ケッチンの発生を適切に回避するとともに、エンジン１０の始動性の低下を抑制するための処理である。つまり、点火タイミングを圧縮上死点以降に設定すると、混合気の燃焼に伴い発生するエネルギをクランク軸２６を正回転させる方向のエネルギとして適切に用いることが可能となる。そしてＡＴＤＣ点火処理をケッチンが発生するおそれのある状況に限って行うことで、エンジン１０の始動性の低下を抑制することが可能となる。
【００６０】
なお、上記ＡＴＤＣ点火処理によって点火タイミングを過度に遅角させないことが望ましい。これは、点火タイミングが過度に遅角されることで混合気の燃焼を適切に開始することができなくなるおそれがあるためである。
【００６１】
一方、上記ステップＳ１８においてケッチン発生率Ｐが第２の閾値ＰＫ２以上であると判断された場合には、ステップＳ２２に進み、点火プラグ２２による点火を禁止する処理（点火禁止処理）を行う。この処理は、ケッチンの発生を回避するための処理である。つまり、点火プラグ２２の点火を禁止すると、混合気の燃焼が行われず、混合気の燃焼に伴いクランク軸２６に作用するクランク軸２６を逆回転させる方向の回転力を大きく低下させることが可能となる。
【００６２】
なお、ステップＳ１６、Ｓ２０、Ｓ２２の処理が完了する場合には、この一連の処理を一旦終了する。
【００６３】
図６に、本実施形態にかかるケッチン発生回避処理の一例を示す。
【００６４】
まず、行程判別が完了していない場合におけるケッチン発生回避処理について説明する。詳しくは、図６（ａ−１）に吸気圧の推移を示し、図６（ｂ−１）にクランク信号の推移を示し、図６（ｃ−１）にイグニッションスイッチの操作状態の推移を示し、図６（ｄ−１）に圧縮比Ｒｃの取得タイミングの推移を示し、図６（ｅ−１）に油温ＴＨＯの検出タイミングの推移を示し、図６（ｆ−１）にスロットル開度ＶＴＡの検出タイミングの推移を示し、図６（ｇ−１）にエンジン回転速度ＮＥの検出タイミングの推移を示し、図６（ｈ−１）にケッチン発生率Ｐの算出タイミングの推移を示す。なお、クランク信号付近に付された番号は、欠け歯部が最初に検出された後、最初に検出されるクランク信号の立ち下がりエッジを基準（０番）として、クランク信号の立ち下がりエッジ毎にカウントした値である。
【００６５】
図示されるように、行程判別が完了していない場合、吸気行程であると想定される期間内の所定のタイミング（時刻ｔ２）又は膨張行程であると想定される期間内の所定のタイミング（時刻ｔ５）にスロットル開度ＶＴＡが検出され、上死点直前となるタイミング（時刻ｔ１，ｔ４）においてエンジン回転速度ＮＥが検出される。すなわち、スロットル開度ＶＴＡ及びエンジン回転速度ＮＥは、３６０°ＣＡ毎に検出される。ここで上死点直前となるタイミング（時刻ｔ１，ｔ４）においてケッチン発生率Ｐが予測され、予測されたケッチン発生率Ｐに基づき、通常の点火制御処理、ＡＴＤＣ点火処理及び点火禁止処理のうちいずれを実行するかが決定される。
【００６６】
次に、行程判別が完了した場合におけるケッチン発生回避処理について説明する。なお、図６（ａ−２）〜図６（ｈ−２）は、図６（ａ−１）〜図６（ｈ−１）に対応している。
【００６７】
図示されるように、行程判別が完了した場合には、吸気バルブ２０の開弁期間内の所定のタイミング（時刻ｔ１）においてスロットル開度ＶＴＡが検出され、圧縮上死点直前となるタイミング（時刻ｔ２）においてエンジン回転速度ＮＥが検出される。すなわち、スロットル開度ＶＴＡ及びエンジン回転速度ＮＥは、７２０°ＣＡ毎に検出される。ここで圧縮上死点直前となるタイミング（時刻ｔ２）においてケッチン発生率Ｐが予測され、予測されたケッチン発生率Ｐに基づき、通常の点火制御処理、ＡＴＤＣ点火処理及び点火禁止処理のうちいずれを実行するかが決定される。
【００６８】
このように、本実施形態では、ロジスティック回帰方程式に基づきケッチン発生率Ｐを予測し、予測されたケッチン発生率Ｐに基づき点火制御処理を変更するケッチン発生回避処理を行うことで、エンジン１０の始動時におけるケッチンの発生を好適に回避するとともに、エンジン１０の始動性の低下を好適に抑制することができる。
【００６９】
以上詳述した本実施形態によれば、以下の効果が得られるようになる。
【００７０】
（１）ロジスティック回帰方程式に基づくケッチン発生率Ｐの予測精度を表すパラメータである自由度二重調整寄与率が規定値以上であること及び自由度二重調整寄与率が最大となるように説明変数として用いる複数のパラメータを選択することを条件として、実機試験によって上式（２）で表されるロジスティック回帰方程式を設定した。そして、上記ロジスティック回帰方程式を用いてケッチン発生率Ｐを予測した。これにより、ケッチンの発生を検出するためのセンサ等を備えることなくケッチンが発生するおそれのある状況を適切に把握することができ、ひいては部品数や車両重量の増大を好適に抑制したり、ケッチンが発生するおそれのある状況を制御ロジックの変更のみで把握したりすることができる。更に、ケッチン発生率Ｐについてロジスティック回帰方程式を用いてモデル化するため、ケッチン発生率Ｐの予測に用いる回帰方程式を決定するための作業工数を低減することなどもできる。
【００７１】
（２）ケッチン発生率Ｐを予測するために用いるエンジン回転速度ＮＥとして、圧縮上死点直前（又は上死点直前）のタイミングで検出されるエンジン回転速度ＮＥを用いた。これにより、ケッチン発生率Ｐの予測精度を好適に向上させることができる。
【００７２】
（３）予測されたケッチン発生率Ｐが第１の閾値ＰＫ１以上であって且つ第２の閾値ＰＫ２未満であると判断された場合、ＡＴＤＣ点火処理を行った。これにより、ケッチンの発生を好適に回避するとともに、エンジン１０の始動性の低下を好適に抑制することができる。
【００７３】
（４）予測されたケッチン発生率Ｐが第２の閾値ＰＫ２以上であると判断された場合、点火禁止処理を行った。これにより、ケッチンの発生を好適に回避することができる。
【００７４】
（その他の実施形態）
なお、上記実施形態は、以下のように変更して実施してもよい。
【００７５】
・上記実施形態において、上式（１）の右辺にエンジン回転速度ＮＥ、スロットル開度ＶＴＡ、油温ＴＨＯ及び圧縮比Ｒｃについての交互作用項を加えてもよい。具体的には例えば、エンジン回転速度ＮＥの低下度合いに対するケッチン発生率Ｐの上昇度合いが、スロットル開度ＶＴＡが大きいほど高くなるならば、エンジン回転速度ＮＥ及びスロットル開度ＶＴＡについての交互作用項「β５×Ｘ１×Ｘ２（β５：偏回帰係数）」を加えればよい。また例えば、エンジン回転速度の低下度合いに対するケッチン発生率Ｐの上昇度合いが、油温ＴＨＯが低いほど高くなるならば、エンジン回転速度ＮＥ及び油温ＴＨＯについての交互作用項「β６×Ｘ１×Ｘ３（β６：偏回帰係数）」を加えればよい。
【００７６】
・上記実施形態において、圧縮比についての偏回帰係数と説明変数との乗算値「β４×Ｘ４」を切片β０に含めてもよい。
【００７７】
・ロジスティック回帰方程式で用いる説明変数としては、上記実施形態に例示したものに限らない。例えば、吸気量と相関を有するパラメータであるスロットル開度に代えて、吸気圧を用いてもよい。また例えば、燃焼室１８に供給される吸気量を直接検出する吸気量検出手段（例えば、吸気通路１２上に設けられるエアフローメータ）を備え、スロットル開度ＶＴＡに代えて吸気量検出手段によって検出された吸気量を用いてもよい。具体的には、吸気バルブ２０の開弁期間における吸気量の検出値の積分演算値や、上記開弁期間の所定タイミング（吸気量の検出値が最大となるタイミング）における吸気量の検出値を用いればよい。
【００７８】
・上記実施形態において、燃焼室１８に供給される吸気温を検出する吸気温センサを備え、上記ロジスティック回帰方程式の説明変数として吸気温を追加してもよい。これは、吸気温が高くなると、空気密度が低下することで吸気量が減少し、ケッチン発生率Ｐが低下することに基づくものである。この場合、吸気温についての偏回帰係数の符号は負になると考えられる。
【００７９】
・上記実施形態では、予測されたケッチン発生率Ｐに応じて、ＡＴＤＣ点火処理又は点火禁止処理のいずれかを使い分ける制御ロジックとしたがこれに限らない。例えば、予測されたケッチン発生率Ｐが所定の閾値以上になると判断された場合、ＡＴＤＣ点火処理又は点火禁止処理のうちいずれか一方のみを行う制御ロジックとしてもよい。
【００８０】
・上記実施形態では、所定時間毎に油温ＴＨＯを検出したがこれに限らず、例えば、エンジン回転速度ＮＥの検出タイミングと同様に、圧縮上死点直前のタイミング毎に検出するようにしてもよい。
【００８１】
・上記実施形態では、自由度二重調整寄与率Ｒａｄ＾２が最大となるようにケッチン発生率に影響を及ぼす複数のパラメータを選択したがこれに限らない。例えば、自由度二重調整寄与率Ｒａｄ＾２が規定値（０．７）以上となるならば、自由度二重調整寄与率Ｒａｄ＾２が最大になるとの条件を除いて上記パラメータを選択してもよい。
【００８２】
・上記実施形態では、ロジスティック回帰方程式の切片及び偏回帰係数を、予め実験により適合された値（固定値）としたがこれに限らない。例えば、実際にエンジン１０が使用される状況下において、エンジン回転速度ＮＥ等の各説明変数の値に対応するケッチンの発生の有無に関する情報を記憶させる手段（ＥＣＵ３８の不揮発性メモリ）と、記憶された上記各説明変数の値に対応するケッチンの発生の有無に関する情報に基づき、ケッチン発生率Ｐを算出する手段とを備え、算出されたケッチン発生率Ｐに基づき上記偏回帰係数等を更新する処理を行ってもよい。これにより、エンジン１０の経年劣化等の影響を加味することでケッチン発生率Ｐの予測精度を維持することなどが期待できる。
【００８３】
なおこの場合、圧縮比についての説明変数を含めない構成とすることが望ましい。こうした構成によれば、車両の仕様毎に圧縮比の影響を適切に反映した切片β０が自動的に算出されることとなる。
【００８４】
・重回帰方程式に基づくケッチン発生率Ｐの予測手法としては、上記実施形態に例示したものに限らない。例えば、複数の説明変数のＮ（Ｎ：正の整数）次多項式として表される重回帰方程式によってケッチン発生率Ｐを予測することが可能ならば、上記重回帰方程式に基づきケッチン発生率Ｐを予測してもよい。ただしこの場合、ケッチン発生率Ｐを算出するための演算式が過度に複雑となることでＥＣＵ３８の演算負荷が過度に増大しないように重回帰方程式を決定するのが望ましい。また例えば、回帰方程式を用いたケッチン発生率Ｐの予測に代えて、エンジン回転速度ＮＥ、スロットル開度ＶＴＡ及び油温ＴＨＯを入力として、エンジン回転速度ＮＥ、スロットル開度ＶＴＡ及び油温ＴＨＯと関係付けられたケッチン発生率Ｐが規定されたマップを用いてケッチン発生率Ｐを予測してもよい。この場合、予め実機試験等によって得られるケッチン発生率Ｐ等の情報を重回帰分析して得られる重回帰方程式に基づき上記マップを作成すればよい。
【００８５】
・重回帰方程式に基づき予測される内燃機関の動作に関する異常発生率としては、ケッチン発生率に限らない。例えば、ノッキングが発生する確率（ノッキング発生率）としてもよい。この場合、説明変数として、燃料性状（例えば、アルコール濃度や燃料のオクタン価）や、吸気量、吸気温、冷却水温、更には圧縮比を用いればよい。詳しくは、アルコールと石油との混合割合が任意の燃料が使用されるエンジンにおいては、アルコール濃度が低いほど、燃料の燃焼に伴う熱発生量が多くなることに起因してノッキング発生率が高くなる。また、圧縮比が高かったり、吸気量が多かったりするほど、燃焼室１８での圧縮による混合気の温度の上昇度合いが大きくなることに起因してノッキング発生率が高くなる。更に、吸気温や水温が高いほど、燃焼室１８内の温度が高くなることに起因してノッキング発生率が高くなる。このため、上記各パラメータは、ノッキング発生率に影響を及ぼすパラメータとなる。
【００８６】
ここで、予測されたノッキング発生率に基づき、内燃機関の燃焼状態を制御するためのアクチュエータの操作状態を変更する処理としては、例えば点火プラグ２２の点火タイミングを遅角補正する処理とすればよい。具体的には、予測されたノッキング発生率が所定の閾値以上になると判断された場合、点火タイミングを所定角度遅角させる処理を行えばよい。これにより、例えば燃焼室１８内の圧力（筒内圧）を検出する筒内圧センサやノックセンサ等、ノッキングを検出する手段を備えることなく、ノッキングが発生するおそれのある状況を把握することができ、ノッキングを回避することができる。
【００８７】
また、上記異常発生率を、失火が発生する確率（失火発生率）としてもよい。この場合、説明変数として、圧縮比、点火タイミングについての情報、混合気の空燃比などを用いればよい。具体的には、圧縮比が高いほど、点火タイミングにおける筒内圧が高くなり、失火発生率が高くなる。これは、筒内圧が高いほど、点火プラグ２２の電極間に放電火花を発生させるために要求される上記電極間に印加される電圧（要求電圧）が高くなるため、筒内圧が高いほど、要求電圧が印加電圧を超えるおそれが大きくなり、失火発生率が高くなる。また、上記点火タイミングの情報として、機関運転状態毎に予め適合された燃焼状態を良好なものとする点火タイミング（基準タイミング）と実際の点火タイミングとの偏差を用いればよい。これは、上記偏差が大きくなるほど失火発生率が高くなることに鑑みたものである。また、混合気の空燃比が過度にリッチ側又はリーン側になる等、混合気の空燃比と基準となる空燃比（例えば理論空燃比）とのずれが大きくなるほど、失火発生率が高くなる。したがって、上記各パラメータは、失火発生率に影響を及ぼすパラメータになると考えられる。
【００８８】
ここで、予測された失火発生率に基づき、内燃機関の燃焼状態を制御するためのアクチュエータの操作状態を変更する処理としては、例えば点火プラグ２２の点火タイミングや燃料噴射弁１９からの燃料噴射量を補正する処理とすればよい。これにより、例えば失火を検出するセンサ（例えば筒内圧センサ）を備えることなく、失火が発生するおそれのある状況を把握することができ、失火を回避することができる。
【００８９】
なお、説明変数にエンジン回転速度を加えることで失火発生率を推定してもよい。具体的には、エンジン回転速度についての情報として、膨張行程の開始タイミングのエンジン回転速度から膨張行程の終了タイミングのエンジン回転速度を減算した値（回転速度差）を用いればよい。つまり、失火が発生すると、混合気の燃焼によってクランク軸２６に回転力が付与されないことから、膨張行程の開始タイミングのエンジン回転速度と比較して膨張行程の終了タイミングのエンジン回転速度が低くなると考えられる。このため、上記回転速度差が負の値となる場合、失火が発生している蓋然性が高いと考えられ、上記回転速度差が正の値となる場合には、失火が発生している蓋然性が低いと考えられる。したがって、上記回転速度差は、失火発生率に影響を及ぼすパラメータとなる。このため、上記回転速度差に加えて、点火タイミングについての情報や、混合気の空燃比を用いることで、失火発生率の推定精度を向上させることが期待できる。
【００９０】
また、エンジン１０の動作状態の推定対象としては、失火発生率に限らない。例えば、１燃焼サイクル（７２０°ＣＡ）のうち、ロータ２６ａの欠け歯部によって分割される期間の一方を表行程（例えば、排気行程及び吸気行程）とし、他方を裏行程（例えば、圧縮行程及び膨張行程）として定義し、現在の行程が表行程である確率（表行程率）を推定してもよい。この場合、説明変数として、吸気圧、バッテリ電圧及びエンジン回転速度についての情報等を用いればよい。具体的には例えば、吸気圧についての情報として、クランキングが開始されてから最初の欠け歯部が検出されるタイミングの吸気圧を用いればよい。つまり、欠け歯部が検出されるタイミングは、吸気行程から圧縮行程へと移行する期間内又は膨張行程から排気行程へと移行する期間内となり、吸気行程から圧縮行程へと移行する期間内の吸気圧は通常、膨張行程から排気行程へと移行する期間の吸気圧よりも低くなる。このため、上記最初の欠け歯部が検出されるタイミングの吸気圧が低いほど、上記最初の欠け歯部が検出されるタイミングが吸気行程から圧縮行程へと移行する期間内である確率が高くなる。したがって、上記最初の欠け歯部が検出されるタイミングの吸気圧は、表行程率に影響を及ぼすパラメータとなる。また、バッテリ電圧についての情報として、上記最初の欠け歯部が検出されるタイミングから所定期間内のバッテリ電圧の変動量を用いればよい。つまり、クランキング時の圧縮行程において燃焼室１８での圧縮に伴いクランク軸に作用する回転力が増大することに起因して、圧縮行程におけるバッテリ電圧の変動量は通常、他の行程におけるバッテリ電圧の変動量よりも大きくなる。このため、バッテリ電圧の変動量が大きいほど、上記最初の欠け歯部が検出されるタイミングが圧縮行程である確率が高くなる。したがって、バッテリ電圧の変動量は、表行程率に影響を及ぼすパラメータとなる。更に、エンジン回転速度についての情報として、上記最初の欠け歯部が検出されるタイミングから所定期間内のエンジン回転速度の変動量を用いればよい。つまり、クランキング時の圧縮行程において燃焼室１８での圧縮に伴いクランク軸に作用する回転力が増大することに起因して、圧縮行程におけるエンジン回転速度の変動量は通常、他の行程におけるエンジン回転速度の変動量よりも小さくなる。このため、エンジン回転速度の変動量が小さいほど、上記最初の欠け歯部が検出されるタイミングが圧縮行程である確率が高くなる。したがって、エンジン回転速度の変動量は、表行程率に影響を及ぼすパラメータとなる。
【００９１】
ここで、推定された表行程率に基づき行程判別を行えばよい。具体的には、クランキングが開始されてから欠け歯部が最初に検出されたことと、予測された表行程率とに基づき行程判別を行えばよい。予測された表行程率に基づく行程判別によれば、例えば上記実施形態で示した行程判別よりも早期に行程判別を完了することができ、その後のエンジン１０の始動制御や燃焼制御等を早期に実行することなどが期待できる。
【００９２】
・エンジン１０を始動させる手法としては、上記実施形態に例示したもの（スタータ式の始動手法）に限らない。例えば、クランク軸２６と機械的に連結されて且つクランク軸２６に初期回転を付与すべくユーザによって操作されるキックペダルを備え、キックペダルの操作によってエンジン１０を始動させるキック式のものであってもよい。この場合、ケッチンが発生することで、キックペダルが逆転し、キックペダルがユーザに接触する等、不都合が生じるおそれがある。このため、上記ケッチン発生回避処理を行うことで、上記不都合が生じる事態の発生を回避することができる。
【００９３】
・本願発明が適用される内燃機関としては、空冷式のものに限らず、水冷式のものであってもよい。この場合、エンジンを冷却する冷却水の水温を検出する水温センサを備え、説明変数として油温ＴＨＯに代えて、水温センサの出力値に基づく冷却水温を用いてもよい。また、内燃機関としては、単気筒のものに限らず複数の気筒を有するものであってもよく、火花点火式のものに限らず圧縮着火式のものであってもよい。
【００９４】
・本願発明が適用される車両としては、四輪バギーに限らない。例えば、自動二輪車（バイク）に適用してもよい。
【符号の説明】
【００９５】
１０…エンジン、１４ａ…スロットルセンサ、１６…吸気圧センサ、１８…燃焼室、２２…点火プラグ、３２…スタータ、３４…クランク角度センサ、３６…油温センサ、３８…ＥＣＵ（内燃機関の制御装置の一実施形態）。
【図１】
2012002204.tif 000003
【図２】
2012002204.tif 000004
【図３】
2012002204.tif 000005
【図４】
2012002204.tif 000006
【図５】
2012002204.tif 000007
【図６】
2012002204.tif 000008