JP2007211746A

JP2007211746A - 内燃機関の制御装置

Info

Publication number: JP2007211746A
Application number: JP2006035377A
Authority: JP
Inventors: Satoshi Tanaka; 聡田中
Original assignee: Toyota Motor Corp
Current assignee: Toyota Motor Corp
Priority date: 2006-02-13
Filing date: 2006-02-13
Publication date: 2007-08-23
Anticipated expiration: 2026-02-13
Also published as: JP4605040B2

Abstract

【課題】内燃機関の制御に用いられるマップ値の学習を効率良く且つ高い精度で行うことができる内燃機関の制御装置を提供する。
【解決手段】機関運転領域毎にマップ値が定められていると共に機関運転状態が特定の機関運転領域にある場合にはその機関運転領域に対応するマップ値を用いて内燃機関を制御する内燃機関の制御装置において、運転パラメータの実測値を検出する実測値検出手段と、機関運転状態が特定の機関運転領域にある場合に特定の機関運転領域に対応するマップ値を用いて運転パラメータの推定値を算出する推定値算出手段と、運転パラメータの実測値と推定値とに基づいて特定の機関運転領域に対応するマップ値を修正するマップ値修正手段とを具備し、マップ値修正手段は、機関運転状態が最長学習時間よりも長く特定の機関運転領域にある場合には最長学習時間が経過するまでの運転パラメータの実測値と推定値に基づいてマップ値を修正する。
【選択図】図２１

Description

本発明は、内燃機関の制御装置に関する。

内燃機関の運転を適切に行うために内燃機関の燃焼室内の吸気ガスの空燃比を目標空燃比に一致させる必要がある。このため多くの内燃機関では内燃機関の運転に関するパラメータ（以下、「運転パラメータ」と称す）の値を検出すると共に検出した運転パラメータの値に基づいて吸気ガスの空燃比が目標空燃比となるような燃料噴射量を算出し、算出された量の燃料を燃料噴射弁から噴射するようにしている。

このような燃料噴射量は、一般的に、実際に検出された運転パラメータ（例えば、エアフロメータの出力値、機関回転数、スロットル開度等）の値や、これら運転パラメータの値に応じて定まるマップ値を用いて数値計算によって算出される。しかしながら、使用による特性の変化や製作公差によるばらつきなどにより斯かる数値計算のみによっては最適な燃料噴射量を正確に算出することができない場合がある。そこで、内燃機関の制御装置の多くは斯かる数値計算に用いられるマップ値等の学習を行い、数値計算によって常に最適な燃料噴射量を算出することができるようにしている。

特許文献１に記載の制御装置では、吸入空気量と機関回転数から定まる基本燃料噴射量を算出すると共に、この基本燃料噴射量に空燃比フィードバック補正係数及び所定のマップ値を乗算することで最適な燃料噴射量を算出することとしている。このマップ値は機関運転領域毎に設定されており、内燃機関の運転に伴って学習せしめられる。空燃比フィードバックのみでは実際の空燃比を目標空燃比に到達させるのに時間がかかるところ、学習されたマップ値を用いることにより実際の空燃比を目標空燃比に迅速に到達させることができるようにしている。

また、特許文献１に記載の制御装置では、学習初期においては、機関運転領域毎にマップ値を定めずに、全ての機関運転領域に対して一つのマップ値となるように学習が行われる。その後、所定時間経過後に複数の機関運転領域に分割すると共に各機関運転領域におけるマップ値の初期値を上述したようにして学習された値とし、その後各機関運転領域毎にマップ値の学習が行われる。

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

ところで、例えば内燃機関が定常運転を行っているような場合、機関運転状態は長期間に亘って同一の機関運転領域内にあることになる。通常、上述したような学習は機関運転状態が特定の機関運転領域から離脱してから、機関運転状態がその特定の機関運転領域内にあったときの運転パラメータの検出値等に基づいて行われる。したがって、長期間に亘って続く定常運転中には学習が行われず、よってその定常運転中には学習結果を反映させることができない。すなわち、内燃機関の制御に用いられるマップ値の学習を効率良く行うことができない。

また、各機関運転領域について一つのマップ値が設定されているが、実際には或る機関運転領域内の各機関運転状態毎に最適なマップ値が存在し、或る機関運転領域に対応するマップ値としてはこの機関運転領域の中心となる機関運転状態について最適な値とするのが好ましい。しかしながら、特許文献１に記載の制御装置では、このように或る機関運転領域の中心となる機関運転状態について最適な値となるようにマップ値の学習を行ってはおらず、よって最適なマップ値を高い精度で算出することができない。

したがって本発明の目的は、内燃機関の制御に用いられるマップ値の学習を効率良く且つ高い精度で行うことができる内燃機関の制御装置を提供することにある。

上記課題を解決するために、第１の発明では、機関運転領域毎にマップ値が定められていると共に、実際の機関運転状態が特定の機関運転領域にある場合にはその機関運転領域に対応するマップ値を用いて内燃機関を制御する、内燃機関の制御装置において、実際の機関運転状態が特定の機関運転領域にある場合における運転パラメータの実測値を検出する実測値検出手段と、実際の機関運転状態が上記特定の機関運転領域にある場合に、該特定の機関運転領域に対応するマップ値を用いて運転パラメータの推定値を算出する推定値算出手段と、実際の機関運転状態が上記特定の機関運転領域にある場合において上記実測値検出手段によって検出された運転パラメータの実測値と上記推定値算出手段によって算出された推定値とに基づいて上記特定の機関運転領域に対応するマップ値を修正するマップ値修正手段とを具備し、上記マップ値修正手段は、実際の機関運転状態が上記特定の機関運転領域に侵入してから所定の最長学習時間よりも長く該特定の機関運転領域にある場合には該最長学習時間が経過するまでの上記運転パラメータの実測値と推定値に基づいてマップ値を修正する。
第１の発明によれば、実際の機関運転状態が最長学習時間よりも長く特定の機関運転領域内にある場合には、最長学習時間までのデータに基づいてマップ値の修正が行われる。このため、内燃機関が長期間に亘って定常運転を行っているような場合においても、実際の機関運転状態がその特定の機関運転領域から離脱しなくてもマップ値の修正を行うことができ、よって効率的にマップ値の修正を行うことができる。

第２の発明では、第１の発明において、上記最長学習時間は機関運転状態に応じて定められる。

上記課題を解決するために、第３の発明では、機関運転領域毎にマップ値が定められていると共に、実際の機関運転状態が特定の機関運転領域にある場合にはその機関運転領域に対応するマップ値を用いて内燃機関を制御する、内燃機関の制御装置において、実際の機関運転状態が特定の機関運転領域にある場合における運転パラメータの実測値を検出する実測値検出手段と、実際の機関運転状態が上記特定の機関運転領域にある場合に、該特定の機関運転領域に対応するマップ値を用いて上記運転パラメータの推定値を算出する推定値算出手段と、実際の機関運転状態が上記特定の機関運転領域にある場合において上記実測値検出手段によって検出された運転パラメータの実測値と上記推定値算出手段によって算出された推定値とに基づいて上記特定の機関運転領域に対応するマップ値を修正するマップ値修正手段とを具備し、上記マップ値修正手段は、実際の機関運転状態が上記特定の機関運転領域に侵入してから所定の最長学習時間よりも長く該特定の機関運転領域にある場合には、実際の機関運転状態が上記特定の機関運転領域にある場合における上記運転パラメータの実測値と推定値とのうち実際の機関運転状態が該特定の機関運転領域の中心領域にある場合における上記運転パラメータの実測値と推定値とを選択し、選択された上記運転パラメータの実測値と推定値とに基づいて当該機関運転領域に対応するマップ値を修正する。
第３の発明によれば、実際の機関運転状態が中心領域にある場合におけるデータに基づいてマップ値の修正が行われる。したがって、或る機関運転領域の中心となる機関運転状態について最適な値となるようにマップ値の修正が行われており、よって最適なマップ値を高い精度で算出することができる。
なお、或る機関運転領域の「中心領域」とは、その機関運転領域に含まれる領域であってその機関運転領域の中心近傍の領域を意味する。

第４の発明では、第１〜第３のいずれか一つの発明において、上記マップ値修正手段は、実際の機関運転状態が上記特定の機関運転領域に侵入してから上記最長学習時間よりも短い最短学習時間を経過することなく該特定の機関運転領域から離脱した場合には、マップ値の修正を行わない。

第５の発明では、第１〜第４のいずれか一つの発明において、上記マップ値が流量係数として用いられるスロットルモデル計算式であってスロットル弁開度と吸気管内圧力とを用いてスロットル通過空気流量の推定値を算出するスロットルモデル計算式によって算出されたスロットル通過空気流量の推定値を用いて内燃機関が制御され、上記実測値検出手段は所定の遅れをもってスロットル通過空気流量に対応する出力値を出力するエアフロメータであり、上記推定値算出手段はスロットル弁開度と吸気管内圧力とを用いて上記スロットルモデル計算式によって算出されるスロットル通過空気流量の推定値を用いてエアフロメータによって出力されるであろう出力値をエアフロメータの出力値の予想値として算出する予想値算出手段であり、上記マップ値修正手段は、エアフロメータの出力値とその予想値とに基づいて上記スロットルモデル計算式に用いられるマップ値を修正する。

第６の発明では、第１〜第４のいずれか一つの発明において、上記マップ値および吸気管内圧力を用いて筒内吸入空気流量の推定値を算出する吸気弁モデル計算式によって算出された筒内吸入空気流量の推定値を用いて内燃機関が制御され、上記実測値検出手段は、実際の吸気管内圧力を検出する吸気管内圧力検出手段であり、上記推定値算出手段は、筒内吸入空気流量を用いて吸気管モデル計算式によって吸気管内圧力の推定値を算出する吸気管内圧力推定値算出手段であり、上記マップ値修正手段は、上記吸気管内圧力推定手段によって算出された吸気管内圧力の推定値と上記吸気管内圧力検出手段によって検出された実際の吸気管内圧力とに基づいて上記吸気弁モデル計算式を修正する。

本発明によれば、内燃機関の制御に用いられるマップ値の学習を効率良く且つ高い精度で行うことができる内燃機関の制御装置が提供される。

以下、図面を参照して本発明の実施形態について説明する。図１に概略的に示した機関本体１は筒内噴射型火花点火式内燃機関を示す。しかしながら、本発明をポート噴射型火花点火式内燃機関等の別の火花点火式内燃機関や圧縮自着火式内燃機関に適用してもよい。

図１に示したように、本発明の実施形態では機関本体１はシリンダブロック２と、シリンダブロック２内で往復動するピストン３と、シリンダブロック２上に固定されたシリンダヘッド４とを具備する。ピストン３とシリンダヘッド４との間には燃焼室５が形成される。シリンダヘッド４には各気筒毎に吸気弁６と、吸気ポート７と、排気弁８と、排気ポート９とが配置される。さらに、図１に示したようにシリンダヘッド４の内壁面の中央部には点火プラグ１０が配置され、シリンダヘッド４内壁面周辺部には燃料噴射弁１１が配置される。また、ピストン３の頂面には燃料噴射弁１１の下方から点火プラグ１０の下方まで延びるキャビティ１２が形成されている。さらに、シリンダヘッド４には、吸気弁６の位相角およびバルブリフト量を連続的に変更可能な吸気弁制御装置１３が設けられている。

各気筒の吸気ポート７は吸気枝管１４を介してサージタンク１５に連結され、サージタンク１５は吸気管１６を介してエアクリーナ１７に連結される。吸気管１６内にはステップモータ１８によって駆動されるスロットル弁１９が配置される。一方、各気筒の排気ポート９は排気管２０に連結され、この排気管２０は排気浄化触媒２１を内蔵したケーシング２２に連結される。なお、以下の説明では、スロットル弁１９から吸気弁６までの吸気枝管１４、サージタンク１５、吸気管１６等の部分を吸気管部分２３と称す。

電子制御ユニット（ＥＣＵ）３１はディジタルコンピュータからなり、双方向性バス３２を介して相互に接続されたＲＡＭ（ランダムアクセスメモリ）３３、ＲＯＭ（リードオンリメモリ）３４、ＣＰＵ（マイクロプロセッサ）３５、入力ポート３６および出力ポート３７を具備する。サージタンク１３には、吸気管部分２３内の空気（吸気ガス）の圧力を検出するための吸気管内圧力センサ４０が設けられており、吸気管内圧力センサ４０は吸気管部分２３内の圧力に比例した出力電圧を発生し、この出力電圧が対応するＡＤ変換器３８を介して入力ポート３６に入力される。

また、スロットル弁１９の上流側の吸気管１６には吸気管１６内を通過して流れる吸入空気流量を検出するエアフロメータ４１が設けられており、本実施形態におけるエアフロメータ４１の構成については後述する。さらに、エアクリーナ１７近傍には吸気温度を検出する吸気温度センサ４２と、大気圧を検出する大気圧センサ４３とが設けられる。スロットル弁１９にはスロットル弁１９の開度を検出するスロットル弁開度センサ４４が設けられており、スロットル弁開度センサ４４はスロットル弁開度に対応する出力信号を発生させる。これらエアフロメータ４１、吸気温度センサ４２、大気圧センサ４３およびスロットル弁開度センサ４４はそれぞれ吸入空気流量（質量流量）、吸気温度（大気温度）、大気圧およびスロットル弁開度に対応する出力信号を発生し、この出力信号が対応するＡＤ変換器３８を介して入力ポート３６に入力される。

また、アクセルペダル４５にはアクセルペダル４５の踏込み量に比例した出力電圧を発生する負荷センサ４６が接続され、負荷センサ４６の出力電圧は対応するＡＤ変換器３８を介して入力ポート３６に入力される。また、クランク角センサ４７は例えばクランクシャフトが１５度回転する毎に出力パルスを発生し、この出力パルスが入力ポート３６に入力される。ＣＰＵ３５ではこのクランク角センサ４５の出力パルスから機関回転数が計算される。

一方、出力ポート３７は対応する駆動回路３９を介して点火プラグ１０、燃料噴射弁１１、吸気弁制御装置１３およびステップモータ１８に接続される。

図２は本実施形態で用いられるエアフロメータ４１の概略斜視図を、図３は図２に示したエアフロメータ４１の熱線計量部の拡大斜視図である。本実施形態では、エアフロメータ４１は、熱線式流量計とされており、図２に示したように吸気管１６内を流れる空気の一部をバイパスさせるバイパス通路と、このバイパス通路にバイパスされた吸入空気の質量流量を計測する熱線計量部４１ａと、計測された質量流量に応じた電圧を出力する信号処理部４１ｂとを有する。熱線計量部４１ａは図３に示したように白金熱線から成る吸気温度計測用抵抗４１ａ１と、この吸気温度計測用抵抗４１ａ１を信号処理部４１ｂに連結して保持するサポート部４１ａ２と、加熱用抵抗（ボビン部）４１ａ３と、この加熱用抵抗４１ａ３を信号処理部４１ｂに連結して保持するサポート部４１ａ４とを備える。信号処理部４１ｂは、吸気温度計測用抵抗４１ａ１と加熱用抵抗４１ａ３とで構成されるブリッジ回路を有し、このブリッジ回路により吸気温度計測用抵抗４１ａ１と加熱用抵抗４１ａ３との温度差を常に一定に維持するように加熱用抵抗４１ａ３に供給する電力を調整すると共に、この供給する電力を電圧に変換して出力するようになっている。エアフロメータ４１の出力電圧Ｖｇとエアフロメータ４１が配置された吸気管１６内を通過する空気の流量（以下、「エアフロ通過空気流量」と称す）との関係は図４に示した通りである。

ところで、内燃機関の制御装置では、内燃機関の燃焼室５において燃焼される混合気の空燃比を目標空燃比にするために、吸気弁６が閉じたときに燃焼室５内に充填されている空気の量（以下、「筒内充填空気量Ｍｃ」と称す）を推定し、推定された筒内充填空気量Ｍｃに基づいて混合気の空燃比が目標空燃比となるように燃料噴射弁１１によって内燃機関の燃焼室５（または吸気通路）に噴射する燃料の量（以下、「燃料噴射量」と称す）を定めている。したがって、内燃機関の燃焼室５において燃焼される混合気の空燃比を正確に目標空燃比とするためには、筒内充填空気量Ｍｃを正確に推定する必要がある。

通常、筒内充填空気量Ｍｃは、エアフロメータ等の多数のセンサ、およびこれらセンサからの出力値を引数とした多数のマップから推定される。ところが、このようにマップを用いて筒内充填空気量Ｍｃを推定する場合、推定される筒内充填空気量Ｍｃの値をより正確なものにするためには、必要なマップの数およびその引数の数が多くなる。このようにマップの数が多くなると、マップを保存するためのＥＣＵのＲＯＭを記憶容量の大きいものにしなければならず、内燃機関の制御装置の製造コストが高くなってしまう。さらに、各マップを作成するにはマップが用いられる内燃機関の形式毎に適合作業を行わなければならないが、この適合作業における測定点はマップの数およびその引数の数に応じて増大するため、マップの数およびその引数の数が多くなると適合作業の工数も増大してしまう。

