JP6012892B1 - 内燃機関の制御装置及びその制御方法 - Google Patents

Abstract

【課題】クランク軸角度の検出誤差を精度良く補正することができる内燃機関の制御装置及びその制御方法を提供する。【解決手段】特定クランク角センサ６により角度間隔Δθｄ及び時間間隔ΔＴｄを検出する角度情報検出部５１と、角度間隔Δθｄ又は時間間隔ΔＴｄを補正値Ｋｃにより補正する角度情報補正部５２と、補正後の角度間隔Δθｄｃ及び時間間隔ΔＴｄｃに基づいて、クランク角速度ωｄ、クランク角加速度αｄ、及び角加速度変化量Δαｄを算出する角度情報算出部５３と、角加速度変化量Δαｄがゼロに近づくように補正値Ｋｃを変化させる補正値変化部５４と、を備えた内燃機関の制御装置５０。【選択図】図３

Description

本発明は、クランク軸と一体回転する回転部材に、予め定められた複数のクランク軸角度に設けられた複数の被検出部と、非回転部材に固定され、前記被検出部を検出する特定クランク角センサと、を備えた内燃機関の制御装置及びその制御方法に関するものである。

上記のような制御装置に関して、例えば下記の特許文献１及び特許文献２に記載された技術が既に知られている。特許文献１の技術では、クランク角センサの出力信号に基づいて、クランク角速度及びクランク角加速度を算出し、クランク角速度及びクランク角加速度に基づいて、燃焼により生じたガス圧トルクを算出し、ガス圧トルクによる仕事量を算出するように構成されている。

特許文献２の技術では、クランク角センサの出力信号に基づいてクランク軸の瞬時回転速度を検出し、瞬時回転速度に基づいて計測パラメータを算出し、瞬時回転速度の平均速度に対応する理想パラメータを記憶手段から取得し、理想パラメータに対する計測パラメータの誤差を学習するように構成されている。

特開２００９−２７５６１８号公報 特開２０１３−８７７２４号公報

ところで、被検出部に製造誤差があると、クランク角速度及びクランク角加速度に検出誤差が生じる。しかしながら、特許文献１には、どのようにして、この検出誤差に対応するか、具体的な方法は開示されていない。そのため、特許文献１の技術では、クランク角速度及びクランク角加速度に基づいて行う制御の精度が悪化するおそれがあった。

特許文献２の技術において、理想パラメータは、運転状態や外部負荷に応じて変化すると考えられるため、理想パラメータを基準とした学習では、誤学習が生じるおそれがある。

そこで、クランク軸角度の検出誤差を精度良く補正することができる内燃機関の制御装置及びその制御方法が求められる。

本発明に係る内燃機関の制御装置は、クランク軸と一体回転する回転部材に、予め定められた複数のクランク軸角度に設けられた複数の被検出部と、非回転部材に固定され、前記被検出部を検出する特定クランク角センサと、を備えた内燃機関の制御装置であって、前記特定クランク角センサの出力信号に基づいてクランク軸角度を検出すると共に前記クランク軸角度を検出した検出時刻を検出し、検出した前記クランク軸角度である検出角度及び前記検出時刻に基づいて、前記検出角度の間の角度区間に対応する角度間隔及び時間間隔を算出する角度情報検出部と、前記角度区間のそれぞれの前記角度間隔又は前記時間間隔を、前記角度区間のそれぞれに対応して１つずつ設けた補正値により補正する角度情報補正部と、前記角度区間のそれぞれの前記補正値による補正後の前記角度間隔及び前記
時間間隔に基づいて、前記検出角度又は前記角度区間のそれぞれに対応する、前記クランク軸角度の時間変化率であるクランク角速度、前記クランク角速度の時間変化率であるクランク角加速度、及び前記クランク角加速度の時間変化量である角加速度変化量を算出する角度情報算出部と、前記検出角度又は前記角度区間のそれぞれの前記角加速度変化量がゼロに近づくように、前記角度区間のそれぞれの前記補正値を変化させる補正値変化部と、を備えるものである。

また、本発明に係る内燃機関の制御方法は、クランク軸と一体回転する回転部材に、予め定められた複数のクランク軸角度に設けられた複数の被検出部と、非回転部材に固定され、前記被検出部を検出する特定クランク角センサと、を備えた内燃機関の制御方法であって、前記特定クランク角センサの出力信号に基づいてクランク軸角度を検出すると共に前記クランク軸角度を検出した検出時刻を検出し、検出した前記クランク軸角度である検出角度及び前記検出時刻に基づいて、前記検出角度の間の角度区間に対応する角度間隔及び時間間隔を算出する角度情報検出ステップと、前記角度区間のそれぞれの前記角度間隔又は前記時間間隔を、前記角度区間のそれぞれに対応して１つずつ設けた補正値により補正する角度情報補正ステップと、前記角度区間のそれぞれの前記補正値による補正後の前記角度間隔及び前記時間間隔に基づいて、前記検出角度又は前記角度区間のそれぞれに対応する、前記クランク軸角度の時間変化率であるクランク角速度、前記クランク角速度の時間変化率であるクランク角加速度、及び前記クランク角加速度の時間変化量である角加速度変化量を算出する角度情報算出ステップと、前記検出角度又は前記角度区間のそれぞれの前記角加速度変化量がゼロに近づくように、前記角度区間のそれぞれの前記補正値を変化させる補正値変化ステップと、を実行するものである。

本発明に係る内燃機関の制御装置及びその制御方法によれば、クランク角加速度の時間変化量である角加速度変化量がゼロに近づくように、補正値を変化させることにより、クランク軸角度の検出誤差を精度良く補正することができる。

本発明の実施の形態１に係る内燃機関及び制御装置の概略構成図である。 本発明の実施の形態１に係る内燃機関及び制御装置の概略構成図である。 本発明の実施の形態１に係る制御装置のブロック図である。 本発明の実施の形態１に係る制御装置のハードウェア構成図である。 本発明の実施の形態１に係る角度情報検出処理を説明するためのタイムチャートである。 本発明の実施の形態１に係る、記憶装置に記憶される補正値を説明するための図である。 本発明の実施の形態１に係る角度情報算出処理を説明するためのタイムチャートである。 本発明の実施の形態１に係る角加速度変化量がゼロより大きい場合の補正値の変化を説明するためのタイムチャートである。 本発明の実施の形態１に係る角加速度変化量がゼロより小さい場合の補正値の変化を説明するためのタイムチャートである。 本発明の実施の形態１に係る異なる運転状態における時間間隔の挙動を示すチャートである。 本発明の実施の形態１に係る未燃焼及び低負荷の運転状態における時間間隔の挙動を示すチャートである。 本発明の実施の形態１に係る、記憶装置に記憶される角度間隔及び時間間隔を説明するための図である。 本発明の実施の形態１に係る補正値変化処理を示すフローチャートである。 本発明の実施の形態１に係る補正なし及び補正ありの時間間隔の挙動を示すチャートである。 本発明の実施の形態１に係る補正なし及び補正ありの時間間隔のスペクトル解析結果を示す図である。 補正なしの比較例の場合の結果と、実施の形態１に係る補正ありの場合の結果を比較する図である。 本発明の実施の形態１に係る制御装置の概略的な処理を示すフローチャートである。 本発明のその他の実施の形態に係る角度情報算出処理を説明するためのタイムチャートである。

１．実施の形態１
実施の形態１に係る内燃機関１の制御装置５０（以下、単に制御装置５０と称す）について図面を参照して説明する。図１及び図２は、本実施の形態に係る内燃機関１及び制御装置５０の概略構成図であり、図３は、本実施の形態に係る制御装置５０のブロック図である。内燃機関１及び制御装置５０は、車両に搭載され、内燃機関１は、車両（車輪）の駆動力源となる。

１−１．内燃機関１の構成
まず、内燃機関１の構成について説明する。図１に示すように、内燃機関１は、空気と燃料の混合気を燃焼する気筒７を備えている。内燃機関１は、気筒７に空気を供給する吸気路２３と、気筒７で燃焼した排気ガスを排出する排気路１７とを備えている。内燃機関１は、ガソリンエンジンとされている。内燃機関１は、吸気路２３を開閉するスロットルバルブ４を備えている。スロットルバルブ４は、制御装置５０により制御される電気モータにより開閉駆動される電子制御式スロットルバルブとされている。スロットルバルブ４には、スロットルバルブ４の開度に応じた電気信号を出力するスロットル開度センサ１９が設けられている。

スロットルバルブ４の上流側の吸気路２３には、吸気路２３に吸入される吸入空気量に応じた電気信号を出力するエアフローセンサ３が設けられている。内燃機関１は、排気ガス還流装置２０を備えている。排気ガス還流装置２０は、排気路１７から吸気マニホールド１２に排気ガスを還流するＥＧＲ流路２１と、ＥＧＲ流路２１を開閉するＥＧＲバルブ２２と、を有している。吸気マニホールド１２は、スロットルバルブ４の下流側の吸気路２３の部分である。ＥＧＲバルブ２２は、制御装置５０により制御される電気モータにより開閉駆動される電子制御式ＥＧＲバルブとされている。

吸気マニホールド１２には、吸気マニホールド１２内の圧力に応じた電気信号を出力するマニホールド圧センサ８が設けられている。吸気マニホールド１２の下流側の部分には、燃料を噴射するインジェクタ１３が設けられている。なお、インジェクタ１３は、気筒７内に直接燃料を噴射するように設けられてもよい。内燃機関１には、大気圧に応じた電気信号を出力する大気圧センサ３３が設けられている。

気筒７の頂部には、空気と燃料の混合気に点火する点火プラグと、点火プラグに点火エネルギーを供給する点火コイル１６と、が設けられている。また、気筒７の頂部には、吸気路２３から気筒７内に吸入される吸入空気量を調節する吸気バルブ１４と、シリンダ内から排気路１７に排出される排気ガス量を調節する排気バルブ１５と、が設けられている。

図２に示すように、内燃機関１は、複数の気筒７（本例では３つ）を備えている。各気筒７内には、ピストン５が備えられている。各気筒７のピストン５は、コンロッド９及びクランク３２を介してクランク軸２に接続されている。クランク軸２は、ピストン５の往復運動によって回転駆動される。各気筒７で発生した燃焼ガス圧は、ピストン５の頂面を押圧し、コンロッド９及びクランク３２を介してクランク軸２を回転駆動する。クランク軸２は、車輪に駆動力を伝達する動力伝達機構に連結されている。動力伝達機構は、変速装置、ディファレンシャルギヤ等から構成される。

内燃機関１は、クランク軸２と一体回転する信号板１０を備えている。信号板１０は、予め定められた複数のクランク軸角度に複数の歯を設けている。本実施の形態では、信号板１０は、１０ｄｅｇ間隔で歯が並べられている。信号板１０の歯には、一部の歯が欠けた欠け歯部分が設けられている。内燃機関１は、エンジンブロック２４に固定され、信号板１０の歯を検出する第１クランク角センサ１１を備えている。

