(57)【要約】【課題】筒内相対圧力センサの出力に基いて、シリンダ内のガスの状態を精度良く取得し得る内燃機関の筒内ガス状態取得装置を提供すること。【解決手段】この筒内ガス状態取得装置は、筒内相対圧力センサから得られる圧縮行程中のＮ個の筒内相対圧力検出値Ｐｍｅａｓ（Ｈ）　（Ｈ＝１，・・・，Ｎ）に対し、圧縮行程開始時点での筒内ガス状態を仮設定することで得られる同時点での仮の筒内絶対圧力値Ｐ０から筒内相対圧力検出値Ｐｍｅａｓ（１）を減算して得られる仮の圧力補正値ΔＰを加算して、Ｎ個の仮の筒内絶対圧力設定値Ｐｃ（Ｈ）を設定する。一方、前記仮設定された前記時点での筒内ガス状態を初期条件としてエネルギー保存則に基くシリンダモデルによりＮ個の仮の筒内絶対圧力計算値Ｐｃｍ（Ｈ）を計算する。そして、Ｎ個の仮の筒内絶対圧力設定値Ｐｃ（Ｈ）とＮ個の仮の筒内絶対圧力計算値Ｐｃｍ（Ｈ）の偏差が最小となるように仮設定した前記時点での筒内ガス状態を同定する。【選択図】　図４
2004108348.tif 000002
【特許請求の範囲】
【請求項１】
シリンダ内の圧力を所定の基準圧力からの相対圧力である筒内相対圧力検出値として検出する筒内相対圧力センサを備え、前記筒内相対圧力検出値に基いて前記シリンダ内に吸入されているガスの状態を取得する内燃機関の筒内ガス状態取得装置であって、
圧縮行程時に前記シリンダ内に吸入されている筒内吸入ガス量、同圧縮行程の所定の時点における同シリンダ内の筒内ガス温度、及び同所定の時点における同シリンダ内の筒内絶対圧力のうちの２つの値を仮設定することで同所定の時点における同シリンダ内に吸入されているガスの状態を仮設定する仮設定手段と、
前記仮設定された前記所定の時点におけるシリンダ内に吸入されているガスの状態から得られる同所定の時点におけるシリンダ内の仮の筒内絶対圧力の値と前記所定の時点における筒内相対圧力検出値との比較により仮の圧力補正値を算出し、前記所定の時点を含む前記圧縮行程の所定の区間内における複数の筒内相対圧力検出値と前記仮の圧力補正値とから同所定の区間内における同シリンダ内の複数の仮の筒内絶対圧力設定値を順次設定する第１計算手段と、
前記仮設定手段により仮設定された前記所定の時点におけるシリンダ内のガスの状態を初期条件として、エネルギー保存則に基いて求められた同シリンダについてのモデルを使用して前記所定の区間内におけるシリンダ内の複数の仮の筒内絶対圧力計算値を順次計算する第２計算手段と、
前記所定の区間内における前記複数の仮の筒内絶対圧力設定値と前記複数の仮の筒内絶対圧力計算値との偏差の程度を示す偏差指標値を計算する偏差指標値計算手段と、
前記偏差指標値が示す前記偏差の程度が小さくなるように前記仮設定する２つの値を逐次補正して、前記仮設定手段、前記第１計算手段、前記第２計算手段、及び前記偏差指標値計算手段が行う一連の各処理を、同偏差の程度が略最小になるまで繰返し実行させるとともに、同偏差の程度が略最小になった時点で前記仮設定されている前記２つの値をそれぞれ真の値として設定し、前記所定の時点におけるシリンダ内に吸入されているガスの真の状態を取得する筒内ガス状態取得手段と、
を備えた内燃機関の筒内ガス状態取得装置。
【請求項２】
請求項１に記載の内燃機関の筒内ガス状態取得装置であって、
前記筒内ガス状態取得手段は、前記所定の時点におけるシリンダ内の真の筒内絶対圧力を取得するように構成され、
前記真の筒内絶対圧力の値と前記所定の時点における筒内相対圧力検出値との比較により真の圧力補正値を取得する圧力補正値取得手段を備えた内燃機関の筒内ガス状態取得装置。
【請求項３】
請求項１又は請求項２に記載の内燃機関の筒内ガス状態取得装置において、
前記所定の時点は、吸気弁が閉弁する圧縮行程開始時点に設定されるとともに、
前記所定の区間は、その始期が前記圧縮行程開始時点に、その終期が圧縮行程終了時点になるように設定された内燃機関の筒内ガス状態取得装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、内燃機関のシリンダ内の圧力を相対圧力として検出する筒内相対圧力センサの出力値に基いて圧縮行程時にシリンダ内に吸入されている筒内吸入ガス量等のシリンダ内に吸入されているガスの状態を取得する内燃機関の筒内ガス状態取得装置に関する。
【０００２】
【従来の技術】
内燃機関により燃焼される混合気の空燃比を所定の値とするためには、同内燃機関のシリンダ（気筒、燃焼室）内に吸入される空気の量（以下、「筒内吸入空気量」と称呼する。）を精度良く求める必要がある。このため、例えば、下記特許文献１に開示された内燃機関の筒内吸入空気量検出装置は、シリンダ内の筒内絶対圧力を検出する筒内絶対圧力センサと、シリンダ内に吸入される吸入空気の温度を検出する吸気温度センサと、シリンダ内から排出される排気ガスの温度を検出する排気温度センサとを備え、吸気温度センサの出力及び排気温度センサの出力等に基いて圧縮行程開始時点におけるシリンダ内に吸入されているガスの筒内ガス温度を推定するとともに、同推定された圧縮行程開始時点での筒内ガス温度と、筒内絶対圧力センサにより得られる圧縮行程開始時点での筒内絶対圧力、及び状態方程式等を使用して筒内吸入空気量を検出するようになっている。
【０００３】
ところで、上記したような絶対圧力を検出できる絶対圧力センサは、一般に、比較的狭い圧力範囲内では精度良く絶対圧力を検出できるものの、シリンダ内の圧力等のように比較的広い範囲内で変動する絶対圧力を同広い範囲の全域に渡って精度良く検出することができない。従って、シリンダ内の圧力を広範囲に渡って精度良く検出するためには、一般に、広範囲に渡って精度良く相対圧力（基準圧力からの相対圧力）を検出できる相対圧力センサを使用することが好適である。ところが、シリンダ内の圧力をかかる相対圧力センサにより検出する場合、同相対圧力センサにより得られるシリンダ内の筒内相対圧力の値を、上記筒内絶対圧力の値のように直接状態方程式に適用して筒内吸入空気量を求めることができない。
【０００４】
このため、下記特許文献２に開示された内燃機関の筒内吸入空気量検出装置は、圧縮行程時、シリンダ内のガスの状態はポリトロープ指数が一定のポリトロープ変化（可逆変化）すると仮定した上で、圧縮行程時における異なる２点での各（相対）圧力を筒内相対圧力センサにより検出して各圧力間の圧力差を求め、この圧力差の値と、エンジン回転速度と、同圧力差の値及びエンジン回転速度と筒内吸入空気量との関係を規定する実験的に求められたテーブルと、に基いて筒内吸入空気量を検出するようになっている。
【０００５】
【特許文献１】
特開平１１−１６６４４７号公報
【特許文献２】
特開平２−２３８１４９号公報
【０００６】
【発明が解決しようとする課題】
しかしながら、実際の内燃機関の運転状態は時々刻々と変化するので、圧縮行程時におけるシリンダと同シリンダ外部との間で伝達される熱量も時々刻々と変化し、その結果、実際の圧縮行程時におけるシリンダ内のガスの状態は、常にポリトロープ指数が一定のポリトロープ変化するとは限らない（非可逆変化する場合もある）。従って、上記特許文献２に開示された装置では、筒内吸入空気量が常に精度良く検出されるとは限らない。
【０００７】
さらには、上記特許文献２に開示された装置では、圧縮行程時における２点のみの筒内相対圧力値に基いて筒内吸入空気量を推定しているので、検出点数の少なさに起因する検出精度の低下により、上記圧力差の検出精度が低下して筒内吸入空気量が精度良く検出されない場合もある。
【０００８】
また、内燃機関の異常燃焼（例えば、ノッキング、ミスファイア等）の発生を検出するためには、シリンダ内の筒内絶対圧力を精度良く求める必要がある。ここで、シリンダ内の筒内吸入空気量を精度良く検出できれば、圧縮行程時におけるシリンダ内のガスに状態方程式を適用することで、シリンダ内の筒内絶対圧力をも精度良く計算でき、この計算された筒内絶対圧力値と、筒内相対圧力センサにより得られる筒内相対圧力値との差を圧力補正値（校正値）として設定することで、その後において、筒内相対圧力センサの出力値と同圧力補正値とに基いて筒内絶対圧力を精度良く求めることができるようになる。
【０００９】
ところが、上述したように、上記特許文献２に開示された装置では、筒内吸入空気量が精度良く検出されない場合もあるので、筒内絶対圧力をも精度良く求めることができない。即ち、上記特許文献２に開示された装置では、筒内相対圧力センサの出力値に基いて筒内吸入空気量（ガス量）、筒内絶対圧力等の圧縮行程時にシリンダ内に吸入されているガスの状態を精度良く取得することができないという問題がある。
【００１０】
従って、本発明の目的は、シリンダ内の相対圧力を検出できる筒内相対圧力センサの出力に基いて、筒内吸入ガス量（空気量）、筒内絶対圧力等のシリンダ内に吸入されているガスの状態を精度良く取得することができる内燃機関の筒内ガス状態取得装置を提供することにある。
【００１１】
【本発明の概要】
本発明の特徴は、シリンダ内の圧力を所定の基準圧力からの相対圧力である筒内相対圧力検出値として検出する筒内相対圧力センサを備え、前記筒内相対圧力検出値に基いて前記シリンダ内に吸入されているガスの状態を取得する内燃機関の筒内ガス状態取得装置が、圧縮行程時に前記シリンダ内に吸入されている筒内吸入ガス量、同圧縮行程の所定の時点における同シリンダ内の筒内ガス温度、及び同所定の時点における同シリンダ内の筒内絶対圧力のうちの２つの値を仮設定することで同所定の時点における同シリンダ内に吸入されているガスの状態を仮設定する仮設定手段と、前記仮設定された前記所定の時点におけるシリンダ内に吸入されているガスの状態から得られる同所定の時点におけるシリンダ内の仮の筒内絶対圧力の値と前記所定の時点における筒内相対圧力検出値との比較により仮の圧力補正値を算出し、前記所定の時点を含む前記圧縮行程の所定の区間内における複数の筒内相対圧力検出値と前記仮の圧力補正値とから同所定の区間内における同シリンダ内の複数の仮の筒内絶対圧力設定値を順次設定する第１計算手段と、前記仮設定手段により仮設定された前記所定の時点におけるシリンダ内のガスの状態を初期条件として、エネルギー保存則に基いて求められた同シリンダについてのモデルを使用して前記所定の区間内におけるシリンダ内の複数の仮の筒内絶対圧力計算値を順次計算する第２計算手段と、前記所定の区間内における前記複数の仮の筒内絶対圧力設定値と前記複数の仮の筒内絶対圧力計算値との偏差の程度を示す偏差指標値を計算する偏差指標値計算手段と、前記偏差指標値が示す前記偏差の程度が小さくなるように前記仮設定する２つの値を逐次補正して、前記仮設定手段、前記第１計算手段、前記第２計算手段、及び前記偏差指標値計算手段が行う一連の各処理を、同偏差の程度が略最小になるまで繰返し実行させるとともに、同偏差の程度が略最小になった時点で前記仮設定されている前記２つの値をそれぞれ真の値として設定し、前記所定の時点におけるシリンダ内に吸入されているガスの真の状態を取得する筒内ガス状態取得手段とを備えたことにある。
【００１２】
ここにおいて、前記仮設定手段は、圧縮行程時にシリンダ内に吸入されている筒内吸入ガス量、同圧縮行程の所定の時点における同シリンダ内の筒内ガス温度、及び同所定の時点における同シリンダ内の筒内絶対圧力（以下、「３つの値」と称呼する。）のうちの２つの値を仮設定することで、同仮設定された２つの値と、前記所定の時点での既知のピストン位置から求まる同所定の時点での既知のシリンダ容積と、圧縮行程時にシリンダ内に密閉されているガスに適用され得る状態方程式とに基いて残りの一つの値を算出し、仮設定することができる。
【００１３】
従って、前記第１計算手段が仮の圧力補正値を算出する際に使用する所定の時点におけるシリンダ内の仮の筒内絶対圧力の値は、所定の時点におけるシリンダ内の筒内絶対圧力の値が上記仮設定された２つの値のうちの一つとして仮設定されている場合には、その仮設定された値であり、所定の時点におけるシリンダ内の筒内絶対圧力の値が上記仮設定された２つの値に含まれていない場合には、同仮設定された２つの値と上記状態方程式等に基いて算出され、仮設定された値である。
【００１４】
第１計算手段は、このようにして仮設定された所定の時点における仮の筒内絶対圧力の値と、筒内相対圧力センサにより得られる同所定の時点における筒内相対圧力検出値との比較により仮の圧力補正値を計算する。そして、第１計算手段は、広範囲に渡って精度良く相対圧力（所定の基準圧力からの相対圧力）を検出できる前記筒内相対圧力センサにより得られる圧縮行程の所定の区間内における複数の筒内相対圧力検出値と、上記のように計算した仮の圧力補正値とに基いて、同所定の区間における複数の仮の筒内絶対圧力設定値を順次計算し設定することにより、同所定の区間内において、同複数の仮の筒内絶対圧力設定値を通る仮の筒内絶対圧力設定値の波形を求める。
【００１５】
また、第２計算手段は、圧縮行程時におけるシリンダと同シリンダ外部との間で伝達される熱量も考慮し得るとともに非可逆変化をも扱うことが可能なエネルギー保存則に基いて求められたシリンダ内についてのモデルを使用して、前記仮設定手段により仮設定された所定の時点におけるシリンダ内のガスの状態を初期条件として、所定の区間における複数の仮の筒内絶対圧力計算値を順次計算することにより、同所定の区間内において、同複数の仮の筒内絶対圧力計算値を通る仮の筒内絶対圧力計算値の波形を求める。
【００１６】
偏差指標値計算手段は、上記のように求められた所定の区間における複数の仮の筒内絶対圧力設定値と複数の仮の筒内絶対圧力計算値との偏差の程度を示す偏差指標値を計算する。この偏差指標値は、例えば、所定の区間内における同一クランク角度における仮の筒内絶対圧力設定値と仮の筒内絶対圧力計算値の差を２乗した値を同所定区間内に渡り積算した値であって、これに限定されない。
【００１７】
そして、筒内ガス状態取得手段は、偏差指標値が略最小になるように、仮設定手段により仮設定される上記２つの値を同定し、同定された２つの値をそれぞれ真の値として設定する。また、この筒内ガス状態取得手段は、上記２つの真の値をそれぞれ設定することで、同設定された２つの真の値と、所定の時点での既知のシリンダ容積と、圧縮行程時にシリンダ内に密閉されているガスに適用され得る状態方程式とに基いて上記３つの値のうちの残りの一つの真の値を算出することができる。具体的には、圧縮行程時にシリンダ内に吸入されている真の筒内吸入ガス量、同圧縮行程の所定の時点における同シリンダ内の真の筒内ガス温度、及び同所定の時点における同シリンダ内の真の筒内絶対圧力をそれぞれ求めることができる。
【００１８】
このようにして筒内ガス状態取得手段により求められた３つの真の値は、広範囲に渡って精度良く得られる筒内相対圧力検出値（実測値）に基いた仮の筒内絶対圧力設定値の波形と、圧縮行程時におけるシリンダと同シリンダ外部との間で伝達される熱量も考慮し得るとともに非可逆変化をも扱うことが可能であって経験則ではなく物理法則に従って表されたモデルにより得られる計算精度の高い仮の筒内絶対圧力計算値の波形とが略一致するように同定された値である。
【００１９】
また、筒内相対圧力検出値（実測値）に基いた筒内絶対圧力設定値の波形は、圧縮行程の所定の区間における多数の検出値（実測値）に基いて求められ得るので、検出点数の少なさに起因する検出精度の低下の影響が発生しにくい波形である。従って、本発明による筒内ガス状態取得装置によれば、筒内相対圧力センサの出力に基いて、上記３つの真の値、即ち、これらの値に基く上記所定の時点におけるシリンダ内のガスの真の状態を精度良く取得することができる。
【００２０】
また、これにより、例えば、流体力学等に基づく式により表される吸気系のモデルと、所定のテーブル（マップ）を用いて筒内吸入空気量を推定するように構成された筒内吸入空気量推定装置を有する内燃機関に対して上記本発明による筒内ガス状態取得装置を適用すれば、通常は、上記モデルとテーブルとに基いて筒内吸入空気量を推定するとともに、所定のタイミング毎に本発明による筒内ガス状態取得装置により真の筒内吸入ガス量を求め、同真の筒内吸入ガス量と推定された筒内吸入空気量との比較により上記テーブルを補正することで、筒内吸入空気量推定装置による筒内吸入空気量の推定精度を向上させることができる。
【００２１】
上記本発明による内燃機関の筒内ガス状態取得装置においては、前記筒内ガス状態取得手段は、前記所定の時点におけるシリンダ内の真の筒内絶対圧力を取得するように構成され、前記真の筒内絶対圧力の値と前記所定の時点における筒内相対圧力検出値との比較により真の圧力補正値を取得する圧力補正値取得手段を備えるように構成されることが好適である。
【００２２】
これによれば、圧力補正値取得手段により真の圧力補正値が取得された後において、筒内相対圧力センサの出力値である筒内相対圧力検出値と同真の圧力補正値とに基いてシリンダ内の筒内絶対圧力を精度良く求めることができるようになる。従って、内燃機関の異常燃焼（例えば、ノッキング、ミスファイア等）の発生を筒内絶対圧力の変化に基いて検出する装置を有する内燃機関に対して本発明による筒内ガス状態取得装置を適用すれば、内燃機関の異常燃焼の発生を精度良く検出できるようになる。
【００２３】
また、上記いずれかの筒内ガス状態取得装置においては、前記所定の時点は、吸気弁が閉弁する圧縮行程開始時点に設定されるとともに、前記所定の区間は、その始期が前記圧縮行程開始時点に、その終期が圧縮行程終了時点になるように設定されることが好適である。ここで、圧縮行程終了時点は、シリンダ内のガスが燃焼を開始して爆発行程に推移する時点である。これによれば、所定の区間の幅を最大とすることができ、筒内絶対圧力値設定値の波形は、より一層多数の検出値（実測値）に基いて求められ得るので、検出点数の少なさに起因する検出精度の低下の影響を最小限に抑制することが可能となる。
【００２４】
【発明の実施の形態】
以下、本発明による内燃機関の筒内ガス状態取得装置を含む燃料噴射量制御装置の実施形態について図面を参照しつつ説明する。図１は、この燃料噴射量制御装置を火花点火式多気筒（例えば、４気筒）内燃機関１０に適用したシステムの概略構成を示している。
【００２５】
内燃機関１０は、シリンダブロック、シリンダブロックロワーケース、及びオイルパン等を含むシリンダブロック部２０と、シリンダブロック部２０の上に固定されるシリンダヘッド部３０と、シリンダブロック部２０にガソリン混合気を供給するための吸気系統４０と、シリンダブロック部２０からの排ガスを外部に放出するための排気系統５０とを含んでいる。
【００２６】
シリンダブロック部２０は、シリンダ２１、ピストン２２、コンロッド２３、及びクランク軸２４を含んでいる。ピストン２２はシリンダ２１内を往復動し、ピストン２２の往復動がコンロッド２３を介してクランク軸２４に伝達され、これにより同クランク軸２４が回転するようになっている。シリンダ２１とピストン２２のヘッドは、シリンダヘッド部３０とともに燃焼室２５を形成している。
【００２７】
シリンダヘッド部３０は、燃焼室２５に連通した吸気ポート３１、吸気ポート３１を開閉する吸気弁３２、吸気弁３２を駆動するインテークカムシャフトを含むとともに同インテークカムシャフトの位相角及び同吸気弁３２のバルブリフト量（最大バルブリフト量）を連続的に変更し得る吸気弁制御装置３３、吸気弁制御装置３３のアクチュエータ３３ａ、燃焼室２５に連通した排気ポート３４、排気ポート３４を開閉する排気弁３５、排気弁３５を駆動するエキゾーストカムシャフト３６、点火プラグ３７、点火プラグ３７に与える高電圧を発生するイグニッションコイルを含むイグナイタ３８、及び燃料を吸気ポート３１内に噴射するインジェクタ（燃料噴射手段）３９を備えている。
【００２８】
吸気系統４０は、吸気ポート３１に連通し同吸気ポート３１とともに吸気通路を形成するインテークマニホールドを含む吸気管４１、吸気管４１の端部に設けられたエアフィルタ４２、吸気管４１内にあって吸気通路の開口断面積を可変とするスロットルバルブ４３、及びスワールコントロールバルブ（以下、「ＳＣＶ」と称呼する。）４４を備えている。スロットルバルブ４３は、ＤＣモータからなるスロットルバルブアクチュエータ４３ａにより吸気管４１内で回転駆動されるようになっている。ＳＣＶ４４は、前記スロットルバルブ４３よりも下流で前記インジェクタ３９よりも上流の位置にて前記吸気管４１に対し回動可能に支持されるとともに、ＤＣモータからなるＳＣＶアクチュエータ４４ａにより回転駆動されるようになっている。
【００２９】
排気系統５０は、排気ポート３４に連通したエキゾーストマニホールド５１、エキゾーストマニホールド５１に接続されたエキゾーストパイプ５２、及びエキゾーストパイプ５２に介装された触媒コンバータ（三元触媒装置）５３を備えている。
【００３０】
一方、このシステムは、熱線式エアフローメータ６１、吸気温センサ６２、大気圧センサ（スロットルバルブ上流圧力センサ）６３、スロットルポジションセンサ６４、ＳＣＶ開度センサ６５、カムポジションセンサ６６、吸気弁リフト量センサ６７、クランクポジションセンサ６８、水温センサ６９、Ｏ２センサ７０、アクセル開度センサ７１、及び筒内相対圧力センサ７２を備えている。
【００３１】
エアフローメータ６１は、吸気管４１内を流れる吸入空気の質量流量に応じた電圧Ｖｇを出力するようになっている。吸気温センサ６２は、エアフローメータ６１内に備えられていて、吸入空気の温度を検出し、吸気温度Ｔａを表す信号を出力するようになっている。大気圧センサ６３は、スロットルバルブ４３の上流の圧力（即ち、大気圧）を検出し、スロットルバルブ上流圧力Ｐａを表す信号を出力するようになっている。スロットルポジションセンサ６４は、スロットルバルブ４３の開度（スロットルバルブ開度）を検出し、スロットルバルブ開度ＴＡを表す信号を出力するようになっている。ＳＣＶ開度センサ６５は、ＳＣＶ４４の開度を検出し、ＳＣＶ開度θｉｖを表す信号を出力するようになっている。
【００３２】
カムポジションセンサ６６は、インテークカムシャフトが９０°回転する毎に（即ち、クランク軸２４が１８０°回転する毎に）一つのパルスを有する信号（Ｇ２信号）を発生するようになっている。吸気弁リフト量センサ６７は、吸気弁３１のリフト量を検出し、吸気弁が全閉のとき「０」の値をとる吸気弁リフト量Ｌを表す信号を出力するようになっている。クランクポジションセンサ６８は、クランク軸２４が１０°回転する毎に幅狭のパルスを有するとともに同クランク軸２４が３６０°回転する毎に幅広のパルスを有する信号を出力するようになっている。この信号は、クランク角θｃａ、及びエンジン回転速度Ｎｅを表す。
【００３３】
水温センサ６９は、内燃機関１０の冷却水の温度を検出し、冷却水温ＴＨＷを表す信号を出力するようになっている。Ｏ２センサ７０は、触媒コンバータ５３に流入する排ガス中の酸素濃度に応じた信号を出力するようになっている。アクセル開度センサ７１は、運転者によって操作されるアクセルペダルの操作量Ａｃｃｐを表す信号を出力するようになっている。筒内相対圧力センサ７２は、シリンダ２１内（燃焼室２５内）の圧力を所定の基準圧力からの相対圧力として検出し、筒内相対圧力検出値Ｐｍｅａｓを表す信号を出力するようになっている。
【００３４】
電気制御装置８０は、互いにバスで接続されたＣＰＵ８１、ＣＰＵ８１が実行するプログラム、テーブル（マップ）、定数等を予め記憶したＲＯＭ８２、ＣＰＵ８１が必要に応じてデータを一時的に格納するＲＡＭ８３、電源が投入された状態でデータを格納するとともに同格納したデータを電源が遮断されている間も保持するバックアップＲＡＭ８４、及びＡＤコンバータを含むインターフェース８５等からなるマイクロコンピュータである。インターフェース８５は、前記センサ６１〜７２と接続され、ＣＰＵ８１にセンサ６１〜７２からの信号を供給するとともに、同ＣＰＵ８１の指示に応じて吸気弁制御装置３３のアクチュエータ３３ａ、イグナイタ３８、インジェクタ３９、スロットルバルブアクチュエータ４３ａ、及びＳＣＶアクチュエータ４４ａに駆動信号を送出するようになっている。
【００３５】
次に、上記のように構成された燃料噴射量制御装置によるシミュレーションモデルを用いた燃料噴射量の決定方法（筒内吸入空気量Ｍｃの推定方法）について説明する。以下に述べる処理は、ＣＰＵ８１がプログラムを実行することによりなされる。
【００３６】
（燃料噴射量ｆｉの決定方法・筒内吸入空気量Ｍｃの推定方法）
この燃料噴射量制御装置（吸入空気量推定装置）は、吸気行程にある気筒の吸気弁３２が閉じる前に同気筒に対して燃料を噴射しなければならないので、吸気弁３２が閉じた時点で（即ち、吸気弁閉時に）同気筒内に吸入されているであろう吸入空気量（筒内（燃焼室内）吸入空気量）を予測する必要がある。一方、吸気弁閉時の吸気管圧力ＰＭＦＷＤは、燃焼室２５に吸入されている筒内吸入空気量Ｍｃと比例関係にある。従って、吸気管圧力ＰＭＦＷＤを予測することができれば、実際の筒内吸入空気量Ｍｃを推定することができる。
【００３７】
そこで、本燃料噴射量制御装置は、吸気弁閉時の吸気管圧力ＰＭＦＷＤを予測・推定し、推定した吸気管圧力ＰＭＦＷＤを一気筒の排気量と空気密度の積で除した値に所定の係数を乗算することにより筒内吸入空気量Ｍｃを求め、下記数１に基づいて燃料噴射量ｆｉを決定する。数１において、Ｋは設定空燃比に応じて変化する係数である。
【００３８】
【数１】
ｆｉ　＝Ｋ・Ｍｃ
【００３９】
以下、吸気弁閉時の吸気管圧力ＰＭＦＷＤの推定方法について、同推定に使用するモデルとともに説明する。図２に示したように、吸気弁閉時の吸気管圧力ＰＭＦＷＤは電子制御スロットルモデルＭ１、スロットルモデルＭ２、吸気弁モデルＭ３、及びインテークマニホールドモデルＭ４により推定される。
【００４０】
（１）電子制御スロットルモデルＭ１
電子制御スロットルモデルＭ１は、現時点までのアクセルペダル操作量Ａｃｃｐに基づいて吸気弁閉時のスロットルバルブ開度ＴＡＳを推定するモデルである。本実施形態においては、スロットルバルブ電子制御ロジックＡ１にて、アクセル開度センサ７１により検出されたアクセルペダル操作量Ａｃｃｐと、図３に示したアクセルペダル操作量Ａｃｃｐと目標スロットルバルブ開度θｒとの関係を規定するテーブルとに基づいて暫定的な目標スロットルバルブ開度θｒ１が求められ、この暫定的な目標スロットルバルブ開度θｒ１を所定時間Ｔ（例えば、６４ｍｓｅｃ）だけ遅延させた値が最終的な目標スロットルバルブ開度θｒとして決定される。そして、スロットルバルブ電子制御ロジックＡ１（電気制御装置８０）は、実際のスロットルバルブ開度ＴＡが目標スロットルバルブ開度θｒとなるようにスロットルバルブアクチュエータ４３ａに対して駆動信号を送出する。
【００４１】
このように、目標スロットルバルブ開度θｒは、現時点から所定時間Ｔだけ前の時点におけるアクセルペダル操作量Ａｃｃｐに応じて決定されるから、現時点から吸気弁閉時までの時間をｔとすると、吸気弁閉時の目標スロットルバルブ開度θｒは、現時点から時間（Ｔ−ｔ）前における暫定的な目標スロットルバルブ開度θｒ１と等しい。また、目標スロットルバルブ開度θｒは、スロットルバルブアクチュエータ４３ａの作動遅れ時間を無視すれば、スロットルバルブ開度ＴＡＳと等しい。このような考えに基づき、電子制御スロットルモデルＭ１は、検出されるエンジン回転速度Ｎｅと、内燃機関１０の運転状態に応じて別途定められる吸気弁の開閉タイミング（進角量）ＶＴ（上記信号Ｎｅと上記Ｇ２信号とにより求めた実際の開閉タイミングＶＴでも良い。）と等に基づいて現時点から吸気弁閉時までの時間ｔを求め、同時間ｔと、現時点から所定時間Ｔだけ前の時点から現時点までのアクセルペダル操作量Ａｃｃｐ（又は、暫定的な目標スロットルバルブ開度θｒ１）の変化の経緯とに基づいて吸気弁閉時のスロットルバルブ開度ＴＡＳを推定する。なお、スロットルバルブアクチュエータ４３ａの作動遅れ時間を考慮に加えて、吸気弁閉時のスロットルバルブ開度ＴＡＳを推定してもよい。
【００４２】
（２）スロットルモデルＭ２
スロットルモデルＭ２は、スロットルバルブ４３を通過する空気量（スロットル通過空気量）ｍｔを、エネルギー保存則、運動量保存則、質量保存則、及び状態方程式に基づいて得られた下記数２及び下記数３に基づいて推定するモデルである。下記数２及び下記数３において、μは流量係数、Ａｔはスロットル開口面積、νはスロットルバルブ４３を通過する空気の流速、Ｐａはスロットルバルブ上流圧力、Ｐｍは吸気管圧力、Ｔａは吸気温度、ρｍは吸気密度、Ｒは気体定数、及びκは比熱比（以下、本吸入空気量推定装置において、κを一定値として扱う。）である。
【００４３】
【数２】
ｍｔ＝μ・Ａｔ・ν・ρｍ＝μ・Ａｔ・｛Ｐａ／（Ｒ・Ｔａ）１／２｝・Φ（Ｐｍ／Ｐａ）
【００４４】
【数３】
2004108348.tif 000003
【００４５】
ここで、上記数２は、ｋ１を所定の係数（＝μ・Ａｔ・｛Ｐａ／（Ｒ・Ｔａ）１／２｝）、ｍｔｓを吸気弁閉時のスロットル通過空気量とするとき下記数４に書き換えられる。また、数４において、内燃機関１０が定常状態にある場合（スロットルバルブ開度が一定である場合）のスロットル通過空気量をｍｔｓＴＡ、及び吸気管圧力をＰｍＴＡとすると、下記数５が得られるので、数４及び数５から係数ｋ１を消去して下記数６を得ることができる。
【００４６】
【数４】
ｍｔｓ＝ｋ１・Φ（Ｐｍ／Ｐａ）
【００４７】
【数５】
ｍｔｓＴＡ＝ｋ１・Φ（ＰｍＴＡ／Ｐａ）
【００４８】
【数６】
ｍｔｓ＝｛ｍｔｓＴＡ／Φ（ＰｍＴＡ／Ｐａ）｝・Φ（Ｐｍ／Ｐａ）
【００４９】
上記数６の右辺における値｛ｍｔｓＴＡ／Φ（ＰｍＴＡ／Ｐａ）｝は、スロットルバルブ開度ＴＡが一定であるときの吸入空気流量（スロットル通過空気量）に関する値であり、スロットルバルブ開度ＴＡ、エンジン回転速度Ｎｅ、吸気弁の開閉タイミングＶＴ、及びスロットルバルブ上流圧力Ｐａが決定されると、実質的に一意に定まる値である。スロットルモデルＭ２は、スロットルバルブ開度ＴＡ、エンジン回転速度Ｎｅ、吸気弁の開閉タイミングＶＴ、及びスロットルバルブ上流圧力Ｐａと、値｛ｍｔｓＴＡ／Φ（ＰｍＴＡ／Ｐａ）｝との関係を規定したテーブルをＲＯＭ８２内に記憶していて、このテーブルと吸気弁閉時の推定スロットルバルブ開度ＴＡＳ、実際のエンジン回転速度Ｎｅ、実際の吸気弁の開閉タイミングＶＴ、及び実際のスロットルバルブ上流圧力Ｐａとに基づいて値｛ｍｔｓＴＡ／Φ（ＰｍＴＡ／Ｐａ）｝を求める。
【００５０】
また、数６の右辺における値Φ（Ｐｍ／Ｐａ）は、上記数３から理解されるように、比熱比κが一定であるとき、吸気管圧力Ｐｍとスロットルバルブ上流圧力Ｐａにより決定される値である。スロットルモデルＭ２は、吸気管圧力Ｐｍ及びスロットルバルブ上流圧力Ｐａと、値Φ（Ｐｍ／Ｐａ）との関係を規定したテーブルをＲＯＭ８２内に記憶していて、このテーブルと、後述するインテークマニホールドモデルＭ４が現時点で既に演算している最新の吸気管圧力Ｐｍ、及び実際のスロットルバルブ上流圧力Ｐａとに基づいて値Φ（Ｐｍ／Ｐａ）を求める。以上により、吸気弁閉時のスロットル通過空気量ｍｔｓが求められる。
【００５１】
（３）吸気弁モデルＭ３
吸気弁モデルＭ３は、吸気管圧力Ｐｍ、吸気管内温度Ｔｍ、及び吸気温度ＴＨＡ等から筒内吸入空気流量ｍｃを推定するモデルである。吸気弁閉時の気筒内圧力は吸気弁３２の上流の圧力、即ち吸気弁閉時の吸気管圧力Ｐｍとみなすことができるので、筒内吸入空気流量ｍｃは吸気弁閉時の吸気管圧力Ｐｍに比例する。そこで、吸気弁モデルＭ３は筒内吸入空気流量ｍｃを、経験則に基づく下記数７にしたがって求める。
【００５２】
【数７】
ｍｃ＝（Ｔａ／Ｔｍ）・（ｃ・Ｐｍ−ｄ）
【００５３】
数７において、値ｃは比例係数、値ｄは筒内に残存していた既燃ガス量である。吸気弁モデルＭ３は、エンジン回転速度Ｎｅ、及び吸気弁の開閉タイミングＶＴと、比例係数ｃ、及び既燃ガス量ｄとの関係をそれぞれ規定するテーブルをＲＯＭ８２内に格納していて、実際のエンジン回転速度Ｎｅと、実際の吸気弁の開閉タイミングＶＴと前記格納しているテーブルとから比例係数ｃ、及び既燃ガス量ｄを求める。また、吸気弁モデルＭ３は、演算時点において、後述するインテークマニホールドモデルＭ４により既に推定されている直前（最新）の吸気弁閉時の吸気管圧力Ｐｍと直前の吸気管内空気温度Ｔｍとを上記数７に適用し、吸気弁閉時の筒内吸入空気流量ｍｃを推定する。
【００５４】
（４）インテークマニホールドモデルＭ４
インテークマニホールドモデルＭ４は、質量保存則とエネルギー保存則とにそれぞれ基づいた下記数８及び下記数９にしたがって、吸気弁閉時の吸気管圧力Ｐｍと、吸気弁閉時の吸気管内温度Ｔｍとを求める。なお、Ｖは吸気管の容積、Ｒは気体定数、ｍｔはスロットル通過空気量、Ｔａはスロットルバルブ通過空気温度（即ち、吸気温度Ｔａ）である。
【００５５】
【数８】
ｄＰｍ／ｄｔ＝κ・（Ｒ／Ｖ）・（ｍｔ・Ｔａ−ｍｃ・Ｔｍ）
【００５６】
【数９】
ｄ（Ｐｍ／Ｔｍ）／ｄｔ＝（Ｒ／Ｖ）・（ｍｔ−ｍｃ）
【００５７】
図２に示したように、インテークマニホールドモデルＭ４は、スロットルモデルＭ２により推定されたスロットル通過空気量ｍｔｓを上記数８，数９におけるスロットル通過空気量ｍｔとして使用し、吸気弁モデルＭ３により推定された吸気弁閉時の筒内吸入空気流量ｍｃを上記数８，数９の筒内吸入空気流量ｍｃとして使用する。このインテークマニホールドモデルＭ４により推定された吸気管圧力Ｐｍが、前記吸気弁閉時の推定吸気管圧力ＰＭＦＷＤとなる。
【００５８】
以上のようにして、本燃料噴射量制御装置は、電子制御スロットルモデルＭ１、スロットルモデルＭ２、吸気弁モデルＭ３、インテークマニホールドモデルＭ４の各シミュレーションモデルを用いて吸気弁閉時の吸気管圧力ＰＭＦＷＤを予測・推定し、推定した吸気管圧力ＰＭＦＷＤに対して所定の計算を行うことにより筒内吸入空気量Ｍｃを求め、上記数１に基づいて燃料噴射量ｆｉを決定する。
【００５９】
（圧縮行程開始時点におけるシリンダ内のガス状態の取得）
本燃料噴射量制御装置は、通常、上述したように各種モデルを使用して、各気筒が吸気行程を迎える度に同吸気行程を迎える特定気筒における筒内吸入空気量Ｍｃを求める一方で、本燃料噴射量制御装置に含まれる本発明による筒内ガス状態取得装置は、所定の条件が成立する度に、以下のようにして、同特定気筒における圧縮行程開始時点におけるシリンダ内のガス状態を取得して、同シリンダ内に吸入されている筒内吸入ガス量（質量）、圧縮行程開始時点における同シリンダ内の筒内絶対圧力等を求める。
【００６０】
（１）筒内相対圧力検出値の取得
先ず、本装置は、クランクポジションセンサ６８により得られる信号とカムポジションセンサ６６により得られる上記Ｇ２信号とにより求められる吸気弁３２の開閉タイミング（進角量）ＶＴに基き、上記特定気筒の吸気弁３２が閉弁する（開状態から閉状態に移行する）圧縮行程開始時点（所定の時点）におけるクランク角θｃａを求めて、その値をクランク角θ（１）として設定する。
【００６１】
そして、本装置は、図４に示すように、クランクポジションセンサ６８により得られるクランク角θｃａがクランク角θ（１）になると、その時点における筒内相対圧力センサ７２により得られる筒内相対圧力検出値Ｐｍｅａｓを筒内相対圧力検出値Ｐｍｅａｓ（１）としてＲＡＭ８３に格納する。以降、本装置は、クランク角θｃａが所定の一定間隔（計算間隔）Δθだけ進む度にクランク角θｃａの値を順にクランク角θ（Ｈ）　（Ｈ＝２，３，・・・）として設定していくとともに、クランク角θｃａがクランク角θ（Ｈ）　（Ｈ＝２，３，・・・）になる度にその時点における上記筒内相対圧力検出値Ｐｍｅａｓを筒内相対圧力検出値Ｐｍｅａｓ（Ｈ）　（Ｈ＝２，３，・・・）として順にＲＡＭ８３に格納していく。
【００６２】
本装置は、このような処理を上記特定気筒の点火プラグ３７が点火されて同特定気筒が爆発行程に移行する圧縮行程終了時点まで継続して行い、同圧縮行程終了時点における上記変数Ｈの値をサンプル数Ｎとして設定する。このようにして、本装置は、所定の時点である圧縮行程開始時点から圧縮行程終了時点までの所定の区間内における複数の（Ｎ個の）筒内相対圧力検出値Ｐｍｅａｓ（Ｈ）　（Ｈ＝１，２，・・・，Ｎ−１，Ｎ）を取得して順次ＲＡＭ８３に記憶する。
【００６３】
（２）筒内吸入ガス量、及び筒内ガス平均温度の仮設定
次に、本装置は、上記複数の筒内相対圧力検出値Ｐｍｅａｓ（Ｈ）　（Ｈ＝１，２，・・・，Ｎ−１，Ｎ）を取得した上記圧縮行程時（以下、この圧縮行程を「特定圧縮行程」と呼ぶ。）に上記特定気筒のシリンダ内に吸入されている筒内吸入ガス（混合気）量（質量）Ｍｃｍ、及び特定圧縮行程開始時点における同シリンダ内の筒内ガス平均温度Ｔｃｍを、それぞれ下記数１０及び下記数１１により仮設定する。
【００６４】
【数１０】
Ｍｃｍ＝ｆｉ＋Ｍｃ
【００６５】
【数１１】
Ｔｃｍ＝Ｔｍ
【００６６】
上記数１０において、ｆｉは上記特定気筒の上記特定圧縮行程開始時点の直前に同特定気筒内に実際に噴射された燃料噴射量であり、Ｍｃは上記各モデルを用いて推定されている同特定気筒の上記特定圧縮行程時における筒内吸入空気量である。上記数１１において、Ｔｍは上記インテークマニホールドモデルＭ４により推定されている上記特定気筒の吸気弁閉時（上記特定圧縮行程開始時点）の吸気管内温度Ｔｍである。このように、上記数１０及び上記数１１を利用して、筒内吸入ガス量Ｍｃｍ、及び筒内ガス平均温度Ｔｃｍをそれぞれ仮設定する手段が、仮設定手段を構成する。
【００６７】
（３）仮の圧力補正値の算出、及び仮の筒内絶対圧力設定値の設定
次に、本装置は、上記のように仮設定した筒内吸入ガス量Ｍｃｍの値及び筒内ガス平均温度Ｔｃｍの値を利用して、状態方程式に基く下記数１２により上記特定気筒の上記特定圧縮行程開始時点における仮の筒内絶対圧力Ｐ０（上記特定圧縮行程開始時点における上記特定気筒のシリンダ内のガスの状態の一つ）を求める（仮設定する）。
【００６８】
【数１２】
Ｐ０＝（Ｍｃｍ・Ｒｇ・Ｔｃｍ）／Ｖｃｍ（１）
【００６９】
上記数１２において、Ｒｇはガス定数であり（上述した気体定数Ｒとは異なる）、Ｖｃｍ（１）は、上記特定圧縮行程開始時点、即ちクランク角θｃａが上記クランク角θ（１）となる時点における既知のピストン２２の位置を利用して求まるシリンダ容積である。
【００７０】
次に、本装置は、上記のようにして求めた仮の筒内絶対圧力Ｐ０の値と、上記特定圧縮行程開始時点における上記筒内相対圧力検出値Ｐｍｅａｓ（１）とから、下記数１３により図４に示す仮の圧力補正値ΔＰを求める。
【００７１】
【数１３】
ΔＰ＝Ｐ０−Ｐｍｅａｓ（１）
【００７２】
そして、本装置は、上記のようにして求めた仮の圧力補正値ΔＰと、上記Ｎ個の筒内相対圧力検出値Ｐｍｅａｓ（Ｈ）　（Ｈ＝１，２，・・・，Ｎ−１，Ｎ）の各々とに基いて、下記数１４により図４に示す上記所定の区間内におけるＮ個の仮の筒内絶対圧力設定値Ｐｃ（Ｈ）（Ｈ＝１，２，・・・，Ｎ−１，Ｎ）を求める。このように、上記数１３を利用して仮の圧力補正値ΔＰを求めるとともに、下記数１４を利用して仮の筒内絶対圧力設定値Ｐｃ（Ｈ）　（Ｈ＝１，２，・・・，Ｎ−１，Ｎ）を求める手段が、第１計算手段を構成する。
【００７３】
【数１４】
Ｐｃ（Ｈ）＝Ｐｍｅａｓ（Ｈ）＋ΔＰ　（Ｈ＝１，２，・・・，Ｎ−１，Ｎ）
【００７４】
（４）仮の筒内絶対圧力計算値の計算
次に、本装置は、上記のようにして仮設定した上記特定圧縮行程開始時点における上記特定気筒のシリンダ内のガスの状態のうちの仮の筒内絶対圧力Ｐ０の値、及び筒内吸入ガス量Ｍｃｍを初期条件として、エネルギー保存則に基づく下記数１５〜下記数１７、及び状態方程式に基いた下記数１８により記述されたシリンダモデルにより、図４に示す上記所定の区間内におけるＮ個の仮の筒内絶対圧力計算値Ｐｃｍ（Ｈ）　（Ｈ＝１，２，・・・，Ｎ−１，Ｎ）を求める。
【００７５】
【数１５】
2004108348.tif 000004
【００７６】
【数１６】
2004108348.tif 000005
【００７７】
【数１７】
2004108348.tif 000006
【００７８】
【数１８】
2004108348.tif 000007
【００７９】
上記数１５はエネルギー保存則に基く基本式であって、上記特定気筒を模式的に表した図５に示すように、同数１５において、Ｐｃｍは仮の筒内絶対圧力計算値、Ｖｃｍはシリンダ容積、Ｑはシリンダ２１外部（シリンダ壁面、吸気ポート等）からシリンダ２１内に伝達される単位時間当たりの熱伝達量、κは比熱比である。上記数１６は、上記数１５における熱伝達量Ｑを表す基本式であって、同数１６において、Ａは筒内表面積、ｈｗは熱伝達率、Ｔｃｍは筒内ガス平均温度、Ｔｗは筒内壁面温度である。
【００８０】
数１７は所謂ウォッシーニ（Ｗｏｓｃｈｎｉ）モデルを使用して上記数１６における熱伝達率ｈｗを表した式であって、同数１７において、ｄｉはシリンダ内径、Ｃ１は定数（本例では、２．２８）、Ｃｗは平均ピストン速度である。上記数１８は上記特定圧縮行程において特定気筒内に密閉されているガス（混合気）に適用される状態方程式に基く式であって、同数１８において、Ｒｇは上述した気体定数Ｒとは異なるガス定数、Ｍｃｍは筒内吸入ガス量（質量）である。
【００８１】
ここで、シリンダ容積Ｖｃｍ、シリンダ容積の時間微分値ｄＶｃｍ／ｄｔ、筒内表面積Ａはクランクポジションセンサ６８により得られるクランク角θｃａに基いて求めることができ、平均ピストン速度Ｃｗはエンジン回転速度Ｎｅに基いて求めることができ、筒内壁面温度Ｔｗは水温センサ６９により得られる冷却水温ＴＨＷにて代用することができる。また、比熱比κは筒内ガス平均温度Ｔｃｍの関数として表すことができる。
【００８２】
従って、上記数１５〜上記数１８により記述されたシリンダモデルは、筒内吸入ガス量Ｍｃｍ、及び時刻ｔにおける仮の筒内絶対圧力計算値Ｐｃｍの値が求められていると、上記数１８により時刻ｔにおける筒内ガス平均温度Ｔｃｍを求めることができ、その結果、上記数１７により時刻ｔにおける熱伝達率ｈｗを求めることができる。そして、上記数１６により時刻ｔにおける熱伝達量Ｑを求めることができ、その結果、理論上、上記数１５により時刻（ｔ＋ｄｔ）における仮の筒内絶対圧力計算値Ｐｃｍを求めることができる。よって、筒内吸入ガス量Ｍｃｍ、及び仮の筒内絶対圧力Ｐ０の値が初期条件（上記特定圧縮行程開始時点における値）として付与されると、上記シリンダモデルは、理論上、上記所定の区間内における仮の筒内絶対圧力計算値Ｐｃｍを順次求めることができる。
【００８３】
実際には、本装置は、以下のようにして、上記所定の区間内におけるＮ個の仮の筒内絶対圧力計算値Ｐｃｍ（Ｈ）　（Ｈ＝１，２，・・・，Ｎ−１，Ｎ）を求める。先ず、上記数１５において、ｄ／ｄｔを（ｄ／ｄθ）・（ｄθ／ｄｔ）＝（ｄ／ｄθ）・ω　（θはクランク角，ωはクランクシャフトの角速度（即ち、エンジン回転速度Ｎｅ））に書き改めると、下記数１９が得られる。
【００８４】
【数１９】
2004108348.tif 000008
【００８５】
上記数１９の両辺に（κ−１）／（ω・Ｖｃｍ）を乗算して整理すると下記数２０が得られる。
【００８６】
【数２０】
2004108348.tif 000009
【００８７】
ここで、上記数２０を、クランク角θに関して計算間隔Δθをもって離散化するため、数２０の左辺のｄＰｃｍ／ｄθ，ｄκ／ｄθを、それぞれ下記数２１，下記数２２を用いて書き改め、さらに、数２０におけるＰｃｍをＰｃｍ（θ）、κをκ（θ）、ＶｃｍをＶｃｍ（θ）、ｄＶｃｍ／ｄθをｄＶｃｍ（θ）／ｄθ、ωをω（θ）、ＱをＱ（θ）にそれぞれ書き改めると、下記数２３が得られる。
【００８８】
【数２１】
2004108348.tif 000010
【００８９】
【数２２】
2004108348.tif 000011
【００９０】
【数２３】
2004108348.tif 000012
【００９１】
また、上記数１６において、ＡをＡ（θ）、ｈｗをｈｗ（θ）、ＴｃｍをＴｃｍ（θ）にそれぞれ書き改めると、下記数２４が得られる。
【００９２】
【数２４】
2004108348.tif 000013
【００９３】
また、上記数１７及び上記数１８を書き改めると、それぞれ、下記数２５及び下記数２６が得られる。
【００９４】
【数２５】
2004108348.tif 000014
【００９５】
【数２６】
2004108348.tif 000015
【００９６】
そして、本装置は、初期条件として、筒内吸入ガス量Ｍｃｍを上記のように仮設定された筒内吸入ガス量Ｍｃｍに設定し、上記仮の筒内絶対圧力Ｐ０の値を、クランク角θｃａが上記クランク角θ（１）となる上記特定圧縮行程開始時点における仮の筒内絶対圧力計算値Ｐｃｍ（１）に設定することにより、上記数２６によりクランク角θｃａがクランク角θ（１）となるときにおける筒内ガス平均温度Ｔｃｍ（１）を求めることができ、その結果、上記数２５によりクランク角θｃａがクランク角θ（１）となるときにおける熱伝達率ｈｗ（１）を求めることができる。また、上記筒内ガス平均温度Ｔｃｍ（１）によりクランク角θｃａがクランク角θ（１）となるときにおける比熱比κ（１）を求めることができる。なお、κ（０）はκ（１）と同一の値として設定される。
【００９７】
また、本装置は、上記数２４によりクランク角θｃａがクランク角θ（１）となるときにおける熱伝達量Ｑ（１）を求めることができ、その結果、上記数２３によりクランク角θｃａがクランク角θ（２）（クランク角θ（１）＋Δθ）となるときにおける仮の筒内絶対圧力計算値Ｐｃｍ（２）を求めることができる。換言すれば、本装置は、仮の筒内絶対圧力計算値Ｐｃｍ（１）が求まれば、上記数２３〜上記数２６を使用することにより仮の筒内絶対圧力計算値Ｐｃｍ（２）を求めることができる。
【００９８】
従って、上記数２３〜上記数２６を繰り返し使用することにより、本装置は、図４に示すように、上記所定の区間内におけるＮ個の仮の筒内絶対圧力計算値Ｐｃｍ（Ｈ）　（Ｈ＝１，２，・・・，Ｎ−１，Ｎ）を順次求めることができる。このように、上記数２３〜上記数２６を利用して仮の筒内絶対圧力計算値Ｐｃｍ（Ｈ）　（Ｈ＝１，２，・・・，Ｎ−１，Ｎ）を求める手段が、第２計算手段を構成する。
【００９９】
（５）偏差指標値の計算
次に、本装置は、上記所定の区間におけるＮ個の仮の筒内絶対圧力設定値Ｐｃ（Ｈ）とＮ個の仮の筒内絶対圧力計算値Ｐｃｍ（Ｈ）　（Ｈ＝１，２，・・・，Ｎ−１，Ｎ）との偏差の程度（図４において多数のドットで示された領域の面積に応じて変化する値）を示す偏差指標値として、下記数２７により表される値Ｅを採用する。このように、下記数２７を利用して偏差指標値Ｅを求める手段が、偏差指標値計算手段を構成する。
【０１００】
【数２７】
2004108348.tif 000016
【０１０１】
（６）筒内吸入ガス量Ｍｃｍ、及び筒内ガス平均温度Ｔｃｍの同定
上述したように偏差指標値Ｅを求めた後、本装置は、上記偏差指標値Ｅが最小になるように、上記した（２）〜（５）の一連の各処理を繰り返し実行することにより、上記数１０及び上記数１１により仮設定した筒内吸入ガス量Ｍｃｍ及び上記特定圧縮行程開始時点での筒内ガス平均温度Ｔｃｍを同定し、同定された筒内吸入ガス量Ｍｃｍの値及び筒内ガス平均温度Ｔｃｍの値を、それぞれ真の値である確定筒内吸入ガス量Ｍｃｆｉｎ及び確定筒内ガス平均温度Ｔｃｆｉｎとして設定する。かかる同定方法の詳細は後述する。
【０１０２】
そして、本装置は、上記数１２において、Ｍｃｍの代わりに確定筒内吸入ガス量Ｍｃｆｉｎの値を、Ｔｃｍの代わりに確定筒内ガス平均温度Ｔｃｆｉｎをそれぞれ適用することにより、上記特定圧縮行程開始時点での筒内絶対圧力の真の値である確定筒内絶対圧力Ｐｃｆｉｎを求め、上記特定圧縮行程開始時点での特定気筒のシリンダ内に吸入されているガスの真の状態を完全に取得する。このようにして、上記特定圧縮行程開始時点での特定気筒のシリンダ内に吸入されているガスの真の状態を取得する手段が、筒内ガス状態取得手段を構成する。
【０１０３】
また、これにより、本装置は、上記確定筒内吸入ガス量Ｍｃｆｉｎから実際の燃料噴射量ｆｉを減じることにより上記特定気筒の上記特定圧縮行程時における真の筒内吸入空気量を求め、この真の筒内吸入空気量の値と、上述した各モデルを用いて推定されている同特定気筒の同特定圧縮行程時における筒内吸入空気量Ｍｃの値との比較により、先に説明した数７における比例係数ｃ及び既燃ガス量ｄをそれぞれ求めるためにＲＯＭ８２に格納されている各テーブルを補正する。
【０１０４】
さらには、本装置は、下記数２８により真の圧力補正値である確定圧力補正値ΔＰｆｉｎを算出する。
【０１０５】
【数２８】
ΔＰｆｉｎ＝Ｐｃｆｉｎ　−　Ｐｍｅａｓ（１）
【０１０６】
従って、本装置は、上記数２８により確定圧力補正値ΔＰｆｉｎを算出した時点以降、筒内相対圧力センサ７２の出力値である筒内相対圧力検出値Ｐｍｅａｓに確定圧力補正値ΔＰｃｆｉｎを加算することにより、シリンダ２１内の筒内絶対圧力を精度良く求めることができるようになり、内燃機関１０の異常燃焼（例えば、ノッキング、ミスファイア等）の発生を精度良く検出できるようになる。このように、上記数２８を利用して確定圧力補正値ΔＰｆｉｎを求める手段が、圧力補正値取得手段を構成する。以上のようにして、上記特定気筒の上記特定圧縮行程開始時点におけるシリンダ内のガス状態が取得される。
【０１０７】
（実際の作動）
次に、以上のように構成された本燃料噴射量制御装置、及び同燃料噴射量制御装置に含まれる本発明による筒内ガス状態取得装置の実際の作動について、電気制御装置８０のＣＰＵ８１が実行するルーチンをフローチャートにより示した図６〜図１２を参照しながら説明する。
【０１０８】
（スロットルバルブ制御）
電気制御装置８０のＣＰＵ８１は、図６にフローチャートにより示したスロットルバルブ開度を制御するためのルーチンを所定時間（１ｍｓｅｃ）の経過毎に実行するようになっている。従って、所定のタイミングとなると、ＣＰＵ８１はステップ６００から処理を開始し、ステップ６０５に進んでアクセルペダル操作量Ａｃｃｐ読み込む。次いで、ＣＰＵ８１はステップ６１０に進み、同ステップ６１０にて図３と同じテーブルを用いることにより上記読み込んだアクセルペダル操作量Ａｃｃｐに基づく暫定的な目標スロットルバルブ開度θｒ１を求める。
【０１０９】
次に、ＣＰＵ８１はステップ６１５に進んで変数Ｉを「６４」に設定し、続くステップ６２０にて記憶値θｒ（Ｉ）にθｒ（Ｉ−１）の値を格納する。現時点では、変数Ｉは「６４」であるから、記憶値θｒ（６４）に記憶値θｒ（６３）の値が格納される。次いで、ＣＰＵ８１はステップ６２５に進み、変数Ｉが「１」と等しくなったか否かを判定する。この場合、変数Ｉの値は「６４」であるから、ＣＰＵ８１はステップ６２５にて「Ｎｏ」と判定してステップ６３０に進み、同ステップ６３０にて変数Ｉの値を「１」だけ減少し、その後上記ステップ６２０に戻る。この結果、ステップ６２０が実行されると、記憶値θｒ（６３）に記憶値θｒ（６２）の値が格納される。このような処理は、変数Ｉの値が「１」となるまで繰り返し実行される。
【０１１０】
その後、ステップ６３０の処理が繰り返されて変数Ｉの値が「１」となると、ＣＰＵ８１はステップ６２５にて「Ｙｅｓ」と判定してステップ６３５に進み、同ステップ６３５にて前記ステップ６１０にて求めた現時点における暫定的な目標スロットルバルブ開度θｒ１を記憶値θｒ（０）に格納する。以上により、現時点からＩｍｓｅｃ前（０ｍｓｅｃ≦Ｉｍｓｅｃ≦６４ｍｓｅｃ，Ｉは整数）の暫定的な目標スロットルバルブ開度θｒ（Ｉ）（Ｉ＝６４，６３，６２，・・・，２，１，０）がＲＡＭ８３内に記憶されることになる。
【０１１１】
次に、ＣＰＵ８１はステップ６４０に進み、同ステップ６４０にて記憶値θｒ（６４）を最終的な目標スロットルバルブ開度θｒとして設定し、続くステップ６４５にて実際のスロットルバルブ開度が目標スロットルバルブ開度θｒと等しくなるように、スロットルバルブアクチュエータ４３ａに対し駆動信号を出力し、その後ステップ６９５にて本ルーチンを一旦終了する。
【０１１２】
以降においても、上記ルーチンの処理は１ｍｓｅｃの経過毎に実行される。この結果、実際のスロットルバルブ開度が、所定時間Ｔ（＝６４ｍｓｅｃ）前のアクセルペダル操作量Ａｃｃｐに基づく目標スロットルバルブ開度θｒと等しくなるように制御される。これにより、上記電子制御スロットルモデルＭ１による吸気弁閉時のスロットルバルブ開度ＴＡＳの推定が可能となる。
【０１１３】
（燃料噴射量ｆｉの計算、噴射指示）
ＣＰＵ８１は、吸気行程を迎える特定気筒のクランク角θｃａが、その気筒の吸気上死点から所定クランク角度だけ前の角度（例えば、ＢＴＤＣ９０°）になると、図７の燃料噴射量ｆｉの計算ルーチンの処理をステップ７００から開始してステップ７０５に進み、図２に示した各モデルに従って別途計算されている吸気弁閉時の吸気管圧力ＰＭＦＷＤに基づく筒内吸入空気量Ｍｃを読み込む。なお、筒内吸入空気量Ｍｃは、所定時間毎に繰り返し実行される図２に示した各モデルに従う図示しないルーチンにより求められている。
【０１１４】
次に、ＣＰＵ８１はステップ７１０に進んで、ステップ７０５にて読み込んだ筒内吸入空気量Ｍｃの値と上記数１の右辺とに基づいて燃料噴射量ｆｉを計算する。そして、ＣＰＵ８１はステップ７１５に進んで、ステップ７１０にて計算された燃料噴射量ｆｉだけ燃料を噴射するように前記特定気筒に対するインジェクタ３９に駆動信号を送出し、ステップ７９５に進んで本ルーチンを一旦終了する。
【０１１５】
（筒内吸入ガス量の計算開始判定）
ＣＰＵ８１は、図８にフローチャートにより示した筒内吸入ガス量を計算するためのルーチンを所定時間の経過毎に実行するようになっている。従って、所定のタイミングとなると、ＣＰＵ８１はステップ８００から処理を開始し、ステップ８０５に進んで、筒内相対圧力取得処理実行中フラグＸＨＡＮの値が「０」であるか否かを判定する。筒内相対圧力取得処理実行中フラグＸＨＡＮは、その値が「１」のとき筒内相対圧力センサ７２の出力に基いて筒内相対圧力検出値Ｐｍｅａｓの取得処理を実行していることを示し、その値が「０」のとき同筒内相対圧力検出値Ｐｍｅａｓの取得処理を実行していないことを示す。
【０１１６】
いま、後述する筒内吸入ガス量等の計算開始条件が成立しておらず、且つ、前記筒内相対圧力検出値Ｐｍｅａｓの取得処理を実行していないとして説明を続けると、筒内相対圧力取得処理実行中フラグＸＨＡＮの値が「０」になっている。従って、ＣＰＵ８１はステップ８０５にて「Ｙｅｓ」と判定してステップ８１０に進み、筒内吸入ガス量等の計算開始条件が成立していて、且つ吸気行程を迎える特定気筒のクランク角θｃａが上記開閉タイミング（進角量）ＶＴに基いて計算される吸気弁閉時の角度より計算間隔Δθ（一定値）だけ前の角度になっているか否かを判定する。
【０１１７】
この筒内吸入ガス量等の計算開始条件は、冷却水温ＴＨＷが所定温度以上であり、図示しない車速センサにより得られた車速が所定の高車速以上であり、スロットル弁開度ＴＡの単位時間あたりの変化量が所定量以下である、機関が定常運転されている場合に成立する。更に、かかる計算開始条件に、前回の筒内吸入ガス量等の計算時点から所定時間以上が経過したこと、前回の筒内吸入ガス量等の計算時点から車両が所定距離以上運転されたこと、前回の筒内吸入ガス量等の計算時点から内燃機関１０が所定時間以上運転されたことの任意の一つ、又は一つ以上を加えても良い。現段階では、上述したように、筒内吸入ガス量等の計算開始条件は成立していないから、ＣＰＵ８１はステップ８１０にて「Ｎｏ」と判定してステップ８９５に進み、本ルーチンを一旦終了する。
【０１１８】
次に、上記筒内吸入ガス量等の計算開始条件が成立したものとして説明を続けると、ＣＰＵ８１はステップ８１０に進んだとき、特定気筒のクランク角θｃａが吸気弁閉時の角度より計算間隔Δθだけ前の角度になっていると同ステップ８１０にて「Ｙｅｓ」と判定してステップ８１５に進み、筒内相対圧力取得処理を開始するため、筒内相対圧力取得処理実行中フラグＸＨＡＮの値を「１」に設定し、続くステップ８２０にて変数Ｈの値を「０」に設定するとともに、続くステップ８２５にて現時点でクランクポジションセンサ６８により得られる前記特定気筒のクランク角θｃａの値を変数θ１に格納した後、ステップ８９５に進んで本ルーチンを一旦終了する。
【０１１９】
以降、ＣＰＵ８１は図８のルーチンをステップ８００から繰り返し実行するが、筒内相対圧力取得処理実行中フラグＸＨＡＮの値が「１」になっていることから、ステップ８０５に進んだとき、同ステップ８０５にて「Ｎｏ」と判定してステップ８３０に進み、前記特定気筒の点火プラグ３７が点火される前であるか否か（現時点が特定気筒の特定圧縮行程終了時点より前であるか否か）を判定するようになる。現時点では、前記特定気筒は吸気弁閉時直前（特定圧縮行程開始直前）の状態にあるから、ＣＰＵ８１はステップ８３０にて「Ｙｅｓ」と判定してステップ８９５に進んで本ルーチンを一旦終了する。これにより、特定気筒の特定圧縮行程が終了するまでステップ８００，８０５，８３０，８９５の処理が繰り返し実行されて、筒内相対圧力取得処理実行中フラグＸＨＡＮの値は「１」に維持される。
【０１２０】
（筒内相対圧力等の取得）
一方、ＣＰＵ８１は図９に示した筒内相対圧力等の取得ルーチンを所定時間の経過毎に繰り返し実行している。従って、所定のタイミングになると、ＣＰＵ８１はステップ９００から処理を開始してステップ９０５に進み、筒内相対圧力取得処理実行中フラグＸＨＡＮの値が「１」であるか否かを判定する。ここで、ＣＰＵ８１は、筒内相対圧力取得処理実行中フラグＸＨＡＮの値が「０」であれば直ちにステップ９９５に進んで本ルーチンを一旦終了するが、現時点では先の図８のステップ８１５の処理により筒内相対圧力取得処理実行中フラグＸＨＡＮの値は「１」になっているので、ステップ９０５にて「Ｙｅｓ」と判定してステップ９１０に進み、現時点でクランクポジションセンサ６８により得られる前記特定気筒のクランク角θｃａの値を変数θ２に格納する。
【０１２１】
次に、ＣＰＵ８１はステップ９１５に進んで、ステップ９１０の処理時点でのクランク角θｃａの値が格納されている変数θ２の値から先の図８のステップ８２５の処理時点でのクランク角θｃａの値が格納されている変数θ１の値を減算した値（角度）が上記計算間隔Δθ以上になっているか否かを判定する。現時点は先の図８のステップ８２５の処理を実行した直後であるから、現時点では変数θ２の値から変数θ１の値を減算した値が上記計算間隔Δθ未満であり、ＣＰＵ８１はステップ９１５にて「Ｎｏ」と判定して直ちにステップ９９５に進んで本ルーチンを一旦終了する。
【０１２２】
以降、ＣＰＵ８１は、変数θ２の値から変数θ１の値を減算した値が上記計算間隔Δθ以上になるまで、上記ステップ９００〜ステップ９１５の処理を繰り返し実行する。そして、ステップ９１０の処理が繰り返されてクランク角θｃａの増大とともに変数θ２の値が増大し、変数θ２の値から変数θ１の値を減算した値が上記計算間隔Δθに到達すると、ＣＰＵ８１はステップ９１５に進んだとき、「Ｙｅｓ」と判定してステップ９２０以降の筒内相対圧力等の取得処理を開始する。なお、この時点は、先の図８のステップ８１０にて「Ｙｅｓ」と判定した時点から計算間隔Δθだけ経過した時点であるので、特定気筒の吸気弁閉時（特定圧縮行程開始時点）に対応している。
【０１２３】
ＣＰＵ８１はステップ９２０に進むと、その時点での変数Ｈの値を「１」だけ増大した値を新たな変数Ｈとして格納する。現時点では、先の図８のステップ８２０の処理により変数Ｈの値は「０」になっているので、この処理により変数Ｈの値は「１」に設定される。次に、ＣＰＵ８１はステップ９２５に進み、現時点でのクランク角θｃａが格納されている変数θ２の値をクランク角θ（Ｈ）に格納する。これにより、上記特定圧縮行程開始時点でのクランク角θｃａの値がクランク角θ（１）に格納される。
【０１２４】
次いで、ＣＰＵ８１はステップ９３０に進み、現時点での筒内相対圧力検出値Ｐｍｅａｓを筒内相対圧力検出値Ｐｍｅａｓ（Ｈ）に格納する。これにより、上記特定圧縮行程開始時点での筒内相対圧力検出値Ｐｍｅａｓが筒内相対圧力検出値Ｐｍｅａｓ（１）に格納される。次に、ＣＰＵ８１はステップ９３５に進み、クランク角θ（Ｈ）の値と、クランク角θｃａの関数ｆ１とに基づいて得られるシリンダ容積の値をシリンダ容積Ｖｃｍ（Ｈ）として設定する。これにより、上記特定圧縮行程開始時点でのシリンダ容積の値がシリンダ容積Ｖｃｍ（１）に格納される。
【０１２５】
次に、ＣＰＵ８１はステップ９４０に進み、クランク角θ（Ｈ）の値と、クランク角θｃａの関数ｆ２とに基づいて得られるシリンダ容積のクランク角についての微分値をシリンダ容積微分値ｄＶｃｍ（Ｈ）／ｄθとして設定する。これにより、上記特定圧縮行程開始時点でのシリンダ容積微分値がシリンダ容積微分値Ｖｃｍ（１）／ｄθに格納される。
【０１２６】
次いで、ＣＰＵ８１はステップ９４５に進み、クランク角θ（Ｈ）の値と、クランク角θｃａの関数ｆ３とに基づいて得られる筒内表面積の値を筒内表面積Ａ（Ｈ）として設定する。これにより、上記特定圧縮行程開始時点での筒内表面積の値が筒内表面積Ａ（１）に格納される。
【０１２７】
次に、ＣＰＵ８１はステップ９５０に進んで、現時点でのエンジン回転速度Ｎｅの値を角速度ω（Ｈ）に格納する。これにより、上記特定圧縮行程開始時点でのエンジン回転速度Ｎｅの値が角速度ω（１）に格納される。そして、ＣＰＵ８１はステップ９５５に進んで、変数θ２の値を変数θ１として格納した後ステップ９９５に進んで本ルーチンを一旦終了する。
【０１２８】
以降、ＣＰＵ８１は、筒内相対圧力取得処理実行中フラグＸＨＡＮの値が「１」である限りにおいて、ステップ９１０，９１５の処理を繰り返し実行する。ここで、先に述べたとおり筒内相対圧力取得処理実行中フラグＸＨＡＮの値は特定圧縮行程が終了するまで「１」に維持されている。従って、現時点から特定圧縮行程が終了するまでの間、ステップ９１０の処理が繰り返されてクランク角θｃａの増大とともに変数θ２の値が増大し、同変数θ２の値から既にステップ９５５にて設定されている変数θ１の最新値を減算した値が上記計算間隔Δθに到達する度ごとに、ステップ９２０以降の筒内相対圧力等の取得処理が実行される。
【０１２９】
この結果、特定圧縮行程開始時点以降、同特定圧縮行程終了時点までの間、クランク角θｃａの値が計算間隔Δθだけ増大する毎に、変数Ｈが「１」ずつ増大するとともに、その各々の時点における各値、即ち、クランク角θ（Ｈ）、筒内相対圧力検出値Ｐｍｅａｓ（Ｈ）、シリンダ容積Ｖｃｍ（Ｈ）、シリンダ容積微分値Ｖｃｍ（Ｈ）／ｄθ、筒内表面積Ａ（Ｈ）、角速度ω（Ｈ）　（Ｈ＝１，２，３，・・・）が順次設定されていく。
【０１３０】
そして、上記特定気筒の点火プラグ３７が点火されて上記特定圧縮行程が終了すると、ＣＰＵ８１は図８のステップ８３０に進んだとき、「Ｎｏ」と判定してステップ８３５に進み、現時点での変数Ｈの値をサンプル数Ｎとして格納する。次いで、ＣＰＵ８１はステップ８４０に進み、現時点にて水温センサ６９により得られる冷却水温ＴＨＷの値を筒内壁面温度Ｔｗとして設定し、続くステップ８４５にて、同ステップ８４５内に記載の式に基き、既に図９のステップ９５０にて算出されているＮ個の角速度ω（Ｈ）　（Ｈ＝１，２，・・・，Ｎ−１，Ｎ）の平均値を平均回転速度Ｎｅｍｅａｎとして設定するとともに、続くステップ８５０にて、平均回転速度Ｎｅｍｅａｎの値と、エンジン回転速度Ｎｅの関数ｆ４とに基いて得られる平均ピストン速度を平均ピストン速度Ｃｗとして格納する。
【０１３１】
次に、ＣＰＵ８１はステップ８５５に進んで、上記特定気筒の上記特定圧縮行程開始時点の直前に同特定気筒内に実際に噴射された燃料噴射量ｆｉの値に、図２に示した各モデルに従う図示しないルーチンにより推定されている同特定気筒の同特定圧縮行程時における筒内吸入空気量Ｍｃの値を加えた値を筒内吸入ガス量Ｍとして仮設定する。
【０１３２】
次いで、ＣＰＵ８１はステップ８６０に進み、図２に示したインテークマニホールドモデルＭ４に従う図示しないルーチンにより推定されている上記特定気筒の吸気弁閉時（上記特定圧縮行程開始時点）の吸気管内温度Ｔｍの値を筒内ガス平均温度Ｔとして仮設定する。そして、ＣＰＵ８１はステップ８６５に進んで、前記筒内吸入ガス量Ｍの値を筒内吸入ガス量Ｍｃｍとして設定するとともに、続くステップ８７０にて前記筒内ガス平均温度Ｔの値を筒内ガス平均温度Ｔｃｍとして設定した後、ステップ８７５に進んで図１０に示した偏差指標値の計算ルーチンを実行する。
【０１３３】
（偏差指標値の計算）
即ち、ＣＰＵ８１はステップ１０００から処理を開始し、ステップ１００５に進んで、先の図８のステップ８６５，８７０にてそれぞれ設定した筒内吸入ガス量Ｍｃｍの値及び筒内ガス平均温度Ｔｃｍの値と、図９のステップ９３５にて既に設定されているシリンダ容積Ｖｃｍ（１）の値と、上記数１２の右辺に相当するステップ１００５内に記載の式とに基いて上記特定圧縮行程開始時点における仮の筒内絶対圧力Ｐ０を算出する。
【０１３４】
次に、ＣＰＵ８１はステップ１０１０に進んで、仮の筒内絶対圧力Ｐ０の値を仮の筒内絶対圧力設定値Ｐｃ（１）として設定するとともに、続くステップ１０１５にて同仮の筒内絶対圧力Ｐ０の値を仮の筒内絶対圧力計算値Ｐｃｍ（１）として設定する。次いで、ＣＰＵ８１は１０２０に進み、筒内絶対圧力設定値Ｐｃ（１）の値と、図９の９３０にて既に設定されている筒内相対圧力検出値Ｐｍｅａｓ（１）と、上記数１３の右辺に相当するステップ１０２０内に記載の式とに基き仮の圧力補正値ΔＰを算出する。
【０１３５】
次に、ＣＰＵ８１はステップ１０２５に進んで、偏差指標値Ｅの値を「０」にクリアするとともに、続くステップ１０３０にて変数Ｉを「１」に設定する。次いで、ステップ１０３５に進んで、図９のステップ９３０にて設定されている筒内相対圧力検出値Ｐｍｅａｓ（Ｉ）の値と、ステップ１０２０にて算出した仮の圧力補正値ΔＰと、上記数１４の右辺に相当するステップ１０３５内に記載の式とに基き仮の筒内絶対圧力設定値Ｐｃ（Ｉ）を算出する。
【０１３６】
次に、ＣＰＵ８１はステップ１０４０に進んで、仮の筒内絶対圧力計算値Ｐｃｍ（Ｉ）の値（現時点では変数Ｉの値は「１」であり、ステップ１０１５にて設定した仮の筒内絶対圧力計算値Ｐｃｍ（１）の値）と、ステップ１００５にて使用した筒内吸入ガス量Ｍｃｍの値と、図９のステップ９３５にて既に設定されているシリンダ容積Ｖｃｍ（Ｉ）の値と、上記数２６の右辺に相当するステップ１０４０内に記載の式とに基き筒内ガス平均温度Ｔｃｍ（Ｉ）を算出する。
【０１３７】
次いで、ＣＰＵ８１はステップ１０４５に進み、前記筒内ガス平均温度Ｔｃｍ（Ｉ）の値と、筒内ガス平均温度Ｔｃｍの関数ｆ５とに基いて得られる比熱比を比熱比κ（Ｉ）として設定し、続くステップ１０５０にて、図８のステップ８５０にて設定した平均ピストン速度Ｃｗの値と、ステップ１０４０にて使用した仮の筒内絶対圧力計算値Ｐｃｍ（Ｉ）の値と、ステップ１０４０にて算出した筒内ガス平均温度Ｔｃｍ（Ｉ）の値と、上記数２５の右辺に相当するステップ１０５０内に記載の式とに基いて熱伝達率ｈｗ（Ｉ）を算出する。
【０１３８】
次に、ＣＰＵ８１はステップ１０５５に進んで、図９のステップ９４５にて設定した筒内表面積Ａ（Ｉ）の値と、前記熱伝達率ｈｗ（Ｉ）の値と、ステップ１０４０にて算出した筒内ガス平均温度Ｔｃｍ（Ｉ）の値と、図８のステップ８４０にて設定した筒内壁面温度Ｔｗの値と、上記数２４の右辺に相当するステップ１０５５内に記載の式とに基いて単位時間当たりの熱伝達量Ｑ（Ｉ）を算出する。
【０１３９】
次いで、ＣＰＵ８１はステップ１０６０に進んで、ステップ１０４５にて設定したκ（Ｉ）の値及びκ（Ｉ−１）の値（現時点では変数Ｉの値は「１」であり、κ（０）の値はκ（１）の値と同一）と、ステップ１０４０にて使用したシリンダ容積Ｖｃｍ（Ｉ）の値と、図９のステップ９４０にて設定したシリンダ容積微分値ｄＶｃｍ（Ｉ）／ｄθと、ステップ１０４０にて使用した仮の筒内絶対圧力計算値Ｐｃｍ（Ｉ）の値と、ステップ１０５５にて算出した熱伝達量Ｑ（Ｉ）の値と、図９のステップ９５０にて設定した角速度ω（Ｉ）の値と、上記数２３の右辺に相当するステップ１０６０内に記載の式とに基いて仮の筒内絶対圧力計算値Ｐｃｍ（Ｉ＋１）を算出する。
【０１４０】
現時点では、変数Ｉの値は「１」であるので、上述したステップ１０３５〜１０６０では、仮の筒内絶対圧力設定値Ｐｃ（１）、筒内ガス平均温度Ｔｃｍ（１）、比熱比κ（１）、熱伝達率ｈｗ（１）、熱伝達量Ｑ（１）、及び仮の筒内絶対圧力計算値Ｐｃｍ（２）が算出・設定される。
【０１４１】
次に、ＣＰＵ８１はステップ１０６５に進み、その時点での偏差指標値Ｅの値に、仮の筒内絶対圧力計算値Ｐｃｍ（Ｉ）から仮の筒内絶対圧力設定値Ｐｃ（Ｉ）を減じた値を二乗した値を加えた値を新たな偏差指標値Ｅとして設定する。現時点では、偏差指標値Ｅの値は「０」であり、変数Ｉの値は「１」であるので、この処理により、偏差指標値Ｅの値は（Ｐｃｍ（１）−Ｐｃ（１））２となる。
【０１４２】
そして、ＣＰＵ８１はステップ１０７０に進んで変数Ｉの値が図８のステップ８３５にて設定したサンプル数Ｎの値と等しくなったか否かを判定する。現時点では変数Ｉの値は「１」であるので、ＣＰＵ８１はステップ１０７０にて「Ｎｏ」と判定してステップ１０７５に進み、変数Ｉの値を「１」だけ増大し、その後ステップ１０３５に戻る。
【０１４３】
この結果、変数Ｉの値は「２」になり、続くステップ１０３５〜１０６０では、仮の筒内絶対圧力設定値Ｐｃ（２）、筒内ガス平均温度Ｔｃｍ（２）、比熱比κ（２）、熱伝達率ｈｗ（２）、熱伝達量Ｑ（２）、及び仮の筒内絶対圧力計算値Ｐｃｍ（３）が算出・設定される。ここで、ステップ１０４０にて使用される仮の筒内絶対圧力計算値Ｐｃｍ（２）としては、直前にステップ１０６０の処理を実行した際に計算した仮の筒内絶対圧力計算値Ｐｃｍ（２）が使用される。
【０１４４】
次に、ＣＰＵ８１は再びステップ１０６５に進み、その時点での偏差指標値Ｅの値に（Ｐｃｍ（２）−Ｐｃ（２））２の値を加えた値を新たな偏差指標値Ｅとして設定する。このような処理は、ステップ１０７５の処理が繰り返されて変数Ｉの値がサンプル数Ｎになるまで繰り返し実行される。これにより、偏差指標値Ｅが上記数２７に示した式に基いた値となる。
【０１４５】
そして、変数Ｉの値がサンプル数Ｎになると、ＣＰＵ８１はステップ１０７０にて「Ｙｅｓ」と判定してステップ１０９５に進み、本ルーチンを終了するとともに、図８のステップ８８０に戻る。このように、図１０の偏差指標値の計算ルーチンは、筒内吸入ガス量Ｍｃｍの値及び筒内ガス平均温度Ｔｃｍの値を設定すると、同設定した筒内吸入ガス量Ｍｃｍの値及び筒内ガス平均温度Ｔｃｍの値に基いた偏差指標値Ｅを計算するルーチンである。
【０１４６】
次いで、ＣＰＵ８１はステップ８８０に進むと、上述した図１０のルーチンにより計算された偏差指標値Ｅの値を偏差指標値前回値Ｅ１に格納するとともに、続くステップ８８５にて筒内相対圧力取得処理実行中フラグＸＨＡＮの値を「０」に設定した後、ステップ８９５に進んで図８のルーチンを一旦終了する。
【０１４７】
以降、ＣＰＵ８１は図８のルーチンをステップ８００から繰り返し実行するが、筒内相対圧力取得処理実行中フラグＸＨＡＮの値が「０」になっているので、ステップ８０５に進んだとき「Ｙｅｓ」と判定してステップ８１０に進むようになる。この時点では、筒内吸入ガス量等の計算を終了した直後であって先に説明した筒内吸入ガス量等の計算開始条件が成立していないので、ＣＰＵ８１はステップ８１０にて「Ｎｏ」と判定してステップ８９５に進み本ルーチンを一旦終了する。この結果、次回の筒内吸入ガス量等の計算開始条件が成立する時点までステップ８００，８０５，８１０，８９５の処理が繰り返し実行されて、筒内相対圧力取得処理実行中フラグＸＨＡＮの値は「０」に維持される。
【０１４８】
（偏差指標値最小化処理）
また、ＣＰＵ８１は図１１及びこれに続く図１２に示した偏差指標値を最小にするためのルーチンを所定時間の経過毎に繰り返し実行するようになっている。従って、所定のタイミングになると、ＣＰＵ８１はステップ１１００から処理を開始し、ステップ１１０２に進んで、筒内相対圧力取得処理実行中フラグＸＨＡＮの値が「１」から「０」に変化したか否かをモニタする。このとき、上記特定気筒の上記特定圧縮行程が終了し、上述した図８のステップ８８５にて筒内相対圧力取得処理実行中フラグＸＨＡＮの値が「１」から「０」に変更されると、ＣＰＵ８１はステップ１１０２にて「Ｙｅｓ」と判定してステップ１１０４に進む。なお、筒内相対圧力取得処理実行中フラグＸＨＡＮの値が変化していなければ、ＣＰＵ８１はステップ１１０２から図１２のステップ１１９５に直接進んで本ルーチンを一旦終了する。
【０１４９】
いま、上記特定気筒の上記特定圧縮行程が終了した直後であるとすると、図８のステップ８８５にて筒内相対圧力取得処理実行中フラグＸＨＡＮの値が「１」から「０」に変更された直後であるから、ＣＰＵ８１はステップ１１０２からステップ１１０４に進み、その時点での筒内吸入ガス量Ｍの値に同定用補正量ΔＭを加えた値を筒内吸入ガス量Ｍ（１）として設定するとともにその時点での筒内ガス平均温度Ｔの値を筒内ガス平均温度Ｔ（１）として設定し、続くステップ１１０６にて、その時点での筒内吸入ガス量Ｍの値から同定用補正量ΔＭを値を減じた値を筒内吸入ガス量Ｍ（２）として設定するとともにその時点での筒内ガス平均温度Ｔの値を筒内ガス平均温度Ｔ（２）として設定し、続くステップ１１０８にて、その時点での筒内吸入ガス量Ｍの値を筒内吸入ガス量Ｍ（３）として設定するとともにその時点での筒内ガス平均温度Ｔの値に同定用補正温度ΔＴを加えた値を筒内ガス平均温度Ｔ（３）として設定し、続くステップ１１１０にて、その時点での筒内吸入ガス量Ｍの値を筒内吸入ガス量Ｍ（４）として設定するとともにその時点での筒内ガス平均温度Ｔの値から同定用補正温度ΔＴを減じた値を筒内ガス平均温度Ｔ（４）として設定する（現時点では、筒内吸入ガス量Ｍの値及び筒内ガス平均温度Ｔの値は、それぞれ図８のステップ８５５及びステップ８６０にて設定された値である。）。
【０１５０】
次に、ＣＰＵ８１はステップ１１１２に進んで変数Ｉを「１」に設定してステップ１１１４に進み、筒内吸入ガス量Ｍ（Ｉ）の値及び筒内ガス平均温度Ｔ（Ｉ）の値をそれぞれ筒内吸入ガス量Ｍｃｍ及び筒内ガス平均温度Ｔｃｍとして設定する。現時点では変数Ｉの値は「１」であるので、筒内吸入ガス量Ｍ（１）の値及び筒内ガス平均温度Ｔ（１）の値がそれぞれ筒内吸入ガス量Ｍｃｍ及び筒内ガス平均温度Ｔｃｍとして設定される。
【０１５１】
次いで、ＣＰＵ８１はステップ１１１６に進んで、先に説明した図１０の偏差指標値の計算ルーチンを実行する。これにより、筒内吸入ガス量Ｍ（１）の値及び筒内ガス平均温度Ｔ（１）の値に基いた偏差指標値Ｅが算出される。次に、ＣＰＵ８１はステップ１１１８に進んで、前記偏差指標値Ｅの値を偏差指標値Ｅ（Ｉ）として設定する。現時点では、変数Ｉの値は「１」であるので前記偏差指標値Ｅが偏差指標値Ｅ（１）として設定される。
【０１５２】
次いで、ＣＰＵ８１はステップ１１２０に進み、変数Ｉの値が「４」と等しいか否かを判定する。現時点では、変数Ｉの値は「１」であるので、ＣＰＵ８１はステップ１１２０にて「Ｎｏ」と判定してステップ１１２２に進み、変数Ｉの値を「１」だけ増大し、その後ステップ１１１４に戻る。この結果、続くステップ１１１４〜１１１８にて、筒内吸入ガス量Ｍ（２）の値及び筒内ガス平均温度Ｔ（２）の値に基いた偏差指標値Ｅが算出され、同偏差指標値Ｅが偏差指標値Ｅ（２）として設定される。このような処理は、ステップ１１２２の処理が繰り返されて変数Ｉの値が「４」になるまで繰り返される。この結果、偏差指標値Ｅ（Ｉ）は筒内吸入ガス量Ｍ（Ｉ）の値及び筒内ガス平均温度Ｔ（Ｉ）に基いて計算される値になるように設定される（Ｉ＝１，２，３，４）。
【０１５３】
そして、変数Ｉの値が「４」になると、ＣＰＵ８１はステップ１１２０にて「Ｙｅｓ」と判定して図１２のステップ１１２４に進み、前記偏差指標値Ｅ（Ｉ）　（Ｉ＝１，２，３，４）のうち最小のものに対応する変数Ｉの値をＩｍｉｎとして設定するとともに、ステップ１１２６に進んで、最小偏差指標値Ｅ（Ｉｍｉｎ）からその時点での偏差指標値前回値Ｅ１（現時点では、図８のステップ８８０にて設定されている値）を減じた値を偏差指標値変化量ΔＥとして設定する。
【０１５４】
次に、ＣＰＵ８１はステップ１１２８に進み、前記偏差指標値変化量ΔＥの絶対値が最小化判定基準値ΔＥｒｅｆ以上であるか否かを判定する。いま、偏差指標値変化量ΔＥの絶対値が最小化判定基準値ΔＥｒｅｆ以上であるとして説明を続けると、ＣＰＵ８１はステップ１１２８にて「Ｙｅｓ」と判定してステップ１１３０に進み、筒内吸入ガス量Ｍ（Ｉｍｉｎ）の値及び筒内ガス平均温度Ｔ（Ｉｍｉｎ）の値を、それぞれ筒内吸入ガス量Ｍ及び筒内ガス平均温度Ｔとして再設定するとともに、続くステップ１１３２にて上記最小偏差指標値Ｅ（Ｉｍｉｎ）を偏差指標値前回値Ｅ１として再設定した後、図１１のステップ１１０４に戻る。
【０１５５】
この結果、続くステップ１１０４〜１１１０にて、筒内吸入ガス量Ｍ（Ｉ）及び筒内ガス平均温度Ｔ（Ｉ）　（Ｉ＝１，２，３，４）の各値が先のステップ１１３０にて設定された筒内吸入ガス量Ｍ及び筒内ガス平均温度Ｔに基いて再設定されて、続くステップ１１１２〜１１２６の処理が実行される。これにより、最小偏差指標値Ｅ（Ｉｍｉｎ）が更新されてより小さい値となり、更新された新たな最小偏差指標値Ｅ（Ｉｍｉｎ）から先のステップ１１３２にて設定された偏差指標値前回値Ｅ１を減じた値が新たな偏差指標値変化量ΔＥとして設定される。
【０１５６】
そして、ステップ１１２８にて再度、前記新たな偏差指標値変化量ΔＥの絶対値が上記最小化判定基準値ΔＥｒｅｆ以上であるか否かが判定され、その結果、新たな偏差指標値変化量ΔＥの絶対値が未だ上記最小化判定基準値ΔＥｒｅｆ以上であれば、同新たな偏差指標値変化量ΔＥの絶対値が上記最小化判定基準値ΔＥｒｅｆ未満になるまで上述した処理が繰り返し実行される。
【０１５７】
その結果、偏差指標値変化量ΔＥの絶対値が上記最小化判定基準値ΔＥｒｅｆ未満になると、その時点での最小偏差指標値Ｅ（Ｉｍｉｎ）が偏差指標値Ｅの最小値であることになる。このとき、ＣＰＵ８１はステップ１１２８にて「Ｎｏ」と判定してステップ１１３４に進み、その時点での筒内吸入ガス量Ｍ（Ｉｍｉｎ）の値及び筒内ガス平均温度Ｔ（Ｉｍｉｎ）の値を、それぞれ確定筒内吸入ガス量Ｍｃｆｉｎ及び確定筒内ガス平均温度Ｔｃｆｉｎとして設定するとともに、続くステップ１１３６にて、同確定筒内吸入ガス量Ｍｃｆｉｎの値及び同確定筒内ガス平均温度Ｔｃｆｉｎの値と、上記数１２の右辺に相当するステップ１１３６内に記載の式とに基き確定筒内絶対圧力Ｐｃｆｉｎを算出する。これにより、筒内吸入ガス量Ｍｃｍ及び筒内ガス平均温度Ｔｃｍの同定が完了する。
【０１５８】
次いで、ＣＰＵ８１はステップ１１３８に進んで、前記確定筒内絶対圧力Ｐｃｆｉｎの値と、図９のステップ９３０にて設定されている筒内相対圧力検出値Ｐｍｅａｓ（１）と、上記数２８の右辺に相当するステップ１１３８内に記載の式とに基いて確定圧力補正値ΔＰｆｉｎを算出する。
【０１５９】
次に、ＣＰＵ８１はステップ１１４０に進んで、上記確定筒内吸入ガス量Ｍｃｆｉｎから実際の燃料噴射量ｆｉを減じることにより上記特定気筒の上記特定圧縮行程時における真の筒内吸入空気量を求め、この真の筒内吸入空気量の値と、上述した各モデルを用いて推定されている同特定気筒の同特定圧縮行程時における筒内吸入空気量Ｍｃの値とを比較することにより、先に説明した吸気弁モデルＭ３を記述する上記数７における比例係数ｃ及び既燃ガス量ｄをそれぞれ求めるためにＲＯＭ８２に格納されている各テーブルを補正する。そして、ＣＰＵ８１はステップ１１９５に進んで本ルーチンを一旦終了する。
【０１６０】
以降、ＣＰＵ８１は図１１及びこれに続く図１２のルーチンを繰り返し実行するが、筒内相対圧力取得処理実行中フラグＸＨＡＮの値が「０」に維持されているので、図１１のステップ１１０２にて「Ｎｏ」と判定して図１２のステップ１１９５に直接進んで本ルーチンを一旦終了するようになる。
【０１６１】
以上説明したように、本発明による筒内ガス状態取得装置の実施形態によれば、筒内相対圧力センサ７２の出力に基いて、上記特定気筒の上記特定圧縮行程開始時点におけるシリンダ内のガスの真の状態を精度良く取得することができた。
【０１６２】
また、上述した筒内吸入ガス量等の計算開始条件が成立する毎に、真の筒内吸入ガス量（Ｍｃｆｉｎ−ｆｉ）を求め、同真の筒内吸入ガス量（Ｍｃｆｉｎ−ｆｉ）と上記各モデルにより推定された筒内吸入空気量Ｍｃとを比較して吸気弁モデルＭ３を記述する上記数７における比例係数ｃ及び既燃ガス量ｄをそれぞれ求めるための各テーブルを補正する。この結果、上記各モデルにより構成される筒内吸入空気量推定装置による筒内吸入空気量Ｍｃの推定精度を向上させることができた。
【０１６３】
また、確定圧力補正値ΔＰｆｉｎが取得された後においては、筒内相対圧力センサ７２の出力値である筒内相対圧力検出値Ｐｍｅａｓと同確定圧力補正値ΔＰｆｉｎとに基いてシリンダ２１内の筒内絶対圧力を精度良く求めることができるようになる。従って、内燃機関の異常燃焼（例えば、ノッキング、ミスファイア等）の発生を精度良く検出できるようになった。
【０１６４】
さらには、最初に仮設定される筒内吸入ガス量Ｍｃｍの値は、上記特定気筒の上記特定圧縮行程開始時点の直前に同特定気筒内に実際に噴射された燃料噴射量ｆｉと上記各モデルを用いて推定されている同特定気筒の上記特定圧縮行程時における筒内吸入空気量Ｍｃとを加算した値であり、最初に仮設定される筒内ガス平均温度Ｔｃｍの値は、上記インテークマニホールドモデルＭ４により推定されている上記特定気筒の吸気弁閉時（上記特定圧縮行程開始時点）の吸気管内温度Ｔｍである。従って、最初に仮設定される筒内吸入ガス量Ｍｃｍの値、及び最初に仮設定される筒内ガス平均温度Ｔｃｍの値は、それぞれ、実際の真の値に近い値である可能性が高いので、偏差指標値Ｅが最小となるように筒内吸入ガス量Ｍｃｍの値及び筒内ガス平均温度Ｔｃｍの値を同定する際に実行される処理の繰り返し回数を少なくすることができた。
【０１６５】
本発明は上記実施形態に限定されることはなく、本発明の範囲内において種々の変形例を採用することができる。例えば、上記実施形態では、仮設定手段は、（特定圧縮行程時にシリンダに吸入されている）筒内吸入ガス量Ｍｃｍと、（特定圧縮行程開始時点における）筒内ガス平均温度Ｔｃｍとを仮設定するように構成されているが、前記筒内吸入ガス量Ｍｃｍ及び前記筒内ガス平均温度Ｔｃｍのいずれか一方と、（特定圧縮行程開始時点における）仮の筒内絶対圧力Ｐ０とを仮設定するように構成してもよい。この場合、前記筒内吸入ガス量Ｍｃｍ及び前記筒内ガス平均温度Ｔｃｍの他方の値は、上記数１２を変形した式に基いて求めることができ、これにより、上記特定圧縮行程開始時点における上記特定気筒のシリンダ内のガスの状態が仮設定される。
【０１６６】
また、上記実施形態では、複数の（サンプル数Ｎ個の）仮の筒内絶対圧力設定値Ｐｃ（Ｈ）　（Ｈ＝１，２，・・・，Ｎ−１，Ｎ）等が設定される圧縮行程の所定の区間の終期を特定圧縮行程終了時点（点火プラグにより点火される時点）に設定しているが、内燃機関１０の運転状態に応じて所定の期間（例えば、アクセルペダルの操作量Ａｃｃｐが「０」となる期間）だけ燃料を噴射しない制御（所謂フューエルカット制御）が実行される場合、前記所定の区間の終期を、通常は特定圧縮行程終了時点に設定し、前記所定の期間内はピストンの位置が圧縮上死点に到達した時点に設定してもよい。
【０１６７】
また、上記実施形態では、エネルギー保存則、及び状態方程式に基くシリンダモデルを記述する上記数１５〜上記数１８のうちの熱伝達率ｈｗ（熱伝達量Ｑ）を求めるための上記数１７において、定数Ｃ１を一定値「２．２８」に設定しているが、値Ｃ１を仮設定手段により仮設定し、筒内ガス状態取得手段により同値Ｃ１の真の値を求めるように（偏差指標値Ｅが最小となるように同値Ｃ１を同定するように）構成してもよい。
【０１６８】
また、上記実施形態では、第１計算手段は、筒内ガス状態取得手段が仮設定手段、第１計算手段、第２計算手段、及び偏差指標値計算手段が行う一連の各処理を繰返し実行させる前に予め取得してある圧縮行程の所定の区間内における複数の筒内相対圧力検出値Ｐｍｅａｓ（Ｈ）　（Ｈ＝１，２，・・・，Ｎ−１，Ｎ）に基いて、複数の仮の筒内絶対圧力設定値Ｐｃ（Ｈ）　（Ｈ＝１，２，・・・，Ｎ−１，Ｎ）を設定するように構成されているが、、前記第１計算手段を、前記筒内ガス状態取得手段が前記一連の各処理を繰返し実行する毎に、前記複数の筒内相対圧力検出値Ｐｍｅａｓ（Ｈ）　（Ｈ＝１，２，・・・，Ｎ−１，Ｎ）を新たに取得して同新たに取得した複数の筒内相対圧力検出値Ｐｍｅａｓ（Ｈ）　（Ｈ＝１，２，・・・，Ｎ−１，Ｎ）に基いて前記複数の仮の筒内絶対圧力設定値Ｐｃ（Ｈ）　（Ｈ＝１，２，・・・，Ｎ−１，Ｎ）を設定するように構成してもよい。
【図面の簡単な説明】
【図１】本発明による筒内ガス状態取得装置を含む燃料噴射量制御装置を火花点火式多気筒内燃機関に適用したシステムの概略構成図である。
【図２】図１に示した電気制御装置が筒内吸入空気量を推定するために採用した各種モデルの接続関係を示した機能ブロック図である。
【図３】図１に示したＣＰＵが参照するアクセルペダル操作量と目標スロットルバルブ開度との関係を規定したテーブルを示す図である。
【図４】図１に示した特定気筒の圧縮行程におけるクランク角に対する、複数の筒内相対圧力検出値、複数の仮の筒内絶対圧力設定値、及び複数の仮の筒内絶対圧力計算値の各々に基く波形の一例を示した図である。
【図５】シリンダモデルを表すために使用する変数を説明するためシリンダ及びその近傍を概念的に示した図である。
【図６】図１に示したＣＰＵが実行するスロットルバルブ開度を制御するためのルーチンを示したフローチャートである。
【図７】図１に示したＣＰＵが実行する燃料噴射量を計算するためのルーチンを示したフローチャートである。
【図８】図１に示したＣＰＵが実行する筒内吸入ガス量の計算を開始するか否かを決定するためのルーチンを示したフローチャートである。
【図９】図１に示したＣＰＵが実行する複数の筒内相対圧力検出値等を取得するためのルーチンを示したフローチャートである。
【図１０】図１に示したＣＰＵが実行する偏差指標値を計算するためのルーチンを示したフローチャートである。
【図１１】図１に示したＣＰＵが実行する偏差指標値を最小化するためのルーチンの前半部を示したフローチャートである。
【図１２】図１に示したＣＰＵが実行する偏差指標値を最小化するためのルーチンの後半部を示したフローチャートである。
【符号の説明】
１０…火花点火式多気筒内燃機関、２０…シリンダブロック部（エンジン本体部）、２１…シリンダ、２５…燃焼室、３１…吸気ポート、３２…吸気弁、３９…インジェクタ、４１…吸気管、４３…スロットルバルブ、６８…クランクポジションセンサ、７１…アクセル開度センサ、７２…筒内相対圧力センサ、８０…電気制御装置、８１…ＣＰＵ。
【図１】
2004108348.tif 000017
【図２】
2004108348.tif 000018
【図３】
2004108348.tif 000019
【図４】
2004108348.tif 000020
【図５】
2004108348.tif 000021
【図６】
2004108348.tif 000022
【図７】
2004108348.tif 000023
【図８】
2004108348.tif 000024
【図９】
2004108348.tif 000025
【図１０】
2004108348.tif 000026
【図１１】
2004108348.tif 000027
【図１２】
2004108348.tif 000028