そこで、マップを用いずに様々なモデルを用いて、数値計算により筒内充填空気量Ｍｃを算出する内燃機関の制御装置が検討されている。このような制御装置では、数値計算を多用することにより必要なマップの数を極力減らすようにしており、これにより適合作業を行う際の工数を大幅に削減しながらも、筒内充填空気量Ｍｃを正確に算出することができる。このような制御装置のうち本願出願人により提案されたものの一つに、図５に示したエアモデルを搭載した制御装置がある。図示したエアモデルは内燃機関に適用される最も単純なモデルであり、以下このエアモデルについて説明する。

エアモデルは、図５に示したようにスロットルモデルＭ１０、吸気管モデルＭ２０、吸気弁モデルＭ３０を備える。スロットルモデルＭ１０には、スロットル弁開度センサ４４によって検出されたスロットル弁１９の開度（スロットル弁開度）θｔと、大気圧センサ４３によって検出された内燃機関周囲の大気圧（または、吸気管１６に吸入される空気の圧力）Ｐａと、吸気温度センサ４２によって検出された内燃機関周囲の大気温度（または、吸気管１６に吸入される空気の温度）Ｔａと、後述する吸気管モデルＭ２０において算出された吸気管部分２３内の圧力（以下、「吸気管内圧力」と称す）Ｐｍとが入力され、これら入力された各パラメータの値を後述するスロットルモデルＭ１０のモデル式に代入することで、単位時間当たりにスロットル弁１９を通過する空気の流量（以下、「スロットル通過空気流量の推定値ｍｔ」と称す）が算出される。スロットルモデルＭ１０において算出されたスロットル通過空気流量の推定値ｍｔは、吸気管モデルＭ２０へ入力される。

吸気管モデルＭ２０には、スロットルモデルＭ１０において算出されたスロットル通過空気流量の推定値ｍｔと、以下で詳述する単位時間当たりに燃焼室５内に流入する空気の流量（以下、「筒内吸入空気流量ｍｃ」と称す。なお、筒内吸入空気流量ｍｃの定義については、吸気弁モデルＭ３０の説明において詳述する）とが入力され、これら入力された各パラメータの値を後述する吸気管モデルＭ２０のモデル式に代入することで、吸気管部分２３内に存在する空気の圧力（以下、「吸気管内圧力Ｐｍ」と称す）と吸気管部分２３内に存在する空気の温度（以下、「吸気管内温度Ｔｍ」と称す）とが算出される。吸気管モデルＭ２０において算出された吸気管内圧力Ｐｍと吸気管内温度Ｔｍは共に吸気弁モデルＭ３０へ入力され、さらに吸気管内圧力ＰｍはスロットルモデルＭ１０にも入力される。

吸気弁モデルＭ３０には、吸気管モデルＭ２０において算出された吸気管内圧力Ｐｍおよび吸気管内温度Ｔｍの他に大気温度Ｔａが入力され、これら入力された各パラメータの値を後述する吸気弁モデルＭ３０のモデル式に代入することで、筒内吸入空気流量ｍｃが算出される。算出された筒内吸入空気流量ｍｃは、筒内充填空気量Ｍｃに変換され、この筒内充填空気量Ｍｃに基づいて燃料噴射弁からの燃料噴射量が決定される。また、吸気弁モデルＭ３０において算出された筒内吸入空気流量ｍｃは吸気管モデルＭ２０に入力される。

図５から分かるように、エアモデルでは或るモデルにおいて算出されたパラメータの値が別のモデルへの入力値として利用されるので、エアモデル全体では、実際に入力される値はスロットル弁開度θｔ、大気圧Ｐａ、および大気温度Ｔａの三つのパラメータのみであり、これら三つのパラメータから筒内充填空気量Ｍｃが算出される。

次に、エアモデルの各モデルＭ１０〜Ｍ３０について説明する。
スロットルモデルＭ１０では、大気圧Ｐａ、大気温度Ｔａ、吸気内圧力Ｐｍ、スロットル弁開度θｔから、下記式（１）に基づいてスロットル通過空気流量の推定値ｍｔが算出される。ここで、式（１）におけるμｔはスロットル弁における流量係数で、スロットル弁開度θｔの関数であり、よって図６に示したようなマップから定まる。また、Ａｔはスロットル弁の開口断面積を示し、スロットル弁開度θｔの関数であり、図７に示したようなマップから定まる。なお、これら流量係数μｔおよび開口断面積Ａｔをまとめたμｔ・Ａｔをスロットル弁開度θｔから一つのマップで求まるようにしてもよい。また、Ｒは気体定数に関する定数であり、実際には気体定数を１mol当たりの気体（空気）の質量Ｍlmolで除算した値である。

また、Φ（Ｐｍ／Ｐａ）は下記式（２）に示した関数であり、この式（２）におけるκは比熱比（一定値とする）である。この関数Φ（Ｐｍ／Ｐａ）は図８に示したようなグラフに表すことができるので、このようなグラフをマップとしてＥＣＵ３１のＲＯＭ３４に保存し、実際には式（２）を用いて計算するのではなくマップからΦ（Ｐｍ／Ｐａ）の値を求めるようにしてもよい。

これらスロットルモデルＭ１０の式（１）および式（２）は、スロットル弁１９上流の気体の圧力を大気圧Ｐａ、スロットル弁１９上流の気体の温度を大気温度Ｔａ、スロットル弁１９の下流の気体の圧力を吸気管内圧力Ｐｍとして、図９に示したようなスロットル弁１９のモデルに対して、質量保存則、エネルギ保存則および運動量保存則を適用し、さらに気体の状態方程式、比熱比の定義式、およびマイヤーの関係式を利用することによって得られる。

吸気管モデルＭ２０では、スロットル通過空気流量の推定値ｍｔ、筒内吸入空気流量ｍｃ、および大気温度Ｔａから、下記式（３）および式（４）に基づいて吸気管内圧力Ｐｍおよび吸気管内温度Ｔｍが算出される。なお、式（３）および式（４）におけるＶｍはスロットル弁１９から吸気弁６までの吸気枝管１４、サージタンク１５、吸気管１６等の部分（吸気管部分２３）の容積に等しい定数である。

ここで、吸気管モデルＭ２０について図１０を参照して説明する。吸気管部分２３の総気体量（総空気量）をＭとすると、総気体量Ｍの時間的変化は、吸気管部分に流入する気体の流量、すなわちスロットル通過空気流量の推定値ｍｔと、吸気管部分から流出する気体の流量、すなわち筒内吸入空気流量ｍｃとの差に等しいため、質量保存則により下記式（５）が得られ、この式（５）および吸気管部分２３における気体の状態方程式（Ｐｍ・Ｖｍ＝Ｍ・Ｒ・Ｔｍ）より、式（３）が得られる。

また、吸気管部分２３の気体のエネルギＭ・Ｃｖ・Ｔｍの時間的変化量は、吸気管部分２３に流入する気体のエネルギと吸気管部分２３から流出する気体のエネルギとの差に等しい。このため、吸気管部分２３に流入する気体の温度を大気温度Ｔａ、吸気管部分２３から流出する気体の温度を吸気管内温度Ｔｍとすると、エネルギ保存則により下記式（６）が得られ、この式（６）および上記気体の状態方程式より、式（４）が得られる。

吸気弁モデルＭ３０では、吸気管内圧力Ｐｍ、吸気管内温度Ｔｍ、および大気温度Ｔａから、下記式（７）に基づいて、筒内吸入空気流量ｍｃが算出される。なお、式（７）におけるａ、ｂは、機関回転数Ｎｅから、さらに吸気弁６の位相角（バルブタイミング）および作用角を変更できる可変動弁機構を備えた内燃機関の場合には吸気弁６の位相角、作用角から定まる値である。

上述した吸気弁モデルＭ３０について図１１を参照して説明する。一般に、吸気弁６が閉じたときに燃焼室５内に充填されている空気の量である筒内充填空気量Ｍｃは、吸気弁６が閉弁するとき（吸気弁閉弁時）に確定し、吸気弁６閉弁時の燃焼室５内の圧力に比例する。また、吸気弁６閉弁時の燃焼室５内の圧力は吸気弁６上流の気体の圧力、すなわち吸気管内圧力Ｐｍと等しいとみなすことができる。したがって、筒内充填空気量Ｍｃは、吸気管内圧力Ｐｍに比例すると近似することができる。

ここで、一定時間（例えば、クランク角７２０°分）当たりに吸気管部分２３から流出する全空気流量を平均化したもの、または一定時間（例えば、クランク角７２０°分）当たりに吸気管部分２３から全ての気筒の燃焼室５に吸入される空気量を上記一定時間で除算したものを筒内吸入空気流量ｍｃ（以下で詳述する）とすると、筒内充填空気量Ｍｃが吸気管内圧力Ｐｍに比例することから、筒内吸入空気流量ｍｃも吸気管内圧力Ｐｍに比例すると考えられる。このことから、理論および経験則に基づいて、上記式（７）が得られる。なお、式（７）における値ａは比例係数であり、値ｂは燃焼室５内に残存していた既燃ガスを表す値（排気弁８閉弁時に燃焼室５内に残る既燃ガス量を後述する時間ΔＴ180°で除算したものと考えられる）である。また、実際の運転では過渡運転時に吸気管内温度Ｔｍが大きく変化する場合があるため、これに対する補正として理論および経験則に基づいて導かれたＴａ／Ｔｍが乗算されている。

ここで、筒内吸入空気流量ｍｃについて、図１２を参照して内燃機関が４気筒である場合について説明する。なお、図１２は横軸がクランクシャフトの回転角度、縦軸が単位時間当たりに吸気管部分２３から燃焼室５に実際に流入する空気の流量である。図１２に示したように、４気筒の内燃機関では、吸気弁６が例えば１番気筒、３番気筒、４番気筒、２番気筒の順に開弁し、各気筒に対応する吸気弁６の開弁量に応じて吸気管部分２３から各気筒の燃焼室５内へ空気が流入する。例えば、吸気管部分２３から各気筒の燃焼室５内に流入する空気の流量の変位は図１２に破線で示した通りであり、これを総合して吸気管部分２３から全気筒の燃焼室５に流入する空気の流量は図１２に実線で示した通りである。また、例えば１番気筒への筒内充填空気量Ｍｃは図１２に斜線で示した通りである。

これに対して、実線で示した吸気管部分２３から全ての気筒の燃焼室５に流入する空気の流量を平均化したものが筒内吸入空気流量ｍｃであり、図中に一点鎖線で示す。そして、この一点鎖線で示した筒内吸入空気流量ｍｃに、４気筒の場合にはクランクシャフトが１８０°（すなわち、４ストローク式内燃機関において１サイクル中にクランクシャフトが回転する角度７２０°を気筒数で割った角度）回転するのにかかる時間ΔＴ180°を乗算したものが筒内充填空気量Ｍｃとなる。したがって、吸気弁モデルＭ３０で算出された筒内吸入空気流量ｍｃにΔＴ180°を乗算することで、筒内充填空気量Ｍｃが算出される（Ｍｃ＝ｍｃ・ΔＴ180°）。より詳細には、筒内充填空気量Ｍｃが吸気弁閉弁時の圧力に比例することを考慮して、吸気弁閉弁時の筒内吸入空気流量ｍｃにΔＴ180°を乗算したものが筒内充填空気量Ｍｃとされる。

次に、上記エアモデルを内燃機関の制御装置に実装して、実際に筒内充填空気量Ｍｃを算出する場合について説明する。筒内充填空気量Ｍｃはエアモデルを用いて、上記式（１）、式（３）、式（４）、および式（７）を解くことにより表される。この場合、ＥＣＵ３１で処理するために、これらの式を離散化する必要がある。時刻ｔ、計算間隔Δｔを用いて式（１）、式（３）、式（４）、および式（７）を離散化すると、それぞれ下記式（８）、式（９）、式（１０）、および式（１１）が得られる。なお、吸気管内温度Ｔｍ（ｔ＋Δｔ）は、式（９）および式（１０）によってそれぞれ算出されたＰｍ／Ｔｍ（ｔ＋Δｔ）およびＰｍ（ｔ＋Δｔ）から、式（１２）によって算出される。

このようにして実装されたエアモデルでは、スロットルモデルＭ１０の式（８）で算出された時刻ｔにおけるスロットル通過空気流量の推定値ｍｔ（ｔ）と、吸気弁モデルＭ３０の式（１１）で算出された時刻ｔにおける筒内吸入空気流量ｍｃ（ｔ）とが、吸気管モデルＭ２０の式（９）および式（１０）に代入され、これにより時刻ｔ＋Δｔにおける吸気管内圧力Ｐｍ（ｔ＋Δｔ）および吸気管内温度Ｔｍ（ｔ＋Δｔ）が算出される。次いで、算出されたＰｍ（ｔ＋Δｔ）およびＴｍ（ｔ＋Δｔ）は、スロットルモデルＭ１０および吸気弁モデルＭ３０の式（８）および式（１１）に代入され、これにより時刻ｔ＋Δｔにおけるスロットル通過空気流量の推定値ｍｔ（ｔ＋Δｔ）および筒内吸入空気流量ｍｃ（ｔ＋Δｔ）が算出される。そして、このような計算を繰り返すことによって、スロットル弁開度θｔ、大気圧Ｐａ、および大気温度Ｔａから、任意の時刻ｔにおける筒内吸入空気流量ｍｃが算出され、算出された筒内吸入空気流量ｍｃに上記時間ΔＴ180°を乗算することで、任意の時刻ｔにおける筒内充填空気量Ｍｃが算出される。

なお、内燃機関の始動時には、すなわち時刻ｔ＝０においては、吸気管内圧力Ｐｍは大気圧と等しい（Ｐｍ（０）＝Ｐａ）とされ、吸気管内温度Ｔｍは大気温度と等しい（Ｔｍ（０）＝Ｔａ）とされて、各モデルＭ１０〜Ｍ３０における計算が開始される。

なお、上記エアモデルでは、大気温度Ｔａおよび大気圧Ｐａが一定であるとしているが、時刻によって変化する値としてもよく、例えば、大気温度を検出するための吸気温度センサによって時刻ｔにおいて検出された値を大気温度Ｔａ（ｔ）、大気圧を検出するための大気圧センサによって時刻ｔにおいて検出された値を大気圧Ｐａ（ｔ）として上記式（８）、式（１０）、および式（１１）に代入するようにしてもよい。

ところで、上記したように図５の装置ではスロットル弁開度θｔ、大気温度Ｔａおよび大気圧Ｐａのみを入力するモデル式を用いて機関始動時からの逐次計算により筒内充填吸気量が算出される。このため、エアモデルの各要素モデルにモデル化誤差が生じていると、算出される筒内充填吸気量は不正確になり燃料噴射量の正確な算出ができなくなる。

一方、エアモデルの各要素モデル（スロットルモデルＭ１０、吸気管モデルＭ２０および吸気弁モデルＭ３０）を考えると、モデル化そのものには誤差がない場合であっても、使用による特性の変化や製作公差によるばらつきなどによりモデル化誤差が生じている場合がある。このため、エアモデルにより求めた筒内充填吸気量等の値をエアフロメータや吸気管内圧力から求めた吸気量と比較し、モデル式から計算した吸気量をエアフロメータなどにより求めた実測値に合致するように修正することが行われている。

しかし、このような従来の修正方法は各要素モデルの集合であるエアモデルの全体を一律に修正するものであるのに対して、実際には各モデルに発生する誤差は一律ではなく、それぞれのモデル毎に相違している。このため、全体を一律に修正する方法では正確なモデルの修正を行うことはできず、例えば内燃機関が定常運転を行っているときに筒内充填吸気量を実測値に合致するように補正すると、内燃機関が過渡運転を行っているときなどでは吸気圧力の計算値が実測値から大きく離れてしまい、過渡運転を行っているときの筒内充填吸気量に大きな誤差が生じるような場合がある。

そこで、本実施形態では、実測値に基づいてモデル修正するものの、エアモデル全体（Ｍ１０〜Ｍ３０）を一律に修正するのではなく、それぞれの要素モデルに応じた補正を別個に行う。

ここで、エアモデルに含まれる各要素モデルのうち、スロットルモデルＭ１０と吸気弁モデルＭ３０とは製作・組み立て公差によるばらつきや、使用による特性変化（流量係数の変化等）が生じやすい。これに対して、吸気管モデルＭ２０は単に容積のみが問題となる比較的単純なモデルであり、ばらつきや使用による特性変化は生じにくいため、当初のモデル化さえ正確にできていれば、運転中にこのモデルを修正する必要はない。したがって、本実施形態ではモデル化誤差が生じやすいスロットルモデルＭ１０と吸気弁モデルＭ３０とを機関運転中に実測値に基づいて個別にモデル式を修正するようにしている。

以下、本実施形態のモデル修正方法について説明する。本実施形態では、現在の実際のエアフロ通過空気流量、スロットル弁開度、吸気管内圧力を、それぞれエアフロメータ４１、スロットル弁開度センサ４４、吸気圧センサ４０で実測し、これらの値を用いて図５の筒内充填吸気量計算とは独立してスロットルモデルＭ１０、吸気弁モデルＭ３０等についての修正操作を行う。

