(57)【要約】【課題】吸気管内の脈動による影響も含めて筒内充填空気量を正確に推定する。【解決手段】筒内充填空気量を推定する際には、吸気管内圧力の脈動によって筒内に充填される空気量（余分空気量）も考慮される。余分空気量は吸気バルブ９の開弁による吸気管内圧力の降下量に基づき求められる。吸気管内圧力は、圧力センサ２２からの出力をＡ／Ｄ変換によりサンプリング周期Ｘ１で量子化して得られる値であり、量子化誤差やノイズによる精度悪化は避けられないものである。こうした精度悪化を抑制すべく、余分空気量の算出に用いられる降下量として、吸気バルブ９の開閉サイクル複数回分の上記降下量を平均化したもの（平均値ΔＰｍdwnave）が採用される。また、サンプリング周期Ｘ１に対する吸気管３内の圧力の脈動周期の長さが所定レベル未満である旨判断されたときには、通常よりも高精度の平均化法により平均値ΔＰｍdwnaveが求められる。【選択図】      図１
2005282563.tif 000002
【特許請求の範囲】
【請求項１】
吸気バルブの開閉に伴い吸気管内に圧力の脈動が生じる内燃機関に適用され、前記吸気管に流入する空気の量を流入空気量として算出し、その流入空気量を用いて筒内充填空気量を推定する内燃機関の筒内充填空気量推定装置において、
前記吸気管の圧力を検出する圧力センサからの出力をＡ／Ｄ変換により所定のサンプリング周期で量子化し、その量子化後の値を前記吸気管の圧力として用いて前記吸気バルブの開弁による前記吸気管内の圧力の降下量を算出する降下量算出手段と、
前記圧力の降下量に基づき、当該圧力の脈動によって筒内に充填される空気の量を余分空気量として算出する余分空気量算出手段と、
前記流入空気量及び前記余分空気量に基づき前記筒内充填空気量を推定する筒内充填空気量推定手段とを備え、
前記降下量算出手段は、前記吸気バルブの開閉サイクル複数回分の前記降下量の平均値を前記余分空気量算出用の降下量とするものであって、前記Ａ／Ｄ変換のサンプリング周期に対する前記圧力の脈動周期の長さが所定レベル未満であるときには、そうでないときより高精度の平均化法によって前記降下量の平均値を求めるものである
ことを特徴とする内燃機関の筒内充填空気量推定装置。
【請求項２】
前記降下量算出手段は、
前記サンプリング周期に対する前記圧力の脈動周期の長さが所定レベル以上であるときには、前記量子化後の圧力を用いて吸気バルブ開弁中の同圧力の最大値と最小値との差を当該吸気バルブの一回の開閉サイクルでの前記降下量として算出するとともにメモリに記憶し、この降下量について開閉サイクル複数回分の平均をとって得られる平均値を前記余分空気量算出用の降下量とするものであって、
前記サンプリング周期に対する前記圧力の脈動周期の長さが所定レベル未満であるときには、サンプリング周期毎に量子化後の圧力と、そのときの時刻とをメモリに記憶し、それら圧力及び時刻についての前記吸気バルブの開閉サイクル複数回分のデータから前記圧力の開閉サイクル一回分の平均波形を求め、その平均波形における圧力の最大値と最小値との差を前記降下量の前記開閉サイクル複数回分の平均値とし前記余分空気量算出用の降下量とするものである
請求項１記載の内燃機関の筒内充填空気量推定装置。
【請求項３】
吸気バルブの開閉に伴い吸気管内に圧力の脈動が生じる内燃機関に適用され、前記吸気管に流入する空気の量を流入空気量として算出し、その流入空気量を用いて筒内充填空気量を推定する内燃機関の筒内充填空気量推定装置において、
前記吸気管の圧力を検出する圧力センサからの出力をＡ／Ｄ変換により所定のサンプリング周期で量子化し、その量子化後の値を前記吸気管の圧力として用いて前記吸気バルブの開弁による前記吸気管内の圧力の降下量を算出する降下量算出手段と、
前記圧力の降下量に基づき、当該圧力の脈動によって筒内に充填される空気の量を余分空気量として算出する余分空気量算出手段と、
前記流入空気量及び前記余分空気量に基づき前記筒内充填空気量を推定する筒内充填空気量推定手段とを備え、
前記降下量算出手段は、前記吸気バルブの開閉サイクル複数回分の前記降下量の平均値を前記余分空気量算出用の降下量とするものであって、前記Ａ／Ｄ変換のサンプリング周期に対する前記圧力の脈動周期の長さに応じて、前記降下量の平均値を求めるための開閉サイクル回数を可変とする
ことを特徴とする内燃機関の筒内充填空気量推定装置。
【請求項４】
吸気バルブの開閉に伴い吸気管内に圧力の脈動が生じる内燃機関に適用され、前記吸気管に流入する空気の量を流入空気量として算出し、その流入空気量を用いて筒内充填空気量を推定する内燃機関の筒内充填空気量推定装置において、
前記吸気管の圧力を検出する圧力センサからの出力をＡ／Ｄ変換により所定のサンプリング周期で量子化し、その量子化後の値を前記吸気管の圧力として用いて前記吸気バルブの開弁による前記吸気管内の圧力の降下量を算出する降下量算出手段と、
前記圧力の降下量に基づき、当該圧力の脈動によって筒内に充填される空気の量を余分空気量として算出する余分空気量算出手段と、
前記流入空気量及び前記余分空気量に基づき前記筒内充填空気量を推定する筒内充填空気量推定手段とを備え、
前記降下量算出手段は、前記吸気バルブの開閉サイクル複数回分の前記降下量の平均値を前記余分空気量算出用の降下量とするものであって、前記Ａ／Ｄ変換のサンプリング周期に対する前記圧力の脈動周期の長さに応じて、当該サンプリング周期を変更する
ことを特徴とする内燃機関の筒内充填空気量推定装置。
【請求項５】
前記降下量算出手段は、前記Ａ／Ｄ変換のサンプリング周期に対する前記圧力の脈動周期の長さが所定レベル未満であるときには、前記サンプリング周期を前記脈動周期の長さが所定レベル以上であるときよりも短くする
請求項４記載の内燃機関の筒内充填空気量推定装置。
【請求項６】
前記降下量算出手段は、前記Ａ／Ｄ変換のサンプリング周期に対する前記圧力の脈動周期の長さが、前記開閉サイクル毎における前記量子化後の吸気管の圧力の推移を同期させる値である旨判断されるとき、前記サンプリング周期をそれまでの値と異なる値に変更する
請求項４又は５記載の内燃機関の筒内充填空気量推定装置。
【請求項７】
吸気バルブの開閉に伴い吸気管内に圧力の脈動が生じる内燃機関に適用され、前記吸気管に流入する空気の量を流入空気量として算出し、その流入空気量を用いて筒内充填空気量を推定する内燃機関の筒内充填空気量推定装置において、
前記吸気管の圧力を検出する圧力センサからの出力をＡ／Ｄ変換により所定のサンプリング周期で量子化し、その量子化後の値を前記吸気管の圧力として用いて前記吸気バルブの開弁による前記吸気管内の圧力の降下量を算出する降下量算出手段と、
前記圧力の降下量に基づき、当該圧力の脈動によって筒内に充填される空気の量を余分空気量として算出する余分空気量算出手段と、
前記流入空気量及び前記余分空気量に基づき前記筒内充填空気量を推定する筒内充填空気量推定手段とを備え、
前記降下量算出手段は、前記吸気バルブの開閉サイクル複数回分の前記降下量の平均値を前記余分空気量算出用の降下量とするものであって、前記Ａ／Ｄ変換のサンプリング周期に対する前記圧力の脈動周期の長さが、前記吸気バルブの開閉サイクル毎における前記量子化後の吸気管の圧力の推移を同期させる値である旨判断されるとき、当該開閉サイクルでの前記吸気管内の圧力の降下量の算出を行わないものである
ことを特徴とする内燃機関の筒内充填空気量推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、内燃機関の筒内充填空気量推定装置に関するものである。
【背景技術】
【０００２】
内燃機関においては、吸気バルブ閉弁時に筒内に存在する空気の量（以下、筒内充填空気量という）に基づき、燃料噴射量制御などの運転制御が行われるため、その筒内充填空気量を正確に推定することが適切な運転制御を行う上で重要になる。
【０００３】
このため、例えば特許文献１では、吸気効率を気筒毎に推定するとともに、それら吸気効率を各々対応する気筒の筒内充填空気量の推定に用いるようにしている。この場合、内燃機関の吸気系でのデポジットの付着や部品の組み付け誤差等により気筒毎に吸気効率に違いが生じ、それに起因して各気筒毎に筒内充填空気量がばらついたとしても、そのばらつきも考慮して筒内充填空気量を推定することができる。
【特許文献１】特開２００１−２３４７９８公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
ところで、内燃機関においては、各気筒の吸気バルブが順々に開弁し、それによって吸気管内の圧力に脈動が生じることとなる。こうした吸気管内の圧力の脈動は、各気筒における筒内充填空気量に密接に関係している。しかし、上記特許文献１では、各気筒の筒内充填空気量を推定する際に上記脈動を考慮していないため、推定される各気筒の筒内充填空気量が上記脈動の影響を受ける分だけ不正確なものとなる。
【０００５】
本発明はこのような実情に鑑みてなされたものであって、その目的は、吸気管内の脈動による影響も含めて筒内充填空気量を正確に推定することのできる内燃機関の筒内充填空気量推定装置を提供することにある。
【課題を解決するための手段】
【０００６】
以下、上記目的を達成するための手段及びその作用効果について記載する。
上記目的を達成するため、請求項１記載の発明では、吸気バルブの開閉に伴い吸気管内に圧力の脈動が生じる内燃機関に適用され、前記吸気管に流入する空気の量を流入空気量として算出し、その流入空気量を用いて筒内充填空気量を推定する内燃機関の筒内充填空気量推定装置において、前記吸気管の圧力を検出する圧力センサからの出力をＡ／Ｄ変換により所定のサンプリング周期で量子化し、その量子化後の値を前記吸気管の圧力として用いて前記吸気バルブの開弁による前記吸気管内の圧力の降下量を算出する降下量算出手段と、前記圧力の降下量に基づき、当該圧力の脈動によって筒内に充填される空気の量を余分空気量として算出する余分空気量算出手段と、前記流入空気量及び前記余分空気量に基づき前記筒内充填空気量を推定する筒内充填空気量推定手段とを備え、前記降下量算出手段は、前記吸気バルブの開閉サイクル複数回分の前記降下量の平均値を前記余分空気量算出用の降下量とするものであって、前記Ａ／Ｄ変換のサンプリング周期に対する前記圧力の脈動周期の長さが所定レベル未満であるときには、そうでないときより高精度の平均化法によって前記降下量の平均値を求めるものとした。
【０００７】
上記構成によれば、吸気バルブの開閉に伴う吸気管内の圧力の脈動によって筒内に充填される空気の量、即ち余分空気量を考慮して筒内充填空気量が推定されるため、その推定される筒内充填空気量を上記脈動による影響も含めた正確なものとすることができる。
【０００８】
また、吸気バルブの開弁による吸気管内の圧力の降下量の算出は、圧力センサからの出力をＡ／Ｄ変換により所定のサンプリング周期で量子化し、量子化後の値を吸気管の圧力として用いることで行われる。ただし、Ａ／Ｄ変換後に取得された圧力は量子化誤差やノイズによって精度が悪化するため、それに伴い余分空気量算出用の降下量の精度も悪化する。このため、吸気バルブの開閉サイクル複数回分の降下量の平均値を余分空気量算出用の降下量として採用し、上述した降下量の精度悪化を抑制するようにしている。
【０００９】
ところで、吸気バルブの開閉態様の変化等によって吸気管内の圧力の脈動周期が変化し、その脈動周期がＡ／Ｄ変換のサンプリング周期に近づくと、吸気バルブの一回の開閉サイクル中におけるＡ／Ｄ変換のサンプリング回数が少なくなる。このサンプリング回数の少なさに起因して、量子化後に取得された圧力の吸気バルブ開弁中の最大値及び最小値が実際の値と大きく食い違うようになり、それら最大値と最小値の差である上記降下量の必要精度確保が困難になる。このため、当該降下量の平均値である余分空気量算出用の降下量を、より高精度の平均化法によって求めることで、余分空気量算出用の降下量の必要精度を確保することも考えられるが、こうした高精度の平均化法の採用は計算負荷の増大やメモリ容量の消費を招く。
【００１０】
そこで、請求項１記載の発明では、上記サンプリング周期に対し上記圧力の脈動周期の長さが所定レベル未満であるときだけ、高精度の平均化法を用いて平均余分空気量算出用の降下量を求めるようにしている。一方、上記サンプリング周期に対し脈動周期の長さが所定レベル以上であるときには、計算負荷やメモリ容量の消費が上記高精度の平均化法よりも少ない平均化法で平均余分空気量算出用の降下量が求められる。これにより、余分空気量算出用の降下量の必要精度を確保しつつ、計算負荷の増大やメモリ容量の消費を最小限にとどめることができるようになる。
【００１１】
請求項２記載の発明では、請求項１記載の発明において、前記降下量算出手段は、前記サンプリング周期に対する前記圧力の脈動周期の長さが所定レベル以上であるときには、前記量子化後の圧力を用いて前記吸気バルブ開弁中の同圧力の最大値と最小値との差を当該吸気バルブの一回の開閉サイクルでの前記降下量として算出するとともにメモリに記憶し、この降下量について前記開閉サイクル複数回分の平均をとって得られる平均値を前記余分空気量算出用の降下量とするものであって、前記サンプリング周期に対する前記圧力の脈動周期の長さが所定レベル未満であるときには、サンプリング周期毎に量子化後の圧力と、そのときの時刻とをメモリに記憶し、それら圧力及び時刻についての前記吸気バルブの開閉サイクル複数回分のデータから前記圧力の前記開閉サイクル一回分の平均波形を求め、その平均波形における圧力の最大値と最小値との差を前記降下量の前記開閉サイクル複数回分の平均値とし前記余分空気量算出用の降下量とするものとした。
【００１２】
上記サンプリング周期に対する上記圧力の脈動周期の長さが所定レベル未満であるときには、降下量の平均値を求めるのに高精度の平均化法が用いられる。この高精度の平均化法では、吸気バルブの開閉サイクル複数回分についてのサンプリング周期毎の量子化後の圧力及びそのときの時刻から、当該圧力の開閉サイクル一回分の圧力センサからの出力の平均波形が求められ、その波形における圧力の最大値と最小値との差を上記降下量の平均値としている。このように平均波形から降下量の平均値を求めることで、その平均値を余分空気量算出用の降下量として精度の高いものとすることができる。ただし、上記平均波形を求めるのには、上記圧力及び時刻に関する多数のデータを用いた複雑な計算が必要になる。このため、計算負荷が増大するとともに、それら多数のデータを記憶するのにメモリ容量を大幅に消費することは避けられない。
【００１３】
このため、上記サンプリング周期に対する上記圧力の脈動周期の長さが所定レベル以上であるときには、降下量の平均値を求めるのに通常精度の平均化法が用いられる。この通常精度の平均化法では、吸気バルブの一回の開閉サイクル毎に、量子化後の圧力を用いて吸気バルブ開弁中の同圧力の最大値と最小値との差を当該圧力の降下量として求め、この降下量について吸気バルブの開閉サイクル複数回分の平均をとって当該降下量の平均値を得るようにしている。この場合、上述した高精度の平均化法のような平均波形を求めることは行われず、一回の開閉サイクル毎の降下量の計算及びその計算結果の記憶、並びに記憶された開閉サイクル複数回分の降下量の平均といった処理が行われるだけであるため、計算負荷の増大やメモリ容量の消費は抑制されるようになる。
【００１４】
以上のように、サンプリング周期に対する圧力の脈動周期の長さが所定レベル未満であるか否かに基づき、上述した高精度の平均化法と通常精度の平均化法とを切り換えて実施することで、余分空気量算出用の降下量の必要精度を確保しつつ、計算負荷の増大やメモリ容量の消費を最小限にとどめることができるようになる。
【００１５】
請求項３記載の発明では、吸気バルブの開閉に伴い吸気管内に圧力の脈動が生じる内燃機関に適用され、前記吸気管に流入する空気の量を流入空気量として算出し、その流入空気量を用いて筒内充填空気量を推定する内燃機関の筒内充填空気量推定装置において、前記吸気管の圧力を検出する圧力センサからの出力をＡ／Ｄ変換により所定のサンプリング周期で量子化し、その量子化後の値を前記吸気管の圧力として用いて前記吸気バルブの開弁による前記吸気管内の圧力の降下量を算出する降下量算出手段と、前記圧力の降下量に基づき、当該圧力の脈動によって筒内に充填される空気の量を余分空気量として算出する余分空気量算出手段と、前記流入空気量及び前記余分空気量に基づき前記筒内充填空気量を推定する筒内充填空気量推定手段とを備え、前記降下量算出手段は、前記吸気バルブの開閉サイクル複数回分の前記降下量の平均値を前記余分空気量算出用の降下量とするものであって、前記Ａ／Ｄ変換のサンプリング周期に対する前記圧力の脈動周期の長さに応じて、前記降下量の平均値を求めるための開閉サイクル回数を可変とするものとした。
【００１６】
上記構成によれば、吸気バルブの開閉に伴う吸気管内の圧力の脈動によって筒内に充填される空気の量、即ち余分空気量を考慮して筒内充填空気量が推定されるため、その推定される筒内充填空気量を上記脈動による影響も含めた正確なものとすることができる。
【００１７】
また、吸気バルブの開弁による吸気管内の圧力の降下量の算出は、圧力センサからの出力をＡ／Ｄ変換により所定のサンプリング周期で量子化し、量子化後の値を吸気管の圧力として用いることで行われる。ただし、Ａ／Ｄ変換後に取得された圧力は量子化誤差やノイズによって精度が悪化するため、それに伴い余分空気量算出用の降下量の精度も悪化する。このため、吸気バルブの開閉サイクル複数回分の降下量の平均値を余分空気量算出用の降下量として採用し、上述した降下量の精度悪化を抑制するようにしている。
【００１８】
ところで、吸気バルブの開閉態様の変化等によって吸気管内の圧力の脈動周期が変化し、その脈動周期がＡ／Ｄ変換のサンプリング周期に近づくと、吸気バルブの一回の開閉サイクル中におけるＡ／Ｄ変換のサンプリング回数が少なくなる。このサンプリング回数の少なさに起因して、量子化後に取得された圧力の吸気バルブ開弁中の最大値及び最小値が実際の値と大きく食い違うようになり、それら最大値と最小値の差である上記降下量の必要精度確保が困難になる。このため、上記降下量の平均値を求めるための開閉サイクル回数を多めの値に設定することで、上記降下量の必要精度を確保することも考えられるが、その開閉サイクル回数を多くするほど上記降下量の平均値を算出するための計算負荷やメモリ容量の消費が増大する。
【００１９】
請求項３記載の発明では、上記Ａ／Ｄ変換のサンプリング周期に対する上記圧力の脈動周期の長さに応じて、降下量の平均値を求めるための開閉サイクル回数を可変とすることができる。ただし、上記開閉サイクル回数が多くなると、平均値を求める際の計算負荷や当該計算に用いられるデータの記憶に起因するメモリ容量の消費が大となることは上述したとおりである。このことを考慮し、サンプリング周期に対する圧力の脈動周期の長さが所定レベル未満になったときだけ、上記開閉サイクル回数を多くすることで、余分空気量算出用の降下量の必要精度を確保しつつ、計算負荷の増大やメモリ容量の消費を最小限にとどめることができる。
【００２０】
請求項４記載の発明では、吸気バルブの開閉に伴い吸気管内に圧力の脈動が生じる内燃機関に適用され、前記吸気管に流入する空気の量を流入空気量として算出し、その流入空気量を用いて筒内充填空気量を推定する内燃機関の筒内充填空気量推定装置において、前記吸気管の圧力を検出する圧力センサからの出力をＡ／Ｄ変換により所定のサンプリング周期で量子化し、その量子化後の値を前記吸気管の圧力として用いて前記吸気バルブの開弁による前記吸気管内の圧力の降下量を算出する降下量算出手段と、前記圧力の降下量に基づき、当該圧力の脈動によって筒内に充填される空気の量を余分空気量として算出する余分空気量算出手段と、前記流入空気量及び前記余分空気量に基づき前記筒内充填空気量を推定する筒内充填空気量推定手段とを備え、前記降下量算出手段は、前記吸気バルブの開閉サイクル複数回分の前記降下量の平均値を前記余分空気量算出用の降下量とするものであって、前記Ａ／Ｄ変換のサンプリング周期に対する前記圧力の脈動周期の長さに応じて、当該サンプリング周期を変更するものとした。
【００２１】
上記構成によれば、吸気バルブの開閉に伴う吸気管内の圧力の脈動によって筒内に充填される空気の量、即ち余分空気量を考慮して筒内充填空気量が推定されるため、その推定される筒内充填空気量を上記脈動による影響も含めた正確なものとすることができる。
【００２２】
また、吸気バルブの開弁による吸気管内の圧力の降下量の算出は、圧力センサからの出力をＡ／Ｄ変換により所定のサンプリング周期で量子化し、量子化後の値を吸気管の圧力として用いることで行われる。ただし、Ａ／Ｄ変換後に取得された圧力は量子化誤差やノイズによって精度が悪化するため、それに伴い余分空気量算出用の降下量の精度も悪化する。このため、吸気バルブの開閉サイクル複数回分の降下量の平均値を余分空気量算出用の降下量として採用し、上述した降下量の精度悪化を抑制するようにしている。
【００２３】
ところで、吸気バルブの開閉態様の変化等による吸気管内の圧力の脈動周期の変化は、余分空気量算出用の降下量の精度に大きな影響を及ぼす。
例えば、脈動周期がＡ／Ｄ変換のサンプリング周期に近づくと、吸気バルブの一回の開閉サイクル中におけるＡ／Ｄ変換のサンプリング回数が少なくなる。このサンプリング回数の少なさに起因して、量子化後に取得された圧力の吸気バルブ開弁中の最大値及び最小値が実際の値と大きく食い違うようになり、それら最大値と最小値の差である上記降下量の必要精度確保が困難になる。
【００２４】
また、所定のサンプリング周期での量子化後の吸気管内の圧力の推移は、吸気バルブの開閉サイクル毎の圧力センサからの出力波形の違いに応じて、開閉サイクル毎に異なる状態（非同期状態）となるのが普通である。しかし、サンプリング周期に対する上記脈動周期の長さによっては、開閉サイクル毎における上記圧力の推移が、開閉サイクル毎の圧力センサの出力波形の違いに関係なく同じ状態（同期状態）となることがある。言い換えれば、いかなるサンプリング周期においても、圧力センサの出力波形の開閉サイクル毎の違いに関係なく、上記圧力の推移が複数の開閉サイクルで同期した状態となる脈動周期が存在するということである。このように上記圧力の推移が複数の開閉サイクルで同期すると、同圧力の降下量を開閉サイクル複数回分を平均したとしても、その平均値が余分空気量算出用の降下量として精度のよい値にならず、同降下量の必要精度確保が困難になる。
【００２５】
請求項４記載の発明では、上記Ａ／Ｄ変換のサンプリング周期に対する上記圧力の脈動周期の長さに応じて、当該サンプリング周期を変更することができるため、上述したような状況が生じないようにサンプリング周期を変更することができる。即ち、吸気バルブの開閉状態に応じて変化する脈動周期がサンプリング周期に過度に近づいた状態とならないよう、また、そのときの脈動周期に対し上記圧力の推移が吸気バルブの複数の開閉サイクルで同期した状態となるサンプリング周期とならないよう、当該サンプリング周期を変更することができる。こうしたサンプリング周期の変更により、上述した状況下でも、余分空気量算出用の降下量の必要精度を確保することができるようになる。
【００２６】
請求項５記載の発明では、請求項４記載の発明において、前記降下量算出手段は、前記Ａ／Ｄ変換のサンプリング周期に対する前記圧力の脈動周期の長さが所定レベル未満であるときには、前記サンプリング周期を前記脈動周期の長さが所定レベル以上であるときよりも短くするものとした。
【００２７】
サンプリング周期を短くすると、吸気バルブの一回の開閉サイクル中のサンプリング回数が多くなるため、そのサンプリング回数の少なさに起因して上記降下量の必要精度確保が困難になるのを抑制することはできるが、そのサンプリング回数が多くなると計算負荷が増大する。上記構成によれば、サンプリング周期に対する圧力の脈動周期の長さが所定レベル未満になったときだけ、サンプリング周期が短くされて一回の開閉サイクル中のサンプリング回数が多くされる。一方、上記サンプリング周期に対し脈動周期の長さが所定レベル以上であるときには、サンプリング周期が短くされることはないため、一回の開閉サイクル中のサンプリング回数が必要以上に多くされることはない。以上により、余分空気量算出用の降下量の必要精度を確保しつつ、上記計算負荷の増大を最小限にとどめることができるようになる。
【００２８】
請求項６記載の発明では、請求項４又は５記載の発明において、前記降下量算出手段は、前記Ａ／Ｄ変換のサンプリング周期に対する前記圧力の脈動周期の長さが、前記開閉サイクル毎における前記量子化後の吸気管の圧力の推移を同期させる値である旨判断されるとき、前記サンプリング周期をそれまでの値と異なる値に変更するものとした。
【００２９】
サンプリング周期と脈動周期との関係によっては、前記吸気バルブの開閉サイクル毎における量子化後の吸気管の圧力の推移が同期した状態になる場合がある。上記構成によれば、このような場合にはサンプリング周期がそれまでとは異なる値に変更され、開閉サイクル毎における量子化後の吸気管の圧力の推移が同期しないようにされる。従って、開閉サイクル毎における上記圧力の推移が同期し、余分空気量算出用の降下量の必要精度確保が困難になるのを抑制することができる。
【００３０】
請求項７記載の発明では、吸気バルブの開閉に伴い吸気管内に圧力の脈動が生じる内燃機関に適用され、前記吸気管に流入する空気の量を流入空気量として算出し、その流入空気量を用いて筒内充填空気量を推定する内燃機関の筒内充填空気量推定装置において、前記吸気管の圧力を検出する圧力センサからの出力をＡ／Ｄ変換により所定のサンプリング周期で量子化し、その量子化後の値を前記吸気管の圧力として用いて前記吸気バルブの開弁による前記吸気管内の圧力の降下量を算出する降下量算出手段と、前記圧力の降下量に基づき、当該圧力の脈動によって筒内に充填される空気の量を余分空気量として算出する余分空気量算出手段と、前記流入空気量及び前記余分空気量に基づき前記筒内充填空気量を推定する筒内充填空気量推定手段とを備え、前記降下量算出手段は、前記吸気バルブの開閉サイクル複数回分の前記降下量の平均値を前記余分空気量算出用の降下量とするものであって、前記Ａ／Ｄ変換のサンプリング周期に対する前記圧力の脈動周期の長さが、前記吸気バルブの開閉サイクル毎における前記量子化後の吸気管の圧力の推移を同期させる値である旨判断されるとき、当該開閉サイクルでの前記吸気管内の圧力の降下量の算出を行わないものとした。
【００３１】
上記構成によれば、吸気バルブの開閉に伴う吸気管内の圧力の脈動によって筒内に充填される空気の量、即ち余分空気量を考慮して筒内充填空気量が推定されるため、その推定される筒内充填空気量を上記脈動による影響も含めた正確なものとすることができる。
【００３２】
また、吸気バルブの開弁による吸気管内の圧力の降下量の算出は、圧力センサからの出力をＡ／Ｄ変換により所定のサンプリング周期で量子化し、量子化後の値を吸気管の圧力として用いることで行われる。ただし、Ａ／Ｄ変換後に取得された圧力は量子化誤差やノイズによって精度が悪化するため、それに伴い余分空気量算出用の降下量の精度も悪化する。このため、吸気バルブの開閉サイクル複数回分の降下量の平均値を余分空気量算出用の降下量として採用し、上述した降下量の精度悪化を抑制するようにしている。
【００３３】
ところで、所定のサンプリング周期での量子化後の吸気管内の圧力の推移は、吸気バルブの開閉サイクル毎の圧力センサからの出力波形の違いに応じて、開閉サイクル毎に異なる状態（非同期状態）となるのが普通である。しかし、サンプリング周期に対する脈動周期の長さによっては、開閉サイクル毎における上記圧力の推移が、開閉サイクル毎の圧力センサの出力波形の違いに関係なく、同じ状態（同期状態）となることがある。言い換えれば、いかなるサンプリング周期においても、圧力センサの出力波形の開閉サイクル毎の違いに関係なく、上記圧力の推移が複数の開閉サイクルで同期した状態となる脈動周期が存在するということである。このように開閉サイクル毎における上記圧力の推移が同期すると、同圧力の降下量を開閉サイクル複数回分を平均したとしても、その平均値が余分空気量算出用の降下量として精度のよい値にはならず、同降下量の必要精度確保が困難になる。
【００３４】
請求項７記載の発明では、上述した状況下では、当該開閉サイクルでの前記吸気管内の圧力の降下量の算出を行わないようにしている。この場合、吸気管内の圧力の降下量が当該開閉サイクルよりも以前の開閉サイクル複数回分について平均され、その平均値が余分空気量算出用の降下量として用いられる。従って、上記圧力の推移が開閉サイクル毎で同期する状況となった開閉サイクルでの吸気管内の圧力の降下量が、余分空気量算出用の降下量を算出するための平均化に用いられることはなく、当該降下量の必要精度確保が困難になるのを抑制することができる。
【発明を実施するための最良の形態】
【００３５】
［第１実施形態］
以下、本発明を自動車に搭載される筒内噴射火花点火式の四気筒エンジンに適用した第１実施形態を図１〜図１０に従って説明する。
【００３６】
図１に示されるエンジン１において、一番〜四番気筒＃１〜＃４（一番気筒＃１のみ図示）の各燃焼室２には、吸気管３を通じて空気が吸入されるとともに燃料噴射弁４から直接燃料が噴射供給される。この空気と燃料とからなる混合気に対し点火プラグ５による点火が行われると、同混合気が燃焼してピストン６が往復移動し、エンジン１の出力軸であるクランクシャフト７が回転する。そして、燃焼後の混合気は排気として各燃焼室２から排気管８に送り出される。
【００３７】
エンジン１の出力調整は、自動車の運転者によるアクセルペダル１６の踏み込み操作を通じて行われる。即ち、アクセルペダル１６が踏み込み操作されると、そのときのアクセル踏み込み量に応じてスロットルバルブ１９の開度が調節され、吸気管３に流入する空気の量が変化する。そして、燃焼室２に吸入される空気の量に対応した量の燃料が燃料噴射弁４から噴射され、こうして燃焼室２に充填される混合気の量が変化してエンジン１の出力が調整されるようになる。
【００３８】
エンジン１において、燃焼室２と吸気管３との間は吸気バルブ９の開閉動作によって連通・遮断され、燃焼室２と排気管８との間は排気バルブ１０の開閉動作によって連通・遮断される。これら吸気バルブ９及び排気バルブ１０は、クランクシャフト７の回転が伝達される吸気カムシャフト１１及び排気カムシャフト１２の回転に伴い開閉動作するようになる。
【００３９】
自動車には、エンジン１の燃料噴射制御及びスロットル開度制御といった運転制御を行う電子制御装置１５が搭載されている。電子制御装置１５は、エンジン１の制御に係る各種演算処理を実行するＣＰＵ、その制御に必要なプログラムやデータの記憶されたＲＯＭ、ＣＰＵの演算結果等が一時記憶されるＲＡＭ、外部との間で入・出力するための入・出力ポート等を備えている。また、電子制御装置１５には、以下に示される各種センサからの検出信号が入力される。
【００４０】
・アクセルペダル１６の踏み込み量（アクセル踏込量）を検出するアクセルポジションセンサ１７。
・吸気管３に流入する空気の流量を検出するエアフローメータ１８。
【００４１】
・吸気管３内におけるスロットルバルブ１９よりも下流側の圧力（吸気管内圧力）を検出する圧力センサ２２。
・吸気管３内の空気の温度を検出する温度センサ２３。
【００４２】
・クランクシャフト７の回転に対応した信号を出力するクランクポジションセンサ２５。
・吸気カムシャフト１１の回転位置に対応した信号を出力するカムポジションセンサ２６。
【００４３】
電子制御装置１５は、各気筒＃１〜＃４の燃焼室２に吸入される空気の量（以下、筒内充填空気量という）を推定し、それら筒内充填空気量に基づき燃料噴射量制御などエンジン１の運転制御を実行する。
【００４４】
こうしたエンジン１の運転制御を適切に行うには各気筒＃１〜＃４の筒内充填空気量を正確に推定する必要があるが、これら筒内充填空気量は吸気管３に流入する空気の量（流入空気量）が関係しているだけでなく、各気筒＃１〜＃４での吸気バルブ９の開閉に起因する吸気管３内での圧力の脈動にも密接に関係している。従って、各気筒＃１〜＃４の筒内充填空気量を正確に推定するためには、それら筒内充填空気量を推定する際に吸気管３内の圧力の脈動を考慮に入れて当該推定を行うことが重要になる。
【００４５】
そこで本実施形態では、吸気管３内の圧力の脈動によって燃焼室２内（筒内）に充填される空気の量である余分空気量を、吸気バルブ９の開弁による吸気管内圧力の降下量に基づき算出する。そして、この余分空気量と上記流入空気量とに基づき筒内充填空気量を推定することで、推定された筒内充填空気量を上記脈動を考慮した正確なものとする。
【００４６】
次に、上記筒内充填空気量の推定手順の概要について、図２及び図３を参照して説明する。
なお、図２は、エンジン１における吸気管部分のモデルである「吸気管モデルＭ１」の基本概念を示している。また、図３（ａ）において、実線ｍｔは吸気管３に流入する空気の流量（単位はｇ／ｓｅｃ）についての時間経過に伴う推移を示しており、実線ｍｃｉは吸気管３からｉ番気筒に流出する空気の流量（単位はｇ／ｓｅｃ）についての時間経過に伴う推移を示している。この図において、面積Ｂは上記流入空気量（単位はｇ）に相当し、面積Ａは余分空気量（単位はｇ）に相当する。更に、図３（ｂ）の実線は吸気管内圧力の時間経過に伴う変化を示している。
【００４７】
まず、図２の吸気管モデルＭ１について考える。エンジン１の吸気管部分について質量保存則を適用すると、吸気管内圧力Ｐｍと、吸気管３に流入する空気の流量（以下、流入流量ｍｔという）と、吸気管３からｉ番気筒に流出する空気の流量（以下、流出流量ｍｃｉという）とには、以下の式（１）の関係が成り立つ。
【００４８】
【数１】
2005282563.tif 000003
ここで、Ｔｍは吸気管内温度、Ｖｍは吸気管容積、Ｒａは気体定数を空気の平均分子数で除算した値である。この式（１）を積分することによって、時刻ｔからΔｔ秒間における吸気管内圧力Ｐｍの変化量ΔＰｍを表す以下の式（２）が得られるようになる。
【００４９】
【数２】
2005282563.tif 000004
この式（２）から、吸気管内圧力Ｐｍについては、吸気管３への流入流量ｍｔが流出流量ｍｃｉよりも多ければ上昇し、小さければ下降し、等しければ一定となることが分かる。更に、Δｔ秒間における吸気管内圧力Ｐｍの変化量ΔＰｍは、吸気管３内の空気量の変化分に相当することが分かる。なお、図３（ａ）に示されるように、流出流量ｍｃｉについては吸気バルブ９の開閉に依存して間欠的且つ大きく変動するのに対し、流入流量ｍｔについては吸気管３内の空気が流入空気の緩衝の役割を果たすことから緩やかな変化を示す。このため、流出流量ｍｃｉと流入流量ｍｔとの大小関係は反転を繰り返す。これは、上記式（２）の右辺の括弧内の値が一定周期で正負反転を繰り返すこと、即ち吸気管内圧力Ｐｍが一定周期で上昇・下降を繰り返すことを意味し、吸気管内圧力Ｐｍの脈動を表す。
【００５０】
ここで、各気筒＃１〜＃４の吸気バルブ９の開弁期間は重複しないと仮定する。この場合、ｉ番気筒への空気の吸入に際して、吸気管内圧力Ｐｍが最大値Ｐｍmax をとるのは、ｉ番気筒の流出流量ｍｃｉが増大する過程で、吸気管内圧力Ｐｍの時間微分値が「０」となるとき（ｄＰｍ／ｄｔ＝０）、即ち流出流量ｍｃｉが流入流量ｍｔと等しくなるとき（ｍｃｉ＝ｍｔ）である。なお、このときの時刻を最大時刻Ｔmax とする。一方、ｉ番気筒への空気の吸入に際して、吸気管内圧力Ｐｍが最小値Ｐｍmin をとるのは、流出流量ｍｃｉが減少する過程で、吸気管内圧力Ｐｍの時間微分値が「０」となるとき、即ち流出流量ｍｃｉが流入流量ｍｔと等しくなるときである。なお、このときの時刻を最小時刻Ｔmin とする。
【００５１】
ｉ番気筒への空気の吸入によって生じる吸気管内圧力Ｐｍの降下量ΔＰｍdwn 、即ち吸気管内圧力Ｐｍの最大値Ｐｍmax と最小値Ｐｍmin との差分は、以下の式（３）のように表すことができる。
【００５２】
【数３】
2005282563.tif 000005
なお、式（３）の積分項は図３（ａ）の面積Ａに相当する。従って、降下量ΔＰｍdwn は面積Ａに比例することが分かる。この式（３）は、各気筒＃１〜＃４における吸気バルブ９の開弁期間の重複はない条件のもとでは、以下の式（４）のように変形することができる。
【００５３】
【数４】
2005282563.tif 000006
この式（４）は筒内充填空気量ＭＣ（単位はｇ）を算出するためのものである。式（４）の右辺の積分項は、図３（ａ）の面積Ｂに相当し、最大時刻Ｔmax から最小時刻Ｔmin までに吸気管３内に流入した空気の量（流入空気量）を表している。また、式（４）の右辺の「ΔＰｍdwn ／（Ｒａ・Ｔｍ／Ｖｍ）」という項については、「Ｖｍ／Ｔｍ」を定数ｋに置き換えて「ｋ・ΔＰｍdwn ／Ｔｍ」と表すことができる。この項は、図３（ａ）の面積Ａに相当し、筒内に吸入される上記流入空気量分の空気に加えて、吸気管内圧力Ｐｍの脈動により更に筒内に充填される空気の量（余分空気量）を表している。この余分空気量は、上記吸気管内圧力Ｐｍの降下量ΔＰｍdwn 等に基づき算出される。従って、式（４）に基づき算出される筒内充填空気量ＭＣは、上記最大時刻Ｔmax から最小時刻Ｔmin までに吸気管３内に流入した空気の量（流入空気量）と上記余分空気量を加算した値となっている。
【００５４】
ただし、厳密に言えば、実際の筒内充填空気量は図３（ａ）の面積Ａ及び面積Ｂに面積Ｃを加えた値に相当する。式（４）によって算出された筒内充填空気量ＭＣは、面積Ｃに相当する空気量を無視した近似値となっている。この面積Ｃに相当する空気量は、各気筒＃１〜＃４における吸気バルブ９の開弁期間の重複がない条件のもとでは無視できるほど微少であるが、当該閉弁期間の重複が生じる場合には無視できないほど大きくなる。このため、各気筒＃１〜＃４における吸気バルブ９の開弁期間が重複する場合も含めて筒内充填空気量ＭＣを正確に求めるべく、式（４）の右辺の積分項を面積Ｂ及び面積Ｃを求める項に置き換えて以下の式（５）とする。
【００５５】
【数５】
2005282563.tif 000007
式（５）の右辺の「ｍｔ（Δｔdwn ＋ｉvo−ｉvc）／２」という項は、図３（ａ）の面積Ｂ及び面積Ｃに対応する部分をほぼ台形と見なし、その台形の面積を面積Ｂと面積Ｃの和として算出するものである。ここで、Δｔdwn は最大時刻Ｔmax から最小時刻Ｔmin までの時間である。また、ｉvoはｉ番気筒での吸気バルブ９の開弁時刻であり、ｉvcは当該ｉ番気筒での吸気バルブ９の閉弁時刻である。この式（５）によって算出された筒内充填空気量ＭＣは、上記式（４）を用いて算出した場合に比べてより正確な値となる。
【００５６】
次に、式（５）による筒内充填空気量ＭＣ（余分空気量）の算出に用いられる降下量ΔＰｍdwn について、図４のタイムチャートを参照して詳しく説明する。
図４（ａ）において、実線は、圧力センサ２２の出力波形であって吸気管３内の実際の圧力の推移を示しており、各気筒＃１〜＃４の吸気バルブ９の閉弁期間に対応して上昇するとともに開弁期間に対応して下降するものである。電子制御装置１５は、圧力センサ２２からの出力をＡ／Ｄ変換により所定のサンプリング周期Ｘ１で量子化し、その量子化後の値を吸気管内圧力Ｐｍとする。この吸気管内圧力Ｐｍは、吸気管３内の実際の圧力に対応して推移することから、吸気バルブ９の開弁中に降下するようになる。そして、吸気バルブ９の開弁中の吸気管内圧力Ｐｍの最大値が上述した最大値Ｐｍmax とされ、吸気バルブ９の開弁中の吸気管内圧力Ｐｍの最小値が上述した最小値Ｐｍmin とされる。更に、それら最大値Ｐｍmax と最小値Ｐｍmin との差が吸気バルブ９の開弁による吸気管内圧力Ｐｍの降下量ΔＰｍdwn とされる。
【００５７】
ただし、吸気管内圧力Ｐｍは、量子化誤差やノイズによって精度が悪化するため、それに伴い余分空気量算出用の降下量ΔＰｍdwn の精度も悪化する。このため、同一気筒の吸気バルブ９の開閉サイクル複数回分の降下量ΔＰｍdwn を平均し、その平均値を余分空気量算出用の降下量ΔＰｍdwn として採用することで、上述した精度悪化という不具合の抑制が図られる。なお、ここでの吸気バルブ９の一回の開閉サイクルとは、吸気バルブ９が開弁状態から閉弁状態となった時点から次回に開弁状態から閉弁状態になるまでの期間を表している。言い換えれば、図４の吸気管３内の実際の圧力（実線）が最小値になった時点から次回の最小値になる時点までが吸気バルブ９の一回の開閉サイクルということになる。
【００５８】
ところで、吸気管３内の実際の圧力の脈動周期は、吸気バルブ９の開閉態様の変化等によって変化する。例えば、エンジン低回転時に吸気バルブ９の開閉周期が長くなると、上記圧力の脈動周期は長くなる。また、エンジン１が運転状態に応じて稼働気筒数を可変とすることの可能なものであって、稼働気筒数低減時に休止気筒の吸気バルブ９の開閉駆動を停止する場合には、その開閉駆動の停止された吸気バルブ９の分だけ上記圧力の脈動周期が長くなる。こうした脈動周期がＡ／Ｄ変換のサンプリング周期Ｘ１に対し十分に長い場合には、吸気バルブ９の一回の開閉サイクル中におけるサンプリング回数が多くなるため、吸気管内圧力Ｐｍの最大値Ｐｍmax 及び最小値Ｐｍmin が実際の値と大きく食い違うことはない。
【００５９】
しかし、エンジン高回転時や全気筒稼働時など上記脈動周期が短くなってサンプリング周期Ｘ１に近づくと、図４（ｂ）に示されるように吸気バルブ９の一回の開閉サイクル中におけるサンプリング回数が少なくなるため、吸気管内圧力Ｐｍの最大値Ｐｍmax 及び最小値Ｐｍmin が実際の値と大きく食い違うようになる。その結果、それら最大値Ｐｍmax 及び最小値Ｐｍmin の差である降下量ΔＰｍdwn の必要精度確保が困難になる。
【００６０】
このため、サンプリング周期Ｘ１に対し上記脈動周期の長さが所定レベル未満になる旨判断されると、即ち上記脈動周期の長さが降下量ΔＰｍdwn の必要精度確保が困難になるほど近づくと判断されると、上記降下量ΔＰｍdwn の平均値である余分空気量算出用の降下量ΔＰｍdwn の必要精度を確保するための対策が講じられる。こうした対策として、余分空気量算出用の降下量ΔＰｍdwn を求める際、通常よりも高精度の平均化法によって同降下量ΔＰｍdwn を求めることが行われ、これにより当該降下量ΔＰｍdwn の必要精度が確保される。
【００６１】
次に、筒内充填空気量ＭＣの具体的な推定手順について、筒内充填空気量推定ルーチンを示す図５のフローチャートを参照して説明する。この筒内充填空気量推定ルーチンは、電子制御装置１５を通じて例えば所定時間毎及び各気筒毎に時間割り込みにて実行される。なお、筒内充填空気量算出ルーチンの実行間隔（上記所定時間）については、対応する気筒での吸気バルブ９の開弁期間に対して十分に短い時間とされている。
【００６２】
筒内充填空気量推定ルーチンにおいては、まず吸気バルブ９の閉弁から開弁への切り換え直後であるか否かが判断され（Ｓ１０１）、ここで肯定判定であれば吸気バルブ９の開弁時刻ｉvoが取得される（Ｓ１０２）。続いて、吸気バルブ９の開弁から閉弁への切り換え直後であるか否かが判断され（Ｓ１０３）、ここで肯定判定であれば吸気バルブ９の閉弁時刻ｉvcが取得される（Ｓ１０４）。更に、吸気バルブ９の開弁中における流入流量ｍｔの平均値ｍｔave 、及び吸気管内温度Ｔｍの平均値Ｔｍave が算出される（Ｓ１０５）。
【００６３】
続くステップＳ１０６〜Ｓ１０９の処理は、サンプリング周期Ｘ１に対する吸気管３内の圧力の脈動周期の長さに応じて、余分空気量算出用の降下量ΔＰｍdwn を通常精度の平均化法、或いはそれよりも高精度の平均化法によって求めるためのものである。
【００６４】
この一連の処理において、ステップＳ１０６，Ｓ１０７では、サンプリング周期Ｘ１に対する上記脈動周期の長さが上記所定レベル未満であるか否かが判断される。同判断は、クランクポジションセンサ２５からの検出信号に基づき求められるエンジン回転速度が所定値α未満であるか否か（Ｓ１０６）、及び気筒休止が行われているか否か（Ｓ１０７）に基づいて行われることとなる。
【００６５】
即ち、エンジン回転速度が所定値α未満であり、且つ気筒休止が行われていれば（Ｓ１０６、Ｓ１０７で共にＹＥＳ）、上記脈動周期がサンプリング周期Ｘ１から離れた状態であり、サンプリング周期Ｘ１に対する脈動周期の長さが上記所定レベル以上である旨判断され、ステップＳ１０８に進む。なお、ここでの所定レベルとは、サンプリング周期Ｘ１に対する上記脈動周期の長さが、余分空気量算出用の降下量ΔＰｍdwn の必要精度確保を困難とする長さになるレベルである。また、上記所定値αは、こうしたレベルであるか否かの判断を的確に行うことの可能な値に設定されている。
【００６６】
ステップＳ１０８では、余分空気量算出用の降下量ΔＰｍdwn を通常精度の平均化法で求めるとともに、吸気管内圧力Ｐｍの降下量ΔＰｍdwn 分の降下時間ΔＴdwn を算出するための通常精度算出処理が行われる。その後、式（５）に基づき筒内充填空気量ＭＣが算出される（Ｓ１１０）。ここでは、開弁時刻ｉvo 、閉弁時刻ｉvc 、降下量ΔＰｍdwn 、及び降下時間ΔＴdwn が式（５）に代入されるとともに、ステップＳ１０５で算出された平均値ｍｔave 及び平均値Ｔｍave が各々式（５）の流入流量ｍｔ及び吸気管内温度Ｔｍとして代入される。これにより筒内充填空気量ＭＣが算出されることとなる。
【００６７】
一方、ステップＳ１０６，Ｓ１０７の処理で、エンジン回転速度が所定値α未満でないか、或いは気筒休止が行われていなければ（Ｓ１０６、Ｓ１０７のいずれかでＮＯ）、上記脈動周期がサンプリング周期Ｘ１に近づいた状態であり、サンプリング周期Ｘ１に対する脈動周期の長さが上記所定レベル未満である旨判断され、ステップＳ１０９に進む。
【００６８】
ステップＳ１０９では、余分空気量算出用の降下量ΔＰｍdwn を上記通常精度の平均化法よりも高精度の平均化法で求めるとともに、吸気管内圧力Ｐｍの降下量ΔＰｍdwn 分の降下時間ΔＴdwn を算出するための高精度算出処理が行われる。その後、上記ステップＳ１１０の処理が実行される。
【００６９】
次に、筒内充填空気量推定ルーチンにおけるステップＳ１０８の通常精度算出処理について、通常精度算出ルーチンを示す図６のフローチャートを参照して詳しく説明する。同ルーチンは、電子制御装置１５を通じて、筒内充填空気量推定ルーチン（図５）のステップＳ１０８に進む毎に実行される。
【００７０】
通常精度算出ルーチンにおいては、まず吸気バルブ９が開弁中であるか否かが判断される（Ｓ２０１）。ここで肯定判定であれば、吸気管内圧力Ｐｍの最大値Ｐｍmax 及び最大値Ｐｍmax が得られたときの最大時刻Ｔmax の更新が行われるとともに（Ｓ２０２）、吸気管内圧力Ｐｍの最小値Ｐｍmin 及び最小値Ｐｍmin が得られたときの最小時刻Ｔmin の更新が行われる（Ｓ２０３）。即ち、吸気管内圧力Ｐｍが吸気バルブ９の開弁から現在までで最大の値になると、その値が新たな最大値Ｐｍmax として記憶されるとともに、そのときの時刻が新たな最大時刻Ｔmax として記憶される。一方、吸気管内圧力Ｐｍが吸気バルブ９の開弁後から現在までで最小の値になると、その値が新たな最小値Ｐｍmin として記憶されるとともに、そのときの時刻が新たな最小時刻Ｔmin として記憶される。
【００７１】
一方、吸気バルブ９の開弁から閉弁への切り換え直後には（Ｓ２０４：ＹＥＳ）、吸気バルブ９の開弁中における吸気管内圧力Ｐｍの降下量ΔＰｍdwn が、上記記憶された最大値Ｐｍmax と最小値Ｐｍmin との差をとることによって算出され、ＲＡＭに記憶される（Ｓ２０５）。続いて、ステップＳ２０６では、吸気バルブ９の開閉サイクル複数回分、例えばｎ回分の降下量ΔＰｍdwn の平均値ΔＰｍdwnaveが、余分空気量算出用の降下量ΔＰｍdwn として用いるために算出される。具体的には、今回の開閉サイクルよりもｎ回前の開閉サイクルから今回の開閉サイクルまでの開閉サイクル複数回分の降下量ΔＰｍdwn がＲＡＭから読み出され、それら降下量ΔＰｍdwn の和を「ｎ」で除算することにより、上記平均値ΔＰｍdwnaveが算出されることとなる。更に、ステップＳ２０７では、今回の開閉サイクルでの吸気管内圧力Ｐｍの降下量ΔＰｍdwn 分の降下時間ΔＴdwn が、ステップＳ２０２，Ｓ２０３で記憶された最小時刻Ｔmin と最大時刻Ｔmax との差をとることによって算出される。
【００７２】
以上のように、余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）を通常精度の平均化法によって求める場合、Ａ／Ｄ変換のサンプリング周期Ｘ１に対し吸気管３内の圧力の脈動周期が近くなると、上述したように一回の開閉サイクルで求められる降下量ΔＰｍdwn の精度が悪化し、上記平均値ΔＰｍdwnaveの必要精度確保が困難になる。このため、サンプリング周期Ｘ１に対する脈動周期の長さが上記所定レベル未満である場合には、筒内充填空気量算出ルーチンにおけるステップＳ１０９（図５）の高精度算出処理を通じて、上記平均値ΔＰｍdwnaveがより高精度の平均化法によって求められる。
【００７３】
次に、上記ステップ１０９の高精度算出処理について、高精度算出ルーチンを示す図７のフローチャートを参照して詳しく説明する。同ルーチンは、電子制御装置１５を通じて、筒内充填空気量推定ルーチンのステップＳ１０９に進む毎に実行される。
【００７４】
高精度算出ルーチンにおいては、まず吸気バルブ９の各開閉サイクルについて、サンプリング周期Ｘ１毎の吸気管内圧力Ｐｍ、及びそのときの時刻ＴがＲＡＭに記憶される（Ｓ３０１）。仮に、吸気管内圧力Ｐｍの開閉サイクル複数回分に亘る推移が図８に示すようになるとすると、サンプリング周期Ｘ１毎のタイミングＴ１，Ｔ２，Ｔ３・・・・・・での吸気管内圧力Ｐｍ及び時刻Ｔが上記ステップＳ３０１の処理によってＲＡＭに記憶されることとなる。
【００７５】
続いて、吸気バルブ９の開弁から閉弁への切り換え直後に（Ｓ３０２：ＹＥＳ）、開閉サイクル複数回（例えばｎ回）分の吸気管内圧力Ｐｍ及び時刻Ｔに関するデータから、吸気バルブ９の一回分の開閉サイクル中の圧力センサ２２からの出力（圧力）の平均波形が求められる（Ｓ３０３）。ちなみに、一回の開閉サイクルから得られる吸気管内圧力Ｐｍ及び時刻Ｔのデータに基づき当該圧力Ｐｍの推移を求め、その推移の開閉サイクル複数回分を互いに重ねると、図９に示されるようになる。なお、図９は、例えば一回目〜三回目の開閉サイクルにおける吸気管内圧力Ｐｍの推移を、それぞれ実線、破線、及び一点鎖線で示している。これら吸気管内圧力Ｐｍ及び時刻Ｔについての開閉サイクル複数回分のデータから、吸気バルブ９の一回分の開閉サイクル中の圧力センサ２２からの出力（圧力）の平均波形として、例えば図１０に示されるような波形が得られるようになる。  その後、ステップＳ３０４において、上記平均波形から吸気バルブ９の開弁期間における吸気管内圧力Ｐｍの最大値Ｐｍmax 及びそのときの時刻Ｔである最大時刻Ｔmax が算出される。更に、ステップＳ３０５において、上記平均波形から吸気バルブ９の開弁期間における吸気管内圧力Ｐｍの最小値Ｐｍmin 及びそのときの時刻Ｔである最小時刻Ｔmin が算出される。そして、上記最大値Ｐｍmax と最小値Ｐｍmin との差が吸気バルブ９の開閉サイクル複数回分の降下量ΔＰｍdwn の平均値ΔＰｍdwnaveとして算出される（Ｓ３０６）。この平均値ΔＰｍdwnaveは余分空気量算出用の降下量ΔＰｍdwn として精度の高い値となり、式（５）に基づく筒内充填空気量ＭＣの算出に用いられる。また、上記最小時刻Ｔmin と最大時刻Ｔmax との差が吸気管内圧力Ｐｍの降下量ΔＰｍdwn （平均値ΔＰｍdwnave）分の降下時間ΔＴdwn として算出される（Ｓ３０７）。
【００７６】
以上詳述した本実施形態によれば、以下に示す効果が得られるようになる。
（１）各気筒＃１〜＃４の筒内充填空気量ＭＣを推定する際には、吸気管内圧力Ｐｍの脈動によって筒内に充填される空気量（余分空気量）も考慮されるため、それら推定される筒内充填空気量ＭＣを上記脈動をも考慮した正確な値とすることができる。
【００７７】
（２）余分空気量は吸気バルブ９の開弁による吸気管内圧力Ｐｍの降下量ΔＰｍdwn に基づき求められる。また、吸気管内圧力Ｐｍは、圧力センサ２２からの出力をＡ／Ｄ変換によりサンプリング周期Ｘ１で量子化して得られる値であり、量子化誤差やノイズによる精度悪化は避けられないものである。しかし、余分空気量の算出に用いられる降下量ΔＰｍdwn としては、吸気バルブ９の開閉サイクル複数回分の上記降下量ΔＰｍdwn を平均化したもの（平均値ΔＰｍdwnave）が採用される。この平均値ΔＰｍdwnaveを余分空気量算出用の降下量ΔＰｍdwn として用いることで、降下量ΔＰｍdwn の上記量子化誤差やノイズによる精度悪化を抑制することができる。
【００７８】
（３）サンプリング周期Ｘ１に対する吸気管３内の圧力の脈動周期の長さが所定レベル未満である旨判断されると、高精度の平均化法により余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）が求められる。この高精度の平均化法では、吸気バルブ９の開閉サイクル複数回分の吸気管内圧力Ｐｍ及び時刻Ｔに関するデータから、吸気バルブ９の一回分の開閉サイクル中の圧力センサ２２からの出力（圧力）の平均波形が求められ、その波形における最大値Ｐｍmax と最小値Ｐｍmin との差を上記平均値ΔＰｍdwnaveとしている。このように平均波形から平均値ΔＰｍdwnaveを求めることで、その平均値ΔＰｍdwnaveを余分空気量算出用の降下量ΔＰｍdwn として精度の高いものとすることができる。
【００７９】
ただし、上記平均波形を求めるには、吸気管内圧力Ｐｍ及び時刻Ｔに関する多数のデータ（開閉サイクル複数回分）を用いた複雑な計算が必要になる。このため、電子制御装置１５におけるＣＰＵの計算負荷が増大するとともに、上記多数のデータを記憶するのにＲＡＭのメモリ容量を大幅に消費することは避けられない。
【００８０】
このため、サンプリング周期Ｘ１に対する吸気管３内の圧力の脈動周期の長さが所定レベル以上である旨判断される場合には、通常精度の平均化法により余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）が算出される。この通常精度の平均化法では、吸気バルブ９の一回の開閉サイクル毎に、同バルブ９の開閉弁中における吸気管内圧力Ｐｍの最大値Ｐｍmax と最小値Ｐｍmin との差を吸気管内圧力Ｐｍの降下量ΔＰｍdwn として求める。そして、この降下量ΔＰｍdwn について、吸気バルブ９の開閉サイクル複数回分の平均をとって降下量ΔＰｍdwn の平均値ΔＰｍdwnaveを求め、この平均値ΔＰｍdwnaveを余分空気量算出用の降下量ΔＰｍdwn としている。この場合、余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）を算出するのに、上述した高精度の平均化法のような平均波形を求めることは行われず、一回の開閉サイクル毎の降下量ΔＰｍdwn の計算及びその計算結果の記憶、並びに記憶された開閉サイクル複数回分の降下量ΔＰｍdwn の平均といった処理が行われるだけである。このため、電子制御装置１５におけるＣＰＵの計算負荷の増大やＲＡＭのメモリ容量の消費は抑制されるようになる。
【００８１】
以上のように、サンプリング周期Ｘ１に対する圧力の脈動周期の長さが所定レベル未満であるか否かに基づき、上述した高精度の平均化法と通常精度の平均化法とを切り換えて実施することで、余分空気量算出用の降下量ΔＰｍdwn の必要精度を確保しつつ、計算負荷の増大やメモリ容量の消費を最小限にとどめることができる。
【００８２】
［第２実施形態］
次に、本発明の第２実施形態を図１１及び図１２に従って説明する。
この実施形態では、サンプリング周期Ｘ１に対する吸気管３の圧力の脈動周期の長さに応じて、余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）を得るための開閉サイクル回数を可変としている。そして、サンプリング周期Ｘ１に対する脈動周期の長さが所定レベル未満である旨判断されたときには、平均値ΔＰｍdwnaveを得るための上記開閉サイクル回数を多くして当該平均値ΔＰｍdwnaveの必要精度を確保するようにしている。
【００８３】
図１１及び図１２は、本実施形態の筒内充填空気量推定ルーチンを示すフローチャートである。この筒内充填空気量推定ルーチンにおいては、第１実施形態の筒内充填空気量推定ルーチン（図５）におけるステップＳ１０８，Ｓ１０９に相当する処理（ステップＳ４０８〜Ｓ４１０）のみが同ルーチンと異なっている。  図１１及び図１２の筒内充填空気量推定ルーチンにおけるステップＳ４０１〜Ｓ４０５（図１１）では、開弁時刻ｉvo及び閉弁時刻ｉvcが取得されるとともに、流入流量ｍｔの平均値ｍｔave 及び吸気管内温度Ｔｍの平均値Ｔｍave が算出される。続いて図１２のステップＳ４０６，Ｓ４０７では、サンプリング周期Ｘ１に対し吸気管３内の圧力の脈動周期の長さが所定レベル未満であるか否かが判断される。
【００８４】
ここで、脈動周期の長さが所定レベル以上である旨判断されると（Ｓ４０６、Ｓ４０７で共にＹＥＳ）、平均値ΔＰｍdwnave を得るための開閉サイクル回数が、例えば第１実
施形態における通常精度の平均化法と同じｎ回に設定される（Ｓ４０８）。また、脈動周期の長さが所定レベル未満である旨判断されると（Ｓ４０６とＳ４０７とのいずれかでＮＯ）、上記開閉サイクル回数がｎ回よりも多い回数、例えば（ｎ＋δ）回に設定される（Ｓ４０９）。そして、上記のように設定された開閉サイクルの回数分の降下量ΔＰｍdwn の平均値ΔＰｍdwnaveが求められるとともに、今回の開閉サイクルでの降下時間ΔＴdwn が算出される（Ｓ４１０）。こうした平均値ΔＰｍdwnave及び降下時間ΔＴdwn の算出は、例えば第１実施形態の通常精度算出処理と同じ手順で行われる。なお、この処理によって算出された平均値ΔＰｍdwnaveは余分空気量算出用の降下量ΔＰｍdwn として用いられる。
【００８５】
そして、上記開弁時刻ｉvo、閉弁時刻ｉvc、流入流量ｍｔの平均値ｍｔave 、吸気管内温度Ｔｍの平均値Ｔｍave 、降下量ΔＰｍdwn （平均値ΔＰｍdwnave）、及び降下時間ΔＴdwn に基づき、式（５）を用いて筒内充填空気量ＭＣが算出されることとなる（Ｓ４１１）。
【００８６】
本実施形態によれば、第１実施形態の（１）及び（２）の効果に加え、以下に示す効果が得られるようになる。
（４）余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）を得るための開閉サイクル回数は、サンプリング周期Ｘ１に対する吸気管３の圧力の脈動周期の長さに応じて可変とされる。具体的には、上記開閉サイクル回数については、脈動周期の長さが所定レベル未満であると判断されたときには、上記脈動周期の長さが所定レベル以上であると判断されたときの回数（ｎ回）よりも多い回数（（ｎ＋δ）回）とされる。上記開閉サイクル回数を多くすれば、平均値ΔＰｍdwnaveにおける余分空気量算出用の降下量ΔＰｍdwn としての必要精度を確保し易くはなるが、平均値ΔＰｍdwnaveを算出する際の計算負荷の増大や、その算出に必要とされるデータを記憶することによるメモリ容量の消費の増大は避けられない。しかし、上記のように開閉サイクル回数を脈動周期の長さが所定レベル未満になったときのみ多くすることで、余分空気量算出用の降下量ΔＰｍdwn の必要精度を確保しつつ、計算負荷の増大やメモリ容量の消費を最小限にとどめることができる。
【００８７】
［第３実施形態］
次に、本発明の第３実施形態を図１３〜図１７に従って説明する。
この実施形態では、サンプリング周期Ｘ１に対する吸気管３内の圧力の脈動周期の長さに応じて当該サンプリング周期を変更することにより、余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）の必要精度を確保するようにしている。
【００８８】
図１３及び図１４は、本実施形態の筒内充填空気量推定ルーチンを示すフローチャートである。この筒内充填空気量推定ルーチンにおいては、第１実施形態の筒内充填空気量推定ルーチン（図５）におけるステップＳ１０８，Ｓ１０９に相当する処理（ステップＳ５０８〜Ｓ５１２）のみが同ルーチンと異なっている。  図１３及び図１４の筒内充填空気量推定ルーチンにおけるステップＳ５０１〜Ｓ５０５（図１３）では、開弁時刻ｉvo及び閉弁時刻ｉvcが取得されるとともに、流入流量ｍｔの平均値ｍｔave 及び吸気管内温度Ｔｍの平均値Ｔｍave が算出される。続いて図１４のステップＳ５０６，Ｓ５０７では、サンプリング周期Ｘ１に対し吸気管３内の圧力の脈動周期の長さが所定レベル未満であるか否かが判断される。
【００８９】
ここで、脈動周期の長さが所定レベル以上である旨判断されると（Ｓ５０６、Ｓ５０７で共にＹＥＳ）、サンプリング周期が例えば第１実施形態と同じ値（Ｘ１）に設定される（Ｓ５０８）。また、脈動周期の長さが所定レベル未満である旨判断されると（Ｓ５０６とＳ５０７とのいずれかでＮＯ）、サンプリング周期が第１実施形態の値よりも短い値、例えば（Ｘ１−β）に設定される（Ｓ５０９）。このようにサンプリング周期を短くすることで、脈動周期の長さが所定レベル未満であるとき、吸気バルブ９の一回の開閉サイク
ル中におけるサンプリング回数が少なくなるのを抑制することができる。従って、そのサンプリング回数の少なさに起因して余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）の必要精度確保が困難になるのを抑制することができる。
【００９０】
ところで、圧力センサ２２からの出力を所定のサンプリング周期で量子化することで得られる吸気管内圧力Ｐｍの推移は、吸気バルブ９の開閉サイクル毎の圧力センサ２２の出力波形の違いに応じて異なる状態（非同期状態）となるのが普通である。しかし、サンプリング周期に対する脈動周期の長さによっては、開閉サイクル毎における吸気管内圧力Ｐｍの推移が、開閉サイクル毎の圧力センサ２２からの出力波形の違いに関係なく、例えば図１５に示されるように複数回（図では三回）の開閉サイクルで同期した状態となることがある。言い換えれば、いかなるサンプリング周期においても、圧力センサ２２の出力波形の開閉サイクル毎の違いに関係なく、吸気管内圧力Ｐｍの推移が複数の開閉サイクルで同期した状態となる脈動周期が存在するということである。このように吸気管内圧力Ｐｍの推移が吸気バルブ９の複数の開閉サイクルで同期した状態になると、降下量ΔＰｍdwn を開閉サイクル複数回分を平均して平均値ΔＰｍdwnaveを求め、その平均値ΔＰｍdwnaveを余分空気量算出用の降下量ΔＰｍdwn として用いても、その降下量ΔＰｍdwn の必要精度確保が困難になる。
【００９１】
ここで、吸気管内圧力Ｐｍの推移が複数の開閉サイクルで同期した状態になる場合と、非同期の状態になる場合とで、余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）にどのような違いが生じるかについて、図１６及び図１７を併せ参照して説明する。なお、これらの図は、平均値ΔＰｍdwnaveを得るための開閉サイクル回数を増大させるにつれて、平均値ΔＰｍdwnaveがどのように変化してゆくかを示している。そして、図１６は吸気管内圧力Ｐｍの推移が複数の開閉サイクルで同期した状態になる場合に対応しており、図１７は吸気管内圧力Ｐｍの推移が複数の開閉サイクルで非同期の状態になる場合に対応している。
【００９２】
吸気管内圧力Ｐｍの推移が複数の開閉サイクルで非同期の状態になる場合（図１７）、平均値ΔＰｍdwnaveを得るための開閉サイクル回数を増大させるにつれて、当該平均値ΔＰｍdwnaveが所定の値に収束してゆくことになる。これに対し、吸気管内圧力Ｐｍの推移が複数の開閉サイクルで同期した状態になる場合には（図１６）、平均値ΔＰｍdwnaveを得るための開閉サイクル回数を増大させたとしても、当該平均値ΔＰｍdwnaveが変動して所定の値に収束しなくなる。このことは、上記開閉サイクル回数を増大させても、平均値ΔＰｍdwnaveの精度はよくならないことを意味している。従って、その平均値ΔＰｍdwnaveを余分空気量算出用の降下量ΔＰｍdwn として採用しても、当該降下量ΔＰｍdwn の必要精度を確保することが困難になる。
【００９３】
筒内充填空気量推定ルーチンのステップＳ５１０（図１４）では、上述したように可変設定された現在のサンプリング周期に対する脈動周期の長さが上記のような状態、即ち吸気管内圧力Ｐｍの推移が複数の開閉サイクルで同期した状態となる値であるか否かが判断される。なお、同判断は、例えばエンジン回転速度が、上述した同期状態を生じさせる脈動周期となる所定範囲内の値であるか否かに基づいて行われる。この所定範囲は、現在のサンプリング周期及び気筒休止状態に基づき、上記判断を行うのに適した範囲となるよう可変設定されるものである。
【００９４】
上記ステップＳ５１０で肯定判定がなされると、サンプリング周期が例えば所定値γだけ変更される（Ｓ５１１）。これにより、吸気管内圧力Ｐｍの推移が複数の開閉サイクル毎で同期した状態となることは抑制され、余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）の必要精度確保が可能になる。また、このときの所定値γ分のサンプリング周期の変更は、サンプリング周期を短くする方向と長くする方向とのどちらの方向について行ってもよいが、一回の開閉サイクル中のサンプリング回数を極力多くとるという観点ではサンプリング周期を短くする方向について行うことが好ましい。上記所定値γについては、例えば、吸気管内圧力Ｐｍの推移が複数の開閉サイクル毎で同期した状態から脱するのに必要なサンプリング周期の必要最低の変更量に相当する値に設定される。
【００９５】
ステップＳ５１２では、吸気バルブ９の開閉サイクル複数回（例えばｎ回）分の降下量ΔＰｍdwn の平均値ΔＰｍdwnaveが求められるとともに、今回の開閉サイクルでの降下時間ΔＴdwn が算出される。こうした平均値ΔＰｍdwnave及び降下時間ΔＴdwn の算出は、例えば第１実施形態の通常精度算出処理と同じ手順で行われる。なお、この処理によって算出された平均値ΔＰｍdwnave は余分空気量算出用の降下量ΔＰｍdwn として用いられ
る。
【００９６】
そして、上記開弁時刻ｉvo、閉弁時刻ｉvc、流入流量ｍｔの平均値ｍｔave 、吸気管内温度Ｔｍの平均値Ｔｍave 、降下量ΔＰｍdwn （平均値ΔＰｍdwnave）、及び降下時間ΔＴdwn に基づき、式（５）を用いて筒内充填空気量ＭＣが算出されることとなる（Ｓ５１３）。
【００９７】
本実施形態によれば、第１実施形態の（１）及び（２）の効果に加え、以下に示す効果が得られるようになる。
（５）サンプリング周期Ｘ１に対する吸気管３内の圧力の脈動周期の長さに応じて当該サンプリング周期が変更される。具体的には、サンプリング周期については、脈動周期の長さが所定レベル未満であると判断されたときには、上記脈動周期の長さが所定レベル以上であると判断されたときの値（Ｘ１）よりも短い値（Ｘ１−β）に設定される。このようにサンプリング周期を短くすることで、脈動周期がサンプリング周期に過度に近い状態とならないようにされ、一回の開閉サイクル中のサンプリング回数を多くとることができる。このため、一回の開閉サイクル中のサンプリング回数の少なさに起因して、余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）の必要精度を確保しにくくなることは抑制される。ただし、一回の開閉サイクル中のサンプリング回数を多くすると、平均値ΔＰｍdwnaveを計算するためのＣＰＵの計算負荷が増大することは避けられない。このことを考慮し、サンプリング周期が短くされて上記サンプリング回数が多くされるのは、サンプリング周期に対する脈動周期の長さが所定レベル未満であると判断されたときだけとしている。これにより、平均値ΔＰｍdwnaveの必要精度を確保しつつ、上記計算負荷の増大を最小限にとどめることができるようになる。
【００９８】
（６）上記のように可変設定されるサンプリング周期に対し脈動周期の長さが、吸気管内圧力Ｐｍの推移が複数の開閉サイクルで同期した状態となる値である旨判断されると、サンプリング周期が所定値γ分だけ異なる値に変更される。これにより、吸気管内圧力Ｐｍの推移が複数の開閉サイクル毎で同期した状態とならないようにされ、その状態となることによって余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）の必要精度確保が困難になるのを抑制することができる。
【００９９】
［第４実施形態］
次に、本発明の第４実施形態を図１８及び図１９に従って説明する。
この実施形態は、吸気管内圧力Ｐｍの推移が複数の開閉サイクル毎で同期した状態となるとき、第３実施形態のようにサンプリング周期を所定値γだけ変更する代わりに、今回の開閉サイクルでの吸気管内圧力Ｐｍの降下量ΔＰｍdwn の算出を行わないようにしたものである。なお、このときの余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）の算出は、ＲＡＭに記憶された直近の開閉サイクルｎ回分、即ち今回の開閉サイクルを除いた直近の開閉サイクルｎ回分の降下量ΔＰｍdwn のデータに基づき行われる。
【０１００】
図１８は、本実施形態の筒内充填空気量推定ルーチンを示すフローチャートである。この筒内充填空気量推定ルーチンにおいては、第３実施形態の筒内充填空気量推定ルーチンにおけるステップＳ５０６〜Ｓ５１１（図１４）に相当する処理（ステップＳ６０６〜Ｓ６０８）のみが同ルーチンと異なっている。
【０１０１】
図１８の筒内充填空気量推定ルーチンにおけるステップＳ６０１〜Ｓ６０５では、開弁時刻ｉvo及び閉弁時刻ｉvcが取得されるとともに、流入流量ｍｔの平均値ｍｔave 及び吸気管内温度Ｔｍの平均値Ｔｍave が算出される。続いて、ステップＳ６０６において、サンプリング周期Ｘ１に対する脈動周期の長さが、吸気管内圧力Ｐｍの推移が複数の開閉サイクルで同期した状態となる値であるか否かが判断される。ここで肯定判定がなされると、上述した同期状態であるか否かを判断するためのフラグＦが「１（同期状態）」に設定される（Ｓ６０７）。一方、否定判定であれば、フラグＦが「０（非同期状態）」に設定される（Ｓ６０８）。
【０１０２】
続いてステップＳ６０９では、余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）を求めるとともに、今回の開閉サイクルでの降下時間ΔＴdwn を算出する算出処理が実行される。そして、上記開弁時刻ｉvo、閉弁時刻ｉvc、流入流量ｍｔの平均値ｍｔave 、吸気管内温度Ｔｍの平均値Ｔｍave 、降下量ΔＰｍdwn （平均値ΔＰｍdwnave）、及び降下時間ΔＴdwn に基づき、式（５）を用いて筒内充填空気量ＭＣが算出されることとなる（Ｓ６１０）。
【０１０３】
次に、上記ステップＳ６０９の算出処理について、算出処理ルーチンを示す図１９のフローチャートを参照して詳しく説明する。同ルーチンは、電子制御装置１５を通じて、筒内充填空気量推定ルーチン（図１８）のステップＳ６０９に進む毎に実行される。
【０１０４】
この算出処理ルーチンは、第１実施形態における通常精度算出ルーチン（図６）において、ステップＳ２０５，Ｓ２０６に相当する処理（Ｓ７０５〜Ｓ７０７）のみが同ルーチンと異なっている。
【０１０５】
図１９の算出処理ルーチンにおいては、吸気バルブ９の開弁中（ステップＳ７０１：ＹＥＳ）、最大値Ｐｍmax 及び最大時刻Ｔmax の更新が行われるとともに（Ｓ７０２）、最小値Ｐｍmin 及び最小時刻Ｔmin の更新が行われる（Ｓ７０３）。また、吸気バルブ９の開弁から閉弁への切り換え直後には（Ｓ７０４：ＹＥＳ）、フラグＦが「０（非同期状態）」であるか否かが判断される（Ｓ７０５）。
【０１０６】
ステップＳ７０５で肯定判定であれば、吸気バルブ９の開弁中における吸気管内圧力Ｐｍの降下量ΔＰｍdwn が、上記記憶された最大値Ｐｍmax と最小値Ｐｍmin との差をとることによって算出され、ＲＡＭに記憶される（Ｓ７０６）。続いて、今回の吸気バルブ９の開閉サイクルを含めて直近の開閉サイクル複数回分、例えばｎ回分の降下量ΔＰｍdwn から平均値ΔＰｍdwnaveが、余分空気量算出用の降下量ΔＰｍdwn として用いるために算出される（Ｓ７０７）。
【０１０７】
一方、ステップＳ７０５で否定判定であって、フラグＦが「１（同期状態）」であれば、ステップＳ７０６の処理がスキップされ、今回の開閉サイクルでの吸気管内圧力Ｐｍの降下量ΔＰｍdwn の算出、及びＲＡＭへの記憶が行われなくなる。この場合、続くステップＳ７０７の処理では、今回の吸気バルブ９の開閉サイクルを除いた直近の開閉サイクル複数回（ｎ回）分の降下量ΔＰｍdwn から平均値ΔＰｍdwnaveが算出される。
【０１０８】
その後、ステップＳ７０８では、今回の開閉サイクルでの吸気管内圧力Ｐｍの降下量ΔＰｍdwn 分の降下時間ΔＴdwn が、ステップＳ７０２，Ｓ７０３で記憶された最小時刻Ｔmin と最大時刻Ｔmax との差をとることによって算出される。
【０１０９】
本実施形態によれば、第１実施形態の（１）及び（２）の効果に加え、以下に示す効果が得られるようになる。
（７）吸気管内圧力Ｐｍの推移が複数の開閉サイクルで同期した状態となる場合、余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）の必要精度確保が困難になる。このため、サンプリング周期に対し脈動周期の長さが、上記のような状態となる値である旨判断されると、そのときの開閉サイクルでの吸気管内圧力Ｐｍの降下量ΔＰｍdwn の算出が行われなくなる。更に、今回の開閉サイクルを除いた直近の開閉サイクル複数回分の吸気管内圧力Ｐｍの降下量ΔＰｍdwn から平均値ΔＰｍdwnaveが求められるようになる。同平均値ΔＰｍdwnaveについては、吸気管内圧力Ｐｍの推移が複数の開閉サイクルで同期した状態となった開閉サイクルでの吸気管内圧力Ｐｍの降下量ΔＰｍdwn が用いられることなく算出される。従って、上記平均値ΔＰｍdwnaveを余分空気量算出用の降下量ΔＰｍdwn として採用することで、その降下量ΔＰｍdwn の必要精度確保が困難になるのを抑制することができる。
【０１１０】
なお、上記各実施形態は、例えば以下のように変更することもできる。
・第２〜第４実施形態において、平均値ΔＰｍdwnaveを第１実施形態に記載した高精度の平均化法を用いて求めてもよい。
【０１１１】
・第３実施形態の所定値γとして、より大きい値を採用してもよい。
・気筒休止が行われないエンジンに本発明を適用してもよい。この場合、筒内充填空気量算出ルーチンにおける気筒休止が行われているか否かの判断（Ｓ１０７、Ｓ４０７、Ｓ５０７）については省略される。
【０１１２】
・第１実施形態での高精度な平均化法については、同実施形態に記載された方法に限定されるものではなく、通常精度の平均化法よりも高い精度で平均値ΔＰｍdwnaveを求めることの可能な他の方法を採用することもできる。
【０１１３】
・第３実施形態において、吸気管内圧力Ｐｍの推移が複数の開閉サイクル毎で同期した状態になるか否かの判断に係る処理（Ｓ５１０）、及び同期状態となったときに所定値γ分のサンプリング周期の変更を行う処理（Ｓ５１１）を省略してもよい。この場合でも、上記（５）に準ずる効果は得られるようになる。
【０１１４】
・第３実施形態において、サンプリング周期Ｘ１に対する脈動周期の長さが所定レベル未満であるか否かの判断に係る処理（Ｓ５０６、Ｓ５０７）、及び所定レベル未満である旨の判断に基づきサンプリング周期を短くする処理（Ｓ５０９）を省略してもよい。この場合でも、上記（６）に準ずる効果は得られるようになる。
【図面の簡単な説明】
【０１１５】
【図１】第１実施形態の筒内充填空気量推定装置が適用される直列四気筒エンジン全体を示す略図。
【図２】吸気管モデルの基本概念を示す図。
【図３】（ａ）及び（ｂ）は、時間経過に対する流入流量、流出流量、及び吸気管内圧力の変化を示すタイムチャート。
【図４】（ａ）及び（ｂ）は、圧力センサからの出力（圧力）の推移、及び当該出力をＡ／Ｄ変換によりサンプリング周期Ｘ１で量子化して得られる吸気管内圧力Ｐｍの推移を示すタイムチャート。
【図５】第１実施形態での筒内充填空気量の算出手順を示すフローチャート。
【図６】第１実施形態において、余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）を通常精度で算出する手順、及び吸気管内圧力Ｐｍの降下量ΔＰｍdwn 分の降下時間ΔＴdwn を算出する手順を示すフローチャート。
【図７】第１実施形態において、余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）を高精度で算出する手順、及び吸気管内圧力Ｐｍの降下量ΔＰｍdwn 分の降下時間ΔＴdwn を算出する手順を示すフローチャート。
【図８】吸気バルブの開閉サイクル複数回分の吸気管内圧力Ｐｍの推移を模式的に示すタイムチャート。
【図９】開閉サイクル一回分の吸気管内圧力Ｐｍの推移を、開閉サイクル複数回（三回）分を重ねた状態を模式的に示したタイムチャート。
【図１０】開閉サイクル複数回分の吸気管内圧力Ｐｍ及び時刻Ｔに関するデータから得られる開閉サイクル一回分の圧力センサ２２からの出力（圧力）の平均波形を示す波形図。
【図１１】第２実施形態における余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）を算出する手順、及び吸気管内圧力Ｐｍの降下量ΔＰｍdwn 分の降下時間ΔＴdwn を算出する手順を示すフローチャート。
【図１２】第２実施形態における余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）を算出する手順、及び吸気管内圧力Ｐｍの降下量ΔＰｍdwn 分の降下時間ΔＴdwn を算出する手順を示すフローチャート。
【図１３】第３実施形態における余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）を算出する手順、及び吸気管内圧力Ｐｍの降下量ΔＰｍdwn 分の降下時間ΔＴdwn を算出する手順を示すフローチャート。
【図１４】第３実施形態における余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）を算出する手順、及び吸気管内圧力Ｐｍの降下量ΔＰｍdwn 分の降下時間ΔＴdwn を算出する手順を示すフローチャート。
【図１５】吸気管内圧力Ｐｍの推移が複数の開閉サイクル毎で同期した状態を模式的示すタイムチャート。
【図１６】吸気管内圧力Ｐｍの推移が複数の開閉サイクル毎で非同期となった状態のもとで、平均値ΔＰｍdwnaveを得るための開閉サイクル回数を増大させるにつれて、平均値ΔＰｍdwnaveがどのように変化してゆくかを示すグラフ。
【図１７】吸気管内圧力Ｐｍの推移が複数の開閉サイクル毎で同期した状態のもとで、平均値ΔＰｍdwnaveを得るための開閉サイクル回数を増大させるにつれて、平均値ΔＰｍdwnaveがどのように変化してゆくかを示すグラフ。
【図１８】第４実施形態での筒内充填空気量の算出手順を示すフローチャート。
【図１９】第４実施形態における余分空気量算出用の降下量ΔＰｍdwn （平均値ΔＰｍdwnave）を算出する手順、及び吸気管内圧力Ｐｍの降下量ΔＰｍdwn 分の降下時間ΔＴdwn を算出する手順を示すフローチャート。
【符号の説明】
【０１１６】
１…エンジン、２…燃焼室、３…吸気管、４…燃料噴射弁、５…点火プラグ、６…ピストン、７…クランクシャフト、８…排気管、９…吸気バルブ、１０…排気バルブ、１１…吸気カムシャフト、１２…排気カムシャフト、１５…電子制御装置（降下量算出手段、余分空気量算出手段、筒内充填空気量推定手段）、１６…アクセルペダル、１９…スロットルバルブ、１７…アクセルポジションセンサ、１８…エアフローメータ、２２…圧力センサ、２３…温度センサ、２５…クランクポジションセンサ、２６…カムポジションセンサ。
【図１】
2005282563.tif 000008
【図２】
2005282563.tif 000009
【図３】
2005282563.tif 000010
【図４】
2005282563.tif 000011
【図５】
2005282563.tif 000012
【図６】
2005282563.tif 000013
【図７】
2005282563.tif 000014
【図８】
2005282563.tif 000015
【図９】
2005282563.tif 000016
【図１０】
2005282563.tif 000017
【図１１】
2005282563.tif 000018
【図１２】
2005282563.tif 000019
【図１３】
2005282563.tif 000020
【図１４】
2005282563.tif 000021
【図１５】
2005282563.tif 000022
【図１６】
2005282563.tif 000023
【図１７】
2005282563.tif 000024
【図１８】
2005282563.tif 000025
【図１９】
2005282563.tif 000026