内燃機関１は、クランク軸２とチェーン２８で連結されたカム軸２９を備えている。カム軸２９は、吸気バルブ１４及び排気バルブ１５を開閉駆動する。クランク軸２が２回転する間に、カム軸２９は１回転する。内燃機関１は、カム軸２９と一体回転するカム用の信号板３１を備えている。カム用の信号板３１は、予め定められた複数のカム軸角度に複数の歯を設けている。内燃機関１は、エンジンブロック２４に固定され、カム用の信号板３１の歯を検出するカム角センサ３０を備えている。

制御装置５０は、第１クランク角センサ１１及びカム角センサ３０の２種類の出力信号に基づいて、各ピストン５の上死点を基準としたクランク軸角度を検出すると共に、各気筒７の行程を判別する。なお、内燃機関１は、吸入行程、圧縮行程、膨張行程、及び排気行程の４行程機関とされている。

内燃機関１は、クランク軸２と一体回転するフライホイール２７を備えている。フライホイール２７の外周部は、リングギア２５とされており、リングギア２５は、予め定められた複数のクランク軸角度に複数の歯を設けている。リングギア２５の歯は、周方向に等角度間隔で設けられている。本例では４ｄｅｇ間隔で、９０個の歯が設けられている。リングギア２５の歯には欠け歯部分は設けられていない。内燃機関１は、エンジンブロック２４に固定され、信号板１０の歯を検出する第２クランク角センサ６を備えている。第２クランク角センサ６は、リングギア２５の径方向外側に、リングギア２５と間隔を空けて対向配置されている。

第１クランク角センサ１１、カム角センサ３０、及び第２クランク角センサ６は、クランク軸２の回転による、各センサと歯の距離の変化に応じた電気信号を出力する。各角センサ１１、３０、６の出力信号は、センサと歯の距離が近い場合と、遠い場合とで信号がオンオフする矩形波となる。各角センサ１１、３０、６には、例えば、電磁ピックアップ式のセンサが用いられる。

フライホイール２７（リングギア２５）は、信号板１０の歯数よりも多い歯数を有しており、また、欠け歯部分もないため、高分解能の角度検出を期待できる。また、フライホイール２７は、信号板１０の質量よりも大きい質量を有しており、高周波振動が抑制されるため、高精度の角度検出を期待できる。

本実施の形態では、第２クランク角センサ６が、本発明における「特定クランク角センサ」に相当し、フライホイール２７が、本発明における「回転部材」に相当し、フライホイール２７に設けられたリングギア２５の歯が、本発明における「被検出部」に相当し、エンジンブロック２４が、本発明における「非回転部材」に相当する。

１−２．制御装置５０の構成
次に、制御装置５０について説明する。
制御装置５０は、内燃機関１を制御対象とする制御装置である。図３に示すように、制御装置５０は、角度情報検出部５１、角度情報補正部５２、角度情報算出部５３、及び補正値変化部５４、筒内圧力推定部５５、燃焼パラメータ算出部５６、及び燃焼制御部５７等の制御部を備えている。制御装置５０の各制御部５１〜５７等は、制御装置５０が備えた処理回路により実現される。具体的には、制御装置５０は、図４に示すように、処理回路として、ＣＰＵ（Central Processing Unit）等の演算処理装置９０（コンピュータ）
、演算処理装置９０とデータのやり取りをする記憶装置９１、演算処理装置９０に外部の信号を入力する入力回路９２、及び演算処理装置９０から外部に信号を出力する出力回路９３等を備えている。

記憶装置９１として、演算処理装置９０からデータを読み出し及び書き込みが可能に構成されたＲＡＭ（Random Access Memory）や、演算処理装置９０からデータを読み出し可能に構成されたＲＯＭ（Read Only Memory）等が備えられている。入力回路９２は、各種のセンサやスイッチが接続され、これらセンサやスイッチの出力信号を演算処理装置９０に入力するＡ／Ｄ変換器等を備えている。出力回路９３は、電気負荷が接続され、これら電気負荷に演算処理装置９０から制御信号を出力する駆動回路等を備えている。

そして、制御装置５０が備える各制御部５１〜５７等の各機能は、演算処理装置９０が、ＲＯＭ等の記憶装置９１に記憶されたソフトウェア（プログラム）を実行し、記憶装置９１、入力回路９２、及び出力回路９３等の制御装置５０の他のハードウェアと協働することにより実現される。なお、各制御部５１〜５７等が用いるテーブル、判定値等の設定データは、ソフトウェア（プログラム）の一部として、ＲＯＭ等の記憶装置９１に記憶されている。また、各制御部５１〜５７等が算出した補正値Ｋｃ、角度間隔Δθｄ、時間間隔ΔＴｄ、各算出値、及び各検出値等のデータは、ＲＡＭ等の書き換え可能な記憶装置９１に記憶される。

本実施の形態では、入力回路９２には、第１クランク角センサ１１、カム角センサ３０、第２クランク角センサ６、エアフローセンサ３、スロットル開度センサ１９、マニホールド圧センサ８、大気圧センサ３３、空燃比センサ１８、及びアクセルポジションセンサ２６等が接続されている。出力回路９３には、スロットルバルブ４（電気モータ）、ＥＧＲバルブ２２（電気モータ）、インジェクタ１３、及び点火コイル１６等が接続されている。なお、制御装置５０には、図示していない各種のセンサ、スイッチ、及びアクチュエータ等が接続されている。制御装置５０は、各種センサの出力信号に基づいて、吸入空気量、吸気マニホールド１２内の圧力、大気圧、空燃比、及びアクセル開度等の内燃機関１の運転状態を検出する。

制御装置５０は、基本的な制御として、入力された各種センサの出力信号等に基づいて、燃料噴射量、点火時期等を算出し、インジェクタ１３及び点火コイル１６等を駆動制御する。制御装置５０は、アクセルポジションセンサ２６の出力信号等に基づいて、運転者が要求している内燃機関１の出力トルクを算出し、当該要求出力トルクを実現する吸入空気量となるように、スロットルバルブ４等を制御する。具体的には、制御装置５０は、目標スロットル開度を算出し、スロットル開度センサ１９の出力信号に基づき検出したスロットル開度が、目標スロットル開度に近づくように、スロットルバルブ４の電気モータを駆動制御する。また、制御装置５０は、入力された各種センサの出力信号等に基づいて、ＥＧＲバルブ２２の目標開度を算出し、ＥＧＲバルブ２２の電気モータを駆動制御する。

＜角度情報検出部５１＞
角度情報検出部５１は、図５に示すように、特定クランク角センサとされた第２クランク角センサ６の出力信号に基づいてクランク軸角度θｄを検出すると共にクランク軸角度θｄを検出した検出時刻Ｔｄを検出する。そして、角度情報検出部５１は、検出したクランク軸角度θｄである検出角度θｄ及び検出時刻Ｔｄに基づいて、検出角度θｄの間の角度区間Ｓｄに対応する角度間隔Δθｄ及び時間間隔ΔＴｄを算出する。

本実施の形態では、角度情報検出部５１は、第２クランク角センサ６の出力信号（矩形波）の立下りエッジ（又は立上りエッジ）を検出した時のクランク軸角度θｄを判定するように構成されている。角度情報検出部５１は、基点角度（例えば、第１気筒７のピストン５の上死点である０deg）に対応する立下りエッジである基点立下りエッジを判定し、基点立下りエッジを基点にカウントアップした立下りエッジの番号ｎ（以下、角度識別番号ｎと称す）に対応するクランク軸角度θｄを判定する。例えば、角度情報検出部５１は、基点立下りエッジを検出した時に、クランク軸角度θｄを基点角度（例えば、０deg）に設定すると共に角度識別番号ｎを０に設定する。そして、角度情報検出部５１は、立下りエッジを検出する毎に、クランク軸角度θｄを、予め設定された角度間隔Δθｄ（本例では４deg）ずつ増加させると共に角度識別番号ｎを１つずつ増加させる。或いは、角度情報検出部５１は、角度識別番号ｎとクランク軸角度θｄとの関係が予め設定された角度テーブルを用い、今回の角度識別番号ｎに対応するクランク軸角度θｄを読み出すように構成されてもよい。角度情報検出部５１は、クランク軸角度θｄ（検出角度θｄ）を角度識別番号ｎに対応付ける。角度識別番号ｎは、最大番号（本例では９０）の後、１に戻る。角度識別番号ｎ＝１の前回の角度識別番号ｎは９０になり、角度識別番号ｎ＝９０の次回の角度識別番号ｎは１になる。

本実施の形態では、角度情報検出部５１は、後述する、第１クランク角センサ１１及びカム角センサ３０に基づいて検出した参照クランク軸角度θｒを参照して、第２クランク角センサ６の基点立下りエッジを判定する。例えば、角度情報検出部５１は、第２クランク角センサ６の立下りエッジを検出した時の参照クランク軸角度θｒが、基点角度に最も近い立下りエッジを、基点立下りエッジと判定する。

また、角度情報検出部５１は、第１クランク角センサ１１及びカム角センサ３０に基づいて判別した各気筒７の行程を参照して、クランク軸角度θｄに対応する各気筒７の行程を判定する。

角度情報検出部５１は、第２クランク角センサ６の出力信号（矩形波）の立下りエッジを検出した時の検出時刻Ｔｄを検出し、検出時刻Ｔｄを角度識別番号ｎに対応付ける。具体的には、角度情報検出部５１は、演算処理装置９０が備えたタイマー機能を用いて、検出時刻Ｔｄを検出する。

角度情報検出部５１は、図５に示すように、立下りエッジを検出した時に、今回の角度識別番号（ｎ）に対応する検出角度θｄ（ｎ）と、前回の角度識別番号（ｎ−１）に対応する検出角度θｄ（ｎ−１）との間の角度区間を、今回の角度識別番号（ｎ）に対応する角度区間Ｓｄ（ｎ）に設定する。

また、角度情報検出部５１は、式（１）に示すように、立下りエッジを検出した時に、今回の角度識別番号（ｎ）に対応する検出角度θｄ（ｎ）と、前回の角度識別番号（ｎ−１）に対応する検出角度θｄ（ｎ−１）との偏差を算出して、今回の角度識別番号（ｎ）（今回の角度区間Ｓｄ（ｎ））に対応する角度間隔Δθｄ（ｎ）に設定する。
Δθｄ（ｎ）＝θｄ（ｎ）−θｄ（ｎ−１） ・・・（１）
本実施の形態では、リングギア２５の歯の角度間隔は、全て等しくされているので、角度情報検出部５１は、全ての角度識別番号ｎの角度間隔Δθｄを、予め設定された角度（本例では４deg）に設定する。

また、角度情報検出部５１は、式（２）に示すように、立下りエッジを検出した時に、今回の角度識別番号（ｎ）に対応する検出時刻Ｔｄ（ｎ）と、前回の角度識別番号（ｎ−１）に対応する検出時刻Ｔｄ（ｎ−１）との偏差を算出して、今回の角度識別番号（ｎ）（今回の角度区間Ｓｄ（ｎ））に対応する時間間隔ΔＴｄ（ｎ）に設定する。
ΔＴｄ（ｎ）＝Ｔｄ（ｎ）−Ｔｄ（ｎ−１） ・・・（２）