まず、スロットルモデルＭ１０の修正操作について説明する。図１３は、スロットルモデルＭ１０の修正操作を説明するブロック図である。図１３に示したように、スロットルモデルの修正操作を行う際には、エアフロメータモデルＭ４０が用いられる。エアフロメータモデルＭ４０には、スロットルモデルＭ１０によって算出されたスロットル通過空気流量（以下、「スロットル通過空気流量の推定値」と称す）ｍｔが入力され、この入力されたパラメータの値を後述するエアフロメータモデルＭ４０のモデル式に代入することで、実際にスロットル弁１９を通過する空気の流量（以下、「実際のスロットル通過空気流量」と称す）がこの値のときにエアフロメータ６１が出力するであろう出力値が算出される。

すなわち、実際のエアフロメータ４１の出力は、実際のエアフロ通過空気流量（エアフロ通過空気流量とスロットル通過空気流量はほぼ同一であると考えられるため、以下ではスロットル通過空気流量として説明する）に対して固有の応答特性に基づく応答遅れを有しており、したがってそのときの実際のスロットル通過空気流量のみに対応した値となっているわけではない。エアフロメータモデルＭ４０は、この応答特性をシュミレートしたモデルであり、上記エアフロメータ４１の応答遅れを考慮した上で、スロットル弁１９をスロットル通過空気流量の推定値ｍｔだけ空気が流れている場合にエアフロメータ４１が出力するであろう出力値（以下、「予想出力値」と称す）を算出する。

ここで、エアフロメータモデルＭ４０について具体的に説明する。エアフロメータ４１では、上述したように吸気温度計測用抵抗４１ａ１と加熱用抵抗（ボビン部）４１ａ３との温度差を常に一定に維持するように加熱用抵抗４１ａ３に供給する電力を調整するとともに、その供給電力に基づいてエアフロ通過空気流量を算出するようにしている。ここで、この供給電力はボビン部４１ａ３およびサポート部４１ａ４から吸気管１６内を通過する空気への放熱量を示すことから、エアフロメータ４１はボビン部４１ａ３およびサポート部４１ａ４からの放熱量に基づいてエアフロ通過空気流量を算出すると言い換えることができる。

ここで、ボビン部４１ａ３は、より詳細には円筒状のセラミックスボビンに白金熱線を巻回し、その外周にガラスをコーティングすることにより形成される。このため、白金熱線から周囲の空気への放熱は白金熱線と周囲の空気との間にガラス層が介在することにより遅れてしまう。したがって、吸気管１６内を通過する空気の流量が急激に増大したような場合であっても、ボビン部４１ａ３からの単位時間当たりの放熱量（以下、単に「放熱量」と称す）は直ぐには増大せず、或る程度の遅れをもって増大することになる。換言すると、ボビン部４１ａ３からの放熱量は実際のエアフロ通過空気流量に対して応答遅れが存在する。また、同様なことがサポート部４１ａ４からの放熱量についても言え、サポート部４１ａ４からの放熱量は実際のエアフロ通過空気流量に対して応答遅れが存在する。

この応答遅れは、一次遅れに近似することができることがわかっており、ボビン部４１ａ３の放熱量の応答遅れは次の式（１３）で表される。ここで、式（１３）におけるω_bはエアフロメータ４１の出力値から換算したボビン部４１ａ３の放熱量、すなわち放熱遅れが生じた結果ボビン部４１ａ３の白金熱線から実際に放熱される放熱量（以下、「遅れ放熱量」と称す）を示している。また、式（１３）におけるＷ_bは、応答遅れを補償した放熱量、すなわち放熱遅れが生じないと仮定した場合におけるボビン部４１ａ３の白金熱線から放熱される放熱量（以下、「完全放熱量」と称す）を示している。すなわち、完全放熱量Ｗ_bは、内燃機関が定常運転を行っているときにおける放熱量に等しく、基本的にエアフロメータ４１近傍の吸気管１６内を通過する空気流量のみの関数である。エアフロメータ４１近傍の吸気管１６内を通過する空気流量はスロットル通過空気流とほぼ等しいため、完全放熱量Ｗ_bはスロットル通過空気流量の関数と考えることができ、スロットル通過空気流量とボビン部４１ａ３からの完全放熱量Ｗ_bとの関係は図１４のように表すことができる。さらに、式（１３）におけるτ_bは、ボビン部４１ａ３からの放熱における一次遅れの時定数であり、その算出方法については後述する。

同様に、サポート部４１ａ４の放熱量の応答遅れは式（１４）で表される。サポート部４１ａ４からの完全放熱量Ｗ_bもスロットル通過空気流量の関数と考えることができ、スロットル通過空気流量とサポート部４１ａ４からの完全放熱量Ｗ_sとの関係は図１４のように表すことができる。また、式（１４）におけるτ_sは、サポート部４１ａ４からの放熱における一次遅れの時定数である。

式（１３）および式（１４）によりボビン部４１ａ３およびサポート部４１ａ４からの遅れ放熱量ω_b、ω_sが算出される。エアフロメータ４１では、放熱遅れが生じた結果ボビン部４１ａ３およびサポート部４１ａ４から実際に放熱される放熱量、すなわち遅れ放熱量に応じて出力電圧が変化するため、式（１３）および式（１４）によって算出された遅れ放熱量の和（ω_b＋ω_s）に基づいてエアフロメータ４１が出力するであろう出力値（予想出力値）ＡＦＭｍｔが算出される。本実施形態では、エアフロメータ４１における遅れ放熱量の和（ω_b＋ω_s）とエアフロメータ４１の予想出力値ＡＦＭｍｔの関係を予め実験的にまたは計算によって求め、図１５に示したようなマップとしてＥＣＵ３１のＲＯＭ３４に記憶させておく。そして、式（１３）および式（１４）から算出された遅れ放熱量の和ω_b＋ω_sに基づいて上記マップを用いてエアフロメータ４１の予想出力値ＡＦＭｍｔが算出される。なお、上記マップを用いた遅れ放熱量の和ω_b＋ω_sに基づくエアフロメータ４１の予想出力値ＡＦＭｍｔの算出は下記式（１５）のように表すことができる。

なお、ボビン部４１ａ３の一次遅れの時定数τ_bは下記式（１６）によって、またサポート部４１ａ４の一次遅れの時定数τ_sは下記式（１７）によって算出される。式（１６）および式（１７）において、ｕは、エアフロメータ４１の検出部における流路、すなわちエアフロメータ４１のバイパス通路における単位断面積当たりの空気流量である。単位断面積当たりの空気流量ｕはエアフロメータ４１の出力値ＡＦＭに基づいてマップによりまたは所定の計算式により算出される。また、ｋ_b、ｋ_s、ｍ_b、ｍ_sは予め実験または計算によって求められる定数であり、ｋ_b、ｍ_bがボビン部４１ａ３についての定数、ｋ_s、ｍ_sがサポート部４１ａ４についての定数をそれぞれ示している。ボビン部４１ａ３とサポート部４１ａ４とでは応答遅れの度合いが異なるので、ボビン部４１ａ３とサポート部４１ａ４とを分離して時定数を設定することによってスロットル通過空気流量から予想出力値を算出するための算出精度を向上させることとしている。

次に、上記エアフロメータモデルＭ４０を内燃機関の制御装置に実装して、実際にスロットルモデルＭ１０によって算出されたスロットル通過空気流量の推定値ｍｔからエアフロメータ４１の予想出力値ＡＦＭｍｔを算出する場合について説明する。エアフロメータ４１の予想出力値ＡＦＭｍｔはエアフロメータモデルＭ４０を用いて、上記式（１３）、式（１４）を解くことにより表される。この場合、ＥＣＵ３１で処理するために、これらの式を離散化する必要がある。時刻ｔ、計算間隔Δｔを用いて式（１３）、式（１４）を離散化すると、それぞれ下記式（１８）、式（１９）が得られる。

このようにして実装されたエアフロメータモデルＭ４０では、時刻ｔにおけるボビン部４１ａ３からの遅れ放熱量ω_b（ｔ）、サポート部４１ａ４からの遅れ放熱量ω_s（ｔ）およびスロットルモデルＭ１０の式（８）で算出された時刻ｔにおけるスロットル通過空気流量の推定値ｍｔ（ｔ）が式（１８）および式（１９）に代入され、これにより時刻ｔ＋Δｔにおけるボビン部４１ａ３からの遅れ放熱量ω_b（ｔ＋Δｔ）およびサポート部４１ａ４からの遅れ放熱量ω_s（ｔ＋Δｔ）が算出される。そして、これら遅れ放熱量ω_b（ｔ＋Δｔ）およびω_s（ｔ＋Δｔ）を用いて式（１５）により時刻ｔ＋Δｔにおけるエアフロメータ４１の予想出力値ＡＦＭｍｔ（ｔ＋Δｔ）が算出される。

再びスロットルモデルの修正操作について話を戻すと、本実施形態では、吸気圧センサ４０とスロットル弁開度センサ４４とで実測した現在の吸気管内圧力Ｐｍとスロットル弁開度θｔ（および大気圧Ｐａ、大気温度Ｔａ）とを用いてスロットルモデルＭ１０（図５）によって算出されるスロットル通過空気流量の推定値ｍｔが、エアフロメータ４１の出力から算出した現在のスロットル通過空気流量に一致するようにスロットルモデルの式（１）を修正する。

すなわち、図１３において、実測値Ｐｍとスロットル弁開度θｔとを用いてスロットルモデルＭ１０の式（１）によって算出されたスロットル通過空気流量の推定値ｍｔをエアフロメータモデルＭ４０に入力すると、実際のスロットル通過空気流量が上記スロットル通過空気流量の推定値ｍｔとなっていたと仮定した場合におけるエアフロメータ４１の予想出力値ＡＦＭｍｔに変換される。

したがって、スロットルモデルＭ１０にモデル化誤差が含まれていない場合、すなわち実際のスロットル通過空気流量とスロットルモデルＭ１０によって算出されたスロットル通過空気流量の推定値ｍｔとが一致している場合には、実際のエアフロメータ４１の出力値ＡＦＭとエアフロメータモデルＭ４０の予想出力値ＡＦＭｍｔとが一致するはずであり、一致していない場合には実際のスロットル通過空気流量とスロットル通過空気流量の推定値ｍｔとが一致していないと考えられ、よってスロットルモデルＭ１０にモデル化誤差が生じていると考えられる。

したがって、出力値ＡＦＭと予想出力値ＡＦＭｍｔとが一致しない場合には、スロットルモデルＭ１０に生じているモデル化誤差を補償するために、予想出力値ＡＦＭｍｔがエアフロメータ４１の出力値（実測値）ＡＦＭに一致するようにスロットルモデルＭ１０の計算式を修正する必要がある。

ここで、本実施形態では、上述したようにスロットルモデルＭ１０の計算式として、下記式（２０）を用いている。

上記式（２０）のパラメータのうち流量係数μｔは製品毎のばらつきや使用による特性変化の影響を受けやすい。したがって、本実施形態では、スロットルモデルＭ１０の式（２０）により算出されたスロットル通過空気流量の推定値ｍｔをエアフロメータモデルＭ４０の式（１３）および式（１４）に代入して算出された予想出力値ＡＦＭｍｔがエアフロメータ４１の出力値ＡＦＭに一致するように流量係数μｔの値を修正することとしている。これにより、スロットルモデルＭ１０が修正され、スロットルモデルＭ１０を用いた計算結果をエアフロメータ４１の出力値と一致させることができるようになり、よってスロットルモデルＭ１０によって算出されたスロットル通過空気流量の推定値ｍｔを実際のスロットル通過空気流量に一致させることができるようになる。なお、当然ながらスロットルモデルＭ１０の上記修正結果は図５筒内充填吸気量Ｍｃの計算にも同時に反映され、筒内充填吸気量Ｍｃの推定精度が向上するようになる。

このような流量係数μｔの値の修正方法としては、例えば、下記式（２１）によるものが考えられる。なお、式（２１）におけるμｔ’は修正後の流量係数である。

ところが、式（２１）により流量係数μｔの値を修正すると、内燃機関が定常運転を行っているときには適切に修正することができるが、内燃機関が過渡運転を行っているときには適切に修正することができない。以下、このことについて図１６を参照して簡単に説明する。

図１６は、スロットル弁開度θｔと空気流量に関する各種パラメータの値とのタイムチャートである。図中、時刻ｔ₀はスロットル弁開度θｔがθｔ₀となったとき、時刻ｔ₁はスロットル弁開度θｔがθｔ₁となったときを示している。また、図１６は、実際のスロットル通過空気流量（図１６中の実線）と、スロットルモデルＭ１０によって算出されたスロットル通過空気流量の推定値ｍｔ（図１６中の一点鎖線）とが時刻ｔ₀前まではずれているが、時刻ｔ₀以降はほぼ一致している場合のタイムチャートを示している。すなわち、図１６は、時刻ｔ₀前、すなわちスロットル弁開度θｔがθｔ₀よりも小さいときには、流量係数μｔの値に誤差が生じていてスロットルモデルＭ１０によって算出されたスロットル通過空気流量の推定値ｍｔが実際のスロットル通過空気流量からずれてしまっていると共に、時刻ｔ₀以降、すなわちスロットル弁開度θｔがθｔ₀以上であるときには、流量係数μｔの値に誤差がほとんどなく、スロットルモデルＭ１０によって算出されたスロットル通過空気流量の推定値ｍｔが実際のスロットル通過空気流量にほぼ一致している場合のタイムチャートを示している。

上述したように、エアフロメータ４１には放熱遅れに基づく応答遅れがあるため、図１６に示したように実際のスロットル通過空気流量が変化している場合に、エアフロメータ４１の出力値ＡＦＭ（図１６中の破線）には実際のスロットル通過空気流量に対して遅れが生じる。特に、エアフロメータ４１の出力値ＡＦＭを実際のスロットル通過空気流量に対して一次遅れで近似することができることを考えると、エアフロメータ４１の現在の出力値ＡＦＭは現在の実際のスロットル通過空気流量のみに応じて変化するわけではなく、過去の実際のスロットル通過空気流量の履歴にも応じて変化する。

また、上述したようにエアフロメータモデルＭ４０ではエアフロメータ４１の放熱遅れに基づく応答遅れが考慮されているため、図１６に示したように実際のスロットル通過空気流量が変化している場合には、エアフロメータモデルＭ４０によって算出される予想出力値ＡＦＭｍｔにもスロットルモデルＭ１０によって算出されたスロットル通過空気流量の推定値ｍｔに対して遅れが生じている。

このため、図１６に示したように、時刻ｔ₀以降においてスロットルモデルＭ１０によって算出されたスロットル通過空気流量の推定値ｍｔが実際のスロットル通過空気流量にほぼ一致しているにも関わらず、時刻ｔ₀以降においてエアフロメータ４１の出力値ＡＦＭとエアフロメータモデルＭ４０によって算出された予想出力値ＡＦＭｍｔとは一致していない。したがって、例えば、出力値ＡＦＭと予想出力値ＡＦＭｍｔとを用いて式（２１）によってスロットルモデルＭ１０の流量係数μｔを修正しようとすると、時刻ｔ₀以降において流量係数μｔには誤差が無いにも関わらず、流量係数μｔが修正されてしまうことになる。

或いは、流量係数μｔの値の修正方法としては、例えば、下記式（２２）によるものが考えられる。このように、積分を用いることにより、エアフロメータ４１の出力値ＡＦＭに存在するノイズの影響を抑制することができる。

ところが、式（２２）により流量係数μｔの値を修正する場合であっても、内燃機関が定常運転を行っているときには適切に修正することができるが、内燃機関が過渡運転を行っているときには適切に修正することができない。以下、このことについて簡単に説明する。

ここで、流量係数μｔのマップは一定の機関運転領域に区切られており、流量係数μｔのマップ値は機関運転領域毎に設定される。具体的には、流量係数μｔのマップは図１７に示したようにスロットル弁開度θｔを複数のマップ領域に分割し、スロットル弁開度θｔの各マップ領域毎に一つの流量係数μｔのマップ値が設定される。したがって、流量係数μｔを修正する場合にはスロットル弁開度θｔの各マップ領域毎にそのマップ領域に対応する流量係数μｔの値を修正することになる。

