(57)【特許請求の範囲】
【請求項１】
燃焼室内に導入される吸入空気の量を吸気通路のスロットル弁により調節する内燃機関の吸入空気量予測装置において、
前記吸気通路内の吸入空気に関する流体モデルと同吸入空気の現時点での温度情報とに基づいて前記吸気通路の前記スロットル弁よりも下流側の部分に存在する吸入空気の現時点から所定期間後の温度を推定し該推定結果を反映しつつ現時点から所定期間後に前記下流側部分から前記燃焼室内に導入される吸入空気の量を予測する予測手段を備える
ことを特徴とする内燃機関の吸入空気量予測装置。
【請求項２】
請求項１に記載の内燃機関の吸入空気量予測装置において、
前記予測手段は、前記吸気通路内の吸入空気に関する流体モデルと同吸入空気の現時点での温度情報とに基づいて、前記吸気通路の前記スロットル弁よりも下流側の部分に存在する吸入空気の現時点から所定期間後の温度に加えさらに同吸入空気の現時点から所定期間後の圧力を推定し、この温度及び圧力の推定結果を反映しつつ現時点から所定期間後に前記下流側部分から前記燃焼室内に導入される吸入空気の量を予測する
ことを特徴とする内燃機関の吸入空気量予測装置。
【請求項３】
請求項１または２に記載の内燃機関の吸入空気量予測装置において、
前記予測手段は、前記下流側部分の吸入空気に関する熱エネルギ保存則並びに質量保存則に基づき構築される流体モデルと、前記下流側部分に流入する吸入空気の現時点での量及び温度、並びに前記下流側部分の吸入空気の現時点での圧力及び温度とに基づいて前記下流側部分の吸入空気の現時点から所定期間後の圧力及び温度を推定し、該推定される圧力及び温度に基づいて前記吸入空気量を予測するものである
ことを特徴とする内燃機関の吸入空気量予測装置。
【請求項４】
請求項３に記載の内燃機関の吸入空気量予測装置において、
前記予測手段は、前記下流側部分に流入する吸入空気の現時点での量をｍt(0)及び温度をＴha、前記下流側部分の吸入空気の現時点での圧力をＰm(0)及び温度をＴm(0)、前記下流側部分から燃焼室に導入される吸入空気の量をｍc(i)、空気の比熱比をκ、気体定数をＲ、前記下流側部分の容積をＶとしたとき、前記各保存則から得られる離散式（１），（２）、並びに前記吸入空気量ｍc(i)と前記下流側部分の吸入空気の圧力Ｐm(i)及び温度Ｔm(i)との関係式（３）
【数１】
0004075233.tif 000032
に基づく演算を所定回数ｎだけ繰り返す処理を実行することにより、前記下流側部分の吸入空気の現時点から所定期間tfwd後の圧力Ｐm(n)及び温度Ｔm(n)を推定するものである
ことを特徴とする内燃機関の吸入空気量予測装置。
【請求項５】
請求項４に記載の内燃機関の吸入空気量予測装置において、
前記予測手段は、前記一連の繰返処理を所定の割込周期毎に実行するものであって、今回の割込タイミングにおいて、次回の割込タイミングでの繰返処理に用いられる前記下流側部分の吸入空気の現時点での圧力Ｐm(0)及び温度Ｔm(0)の少なくとも一方を算出するものである
ことを特徴とする内燃機関の吸入空気量予測装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
この発明は、機関燃焼室内に導入される吸入空気の量を予測する内燃機関の吸入空気量予測装置に関する。
【０００２】
【従来の技術】
吸気通路を通じて内燃機関の燃焼室に導入される空気の量（吸入空気量）は、機関回転数とともに機関運転状態を最も反映するものとして、燃料噴射制御や点火時期制御といった内燃機関の制御量を決定する際に用いられており、通常、これは吸気通路に設けられたエアフローメータや吸気圧センサからの検出信号に基づいて求められる。
【０００３】
ところで、吸入空気量に基づいて内燃機関の制御量を算出するに際し、エアフローメータ等の検出信号から求められる現時点での吸入空気量をそのまま用いるようにすると、過渡運転時にあってはその制御量が機関制御に実際に反映される時点とは異なる吸入空気量に基づいて同制御量が決定されるようになる。このように制御量の決定に用いた吸入空気量と、その制御量が機関制御に反映されるときの吸入空気量とが異なっていると、実際の機関運転状態に即した制御を行うことができず、その制御精度を高めるにしても自ずと限界が生じるようになる。
【０００４】
そこで、特開平２−４２１６０号公報に記載されるように、現時点から所定期間後の吸入空気量を予測し、その予測値に基づいて制御量を決定するようにした装置が従来より提案されている。こうした予測を行うことで、制御量が機関制御に反映される際の吸入空気量を予め見越して同制御量を決定することができるようになり、制御精度の向上が図られるようになる。
【０００５】
また、吸入空気量は、機関回転数やスロットル弁の開度の他、吸気温によっても変化する。例えば、吸気温が低くなるほど吸入空気の密度が大きくなるため、吸入空気量も増大するようになる。従って、上記のように吸入空気量を予測する際には、こうした吸気温の影響を考慮することも精度良い予測結果を得る上で重要になる。
【０００６】
この点、上記装置では、吸気通路のスロットル弁よりも上流側の位置での吸気温を検出し、その吸気温（以下、吸気通路の他の位置での吸気温と特に区別する必要がある場合には、これを「上流側吸気温」という）を吸入空気量の予測に反映させるようにしている。
【０００７】
【発明が解決しようとする課題】
しかしながら、上記従来の装置は、吸入空気量の予測に吸気温を反映させているとはいえ、あくまでもそれは、吸気通路内において吸気温の位置的な変化及びその時間的な変化が予測期間中に生じないことを前提条件としてその予測を行うものでしかない。
【０００８】
確かに、スロットル弁の開度が一定に保持される定常運転時にあっては、こうした吸気通路内における吸気温の位置的な変化は殆ど無く、また時間的な吸気温の変化も無視できる範囲にあるため、上記前提条件のもとで吸入空気量の予測をしたとしても、即ち前述した上流側吸気温のみに基づき吸入空気量の予測をしたとしても、その予測結果に大きな誤差が生じることはない。
【０００９】
ところが、過渡運転時においてスロットル弁の開度が急激に変化する場合にあっては、吸気温が吸気通路内の位置に応じて異なるようになり、またその時間的な変化も大きく無視できないものとなる。従って、上記前提条件はこうした吸気通路内における吸入空気の実際の温度状態に沿わないものとなる。
【００１０】
例えば、スロットル弁の開度が増大することにより、同弁の下流側部分に流入する吸入空気の量が急激に増大した場合、同下流側部分の吸入空気が圧縮された状態になり、その温度（以下、「下流側吸気温」という）は前述した上流側吸気温よりも一時的に上昇するようになる。一方、スロットル弁の開度が減少することにより、前記下流側部分に流入する吸入空気の量が急激に減少した場合には逆に、同下流側部分の吸入空気が膨張した状態になるため、下流側吸気温が上流側吸気温よりも一時的に低下するようになる。
【００１１】
このようにスロットル弁の開度が変化する過渡運転時にあっては、一時的にせよ、下流側吸気温が時間的に大きく変動して上流側吸気温と異なる温度になり得る。その結果、吸気通路内において吸気温の位置的及び時間的な変化が生じないことを前提条件として予測を行う従来の装置にあっては、こうした過渡運転時における吸気温のいわば動的な挙動による影響が吸入空気量の予測に反映されることはなく、従ってその予測結果の誤差も無視できないものとなる。
【００１２】
この発明はこうした実情に鑑みてなされたものであり、その目的は吸気通路を通じて燃焼室内に導入される吸入空気の量を、過渡運転時であっても精度良く予測することのできる内燃機関の吸入空気量予測装置を提供することにある。
【００１３】
【課題を解決するための手段】
上記目的を達成するための手段及びその作用効果について以下に記載する。
請求項１に記載した発明は、燃焼室内に導入される吸入空気の量を吸気通路のスロットル弁により調節する内燃機関の吸入空気量予測装置において、前記吸気通路内の吸入空気に関する流体モデルと同吸入空気の現時点での温度情報とに基づいて前記吸気通路の前記スロットル弁よりも下流側の部分に存在する吸入空気の現時点から所定期間後の温度を推定し該推定結果を反映しつつ現時点から所定期間後に前記下流側部分から前記燃焼室内に導入される吸入空気の量を予測する予測手段を備えるようにしている。
請求項２に記載した発明は、請求項１に記載の内燃機関の吸入空気量予測装置において、前記予測手段は、前記吸気通路内の吸入空気に関する流体モデルと同吸入空気の現時点での温度情報とに基づいて、前記吸気通路の前記スロットル弁よりも下流側の部分に存在する吸入空気の現時点から所定期間後の温度に加えさらに同吸入空気の現時点から所定期間後の圧力を推定し、この温度及び圧力の推定結果を反映しつつ現時点から所定期間後に前記下流側部分から前記燃焼室内に導入される吸入空気の量を予測するようにしている。
【００１４】
上記請求項１または２に記載した発明において、予測手段は、吸気通路内の吸入空気に関する流体モデルと同吸入空気の現時点での温度情報（温度や温度変化の態様等、或いはこれらと相関を有する物理量）とに基づいて前記下流側部分の吸入空気の温度を推定し、その推定結果を吸入空気量の予測に反映させる。従って、この下流側部分の温度が現時点から過渡的に変化する場合であっても、その過渡的な変化を把握した上で吸入空気量を予測することができ、極めて精度の良い予測結果を得ることができるようになる。
【００１５】
また、吸気通路を通じて燃焼室に導入される吸入空気の量は、上記下流側部分の吸入空気の圧力が高くなるほど、またその温度が低くなるほど多くなる傾向があるため、これを上記圧力及び温度に基づいて求めることができる。従って、下流側部分の吸入空気の現時点から所定期間後の圧力及び温度を推定することにより、それら推定される圧力及び温度に基づいて現時点から所定期間後に燃焼室に導入される吸入空気の量を予測することができるようになる。
【００１６】
更に、この下流側部分の吸入空気の現時点から所定期間後の圧力及び温度は、これらを、下流側部分の吸入空気についての熱エネルギ保存則並びに質量保存則に基づいて、下流側部分に流入する吸入空気の現時点での量及び温度、並びにその下流側部分の現時点での吸入空気の圧力及び温度といった物理量と関連付けることにより推定することができる。
【００１７】
従って、請求項１または２に記載した発明をより具体化した構成として、請求項３に記載した発明のように、上記吸入空気量予測装置において、
・前記予測手段は、前記下流側部分の吸入空気に関する熱エネルギ保存則並びに質量保存則に基づき構築される流体モデルと、前記下流側部分に流入する吸入空気の現時点での量及び温度、並びに前記下流側部分の吸入空気の現時点での圧力及び温度とに基づいて前記下流側部分の吸入空気の現時点から所定期間後の圧力及び温度を推定し、該推定される圧力及び温度に基づいて前記吸入空気量を予測するものである、
といった構成を採用することができる。
【００１８】
また、請求項４に記載した発明は、請求項３に記載の内燃機関の吸入空気量予測装置において、前記予測手段は、前記下流側部分に流入する吸入空気の現時点での量をｍt(0)及び温度をＴha、前記下流側部分の吸入空気の現時点での圧力をＰm(0)及び温度をＴm(0)、前記下流側部分から燃焼室に導入される吸入空気の量をｍc(i)、空気の比熱比をκ、気体定数をＲ、前記下流側部分の容積をＶとしたとき、前記各保存則から得られる離散式（１），（２）、並びに前記吸入空気量ｍc(i)と前記下流側部分の吸入空気の圧力Ｐm(i)及び温度Ｔm(i)との関係式（３）
【００１９】
【数２】
0004075233.tif 000002
に基づく演算を所定回数ｎだけ繰り返す処理を実行することにより、前記下流側部分の吸入空気の現時点から所定期間tfwd後の圧力Ｐm(n)及び温度Ｔm(n)を推定するものであるとしている。
【００２０】
上記構成によれば、前記流体モデルが上記各式（１），（２）に示されるように離散化されて表現されているため、前記下流側部分の吸入空気の現時点から所定期間後の圧力Ｐm(n)及び温度Ｔm(n)を比較的簡単なアルゴリズムに基づいて推定することができるようになる。
【００２１】
また、こうした演算処理を行うにあたり、請求項５に記載した発明のように、請求項４に記載の内燃機関の吸入空気量予測装置において、
・前記予測手段は、前記一連の繰返処理を所定の割込周期毎に実行するものであって、今回の割込タイミングにおいて、次回の割込タイミングでの繰返処理に用いられる前記下流側部分の吸入空気の現時点での圧力Ｐm(0)及び温度Ｔm(0)の少なくとも一方を算出するものである、
といった構成を採用すれば、上記吸入空気の圧力Ｐm(0)、或いはその温度Ｔm(0)をセンサ等を用いずとも上記繰返処理を通じて逐次求めることができ、吸入空気量予測装置の構成の簡略化を図ることができるようになる。
【００２２】
また、一般に、汎用の温度センサはその応答性が低いため、同センサを用いて吸入空気の温度Ｔm(0)を検出するようにした場合、同温度Ｔm(0)の変化を精度良く検出するのは難しく、その検出結果に誤差が生じるのは避けきれない。
【００２３】
この点、上記構成において、特に吸入空気量の温度Ｔm(0)を算出するようにすれば、構成の簡略化を図ることができるとともに、こうしたセンサの検出誤差に起因する予測精度の低下を回避することができるようになる。
【００２４】
【発明の実施の形態】
［第１の実施形態］
以下、本発明の第１の実施形態について説明する。
【００２５】
本実施形態では、吸入空気量を検出するセンサとしてエアフローメータを備える内燃機関に本発明に係る吸入空気量予測装置を適用するようにしている。はじめに、図１を参照して、この内燃機関１０の概要と同予測装置の構成について説明する。
【００２６】
内燃機関１０の吸気通路２０は、スロットル弁２４が設けられたスロットルボディ２１の内部通路、このスロットルボディ２１の下流側に接続されたサージタンク２２、及び同サージタンク２２と燃焼室２９とを接続する吸気ポート２３等により構成されている。吸気通路２０を流れる吸入空気は、スロットル弁２４によって調量された後、吸気弁２６の開弁時に燃焼室２９内に導入される。
【００２７】
吸入空気量予測装置は、現時点から所定期間後に燃焼室２９内に導入される吸入空気の量を予測するものであり、大きくは、演算装置３０とエアフローメータ４１をはじめとする各種センサ４１〜４５とによって構成されている。
【００２８】
演算装置３０は、演算処理を実行する演算部３１と、こうした演算処理に係るプログラムやその実行に際して用いられる各種関数データが記憶されたメモリ３２とを備えている。この演算装置３０は、例えば、内燃機関１０の各種制御を実行する電子制御装置によって構成される。
【００２９】
センサ４１〜４５には、上記エアフローメータ４１の他、吸気温センサ４２、スロットルセンサ４３、カム角センサ４４、及びクランク角センサ４５等が含まれる。エアフローメータ４１、吸気温センサ４２はいずれも、吸気通路２０においてスロットル弁２４より上流側の部分に設けられており、同弁２４を介してサージタンク２２内に流入する吸入空気の量とその温度（以下、「上流側吸気温Ｔha」という）をそれぞれ検出する。
【００３０】
また、本実施形態では、上記エアフローメータ４１として、熱線からなる検出部４１０を備えた熱線式のものが用いられている。図２は、この検出部４１０の断面構造を示している。同図に示すように、この検出部４１０は、例えば白金等からなる熱線４１２と、同熱線４１２の周囲を被覆するガラス層４１４とによって構成されている。
【００３１】
スロットルセンサ４３は、スロットル弁２４の近傍に設けられており、同弁２４の開度（スロットル開度ＴＡ）を検出する。カム角センサ４４は、吸気弁２６を開閉駆動するカムシャフト（図示略）の近傍に、クランク角センサ４５は、このカムシャフトを回転駆動するクランクシャフト（図示略）の近傍にそれぞれ設けられている。演算装置３０は、これらカム角センサ４４及びクランク角センサ４５の出力信号を取り込み、クランクシャフトの回転速度（機関回転数ＮＥ）及びその回転位相（クランク角ＣＡ）並びに吸気弁２６の開閉時期（バルブタイミングＶＴ）をそれぞれ検出する。尚、本実施形態の内燃機関１０には、上記吸気弁２６の開閉時期を機関回転数ＮＥや機関負荷等の機関運転状態に基づいて変更する機構（図示略）が設けられている。従って、上記バルブタイミングＶＴも機関運転状態に応じて適宜変更されることとなる。
【００３２】
次に、こうした本実施形態の装置による吸入空気量の予測処理について説明する。
はじめにこの予測処理の前提となる原理について説明する。
【００３３】
まず、吸気通路２０から燃焼室２９までの空間を、スロットル弁２４より上流側の部分（以下、「上流側部分」という）２７、同スロットル弁２４より下流側で吸気弁２６よりも上流側の部分（以下、「下流側部分」という）２８、及び燃焼室２９といった３つの部分に分け、これら各部分２７〜２９のうち下流側部分２８に存在する吸入空気についての流体モデルを構築する。
【００３４】
即ち、この下流側部分２８に存在する吸入空気の質量を「Ｍ」、同下流側部分２８から燃焼室２９内に時間当たりに流出する吸入空気の量（以下、「筒内吸入空気量」という）を「ｍc」、上流側部分２７から下流側部分２８に流入する吸入空気の量（以下、「スロットル弁通過量」という）を「ｍt」とすると、上記質量Ｍの時間当たりの変化量ｄ（Ｍ）／ｄｔは質量保存則に基づいて次式（５）のように表すことができる。
【００３５】
【数３】
0004075233.tif 000003
一方、下流側部分２８に存在する吸入空気の温度（以下、「下流側吸気温」という）を「Ｔm」、その等容比熱及び等圧比熱をそれぞれ「Ｃv」、「Ｃp」とすると、同吸入空気の有する熱エネルギは「Ｃv・Ｍ・Ｔm」となる。また、上流側部分２７から下流側部分２８に流入する吸入空気の有する熱エネルギが「Ｃp・ｍt・Ｔha」、下流側部分２８から燃焼室２９に流出する吸入空気の有する熱エネルギが「Ｃp・ｍc・Ｔm」であることから、上記熱エネルギＣv・Ｍ・Ｔmの時間当たりの変化量ｄ（Ｃv・Ｍ・Ｔm）／ｄｔは熱エネルギ保存則に基づいて次式（６）のように表すことができる。
【００３６】
【数４】
0004075233.tif 000004
また、下流側部分２８に存在する吸入空気についての状態方程式は、下流側部分２８の容積を「Ｖ」、同部分２８での吸入空気の圧力（以下、「下流側吸気圧」という）を「Ｐm」、気体定数を「Ｒ」とすると、次の式（７）のようになる。
【００３７】
【数５】
0004075233.tif 000005
ここで、空気の比熱比を「κ」（＝Ｃp／Ｃv）とおくと、上記式（７）を用いて式（５）は以下の式（９）、式（６）は以下の式（８）のように表すことができる。
【００３８】
【数６】
0004075233.tif 000006
更に、これら各式（８），（９）の左辺は、離散間隔を「Δｔ」とすると次式（１０），（１１）のように表すことができる。
【００３９】
【数７】
0004075233.tif 000007
従って、これら各式（８）〜（１１）から以下の離散式（１２），（１３）が導き出される。尚、この式（１２），（１３）において、括弧内の「ｉ」は今回値、「ｉ−１」は前回値をそれぞれ示す。
【００４０】
【数８】
0004075233.tif 000008
一方、筒内吸入空気量ｍcは、次の実験式（１４）を用いて求めることができる。
【００４１】
【数９】
0004075233.tif 000009
上式（１４）において、「ａ」、「ｂ」はいずれも定数であり、下流側部分２８から燃焼室２９までの部分における吸入空気の流れに影響を与える機関制御パラメータ、即ちここでは、機関回転数ＮＥ及びバルブタイミングＶＴに基づいて設定される。これら定数ａ或いは定数ｂと機関回転数ＮＥ及びバルブタイミングＶＴとの関係は予め実験により定められ、演算装置３０のメモリ３２に関数マップとして記憶されている。また、「Ｔo」は、こうした実験により各定数ａ，ｂを設定したときの吸入空気の基準温度（定数）である。
【００４２】
この実験式（１４）を用いると、上記各式（１２），（１３）は、次の各式（１５），（１６）のように変形することができる。
【００４３】
【数１０】
0004075233.tif 000010
また、下流側吸気温Ｔmは、これら各式（１５），（１６）から得られる以下の式（１７）を用いて求めることができる。
【００４４】
【数１１】
0004075233.tif 000011
これら各式（１５）〜（１７）から明らかなように、現時点でのスロットル弁通過量ｍt、下流側吸気圧Ｐm、下流側吸気温Ｔm、及び上流側吸気温Ｔhaに基づいて、現時点から所定期間後、即ちここでは上記離散間隔Δｔ後の下流側吸気圧Ｐm及び下流側吸気温Ｔmを推定することができ、更にこうして推定される下流側吸気圧Ｐm及び下流側吸気温Ｔmを式（１４）に代入することにより、現時点から離散間隔Δｔ後の筒内吸入空気量ｍcを推定することができる。
【００４５】
また、本実施形態の予測処理では、上記スロットル弁通過量ｍtを以下の異なる二つの方法に従って求めるようにしている。以下、このスロットル弁通過量ｍtを求める際の方法について説明する。
【００４６】
まず、エアフローメータ４１の検出信号に基づいてスロットル弁通過量ｍtを求める際の手順について説明する。
エアフローメータ４１は、周知のように、熱線４１２の周囲を吸入空気が通過する際に同熱線４１２から奪われる熱量が吸入空気量に応じて変化するのを利用して同量を検出するものである。このエアフローメータ４１の検出信号と吸入空気量との関係は予め実験によって求められ、演算装置３０のメモリ３２に関数マップとして記憶されている。従って、エアフローメータ４１の検出信号とこの関数マップとに基づいて吸入空気量を求めることができる。
【００４７】
しかしながら、上記ガラス層４１４は、熱線４１２と比較して一般に熱容量が大きいため、その温度Ｔgsは検出部４１０の近傍を通過する吸入空気の量の変化に対して遅れて変化する傾向がある。このため、吸入空気量が変化しない定常運転時にあっては、上記関数マップから求められる吸入空気量は真の吸入空気量と一致するものの、同吸入空気量が変化する過渡運転時にあっては、関数マップから求められる吸入空気量が真の吸入空気量からずれるようになる。
【００４８】
そこで、本実施形態では、こうしたガラス層４１４の温度変化についての応答遅れを予め見越して真の吸入空気量を求めるために、同層４１４に関する熱平衡モデルを構築し、その逆モデルに基づいて吸入空気量を算出するようにしている。尚、以下の説明では、真の吸入空気量をスロットル弁通過量ｍtとし、エアフローメータ４１の検出信号と上記関数マップとに基づいて得られる吸入空気量をマップ演算値ＧＡとして区別する。
【００４９】
まず、熱線４１２の温度を「Ｔht」とすると（図２参照）、熱線４１２からガラス層４１４に伝達される熱量と、同層４１４から吸入空気に伝達される熱量との釣り合いから、ガラス層４１４の温度Ｔgsの変化量ｄ（Ｔgs）／ｄｔは、以下の式（１８）のように表すことができる。
【００５０】
【数１２】
0004075233.tif 000012
ここで、「Ａ」，「Ｂ」，「Ｃ」，「Ｄ」は、熱線４１２の断面積、長さ、及びその抵抗率や、ガラス層４１４と熱線４１２との間の熱伝達率、ガラス層４１４と吸入空気との間の熱伝達率等々に応じて決定される定数である。
【００５１】
また、定常運転時には、ガラス層４１４と、熱線４１２及び吸入空気との間の熱の授受が無くなるため、次式（１９），（２０）に示すように、式（１８）の左辺、即ちガラス層４１４の温度変化量ｄ（Ｔgs）／ｄｔは「０」になり、また、マップ演算値ＧＡとスロットル弁通過量ｍtとは等しくなる。
【００５２】
【数１３】
0004075233.tif 000013
従って、上記各式（１９），（２０）から次の式（２１）が得られる。
【００５３】
【数１４】
0004075233.tif 000014
更に、この式（２１）を用いて式（１８）を変形すると、ガラス層４１４の温度Ｔgsが消去された次式（２２）が得られる。
【００５４】
【数１５】
0004075233.tif 000015
また、この式（２２）の左辺は、離散間隔を「Δｔ」とすると次式（２３）のように表すことができる。
【００５５】
【数１６】
0004075233.tif 000016
従って、上記各式（２２），（２３）から以下の離散式（２４）が導き出される。尚、この式（２４）において、「α」、「β」は上記定数Ａ〜Ｄによって定まる定数であり、また、括弧内の「ｊ」は今回値、「ｊ−１」は前回値をそれぞれ示している。
【００５６】
【数１７】
0004075233.tif 000017
この式（２４）から明らかなように、現時点でのスロットル弁通過量ｍt(j)は、現時点でのマップ演算値ＧＡ(j)と、離散間隔Δｔ前のマップ演算値ＧＡ(j-1)とに基づいて推定することができる。
【００５７】
次に、スロットル開度ＴＡ等の機関運転状態に基づいてスロットル弁通過量ｍtを求める際の手順について説明する。
スロットル弁２４の開口部分を一種のオリフィスと仮定すると、同オリフィスを通過する流体の量（スロットル弁通過量ｍt）は、周知のように次の各式（２５），（２６）によって表される。
【００５８】
【数１８】
0004075233.tif 000018
ここで、「μ」はオリフィス（スロットル弁２４の開口部分）の流量係数、「Ｐa」は大気圧（一定値）である。また、「Ａt」は、スロットル弁２４の開口面積であり、これはスロットル開度ＴＡから一義的に求めることができる。
【００５９】
ここで、スロットル開度ＴＡが一定に保持される定常運転時を想定すると、式（２５）は次式（２７）のようになる。
【００６０】
【数１９】
0004075233.tif 000019
ここで、「ｍtTA」、「ＰmTA」は、定常運転時におけるスロットル弁通過量ｍt、下流側吸気圧Ｐmをそれぞれ示している。これら定常運転時におけるスロットル弁通過量ｍtTA及び下流側吸気圧ＰmTAは、上流側部分２７から燃焼室２９までの部分における吸入空気の流れに影響を与える機関制御パラメータ、即ちここでは、スロットル開度ＴＡ、機関回転数ＮＥ、及びバルブタイミングＶＴに基づいて決定することができる。これら定常運転時のスロットル弁通過量ｍtTA或いは下流側吸気圧ＰmTAと、スロットル開度ＴＡ、機関回転数ＮＥ、及びバルブタイミングＶＴとの関係は、予め実験によって求められ、演算装置３０のメモリ３２に関数マップとして記憶されている。
【００６１】
この式（２７）を用いると、式（２５）は次式（２８）のように表すことができる。
【００６２】
【数２０】
0004075233.tif 000020
更に、定常運転時にあっては、スロットル弁通過量ｍtTAと筒内吸入空気量ｍcとが等しくなり（ｍtTA＝ｍc）、また、下流側吸気温Ｔmと上流側吸気温Ｔhaとが等しくなるため（Ｔm＝Ｔha）、上記式（１４）は次式（２９）のように表すことができる。
【００６３】
【数２１】
0004075233.tif 000021
そして、この式（２９）を用いると、上記式（２８）は更に次式（３０）のように表すことができる。
【００６４】
【数２２】
0004075233.tif 000022
この式（３０）から明らかなように、スロットル弁通過量ｍtは、現時点での上流側吸気温Ｔha及び下流側吸気圧Ｐm、大気圧Ｐa、並びに定常運転時であると仮定したときの下流側吸気圧Ｐmの値ＰmTA（ＴＡ，ＮＥ，ＶＴ）に基づいて推定することができる。
【００６５】
本実施形態では、このように各別の方法に基づいて求められるスロットル弁通過量ｍtと、前述した各式（１４）〜（１７）とに基づいて現時点から所定期間後の筒内吸入空気量ｍcを予測するようにしている。尚、以下の説明では、スロットル開度ＴＡ等の機関運転状態に基づいて求められるスロットル弁通過量ｍtをスロットル弁通過量推定値ｍt1、エアフローメータ４１の検出信号に基づいて求められるスロットル弁通過量ｍtをスロットル弁通過量測定値ｍt2として区別する。
【００６６】
次に、本実施形態の装置による予測処理の詳細について図３及び図４を併せ参照して説明する。
図３及び図４は、筒内吸入空気量ｍcを予測する際の手順の一例を示すフローチャートである。この一連の処理は、演算装置３０により所定時間毎の割込周期（例えば「８msec.」毎）で実行される。
【００６７】
この処理に際しては、まず、吸入空気量を予測する際の予測期間tfwdが設定される（ステップ１００）。
この予測期間tfwdは、前回の燃料噴射での燃料噴射量、換言すれば燃料噴射時間に応じて設定されるものである。従って、ここでは、現時点から前回の燃料噴射に要した時間が経過した後の筒内吸入空気量ｍcが予測されることとなる。因みに、この燃料噴射時間は、この処理とは別の燃料噴射に係る処理を通じて算出されるものである。また、この燃料噴射に係る処理では、本処理を通じて算出される最新の筒内吸入空気量ｍcに基づいて燃料噴射量が算出される。
【００６８】
次に、現時点から予測期間tfwd後の下流側吸気圧Ｐm、下流側吸気温Ｔm、及びそれらの比Ｐm／Ｔmが前記スロットル弁通過量推定値ｍt1に基づいて算出される（ステップ２００）。尚、以下では、スロットル弁通過量推定値ｍt1に基づいて算出される予測期間tfwd後の下流側吸気圧Ｐm、下流側吸気温Ｔm、及びそれらの比Ｐm／Ｔmをそれぞれ、「Ｐm1tfwd」、「Ｔm1tfwd」、「Ｐm1tfwd／Ｔm1tfwd」と、また、その算出途中での下流側吸気圧Ｐm、下流側吸気温Ｔm、及びそれらの比Ｐm／Ｔmを、「Ｐm1(i)」、「Ｔm1(i)」、「Ｐm1／Ｔm1(i)」とそれぞれ表記する。
【００６９】
この下流側吸気圧Ｐm1tfwd及び下流側吸気温Ｔm1tfwdの算出に際しては、上記各式（１５）〜（１７）及び式（３０）に基づく演算を所定回数だけ繰り返す処理が行われる。
【００７０】
即ち、図４に示すように、まず、この処理の繰返回数を表すカウンタ値ｉ（初期値：「０」）がインクリメントされ（ステップ２０２）、次に予測期間tfwdが所定時間Δt1以上であるか否かが判断される（ステップ２０４）。ここで、上記所定時間Δt1は、本ルーチンの割込周期と等しく（例えば「８msec.」）設定されている。そして、予測期間tfwdがこの所定時間Δt1以上であると判断された場合（ステップ２０４：ＹＥＳ）、上記各式（１５），（１６）の離散間隔Δｔが所定時間Δt1と等しく設定される（ステップ２０６）。
【００７１】
次に、式（３０）に基づいてスロットル弁通過量推定値ｍt1(i)が算出され、更にその後、各式（１５）〜（１７）に基づいて、下流側吸気圧Ｐm1(i)、下流側吸気温Ｔm1(i)、及びそれらの比Ｐm1(i)／Ｔm1(i)がそれぞれ算出される（ステップ２０８）。そして、現在の予測期間tfwdから所定時間Δt1が減算され、その減算値（tfwd−Δt1）が新たな予測期間tfwdとして設定される（ステップ２１０）。
【００７２】
次に、この新たな予測期間tfwdが「０」以下であるか否かが判断され（ステップ２１２）、同予測期間tfwdが「０」より長いと判断された場合（ステップ２１２：ＮＯ）には、ステップ２０２以降の処理が再度実行される。この際、ステップ２０４において予測期間tfwdが所定時間Δt1未満であると判断されると（ステップ２０４：ＮＯ）、離散間隔Δｔが予測期間tfwdと等しく設定されて（ステップ２０７）、ステップ２０８以降の処理が実行される。
【００７３】
そして、ステップ２０２〜２１２の処理が繰り返されることにより、予測期間tfwdが「０」以下になると（ステップ２１２：ＹＥＳ）、下流側吸気圧Ｐm1(i)、下流側吸気温Ｔm1(i)、及び比それらのＰm1／Ｔm1(i)が、予測期間tfwd後の下流側吸気圧Ｐm1tfwd、下流側吸気温Ｔm1tfwd、及び比Ｐm1tfwd／Ｔm1tfwdとしてそれぞれ設定される（ステップ２１４）。
【００７４】
例えば、本ルーチンの割込周期が「８msec.」であり、先のステップ１００の処理を通じて予測期間tfwdが「３０msec.」に設定された場合、離散間隔Δｔは上記ステップ２０４，２０６，２０７の処理を通じて、「８msec.」→「８msec.」→「８msec.」→「６msec.」といった態様で順に設定され、ステップ２０２〜ステップ２１２までの処理が「４回」繰り返されることとなる。そして、カウンタ値ｉが「４」であるときの下流側吸気圧Ｐm1（４）、下流側吸気温Ｔm1（４）、及びそれらの比Ｐm1／Ｔm1（４）が、現時点から予測期間tfwd後の下流側吸気圧Ｐm1tfwd、下流側吸気温Ｔm1tfwd、及び比Ｐm1tfwd／Ｔm1tfwdとしてそれぞれ設定されることとなる。
【００７５】
その後、カウンタ値ｉが「０」にリセットされるとともに、離散間隔Δｔが上記所定時間Δt1と等しく設定された後（ステップ２１６）、処理は図３に示すステップ３００に戻される。
【００７６】
ステップ３００では、エアフローメータ４１の検出信号に基づいて現時点でのマップ演算値ＧＡ(j)（以下、「ＧＡ(0)」と表記する）が算出される。次に、この現時点でのマップ演算値ＧＡ(0)と前回の割込タイミングにおいて算出されたマップ演算値ＧＡ(j-1)（以下、「ＧＡ(-1)」と表記する）とに基づき、先の式（２４）を用いて現時点でのスロットル弁通過量測定値ｍt2(0)が算出される（ステップ４００）。
【００７７】
そして、このスロットル弁通過量測定値ｍt2(0)に基づき、先の各式（１５）〜（１７）を用いて離散間隔Δｔ後の下流側吸気圧Ｐm(i)、下流側吸気温Ｔm(i)、及びそれらの比Ｐm／Ｔm(i)がそれぞれ算出される（ステップ５００）。尚、このステップ５００において算出される下流側吸気圧Ｐm(i)、下流側吸気温Ｔm(i)、及びそれらの比Ｐm／Ｔm(i)のうち、今回の割込タイミングで算出される値をそれぞれ「Ｐm2(1)」、「Ｔm2(1)」、「Ｐm2／Ｔm2(1)」と、前回の割込タイミングで算出された値をそれぞれ「Ｐm2(0)」、「Ｔm2(0)」、「Ｐm2／Ｔm2(0)」と表記する。
【００７８】
次に、以下の各式（３１），（３２），（３３）に基づいて現時点から予測期間tfwd後の下流側吸気圧Ｐmtfwd、下流側吸気温Ｔmtfwd、及びそれらの比Ｐmtfwd／Ｔmtfwdが算出されるとともに、上記式（１４）から得られる式（３４）に基づいて予測期間tfwd後における筒内吸入空気量ｍctfwdが算出される（ステップ６００）。
【００７９】
【数２３】
0004075233.tif 000023
こうして筒内吸入空気量ｍctfwdが算出された後、次回の割込タイミングにおける演算処理の初期値がそれぞれ設定される（ステップ７００）。
【００８０】
図５は、こうした初期値の設定態様を示すタイミングチャートである。同図では、時刻ｔ１において機関運転が開始されるとともに、本ルーチンが開始されるものとする。
【００８１】
同図に示すように、２回目以降の割込タイミング（時刻ｔ２，ｔ３）における初期値、即ち下流側吸気圧Ｐm1(0)，Ｐm2(0)、下流側吸気温Ｔm1(0)，Ｔm2(0)、これらの比Ｐm1／Ｔm1(0)，Ｐm2／Ｔm2(0)、並びにマップ演算値ＧＡ(-1)は、前回の割込タイミングにおいて離散間隔Δｔ後の値（Ｐm1（１）等）或いは現時点での値（ＧＡ(0)）として算出された値を用いてそれぞれ以下のように設定される。
【００８２】
【数２４】
0004075233.tif 000024
一方、最初の割込タイミング（時刻ｔ１）における初期値については、前回の割込タイミングでの算出値が存在しないことから、以下のように設定される。
【００８３】
【数２５】
0004075233.tif 000025
このようにして次回の割込タイミングにおける各初期値を設定した後、本ルーチンの処理は一旦終了される。
【００８４】
以上説明したように、本実施形態では、下流側部分２８の吸入空気に関する熱エネルギ保存則並びに質量保存則に基づいて流体モデルを構築し、その流体モデルと、スロットル弁通過量ｍt及び上流側吸気温Ｔha、並びに下流側吸気圧Ｐm及び下流側吸気温Ｔmとに基づいて現時点から予測期間tfwd後の下流側吸気圧Ｐmtfwd及び下流側吸気温Ｔmtfwdを推定し、更にそれらに基づいて現時点から予測期間tfwd後の筒内吸入空気量ｍctfwdを予測するようにしている。
【００８５】
（１）従って、下流側吸気温Ｔmが現時点から過渡的に変化する場合であっても、その過渡的な変化を把握した上で筒内吸入空気量ｍcを予測することができ、極めて精度の良い予測結果を得ることができるようになる。
【００８６】
（２）また、こうした吸入空気の流体モデルを構築する際に、下流側部分２８の吸入空気に関する質量保存則及び熱エネルギ保存則を先の式（１２），（１３）に示されるような離散式によって表現し、これら離散式（１２），（１３）に基づく演算を所定回数繰り返す処理を実行することにより、現時点から予測期間tfwd後の下流側吸気圧Ｐmtfwd及び下流側吸気温Ｔmtfwdを推定するようにしている。従って、これら筒内吸入空気量ｍctfwdの予測に必要な下流側吸気圧Ｐmtfwd及び下流側吸気温Ｔmtfwdを比較的簡単なアルゴリズムに基づいて推定することができるようになる。
【００８７】
（３）更に、上記繰返処理を所定の割込周期で実行する際に、次回の割込タイミングで用いられる現時点での下流側吸気圧Ｐm(0)及び下流側吸気温Ｔm(0)を今回の割込タイミングにおいて算出するようにしている。このため、下流側部分２８内の圧力や温度を検出するためのセンサ等を用いずとも、これら下流側吸気圧Ｐm(0)及び下流側吸気温Ｔm(0)を上記繰返処理を通じて逐次求めることができ、吸入空気量予測装置の構成の簡略化を図ることができるようになる。
【００８８】
（４）特に、現時点での下流側吸気温Ｔm(0)を推定することにより、一般に応答性の低い汎用の温度センサを用いる必要が無くなり、こうしたセンサの検出誤差に起因する予測精度の低下を回避することができるようになる。
【００８９】
また、本実施形態では、スロットル弁通過量ｍtとして、スロットル開度ＴＡ等の機関運転状態から推定されるスロットル弁通過量推定値ｍt1と、エアフローメータ４１の検出信号に基づいて算出されるスロットル弁通過量測定値ｍt２とをそれぞれ求めるようにしている。
【００９０】
ここで、これらスロットル弁通過量推定値ｍt1及びスロットル弁通過量測定値ｍt２を用いて求められる下流側吸気圧Ｐm1(i)，Ｐm2(0)，Ｐm2(1)及び下流側吸気温Ｔm1(i)，Ｔm2(0)，Ｔm2(1)は、実際の吸気圧や吸気温と比較した場合、以下のような傾向があることが本発明者らの実験によって確認されている。
【００９１】
即ち、スロットル弁通過量推定値ｍt1を用いて求められる下流側吸気圧Ｐm1(i)、下流側吸気温Ｔm1(i)、及びそれらの比Ｐm1／Ｔm1(i)に関しては、特に、その時間的変化が吸気圧や吸気温の実際の時間的変化を正確に反映する傾向がある。一方、スロットル弁通過量測定値ｍt２を用いて求められる下流側吸気圧Ｐm2(0)，Ｐm2(1)、下流側吸気温Ｔm2(0)，Ｔm2(1)、及びそれらの比Ｐm2／Ｔm2(0)，Ｐm2／Ｔm2(1)に関しては、定常運転時において吸気圧や吸気温、或いはそれらの比の実際の値とよく一致する傾向がある。
【００９２】
このため、先の各式（３１）〜（３３）を用いて最終的な予測期間tfwd後の下流側吸気圧Ｐmtfwd及び下流側吸気温Ｔmtfwdを算出することにより、過渡運転時にあっては、実際の吸気圧及び吸気温の時間的変化を正確にその算出に反映させることができるとともに、定常運転時、即ち、
【００９３】
【数２６】
0004075233.tif 000026
が満たされるときには、
【００９４】
【数２７】
0004075233.tif 000027
となることから、下流側吸気圧Ｐmtfwd及び下流側吸気温Ｔmtfwdがスロットル弁通過量測定値ｍt２のみに基づいて求められるようになる。
【００９５】
（５）従って、本実施形態によれば、定常運転時であれ過渡運転時であれ、現時点から予測期間tfwd後の下流側吸気圧Ｐmtfwd及び下流側吸気温Ｔmtfwdを正確に推定することができ、ひいては筒内吸入空気量ｍctfwdをより正確に予測することができるようになる。
【００９６】
（６）更に、本実施形態では、エアフローメータ４１のガラス層４１４に関する熱平衡モデルを構築し、その逆モデルに基づいて上記スロットル弁通過量測定値ｍt２を求めるようにしているため、同層４１４の温度変化についての応答遅れを予め見越してこれを求めることができるようになる。従って、より正確なスロットル弁通過量測定値ｍt２を求めることができ、これに基づいて下流側吸気圧Ｐmtfwdや下流側吸気温Ｔmtfwdの推定を更に正確に行うことができるようになる。
【００９７】
［第２の実施形態］
以下、本発明の第２の実施形態について第１の実施形態との相違点を中心に説明する。
【００９８】
本実施形態では、吸入空気量を検出するセンサとして吸気圧センサを備える内燃機関に本発明に係る吸入空気量予測装置を適用するようにしている点が第１の実施形態と相違している。この吸気圧センサ４６は、図６に示すように、サージタンク２２に取り付けられており、その内部の吸入空気の圧力、即ち下流側吸気圧Ｐmを検出する。
【００９９】
次に、本実施形態の装置による予測処理の詳細について説明する。
図７は、筒内吸入空気量ｍcを予測する際の手順の一例を示すフローチャートである。この一連の処理は、演算装置３０により所定時間毎の割込周期（例えば「８msec.」毎）で実行される。尚、同図７に示すステップ１００及びステップ２００では、先の図３及び図４に示すステップ１００及びステップ２００〜２１６と同様の処理が実行されるため、その処理内容の説明は省略する。
【０１００】
これらステップ１００，２００の各処理が実行された後、ステップ３５０では、吸気圧センサ４６の検出信号に基づいて現時点での下流側吸気圧Ｐm(i)（以下、「Ｐm3(0)」と表記する）が算出される。次に、この今回の割込タイミングにおいて検出された下流側吸気圧Ｐm3(0)及び前回の割込タイミングにおいて検出された下流側吸気圧Ｐm(i-1)（以下、「Ｐm3(-1)」と表記する）と、先の式（１５）から得られる次式（３５）とに基づいて、現時点から離散間隔Δｔ前での（即ち前回の割込タイミングでの）スロットル弁通過量ｍt(i-1)（以下、「スロットル弁通過量測定値ｍt3(-1)」という）が算出される（ステップ４５０）。
【０１０１】
【数２８】
0004075233.tif 000028
そして、更にこのスロットル弁通過量測定値ｍt3(-1)と、先の式（１６），（１７）から得られる次の式（３６），（３７）とに基づいて、現時点での下流側吸気温Ｔm(0)（以下、「Ｔm3(0)」と表記する）、並びに下流側吸気圧Ｐm3(0)と下流側吸気温Ｔm3(0)との比Ｐm3／Ｔm3(0)がそれぞれ算出される（ステップ５５０）。
【０１０２】
【数２９】
0004075233.tif 000029
尚、上記ステップ４５０及びステップ５５０において、式（３５）及び（３６）に基づいて下流側吸気温Ｔm3(0)等を算出する場合、上流側吸気温Ｔhaは前回の割込タイミングで読み込まれた値が用いられ、また、各定数ａ，ｂについては前回の割込タイミングでの機関回転数ＮＥ及びバルブタイミングＶＴに基づいて設定された値が用いられる。また、今回が最初の割込タイミングである場合には、上流側吸気温Ｔhaは現時点での値、各定数ａ，ｂは機関停止時の機関回転数ＮＥ（＝「０」）及びバルブタイミングＶＴに基づいて設定される値が用いられる。
【０１０３】
次に、先の各式（３１），（３２）において、「Ｐm2(0)」、「Ｐm2／Ｔm2(0)」をそれぞれ「Ｐm3(0)」、「Ｐm3／Ｔm3(0)」に置換した式と、各式（３３），（３４）とを用いて、現時点から予測期間tfwd後の下流側吸気圧Ｐmtfwd、下流側吸気温Ｔmtfwd、及びそれらの比Ｐmtfwd／Ｔmtfwdが算出されるとともに、予測期間tfwd後における筒内吸入空気量ｍctfwdが算出される（ステップ６５０）。
【０１０４】
こうして筒内吸入空気量ｍctfwdが算出された後、次回の割込タイミングにおける演算処理の初期値がそれぞれ設定される（ステップ７５０）。
図８は、こうした初期値の設定態様を示すタイミングチャートである。同図では、時刻ｔ１において機関運転が開始されるとともに、本ルーチンが開始されるものとする。
【０１０５】
同図に示すように、２回目以降の割込タイミング（時刻ｔ２，ｔ３）における初期値、即ち下流側吸気圧Ｐm1(0)，Ｐm3(-1)、下流側吸気温Ｔm1(0)，Ｔm3(-1)、並びにこれらの比Ｐm1／Ｔm1(0)，Ｐm3／Ｔm3(-1)は、前回の割込タイミングにおいて離散間隔Δｔ後の値（Ｐm1（１）等）或いは現時点での値として求められた値（Ｐm3(0)，Ｔm3(0)等）を用いてそれぞれ以下のように設定される。
【０１０６】
【数３０】
0004075233.tif 000030
一方、最初の割込タイミング（時刻ｔ１）における初期値については、前回の割込タイミングでの算出値或いは検出値が存在しないことから、以下のように設定される。
【０１０７】
【数３１】
0004075233.tif 000031
このようにして次回の割込タイミングにおける各初期値を設定した後、本ルーチンの処理は一旦終了される。
【０１０８】
以上説明した第２の実施形態においても、第１の実施形態に記載した（１）、（２）、及び（４）に関しては略同等の作用効果を得ることができる。
また第１の実施形態に記載した（３）の作用効果については、下流側吸気圧Ｐmを検出するための吸気圧センサ４６が別途必要になるものの、逆にエアフローメータ４１は不要になり、また、少なくとも下流側吸気温Ｔmを検出するセンサを必要としないため、この点では吸入空気量予測装置の構成の簡略化を図ることができるようになる。
【０１０９】
更に、吸気圧センサ４６の検出信号に基づいて求められる下流側吸気圧Ｐm3(0)，Ｐm3（１）、下流側吸気温Ｔm3(0)，Ｔm3（１）、及びそれらの比Ｐm3／Ｔm3(0)，Ｐm3／Ｔm3（１）は、定常運転時において吸気圧や吸気温、或いはそれらの比の実際の値とよく一致する傾向がある。従って、本実施形態によっても第１の実施形態において記載した（５）と同等の作用効果を奏することができる。
【０１１０】
以上説明した各実施形態は、以下のように構成を変更して実施することもできる。
・上記各実施形態では、スロットル弁２４を通過する吸入空気の量として、スロットル開度ＴＡ等の機関運転状態に基づく推定値と、エアフローメータ４１や吸気圧センサ４６等のセンサの検出信号に基づく実測値とをそれぞれ求め、これら各値を用いて筒内吸入空気量ｍcを予測するようにしたが、これら推定値及び実測値の一方に基づいて同予測を行うようにしてもよい。
【０１１１】
・上記各実施形態では、吸気弁２６のバルブタイミングを変更する機構を備えた内燃機関１０を想定して筒内吸入空気量ｍcの予測を行うようにしたが、同機構を備えていない内燃機関においても同様の手順に従って筒内吸入空気量ｍcの予測を行うことができる。
【０１１２】
・また、例えば吸気通路２０の吸入空気の流れを制御する機構、例えば燃焼室２９に発生するスワールの強度を制御するスワール弁等が同吸気通路２０内に設けられる内燃機関において上記予測を行う場合には、上記定数ａ，ｂ及び定常運転時における下流側吸気圧ＰmTA及びスロットル弁通過量ｍtTAを上記機構の制御量（例えばスワール弁の開度等）の関数として設定するようにする。
【０１１３】
・上記各実施形態の予測処理においては、筒内吸入空気量ｍcを求める際に、同空気量ｍcを下流側吸気温Ｔm及び下流側吸気圧Ｐmの関数とする実験式（１４）を用いるようにしたが、同空気量ｍcを求める方法はこれに限定されるものではない。
【０１１４】
・上記各実施形態では、吸入空気の流体モデルとして得られる上記各式（８），（９）を離散化し、繰返演算を実行することにより、下流側吸気圧Ｐm及び下流側吸気温Ｔmを算出するようにしたが、これら各式（８），（９）を解析的に解いた演算式に基づいて、これら下流側吸気圧Ｐm及び下流側吸気温Ｔmを算出するようにしてもよい。
【０１１５】
・上記各実施形態では、燃料噴射量の算出に用いられる筒内吸入空気量ｍcを予測すべく、上記予測期間tfwdを同燃料噴射量の算出に適した値に設定するようにしたが、この予測処理は予測期間tfwdを適宜設定することにより、点火時期制御等、その他の機関制御にも適用することができる。
【０１１６】
・上記各実施形態では、図３や図７に示す予測処理ルーチンの割込周期と、離散間隔Δｔとを等しく設定するようにしたが、例えば、この離散間隔Δｔを次回の割込タイミングでの初期値を決定できる範囲、即ち割込周期の「１／２」倍等、１／ｎ倍（ｎ：整数）となるように変更することもできる。
【図面の簡単な説明】
【図１】第１の実施形態における吸入空気量予測装置の概略構成図。
【図２】エアフローメータの検出部の断面図。
【図３】第１の実施形態の装置による吸入空気量の予測手順を示すフローチャート。
【図４】同じく吸入空気量を予測する際の手順を示すフローチャート。
【図５】第１の実施形態の装置による各割込タイミングでの初期値の設定手順を示すタイミングチャート。
【図６】第２の実施形態における吸入空気量予測装置の概略構成図。
【図７】第２の実施形態の装置による吸入空気量の予測手順を示すフローチャート。
【図８】第２の実施形態の装置による各割込タイミングでの初期値の設定手順を示すタイミングチャート。
【符号の説明】
１０…内燃機関、２０…吸気通路、２１…スロットルボディ、２２…サージタンク、２３…吸気ポート、２４…スロットル弁、２７…上流側部分、２８…下流側部分、２９…燃焼室、２６…吸気弁、３０…演算装置、３１…演算部、３２…メモリ、４１…エアフローメータ、４２…吸気温センサ、４３…スロットルセンサ、４４…カム角センサ、４５…クランク角センサ、４６…吸気圧センサ、４１０…検出部、４１２…熱線。
【図１】
0004075233.tif 000033
【図２】
0004075233.tif 000034
【図３】
0004075233.tif 000035
【図４】
0004075233.tif 000036
【図５】
0004075233.tif 000037
【図６】
0004075233.tif 000038
【図７】
0004075233.tif 000039
【図８】
0004075233.tif 000040