角度情報検出部５１は、第１クランク角センサ１１及びカム角センサ３０の２種類の出力信号に基づいて、第１気筒７のピストン５の上死点を基準とした参照クランク軸角度θｒを検出すると共に、各気筒７の行程を判別する。例えば、角度情報検出部５１は、第１クランク角センサ１１の出力信号（矩形波）の立下りエッジの時間間隔から、信号板１０の欠け歯部分の直後の立下りエッジを判定する。そして、角度情報検出部５１は、欠け歯部分の直後の立下りエッジを基準にした各立下りエッジと、上死点を基準にした参照クランク軸角度θｒと対応関係を判定し、各立下りエッジを検出した時の、上死点を基準とした参照クランク軸角度θｒを算出する。また、角度情報検出部５１は、第１クランク角センサ１１の出力信号（矩形波）における欠け歯部分の位置と、カム角センサ３０の出力信号（矩形波）との関係から、各気筒７の行程を判別する。

＜角度情報補正部５２＞
角度情報補正部５２は、角度区間Ｓｄのそれぞれの角度間隔Δθｄ又は時間間隔ΔＴｄを、角度区間Ｓｄのそれぞれに対応して１つずつ設けた補正値Ｋｃにより補正する。

本実施の形態では、角度情報補正部５２は、各角度識別番号ｎの角度区間Ｓｄ（ｎ）に１つずつ補正値Ｋｃ（ｎ）を設けている。本例では、角度識別番号ｎ及び角度区間Ｓｄは９０設けられているので、補正値Ｋｃも９０設けられている。各補正値Ｋｃは、図６に示すように、各角度識別番号ｎに対応付けられて、制御装置５０のＲＡＭ等の書き換え可能な記憶装置９１に記憶される。

角度情報補正部５２は、式（３）に示すように、今回の角度識別番号（ｎ）に対応する角度間隔Δθｄ（ｎ）又は時間間隔ΔＴｄ（ｎ）に、今回の角度識別番号（ｎ）に対応する補正値Ｋｃ（ｎ）を乗算して、今回の角度識別番号（ｎ）に対応する補正後の角度間隔Δθｄｃ（ｎ）又は時間間隔ΔＴｄｃ（ｎ）を算出するように構成されている。
Δθｄｃ（ｎ）＝Ｋｃ（ｎ）×Δθｄ（ｎ）
又は ・・・（３）
ΔＴｄｃ（ｎ）＝Ｋｃ（ｎ）×ΔＴｄ（ｎ）

本実施の形態では、補正値Ｋｃにより時間間隔ΔＴｄが補正される場合について説明する。なお、補正値Ｋｃにより補正されていない角度間隔Δθｄも、説明の便宜上、補正後の角度間隔Δθｄｃと称す。

＜角度情報算出部５３＞
角度情報算出部５３は、角度区間Ｓｄのそれぞれの補正値Ｋｃによる補正後の角度間隔Δθｄｃ及び時間間隔ΔＴｄｃに基づいて、検出角度θｄ又は角度区間Ｓｄのそれぞれに対応する、クランク軸角度θｄの時間変化率であるクランク角速度ωｄ、クランク角速度ωｄの時間変化率であるクランク角加速度αｄ、及びクランク角加速度αｄの時間変化量である角加速度変化量Δαｄを算出する。

本実施の形態では、図７に示すように、角度情報算出部５３は、処理対象とする角度区間Ｓｄ（ｎ）に対応する補正後の角度間隔Δθｄｃ（ｎ）及び時間間隔ΔＴｄｃ（ｎ）に基づいて、処理対象の角度区間Ｓｄ（ｎ）に対応するクランク角速度ωｄ（ｎ）を算出する。具体的には、角度情報算出部５３は、式（４）に示すように、処理対象の角度区間Ｓｄ（ｎ）に対応する補正後の角度間隔Δθｄｃ（ｎ）を補正後の時間間隔ΔＴｄｃ（ｎ）で除算して、クランク角速度ωｄ（ｎ）を算出する。
ωｄ（ｎ）＝Δθｄｃ（ｎ）／ΔＴｄｃ（ｎ） ・・・（４）

角度情報算出部５３は、処理対象とする検出角度θｄ（ｎ）の直前１つの角度区間Ｓｄ（ｎ）に対応するクランク角速度ωｄ（ｎ）及び補正後の時間間隔ΔＴｄｃ（ｎ）、並びに処理対象の検出角度θｄ（ｎ）の直後１つの角度区間Ｓｄ（ｎ＋１）に対応するクランク角速度ωｄ（ｎ＋１）及び補正後の時間間隔ΔＴｄｃ（ｎ＋１）に基づいて、処理対象の検出角度θｄ（ｎ）に対応するクランク角加速度αｄ（ｎ）を算出する。具体的には、角度情報算出部５３は、式（５）に示すように、直後のクランク角速度ωｄ（ｎ＋１）から直前のクランク角速度ωｄ（ｎ）を減算した減算値を、直後の補正後の時間間隔ΔＴｄｃ（ｎ＋１）と直前の補正後の時間間隔ΔＴｄｃ（ｎ）の平均値で除算して、クランク角加速度αｄ（ｎ）を算出する。
αｄ（ｎ）＝｛ωｄ（ｎ＋１）−ωｄ（ｎ）｝
／｛ΔＴｄｃ（ｎ＋１）＋ΔＴｄｃ（ｎ）｝×２ ・・・（５）

角度情報算出部５３は、処理対象とする検出角度θｄ（ｎ）の直前１つの検出角度θｄ（ｎ−１）に対応するクランク角加速度αｄ（ｎ−１）、及び処理対象の検出角度θｄ（ｎ）の直後１つの検出角度θｄ（ｎ＋１）に対応するクランク角加速度αｄ（ｎ＋１）に基づいて、処理対象の検出角度θｄ（ｎ）に対応する角加速度変化量Δαｄ（ｎ）を算出する。具体的には、角度情報算出部５３は、式（６）に示すように、直後のクランク角加速度αｄ（ｎ＋１）から直前のクランク角加速度αｄ（ｎ−１）を減算して、角加速度変化量Δαｄ（ｎ）を算出する。このように、中心差分により角加速度変化量Δαｄを算出するように構成されている。
Δαｄ（ｎ）＝αｄ（ｎ＋１）−αｄ（ｎ−１） ・・・（６）

以上の処理をまとめると、１つの角加速度変化量Δαｄ（ｎ）を算出するために以下のような処理が実行される。角度情報算出部５３は、角加速度変化量Δαｄ（ｎ）を算出する検出角度θｄである対象検出角度θｄ（ｎ）の直前２つの角度区間Ｓｄ（ｎ−１）、Ｓｄ（ｎ）及び直後２つの角度区間Ｓｄ（ｎ＋１）、Ｓｄ（ｎ＋２）のそれぞれの補正後の角度間隔Δθｄｃ及び時間間隔ΔＴｄｃに基づいて、当該直前２つ及び直後２つの角度区間Ｓｄのそれぞれのクランク角速度ωｄ（ｎ−１）、ωｄ（ｎ）、ωｄ（ｎ＋１）、ωｄ（ｎ＋２）を算出する。

そして、角度情報算出部５３は、対象検出角度θｄ（ｎ）の直前２つの角度区間Ｓｄのそれぞれのクランク角速度ωｄ（ｎ−１）、ωｄ（ｎ）及び時間間隔ΔＴｄ（ｎ−１）、ΔＴｄ（ｎ）に基づいて、当該直前２つの角度区間Ｓｄに対応する１つのクランク角加速度αｄ（ｎ−１）を算出する。また、角度情報算出部５３は、対象検出角度θｄ（ｎ）の直後２つの角度区間Ｓｄのそれぞれのクランク角速度ωｄ（ｎ＋１）、ωｄ（ｎ＋２）及び時間間隔ΔＴｄ（ｎ＋１）、ΔＴｄ（ｎ＋２）に基づいて、当該直後２つの角度区間Ｓｄに対応する１つのクランク角加速度αｄ（ｎ＋１）を算出する。そして、角度情報算出部５３は、対象検出角度θｄ（ｎ）の直前２つの角度区間Ｓｄに対応するクランク角加速度αｄ（ｎ−１）及び直後２つの角度区間Ｓｄに対応するクランク角加速度αｄ（ｎ＋１）に基づいて、対象検出角度θｄ（ｎ）に対応する１つの角加速度変化量Δαｄ（ｎ）を算出する。

本実施の形態では、角度情報算出部５３は、角加速度変化量Δαｄとして、クランク角加速度αｄの時間変化率であるクランク角躍度δｄを算出するように構成されている。具体的には、角度情報算出部５３は、式（７）に示すように、対象検出角度θｄ（ｎ）の直後のクランク角加速度αｄ（ｎ＋１）から直前のクランク角加速度αｄ（ｎ−１）を減算した減算値を、対象検出角度θｄ（ｎ）の直後の補正後の時間間隔ΔＴｄｃ（ｎ＋１）と直前の補正後の時間間隔ΔＴｄｃ（ｎ）の合計値で除算して、対象検出角度θｄ（ｎ）のクランク角躍度δｄ（ｎ）を算出する。
δｄ（ｎ）＝｛αｄ（ｎ＋１）−αｄ（ｎ−１）｝
／｛ΔＴｄｃ（ｎ＋１）＋ΔＴｄｃ（ｎ）｝ ・・・（７）

以上のように、角加速度変化量Δαｄ又はクランク角躍度δｄを算出するためには、対象検出角度θｄ（ｎ）の直後２つの角度区間Ｓｄ（ｎ＋１）、Ｓｄ（ｎ＋２）の角度間隔Δθｄ及び時間間隔ΔＴｄが必要になる。そのため、角度情報算出部５３は、最新の検出角度θｄに対して、少なくとも２つ遅れの検出角度θｄを対象検出角度に設定して、角加速度変化量Δαｄ又はクランク角躍度δｄの算出を行う。

＜補正値変化部５４＞
補正値変化部５４は、検出角度θｄ又は角度区間Ｓｄ（本例では検出角度θｄ）のそれぞれの角加速度変化量Δαｄ（又はクランク角躍度δｄ）がゼロに近づくように、角度区間Ｓｄのそれぞれの補正値Ｋｃを変化させる。

図８の例に示すように、対象検出角度θｄ（ｎ）の角加速度変化量Δαｄ（ｎ）がゼロより大きい場合は、対象検出角度θｄ（ｎ）の直後のクランク角加速度αｄ（ｎ＋１）が、直前のクランク角加速度αｄ（ｎ−１）よりも大きく、クランク角加速度αｄが増加中である。そのため、角加速度変化量Δαｄ（ｎ）をゼロに近づけるためには、直前のクランク角加速度αｄ（ｎ−１）を増加させるか、直後のクランク角加速度αｄ（ｎ＋１）を減少させればよい。