図１６に示した例を参照すると、スロットル弁開度θｔの一つのマップ領域がθｔ₀からθｔ₁である場合、スロットル弁開度θｔがそのマップ領域内にある期間は時刻ｔ₀から時刻ｔ₁までとなっている。したがって、このスロットル弁開度θｔがこのマップ領域にある場合における流量係数μｔの値の誤差率を、時刻ｔ₀から時刻ｔ₁までのエアフロメータ４１の出力値ＡＦＭの平均値とエアフロメータモデルＭ４０によって算出された予想出力値ＡＦＭｍｔの平均値との比に基づいて算出すると、比較的正確に流量係数μｔの値の誤差を補償することができると考えられ、斯かる考えに基づいて式（２２）による流量係数μｔの値の修正が提案されている。より詳細には、式（２２）では、エアフロメータ４１の出力値ＡＦＭの平均値とその予想出力値ＡＦＭｍｔの平均値との比は、エアフロメータ４１の出力値ＡＦＭの積分値とその予想出力値ＡＦＭｍｔの積分値との比に等しいため、式（２２）では現在マップに保存されている流量係数μｔ（θｔ）に、エアフロメータ４１の出力値ＡＦＭの積分値（図１６中の面積Ｓ₁＋Ｓ₂）を予想出力値ＡＦＭｍｔの積分値（図１６中の面積Ｓ₁）で除算したものを乗算することによって（θｔ・（Ｓ₁＋Ｓ₂）／Ｓ₁）、流量係数μｔ（θｔ）を修正した値μｔ’（θｔ）を算出している。

ところが、上述したように、図１６に示した例では、時刻ｔ₀以降は実際のスロットル通過空気流量とスロットルモデルＭ４０により算出されたスロットル通過空気流量の推定値ｍｔとがほぼ一致しているにも関わらず、エアフロメータ４１に応答遅れが存在するため、エアフロメータ４１の出力値ＡＦＭの時刻ｔ₀から時刻ｔ₁までの積分値（Ｓ₁＋Ｓ₂）とエアフロメータモデルＭ４０によって算出された予想出力値ＡＦＭｍｔの時刻ｔ₀から時刻ｔ₁までの積分値（Ｓ₁）とが一致しない。このため、スロットル弁開度θｔがθｔ₁からθｔ₂であるマップ領域における流量係数μｔには誤差が無いにも関わらず流量係数μｔが修正されてしまうことになる。

このように、式（２１）および式（２２）により流量係数μｔの値を修正すると、内燃機関が定常運転を行っているときには的確に修正することができるが、内燃機関が過渡運転を行っているときには的確に修正することができない。そこで、本発明では、内燃機関が過渡運転を行っているときにおいても的確に流量係数μｔの値を修正することができるような内燃機関の吸入空気量推定装置を提供する。

図１８は、図１６と同様に、スロットル弁開度θｔと空気流量に関する各種パラメータの値とのタイムチャートである。以下では、図１８を参照して、スロットル弁開度θｔがθｔ₀からθｔ₁であるマップ領域Ａにおける、流量係数μｔの修正について考える。式（１３）および式（１４）に示したように、エアフロメータモデルＭ４０では、スロットルモデルＭ１０によって算出されたスロットル通過空気流量の推定値ｍｔに対して一次遅れが存在するとしてエアフロメータ４１の予想出力値ＡＦＭｍｔを算出している。このため、エアフロメータモデルＭ４０の予想出力値ＡＦＭｍｔは、スロットルモデルＭ１０によって算出された現在のスロットル通過空気流量の推定値ｍｔだけでなく所定時間前のスロットル通過空気流量の推定値ｍｔに基づいて算出されることになる。

例えば、図１８に示したように時刻ｔ₀においてスロットル弁開度θｔがマップ領域Ａに侵入したとき、時刻ｔ₀以降においてエアフロメータモデルＭ４０によって算出された予想出力値ＡＦＭｍｔは、時刻ｔ₀以降においてスロットルモデルＭ１０によって算出されたスロットル通過空気流量の推定値ｍｔだけでなく、時刻ｔ₀以前にスロットルモデルＭ１０によって算出されたスロットル通過空気流量の推定値ｍｔを用いて算出されることになる。すなわち、エアフロメータモデルＭ４０によって算出される予想出力値ＡＦＭｍｔは、スロットル弁開度θｔがマップ領域Ａ以外のマップ領域にあるときのスロットル通過空気流量の推定値ｍｔに基づいて算出されることになる。

ここで、上述したように図１８に示した例では、スロットル弁開度θｔがθｔ₀からθｔ₁であるマップ領域Ａに対応する流量係数μｔの修正を行おうとしており、流量係数μｔの修正を行うに当たって実際のスロットル通過空気流量とスロットルモデルＭ１０によって算出されたスロットル通過空気流量の推定値ｍｔとの比をエアフロメータ４１の出力値ＡＦＭとエアフロメータモデルＭ４０によって算出された予想出力値ＡＦＭｍｔとの比から算出しようとしている。ところが、スロットル弁開度θｔがマップ領域Ａにある場合においてエアフロメータモデルＭ４０によって算出される予想出力値ＡＦＭｍｔは、マップ領域Ａ以外のマップ領域におけるスロットル通過空気流量の推定値ｍｔの影響を受けており、よってスロットル弁開度θｔがマップ領域Ａにある場合におけるエアフロメータ４１の出力値ＡＦＭとエアフロメータモデルＭ４０によって算出される予想出力値ＡＦＭｍｔとの比は、スロットル弁開度θｔがマップ領域Ａにある場合における実際のスロットル通過空気流量とスロットル通過空気流量の推定値ｍｔとの比を表したものではないものとなっている。

そこで、本実施形態では、スロットル弁開度θｔがマップ領域Ａにある場合における流量係数μｔを修正するときには、エアフロメータモデルＭ４０によって算出された予想出力値ＡＦＭｍｔを直接利用せずに、この予想出力値ＡＦＭｍｔからスロットル弁開度θｔがマップ領域Ａに侵入する前（図１８に示した例では、スロットル弁開度θｔがθｔ₀よりも小さいとき）のスロットル通過空気流量の推定値ｍｔの影響を除去した値（以下、「修正出力値」と称す）ＡＦＭｍｔｅを用いることとする。

本実施形態では、修正出力値ＡＦＭｍｔｅの算出は基本的にエアフロメータモデルＭ４０における予想出力値ＡＦＭｍｔの算出と同様に行われる。すなわち、下記式（２３）および式（２４）によりボビン部４１ａ３およびサポート部４１ａ４の遅れ放熱量ω_be、ω_seが算出され、算出された遅れ放熱量ω_be、ω_seに基づいて式（２５）により図１５に示したマップと同様なマップから修正出力値ＡＦＭｍｔｅが算出される。式（２３）および式（２４）は、ＥＣＵ３１で処理するために時刻ｔ、計算間隔Δｔを用いて下記式（２６）および式（２７）のように離散化されて用いられる。

ただし、エアフロメータモデルＭ４０における予想出力値ＡＦＭｍｔの算出と異なり、修正出力値ＡＦＭｍｔｅの算出はマップ領域毎に行われる。すなわち、スロットル弁開度θｔが或るマップ領域に侵入したとき、図１８に示した例ではスロットル弁開度θｔがマップ領域Ａに侵入したとき（すなわちスロットル弁開度θｔがθｔ₀となったとき）にそのマップ領域における修正出力値ＡＦＭｍｔｅの算出が開始されると共に、スロットル弁開度θｔがそのマップ領域から離脱したとき、図１８に示した例ではスロットル弁開度θｔがマップ領域Ａから離脱したとき（すなわちスロットル弁開度θｔがθｔ₁となったとき）にそのマップ領域における修正出力値ＡＦＭｍｔｅの算出が終了せしめられる。

そして、修正出力値ＡＦＭｍｔｅの算出にあたっては、スロットル弁開度θｔがマップ領域Ａに侵入する前にスロットルモデルＭ１０によって算出されたスロットル通過空気流量の推定値ｍｔ（以下、「過去のスロットル通過空気流量の推定値ｍｔ」と称す）の影響を除去するために、スロットル弁開度θｔがマップ領域Ａに侵入したときの修正出力値ＡＦＭｍｔｅの値（以下、「初期値」と称す）を、エアフロメータ４１の出力値ＡＦＭとほぼ同様な値になるようにしている。

すなわち、上述したように修正出力値ＡＦＭｍｔｅはＥＣＵ３１において上記離散化された式（２６）および式（２７）に基づいて算出されるが、これら式（２６）および式（２７）には初期値以外に過去の（スロットル弁開度θｔがマップ領域Ａに侵入する前の）スロットル通過空気流量の推定値ｍｔが影響するような項は存在しない。逆に言うと、初期値を過去のスロットル通過空気流量の推定値ｍｔの影響を受けていない値に設定しさえすれば、修正出力値ＡＦＭｍｔｅに対する過去のスロットル通過空気流量の推定値ｍｔの影響を最小限に抑制することができる。

本実施形態では、スロットル弁開度θｔがマップ領域Ａに侵入したときの修正出力値ＡＦＭｍｔｅの初期値をそのときのエアフロメータ４１の出力値ＡＦＭとほぼ同様な値としているため、修正出力値ＡＦＭｍｔｅの初期値は過去のスロットル通過空気流量の推定値ｍｔに基づいて定まる値ではなく、スロットル弁開度θｔがマップ領域Ａに侵入したときのエアフロメータ４１の出力値ＡＦＭに基づいて定まる値であり、よって修正出力値ＡＦＭｍｔｅから過去のスロットル通過空気流量の推定値ｍｔの影響が除去される。

このように、修正出力値ＡＦＭｍｔｅの初期値をエアフロメータ４１の出力値ＡＦＭとほぼ同様な値として式（２３）および式（２４）によってその後の修正出力値ＡＦＭｍｔを算出するためには、時刻ｔ₀におけるボビン部４１ａ３の遅れ放熱量およびサポート部４１ａ４の遅れ放熱量を時刻ｔ₀におけるボビン部４１ａ３からの実際の放熱量およびサポート部４１ａ４からの実際の放熱量に一致させることが必要である。ところが、時刻ｔ₀におけるエアフロメータ４１の出力値ＡＦＭから時刻ｔ₀におけるボビン部４１ａ３およびサポート部４１ａ４からの実際の放熱量を算出することは困難である。すなわち、図１５に示したマップ等によってエアフロメータ４１の出力値ＡＦＭからボビン部４１ａ３およびサポート部４１ａ４からの実際の放熱量を合計した値を求めることはできるが、ボビン部４１ａ３およびサポー部４１ａ４それぞれからの実際の放熱量を求めることはできない。

そこで、本実施形態では、エアフロメータ４１のボビン部４１ａ３およびサポート部４１ａ４それぞれからの実際の放熱量を比較的正確に算出することができるオブザーバを用いるようにしている。このオブザーバは、時刻やスロットル弁開度θｔのマップ領域とは無関係に、常にエアフロメータ４１のボビン部４１ａ３およびサポート部４１ａ４それぞれからの実際の放熱量の近似値を算出している。

具体的には、オブザーバは、下記式（２８）および式（２９）によってボビン部４１ａ３からの実際の放熱量の近似値（以下、「近似放熱量」と称す）ω_boおよびサポート部４１ａ４からの実際の放熱量の近似値（以下、「近似放熱量」と称す）ω_soをそれぞれ算出している。式（２８）および式（２９）において、Ｋ_Iは積分制御等に用いられる積分係数と同様な係数である。また、ＡＦＭｍｔｏはオブザーバによって算出された近似放熱量ω_bo、ω_soを合計した値に基づいて下記式（３０）により（すなわち、図１５のマップにより）算出された値である。

したがって、オブザーバでは、エアフロメータ４１の出力値ＡＦＭとオブザーバによって算出された近似放熱量ω_bo、ω_soに基づいて算出されるエアフロメータ４１の出力値（以下、「近似出力値」と称す）ＡＦＭｍｔｏとの偏差（ＡＦＭ−ＡＦＭｍｔｏ）を時間的に累積する。この出力値ＡＦＭと近似出力値ＡＦＭｍｔｏとの偏差の累積値は、実際のスロットル通過空気流量とスロットルモデルＭ１０によって算出されたスロットル通過空気流量の推定値ｍｔとの誤差に基づいて変化し、この誤差の絶対値が大きいほど上記偏差の累積値の絶対値が大きくなる。したがって、スロットルモデルＭ１０によって算出されたスロットル通過空気流量の推定値ｍｔに、出力値ＡＦＭとオブザーバによって算出された近似出力値ＡＦＭｍｔｏとの偏差の累積値に係数Ｋ_Iを乗じた値を加算することで、実際のスロットル通過空気流量とスロットルモデルＭ１０によって算出されたスロットル通過空気流量の推定値ｍｔとの誤差を或る程度補償することができる。これにより、オブザーバでは、実際のスロットル通過空気流量とほぼ等しいスロットル通過空気流量に基づいて近似放熱量ω_bo、ω_soが算出されるため、オブザーバによって算出された近似放熱量ω_bo、ω_soはボビン部４１ａ３およびサポート部４１ａ４からの実際の放熱量とほぼ等しい値となっている。

なお、上記オブザーバを内燃機関の制御装置に実装して、スロットル通過空気流量の推定値ｍｔおよびエアフロメータ４１の出力値ＡＦＭから近似放熱量ω_bo、ω_soを算出するためには、式（２８）および式（２９）を離散化する必要がある。時刻ｔ、計算間隔Δｔを用いて式（２８）および式（２９）を離散化すると、それぞれ下記式（３１）および式（３２）が得られる。なお、式（３１）および式（３２）中のΔＡＦＭＩＮＴ（ｔ）は、時刻ｔにおけるＡＦＭ（ｔ）−ＡＦＭｍｔｏ（ｔ）の累積値であり、下記式（３３）により表される。

このようにして実装されたオブザーバでは、時刻ｔにおけるボビン部４１ａ３およびサポート部４１ａ４からの近似放熱量ω_bo（ｔ）およびω_so（ｔ）、スロットルモデルＭ１０の式（８）によって算出された時刻ｔにおけるスロットル通過空気流量の推定値ｍｔ（ｔ）、時刻ｔにおけるエアフロメータ４１の出力値ＡＦＭ（ｔ）、および時刻ｔにおいてオブザーバによって算出された近似出力値ＡＦＭｍｔｏ（ｔ）が式（３１）、式（３２）および式（３３）に代入され、これにより時刻ｔ＋Δｔにおける近似放熱量ω_bo（ｔ＋Δｔ）およびω_so（ｔ＋Δｔ）が算出される。そして、これら近似放熱量ω_bo（ｔ＋Δｔ）およびω_so（ｔ＋Δｔ）に基づいて式（３０）により時刻ｔ＋Δｔにおける近似出力値ＡＦＭｍｔｏ（ｔ＋Δｔ）が算出され、この時刻ｔ＋Δｔにおける近似出力値ＡＦＭｍｔｏ（ｔ＋Δｔ）はその後の近似放熱量ω_boおよびω_soの算出に利用される。そして、このような計算を繰り返すことによって、スロットル通過空気流量の推定値ｍｔ、エアフロメータ４１の出力値ＡＦＭから、任意の時刻ｔにおける近似放熱量ω_bo、ω_soが算出される。

図１９は、図１８の時刻ｔ₀〜ｔ₁における流量に関するパラメータの推移を拡大して示したタイムチャートである。図１９から分かるように、オブザーバによって算出された近似出力値ＡＦＭｍｔｏ（図１９中の短い破線）は、エアフロメータ４１の出力値ＡＦＭとほぼ同様な値となっている。上述したように修正出力値ＡＦＭｍｔｅの算出にあたっては、スロットル弁開度θｔがマップ領域Ａに侵入したときに、すなわち時刻ｔ₀においてオブザーバによって算出されたボビン部４１ａ３の近似放熱量ω_boおよびサポート部４１ａ４の近似放熱量ω_soが用いられる。したがって、修正出力値ＡＦＭｍｔｅの算出においては、遅れ放熱量ω_be、ω_seの初期値、すなわち時刻ｔ₀における遅れ放熱量ω_be（ｔ₀）、ω_se（ｔ₀）はオブザーバによって算出された時刻ｔ₀における近似放熱量ω_bo（ｔ₀）、ω_se（ｔ₀）とされる（ω_be（ｔ₀）＝ω_bo（ｔ₀）、ω_se（ｔ₀）＝ω_so（ｔ₀））。このため、時刻ｔ₀における修正出力値ＡＦＭｍｔｅは、時刻ｔ₀におけるオブザーバの近似出力値ＡＦＭｍｔｏと同一の値となる。

そして、修正出力値ＡＦＭｍｔｅの算出に当たっては、遅れ放熱量ω_be、ω_seの初期値、すなわち時刻ｔ₀における遅れ放熱量ω_be（ｔ₀）、ω_se（ｔ₀）が算出された後に、上述した式（２６）および式（２７）によって時刻ｔ₀以降における遅れ放熱量ω_be（ｔ）、ω_se（ｔ）が算出され、そしてこれら時刻ｔ₀以降の時刻ｔにおける遅れ放熱量ω_be（ｔ）、ω_se（ｔ）に基づいて上記式（２５）を用いて時刻ｔ₀以降の時刻ｔにおける修正出力値ＡＦＭｍｔｅ（ｔ）が算出される。このようにして算出された修正出力値ＡＦＭｍｔｅは図１９に点線で示したように推移する。すなわち、修正出力値ＡＦＭｍｔｅは時刻ｔ₀においてはオブザーバによって算出された近似出力値ＡＦＭｍｔｏと同一の値となっており、その後徐々にエアフロメータモデルＭ４０によって算出された予想出力値ＡＦＭｍｔに近づくように推移する。

このようにして算出された修正出力値ＡＦＭｍｔｅは、上述したようにスロットル弁開度θｔがマップ領域Ａに侵入する前のスロットル通過空気流量の推定値ｍｔの影響が除去された値となっている。このため、マップ領域Ａにおけるエアフロメータ４１の出力値ＡＦＭと修正出力値ＡＦＭｍｔｅとの比は、マップ領域Ａにおける実際のスロットル通過空気流量とスロットルモデルＭ１０によって算出されたスロットル通過空気流量との比を比較的正確に示しており、すなわちマップ領域ＡにおけるスロットルモデルＭ１０の流量係数μｔの誤差率を示している。したがって、上記式（２２）と同様に、エアフロメータ４１の出力値ＡＦＭと修正出力値ＡＦＭｍｔｅとに基づいて下記式（３４）により誤差率Ｅｒｒが算出される。

すなわち、誤差率Ｅｒｒは、スロットル弁開度θｔが特定のマップ領域（図１８、１９に示した例ではマップ領域Ａ）にある間に亘ってエアフロメータ４１の出力値ＡＦＭを積分した値、すなわち時刻ｔ₀からｔ₁に亘って出力値ＡＦＭを積分した値（図１９中の面積Ｓ₃＋Ｓ₄）を、スロットル弁開度θｔが上記特定のマップ領域（図１８、１９に示した例ではマップ領域Ａ）にある間に亘って修正出力値ＡＦＭｍｔｅを積分した値、すなわち修正出力値ＡＦＭｍｔｅを時刻ｔ₀からｔ₁に亘って積分した値（図１９中の面積Ｓ₃）で除算した値（（Ｓ₃＋Ｓ₄）／Ｓ₃）に等しい。このようにして算出された誤差率Ｅｒｒを上記特定のマップ領域（図１８、１９に示した例ではマップ領域Ａ）の流量係数μｔのマップ値に乗算することで、流量係数μｔのマップ値が的確に修正される。

なお、上記実施形態では、スロットル弁開度θｔがθｔ₀〜θｔ₁であるマップ領域Ａに対応する流量係数μｔのマップ値を修正する場合について示したが、他のマップ領域においても同様にそのマップ領域に対応する流量係数μｔのマップ値が修正せしめられる。これにより、内燃機関の過渡運転中であっても各マップ領域に対応する流量係数μｔのマップ値が的確に修正せしめられる。

ところで、上述したように、修正出力値ＡＦＭｍｔｅを用いることでスロットル弁開度θｔがマップ領域Ａに侵入する前のスロットル通過空気流量の推定値ｍｔの影響が除去される。ところが、エアフロメータ４１の出力値ＡＦＭもスロットル弁開度θｔがマップ領域Ａに侵入する前の実際のスロットル通過空気流量の影響を受けており、上述したように修正出力値ＡＦＭｍｔｅを用いたとしてもエアフロメータ４１の出力値ＡＦＭから過去の実際のスロットル通過空気流量の影響を除去することができない。このため、上述したように、エアフロメータ４１の出力値ＡＦＭと修正出力値ＡＦＭｍｔｅとに基づいて算出されたスロットルモデルＭ１０の流量係数μｔの誤差率は、実際の流量係数μｔの誤差率に対して僅かに異なった値となる。

そこで、上記実施形態の変更例では、エアフロメータ４１の出力値ＡＦＭから過去の実際のスロットル通過空気流量の影響を除去することとしている。以下、エアフロメータ４１の出力値ＡＦＭから過去の実際のスロットル通過空気流量の影響を除去してスロットルモデルＭ１０の流量係数μｔの誤差率を正確に算出する方法について説明する。

ところで、一般に、エアフロメータ４１等のセンサの応答は、初期値応答とゼロ応答とに分離することができる。ここで、初期値応答とは、或る時刻ｔ以降の状態量がゼロであると仮定した場合における時刻ｔ以降のセンサの出力値の推移であり、ゼロ応答とは、或る時刻ｔまで状態量がゼロであると仮定した場合における時刻ｔ以降のセンサの出力値の推移である。

図２０は、エアフロメータ４１の出力値ＡＦＭについての初期値応答曲線とゼロ応答曲線とを示す図である。ここで、Φ₁は時刻０を基準とした初期値応答曲線（破線）を示しており、時刻０以降に実際のスロットル通過空気流量がゼロであると仮定した場合におけるエアフロメータ４１の出力値の推移を示している。また、Φ₂は時刻０を基準としたゼロ応答曲線（一点鎖線）を示しており、時刻０までの実際のスロットル通過空気流量がゼロであると仮定した場合、すなわち時刻０からエアフロメータ４１での計測が開始された場合におけるエアフロメータ４１の出力値の推移を示している。

すなわち、初期値応答曲線は時刻０におけるエアフロメータ４１の出力値に応じて、すなわち時刻０よりも前の実際のスロットル通過空気流量に応じて変化するのに対して、ゼロ応答曲線は時刻０以降における実際のスロットル通過空気流量に応じて変化する。そして、時刻０以降の各時刻ｔにおけるエアフロメータ４１の出力値ＡＦＭ（ｔ）は、初期値応答曲線の時刻ｔにおける値Φ₁（ｔ）とゼロ応答曲線の時刻ｔにおける値Φ₂（ｔ）とを合計した値に等しい（ＡＦＭ（ｔ）＝Φ₁（ｔ）＋Φ₂（ｔ））。そして、このようなエアフロメータ４１の出力値ＡＦＭと、初期値応答曲線の各時刻における値Φ₁と、ゼロ応答曲線の各時刻における値Φ₂との関係は、基準となる時刻（図２０に示した例では時刻０）を如何なる時刻にしても成立する。

上述したようにゼロ応答曲線は所定時刻以降における実際のスロットル通過空気流量に応じて変化するため、特定のマップ領域に侵入した時刻を基準としたエアフロメータ４１の出力値ＡＦＭのゼロ応答曲線からは、その特定のマップ領域に侵入する前の実際のスロットル通過空気流量の影響が除去されている。したがって、斯かるゼロ応答曲線を求め、これと修正出力値ＡＦＭｍｔｅについてのゼロ応答曲線（正確には、修正出力値ＡＦＭｍｔｅについてのゼロ応答曲線に対応する曲線）とを後述するように比較することにより、上記特定のマップ領域に侵入する前の実際のスロットル通過空気流量の影響を除去して、流量係数μの誤差率を的確に算出することができる。

ただし、時刻ｔ₀を基準としたゼロ応答を算出するのは困難であることから、本実施形態では、エアフロメータ４１の出力値ＡＦＭから時刻ｔ₀を基準とした初期値応答曲線の各時刻における値Φ₁を減算することでゼロ応答を算出することとしている（Φ₂（ｔ）＝ＡＦＭ（ｔ）−Φ₁（ｔ））。

ここで、時刻ｔ₀を基準とした初期値応答曲線の各時刻における値Φ₁の算出方法について説明する。上述したように初期値応答曲線の各時刻における値Φ₁は、時刻ｔ₀以降に実際のスロットル通過空気流量がゼロであると仮定した場合におけるエアフロメータ４１の出力値の推移を示しているため、上記式（１３）および式（１４）のスロットル通過空気流量の推定値ｍｔをゼロとすることにより求めることができる。すなわち、初期値応答曲線についての放熱量ω_bi、ω_siは下記式（３５）および式（３６）により表される。この場合、時刻ｔ₀における放熱量ω_bi（ｔ₀）、ω_si（ｔ₀）は、時刻ｔ₀における実際の放熱量の近似値である近似放熱量ω_bo（ｔ₀）、ω_so（ｔ₀）が用いられる。

そして、式（３５）、（３６）により算出された放熱量ω_bi、ω_siに基づいて、式（３７）により図１５に示したマップと同様なマップから初期値応答曲線の各時刻における値Φ₁が算出される。なお、式（３５）および式（３６）は、ＥＣＵ３１で処理するために離散化されて用いられる。

このようにして算出された時刻ｔ₀を基準とした初期値応答曲線の各時刻における値Φ₁をエアフロメータ４１の出力値ＡＦＭから減算することにより、エアフロメータ４１の出力値についての時刻ｔ₀を基準としたゼロ応答曲線の各時刻における値Φ₂が算出される。

また、修正出力値ＡＦＭｍｔｅも時刻ｔ₀において近似放熱量ω_bo、ω_soを利用しているため、修正出力値ＡＦＭｍｔｅから上述した初期値応答曲線の各時刻における値Φ₁を減算することにより、修正出力値ＡＦＭｍｔｅについてのゼロ応答曲線の各時刻における値が算出される。このため、本実施形態では、スロットルモデルＭ１０の流量係数μｔの誤差率Ｅｒｒを算出するにあたって、エアフロメータ４１の出力値ＡＦＭと、修正出力値ＡＦＭｍｔｅと、エアフロメータ４１の出力値についての時刻ｔ₀を基準とした初期値応答曲線の各時刻における値Φ₁とを用いて下記式（３８）により算出される。

このように、流量係数μｔの誤差率を算出するにあたって、修正出力値ＡＦＭｍｔｅを用いることにより、エアフロメータモデルＭ４０によって算出された予想出力値ＡＦＭｍｔに存在する過去のスロットル通過空気流量の推定値ｍｔの影響が除去され、すなわち所定時間前のスロットル通過空気流量の影響で予想出力値ＡＦＭｍｔに含まれている成分を予想出力値ＡＦＭｍｔから除去される。さらに、エアフロメータ４１の出力値ＡＦＭおよび修正出力値ＡＦＭｍｔｅから初期値応答曲線の各時刻における値Φ₁を減算することにより、エアフロメータ４１の出力値ＡＦＭ等に存在する過去の実際のスロットル通過空気流量の影響が除去され、すなわち所定時間前のスロットル通過空気流量の影響でエアフロメータの出力値ＡＦＭに含まれている成分をこのエアフロメータの出力値から除去することができる。

このため、上記式（３８）により算出された誤差率Ｅｒｒはスロットル弁開度θｔが所定のマップ領域内にあるときの実際のスロットル通過空気流量およびスロットルモデルＭ１０によって算出されたスロットル通過空気流量の推定値に基づいて算出されることになるため、スロットル弁開度θｔが所定のマップ領域内にあるときの流量係数の誤差率を正確に算出することができる。

ところで、上記説明では、スロットル弁開度θｔが或るマップ領域Ａ（θｔ₀〜θｔ₁）内にあるときにおける流量係数μｔのマップ値を修正するのに、スロットル弁開度θｔがθｔ₀になったとき（時刻ｔ₀）からθｔ₁になったとき（時刻ｔ₁）まで、すなわちスロットル弁開度θｔがそのマップ領域Ａに侵入してから離脱するまでのエアフロメータ４１の出力値ＡＦＭの積分値と修正出力値ＡＦＭｍｔｅの積分値とを用いて誤差率Ｅｒｒを算出することとしている。

ところが、定常運転時等には、スロットル弁開度θｔがほとんど変化せず、よってスロットル弁開度θｔが或る特定のマップ領域内にある時間が長くなる。ここで、例えば、流量係数μｔのマップ値に誤差が生じている場合にはエアモデル（Ｍ１０〜Ｍ３０）によって算出される筒内充填空気量Ｍｃにも誤差が生じているため、流量係数μｔのマップ値に生じている誤差を迅速に補償する必要がある。しかしながら、定常運転時等にはスロットル弁開度θｔが特定のマップ領域からはなかなか離脱しないため、上述したようにスロットル弁開度θｔが特定のマップ領域内に侵入してから離脱するまでの積分値（平均値）を用いて流量係数μｔの誤差率Ｅｒｒを算出するようにしていると、この特定のマップ領域に対応する流量係数μｔのマップ値を迅速に補償することができなくなってしまう。すなわち、流量係数μｔのマップ値に誤差が生じている場合であっても長期間継続している定常運転時等には流量係数μｔを補償することができない。

そこで、本発明の実施形態では、スロットル弁開度θｔが或るマップ領域内に侵入してから後述する最長学習時間が経過してもまだそのマップ領域内にある場合には、スロットル弁開度θｔがこのマップ領域に侵入してから上記最長学習時間が経過するまでのエアフロメータ４１の出力値ＡＦＭの積分値と修正出力値ＡＦＭｍｔｅの積分値とを用いて流量係数μｔの誤差率Ｅｒｒを算出することとしている。以下、このことについて詳細に説明する。

図２１はスロットル弁開度と空気流量に関する各種パラメータの値とのタイムチャートであって図１９と同様な図である。図２１は、スロットル弁開度及び空気流量に関する各種パラメータの値の推移の例を示している。図２１に示した例では、時刻ｔ₀においてスロットル弁開度θｔがθｔ₀となって或る特定のマップ領域Ａ（θｔ₀〜θｔ₁）に侵入していると共に、時刻ｔ₁においてスロットル開度θｔがθｔ₁となってその特定のマップ領域Ａから離脱している。

図示した例では、スロットル弁開度θｔの単位時間当たりの変化量が非常に小さく、よってスロットル弁開度θｔが一つのマップ領域Ａ内にある時間が長い。このため、時刻ｔ₀から時刻ｔ₁までの時間は最長学習時間Ｔｍａｘよりも長いものとなっている。このような場合、本発明の実施形態では、流量係数μｔの誤差率Ｅｒｒの算出を、スロットル弁開度θｔが或る特定のマップ領域Ａから離脱する（時刻ｔ₁）まで待たずに、スロットル弁開度θｔがその特定のマップ領域Ａ内に侵入したとき（時刻ｔ₀）から最長学習時間Ｔｍａｘが経過した時点（時刻ｔ₂）で行うようにしている。

また、最長学習時間Ｔｍａｘが経過した時点で行われる流量係数μｔの誤差率Ｅｒｒの算出は、下記式（３９）によって行われる。すなわち、上述したように、流量係数μｔの誤差率Ｅｒｒの算出にあたっては、エアフロメータ４１の出力値ＡＦＭの積分値および修正出力値ＡＦＭｍｔｅの積分値が用いられるが、これら積分値としてはスロットル弁開度θｔが或る特定のマップ領域に侵入してから離脱するまで積分（すなわち時刻ｔ₀から時刻ｔ₁までの積分）したものではなく、スロットル弁開度θｔがその特定のマップ領域に侵入してから最長学習時間Ｔｍａｘが経過するまで積分（すなわち時刻ｔ₀から時刻ｔ₂まで積分）したものが用いられる。したがって、図２１に示した例では、流量係数μｔの誤差率Ｅｒｒの算出にあたって、エアフロメータ４１の出力値ＡＦＭの積分値として図中の面積Ｓ₅＋Ｓ₆が用いられ、修正出力値ＡＦＭｍｔｅの積分値として図中の面積Ｓ₅が用いられる。

ここで、最長学習時間とは、エアフロメータ４１の出力値ＡＦＭおよび修正出力値ＡＦＭｍｔｅの積分値に基づいて流量係数μｔの誤差率Ｅｒｒを算出するにあたって、高い精度で流量係数μｔの誤差率Ｅｒｒを算出するのに十分な積分時間を意味する。換言すると、最長学習時間とは、高い精度で流量係数μｔの誤差率Ｅｒｒを算出するのに十分な数の出力値ＡＦＭおよび修正出力値ＡＦＭｍｔｅのデータを取得するのにかかる時間を意味する。したがって、スロットル弁開度θｔが特定のマップ領域に侵入してから最長学習時間Ｔｍａｘが経過するまでのエアフロメータ４１の出力値ＡＦＭおよび修正出力値ＡＦＭｍｔｅの積分値に基づいて算出された流量係数μｔの誤差率Ｅｒｒは、それ以上の長時間に亘って積分されたエアフロメータ４１の出力値ＡＦＭおよび修正出力値ＡＦＭｍｔｅの積分値に基づいて算出された流量係数μｔの誤差率Ｅｒｒとほぼ同一の値となる。

また、最長学習時間は、機関運転状態に応じて変更せしめられる。例えば、スロットル弁開度が小さい場合、スロットル弁を通過する空気の速度が音速に近づくため、吸気弁の開閉に伴って吸気管内圧力が変動してもスロットル弁を通過する空気の流速の変動は小さいものとなり、よってスロットル弁通過空気流量の変動も小さいものとなる。このため、エアフロメータ４１の出力値ＡＦＭおよび修正出力値ＡＦＭｍｔｅの積分時間が短くても、これら値の積分値から流量係数μｔの誤差率Ｅｒｒを比較的正確に算出することができる。このため、スロットル弁開度が小さい場合には最長学習時間は短いものとされ、逆にスロットル弁開度が大きい場合には最長学習時間は長いものとされる。