直前のクランク角加速度αｄ（ｎ−１）を増加させるためには、直前のクランク角速度ωｄ（ｎ）を増加させればよい。そのためには、直前の補正後の角度間隔Δθｄｃ（ｎ）が増加するように、又は直前の補正後の時間間隔ΔＴｄｃ（ｎ）が減少するように、直前の補正値Ｋｃ（ｎ）を変化させればよい。一方、直後のクランク角加速度αｄ（ｎ＋１）を減少させるためには、直後のクランク角速度ωｄ（ｎ＋１）を減少させればよい。そのためには、直後の補正後の角度間隔Δθｄｃ（ｎ＋１）が減少するように、又は直後の補正後の時間間隔ΔＴｄｃ（ｎ＋１）が増加するように、直後の補正値Ｋｃ（ｎ＋１）を変化させればよい。

一方、図９の例に示すように、対象検出角度θｄ（ｎ）の角加速度変化量Δαｄ（ｎ）がゼロより小さい場合は、対象検出角度θｄ（ｎ）の直後のクランク角加速度αｄ（ｎ＋１）が、直前のクランク角加速度αｄ（ｎ−１）よりも小さく、クランク角加速度αｄが減少中である。そのため、角加速度変化量Δαｄ（ｎ）をゼロに近づけるためには、直前のクランク角加速度αｄ（ｎ−１）を減少させるか、直後のクランク角加速度αｄ（ｎ＋１）を増加させればよい。

直前のクランク角加速度αｄ（ｎ−１）を減少させるためには、直前のクランク角速度ωｄ（ｎ）を減少させればよい。そのためには、直前の補正後の角度間隔Δθｄｃ（ｎ）が減少するように、又は直前の補正後の時間間隔ΔＴｄｃ（ｎ）が増加するように、直前の補正値Ｋｃ（ｎ）を変化させればよい。一方、直後のクランク角加速度αｄ（ｎ＋１）を増加させるためには、直後のクランク角速度ωｄ（ｎ＋１）を増加させればよい。そのためには、直後の補正後の角度間隔Δθｄｃ（ｎ＋１）が増加するように、又は直後の補正後の時間間隔ΔＴｄｃ（ｎ＋１）が減少するように、直後の補正値Ｋｃ（ｎ＋１）を変化させればよい。

従って、補正値変化部５４は、以下のように構成されている。まず、補正値Ｋｃにより角度間隔Δθｄが補正されるように構成されている場合について説明する。補正値変化部５４は、対象検出角度θｄ（ｎ）の角加速度変化量Δαｄ（ｎ）が正の値である場合は、対象検出角度θｄ（ｎ）の直前１つの角度区間Ｓｄ（ｎ）の補正後の角度間隔Δθｄ（ｎ）が増加するように、当該角度区間Ｓｄ（ｎ）の補正値Ｋｃ（ｎ）を変化させる（本例では増加させる）増加補正、及び対象検出角度θｄ（ｎ）の直後１つの角度区間Ｓｄ（ｎ＋１）の補正後の角度間隔Δθｄ（ｎ＋１）が減少するように、当該角度区間Ｓｄ（ｎ＋１）の補正値Ｋｃ（ｎ＋１）を変化させる（本例では減少させる）減少補正の一方又は双方を実行する。一方、補正値変化部５４は、対象検出角度θｄ（ｎ）の角加速度変化量Δαｄ（ｎ）が負の値である場合は、対象検出角度θｄ（ｎ）の直前１つの角度区間Ｓｄ（ｎ）の補正後の角度間隔Δθｄ（ｎ）が減少するように、当該角度区間Ｓｄ（ｎ）の補正値Ｋｃ（ｎ）を変化させる（本例では減少させる）減少補正、及び対象検出角度θｄ（ｎ）の直後１つの角度区間Ｓｄ（ｎ＋１）の補正後の角度間隔Δθｄ（ｎ＋１）が増加するように、当該角度区間Ｓｄ（ｎ＋１）の補正値Ｋｃ（ｎ＋１）を変化させる（本例では増加させる）増加補正の一方又は双方を実行する。

次に、補正値Ｋｃにより時間間隔ΔＴｄが補正されるように構成されている場合について説明する。補正値変化部５４は、対象検出角度θｄ（ｎ）の角加速度変化量Δαｄ（ｎ）が正の値である場合は、対象検出角度θｄ（ｎ）の直前１つの角度区間Ｓｄ（ｎ）の補正後の時間間隔ΔＴｄ（ｎ）が減少するように、当該角度区間Ｓｄ（ｎ）の補正値Ｋｃ（ｎ）を変化させる（本例では減少させる）減少補正、及び対象検出角度θｄ（ｎ）の直後１つの角度区間Ｓｄ（ｎ＋１）の補正後の時間間隔ΔＴｄ（ｎ＋１）が増加するように、当該角度区間Ｓｄ（ｎ＋１）の補正値Ｋｃ（ｎ＋１）を変化させる（本例では増加させる）増加補正の一方又は双方を実行する。一方、補正値変化部５４は、対象検出角度θｄ（ｎ）の角加速度変化量Δαｄ（ｎ）が負の値である場合は、対象検出角度θｄ（ｎ）の直前１つの角度区間Ｓｄ（ｎ）の補正後の時間間隔ΔＴｄ（ｎ）が増加するように、当該角度区間Ｓｄ（ｎ）の補正値Ｋｃ（ｎ）を変化させる（本例では増加させる）増加補正、及び対象検出角度θｄ（ｎ）の直後１つの角度区間Ｓｄ（ｎ＋１）の補正後の時間間隔ΔＴｄ（ｎ＋１）が減少するように、当該角度区間Ｓｄ（ｎ＋１）の補正値Ｋｃ（ｎ＋１）を変化させる（本例では減少させる）減少補正の一方又は双方を実行する。

補正値変化部５４は、各角度区間Ｓｄの補正値Ｋｃを、予め設定された変化量ΔＫｃだけ増加又は減少させるように構成されている。なお、補正値変化部５４は、各角度区間Ｓｄの角加速度変化量Δαｄの絶対値に応じて、変化量ΔＫｃを変化させるように構成されてもよい。この場合は、補正値変化部５４は、各角度区間Ｓｄの角加速度変化量Δαｄの絶対値が大きくなに従って、変化量ΔＫｃを大きくする。また、補正値変化部５４は、各角度区間Ｓｄの角加速度変化量Δαｄの絶対値が、予め設定された不感帯判定値以下の場合は、補正値Ｋｃを変化させない不感帯処理を行うように構成されてもよい。

＜補正値Ｋｃの変化許可条件＞
図１０に示すように、内燃機関１に燃料が供給されており、燃焼している場合（一点鎖線）は、燃焼による気筒７内の圧力変動が大きいため、時間間隔ΔＴｄの時間変化が大きくなり、クランク角速度の時間変化が大きくなる。燃料供給が停止され、燃焼停止しており、内燃機関１の負荷が中程度である場合（破線）は、燃焼時よりも気筒７内の圧力変動が減少するが、気筒７内に吸入される空気量が多くなり、気筒７内の空気の圧力変動幅が大きくなる。そのため、時間間隔及びクランク角速度の時間変化は、燃焼時よりも減少するが十分に減少していない。

燃料供給が停止され、内燃機関１の負荷が低い場合（実線）は、気筒７内に吸入される空気量が少なくなり、気筒７内の空気の圧力変動幅が大幅に減少する。そのため、時間間隔及びクランク角速度の時間変化が大幅に減少し、クランク軸２の回転は、等速円運動に近づく。しかし、気筒７内には空気が残存しており、気筒７内の空気の圧力変動が多少生じる。図１１に、図１０の縦軸を拡大して示すように、時間間隔及びクランク角速度の時間変化は、減少するもののゼロにはならない。なお、図１１には、補正なしの時間間隔ΔＴｄを実線で示し、補正値Ｋｃによる補正の目標となる時間間隔を破線で示している。よって、燃焼停止及び低負荷の状態でも、本実施の形態とは異なり、クランク角速度ωｄの時間変化率であるクランク角加速度αｄがゼロに近づくように、補正値Ｋｃを変化させるように構成すると、誤補正が生じる。

一方、時間間隔及びクランク角速度は、一定の傾きで増加又は減少する挙動に近づき、クランク角加速度αｄの時間変化がゼロに近づく。よって、本実施の形態のように、クランク角加速度αｄの時間変化量である角加速度変化量Δαｄがゼロに近づくように、補正値Ｋｃを変化させることにより、適切に補正を行うことができる。

本実施の形態では、補正値変化部５４は、内燃機関１への燃料供給が停止されており、内燃機関１の負荷が予め設定された判定負荷よりも低い変化許可条件が成立した場合に、補正値Ｋｃの変化を許可するように構成されている。

補正値変化部５４は、スロットル開度センサ１９の出力信号に基づいて検出したスロットルバルブ４の開度が予め設定された判定開度よりも低い場合に成立する条件、エアフローセンサ３の出力信号に基づいて検出した内燃機関１（気筒７）に吸入される吸入空気量が予め設定された判定空気量よりも低い場合に成立する条件、及びマニホールド圧センサ８の出力信号に基づいて検出した吸気マニホールド１２内の圧力が予め設定された判定圧力よりも低い場合に成立する条件、内燃機関１の充填効率が予め設定された判定効率よりも低い場合に成立する条件のいずれか１つ以上の条件を論理積により組み合わせた条件が成立した場合に、内燃機関１の負荷が判定負荷よりも低いと判定するように構成されている。

例えば、判定開度は、スロットルバルブ４の開度の変化幅の５％以下の開度に予め設定される。判定空気量は、吸入空気量の変化幅の５％以下の空気量に予め設定される。判定圧力は、吸気マニホールド１２内の圧力の変化幅の５％以下の圧力に予め設定される。判定効率は、充填効率の変化幅の５％以下の充填効率に予め設定される。

＜補正値Ｋｃの変化完了判定＞
補正値変化部５４は、クランク軸２が予め設定された判定角度（本例では、３６０deg）回転する毎に、判定角度間の角加速度変化量Δαｄ（又はクランク角躍度δｄ）についての、当該角加速度変化量Δαｄの平均値に対するばらつき度合を算出する。本実施の形態では、補正値変化部５４は、ばらつき度合として、判定角度間の角加速度変化量Δαｄ（又はクランク角躍度δｄ）の標準偏差σを算出するように構成されている。そして、補正値変化部５４は、前回算出した標準偏差σと今回算出した標準偏差σとの変化量の絶対値｜Δσ｜が、予め設定された判定変化量Ｘσ未満になった場合に、補正値Ｋｃの変化が完了したと判定し、補正値Ｋｃの変化を停止し、補正値Ｋｃを保持する。

各角度区間Ｓｄの補正値Ｋｃが適切に変化されていくと、各検出角度θｄの角加速度変化量Δαｄがゼロに近づいていく。そして、各検出角度θｄの角加速度変化量Δαｄがゼロに近づいていくと、角加速度変化量Δαｄの標準偏差σの変化量の絶対値｜Δσ｜が次第に減少していく。上記の構成によれば、標準偏差σの変化量の絶対値｜Δσ｜が、判定変化量Ｘσ以下まで減少した場合に、補正値Ｋｃの適切な変化が完了したと判定して、補正値Ｋｃの変化処理を停止することができる。よって、制御装置５０の処理負荷を軽減す
ることができる。また、補正値Ｋｃの適切な変化が完了した後も、補正値Ｋｃの変化処理を継続することにより、何らかの要因で、補正値Ｋｃの誤変化が生じることを防止できる。

本実施の形態では、補正値変化部５４は、式（８）を用い、判定角度分のＮ個（本例では９０個）の角加速度変化量Δαｄのデータに基づいて、標準偏差σを算出する。

また、補正値変化部５４は、式（９）に示すように、標準偏差σをＭ回算出する毎に、Ｍ個の標準偏差σの平均値σａｖｅを算出する。そして、補正値変化部５４は、前回算出した標準偏差の平均値σａｖｅ（ｋ−１）と今回算出した標準偏差の平均値σａｖｅ（ｋ）との偏差を、今回算出した標準偏差の平均値σａｖｅ（ｋ）で除算した値の絶対値を、標準偏差の変化量の絶対値｜Δσ｜として算出するように構成されている。よって、無次元化された標準偏差の変化量が算出される。

＜記憶データを用いた補正値Ｋｃの変化＞
本実施の形態では、角度情報検出部５１は、内燃機関１への燃料供給が停止されており、内燃機関１の負荷が予め設定された判定負荷よりも低い変化許可条件が成立した場合に検出した、角度区間Ｓｄのそれぞれの角度間隔Δθｄ及び時間間隔ΔＴｄを記憶する。そして、角度情報補正部５２、角度情報算出部５３、及び補正値変化部５４は、記憶された角度区間Ｓｄのそれぞれの角度間隔Δθｄ及び時間間隔ΔＴｄに対して処理を行って、補正値Ｋｃを変化させる。

この構成によれば、検出した角度間隔Δθｄ及び時間間隔ΔＴｄに対してリアルタイムに処理を行って、補正値Ｋｃを変化させる必要がなくなる。そのため、補正値Ｋｃを変化させる処理を、時間的に分散させて実行することができる。よって、制御装置５０の処理負荷を軽減することができると共に、制御装置５０に求められる処理性能を低減することができる。また、燃焼停止及び低負荷の好条件で記憶した角度間隔Δθｄ及び時間間隔ΔＴｄを用いることができるので、補正値Ｋｃの変化精度を向上させることができる。

角度情報検出部５１は、変化許可条件が成立した場合に、予め設定された記憶回転数分の角度間隔Δθｄ及び時間間隔ΔＴｄを、制御装置５０のＲＡＭ等の書き換え可能な記憶装置９１に記憶する。角度情報検出部５１は、クランク軸角度θｄ及び検出時刻Ｔｄを検出し、角度識別番号ｎ、角度間隔Δθｄ及び時間間隔ΔＴｄを算出する毎に、記憶番号ｉをカウントアップし、図１２に示すように、カウントアップした記憶番号ｉと対応付けて、角度識別番号ｎ、角度間隔Δθｄ及び時間間隔ΔＴｄを記憶装置９１に記憶する。

角度情報補正部５２は、記憶された角度間隔Δθｄ又は時間間隔ΔＴｄに対して補正値Ｋｃにより補正する。角度情報算出部５３は、補正後の角度間隔Δθｄｃ及び時間間隔ΔＴｄｃに基づいて、クランク角速度ωｄ、クランク角加速度αｄ、及び角加速度変化量Δαｄ（又はクランク角躍度δｄ）を算出する。補正値変化部５４は、角加速度変化量Δαｄがゼロに近づくように、補正値Ｋｃを変化させる。

角度情報補正部５２、角度情報算出部５３、及び補正値変化部５４は、記憶された記憶回転数分の角度間隔Δθｄ及び時間間隔ΔＴｄを繰り返し用いて補正値Ｋｃを変化させる処理を行う。好条件で記憶したデータを繰り返し用いて、補正値Ｋｃを高精度に変化させることができる。

＜補正値Ｋｃの変化処理のフローチャート＞
次に、図１３に示すフローチャートを用いて、本実施の形態に係る補正値Ｋｃの変化処理を説明する。
まず、ステップＳ０１で、補正値変化部５４は、内燃機関１への燃料供給が停止されており、内燃機関１の負荷が予め設定された判定負荷よりも低い変化許可条件が成立しているか否かを判定する。角度情報検出部５１は、変化許可条件が成立した場合（ステップＳ０１：Ｙｅｓ）は、上記のように、検出した角度間隔Δθｄ及び時間間隔ΔＴｄを記憶する。本例では、記憶回転数分のデータが記憶される。

ステップＳ０３で、補正値変化部５４は、記憶番号ｉを２に設定する。これは、本実施の形態では、角加速度変化量Δαｄの算出に、１つ前の角度間隔Δθｄ（ｉ−１）及び時間間隔ΔＴｄ（ｉ−１）のデータが必要なためである。そして、ステップＳ０４で、角度情報補正部５２は、記憶番号ｉに対応する角度間隔Δθｄ又は時間間隔ΔＴｄに対して、記憶番号ｉに対応する角度識別番号ｎの補正値Ｋｃにより補正する。本実施の形態では、角度情報補正部５２は、角加速度変化量Δαｄを算出するために、１つ前、今回、１つ後、及び２つ後の記憶番号（ｉ−１）、（ｉ）、（ｉ＋１）、（ｉ＋２）に対応する角度間隔Δθｄ又は時間間隔ΔＴｄに対して、補正値Ｋｃにより補正する。

ステップＳ０５で、角度情報算出部５３は、記憶番号ｉに対応する補正後の角度間隔Δθｄｃ及び時間間隔ΔＴｄｃに基づいて、クランク角速度ωｄ、クランク角加速度αｄ、及び角加速度変化量Δαｄ（又はクランク角躍度δｄ）を算出する。本実施の形態では、角度情報算出部５３は、１つ前、今回、１つ後、及び２つ後の記憶番号（ｉ−１）、（ｉ）、（ｉ＋１）、（ｉ＋２）に対応する補正後の角度間隔Δθｄｃ及び時間間隔ΔＴｄｃを用いて、記憶番号ｉに対応する角加速度変化量Δαｄ（ｉ）を算出する。

ステップＳ０６で、補正値変化部５４は、記憶番号ｉに対応する角加速度変化量Δαｄ（ｉ）（又はクランク角躍度δｄ）がゼロに近づくように、記憶番号ｉに対応する角度識別番号ｎの補正値Ｋｃを変化させる。

ステップＳ０７で、補正値変化部５４は、標準偏差σを算出するための、式（８）における、角加速度変化量Δαｄ（又はクランク角躍度δｄ）のΣの計算を行う。補正値変化部５４は、判定角度分のＮ個（本例では９０個）のΣの算出を行った場合は、式（８）の標準偏差σの算出を行う。

ステップＳ０８で、補正値変化部５４は、記憶番号ｉが、記憶回転数分の記憶データ数以上になったか否かを判定する。補正値変化部５４は、記憶番号ｉが記憶データ数未満の場合（ステップＳ０８：Ｎｏ）は、ステップＳ０９で、記憶番号ｉに１を足して、ステップＳ０４に戻り、次の記憶番号ｉの記憶データを用いたステップＳ０４からステップＳ０７の処理を繰り返し行う。

一方、補正値変化部５４は、記憶番号ｉが記憶データ数以上の場合（ステップＳ０８：Ｙｅｓ）は、ステップＳ１０で、式（９）を用い、記憶番号ｉが２から記憶データ数まで増加する間に算出したＭ個の標準偏差σの平均値σａｖｅを算出する。そして、補正値変化部５４は、式（９）を用い、前回算出した標準偏差の平均値σａｖｅ（ｋ−１）及び今回算出した標準偏差の平均値σａｖｅ（ｋ）に基づいて、標準偏差の変化量の絶対値｜Δσ｜を算出する。

ステップＳ１１で、補正値変化部５４は、標準偏差の変化量の絶対値｜Δσ｜が、判定変化量Ｘσ未満になったか否かを判定する。補正値変化部５４は、標準偏差の変化量の絶対値｜Δσ｜が、判定変化量Ｘσ以上である場合（ステップＳ１１：Ｎｏ）は、ステップＳ０３に戻り、記憶された記憶回転数分の角度間隔Δθｄ及び時間間隔ΔＴｄを用いて補正値Ｋｃを変化させる処理を繰り返し行う。一方、補正値変化部５４は、標準偏差の変化量の絶対値｜Δσ｜が、判定変化量Ｘσ未満である場合（ステップＳ１１：Ｙｅｓ）は、ステップＳ１２に進み、補正値Ｋｃの変化が完了したと判定し、補正値Ｋｃの変化処理を終了する。

＜変化完了後の補正値Ｋｃによる補正結果＞
図１４に、変化完了後の補正値Ｋｃによる補正結果を実線で示す。破線で示す、補正値Ｋｃによる補正前の時間間隔ΔＴｄには、リングギア２５の歯の公差により、高振幅の振動成分が重畳している。変化完了後の補正値Ｋｃによる補正後の時間間隔ΔＴｄｃは、補正前よりも振幅が低減しており、図１１に示した補正目標に近づいている。よって、補正値Ｋｃによりリングギア２５の歯の公差による時間間隔ΔＴｄｃの変動を良好に補正できている。

図１５に、図１４の補正前の時間間隔ΔＴｄ、及び変化完了後の補正値Ｋｃによる補正後の時間間隔ΔＴｄｃのスペクトル解析結果を示す。スペクトル解析には、４ｄｅｇに対応する時間間隔ΔＴｄを１点のデータとし、６４点のデータを用いた高速フーリエ変換が用いられている。横軸の周波数の次数は、スペクトル解析結果の標本位置を表している。１次の周波数は、クランク軸２の回転周波数の約０．７倍（≒６４／９０）に相当し、ｎ次の周波数は、クランク軸２の回転周波数のｎ×０．７倍に相当する。補正により、２次以上の周波数成分が大幅に低減している。補正前では、８．５次の低周波数の周波数成分が大きくなっている。そのため、時間間隔ΔＴｄの振動成分を減少させるローパスフィルタ処理を行う場合に、ローパスフィルタ処理のカットオフ周波数を、８．５次以下に設定する必要があり、８．５次以上の高周波数側に設定することができない。一方、補正後の時間間隔ΔＴｄｃでは、２８次付近の高周波数まで周波数成分を大幅に低減できている。よって、補正後の時間間隔ΔＴｄｃの振動成分を減少させるローパスフィルタ処理を行う場合に、ローパスフィルタ処理のカットオフ周波数を、大幅に高めることができる。従って、補正によって、高周波数帯の時間間隔ΔＴｄの変化を精度よく検出することができ、後述する筒内圧力等の推定精度を向上することができる。

＜筒内圧力推定部５５＞
筒内圧力推定部５５は、内燃機関１のピストン、コンロッド及びクランクを含むクランク軸２の回転系の運動方程式を用い、クランク軸角度θｄ及びクランク角加速度αｄに基づいて、燃焼により生じた燃焼ガス圧トルクＴｂを算出し、燃焼ガス圧トルクＴｂ及びクランク軸角度θｄに基づいて燃焼している気筒ｂの筒内圧力Ｐｃｙｌｂを推定する。