このように、内燃機関の定常運転時等、スロットル弁開度θｔが長時間に亘って或るマップ領域内にある場合に、最長学習時間Ｔｍａｘが経過した時点で流量係数μｔの誤差率Ｅｒｒを算出し、算出された誤差率Ｅｒｒに基づいてそのマップ領域に対応する流量係数μｔのマップ値を修正することにより、スロットル弁開度θｔがそのマップ領域内にある間にそのマップ領域に対応する流量係数μｔのマップ値を修正することができるようになる。このため、流量係数μｔのマップ値を効率的に修正することができる。さらに、最長学習時間が高い精度で流量係数の誤差率を算出するのに十分な積分時間であるため、最長学習時間でエアフロメータ４１の出力値ＡＦＭや修正出力値ＡＦＭｍｔｅの検出、算出を中止しても算出された流量係数μｔの誤差率Ｅｒｒは比較的正確に実際の流量係数の誤差率を示している。したがって、本発明の実施形態によれば、効率的に且つ比較的高い精度で流量係数μｔの誤差を修正することができる。

一方、流量係数μｔの誤差率Ｅｒｒを比較的高い精度で算出するにあたり、エアフロメータ４１の出力値に存在するノイズ等の影響を除去するためにはある程度の数のエアフロメータ４１の出力値ＡＦＭのデータが必要である。また、上述したようにスロットル弁開度θｔが或るマップ領域に侵入したときには、修正出力値ＡＦＭｍｔｅとしてオブザーバによって算出された近似出力値ＡＦＭｍｔｏを用いているため、スロットル弁開度θｔがそのマップ領域に侵入した直後においては修正出力値ＡＦＭｍｔｅの正確性は低い。したがって、流量係数μｔの誤差率Ｅｒｒを比較的高い精度で算出するためには、スロットル弁開度θｔが或るマップ領域内の開度となっていてエアフロメータ４１の出力値ＡＦＭおよび修正出力値ＡＦＭｍｔｅのデータを取得することができる期間がある程度必要である。

そこで、本実施形態では、スロットル弁開度θｔが或るマップ領域に侵入してから最短学習時間が経過せずにそのマップ領域から離脱した場合には流量係数μｔの誤差率Ｅｒｒの修正を行わないこととしている。したがって、上述したように、本発明の実施形態によれば、内燃機関が過渡運転を行っている場合であっても基本的に流量係数μｔのマップ値の修正は行われるが、スロットル弁開度の変化速度があまりにも速い場合には流量係数μｔのマップ値の修正は行われない。

ここで、最短学習時間とは、エアフロメータ４１の出力値ＡＦＭおよび修正出力値ＡＦＭｍｔｅの積分値に基づいて流量係数μｔの誤差率Ｅｒｒを算出するにあたって、ある程度高い精度で流量係数μｔの誤差率Ｅｒｒを算出するのに最低限必要な積分時間を意味する。換言すると、最短学習時間とは、或る程度高い精度で流量係数μｔの誤差率Ｅｒｒを算出するのに最低限必要な数の出力値ＡＦＭおよび修正出力値ＡＦＭｍｔｅのデータを取得するのにかかる時間を意味する。したがって、スロットル弁開度θｔが或るマップ領域に侵入してから最短学習時間Ｔｍｉｎが経過するまでのエアフロメータ４１の出力値ＡＦＭおよび修正出力値ＡＦＭｍｔｅの積分値に基づいて算出された流量係数μｔの誤差率Ｅｒｒは、実際の流量係数μｔの誤差率Ｅｒｒにほぼ等しい。

なお、最短学習時間Ｔｍｉｎも、上述した最長学習時間Ｔｍａｘと同様に、機関運転状態に応じて変更せしめられる。例えば、スロットル弁開度が小さい場合には最短学習時間Ｔｍｉｎは短いものとされ、逆にスロットル弁開度が大きい場合には最短学習時間Ｔｍｉｎは長いものとされる。

このように、最短学習時間Ｔｍｉｎが経過しなかった場合には流量係数μｔの誤差率Ｅｒｒを修正しないようにすることにより、エアフロメータ４１の出力値ＡＦＭに存在するノイズ等の影響や、スロットル弁開度θｔが或るマップ領域に侵入した直後における修正出力値ＡＦＭｍｔｅの正確性の低さの影響を受けて流量係数μｔの誤差率Ｅｒｒの算出精度が低下してしまうのが防止せしめられる。

図２２は、流量係数μｔの修正制御の制御ルーチンを示すフローチャートである。図２２に示した制御ルーチンは、一定時間間隔の割り込みによって実行される。

まず、ステップ１０１において、侵入フラグＸが０であるか否かが判定される。この侵入フラグＸは、スロットル開度θｔが或る特定のマップ領域内にあるか否かを示すフラグであり、スロットル開度θｔがその特定のマップ領域内にある場合に１とされ、それ以外の場合に０とされる。ステップ１０１においてスロットル開度θｔが特定のマップ領域内にない（Ｘ＝０）と判定された場合には、ステップ１０２へと進む。ステップ１０２では、スロットル開度θｔが特定のマップ領域に侵入したか否かが判定される。スロットル開度θｔが特定のマップ領域に侵入していないと判定された場合には制御ルーチンが終了せしめられる。

一方、ステップ１０２において、スロットル開度θｔが特定のマップ領域に侵入したと判定された場合には、ステップ１０３へと進む。ステップ１０３では、そのときのスロットル弁開度θｔ等の機関運転状態に基づいて最長学習時間Ｔｍａｘ及び最短学習時間Ｔｍｉｎが算出され、次いでステップ１０４では時間カウンタＴｃｏｕｎｔがゼロにリセットされる。ここで、時間カウンタＴｃｏｕｎｔは、スロットル弁開度θｔが特定のマップ領域に侵入してから経過した時間を表すカウンタである。そして、ステップ１０５では侵入フラグＸが１にセットされる。

このように侵入フラグＸが１にセットされた後の制御ルーチンでは、まずステップ１０１においてスロットル弁開度θｔが特定のマップ領域内にない（Ｘ＝１）と判定され、ステップ１０６へと進む。ステップ１０６では、時間カウンタＴｃｏｕｎｔの値がＴｃｏｕｎｔ＋１とされる。次いで、ステップ１０７ではエアフロメータ４１の出力値ＡＦＭが検出されると共に、修正出力値ＡＦＭｍｔｅが算出される。

ステップ１０８では、時間カウンタＴｃｏｕｎｔが最長学習時間Ｔｍａｘよりも大きいか否か、すなわちスロットル弁開度θｔが特定のマップ領域に侵入してから最長学習時間Ｔｍａｘが経過したか否かが判定される。スロットル弁開度θｔが特定のマップ領域に侵入してから最長学習時間Ｔｍａｘが経過していない（Ｔｃｏｕｎｔ≦Ｔｍａｘ）と判定された場合にはステップ１０９へと進む。ステップ１０９では、スロットル弁開度θｔが特定のマップ領域から離脱したか否かが判定される。スロットル弁開度θｔが特定のマップ領域から離脱していないと判定された場合には制御ルーチンが終了せしめられる。

一方、ステップ１０９において、スロットル弁開度θｔが特定のマップ領域から離脱したと判定された場合には、ステップ１１０へと進む。ステップ１１０では、時間カウンタＴｃｏｕｎｔが最短学習時間Ｔｍｉｎ以上であるか否か、すなわちスロットル弁開度θｔが特定のマップ領域に侵入してから最短学習時間Ｔｍｉｎが経過したか否かが判定される。ステップ１１０においてスロットル弁開度θｔが特定のマップ領域に侵入してから最短学習時間Ｔｍｉｎが経過した（Ｔｃｏｕｎｔ≧Ｔｍｉｎ）と判定された場合にはステップ１１１へと進む。ステップ１１１では、スロットル弁開度θｔが特定のマップ領域に侵入してからステップ１０７で算出されてきたエアフロメータ４１の出力値ＡＦＭおよび修正出力値ＡＦＭｍｔｅに基づいて流量係数μｔの誤差率Ｅｒｒが算出され、次いでステップ１１２において侵入フラグＸが０にリセットされる。一方、ステップ１１０において、スロットル弁開度θｔが特定のマップ領域に侵入してから最短学習時間Ｔｍｉｎが経過していない（Ｔｃｏｕｎｔ＜Ｔｍｉｎ）と判定された場合には、流量係数μｔの誤差率Ｅｒｒの算出が行われることなくステップ１１３において侵入フラグＸが０にリセットされる。

一方、ステップ１０８においてスロットル弁開度θｔが特定のマップ領域に侵入してから最長学習時間Ｔｍａｘが経過した（Ｔｃｏｕｎｔ＞Ｔｍａｘ）と判定された場合にはステップ１１１へと進む。ステップ１１１では、スロットル弁開度θｔが特定のマップ領域に侵入してから最長学習時間Ｔｍａｘが経過するまでの間にステップ１０７で算出されてきたエアフロメータ４１の出力値ＡＦＭおよび修正出力値ＡＦＭｍｔｅに基づいて流量係数μｔの誤差率Ｅｒｒが算出され、次いでステップ１１２において侵入フラグＸが０にリセットされる。

ところで、上記実施形態では、スロットル開度θｔを複数のマップ領域に分けて、スロットル開度θｔの各マップ領域に対応する流量係数μｔを設定するようにしているが、実際には各スロットル開度θｔの値に対して最適な流量係数μｔの値が存在する。したがって、スロットル開度θｔの或るマップ領域について考えてみると、そのマップ領域内の複数のスロットル開度θｔの値それぞれに対して最適な流量係数μｔの値が存在することになる。したがって、各マップ領域内のスロットル弁開度θｔに対して最適な流量係数μｔは、一つの値に定まるものではなく、或る程度の幅をもったものとなる。

このような状況を考慮すると、或るマップ領域に対応する流量係数として一つの値を設定する場合には、そのマップ領域に対応する流量係数のマップ値としてこのマップ領域の中心のスロットル開度θｔにおいて最適な流量係数μｔの値を用いることが好ましい。このように各マップ領域に対応する流量係数μｔのマップ値を設定することにより、流量係数μｔのマップ値をそのマップ領域内の如何なるスロットル弁開度θｔに対しても比較的適した値とすることができる。

そこで、本発明の第二実施形態では、スロットル弁開度θｔが或るマップ領域に侵入してから上記最長学習時間Ｔｍａｘよりも長い時間の経過後にスロットル弁開度θｔがそのマップ領域から離脱した場合、スロットル弁開度θｔが下記の中心領域にあったときのエアフロメータ４１の出力値ＡＦＭおよび修正出力値ＡＦＭｍｔｅに基づいて流量係数μｔの誤差率Ｅｒｒを算出するようにしている。以下、このことについて詳細に説明する。

図２３はスロットル弁開度と空気流量に関する各種パラメータの値とのタイムチャートであって、図１９および図２１と同様な図である。図２３に示した例では、スロットル弁開度および空気流量に関する各種パラメータの値は図２１に示した例と同様に推移している。また、図２３に示した例では、図２１と同様に、時刻ｔ₀においてスロットル弁開度θｔがθｔ₀となって或る特定のマップ領域Ａ（θｔ₀〜θｔ₁）に侵入していると共に、時刻ｔ₁においてスロットル弁開度θｔがθｔ₁となってその特定のマップ領域Ａから離脱している。

図２３から分かるように、図示した例では、時刻ｔ₀から時刻ｔ₁までの時間は最長学習時間Ｔｍａｘよりも長いものとなっている。このような場合、本実施形態では、流量係数μｔの誤差率Ｅｒｒの算出を、スロットル弁開度θｔが中心領域Ａ’にあったときのエアフロメータ４１の出力値ＡＦＭおよび修正出力値ＡＦＭｍｔｅに基づいて行うこととしている。

ここで、中心領域とは、或るマップ領域に含まれるスロットル弁開度の領域であって、そのマップ領域の中心近傍の領域を意味する。図２３に示した例で説明すると、マップ領域Ａの中心領域Ａ’の範囲（θｔ₃〜θｔ₄）はそのマップ領域Ａの範囲（θｔ₀〜θｔ₁）に含まれている。また、この中心領域Ａ’は、マップ領域Ａの中心、すなわちスロットル弁開度θｔ₀およびθｔ₁の中間のスロットル弁開度θｔ₂（＝（θｔ₀＋θｔ₁）／２）の近傍の領域となっており、より詳細にはマップ領域Ａの中心のスロットル弁開度θｔ₂を中心とした領域となっている。

また、本実施形態では、中心領域Ａ’は、スロットル弁開度θｔが中心領域Ａ’内にある期間が最長学習時間Ｔｍａｘと等しくなるように定められる。図２３に示した例で説明すると、中心領域Ａ’の下限スロットル弁開度θｔ₃と上限スロットル弁開度θｔ₄は、スロットル弁開度θｔが下限スロットル弁開度θｔ₃と上限スロットル弁開度θｔ₄との間にある期間（時刻ｔ₃〜時刻ｔ₄）が最長学習時間Ｔｍａｘと等しくなるように定められる。なお、図示した例では、スロットル弁開度θｔが中心領域Ａ’内にある期間は時刻ｔ₃から時刻ｔ₄まで連続しているが、斯かる期間は連続していない複数の期間から構成されてもよい。

したがって、図示した例では、流量係数μｔの誤差率Ｅｒｒの算出は、下記式（４０）によって行われる。すなわち、上述したように、流量係数μｔの誤差率Ｅｒｒの算出にあたっては、エアフロメータ４１の出力値ＡＦＭの積分値および修正出力値ＡＦＭｍｔｅの積分値が用いられるが、本実施形態ではこれら積分値としてスロットル弁開度θｔが下限スロットル弁開度θｔ₃となった時刻ｔ₃からスロットル弁開度θｔが上限スロットル弁開度θｔ₄となった時刻ｔ₄まで積分した値が用いられる。したがって、図２３に示した例では、流量係数μｔの誤差率Ｅｒｒの算出にあたっては、エアフロメータ４１の出力値ＡＦＭの積分値として図中の面積Ｓ₇＋Ｓ₈が用いられ、修正出力値ＡＦＭｍｔｅの積分値として図中の面積Ｓ₇が用いられる。

なお、上記第一実施形態では、修正出力値ＡＦＭｍｔｅの初期値をスロットル開度θｔがマップ領域Ａに侵入したとき（すなわち、時刻ｔ₀）においてオブザーバによって算出されたボビン部４１ａ３の近似放熱量ω_boおよびサポート部４１ａ４の近似放熱量ω_soに基づいて算出すると共に、その後式（２５）、式（２６）および式（２７）を用いて時刻ｔ₀以降の修正出力値ＡＦＭｍｔｅを算出するようにしている。しかしながら、このように修正出力値ＡＦＭｍｔｅを算出すると、修正出力値ＡＦＭｍｔｅには時刻ｔ₀から時刻ｔ₃までの間のスロットル通過空気流量の推定値ｍｔの影響が含まれた値となってしまう。このため、このように修正出力値ＡＦＭｍｔｅを算出すると、スロットル弁開度θｔが中心領域Ａ’にある場合における流量係数μｔの誤差率Ｅｒｒを的確に算出することができない。

そこで、本実施形態では、図２３に示したように修正出力値ＡＦＭｍｔｅの初期値をスロットル開度θｔが中心領域Ａ’に侵入したとき（すなわち、時刻ｔ₃）においてオブザーバによって算出されたボビン部４１ａ３の近似放熱量ω_boおよびサポート部４１ａ４の近似放熱量ω_soに基づいて算出すると共に、その後式（２５）、式（２６）および式（２７）を用いて時刻ｔ₃以降の修正出力値ＡＦＭｍｔｅを算出するようにしている（このように算出した修正出力値をＡＦＭｍｔｅ’とする）。これにより修正出力値ＡＦＭｍｔｅ’には時刻ｔ₃より前のスロットル通過空気流量の推定値ｍｔの影響が含まれておらず、スロットル開度θｔが中心領域Ａ’にある場合における流量係数μｔの誤差率Ｅｒｒを的確に算出することができるようになる。

このように、内燃機関の定常運転時等、スロットル弁開度θｔが長時間に亘って或る特定のマップ領域内にある場合に、スロットル弁開度θｔが中心領域にあったときのエアフロメータ４１の出力値ＡＦＭおよび修正出力値ＡＦＭｍｔｅに基づいて流量係数μｔの誤差率Ｅｒｒを算出するようにしているため、その特定のマップ領域の中心のスロットル弁開度に対応するように誤差率Ｅｒｒを算出することができる。これにより、この誤差率Ｅｒｒに基づいて修正された流量係数μｔのマップ値は、或る特定のマップ領域の中心のスロットル弁開度に対応するような値となっており、これにより流量係数μｔのマップ値をそのマップ領域内の如何なるスロットル弁開度θｔに対しても比較的適した値とすることができる。

図２４は、流量係数μｔの修正制御の制御ルーチンを示すフローチャートである。図２４に示した制御ルーチンは、一定時間間隔の割り込みによって実行される。ステップ１５１〜ステップ１５７は、図２２のステップ１０１〜ステップ１０７と同様であるため説明を省略する。