本実施の形態では、角度情報補正部５２及び角度情報算出部５３は、補正値Ｋｃの変化処理のために、変化許可条件の成立時に記憶された角度間隔Δθｄ又は時間間隔ΔＴｄに対して行う処理に加えて、リアルタイムに算出された角度間隔Δθｄ又は時間間隔ΔＴｄに対して、リアルタイムに補正値Ｋｃにより補正を行い、リアルタイムにクランク角速度ωｄ及びクランク角加速度αｄを算出するように構成されている。また、角度情報算出部５３は、高周波のノイズ成分を低減するために、補正後の時間間隔ΔＴｄｃ又はクランク角加速度αｄに対してローパスフィルタ処理を行う。

内燃機関１のピストン、コンロッド及びクランクを含むクランク軸２の回転系の運動方程式は、式（１０）で表せる。
ここで、Ｉは、クランク軸２の慣性モーメントであり、Ｐｃｙｌｊは、ｊ番目気筒７の筒内圧力であり、Ｓｐは、ピストン５の頂面の投影面積であり、ｍｐは、ピストン５の質量であり、αｐｊは、ｊ番目気筒７のピストン５の加速度であり、Ｒｊは、ｊ番目気筒７のピストン５に生じた力を、クランク軸２回りのトルクに変換する変換係数であり、Ｔｅｘは、フリクション、補機負荷、及び走行抵抗等の外部からクランク軸２に伝達される外部負荷トルクである。Ｌは、気筒数であり、本実施の形態ではＬ＝３である。また、ｒは、クランクの半径であり、θｄｊは、ｊ番目気筒７のピストン５の上死点を基準にしたクランク軸角度であり、φｊは、ｊ番目気筒７のコンロッドの角度であり、クランク長さとコンロッド長さの比であるコンロッド比とクランク軸角度θｄｊとに基づいて求まる。

筒内圧力推定部５５は、各気筒ｊのクランク軸角度θｄｊに応じて変化するコンロッド９及びクランク３２の幾何学的関係、並びにクランク角加速度αｄに基づいて、各気筒ｊのピストン５の加速度αｐｊを算出する。また、筒内圧力推定部５５は、各気筒ｊのクランク軸角度θｄｊに基づいて、各気筒ｊの変換係数Ｒｊを算出する。

燃焼が行われる圧縮行程の後半及び膨張行程以外の筒内圧力Ｐｃｙｌｊは、吸気マニホールド１２内の圧力、大気圧、クランク軸角度θｄｊに応じた圧力となる。筒内圧力推定部５５は、吸気マニホールド１２内の圧力、大気圧、クランク軸角度θｄｊに基づいて、吸気行程、圧縮行程（後半を除く）又は排気行程となっている各未燃焼気筒ｊの筒内圧力Ｐｃｙｌｕｂｊを推定する。ｂ番目気筒７が圧縮行程の後半及び膨張行程であり、燃焼が行われている場合は、式（１０）を式（１１）のように変形することができる。ここで、Ｐｃｙｌｂは、燃焼気筒ｂの筒内圧力であり、Ｐｃｙｌｕｂｊは、各未燃焼気筒ｊ（ｊ≠ｂ）の筒内圧力である。

燃焼気筒ｂのピストン５が上死点である場合は、式（１１）の右辺の第一項がゼロになるので、式（１１）を外部負荷トルクＴｅｘについて整理すると、式（１２）のようになる。外部負荷トルクＴｅｘは、１サイクルの間、大きく変動しないため、上死点で推定した一定値であると仮定する。

筒内圧力推定部５５は、式（１２）を用い、燃焼気筒ｂのピストン５が上死点である場合の、各未燃焼気筒ｊの筒内圧力Ｐｃｙｌｕｂｊ、ピストン５の加速度αｐｊ、及び変換係数Ｒｊ、並びにクランク角加速度αｄに基づいて外部負荷トルクＴｅｘを推定する。

式（１１）を、燃焼によりクランク軸２に生じた燃焼ガス圧トルクＴｂに相当する「Ｐｃｙｌｂ・Ｓｐ・Ｒｂ」について整理すると式（１３）を得る。

筒内圧力推定部５５は、式（１３）に示すクランク軸２の回転系の運動方程式を用い、クランク軸角度θｄ及びクランク角加速度αｄに基づいて、燃焼ガス圧トルクＴｂを推定する。この際、筒内圧力推定部５５は、上記のように、燃焼気筒ｂのピストンの加速度αｐｂ及び変換係数Ｒｂ、並びに各未燃焼気筒ｊの筒内圧力Ｐｃｙｌｕｂｊ、ピストン５の加速度αｐｊ及び変換係数Ｒｊ、並びに外部負荷トルクＴｅｘを算出する。

そして、筒内圧力推定部５５は、式（１４）に示すように、燃焼ガス圧トルクＴｂを、ピストン５の投影面積Ｓｐ、及び燃焼気筒ｂの変換係数Ｒｂで除算して、燃焼気筒ｂの筒内圧力Ｐｃｙｌｂを算出する。

＜燃焼パラメータ算出部５６＞
燃焼パラメータ算出部５６は、燃焼気筒ｂの筒内圧力Ｐｃｙｌｂに基づいて、熱発生率及び質量燃焼割合ＭＦＢの一方又は双方の燃焼パラメータを算出する。

本実施の形態では、燃焼パラメータ算出部５６は、式（１５）を用い、単位クランク軸角度当たりの熱発生率ｄＱ／ｄθを算出する。ここで、κは、比熱比であり、Ｖｂは、燃焼気筒ｂのシリンダ容積である。燃焼パラメータ算出部５６は、燃焼気筒ｂのクランク軸角度θｄｂ及びコンロッド９及びクランク３２の幾何学的関係に基づいて、シリンダ容積Ｖｂ及び単位クランク軸角度当たりのシリンダ容積変化率ｄＶｂ／ｄθを算出する。

燃焼パラメータ算出部５６は、式（１６）を用い、熱発生率ｄＱ／ｄθを燃焼開始角度θ０からクランク軸角度θｄｂまで積分した瞬時積分値を、燃焼角度区間全体に亘って熱発生率ｄＱ／ｄθを積分した全積分値Ｑ０で除算して、各クランク軸角度θｄｂの質量燃焼割合ＭＦＢを算出する。

＜燃焼制御部５７＞
燃焼制御部５７は、燃焼パラメータに基づいて、点火時期及びＥＧＲ量の一方又は双方を変化させる燃焼制御を行う。本実施の形態では、燃焼制御部５７は、質量燃焼割合ＭＦＢが０．５（５０％）になるクランク軸角度θｄｂ（燃焼中心角度と称す）を判定し、燃焼中心角度が予め設定された目標角度に近づくように、点火時期及びＥＧＲ量の一方又は双方を変化させる。例えば、燃焼制御部５７は、燃焼中心角度が目標角度よりも遅角側である場合は、点火時期を進角側に変化させる、又はＥＧＲバルブ２２の開度を増加させてＥＧＲ量を増加させる。なお、ＥＧＲ量を増加させれば、燃焼速度が緩やかになり、燃焼中心角度が進角側に変化する。一方、燃焼制御部５７は、燃焼中心角度が目標角度よりも進角側である場合は、点火時期を遅角側に変化させる、又はＥＧＲバルブ２２の開度を減少させてＥＧＲ量を減少させる。

或いは、燃焼制御部５７は、熱発生率ｄＱ／ｄθが最大値になるクランク軸角度θｄｂを判定し、当該クランク軸角度θｄｂが予め設定された目標角度に近づくように、点火時期及びＥＧＲ量の一方又は双方を変化させるように構成されてもよい。

＜筒内圧力、熱発生率及び質量燃焼割合の算出精度の向上＞
図１６を用いて、補正値Ｋｃによる補正を行う本実施の形態の場合と、補正を行わない比較例の場合の、筒内圧力、熱発生率及び質量燃焼割合の算出精度を説明する。図１６の左側の列に、補正値Ｋｃによる補正を行わない比較例の場合の各算出値の挙動を示している。図１６の右側の列に、変化完了後の補正値Ｋｃによる補正を行った本実施の形態の場合の各算出値の挙動を示している。

左側の比較例の場合は、上記のように、時間間隔ΔＴｄに重畳している８．５次の周波数成分を除去するために、カットオフ周波数の低いローパスフィルタ処理が行われている。一方、右側の本実施の形態では、ローパスフィルタ処理のカットオフ周波数が比較例の場合よりも、高くされている（本例では１２次）。そのため、左側の比較例の場合は、実測値と比較すると、波形がなまっており、特に熱発生率及び質量燃焼割合ＭＦＢの算出精度が悪くなっている。一方、右側の実施の形態の場合は、実測値と比較しても、波形のな
まりが小さく、特に熱発生率及び質量燃焼割合ＭＦＢの算出精度が大幅に向上している。よって、本実施の形態では、熱発生率及び質量燃焼割合ＭＦＢの燃焼パラメータを用いた燃焼制御の制御精度を向上させることができる。

＜処理全体の概略フローチャート＞
本実施の形態に係る制御装置５０の概略的な処理の手順（内燃機関１の制御方法）について、図１７に示すフローチャートに基づいて説明する。図１７のフローチャートの処理は、演算処理装置９０が記憶装置９１に記憶されたソフトウェア（プログラム）を実行することにより、例えば一定の演算周期毎に繰り返し実行される。

ステップＳ５１で、角度情報検出部５１は、上記のように、特定クランク角センサ６の出力信号に基づいてクランク軸角度θｄを検出すると共に検出時刻Ｔｄを検出し、検出角度θｄ及び検出時刻Ｔｄに基づいて、角度間隔Δθｄ及び時間間隔ΔＴｄを算出する角度情報検出処理（角度情報検出ステップ）を実行する。

次に、ステップＳ５２で、角度情報補正部５２は、上記のように、角度区間Ｓｄのそれぞれの角度間隔Δθｄ又は時間間隔ΔＴｄを、角度区間Ｓｄのそれぞれに対応して１つずつ設けた補正値Ｋｃにより補正する角度情報補正処理（角度情報補正ステップ）を実行する。

ステップＳ５３で、角度情報算出部５３は、上記のように、補正値Ｋｃによる補正後の角度間隔Δθｄｃ及び時間間隔ΔＴｄｃに基づいて、検出角度θｄ又は角度区間Ｓｄのそれぞれに対応する、クランク角速度ωｄ、クランク角加速度αｄ、及び角加速度変化量Δαｄ（又はクランク角躍度δｄ）を算出する角度情報算出処理（角度情報算出ステップ）を実行する。

そして、ステップＳ５４で、補正値変化部５４は、上記のように、検出角度θｄ又は角度区間Ｓｄのそれぞれの角加速度変化量Δαｄ（又はクランク角躍度δｄ）がゼロに近づくように、角度区間Ｓｄのそれぞれの補正値Ｋｃを変化させる補正値変化処理（補正値変化ステップ）を実行する。

また、ステップＳ５５で、筒内圧力推定部５５は、上記のように、内燃機関１のピストン、コンロッド及びクランクを含むクランク軸２の回転系の運動方程式を用い、クランク軸角度θｄ及びクランク角加速度αｄに基づいて、燃焼ガス圧トルクＴｂを算出し、燃焼ガス圧トルクＴｂ及びクランク軸角度θｄに基づいて燃焼気筒ｂの筒内圧力Ｐｃｙｌｂを推定する筒内圧力推定処理（筒内圧力推定ステップ）を実行する。

ステップＳ５６で、燃焼パラメータ算出部５６は、上記のように、燃焼気筒ｂの筒内圧力Ｐｃｙｌｂに基づいて、熱発生率及び質量燃焼割合ＭＦＢの一方又は双方の燃焼パラメータを算出する燃焼パラメータ算出処理（燃焼パラメータ算出ステップ）を実行する。

そして、ステップＳ５７で、燃焼制御部５７は、上記のように、燃焼パラメータに基づいて、点火時期及びＥＧＲ量の一方又は双方を変化させる燃焼制御処理（燃焼制御ステップ）を実行する。

〔その他の実施の形態〕
最後に、本発明のその他の実施の形態について説明する。なお、以下に説明する各実施の形態の構成は、それぞれ単独で適用されるものに限られず、矛盾が生じない限り、他の実施の形態の構成と組み合わせて適用することも可能である。

（１）上記の実施の形態１においては、第２クランク角センサ６が、本発明における「特定クランク角センサ」に相当し、フライホイール２７が、本発明における「回転部材」に相当し、フライホイール２７に設けられたリングギア２５の歯が、本発明における「被検出部」に相当する場合を例に説明した。しかし、本発明の実施の形態はこれに限定されない。すなわち、第１クランク角センサ１１が、本発明における「特定クランク角センサ」に相当し、信号板１０が、本発明における「回転部材」に相当し、信号板１０に設けられた複数の歯が、本発明における「被検出部」に相当してもよい。

（２）上記の実施の形態１においては、内燃機関１は、ガソリンエンジンとされている場合を例として説明した。しかし、本発明の実施の形態はこれに限定されない。すなわち、内燃機関１は、ディーゼルエンジン、ＨＣＣＩ燃焼（Homogeneous-Charge Compression Ignition Combustion）を行うエンジン等の各種の内燃機関とされてもよい。

（３）上記の実施の形態１においては、角度情報算出部５３は、図７を用いて説明したように、対象検出角度θｄ（ｎ）の直前２つ及び直後２つの角度区間Ｓｄのそれぞれのクランク角速度ωｄを算出し、直前２つの角度区間Ｓｄに対応するクランク角加速度αｄ（ｎ−１）及び直後２つの角度区間Ｓｄに対応する１つのクランク角加速度αｄ（ｎ＋１）を算出し、これら２つのクランク角加速度に基づいて、対象検出角度θｄ（ｎ）に対応する１つの角加速度変化量Δαｄ（ｎ）を算出するように構成されている場合を例に説明した。しかし、本発明の実施の形態はこれに限定されない。すなわち、角度情報算出部５３は、補正後の角度間隔Δθｄｃ及び時間間隔ΔＴｄｃに基づいて、検出角度θｄ又は角度区間Ｓｄのそれぞれに対応する、クランク角速度ωｄ、クランク角加速度αｄ、及び角加速度変化量Δαｄを算出すれば、どのような算出方法を用いてもよい。

例えば、図１８に示すように、角度情報算出部５３は、角加速度変化量Δαｄ（ｎ）を算出する角度区間Ｓｄである対象角度区間Ｓｄ（ｎ）、対象角度区間Ｓｄ（ｎ）の直前１つの角度区間Ｓｄ（ｎ−１）、及び対象角度区間Ｓｄ（ｎ）の直後１つの角度区間Ｓｄ（ｎ＋１）のそれぞれの補正後の角度間隔Δθｄｃ及び時間間隔ΔＴｄｃに基づいて、直前１つの角度区間Ｓｄ（ｎ−１）、対象角度区間Ｓｄ（ｎ）、及び直後１つの角度区間Ｓｄ（ｎ＋１）のそれぞれのクランク角速度ωｄ（ｎ−１）、ωｄ（ｎ）、ωｄ（ｎ＋１）を算出する。そして、角度情報算出部５３は、直前１つの角度区間Ｓｄ（ｎ−１）及び対象角度区間Ｓｄ（ｎ）のそれぞれのクランク角速度ωｄ（ｎ−１）、ωｄ（ｎ）及び時間間隔ΔＴｄ（ｎ−１）、ΔＴｄ（ｎ）に基づいて、当該２つの角度区間Ｓｄに対応する１つのクランク角加速度αｄ（ｎ−１）を算出する。また、角度情報算出部５３は、対象検出角度θｄ（ｎ）及び直後１つの角度区間Ｓｄ（ｎ＋１）のそれぞれのクランク角速度ωｄ（ｎ）、ωｄ（ｎ＋１）及び時間間隔ΔＴｄ（ｎ）、ΔＴｄ（ｎ＋１）に基づいて、当該２つの角度区間Ｓｄに対応する１つのクランク角加速度αｄ（ｎ）を算出する。そして、角度情報算出部５３は、２つのクランク角加速度αｄ（ｎ−１）、αｄ（ｎ）に基づいて、対象角度区間Ｓｄ（ｎ）に対応する１つの角加速度変化量Δαｄ（ｎ）を算出するように構成されてもよい。

（４）上記の実施の形態１においては、角度情報算出部５３は、検出角度θｄのそれぞれに対応する角加速度変化量Δαｄを算出し、補正値変化部５４は、対象検出角度θｄ（ｎ）の角加速度変化量Δαｄ（ｎ）がゼロに近づくように、対象検出角度θｄ（ｎ）の直前１つの対象角度区間Ｓｄ（ｎ）の補正値Ｋｃ（ｎ）及び直後１つの対象角度区間Ｓｄ（ｎ＋１）の補正値Ｋｃ（ｎ＋１）の一方又は双方を変化させるように構成されている場合を例に説明した。しかし、本発明の実施の形態はこれに限定されない。すなわち、角度情報算出部５３は、図１８を用いて説明したように、角度区間Ｓｄのそれぞれに対応する角加速度変化量Δαｄを算出し、補正値変化部５４は、処理対象とする角度区間である対象角度区間Ｓｄ（ｎ）の角加速度変化量Δαｄ（ｎ）がゼロに近づくように、対象角度区間Ｓ
ｄ（ｎ）の補正値Ｋｃ（ｎ）を変化させるように構成されてもよい。

具体的には、補正値変化部５４は、補正値Ｋｃにより角度間隔Δθｄを補正する場合であって、対象角度区間Ｓｄ（ｎ）の角加速度変化量Δαｄ（ｎ）が正の値である場合は、対象角度区間Ｓｄ（ｎ）の補正後の角度間隔Δθｄｃ（ｎ）が増加するように補正値Ｋｃ（ｎ）を変化させる増加補正を実行し、対象角度区間Ｓｄ（ｎ）の角加速度変化量Δαｄ（ｎ）が負の値である場合は、対象角度区間Ｓｄ（ｎ）の補正後の角度間隔Δθｄｃ（ｎ）が減少するように補正値Ｋｃ（ｎ）を変化させる減少補正を実行する。或いは、補正値変化部５４は、補正値Ｋｃにより時間間隔ΔＴｄを補正する場合であって、対象角度区間Ｓｄ（ｎ）の角加速度変化量Δαｄ（ｎ）が正の値である場合は、対象角度区間Ｓｄ（ｎ）の補正後の時間間隔ΔＴｄｃ（ｎ）が減少するように補正値Ｋｃ（ｎ）を変化させる減少補正を実行し、対象角度区間Ｓｄ（ｎ）の角加速度変化量Δαｄ（ｎ）が負の値である場合は、対象角度区間Ｓｄ（ｎ）の補正後の時間間隔ΔＴｄｃ（ｎ）が増加するように補正値Ｋｃを変化させる増加補正を実行する。

（５）上記の実施の形態１においては、補正値変化部５４は、内燃機関１への燃料供給が停止されており、内燃機関１の負荷が判定負荷よりも低い変化許可条件が成立した場合に、補正値Ｋｃの変化を許可するように構成されている場合を例に説明した。しかし、本発明の実施の形態はこれに限定されない。すなわち、変化許可条件は、実施の形態１とは異なる適切な条件とされてもよく、或いは、変化許可条件が設けられなくてもよい。

（６）上記の実施の形態１においては、角度情報補正部５２、角度情報算出部５３、及び補正値変化部５４は、変化許可条件が成立した場合に記憶された角度区間Ｓｄのそれぞれの角度間隔Δθｄ及び時間間隔ΔＴｄに対して処理を行って、補正値Ｋｃを変化させるように構成されている場合を例に説明した。しかし、本発明の実施の形態はこれに限定されない。すなわち、角度情報補正部５２、角度情報算出部５３、及び補正値変化部５４は、変化許可条件が成立している場合に、リアルタイムに算出された角度区間Ｓｄのそれぞれの角度間隔Δθｄ又は時間間隔ΔＴｄに対して、リアルタイムに処理を行って、補正値Ｋｃを変化させるように構成されてもよい。

（７）上記の実施の形態１においては、補正値変化部５４は、判定角度間の角加速度変化量Δαｄの標準偏差σを算出し、標準偏差σの前回算出値と今回算出値との変化量の絶対値｜Δσ｜が判定変化量未満になった場合に、補正値Ｋｃの変化を停止し、補正値Ｋｃを保持するように構成されている場合を例に説明した。しかし、本発明の実施の形態はこれに限定されない。すなわち、補正値変化部５４は、補正値Ｋｃの変化を停止しないように構成されてもよく、或いは、標準偏差σの変化量の絶対値｜Δσ｜以外の条件、例えば、標準偏差σが予め設定された判定偏差以下になった場合に、補正値Ｋｃの変化を停止し、補正値Ｋｃを保持するように構成されてもよい。

（８）上記の実施の形態１においては、制御装置５０は、補正値Ｋｃによる補正後の角度間隔Δθｄｃ及び時間間隔ΔＴｄｃに基づいて算出されたクランク角速度ωｄ及びクランク角加速度αｄに基づいて、筒内圧力、熱発生率及び質量燃焼割合を算出し、燃焼制御を行うように構成されている場合を例に説明した。しかし、本発明の実施の形態はこれに限定されない。すなわち、制御装置５０は、補正後の角度間隔Δθｄｃ及び時間間隔ΔＴｄｃに基づいて算出されたクランク角速度ωｄ及びクランク角加速度αｄに基づいて、各気筒７の燃焼の失火検出等の他の制御を行うように構成されてもよい。

（９）上記の実施の形態１においては、補正値変化部５４は、ばらつき度合として、判定角度間の角加速度変化量Δαｄの標準偏差σを算出するように構成されている場合を例に説明した。しかし、本発明の実施の形態はこれに限定されない。すなわち、補正値変化部
５４は、ばらつき度合として、判定角度間の角加速度変化量Δαｄの分散σ＾２を算出するように構成されてもよい。補正値変化部５４は、式（８）の標準偏差σを２乗したものに相当する式（１７）を用いて分散σ＾２を算出する。