ステップ１５８では、スロットル弁開度θｔが特定のマップ領域から離脱したか否かが判定される。スロットル弁開度θｔが特定のマップ領域から離脱していないと判定された場合には制御ルーチンが終了せしめられる。一方、ステップ１５８において、スロットル弁開度θｔが特定のマップ領域から離脱したと判定された場合には、ステップ１５９へと進む。ステップ１５９では、時間カウンタＴｃｏｕｎｔが最短学習時間Ｔｍｉｎ以上であるか否か、すなわちスロットル弁開度θｔが特定のマップ領域に侵入してから最短学習時間Ｔｍｉｎが経過しているか否かが判定される。ステップ１５９においてスロットル弁開度θｔが特定のマップ領域に侵入してから最短学習時間Ｔｍｉｎが経過していない（Ｔｃｏｕｎｔ＜Ｔｍｉｎ）と判定された場合には、流量係数μｔの誤差率Ｅｒｒの算出が行われることなくステップ１６０において侵入フラグＸが０にリセットされる。

一方、ステップ１５９においてスロットル弁開度θｔが特定のマップ領域に侵入してから最短学習時間Ｔｍｉｎが経過している（Ｔｃｏｕｎｔ≧Ｔｍｉｎ）と判定された場合にはステップ１６１へと進む。ステップ１６１では、スロットル弁開度θｔが特定のマップ領域に侵入してから最長学習時間Ｔｍａｘが経過しているか否かが判定される。スロットル弁開度θｔが特定のマップ領域に侵入してから最長学習時間Ｔｍａｘが経過していると判定された場合（Ｔｃｏｕｎｔ＞Ｔｍａｘ）には、ステップ１６２へと進む。

ステップ１６２では、スロットル弁開度θｔが中心領域内にあったときのエアフロメータ４１の出力値ＡＦＭおよびスロットル通過空気流量の推定値ｍｔが抽出される。具体的には、スロットル弁開度θｔがマップ領域Ａの中心のスロットル弁開度（図２３に示した例ではθｔ₂）に近いときの出力値ＡＦＭおよびスロットル通過空気流量の推定値ｍｔから抽出していき、出力値ＡＦＭおよびスロットル通過空気流量の推定値ｍｔのデータの数が最大学習時間Ｔｍａｘに対応する数になるまで抽出が行われる。次いで、ステップ１６３では、ステップ１６２で抽出されたスロットル通過空気流量の推定値ｍｔと抽出されたデータのうち最も早く抽出されたデータの抽出時刻におけるオブザーバの出力とに基づいてスロットル開度θｔが中心領域Ａ’に侵入してからの修正出力値ＡＦＭｍｔｅ’が算出される。次いでステップ１６４では、ステップ１６２において抽出されたエアフロメータ４１の出力値ＡＦＭおよびステップ１６３で算出された修正出力値ＡＦＭｍｔｅ’に基づいて式（４０）により流量係数μｔの誤差率Ｅｒｒが算出され、その後ステップ１６０において侵入フラグＸが０にリセットされる。

一方、ステップ１６１において、スロットル弁開度θｔが特定のマップ領域に侵入してから最長学習時間Ｔｍａｘが経過していないと判定された場合（Ｔｃｏｕｎｔ≦Ｔｍａｘ）には、ステップ１６５へと進む。ステップ１６５では、ステップ１５７で算出されてきたエアフロメータ４１の出力値ＡＦＭおよび修正出力値ＡＦＭｍｔｅの全てに基づいて流量係数μｔの誤差率Ｅｒｒが算出され、その後ステップ１６０において侵入フラグＸが０にリセットされる。

次に、吸気弁モデルＭ３０の修正操作について説明する。上述したように、吸気弁モデルＭ３０の式（７）において用いられる値ａは比例係数であり、また値ｂは燃焼室５内に残存していた既燃ガスを表す値である。これら値ａ、ｂは、機関回転数ＮＥ、吸気弁の位相角および作用角（以下では、吸気弁の開弁時期ＶＶＴを例にとって説明する）に応じて変化する。このため、機関回転数ＮＥおよび吸気弁の開弁時期ＶＶＴとこれら値ａ、ｂとの関係を予め実験的にまたは計算により算出し、図２５（ａ）および図２５（ｂ）に示したようにマップとしてＥＣＵ３１に保存すると共に、吸気弁モデルＭ３０において上記式（７）の計算をするにあたってはクランク角センサ４７によって検出された機関回転数ＮＥおよび吸気弁制御装置１３への吸気弁開弁時期ＶＶＴを用いて上記マップによりこれら値ａ、ｂを算出している。なお、図２５（ａ）は値ａのマップを、図２５（ｂ）は値ｂのマップをそれぞれ示しており、図中のｘ軸は機関回転数ＮＥを、ｙ軸は吸気弁開弁時期ＶＶＴ（ｙ軸が大きくなるほど進角）をそれぞれ示している。

ところが、これら値ａ、ｂには、上述したように使用による特性の変化や製作公差によるばらつきなどにより誤差が生じる場合がある。このように、値ａ、ｂに誤差が生じている場合には、吸気弁モデルＭ３０のモデル化誤差を補償するためにこれら値ａ、ｂのマップ値を修正していく必要がある。ここで、本実施形態では、値ａ、ｂは、機関運転領域毎に定められている。図２５に示した例では、値ａ、ｂは、マップ領域（例えば、図２５中の斜線部分が一つのマップ領域を表している）毎に、すなわち一定間隔の機関回転数毎および一定間隔の吸気弁開弁時期毎に定められている。したがって値ａ、ｂの修正は、各マップ領域毎に行われなければならない。

そこで、本実施形態では、吸気弁モデルＭ３０の式（１１）で用いられる値ａ、ｂをマップ領域毎に修正することとしている。以下では、例として、機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴがマップ領域Ｂ（すなわち、機関回転数ＮＥがＮＥ₁からＮＥ₂の間であって且つ吸気弁開弁時期ＶＶＴがＶＶＴ₁からＶＶＴ₂の間であるマップ領域）にある場合における値ａ、ｂのマップ値ａ_ij、ｂ_ijの修正方法について説明する。

図２６は、内燃機関の運転中における吸気管内圧力Ｐｍ、機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴのタイムチャートである。図示した例では、時刻ｔ₀以降において機関回転数ＮＥがＮＥ₀以上となっており、また時刻ｔ₁以降において機関回転数ＮＥがＮＥ₁以上となっている。すなわち、機関回転数ＮＥは時刻ｔ₀から時刻ｔ₁までに亘ってＮＥ₀とＮＥ₁との間の値となっており、それ以外の時刻においてはＮＥ₀とＮＥ₁との間の値となっていない。一方、吸気弁開弁時期ＶＶＴは時刻ｔ₀から時刻ｔ₁までに亘ってＶＶＴ₁とＶＶＴ₂との間の値となっている。

したがって、機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴの値は、時刻ｔ₀以前はマップ領域Ｂ内にはなく、時刻ｔ₀においてマップ領域Ｂに侵入すると共に時刻ｔ₀から時刻ｔ₁までに亘ってマップ領域Ｂ内にある。その後、時刻ｔ₁においてマップ領域Ｂから離脱して、時刻ｔ₁以降はマップ領域Ｂ内にはない。そして、機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴの値がマップ領域Ｂ内にあるとき、すなわち時刻ｔ₀から時刻ｔ₁までの間において、吸気弁モデルＭ３０ではマップ領域Ｂに対応した値ａ_ij、ｂ_ijが用いられる。

ここで、上記式（６）について再び考えると、式（６）の左辺は吸気管部分２３の気体の内部エネルギの時間的変化量を表しており、式（６）の右辺のＣｐ・ｍｔ・Ｔａは吸気管部分２３に流入する気体のエネルギを、Ｃｐ・ｍｃ・Ｔｍは吸気管部分２３から流出した気体のエネルギをそれぞれ表している。ここで、吸気管部分２３における気体の状態方程式（Ｐｍ・Ｖｍ＝Ｍ・Ｒ・Ｔｍ）を考慮すると、式（６）の左辺は下記式（４１）のように表せる。

したがって、時刻ｔ₀から時刻ｔ₁までの間の内部エネルギの変化量は下記式（４２）のように表せる。特に、吸気管内圧力センサ４０により実際の吸気管内圧力Ｐｍｒが検出されているため、式（４２）において実際の吸気管内圧力Ｐｍｒを用いることによって、時刻ｔ₀から時刻ｔ₁までの間の実際の内部エネルギの変化量を正確に算出することができる。なお、式（４２）においてＰｍｒ（ｔ₀）は時刻ｔ₀における実際の吸気管内圧力を、Ｐｍｒ（ｔ₁）は時刻ｔ₁における実際の吸気管内圧力をそれぞれ示している。

一方、吸気管部分２３に流入する気体のエネルギ（すなわち、式（６）の左辺のＣｐ・ｍｔ・Ｔａ）に関しては、上述したようにスロットルモデルＭ１０の流量係数μｔのマップ値が逐次修正されるため、スロットルモデルＭ１０の流量係数μｔは適切な値となっている。したがって、スロットルモデルＭ１０の式（１）の吸気管内圧力Ｐｍとして吸気管内圧力センサ４０によって検出された実際の吸気管内圧力Ｐｍｒを用いると、吸気管部分２３に流入する気体のエネルギを正確に算出することができる。特に、時刻ｔ₀から時刻ｔ₁までに吸気管部分２３に流入する気体のエネルギは下記式（４３）により算出することができる。

このように時刻ｔ₀から時刻ｔ₁までの間の内部エネルギの変化量が上記式（４２）で、時刻ｔ₀から時刻ｔ₁までに吸気管部分２３に流入する気体のエネルギが上記式（４３）でそれぞれ正確に算出されることから、式（４３）から式（４２）を減算した下記式（４４）により時刻ｔ₀から時刻ｔ₁までに吸気管部分２３から流出した気体のエネルギを正確に算出することができる。

ここで、上述したように吸気管部分２３から流出した気体のエネルギをＣｐ・ｍｃ・Ｔｍで表すことができ、またこのｍｃとして吸気弁モデルＭ３０によって算出された筒内吸入空気流量ｍｃ（すなわち、式（７））を用いると、時刻ｔ₀から時刻ｔ₁までに吸気管部分２３から流出した気体のエネルギは下記式（４５）のように表せる。そして、式（４４）及び式（４５）はいずれも吸気管部分２３から流出した気体のエネルギを表しているため、本来であればこれら式によって算出された値は一致するはずである。

ところが、吸気弁モデルＭ３０のモデル式（７）に誤差があると、式（４４）と式（４５）とによって算出された値は一致しない。特に上記式（４５）では吸気管内圧力として吸気管内圧力センサ４０により検出された実際の吸気管内圧力Ｐｍｒが用いられていることを考慮すると、式（４４）と式（４５）とが一致しない場合は式（７）で用いられる値ａ、ｂのマップ値に誤差がある場合であると考えられる。

そこで、本実施形態では、実際の筒内吸入空気流量に対する吸気弁モデルＭ３０によって算出された筒内吸入空気流量ｍｃの誤差率（実際の筒内吸入空気流量／吸気弁モデルＭ３０によって算出された筒内吸入空気流量）をαとする。この誤差率αを式（４５）に乗算して算出された値と式（４４）によって算出された値は等しくなるため、下記式（４６）が成立する。そして、式（４６）を誤差率αについて解くと、誤差率αは下記式（４７）のように表せる。

このようにして算出された誤差率αは、機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴがマップ領域Ｂにある場合における実際の筒内吸入空気流量と吸気弁モデルＭ３０によって算出された筒内吸入空気流量の推定値ｍｃとの誤差率を表している。そして、この筒内吸入空気流量における誤差を補償するためには、本実施形態では、マップ領域Ｂに対応するマップ値ａ_ij、ｂ_ijに誤差率αを乗算したものを新たな値ａ、ｂのマップ領域Ｂにおけるマップ値ａ_ij、ｂ_ijとすることとしている（ａ_ij＝α・ａ_ij、ｂ_ij＝α・ｂ_ij）。このように吸気弁モデルＭ３０の式（７）に用いられる値ａ、ｂのマップ値を修正することにより、吸気弁モデルＭ３０によって算出される筒内吸入空気流量の推定値ｍｃを実際の筒内吸入空気流量とほぼ等しい値とすることができる。

このように、本実施形態では、機関運転状態（例えば、機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴ等）が特定の機関運転領域（例えば、上記マップ領域Ｂ）に侵入してから離脱するまでに吸気管部分２３から流出した気体のエネルギの実際の値（すなわち、吸気管部分２３のエネルギ保存則を用いて算出された値）と、斯かる気体のエネルギの吸気弁モデルＭ３０を用いて算出された値とが求められ、実際の値と吸気弁モデルＭ３０を用いて算出された値との比率に基づいて上記特定の機関運転領域に対応する値ａ、ｂのマップ値を修正することにより、内燃機関が過渡運転を行っている場合であっても値ａ、ｂのマップ値を的確に修正することができる。

ここで、吸気管モデルＭ２０によって算出された吸気管内圧力の推定値Ｐｍは所定時間前に吸気弁モデルＭ３０によって算出された筒内吸入空気流量の推定値ｍｃを用いて算出されている。このため、例えば吸気管内圧力の実測値と吸気弁モデルＭ３０によって算出された筒内吸入空気流量の推定値ｍｃを用いて吸気管モデルＭ２０によって算出された吸気管内圧力の推定値との偏差に基づいて値ａ、ｂのマップ値を修正しようとすると、吸気管内圧力の実測値が現在の実際の筒内吸入空気流量に対応した値であるのに対して吸気管モデルＭ２０によって算出された吸気管内圧力の推定値は吸気弁モデルＭ３０によって算出された過去の筒内吸入空気流量の推定値ｍｃに対応した値となっているため、内燃機関が定常運転を行っているときには比較的的確に吸気弁モデルを修正することができるが、内燃機関が過渡運転を行っているときには的確に吸気弁モデルを修正することができない。

一方、本実施形態では、機関運転状態（例えば、機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴ等）が特定の機関運転領域（例えば、上記マップ領域Ｂ）に侵入してから離脱するまでに吸気管部分２３から流出した気体のエネルギの実際の値（すなわち、吸気管部分２３のエネルギ保存則を用いて算出された値）と、斯かる気体のエネルギの吸気弁モデルＭ３０を用いて算出された値とを比較し、これら値が一致するように値ａ、ｂのマップ値を修正するようにしている。このように機関運転状態が特定の機関運転領域にある間に吸気管部分２３から流出した気体のエネルギの総量同士を比較することにより、すなわち斯かる間の単位時間当たりのエネルギ流出量の積分値同士を比較することにより、内燃機関が過渡運転を行っている場合であっても、この特定の機関運転領域に対応する値ａ、ｂのマップ値を的確に修正することができるようになる。

次に、上記式（４７）を内燃機関の制御装置に実装して、実際にスロットル通過空気流量の推定値ｍｔおよび吸気管内圧力の実測値Ｐｍｒから誤差率αを算出する場合について説明する。ＥＣＵ３１では式（４７）中の積分項をそのまま算出することはできない。このため、式（４７）中のスロットル通過空気流量の推定値ｍｔについての積分項は時刻ｔ、計算間隔Δｔを用いて離散化された下記式（４８）によって算出され、吸気管内圧力の実測値Ｐｍｒについての積分項は時刻ｔ、計算間隔Δｔを用いて離散化された下記式（４９）によって算出される。

そして、スロットル通過空気流量の推定値ｍｔについての積分項に関しては、時刻ｔにおけるスロットル通過空気流量の推定値ｍｔ（ｔ）が式（４８）に代入され、これにより時刻ｔにおけるスロットル通過空気流量の推定値ｍｔについての積分項の値ｍｔｉｎｔ（ｔ）が算出される。また、吸気管内圧力の実測値Ｐｍｒについての積分項に関しては、時刻ｔにおける吸気管内圧力の実測値Ｐｍｒ（ｔ）が式（４９）に代入され、これにより時刻ｔにおける吸気管内圧力の実測値Ｐｍｒについての積分項の値ｈｉｎｔ（ｔ）が算出される。なお、いずれの積分項の算出にあたっても、機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴが対象となるマップ領域に侵入したとき、すなわち時刻ｔ₀におけるスロットル通過空気流量の推定値ｍｔ（ｔ₀）および吸気管内圧力Ｐｍ（ｔ₀）を初期値として上記式（４８）および式（４９）の計算が行われる。

そして、式（４７）による誤差率αの算出は、機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴが対象となるマップ領域から離脱した後、すなわち時刻ｔ₁以降に行われる。式（４８）についての初期値が時刻ｔ₀におけるスロットル通過空気流量の推定値ｍｔ（ｔ₀）であることから、時刻ｔ₁におけるスロットル通過空気流量の推定値ｍｔについての積分項の値ｍｔｉｎｔ（ｔ₁）は時刻ｔ₀から時刻ｔ₁までのｍｔの積分値を表しており、また、式（４９）についての初期値が時刻ｔ₀における吸気管内圧力の実測値Ｐｍｒ（ｔ₀）であることから、時刻ｔ₁における吸気管内圧力の実測値Ｐｍｒについての積分項の値ｈｉｎｔ（ｔ₁）は時刻ｔ₀から時刻ｔ₁までのａ・Ｐｍｒ＋ｂの積分値を表している。このため、上記マップ領域における誤差率αは下記式（５０）により算出せしめられる。

ところで、上記説明では、吸気管内圧力センサ４０によって実測された吸気管内圧力を用いて式（４３）により誤差率αを算出することとしている。ところが、実際の吸気管部分２３内の圧力は吸気脈動により大きく変動しており、吸気管内圧力センサ４０によって実測される吸気管内圧力も図２７に実線で示したように大きく変動する。一方、上述した吸気管モデルＭ２０や吸気弁モデルＭ３０では、吸気管部分２３内に生じる圧力脈動を無視して吸気管部分２３内の実際の圧力の１サイクル（本実施形態では、クランク角７２０°）当たりの平均値（以下、「サイクル間平均値」と称す。）を吸気管内圧力Ｐｍとして用いている。また、吸気管内圧力センサ４０の出力にはノイズが含まれてしまうことも多い。

このため、吸気管内圧力センサ４０によって実測された時刻ｔ₀および時刻ｔ₁における吸気管内圧力Ｐｍ（ｔ₀）およびＰｍ（ｔ₁）を直接用いて誤差率αを算出すると、すなわち機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴが或るマップ領域に侵入した時およびこのマップ領域から離脱した時に吸気管内圧力センサ４０によって実測された吸気管内圧力Ｐｍｒを直接用いて誤差率αを算出すると、算出された誤差率αの値が不正確なものになってしまう場合がある。

そこで、本実施形態では、上述したように吸気管内圧力センサ４０によって実測された吸気管内圧力を直接用いて誤差率αを算出せずに、吸気管内圧力センサ４０によって実測された吸気管内圧力のサイクル間平均値（図２７中の破線）を用いて誤差率αを算出することとしている。すなわち、機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴが或るマップ領域に侵入した時を含むサイクルにおいて吸気管内圧力センサ４０によって実測された吸気管内圧力のサイクル間平均値をこのマップ侵入時の吸気管内圧力Ｐｍｒとして用い、このマップ領域から離脱した時を含むサイクルにおいて吸気管内圧力センサ４０によって実測された吸気管内圧力のサイクル間平均値をこのマップ離脱時の吸気管内圧力Ｐｍｒとして用いて、上記式（４３）により誤差率αを算出する。

このように、式（４３）に用いる吸気管内圧力Ｐｍｒとして吸気管内圧力センサ４０によって検出された吸気管内圧力のサイクル間平均値を用いることにより、式（４３）に入力される吸気管内圧力Ｐｍｒが、吸気管モデルＭ２０や吸気弁モデルＭ３０で用いられている吸気管内圧力Ｐｍと同様に、吸気管部分２３内に生じる圧力脈動の影響を排除した値となっており、よって式（４３）によって正確に誤差率αを算出することができるようになる。また、誤差率αを算出するにあたって吸気管内圧力センサ４０の出力に含まれているノイズの影響を抑制することができる。

次に、上記実施形態の変更例について説明する。上記実施形態では、誤差率αを算出した後、この誤差率αを値ａおよび値ｂに乗算して値ａ、ｂの両マップ値を同じ割合で修正している。しかしながら、値ａおよび値ｂの両マップ値に必ずしも同程度の誤差が生じている場合ばかりではなく、値ａのみの誤差が大きい場合や値ｂのみの誤差が大きい場合も存在する。したがって、値ａおよび値ｂのマップ値を正確に修正するためには、値ａおよび値ｂを同率で修正せずに、値ａおよび値ｂを別々に修正する必要がある。

そこで、本変更例では、上述したように算出された誤差率αに基づいて値ａおよび値ｂを別々に修正することとしている。本変更例では、まず、機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴが或るマップ領域にある場合において、上記式（４３）に基づいて誤差率αを算出すると共に、誤差率αの算出における積分期間、すなわち機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴが上記マップ領域にある期間中における吸気管内圧力Ｐｍの平均値Ｐｍａｖｅを算出する。便宜上、このようにして算出された誤差率および吸気管内圧力の平均値をそれぞれα₁、Ｐｍ₁とする。その後、再び機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴが上記マップ領域にある場合において、上述したように誤差率αと吸気管内圧力の平均値Ｐｍａｖｅを算出する。便宜上、このようにして算出された誤差率および吸気管内圧力の平均値をそれぞれα₂、Ｐｍ₂とする。

図２８は、機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴが上記マップ領域にある場合における吸気管内圧力Ｐｍと筒内吸入空気流量ｍｃとの関係を示した図である。値ａ、ｂのマップ値が修正される前に吸気弁モデルＭ３０によって算出される筒内吸入空気流量ｍｃと吸気管内圧力Ｐｍとの関係は、例えば、図２８に実線で示された関係になっている。

ここで、上述したようにして算出された誤差率α₁および吸気管内圧力Ｐｍ₁を用いて算出されたα₁（ａ・Ｐｍ₁＋ｂ）は吸気管内圧力がＰｍ₁であるときの実際の筒内吸入空気流量を正確に表しており、また、誤差率α₂および吸気管内圧力Ｐｍ₂を用いて算出されたα₂（ａ・Ｐｍ₂＋ｂ）は吸気管内圧力がＰｍ₂であるときの実際の筒内吸入空気流量を正確に表している。このことから、点（Ｐｍ₁、α₁（ａ・Ｐｍ₁＋ｂ））および点（Ｐｍ₂、α₂（ａ・Ｐｍ₂＋ｂ））の二点を通る直線（図２８中の破線）は、吸気管内圧力と実際の筒内吸入空気流量との関係を比較的正確に表していると考えられる。

したがって、値ａ、ｂのマップ値を修正することで吸気管内圧力Ｐｍと吸気弁モデルＭ３０によって算出される筒内吸入空気流量ｍｃとの関係を、上記二点を通る直線とすることで、吸気弁モデルＭ３０によって算出される筒内吸入空気流量の推定値ｍｃを実際の筒内吸入空気流量とほぼ等しい値にすることができる。

そこで、本実施形態では、上述した誤差率α₁、α₂および吸気管内圧力Ｐｍ₁、Ｐｍ₂を用いることで、値ａ、ｂのマップ値を下記式（５１）、（５２）により修正することとしている。なお、式（５１）、（５２）において、値ａ’および値ｂ’はそれぞれ修正後のマップ値を示している。

このように値ａ、ｂのマップ値を修正することで、修正後の値ａ’および値ｂ’のマップ値を用いた吸気弁モデルＭ３０によって算出される筒内吸入空気流量ｍｃは実際の筒内空乳空気流量とほぼ等しい値となっている。したがって、本変更例によれば、値ａおよび値ｂのマップ値を別々に修正することができ、これにより値ａおよび値ｂのマップ値を正確に修正することができるようになる。

また、吸気弁モデルＭ３０の修正操作においても、スロットルモデルＭ１０における流量係数の修正と同様に、機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴが或るマップ領域内に侵入してから上述したような最長学習時間が経過してもまだそのマップ領域内にある場合には、機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴがこのマップ領域に侵入してから最長学習時間が経過するまでの吸気管内圧力の実測値Ｐｍｒを用いて値ａおよび値ｂの誤差率αを算出するようにしている。これにより、値ａおよび値ｂのマップ値の修正を効率的に且つ比較的高い精度で行うことができる。

なお、吸気弁モデルＭ３０の修正操作において用いられる最長学習時間Ｔｍａｘも、スロットルモデルＭ１０の修正操作の場合と同様に機関運転状態に応じて変更せしめられる。例えば、上述したように吸気管モデルＭ２０や吸気弁モデルＭ３０においては、パラメータとして吸気管内圧力Ｐｍの実際の値ではなくそのサイクル間平均値が用いられる。したがって、吸気管内圧力の実測値Ｐｍｒを用いて吸気弁モデルＭ３０の値ａおよび値ｂの誤差率αを算出する場合には、吸気管内圧力の実測値Ｐｍｒの積分値を算出するにあたってサイクル数が多いことが必要とされる。よって、機関回転数が高いときには単位時間当たりのサイクル数が多くなるため最長学習時間が短いものとされ、逆に機関回転数が低いときには単位時間当たりのサイクル数が少なくなるため最長学習時間が長いものとされる。

さらに、吸気弁モデルＭ３０の修正操作においても、スロットルモデルＭ１０における流量係数の修正と同様に、機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴが或るマップ領域内に侵入してから最短学習時間が経過せずにそのマップ領域から離脱した場合には値ａおよび値ｂのマップ値の修正を行わないこととしている。これにより、吸気管内圧力センサ４０の実測値に存在するノイズ等の影響を受けて値ａおよび値ｂの誤差率αの算出精度が低下してしまうのが防止される。

或いは、機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴが或るマップ領域に侵入してから上記最長学習時間Ｔｍａｘよりも長い時間の経過後にそのマップ領域から離脱した場合に、機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴが中心領域にあったときの吸気管内圧力センサ４０の実測値に基づいて値ａおよび値ｂのマップ値を修正するようにしてもよい。この場合、中心領域Ｂ’は、上述したスロットルモデルＭ１０の修正操作の場合と同様に、機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴが中心領域Ｂ’内にある期間が最長学習時間Ｔｍａｘと等しくなるように定められる。図２９に示した例で説明すると、中心領域Ｂ’の下限機関回転数ＮＥ₃および上限機関回転数ＮＥ₄と下限吸気弁開弁時期ＶＶＴ₃および上限吸気弁開弁時期ＶＶＴ₄とは、機関回転数ＮＥが下限機関回転数ＮＥ₃と上限機関回転数ＮＥ₄との間にあって且つ吸気弁開弁時期ＶＶＴが下限吸気弁開弁時期ＶＶＴ₃と上限吸気弁開弁時期ＶＶＴ₄との間にある期間（時刻ｔ₃〜時刻ｔ₄）が最長学習時間Ｔｍａｘと等しくなるように定められる。

したがって、図２９に示した例では、誤差率αの算出は、下記式（５３）によって行われる。式（５３）からわかるように、吸気管内圧力の実測値Ｐｍｒの積分は、マップ領域Ｂに侵入・離脱した時刻ｔ₀から時刻ｔ₁に亘って行われるのではなく、中心領域Ｂ’に侵入・離脱した時刻ｔ₃から時刻ｔ₄に亘って行われる。これにより、この誤差率αに基づいて修正された値ａおよび値ｂのマップ値は、或る特定のマップ領域の中心の機関回転数および吸気弁開弁時期に対応するような値とされ、これにより値ａ、ｂのマップ値をそのマップ領域内の如何なる機関回転数ＮＥおよび吸気弁開弁時期ＶＶＴに対しても比較的適した値とすることができる。

本発明の吸入空気量推定装置を搭載した内燃機関全体を示す図である。図１に示したエアフロメータの概略斜視図である。図２に示したエアフロメータの熱線計量部の拡大斜視図である。エアフロメータの出力電圧とエアフロ通過空気流量との関係を示す図である。エアモデルのブロック図である。スロットル弁開度と流量係数との関係を示す図である。スロットル弁開度と開口断面積との関係を示す図である。関数Φ（Ｐｍ／Ｐａ）を示す図である。スロットルモデルの基本概念を示す図である。吸気管モデルの基本概念を示す図である。吸気弁モデルの基本概念を示す図である。筒内充填ガス量および筒内吸気ガス量の定義に関する図である。スロットルモデルの修正操作を示す図である。スロットル通過空気流量と完全放熱量との関係を示す図である。遅れ放熱量の和とエアフロメータの予想出力値との関係を示す図である。スロットル弁開度と空気流量に関する各種パラメータの値とのタイムチャートである。スロットル弁開度に基づいて定められた流量係数のマップを示すである。スロットル弁開度と空気流量に関する各種パラメータの値とのタイムチャートである。図１８の一部を拡大した図である。初期値応答とゼロ応答とを説明するための図である。スロットル弁開度と空気流量に関する各種パラメータの値とのタイムチャートであって図１９と同様な図である。第一実施形態の流量係数の修正制御の制御ルーチンを示すフローチャートである。スロットル弁開度と空気流量に関する各種パラメータの値とのタイムチャートであって図１９および図２１と同様な図である。第二実施形態の流量係数の修正制御の制御ルーチンを示すフローチャートである。機関回転数および吸気弁開弁時期に基づいて定められた値ａ、ｂのマップを示す図である。吸気管内圧力、機関回転数、吸気弁開弁時期のタイムチャートである。吸気管内圧力のタイムチャートである。吸気管内圧力と筒内吸入空気流量との関係を示した図である。吸気管内圧力、機関回転数、吸気弁開弁時期のタイムチャートであり、図２６と同様な図である。

符号の説明

１機関本体
５燃焼室
６吸気弁
７吸気ポート
８排気弁
１１燃料噴射弁
１４吸気枝管
１５サージタンク
１６吸気管
１９スロットル弁
２２吸気管部分

Claims

機関運転領域毎にマップ値が定められていると共に、実際の機関運転状態が特定の機関運転領域にある場合にはその機関運転領域に対応するマップ値を用いて内燃機関を制御する、内燃機関の制御装置において、
実際の機関運転状態が特定の機関運転領域にある場合における運転パラメータの実測値を検出する実測値検出手段と、
実際の機関運転状態が上記特定の機関運転領域にある場合に、該特定の機関運転領域に対応するマップ値を用いて運転パラメータの推定値を算出する推定値算出手段と、
実際の機関運転状態が上記特定の機関運転領域にある場合において上記実測値検出手段によって検出された運転パラメータの実測値と上記推定値算出手段によって算出された推定値とに基づいて上記特定の機関運転領域に対応するマップ値を修正するマップ値修正手段とを具備し、
上記マップ値修正手段は、実際の機関運転状態が上記特定の機関運転領域に侵入してから所定の最長学習時間よりも長く該特定の機関運転領域にある場合には該最長学習時間が経過するまでの上記運転パラメータの実測値と推定値に基づいてマップ値を修正する、内燃機関の制御装置。
上記最長学習時間は機関運転状態に応じて定められる、請求項１に記載の内燃機関の制御装置。
機関運転領域毎にマップ値が定められていると共に、実際の機関運転状態が特定の機関運転領域にある場合にはその機関運転領域に対応するマップ値を用いて内燃機関を制御する、内燃機関の制御装置において、
実際の機関運転状態が特定の機関運転領域にある場合における運転パラメータの実測値を検出する実測値検出手段と、
実際の機関運転状態が上記特定の機関運転領域にある場合に、該特定の機関運転領域に対応するマップ値を用いて上記運転パラメータの推定値を算出する推定値算出手段と、
実際の機関運転状態が上記特定の機関運転領域にある場合において上記実測値検出手段によって検出された運転パラメータの実測値と上記推定値算出手段によって算出された推定値とに基づいて上記特定の機関運転領域に対応するマップ値を修正するマップ値修正手段とを具備し、
上記マップ値修正手段は、実際の機関運転状態が上記特定の機関運転領域に侵入してから所定の最長学習時間よりも長く該特定の機関運転領域にある場合には、実際の機関運転状態が上記特定の機関運転領域にある場合における上記運転パラメータの実測値と推定値とのうち実際の機関運転状態が該特定の機関運転領域の中心領域にある場合における上記運転パラメータの実測値と推定値とを選択し、選択された上記運転パラメータの実測値と推定値とに基づいて当該機関運転領域に対応するマップ値を修正する、内燃機関の制御装置。
上記マップ値修正手段は、実際の機関運転状態が上記特定の機関運転領域に侵入してから上記最長学習時間よりも短い最短学習時間を経過することなく該特定の機関運転領域から離脱した場合には、マップ値の修正を行わない、請求項１〜３のいずれか１項に記載の内燃機関の制御装置。
上記マップ値が流量係数として用いられるスロットルモデル計算式であってスロットル弁開度と吸気管内圧力とを用いてスロットル通過空気流量の推定値を算出するスロットルモデル計算式によって算出されたスロットル通過空気流量の推定値を用いて内燃機関が制御され、
上記実測値検出手段は所定の遅れをもってスロットル通過空気流量に対応する出力値を出力するエアフロメータであり、
上記推定値算出手段はスロットル弁開度と吸気管内圧力とを用いて上記スロットルモデル計算式によって算出されるスロットル通過空気流量の推定値を用いてエアフロメータによって出力されるであろう出力値をエアフロメータの出力値の予想値として算出する予想値算出手段であり、
上記マップ値修正手段は、エアフロメータの出力値とその予想値とに基づいて上記スロットルモデル計算式に用いられるマップ値を修正する、請求項１〜４のいずれか１項に記載の内燃機関の制御装置。
上記マップ値および吸気管内圧力を用いて筒内吸入空気流量の推定値を算出する吸気弁モデル計算式によって算出された筒内吸入空気流量の推定値を用いて内燃機関が制御され、
上記実測値検出手段は、実際の吸気管内圧力を検出する吸気管内圧力検出手段であり、
上記推定値算出手段は、筒内吸入空気流量を用いて吸気管モデル計算式によって吸気管内圧力の推定値を算出する吸気管内圧力推定値算出手段であり、
上記マップ値修正手段は、上記吸気管内圧力推定手段によって算出された吸気管内圧力の推定値と上記吸気管内圧力検出手段によって検出された実際の吸気管内圧力とに基づいて上記吸気弁モデル計算式を修正する、請求項１〜４のいずれか１項に記載の内燃機関の制御装置。