なお、本発明は、その発明の範囲内において、実施の形態を適宜、変形、省略したりすることが可能である。

この発明は、クランク軸と一体回転する回転部材に、予め定められた複数のクランク軸角度に設けられた複数の被検出部と、非回転部材に固定され、前記被検出部を検出する特定クランク角センサと、を備えた内燃機関の制御装置及び制御方法に好適に利用することができる。

１：内燃機関、２：クランク軸、６：第２クランク角センサ（特定クランク角センサ）、７：気筒、２４：エンジンブロック（非回転部材）、２５：リングギア（被検出部）、２７：フライホイール（回転部材）、５０：制御装置、５１：角度情報検出部、５２：角度情報補正部、５３：角度情報算出部、５４：補正値変化部、５５：筒内圧力推定部、５６：燃焼パラメータ算出部、５７：燃焼制御部、Δαｄ：角加速度変化量、Δθｄ：角度間隔、Δθｄｃ：補正後の角度間隔、ΔＴｄ：時間間隔、ΔＴｄｃ：補正後の時間間隔、αｄ：クランク角加速度、δｄ：クランク角躍度、θｄ：クランク軸角度（検出角度）、ωｄ：クランク角速度、Ｋｃ：補正値、ＭＦＢ：質量燃焼割合、Ｐｃｙｌｂ：燃焼気筒の筒内圧力、Ｔｂ：燃焼ガス圧トルク、Ｔｄ：検出時刻、ｎ：角度識別番号

Claims (12)

1. クランク軸と一体回転する回転部材に、予め定められた複数のクランク軸角度に設けられた複数の被検出部と、非回転部材に固定され、前記被検出部を検出する特定クランク角センサと、を備えた内燃機関の制御装置であって、
   前記特定クランク角センサの出力信号に基づいてクランク軸角度を検出すると共に前記クランク軸角度を検出した検出時刻を検出し、検出した前記クランク軸角度である検出角度及び前記検出時刻に基づいて、前記検出角度の間の角度区間に対応する角度間隔及び時間間隔を算出する角度情報検出部と、
   前記角度区間のそれぞれの前記角度間隔又は前記時間間隔を、前記角度区間のそれぞれに対応して１つずつ設けた補正値により補正する角度情報補正部と、
   前記角度区間のそれぞれの前記補正値による補正後の前記角度間隔及び前記時間間隔に基づいて、前記検出角度又は前記角度区間のそれぞれに対応する、前記クランク軸角度の時間変化率であるクランク角速度、前記クランク角速度の時間変化率であるクランク角加速度、及び前記クランク角加速度の時間変化量である角加速度変化量を算出する角度情報算出部と、
   前記検出角度又は前記角度区間のそれぞれの前記角加速度変化量がゼロに近づくように、前記角度区間のそれぞれの前記補正値を変化させる補正値変化部と、を備える内燃機関の制御装置。
2. 前記補正値変化部は、前記内燃機関への燃料供給が停止されており、前記内燃機関の負荷が予め設定された判定負荷よりも低い場合に、前記補正値の変化を許可する請求項１に記載の内燃機関の制御装置。
3. 前記角度情報検出部は、前記内燃機関への燃料供給が停止されており、前記内燃機関の負荷が予め設定された判定負荷よりも低い場合に検出した前記角度区間のそれぞれの前記角度間隔及び前記時間間隔を記憶し、
   前記角度情報補正部、前記角度情報算出部、及び前記補正値変化部は、記憶された前記角度区間のそれぞれの前記角度間隔及び前記時間間隔に対して処理を行って、前記補正値を変化させる請求項１又は２に記載の内燃機関の制御装置。
4. 前記補正値変化部は、前記内燃機関のスロットルバルブの開度が予め設定された判定開度よりも低い場合に成立する条件、前記内燃機関に吸入される吸入空気量が予め設定された判定空気量よりも低い場合に成立する条件、前記内燃機関の吸気マニホールド内の圧力が予め設定された判定圧力よりも低い場合に成立する条件、前記内燃機関の充填効率が予め設定された判定効率よりも低い場合に成立する条件のいずれか１つ以上の条件を論理積により組み合わせた条件が成立した場合に、前記内燃機関の負荷が前記判定負荷よりも低いと判定する請求項２又は３に記載の内燃機関の制御装置。
5. 前記角度情報算出部は、前記角加速度変化量を算出する前記検出角度である対象検出角度の直前２つの前記角度区間及び直後２つの前記角度区間のそれぞれの補正後の前記角度間隔及び前記時間間隔に基づいて、当該直前２つ及び直後２つの前記角度区間のそれぞれの前記クランク角速度を算出し、前記直前２つの前記角度区間のそれぞれの前記クランク角速度及び前記時間間隔に基づいて、当該直前２つの前記角度区間に対応する１つの前記クランク角加速度を算出し、前記直後２つの前記角度区間のそれぞれの前記クランク角速度及び前記時間間隔に基づいて、当該直後２つの前記角度区間に対応する１つの前記クランク角加速度を算出し、前記直前２つの前記角度区間に対応する前記クランク角加速度及び前記直後２つの前記角度区間に対応する前記クランク角加速度に基づいて、前記対象検出角度に対応する１つの前記角加速度変化量を算出する請求項１から４のいずれか一項に記載の内燃機関の制御装置。
6. 前記角度情報算出部は、前記検出角度のそれぞれに対応する前記角加速度変化量を算出し、
   前記補正値変化部は、前記補正値により前記角度間隔を補正する場合であって、処理対象とする前記検出角度である対象検出角度の前記角加速度変化量が正の値である場合は、前記対象検出角度の直前１つの前記角度区間の補正後の前記角度間隔が増加するように、当該角度区間の前記補正値を変化させる増加補正、及び前記対象検出角度の直後１つの前記角度区間の補正後の前記角度間隔が減少するように、当該角度区間の前記補正値を変化させる減少補正の一方又は双方を実行し、前記対象検出角度の前記角加速度変化量が負の値である場合は、前記対象検出角度の直前１つの前記角度区間の補正後の前記角度間隔が減少するように、当該角度区間の前記補正値を変化させる減少補正、及び前記対象検出角度の直後１つの前記角度区間の補正後の前記角度間隔が増加するように、当該角度区間の前記補正値を変化させる増加補正の一方又は双方を実行し、或いは、
   前記補正値により前記時間間隔を補正する場合であって、前記対象検出角度の前記角加速度変化量が正の値である場合は、前記対象検出角度の直前１つの前記角度区間の補正後の前記時間間隔が減少するように、当該角度区間の前記補正値を変化させる減少補正、及び前記対象検出角度の直後１つの前記角度区間の補正後の前記時間間隔が増加するように、当該角度区間の前記補正値を変化させる増加補正の一方又は双方を実行し、前記対象検出角度の前記角加速度変化量が負の値である場合は、前記対象検出角度の直前１つの前記角度区間の補正後の前記時間間隔が増加するように、当該角度区間の前記補正値を変化させる増加補正、及び前記対象検出角度の直後１つの前記角度区間の補正後の前記時間間隔が減少するように、当該角度区間の前記補正値を変化させる減少補正の一方又は双方を実行する請求項１から５のいずれか一項に記載の内燃機関の制御装置。
7. 前記角度情報算出部は、前記角度区間のそれぞれに対応する前記角加速度変化量を算出し、
   前記補正値変化部は、前記補正値により前記角度間隔を補正する場合であって、処理対象とする前記角度区間である対象角度区間の前記角加速度変化量が正の値である場合は、前記対象角度区間の補正後の前記角度間隔が増加するように前記補正値を変化させる増加補正を実行し、前記対象角度区間の前記角加速度変化量が負の値である場合は、前記対象角度区間の補正後の前記角度間隔が減少するように前記補正値を変化させる減少補正を実行し、或いは、
   前記補正値により前記時間間隔を補正する場合であって、前記対象角度区間の前記角加速度変化量が正の値である場合は、前記対象角度区間の補正後の前記時間間隔が減少するように前記補正値を変化させる減少補正を実行し、前記対象角度区間の前記角加速度変化量が負の値である場合は、前記対象角度区間の補正後の前記時間間隔が増加するように前記補正値を変化させる増加補正を実行する請求項１から５のいずれか一項に記載の内燃機関の制御装置。
8. 前記補正値変化部は、前記クランク軸が予め設定された判定角度回転する毎に、前記判定角度間の前記角加速度変化量についての、当該角加速度変化量の平均値に対するばらつき度合を算出し、前回算出した前記ばらつき度合と今回算出した前記ばらつき度合との変化量の絶対値が、予め設定された判定変化量未満になった場合に、前記補正値の変化を停止し、前記補正値を保持する請求項１から７のいずれか一項に記載の内燃機関の制御装置。
9. 前記角度情報算出部は、前記角加速度変化量として、前記クランク角加速度の時間変化率であるクランク角躍度を算出する請求項１から８のいずれか一項に記載の内燃機関の制御装置。
10. 前記内燃機関のピストン、コンロッド及びクランクを含む前記クランク軸の回転系の運動方程式を用い、前記クランク軸角度及び前記クランク角加速度に基づいて、燃焼により生じたガス圧トルクを算出し、前記ガス圧トルク及び前記クランク軸角度に基づいて燃焼している気筒の筒内圧力を推定する筒内圧力推定部を備える請求項１から９のいずれか一項に記載の内燃機関の制御装置。
11. 前記筒内圧力に基づいて、熱発生率及び質量燃焼割合の一方又は双方の燃焼パラメータを推定する燃焼パラメータ算出部と、
    前記燃焼パラメータに基づいて、点火時期及びＥＧＲ量の一方又は双方を変化させる燃焼制御部と、を備えた請求項１０に記載の内燃機関の制御装置。
12. クランク軸と一体回転する回転部材に、予め定められた複数のクランク軸角度に設けられた複数の被検出部と、非回転部材に固定され、前記被検出部を検出する特定クランク角センサと、を備えた内燃機関の制御方法であって、
    前記特定クランク角センサの出力信号に基づいてクランク軸角度を検出すると共に前記クランク軸角度を検出した検出時刻を検出し、検出した前記クランク軸角度である検出角度及び前記検出時刻に基づいて、前記検出角度の間の角度区間に対応する角度間隔及び時間間隔を算出する角度情報検出ステップと、
    前記角度区間のそれぞれの前記角度間隔又は前記時間間隔を、前記角度区間のそれぞれに対応して１つずつ設けた補正値により補正する角度情報補正ステップと、
    前記角度区間のそれぞれの前記補正値による補正後の前記角度間隔及び前記時間間隔に基づいて、前記検出角度又は前記角度区間のそれぞれに対応する、前記クランク軸角度の時間変化率であるクランク角速度、前記クランク角速度の時間変化率であるクランク角加速度、及び前記クランク角加速度の時間変化量である角加速度変化量を算出する角度情報算出ステップと、
    前記検出角度又は前記角度区間のそれぞれの前記角加速度変化量がゼロに近づくように、前記角度区間のそれぞれの前記補正値を変化させる補正値変化ステップと、を実行する内燃機関の制御方法。