"吸気通路に設けられた吸気絞り弁によって当該吸気絞り弁を通過する空気量が通過空気量として変更される内燃機関において、当該吸気通路内の空気の状態を表す吸気パラメータを算出する内燃機関の吸気パラメータ算出装置であって、前記吸気パラメータの基本値として基本吸気パラメータを算出する基本吸気パラメータ算出手段と、所定のモデリング手法によって導出された、前記吸気絞り弁の上流側おける前記吸気通路内の圧力である上流側圧力と、前記吸気絞り弁の下流側における前記吸気通路内の圧力である下流側圧力と、前記吸気絞り弁の開度によって定まる開度関数値と、前記通過空気量との関係を定義したモデル式を用いて、当該開度関数値の第1の算出値として第1開度関数値を算出する第1開度関数値算出手段と、前記吸気絞り弁の開度と前記開度関数値との相関関係を表す相関関係モデルを用いて、前記開度関数値の第2の算出値として第2開度関数値を算出する第2開度関数値算出手段と、前記算出された第1開度関数値および前記算出された第2開度関数値の一方と他方との比である関数値比を用いて、補正値を算出する補正値算出手段と、当該算出された補正値で前記基本吸気パラメータを補正することにより、前記吸気パラメータを算出する吸気パラメータ算出手段と、を備えることを特徴とする内燃機関の吸気パラメータ算出装置。前記基本吸気パラメータ算出手段は、前記基本吸気パラメータとして、前記通過空気量の基本値である基本通過空気量を算出する基本通過空気量算出手段を有し、前記補正値算出手段は、前記補正値として、当該基本通過空気量を補正するための第1補正値を算出する第1補正値算出手段を有し、前記吸気パラメータ算出手段は、前記算出された基本通過空気量を前記算出された第1補正値で補正することにより、前記吸気パラメータとして、前記通過空気量を算出する通過空気量算出手段を有することを特徴とする請求項1に記載の内燃機関の吸気パラメータ算出装置。前記第1補正値算出手段は、前記関数値比を誤差と見なしたときの、当該誤差と前記吸気絞り弁の開度との関係を定義した誤差モデル式のモデルパラメータを、オンボードで同定するオンボード同定手段を有し、当該オンボードで同定されたモデルパラメータおよび当該誤差モデル式を用いて、前記第1補正値を算出することを特徴とする請求項2に記載の内燃機関の吸気パラメータ算出装置。前記オンボード同定手段は、前記吸気絞り弁の開度を複数の領域に区分したときの当該複数の領域に対応して複数の重みをそれぞれ設定するとともに、当該複数の重みによる重み付け処理を適用した同定演算アルゴリズムにより、前記モデルパラメータの同定値を算出し、当該複数の重みは互いに同じ値に設定されていることを特徴とする請求項3に記載の内燃機関の吸気パラメータ算出装置。前記オンボード同定手段は、前記吸気絞り弁の開度を複数の領域に区分したときの当該複数の領域に対応して複数の重みをそれぞれ設定するとともに、当該複数の重みによる重み付け処理を適用した同定演算アルゴリズムにより、前記モデルパラメータの同定値を算出し、当該複数の重みは、前記吸気絞り弁の開度がより小さい領域の重みほど、より大きくなるように設定されていることを特徴とする請求項3に記載の内燃機関の吸気パラメータ算出装置。前記内燃機関は、車両の動力源として用いられ、前記第1補正値算出手段は、前記内燃機関の始動時点から所定時間が経過していないこと、および前記内燃機関の始動後における前記車両の総走行距離が所定値未満であることの一方の条件が成立しているときに、前記第1補正値の算出を実行することを特徴とする請求項2ないし5のいずれかに記載の内燃機関の吸気パラメータ算出装置。前記上流側圧力の推定値として推定上流側圧力を算出する推定上流側圧力算出手段をさらに備え、前記第1開度関数値算出手段は、当該算出された推定上流側圧力および前記モデル式を用いて、前記第1開度関数値を算出することを特徴とする請求項2ないし6のいずれかに記載の内燃機関の吸気パラメータ算出装置。前記下流側圧力を検出する下流側圧力検出手段をさらに備え、前記推定上流側圧力算出手段は、前記推定上流側圧力を、前記内燃機関を始動する際に前記下流側圧力検出手段によって検出された下流側圧力に設定することを特徴とする請求項7に記載の内燃機関の吸気パラメータ算出装置。前記吸気通路の前記吸気絞り弁よりも上流側に設けられ、前記吸気通路内を流れる空気量である吸入空気量を表す検出信号を出力する吸入空気量検出手段と、当該吸入空気量検出手段の検出結果に基づき、前記吸入空気量を算出する吸入空気量算出手段と、をさらに備え、前記推定上流側圧力算出手段は、前記第1補正値および前記モデル式を用いて、前記通過空気量の推定値である推定通過空気量を算出するとともに、当該推定通過空気量と前記算出された吸入空気量との比較結果に基づいて、前記推定上流側圧力を更新することを特徴とする請求項7または8に記載の内燃機関の吸気パラメータ算出装置。前記第1補正値算出手段は、前記関数値比および前記第1補正値の一方の値の相加平均値を、前記吸気絞り弁の開度の所定間隔ごとの値に対応して算出する相加平均値算出手段と、当該算出された相加平均値を、記憶用値として前記吸気絞り弁の開度の所定間隔ごとの値に対応して記憶する記憶手段と、を有し、前記吸気絞り弁の開度に応じて前記記憶手段から読み出した前記記憶用値を、前記一方の値として用いることを特徴とする請求項2に記載の内燃機関の吸気パラメータ算出装置。前記吸気通路の前記吸気絞り弁よりも上流側に設けられ、前記吸気通路内を流れる空気量である吸入空気量を表す検出信号を出力する吸入空気量検出手段と、前記上流側圧力を表す検出信号を出力する上流側圧力検出手段と、前記下流側圧力を表す検出信号を出力する下流側圧力検出手段と、前記上流側圧力検出手段および前記下流側圧力検出手段の検出結果に基づき、前記下流側圧力と前記上流側圧力との比である圧力比を算出する圧力比算出手段と、前記吸気絞り弁の開度として弁開度を検出する弁開度検出手段と、をさらに備え、前記基本吸気パラメータ算出手段は、前記基本吸気パラメータとして、前記吸入空気量検出手段の検出結果に基づき、前記吸入空気量の基本値である基本吸入空気量を算出する基本吸入空気量算出手段をさらに有し、前記補正値算出手段は、前記補正値として、前記検出された弁開度が所定開度以上でかつ前記算出された圧力比が臨界圧力比未満のときに、前記基本吸入空気量を補正するための第2補正値を算出する第2補正値算出手段をさらに有し、前記吸気パラメータ算出手段は、前記算出された第2補正値で前記基本吸入空気量を補正することにより、前記吸気パラメータとして、前記吸入空気量を算出する吸入空気量算出手段をさらに有し、前記第1開度関数値算出手段は、前記通過空気量として前記吸入空気量を用いて、前記第1開度関数値を算出することを特徴とする請求項2ないし8のいずれかに記載の内燃機関の吸気パラメータ算出装置。前記基本吸気パラメータ算出手段は、前記下流側圧力検出手段の検出結果に基づき、前記基本吸気パラメータとして、前記下流側圧力の基本値である基本下流側圧力を算出する基本下流側圧力算出手段をさらに有し、前記補正値算出手段は、前記圧力比が臨界圧力比以上のときに、前記補正値として、前記基本下流側圧力を補正するための第3補正値を算出する第3補正値算出手段をさらに有し、前記吸気パラメータ算出手段は、前記算出された第3補正値で前記基本下流側圧力を補正することにより、前記吸気パラメータとして、前記下流側圧力を算出する下流側圧力算出手段をさらに有することを特徴とする請求項11に記載の内燃機関の吸気パラメータ算出装置。前記吸気通路の前記吸気絞り弁よりも上流側に設けられ、前記吸気通路内を流れる空気量である吸入空気量を表す検出信号を出力する吸入空気量検出手段と、前記上流側圧力を表す検出信号を出力する上流側圧力検出手段と、前記下流側圧力を表す検出信号を出力する下流側圧力検出手段と、前記上流側圧力検出手段および前記下流側圧力検出手段の検出結果に基づき、前記下流側圧力と前記上流側圧力との比である圧力比を算出する圧力比算出手段と、前記吸気絞り弁の開度として弁開度を検出する弁開度検出手段と、をさらに備え、前記基本吸気パラメータ算出手段は、前記算出された圧力比が臨界圧力比以上のときに、前記基本吸気パラメータとして、前記下流側圧力検出手段の検出結果に基づき、前記下流側圧力の基本値である基本下流側圧力を算出する基本下流側圧力算出手段をさらに有し、前記補正値算出手段は、前記圧力比が臨界圧力比以上のときに、前記補正値として、前記基本下流側圧力を補正するための第3補正値を算出する第3補正値算出手段をさらに有し、前記吸気パラメータ算出手段は、前記圧力比が臨界圧力比以上のときに、前記算出された第3補正値で前記基本下流側圧力を補正することにより、前記吸気パラメータとして、前記下流側圧力を算出する下流側圧力算出手段をさらに有することを特徴とする請求項2に記載の内燃機関の吸気パラメータ算出装置。前記吸気通路の前記吸気絞り弁よりも上流側に設けられ、前記吸気通路内を流れる空気量である吸入空気量を表す検出信号を出力する吸入空気量検出手段と、前記上流側圧力を表す検出信号を出力する上流側圧力検出手段と、前記下流側圧力を表す検出信号を出力する下流側圧力検出手段と、前記上流側圧力検出手段および前記下流側圧力検出手段の検出結果に基づき、前記下流側圧力と前記上流側圧力との比である圧力比を算出する圧力比算出手段と、前記吸気絞り弁の開度として弁開度を検出する弁開度検出手段と、をさらに備え、前記基本吸気パラメータ算出手段は、前記基本吸気パラメータとして、前記吸入空気量検出手段の検出結果に基づき、前記吸入空気量の基本値である基本吸入空気量を算出する基本吸入空気量算出手段を有し、前記補正値算出手段は、前記検出された弁開度が所定開度以上でかつ前記算出された圧力比が臨界圧力比未満のときに、前記補正値として、前記基本吸入空気量を補正するための第2補正値を算出する第2補正値算出手段をさらに有し、前記吸気パラメータ算出手段は、前記算出された第2補正値で前記基本吸入空気量を補正することにより、前記吸気パラメータとして、前記吸入空気量を算出する吸入空気量算出手段を有し、前記第1開度関数値算出手段は、前記通過空気量として前記吸入空気量を用いて、前記第1開度関数値を算出することを特徴とする請求項1に記載の内燃機関の吸気パラメータ算出装置。前記基本吸気パラメータ算出手段は、前記下流側圧力検出手段の検出結果に基づき、前記基本吸気パラメータとして、前記下流側圧力の基本値である基本下流側圧力を算出する基本下流側圧力算出手段をさらに有し、前記補正値算出手段は、前記圧力比が臨界圧力比以上のときに、前記補正値として、前記基本下流側圧力を補正するための第3補正値を算出する第3補正値算出手段をさらに有し、前記吸気パラメータ算出手段は、前記圧力比が臨界圧力比以上のときに、前記算出された第3補正値で前記基本下流側圧力を補正することにより、前記吸気パラメータとして、前記下流側圧力を算出する下流側圧力算出手段をさらに有することを特徴とする請求項14に記載の内燃機関の吸気パラメータ算出装置。前記上流側圧力を表す検出信号を出力する上流側圧力検出手段と、当該上流側圧力検出手段の検出結果に基づき、前記上流側圧力を算出する上流側圧力算出手段と、前記下流側圧力を表す検出信号を出力する下流側圧力検出手段と、前記上流側圧力検出手段および前記下流側圧力検出手段の検出結果に基づき、前記下流側圧力と前記上流側圧力との比である圧力比を算出する圧力比算出手段と、をさらに備え、前記基本吸気パラメータ算出手段は、前記下流側圧力検出手段の検出結果に基づき、前記基本吸気パラメータとして、前記下流側圧力の基本値である基本下流側圧力を算出する基本下流側圧力算出手段をさらに有し、前記補正値算出手段は、前記圧力比が臨界圧力比以上のときに、前記補正値として、前記基本下流側圧力を補正するための第3補正値を算出する第3補正値算出手段をさらに有し、前記吸気パラメータ算出手段は、前記圧力比が臨界圧力比以上のときに、前記算出された第3補正値で前記基本下流側圧力を補正することにより、前記吸気パラメータとして、前記下流側圧力を算出する下流側圧力算出手段をさらに有することを特徴とする請求項1に記載の内燃機関の吸気パラメータ算出装置。吸気通路に設けられた吸気絞り弁によって当該吸気絞り弁を通過する空気量が通過空気量として変更される内燃機関において、当該吸気通路内の空気の状態を表す吸気パラメータを算出する内燃機関の吸気パラメータ算出方法であって、前記吸気パラメータの基本値として基本吸気パラメータを算出し、所定のモデリング手法によって導出された、前記吸気絞り弁の上流側おける前記吸気通路内の圧力である上流側圧力と、前記吸気絞り弁の下流側における前記吸気通路内の圧力である下流側圧力と、前記吸気絞り弁の開度によって定まる開度関数値と、前記通過空気量との関係を定義したモデル式を用いて、当該開度関数値の第1の算出値として第1開度関数値を算出し、前記吸気絞り弁の開度と前記開度関数値との相関関係を表す相関関係モデルを用いて、前記開度関数値の第2の算出値として第2開度関数値を算出し、前記算出された第1開度関数値および前記算出された第2開度関数値の一方と他方との比である関数値比を用いて、補正値を算出し、当該算出された補正値で前記基本吸気パラメータを補正することにより、前記吸気パラメータを算出することを特徴とする内燃機関の吸気パラメータ算出方法。特許文献1の内燃機関のように、スロットル弁などの、吸気通路の開口面積を変更する吸気絞り弁が吸気通路に設けられている場合、吸気圧は、吸気通路の吸気絞り弁の下流側の圧力である関係上、吸気絞り弁の開度変化の影響を受けやすいという特性を有している。これに対して、特許文献1の吸気パラメータ算出装置の場合、要求吸気圧Pを、吸気絞り弁の開度を考慮することなく、気体の状態方程式から導出した吸気系モデルを用いて算出しているものに過ぎないので、算出精度が低いという問題がある。特に、過渡運転状態のように、吸気絞り弁の開度が変化しやすい条件下では、算出精度の低下度合が大きくなってしまう。本発明は、上記課題を解決するためになされたもので、吸気絞り弁が設けられている場合において、吸気パラメータを精度よく算出することができる内燃機関の吸気パラメータ算出装置および吸気パラメータ算出方法を提供することを目的とする。上記目的を達成するために、請求項1に係る発明は、吸気通路6に設けられた吸気絞り弁スロットル弁7aによって吸気絞り弁を通過する空気量が通過空気量として変更される内燃機関3において、吸気通路6内の空気の状態を表す吸気パラメータ通過空気量GAIRTH、吸入空気量GAIR、吸気圧PBAを算出する内燃機関3の吸気パラメータ算出装置1であって、吸気パラメータの基本値として基本吸気パラメータ基本通過空気量GAIRTHN,基本吸入空気量GAIRN,基本吸気圧PBANを算出する基本吸気パラメータ算出手段ECU2、ステップ6,90,110,144,175と、所定のモデリング手法によって導出された、吸気絞り弁の上流側おける吸気通路6内の圧力である上流側圧力大気圧PAと、吸気絞り弁の下流側における吸気通路6内の圧力である下流側圧力吸気圧PBAと、吸気絞り弁の開度によって定まる開度関数値KTHと、通過空気量GAIRTHとの関係を定義したモデル式[式7,13,23]を用いて、開度関数値の第1の算出値として第1開度関数値モデル式値KTHCALを算出する第1開度関数値算出手段ECU2、ステップ11,103,123,154,183と、吸気絞り弁の開度スロットル弁開度THと開度関数値KTHとの相関関係を表す相関関係モデル図2を用いて、開度関数値の第2の算出値として第2開度関数値マップ値KTHを算出する第2開度関数値算出手段ECU2、ステップ11,100,120,151,180と、算出された第1開度関数値および算出された第2開度関数値の一方と他方との比である関数値比関数値誤差KTHERRを用いて、補正値補正係数KTHCOR、補正係数KAFMERR、補正項PBAERRCORを算出する補正値算出手段ECU2、ステップ3,100~107,120~131,142,146,147,172と、算出された補正値で基本吸気パラメータを補正することにより、吸気パラメータ通過空気量GAIRTH、吸入空気量GAIR、吸気圧PBAを算出する吸気パラメータ算出手段ECU2、ステップ6,91,111,144,175と、を備えることを特徴とする。この内燃機関の吸気パラメータ算出装置によれば、第1開度関数値が、吸気絞り弁の上流側おける吸気通路内の圧力である上流側圧力と、吸気絞り弁の下流側における吸気通路内の圧力である下流側圧力と、吸気絞り弁の開度によって定まる開度関数値と、通過空気量との関係を定義したモデル式を用いて算出され、第2開度関数値が、吸気絞り弁の開度と開度関数値との相関関係を表す相関関係モデルを用いて算出され、算出された第1開度関数値および算出された第2開度関数値の一方と他方との比である関数値比を用いて、補正値が算出され、このように算出された補正値で基本吸気パラメータを補正することにより、吸気パラメータが算出される。この場合、第1および第2開度関数値はいずれも、吸気絞り弁の開度によって定まる関数値であるとともに、そのような2つの値の比を用いて算出した補正値で基本吸気パラメータを補正することにより、吸気パラメータが算出されるので、吸気絞り弁の開度状態を反映させながら、吸気パラメータを算出することができる。これに加えて、第1開度関数値および第2開度関数値の一方と他方との比である関数値比は、2つの開度関数値の間の差異、すなわちモデル式と相関関係モデルとの間の誤差を表すものとなるので、補正値を、そのような誤差を補正できるものとして算出することができる。以上により、過渡運転状態のように、吸気絞り弁の開度が変化しやすい条件下でも、吸気パラメータを従来よりも精度よく算出することができ、商品性を向上させることができる。請求項2に係る発明は、請求項1に記載の内燃機関3の吸気パラメータ算出装置1において、基本吸気パラメータ算出手段は、基本吸気パラメータとして、通過空気量GAIRTHの基本値である基本通過空気量GAIRTHNを算出する基本通過空気量算出手段ECU2、ステップ6,144,175を有し、補正値算出手段は、補正値として、基本通過空気量を補正するための第1補正値補正係数KTHCORを算出する第1補正値算出手段ECU2、ステップ3,142,146,147を有し、吸気パラメータ算出手段は、算出された基本通過空気量を算出された第1補正値で補正することにより、吸気パラメータとして、通過空気量GAIRTHを算出する通過空気量算出手段ECU2、ステップ6,144,175を有することを特徴とする。この内燃機関の吸気パラメータ算出装置によれば、基本吸気パラメータとして、通過空気量の基本値である基本通過空気量が算出され、補正値として、基本吸入空気量を補正するための第1補正値が算出され、算出された基本通過空気量を算出された第1補正値で補正することにより、吸気パラメータとして、通過空気量が算出されるので、吸気絞り弁の開度状態を反映させながら、通過空気量を算出することができる。これに加え、第1補正値が、モデル式と相関関係モデルとの間の誤差を補正できるものとして算出されるので、通過空気量を、そのような誤差が補正された値として算出することができる。以上により、通過空気量を精度よく算出することができる。請求項3に係る発明は、請求項2に記載の内燃機関3の吸気パラメータ算出装置1において、第1補正値算出手段は、関数値比を誤差と見なしたときの、誤差と吸気絞り弁の開度との関係を定義した誤差モデル式8のモデルパラメータAを、オンボードで同定するオンボード同定手段ECU2、ステップ17,40~53,157を有し、オンボードで同定されたモデルパラメータAおよび誤差モデル式を用いて、第1補正値を算出することを特徴とする。この内燃機関の吸気パラメータ算出装置によれば、関数値比を誤差と見なしたときの、誤差と吸気絞り弁の開度との関係を定義した誤差モデル式のモデルパラメータが、オンボードで同定され、オンボードで同定されたモデルパラメータおよび誤差モデル式を用いて、第1補正値が算出されるので、吸気絞り弁における経年変化および個体間のばらつきなどに起因して、誤差モデル式が誤差と吸気絞り弁の開度との間の実際の関係から乖離した場合、すなわちモデル化誤差を生じた場合でも、オンボード同定されたモデルパラメータを用いることによって、そのようなモデル化誤差を迅速に補償でき、誤差モデル式を誤差と吸気絞り弁の開度との実際の関係に迅速に適合させることができる。それにより、補正値による補正精度を向上させることができ、通過空気量の算出精度を向上させることができる。請求項4に係る発明は、請求項3に記載の内燃機関3の吸気パラメータ算出装置1において、オンボード同定手段は、吸気絞り弁の開度を複数の領域に区分したときの複数の領域に対応して複数の重みをそれぞれ設定するとともに、複数の重みによる重み付け処理を適用した同定演算アルゴリズム[式16~18]により、モデルパラメータAの同定値を算出し、複数の重みは互いに同じ値に設定されていることを特徴とする。この内燃機関の吸気パラメータ算出装置によれば、吸気絞り弁の開度を複数の領域に区分したときの複数の領域の各々に対応して複数の重みがそれぞれ設定されるとともに、複数の重みによる重み付け処理を適用した同定演算アルゴリズムにより、モデルパラメータの同定値が算出される。そして、複数の重みは、互いに同じ値に設定されているので、オンボード同定演算中、吸気絞り弁の開度が複数の領域のいずれかに一時的に偏り、その領域に対応する演算結果が他の領域の演算結果よりも一時的に非常に多い状態となったときでも、モデルパラメータの同定値を、そのように偏った演算結果が反映されるのを抑制しながら適切に算出することができ、良好な算出精度を確保することができる。それにより、通過空気量の算出精度をさらに向上させることができる。請求項5に係る発明は、請求項3に記載の内燃機関3の吸気パラメータ算出装置において、オンボード同定手段は、吸気絞り弁の開度を複数の領域に区分したときの複数の領域に対応して複数の重みをそれぞれ設定するとともに、複数の重み重み係数K1~KG4による重み付け処理を適用した同定演算アルゴリズム[式18,27,28]により、モデルパラメータAの同定値を算出し、複数の重みは、吸気絞り弁の開度がより小さい領域の重みほど、より大きくなるように設定されていることを特徴とする。この内燃機関の吸気パラメータ算出装置によれば、吸気絞り弁の開度を複数の領域に区分したときの複数の領域の各々に対応して複数の重みがそれぞれ設定されるとともに、複数の重みによる重み付け処理を適用した同定演算アルゴリズムにより、モデルパラメータの同定値が算出される。そして、複数の重みは、吸気絞り弁の開度がより小さい領域の重みほど、より大きくなるように設定されている。この場合、後述するように、通過空気量の算出誤差は、吸気絞り弁の開口面積の、標準品基準品の値に対するずれに起因して発生するものであり、そのずれの影響度合は、吸気絞り弁の開度が小さいほど、より大きくなる。したがって、複数の重みを、吸気絞り弁の開度がより小さい領域の重みほど、より大きくなるように設定することによって、吸気絞り弁の開口面積の、標準品の値に対するずれの影響度合を反映させながら、モデルパラメータの同定値を算出することができ、それにより、第1補正値の算出精度をさらに向上させることができる。その結果、通過空気量の算出精度をさらに向上させることができる。請求項6に係る発明は、請求項2ないし5のいずれかに記載の内燃機関3の吸気パラメータ算出装置1において、内燃機関3は、車両の動力源として用いられ、第1補正値算出手段は、内燃機関3の始動時点から所定時間が経過していないこと、および内燃機関3の始動後における車両の総走行距離DISTが所定値DLEARN未満であることの一方の条件が成立しているときステップ11,20の判別結果がYESのときに、第1補正値の算出を実行することを特徴とする。この内燃機関の吸気パラメータ算出装置によれば、内燃機関の始動時点から所定時間が経過していないこと、および内燃機関の始動後における車両の総走行距離が所定値未満であることの一方の条件が成立しているときに、第1補正値の算出が実行される。この場合、後述するように、通過空気量の算出誤差は、吸気絞り弁の開口面積の、標準品の値に対するずれに起因して発生するものであり、吸気絞り弁の開度が小さい低開度域での発生頻度が高いとともに、内燃機関の始動時点から所定時間が経過していないときや、内燃機関の始動後における車両の総走行距離が所定値未満であるときには、吸気絞り弁の開度が低開度域で推移しやすい状態となる。以上により、第1補正値の算出を、通過空気量の算出誤差の発生しやすい領域でのみ実行することによって、その算出精度を高めることができる。請求項7に係る発明は、請求項2ないし6のいずれかに記載の内燃機関3の吸気パラメータ算出装置において、上流側圧力の推定値として推定上流側圧力推定大気圧HPAを算出する推定上流側圧力算出手段ECU2、ステップ61,64をさらに備え、第1開度関数値算出手段は、算出された推定上流側圧力およびモデル式を用いて、第1開度関数値を算出することを特徴とする。この内燃機関の吸気パラメータ算出装置によれば、算出された推定上流側圧力およびモデル式を用いて、第1開度関数値が算出されるので、上流側圧力の変化に起因して、第1補正値の算出精度が低下するのを回避することができ、良好な算出精度を確保することができる。請求項8に係る発明は、請求項7に記載の内燃機関3の吸気パラメータ算出装置において、下流側圧力吸気圧PBAを検出する下流側圧力検出手段ECU2、吸気圧センサ22をさらに備え、推定上流側圧力算出手段は、推定上流側圧力を、内燃機関3を始動する際に下流側圧力検出手段によって検出された下流側圧力に設定するステップ61ことを特徴とする。この内燃機関の吸気パラメータ算出装置によれば、推定上流側圧力が、内燃機関を始動する際に下流側圧力検出手段によって検出された下流側圧力に設定されるので、第1開度関数値がそのような下流側圧力を用いて算出されることになる。この場合、内燃機関の停止中は、吸気絞り弁の上流側と下流側の圧力がほぼ同一になるので、内燃機関を始動する際に下流側圧力検出手段によって検出された下流側圧力は、上流側圧力を適切に表すものとなる。それにより、上流側圧力を適切に表す下流側圧力を用いて第1開度関数値を精度よく算出することができる。これに加えて、内燃機関の始動後、上流側圧力の変動が生じた場合でも、その影響を回避することができ、それにより、第1開度関数値の算出精度をさらに向上させることができる。請求項9に係る発明は、請求項7または8に記載の内燃機関の吸気パラメータ算出装置において、吸気通路6の吸気絞り弁よりも上流側に設けられ、吸気通路6内を流れる空気量である吸入空気量GAIRを表す検出信号を出力する吸入空気量検出手段エアフローセンサ20と、吸入空気量検出手段の検出結果に基づき、吸入空気量GAIRを算出する吸入空気量算出手段ECU2、ステップ90,91と、をさらに備え、推定上流側圧力算出手段は、第1補正値補正係数KTHCORおよびモデル式を用いて、通過空気量の推定値である推定通過空気量HGAIRTHを算出するとともに、推定通過空気量HGAIRTHと算出された吸入空気量GAIRとの比較結果空気量偏差DGAIRに基づいて、推定上流側圧力推定大気圧HPAを更新するステップ73~83ことを特徴とする。この内燃機関の吸気パラメータ算出装置によれば、吸入空気量検出手段の検出結果に基づいて算出された吸入空気量と、推定通過空気量との比較結果に基づいて、推定上流側圧力が更新される。この推定通過空気量は、第1補正値およびモデル式を用いて算出されるので、その算出結果において、前述したモデル化誤差を補償しながら、高い算出精度を確保することができる。したがって、そのような推定通過空気量と吸入空気量との比較結果に基づいて、推定上流側圧力を更新することによって、推定上流側圧力の更新精度すなわち推定上流側圧力の算出精度を向上させることができる。請求項10に係る発明は、請求項2に記載の内燃機関3の吸気パラメータ算出装置において、第1補正値算出手段は、関数値比および第1補正値の一方の値の相加平均値を、吸気絞り弁の開度の所定間隔1 ゚ごとの値に対応して算出する相加平均値算出手段ECU2、ステップ186と、算出された相加平均値を、記憶用値として吸気絞り弁の開度の所定間隔ごとの値に対応して記憶する記憶手段ECU2、ステップ186と、を有し、吸気絞り弁の開度に応じて記憶手段から読み出した記憶用値を、一方の値として用いることを特徴とする。この内燃機関の吸気パラメータ算出装置によれば、関数値比および第1補正値の一方の値の相加平均値が、吸気絞り弁の開度の所定間隔ごとの値に対応して算出され、算出された相加平均値が記憶用値として吸気絞り弁の開度の所定間隔ごとの値に対応して記憶されるとともに、吸気絞り弁の開度に応じて記憶手段から読み出した記憶用値が、一方の値として用いられる。この場合、記憶用値は、関数値比または第1補正値の相加平均値であるので、相加平均演算が進行するほど、その算出精度を高めることができる。これに加えて、記憶用値が吸気絞り弁の開度の所定間隔ごとの値に対応して記憶されるので、その所定間隔が小さいほど、記憶用値の算出精度を高めることができる。したがって、そのような記憶用値を関数値比または第1補正値として用いることにより、通過空気量が算出されるので、演算の進行に伴い、通過空気量の算出精度を向上させることができる。これに加えて、所定間隔を小さく設定することによって、その算出精度をさらに向上させることができる。請求項11に係る発明は、請求項2ないし8のいずれかに記載の内燃機関3の吸気パラメータ算出装置において、吸気通路6の吸気絞り弁よりも上流側に設けられ、吸気通路6内を流れる空気量である吸入空気量GAIRを表す検出信号を出力する吸入空気量検出手段エアフローセンサ20と、上流側圧力大気圧PAを表す検出信号を出力する上流側圧力検出手段大気圧センサ24と、下流側圧力吸気圧PBAを表す検出信号を出力する下流側圧力検出手段吸気圧センサ22と、上流側圧力検出手段および下流側圧力検出手段の検出結果に基づき、下流側圧力と上流側圧力との比である圧力比R_Pを算出する圧力比算出手段ECU2と、吸気絞り弁の開度として弁開度スロットル弁開度THを検出する弁開度検出手段ECU2、スロットル弁開度センサ21と、をさらに備え、基本吸気パラメータ算出手段は、基本吸気パラメータとして、吸入空気量検出手段の検出結果に基づき、吸入空気量GAIRの基本値である基本吸入空気量GAIRNを算出する基本吸入空気量算出手段ECU2、ステップ90をさらに有し、補正値算出手段は、補正値として、検出された弁開度が所定開度THB以上でかつ算出された圧力比R_Pが臨界圧力比R_CRIT未満のときに、基本吸入空気量GAIRNを補正するための第2補正値補正係数KAFMERRを算出する第2補正値算出手段ECU2、ステップ100~107をさらに有し、吸気パラメータ算出手段は、算出された第2補正値補正係数KAFMERRで基本吸入空気量GAIRNを補正することにより、吸気パラメータとして、吸入空気量GAIRを算出する吸入空気量算出手段ECU2、ステップ91をさらに有し、第1開度関数値算出手段は、通過空気量として吸入空気量GAIRを用いて、第1開度関数値を算出することを特徴とする。この内燃機関の吸気パラメータ算出装置によれば、基本吸気パラメータとして、吸入空気量検出手段の検出結果に基づき、吸入空気量の基本値である基本吸入空気量が算出され、補正値として、検出された弁開度が所定開度以上でかつ算出された圧力比が臨界圧力比未満のときに、基本吸入空気量を補正するための第2補正値が算出され、算出された第2補正値で基本吸入空気量を補正することにより、吸気パラメータとして、吸入空気量が算出される。この吸気パラメータ算出装置のように、上流側圧力と、下流側圧力と、開度関数値と、通過空気量すなわち吸入空気量との関係を定義したモデル式を用いて、第1開度関数値を算出し、さらに、関数値比を算出した場合、関数値比が表す、モデル式と相関関係モデルとの間の誤差は、弁開度が所定開度以上でかつ算出された圧力比が臨界圧力比未満の領域では、吸入空気量検出手段の検出結果に含まれる誤差に起因して発生することが、本出願人の実験により確認されている。したがって、そのような領域において、吸入空気量検出手段の検出結果に基づいて算出した基本吸入空気量を、モデル式と相関関係モデルとの間の誤差を補正できる第2補正値を用いて補正することにより、吸入空気量を精度よく算出することができる。請求項12に係る発明は、請求項11に記載の内燃機関3の吸気パラメータ算出装置において、基本吸気パラメータ算出手段は、下流側圧力検出手段の検出結果に基づき、基本吸気パラメータとして、下流側圧力の基本値である基本下流側圧力基本吸気圧PBANを算出する基本下流側圧力算出手段ECU2、ステップ110をさらに有し、補正値算出手段は、圧力比R_Pが臨界圧力比R_CRIT以上のときに、補正値として、基本下流側圧力を補正するための第3補正値補正項PBAERRCORを算出する第3補正値算出手段ECU2、ステップ120~131をさらに有し、吸気パラメータ算出手段は、算出された第3補正値補正項PBAERRCORで基本下流側圧力基本吸気圧PBANを補正することにより、吸気パラメータとして、下流側圧力吸気圧PBAを算出する下流側圧力算出手段ECU2、ステップ111をさらに有することを特徴とする。この内燃機関の吸気パラメータ算出装置によれば、下流側圧力検出手段の検出結果に基づき、基本吸気パラメータとして、下流側圧力の基本値である基本下流側圧力が算出され、圧力比が臨界圧力比以上のときに、補正値として、基本下流側圧力を補正するための第3補正値が算出され、算出された第3補正値で基本下流側圧力を補正することにより、吸気パラメータとして、下流側圧力が算出される。この吸気パラメータ算出装置のように、上流側圧力と、下流側圧力と、開度関数値と、通過空気量との関係を定義したモデル式を用いて、第1開度関数値を算出し、さらに、関数値比を算出した場合、関数値比が表す、モデル式と相関関係モデルとの間の誤差は、圧力比が臨界圧力比以上の領域では、下流側圧力検出手段の検出結果に含まれる誤差に起因して発生することが、本出願人の実験により確認されている。したがって、そのような領域において、下流側圧力検出手段の検出結果に基づいて算出した基本下流側圧力を、モデル式と相関関係モデルとの間の誤差を補正できる第3補正値を用いて補正することにより、下流側圧力を精度よく算出することができる。請求項13に係る発明は、請求項2に記載の内燃機関3の吸気パラメータ算出装置において、吸気通路6の吸気絞り弁よりも上流側に設けられ、吸気通路6内を流れる空気量である吸入空気量GAIRを表す検出信号を出力する吸入空気量検出手段エアフローセンサ20と、上流側圧力大気圧PAを表す検出信号を出力する上流側圧力検出手段大気圧センサ24と、下流側圧力吸気圧PBAを表す検出信号を出力する下流側圧力検出手段吸気圧センサ22と、上流側圧力検出手段および下流側圧力検出手段の検出結果に基づき、下流側圧力と上流側圧力との比である圧力比R_Pを算出する圧力比算出手段ECU2と、吸気絞り弁の開度として弁開度スロットル弁開度THを検出する弁開度検出手段ECU2、スロットル弁開度センサ21と、をさらに備え、基本吸気パラメータ算出手段は、算出された圧力比R_Pが臨界圧力比R_CRIT以上のときに、基本吸気パラメータとして、下流側圧力検出手段の検出結果に基づき、下流側圧力の基本値である基本下流側圧力基本吸気圧PBANを算出する基本下流側圧力算出手段ECU2、ステップ110をさらに有し、補正値算出手段は、圧力比R_Pが臨界圧力比R_CRIT以上のときに、補正値として、基本下流側圧力を補正するための第3補正値補正項PBAERRCORを算出する第3補正値算出手段ECU2、ステップ120~131をさらに有し、吸気パラメータ算出手段は、圧力比R_Pが臨界圧力比R_CRIT以上のときに、算出された第3補正値補正項PBAERRCORで基本下流側圧力基本吸気圧PBANを補正することにより、吸気パラメータとして、下流側圧力吸気圧PBAを算出する下流側圧力算出手段ECU2、ステップ111をさらに有することを特徴とする。この内燃機関の吸気パラメータ算出装置によれば、請求項12に係る発明と同様の作用効果を得ることができる。請求項14に係る発明は、請求項1に記載の内燃機関3の吸気パラメータ算出装置において、吸気通路6の吸気絞り弁よりも上流側に設けられ、吸気通路6内を流れる空気量である吸入空気量GAIRを表す検出信号を出力する吸入空気量検出手段エアフローセンサ20と、上流側圧力大気圧PAを表す検出信号を出力する上流側圧力検出手段大気圧センサ24と、下流側圧力吸気圧PBAを表す検出信号を出力する下流側圧力検出手段吸気圧センサ22と、上流側圧力検出手段および下流側圧力検出手段の検出結果に基づき、下流側圧力と上流側圧力との比である圧力比R_Pを算出する圧力比算出手段ECU2と、吸気絞り弁の開度として弁開度スロットル弁開度THを検出する弁開度検出手段ECU2、スロットル弁開度センサ21と、をさらに備え、基本吸気パラメータ算出手段は、基本吸気パラメータとして、吸入空気量検出手段の検出結果に基づき、吸入空気量GAIRの基本値である基本吸入空気量GAIRNを算出する基本吸入空気量算出手段ECU2、ステップ90をさらに有し、補正値算出手段は、補正値として、検出された弁開度が所定開度THB以上でかつ算出された圧力比R_Pが臨界圧力比R_CRIT未満のときに、基本吸入空気量GAIRNを補正するための第2補正値補正係数KAFMERRを算出する第2補正値算出手段ECU2、ステップ100~107をさらに有し、吸気パラメータ算出手段は、算出された第2補正値補正係数KAFMERRで基本吸入空気量GAIRNを補正することにより、吸気パラメータとして、吸入空気量GAIRを算出する吸入空気量算出手段ECU2、ステップ91をさらに有し、第1開度関数値算出手段は、通過空気量として吸入空気量GAIRを用いて、第1開度関数値を算出することを特徴とする。この内燃機関の吸気パラメータ算出装置によれば、請求項11に係る発明と同様の作用効果を得ることができる。請求項15に係る発明は、請求項14に記載の内燃機関3の吸気パラメータ算出装置において、基本吸気パラメータ算出手段は、算出された圧力比R_Pが臨界圧力比R_CRIT以上のときに、基本吸気パラメータとして、下流側圧力検出手段の検出結果に基づき、下流側圧力の基本値である基本下流側圧力基本吸気圧PBANを算出する基本下流側圧力算出手段ECU2、ステップ110をさらに有し、補正値算出手段は、圧力比R_Pが臨界圧力比R_CRIT以上のときに、補正値として、基本下流側圧力を補正するための第3補正値補正項PBAERRCORを算出する第3補正値算出手段ECU2、ステップ120~131をさらに有し、吸気パラメータ算出手段は、圧力比R_Pが臨界圧力比R_CRIT以上のときに、算出された第3補正値補正項PBAERRCORで基本下流側圧力基本吸気圧PBANを補正することにより、吸気パラメータとして、下流側圧力吸気圧PBAを算出する下流側圧力算出手段ECU2、ステップ111をさらに有することを特徴とする。この内燃機関の吸気パラメータ算出装置によれば、請求項12,13と同様の作用効果を得ることができる。請求項16に係る発明は、請求項1に記載の内燃機関3の吸気パラメータ算出装置において、上流側圧力大気圧PAを表す検出信号を出力する上流側圧力検出手段大気圧センサ24と、上流側圧力検出手段の検出結果に基づき、上流側圧力大気圧PAを算出する上流側圧力算出手段ECU2と、下流側圧力吸気圧PBAを表す検出信号を出力する下流側圧力検出手段吸気圧センサ22と、上流側圧力検出手段および下流側圧力検出手段の検出結果に基づき、下流側圧力と上流側圧力との比である圧力比R_Pを算出する圧力比算出手段ECU2と、をさらに備え、基本吸気パラメータ算出手段は、算出された圧力比R_Pが臨界圧力比R_CRIT以上のときに、基本吸気パラメータとして、下流側圧力検出手段の検出結果に基づき、下流側圧力の基本値である基本下流側圧力基本吸気圧PBANを算出する基本下流側圧力算出手段ECU2、ステップ110をさらに有し、をさらに備え、補正値算出手段は、圧力比R_Pが臨界圧力比R_CRIT以上のときに、補正値として、基本下流側圧力を補正するための第3補正値補正項PBAERRCORを算出する第3補正値算出手段ECU2、ステップ120~131をさらに有し、吸気パラメータ算出手段は、圧力比R_Pが臨界圧力比R_CRIT以上のときに、算出された第3補正値補正項PBAERRCORで基本下流側圧力基本吸気圧PBANを補正することにより、吸気パラメータとして、下流側圧力吸気圧PBAを算出する下流側圧力算出手段ECU2、ステップ111をさらに有することを特徴とする。この内燃機関の吸気パラメータ算出装置によれば、請求項12,13,15と同様の作用効果を得ることができる。請求項17に係る発明は、吸気通路6に設けられた吸気絞り弁スロットル弁7aによって吸気絞り弁を通過する空気量が通過空気量として変更される内燃機関3において、吸気通路6内の空気の状態を表す吸気パラメータ通過空気量GAIRTH、吸入空気量GAIR、吸気圧PBAを算出する内燃機関3の吸気パラメータ算出方法であって、吸気パラメータの基本値として基本吸気パラメータ基本通過空気量GAIRTHN,基本吸入空気量GAIRN,基本吸気圧PBANを算出しステップ6,90,110,144,175、所定のモデリング手法によって導出された、吸気絞り弁の上流側おける吸気通路6内の圧力である上流側圧力大気圧PAと、吸気絞り弁の下流側における吸気通路6内の圧力である下流側圧力吸気圧PBAと、吸気絞り弁の開度によって定まる開度関数値KTHと、通過空気量GAIRTHとの関係を定義したモデル式[式7,13,23]を用いて、開度関数値の第1の算出値として第1開度関数値モデル式値KTHCALを算出しステップ11,103,123,154,183、吸気絞り弁の開度スロットル弁開度THと開度関数値KTHとの相関関係を表す相関関係モデル図2を用いて、開度関数値の第2の算出値として第2開度関数値マップ値KTHを算出しステップ11,100,120,151,180、算出された第1開度関数値および算出された第2開度関数値の一方と他方との比である関数値比関数値誤差KTHERRを用いて、補正値補正係数KTHCOR、補正係数KAFMERR、補正項PBAERRCORを算出しステップ3,100~107,120~131,142,146,147,172、算出された補正値で基本吸気パラメータを補正することにより、吸気パラメータ通過空気量GAIRTH、吸入空気量GAIR、吸気圧PBAを算出するステップ6,91,111,144,175ことを特徴とする。この内燃機関の吸気パラメータ算出方法によれば、請求項1の内燃機関の吸気パラメータ算出装置と同様の作用効果を得ることができる。以下、図面を参照しながら、本発明の第1実施形態に係る内燃機関の吸気パラメータ算出装置について説明する。図1に示すように、第1実施形態の吸気パラメータ算出装置1は、ECU2を備えており、このECU2は、後述するように、吸気パラメータとして、通過空気量GAIRTH、吸入空気量GAIRおよび吸気圧PBAを算出する。エンジン3は、図示しない車両に搭載されたガソリンエンジンタイプのものであり、このエンジン3には、燃料噴射弁4および点火プラグ5が気筒ごとに取り付けられている。この燃料噴射弁4は、ECU2に電気的に接続されており、その開弁時間および開弁タイミングがECU2によって制御され、それにより、燃料噴射制御が実行される。また、点火プラグ5も、ECU2に電気的に接続されており、ECU2によって、点火時期に応じたタイミングで燃焼室内の混合気を燃焼させるように、放電状態が制御される。すなわち、点火時期制御が実行される。エンジン3の吸気通路6には、上流側から順に、エアフローセンサ20、スロットル弁機構7、スロットル弁開度センサ21および吸気圧センサ22が設けられている。このエアフローセンサ20は、吸気通路6のエアフローセンサ20付近を流れる空気の流量以下吸入空気量というを検出して、それを表す検出信号をECU2に出力する。ECU2は、後述するように、このエアフローセンサ20の検出信号に基づいて、吸入空気量GAIRを算出する。この吸入空気量GAIRは、質量流量として算出される。なお、本実施形態では、エアフローセンサ20が吸入空気量検出手段に相当する。スロットル弁機構7は、スロットル弁7aおよびこれを開閉駆動するTHアクチュエータ7bなどを備えている。スロットル弁7aは、吸気通路6の途中に回動自在に設けられており、当該回動に伴う開度の変化によりスロットル弁7aを通過する空気量を変化させる。THアクチュエータ7bは、ECU2に接続されたモータにギヤ機構いずれも図示せずを組み合わせたものであり、ECU2からの制御入力信号によって制御されることにより、スロットル弁7aの開度を変化させる。また、スロットル弁開度センサ21は、例えばポテンショメータなどで構成されており、スロットル弁7aの開度以下スロットル弁開度というTHを検出して、それを表す検出信号を、ECU2に出力する。ECU2は、このスロットル弁開度センサ21の検出信号に基づいて、スロットル弁開度THを算出する。このスロットル弁開度THは、角度 ゚として算出される。なお、本実施形態では、スロットル弁7aが吸気絞り弁に、スロットル弁開度THが弁開度に、スロットル弁開度センサ21が弁開度検出手段にそれぞれ相当する。さらに、吸気圧センサ22は、吸気通路6のスロットル弁7aよりも下流側のサージタンクの部分に配置されており、吸気通路6内の圧力以下吸気圧というを検出して、それを表す検出信号をECU2に出力する。ECU2は、後述するように、この吸気圧センサ22の検出信号に基づいて、吸気圧PBAを算出する。この吸気圧PBAは絶対圧として算出される。なお、本実施形態では、吸気圧センサ22が下流側圧力検出手段に相当し、吸気圧PBAが下流側圧力に相当する。一方、ECU2には、吸気温センサ23、大気圧センサ24、クランク角センサ25および4つの車輪速度センサ261つのみ図示がそれぞれ電気的に接続されている。これらの吸気温センサ23および大気圧センサ24はそれぞれ、吸気通路6内の空気の温度以下吸気温というおよび大気の圧力以下大気圧というを検出して、それを表す検出信号をECU2に出力する。ECU2は、これらの吸気温センサ23および大気圧センサ24の検出信号に基づき、吸気温TAおよび大気圧PAをそれぞれ算出する。この吸気温TAは絶対温度として、大気圧PAは絶対圧としてそれぞれ算出される。なお、本実施形態では、大気圧センサ24が上流側圧力検出手段に、大気圧PAが上流側圧力にそれぞれ相当する。また、クランク角センサ25は、マグネットロータおよびMREピックアップで構成されており、図示しないクランクシャフトの回転に伴い、いずれもパルス信号であるCRK信号およびTDC信号をECU2に出力する。このCRK信号は、所定クランク角例えば2 ゚毎に1パルスが出力され、ECU2は、このCRK信号に基づき、エンジン3の回転数以下エンジン回転数というNEを算出する。また、TDC信号は、各気筒のピストンが吸気行程のTDC位置よりも若干、手前の所定のクランク角位置にあることを表す信号であり、所定クランク角毎に1パルスが出力される。さらに、4つの車輪速センサ26の各々は、対応する車輪の回転速度を検出して、それを表す検出信号をECU2に出力する。ECU2は、これらの車輪速センサの検出信号に基づいて、車速VPや、エンジン始動後の総走行距離DISTなどを算出する。一方、ECU2は、CPU、RAM、ROMおよびI/Oインターフェースいずれも図示せずなどからなるマイクロコンピュータで構成されており、以上の各種のセンサ20~26の検出信号などに基づいて、エンジン回転数NEなどの、エンジン3の運転状態を表す各種の運転状態パラメータを算出するとともに、以下に述べるように、通過空気量GAIRTHの算出処理、吸入空気量GAIRの算出処理および吸気圧PBAの算出処理などの各種の算出処理を実行する。なお、本実施形態では、ECU2が、基本吸気パラメータ算出手段、第1開度関数値算出手段、第2開度関数値算出手段、補正値算出手段、吸気パラメータ算出手段、基本通過空気量算出手段、第1補正値算出手段、通過空気量算出手段、オンボード同定手段、推定上流側圧力算出手段、下流側圧力検出手段、圧力比算出手段、弁開度検出手段、基本吸入空気量算出手段、第2補正値算出手段、吸入空気量算出手段、基本下流側圧力算出手段、第3補正値算出手段、および下流側圧力算出手段に相当する。以下、本実施形態における通過空気量GAIRTHの算出手法の原理について説明する。本実施形態のエンジン3の場合、スロットル弁7aが吸気通路6に設けられているので、このスロットル弁7aを通過する空気以下通過空気というを圧縮性流体かつ断熱流であると見なしかつスロットル弁7aをノズルと見なしてモデリングすると、下式1に示すモデル式が得られる。上式1において、uは通過空気の流速を、P1,P2はスロットル弁7aの上流側および下流側圧力を、ρ1はスロットル弁7aの上流側の吸入空気の密度を、κは吸入空気の比熱比をそれぞれ表している。次に、通過空気の流量をG、スロットル弁7aの開口面積をAth、スロットル弁7aの流量係数をCd、空気の温度をT、空気のガス定数をRとし、連続の式および気体の状態方程式を適用して変形すると、下式2が得られる。ここで、上式2の圧力比P2/P1の項を流量関数値Ψとして下式3のように定義し、これを用いて上式2を書き換えると、下式4が得られる。上式4において、開口面積Athおよび流量係数Cdはいずれも、スロットル弁開度THによって決まる関数であるので、値CdAthを開度関数値KTH=CdAthと定義し、これを用いて上式4を書き換えると、下式5が得られる。上式5において、スロットル弁7aの上流側における圧力損失を無視し、上流側圧力P1を大気圧PAに置き換え、流量Gを通過空気量GAIRTHに、温度Tを吸気温TAにそれぞれ置き換えるとともに、流量関数値ΨをFPBAPA、気体定数Rの平方根をRGAS=R12、吸気温TAの平方根をRTTA=TA12とそれぞれ表記すると、下式6に示すモデル式が得られる。なお、下式6のKCは、通過空気量GAIRTHの単位をg/secに変換するための変換係数である。上式6において、流量関数値FPBAPAは、吸気圧PBAと大気圧PAとの比である圧力比R_P=PBA/PAによって決まる値であるので、後述するように、圧力比R_Pに応じて、マップを検索することにより算出される。また、上式6をKTHについて整理するとともに、通過空気量GAIRTHを吸入空気量GAIRに置き換えると、下式7が得られる。ここで、開度関数値KTHは、スロットル弁開度THによって決まる関数値であり、この開度関数値KTHとスロットル弁開度THとの実際の関係を測定し、マップを作成した場合すなわちオフラインで同定した場合、図2に示すものが得られる。上記式7によって算出した開度関数値をモデル式値KTHCALとし、図2のマップ検索により算出した開度関数値をマップ値KTHとするとともに、モデル式値KTHCALとマップ値KTHとの比を値KTHERR=KTHCAL/KTHとした場合、KTHERR=1のときには、KTH=KTHCALとなり、モデル式6のモデル化誤差がないことになる。一方、値KTHERRが値1に対して乖離している場合、その乖離度合が大きいほど、マップ値KTHに対するモデル式値KTHCALの乖離度合が大きいことになり、前述したモデル式6のモデル化誤差が大きいことになる。以上のように、値KTHERRはモデル式6のモデル化誤差を表すものであり、以下の説明では、値KTHERRを関数値誤差KTHERRという。なお、本実施形態では、図2のマップが相関関係モデルに、モデル式値KTHCALが第1開度関数値に、マップ値KTHが第2開度関数値に、関数値誤差KTHERRが関数値比にそれぞれ相当する。この関数値誤差KTHERRは、主に、標準品基準品に対するスロットル弁7aの開口面積のずれに起因して発生する。このずれは、標準品に基づいて設定されたマップ値KTHに対する実製品のずれであり、例えば、スロットルボア径の精度ばらつきや、スロットルボアに付着したカーボンなどの異物に起因して発生するものである。以下の説明では、このずれを第1誤差という。また、関数値誤差KTHERRの算出において、エアフローセンサ20および吸気圧センサ22の検出信号を用いている関係上、関数値誤差KTHERRは、エアフローセンサ20の検出誤差以下第2誤差というや、吸気圧センサ22の検出誤差以下第3誤差というに起因しても発生する。この第1誤差の場合、スロットル弁開度THが大きくなるほどすなわち、スロットル弁7aの開口面積が大きくなるほど、関数値誤差KTHERRに対する影響は小さくなるので、スロットル弁開度THが所定開度THB以上となる領域では、第1誤差を無視することができる。なお、この所定開度THBは、実験的に求めることも可能であるとともに、スロットル弁7aの開口面積の許容誤差に基づいて演算することも可能である。また、圧力比R_Pが臨界圧力比R_CRIT未満の場合、流量関数値FPBAPAは一定値となるので、R_P<R_CRITの領域では、関数値誤差KTHERRは吸気圧センサ22の検出誤差の影響を受けなくなる。したがって、圧力比R_Pが臨界圧力比R_CRITとなるスロットル弁開度THをTH_CRITとすると、TH<TH_CRITの領域すなわちR_P<R_CRITの領域では、第3誤差を無視することができる。その結果、THB≦TH<TH_CRITの領域すなわちTHB≦THかつR_P<R_CRITが成立する領域では、第1誤差および第3誤差を無視することができるので、関数値誤差KTHERRによって第2誤差を特定することが可能となる。以上の理由により、本実施形態の場合、後述するように、TH<THBの領域では、第1誤差を補償補正するために、関数値誤差KTHERRを用いて、基本通過空気量GAIRTHN補正用の補正係数KTHCORが算出され、THB≦THかつR_P<R_CRITの領域では、第2誤差を補償するために、関数値誤差KTHERRを用いて、基本吸入空気量GAIRN補正用の補正係数KAFMERRが算出され、TH_CRIT≦THの領域すなわちR_CRIT≦R_Pの領域では、第3誤差を補償するために、関数値誤差KTHERRを用いて、基本吸気圧PBAN補正用の補正項PBAERRCORが算出される。以下、上述した第1誤差を補償しながら、通過空気量GAIRTHを算出する手法について説明する。まず、関数値誤差KTHERRから値1を減算した値を誤差KTHERRCOR=KTHERR-1とする。この場合、前述したように、値KTHERRが値1に対して乖離している度合がモデル化誤差を表すので、誤差KTHERRCORは、モデル化誤差を表すものとなり、誤差モデル式を、下式8のように定義する。なお、この誤差モデル式8を用いた理由については後述する。この誤差モデル式8のAはモデルパラメータであり、このモデルパラメータAをオンボード同定する演算手法として、最小2乗法を用いた場合、モデルパラメータAの同定演算式として、下式9が得られる。なお、後述するように、モデルパラメータAの実際の演算では、上式9に基づき、スロットル弁開度THのサンプリング領域が偏った場合でも、モデルパラメータAの同定結果における算出精度の低下を回避できるような演算手法が用いられる。次に、同定したモデルパラメータAを用いて、前述した8により、誤差KTHERRCORを算出した後、下式10により、誤差補正係数KTHCORを算出する。上式10に示すように、誤差補正係数KTHCORは、誤差KTHERRCORと値1との和、すなわち関数値誤差KTHERRの逆数に相当する値として算出される。これは以下の理由による。すなわち、関数値誤差KTHERRがモデル式値KTHCALとマップ値KTHとの比である以上、例えば、KTHERR>1の場合、言い換えればモデル式値KTHCALがマップ値KTHを上回っている場合には、その上回っている分を補正するために、関数値誤差KTHERRで通過空気量GAIRTHを除算すればよいことになる。一方、KTHERR<1の場合、すなわちモデル式値KTHCALがマップ値KTHを下回っている場合には、その下回っている分を補正するために、関数値誤差KTHERRで通過空気量GAIRTHを除算すればよいことになる。これに対して、誤差補正係数KTHCORは、後述するように乗算係数として用いられるので、モデル式値KTHCALがマップ値KTHに対して下回っている分または上回っている分を補正するために、関数値誤差KTHERRの逆数に相当する値として算出される。そして、最終的に、誤差補正係数KTHCORを用いて、下式11,12により、通過空気量GAIRTHが算出される。ここで、上式11のGAIRTHNは、基本通過空気量であり、HPAは、後述するように算出される推定大気圧である。この式11は、前述した式7の通過空気量GAIRTHを基本通過空気量GAIRTHNに、大気圧PAを推定大気圧HPAにそれぞれ置き換えたものに相当する。この式11において、大気圧PAに代えて、推定大気圧HPAを用いたのは、大気圧PAの変動を回避しながら、通過空気量GAIRTHの算出精度を向上させるためである。上式12に示すように、通過空気量GAIRTHは、基本通過空気量GAIRTHNを補正係数KTHCORで補正することにより算出される。なお、本実施形態では、基本通過空気量GAIRTHNが基本吸気パラメータに相当し、補正係数KTHCORが補正値および第1補正値に相当する。次に、前述した誤差モデル式8を用いた理由について説明する。まず、開度関数値KTHは、前述したように、開口面積Athおよび流量係数Cdの積であり、スロットル弁7aの半径の2乗値との相関性が高い以上、誤差モデル式としては、誤差を従属変数とし、スロットル弁開度THを独立変数とする2次式を用いた方が、良好な算出精度が得られることになる。その場合、THB≦THの領域では、関数値誤差KTHERRは、前述したように、スロットル弁開度THの算出誤差ではなく、吸入空気量GAIRまたは吸気圧PBAの算出誤差に起因して発生する関係上、吸入空気量GAIRおよび吸気圧PBAの算出誤差がない条件下では、図3に示すように、関数値誤差KTHERR=1となるので、誤差モデル式のモデルパラメータAを同定する必要がないことになる。これに加えて、図3に示すように、誤差モデル式として、関数値誤差KTHERRを従属変数Yとし、スロットル弁開度THを独立変数Xとする2次式を用いた場合、その誤差モデル式は、Y=aX2bXca,b,cはモデルパラメータの形式となることで、3つのモデルパラメータa,b,cを同定する必要が生じる。その結果、同定演算が複雑になり、演算負荷が高くなってしまう。なお、図3の黒丸は、関数値誤差KTHERRの制御タイミングごとの演算結果を示している。これに対して、上記8式を用いた場合、図4に示すように、誤差KTHERRCORを従属変数Yとし、TH-THBを独立変数Xとすると、誤差モデル式がY=AX2の形式となることで、1つのモデルパラメータAのみを同定すればよくなり、それによって、同定演算も容易になり、その演算負荷を低減できることになる。以上の理由により、本実施形態では、誤差モデル式として式8が用いられる。なお、図4の黒丸は、誤差KTHERRCORの制御タイミングごとの演算結果を示している。以下、図5を参照しながら、通過空気量GAIRTHの算出処理について説明する。この算出処理は、前述した算出手法によって通過空気量GAIRTHを算出するものであり、ECU2によって、CRK信号が所定回数発生するごとすなわち所定クランク角ごとに実行される。なお、以下の説明では、算出された各種の値はすべてECU2のRAM内に記憶されるものとするとともに、今回の制御タイミングにおいてRAMから読み出す値すなわち前回の制御タイミングで算出された値を前回値という。同図に示すように、まず、ステップ1図ではS1と略す。以下同じで、スロットル弁開度THが所定開度THB未満であるか否かを判別する。この判別結果がYESのときには、誤差補正係数KTHCORによって通過空気量GAIRTHを補正する必要があると判定して、ステップ2に進み、前述した式8により、誤差KTHERRCORを算出する。次いで、ステップ3で、前述した式10により、誤差補正係数KTHCORを算出する。一方、ステップ1の判別結果がNOのときには、誤差補正係数KTHCORによって通過空気量GAIRTHを補正する必要がないと判定して、ステップ4に進み、誤差補正係数KTHCORを値1に設定する。以上のステップ3または4に続くステップ5で、圧力比R_PHに応じて、図6に示すマップを検索することにより、流量関数値FPBAPAを算出する。この圧力比R_PHは、吸気圧PBAと前述した推定大気圧HPAとの比PBA/HPAに相当する値である。この場合、圧力比R_PHを用いた理由は、前述した推定大気圧HPAを用いた理由と同じであり、推定大気圧HPAは、RAM内に記憶されている、大気圧推定処理で算出された値を用いる。次に、ステップ6に進み、前述した式11,12により、通過空気量GAIRTHを算出する。この場合、吸気温TAの平方根RTTAおよび開度関数値KTHは、RAM内に記憶されている、後述するモデルパラメータ学習処理において算出された値を用いる。以上のように、ステップ6で、通過空気量GAIRTHを算出した後、本処理を終了する。なお、図5の算出処理で算出された通過空気量GAIRTHは、ECU2によって実行される各種の制御処理において用いられる。例えば、燃料噴射制御処理および点火時期制御処理において、過渡運転条件下で、燃料噴射量や点火時期を算出する際、後述するように算出された吸入空気量GAIRを補正するのに用いられる。以下、図7を参照しながら、モデルパラメータ学習処理について説明する。この処理は、関数値誤差KTHERRを反映させながら、誤差モデルのモデルパラメータAの同定演算を行うもの、言い換えればモデルパラメータAの同定値を学習するものであり、ECU2によって、CRK信号が所定回数発生するごとに実行される。同図に示すように、まず、ステップ10で、学習条件フラグF_LEARNが1であるか否かを判別する。この学習条件フラグF_LEARNは、モデルパラメータAの同定値の学習条件が成立しているか否かを表すものであり、その値は、図8に示す手法により設定される。同図8に示すように、まず、ステップ20で、エンジン始動後の車両の総走行距離DISTが所定値DLEARN未満であるか否かを判別する。この判別結果がYESのときには、モデルパラメータAの同定値の学習条件が成立していると判定して、ステップ21に進み、それを表すために、学習条件フラグF_LEARNを1に設定する。その後、本処理を終了する。一方、ステップ20の判別結果がNOのときには、モデルパラメータAの同定値の学習条件が成立していないと判定して、ステップ22に進み、それを表すために、学習条件フラグF_LEARNを0に設定する。その後、本処理を終了する。図7に戻り、ステップ10の判別結果がNOで、モデルパラメータAの同定値の学習条件が成立していないときには、そのまま本処理を終了する。一方、ステップ10の判別結果がYESで、モデルパラメータAの同定値の学習条件が成立しているときには、ステップ11に進み、スロットル弁開度THに応じて、前述した図2のマップを検索することにより、開度関数値のマップ値KTHを算出する。次いで、ステップ12に進み、流量関数値の同定用値FPBAPAiniを算出する。この同定用値FPBAPAiniは、前述した図6の縦軸の流量関数値FPBAPAを同定用値FPBAPAiniに、横軸の圧力比R_PHを圧力比R_PHiにそれぞれ置き換えたマップを、圧力比R_PHiに応じて検索することにより算出される。この圧力比R_PHiは、吸気圧PBAと初期推定大気圧HPAiniとの比PBA/HPAiniとして算出され、この初期推定大気圧HPAiniは、RAM内に記憶されている、大気圧推定処理で算出された値を用いる。この初期推定大気圧HPAiniを用いた理由については後述する。次に、ステップ13で、吸気温TAに応じて、図示しないマップを検索することにより、吸気温TAの平方根RTTAを算出する。その後、ステップ14で、下式13により、開度関数値のモデル式値KTHCALを算出する。この場合、上式13は、前述した式7において、開度関数値KTHをモデル式値KTHCALに、大気圧PAを推定大気圧の初期値HPAiniに、流量関数値FPBAPAを流量関数値の同定用値FPBAPAiniにそれぞれ置き換えたものに相当する。なお、これらの値HPAini,FPBAPAiniを用いた理由については後述する。ステップ14に続くステップ15で、関数値誤差KTHERRを、モデル式値とマップ値との比KTHCAL/KTHに設定する。次いで、ステップ16に進み、スロットル弁開度THが前述した所定開度THB未満であるか否かを判別する。この判別結果がNOのときには、そのまま本処理を終了する。一方、ステップ16の判別結果がYESのときには、スロットル弁開度THが、前述した第1誤差に起因してモデル化誤差が発生する領域にあり、モデルパラメータAの同定演算を行うべきであると判定して、ステップ17に進み、後述するように、モデルパラメータAの同定演算処理を実行する。その後、本処理を終了する。以上のように、このモデルパラメータ学習処理では、学習条件フラグF_LEARN=1のときに、初期推定大気圧HPAiniおよび流量関数値の同定用値FPBAPAiniを用いて、関数値誤差KTHERRが算出される。このような手法によって、関数値誤差KTHERRを算出するのは以下の理由による。すなわち、本実施形態の場合、大気圧PAの推定値を用いて、関数値誤差KTHERRを算出している関係上、この大気圧の推定値に誤差が含まれていると、前述した第1~第3誤差の補償精度修正精度が低下し、ひいては、後述するように、これらの第1~第3誤差を補償しながら算出された3つの値HGAIRTH,GAIR,PBAを用いて算出される推定大気圧HPAの算出精度推定精度も低下してしまう。この場合、そのような大気圧PAの推定誤差の影響を排除するには、吸気圧センサ22の検出信号から算出された吸気圧PBAが大気圧PAの真値と等しくなり、かつ大気圧PAの真値がほとんど変動しない条件下で、関数値誤差KTHERRを算出する必要がある。したがって、それを達成するために、本実施形態では、学習条件フラグF_LEARNを用いて、学習条件の成立を確認するととともに、初期推定大気圧HPAiniおよび流量関数値の同定用値FPBAPAiniを用いて、関数値誤差KTHERRを算出するように構成されている。次に、図9を参照しながら、モデルパラメータAの同定演算処理について説明する。この処理では、以下に述べるように、スロットル弁開度THの領域を、0≦TH<THN1の第1領域と、THN1≦TH<THN2の第2領域と、THN2≦TH<THN3の第3領域と、THN3≦TH<THBの第4領域とに区分けするとともに、これらの4つの領域に対して均等な重み付け処理を施しながら、前述した同定演算式9に基づく演算手法によって、モデルパラメータAが算出される。この場合、THN1~3は、スロットル弁開度THの所定開度であり、0<THN1<THN2<THN3<THBが成立するように設定される。同図に示すように、まず、ステップ40で、誤差KTHERRCORを、関数値誤差KTHERRと値1との偏差KTHERR-1に設定する。次いで、ステップ41に進み、スロットル弁開度THが所定値THN1未満であるか否かを判別する。この判別結果がYESで、スロットル弁開度THが第1領域にあるときには、ステップ44に進み、領域値nを値1に設定する。一方、ステップ41の判別結果がNOで、THN1≦THのときには、ステップ42に進み、スロットル弁開度THが所定値THN2未満であるか否かを判別する。この判別結果がYESで、スロットル弁開度THが第2領域にあるときには、ステップ45に進み、領域値nを値2に設定する。一方、ステップ42の判別結果がNOで、THN2≦THのときには、ステップ43に進み、スロットル弁開度THが所定値THN3未満であるか否かを判別する。この判別結果がYESで、スロットル弁開度THが第3領域にあるときには、ステップ46に進み、領域値nを値3に設定する。一方、ステップ43の判別結果がNOで、スロットル弁開度THが第4領域にあるときには、ステップ47に進み、領域値nを値4に設定する。以上のように、領域値nは、スロットル弁開度THの4つの領域を表す値として算出される。以上のステップ44~47に続くステップ48で、下式14により、第n領域の積分項XXXX[n]を算出する。上式14において、積分項XXXX[n]は、前述した同定演算式9の分母に相当する値であり、値XXXX[n]zは積分項の前回値を表している。また、積分項XXXX[n]の[]内の値nは、上述した領域値であり、この点は以下の説明においても同様である。すなわち、このステップ48では、例えば、領域値n=1のときには第1領域の積分項XXXX[1]が、n=2のときには第2領域の積分項XXXX[2]がそれぞれ算出される。次いで、ステップ49で、下式15により、第n領域の積分項XXY[n]を算出する。上式15において、積分項XXY[n]は、前述した同定演算式9の分子に相当する値であり、値XXY[n]zは積分項の前回値を表している。ステップ49に続くステップ50で、第n領域のサンプリング数SAMPL[n]を、その前回値と値1の和SAMPL[n]z1に設定する。このサンプリング数SAMPL[n]は、第n領域における積分項のサンプリング数すなわち演算結果数を表している。次いで、ステップ51に進み、下式16により、加重平均値XXXXTTLを算出する。上式16を参照すると明らかなように、加重平均値XXXXTTLは、値TH-THB4の相加平均値を領域ごとに算出するとともに、これらの相加平均値に対して、均等重み付けでの加重平均演算を施すことにより算出される。次に、ステップ52で、下式17により、加重平均値XXYTTLを算出する。上式17を参照すると明らかなように、加重平均値XXYTTLは、値KTHERRCORTH-THB2の相加平均値を領域ごとに算出するとともに、これらの相加平均値に対して、均等重み付けでの加重平均演算を施すことにより算出される。ステップ52に続くステップ53で、下式18により、モデルパラメータAを算出した後、本処理を終了する。以上のように、このモデルパラメータAの同定演算処理では、スロットル弁開度THの領域ごとに算出した値KTHERRCORTH-THB2の4つの相加平均値に対して、均等重み付けでの加重平均演算を施すことにより、加重平均値XXYTTLを算出し、領域ごとに算出した値TH-THB4の4つの相加平均値に対して、均等重み付けでの加重平均演算を施すことにより、加重平均値XXXXTTLを算出し、前者を後者で除算することにより、モデルパラメータAが算出される。この演算手法の利点を、図10および図11を参照しながら説明する。図10は、本実施形態の上述した同定演算手法で算出したモデルパラメータAを用いて、誤差KTHERRCORを演算した場合の演算結果例を示しており、図11は、比較のために、前述した式9により算出したモデルパラメータAを用いて、誤差KTHERRCORを演算した場合の演算結果例を示している。なお、両図中の■で示すデータは演算開始直後のものを、▽で示すデータは■で示すデータよりも演算時間が経過したときのものを、×で示すデータは▽で示すデータよりも演算時間がさらに経過したときのものをそれぞれ表している。両図を比較すると明らかなように、図11の場合、演算時間の経過に伴い、スロットル弁開度THのサンプリングデータの多い領域での演算結果が、モデルパラメータAの演算結果に反映される度合が高くなってしまい、その結果、誤差KTHERRCORの演算結果を示す曲線が変化していることが判る。すなわち、モデルパラメータAの同定精度が低下していることが判る。これに対して、図10に示す本実施形態の手法の場合、演算時間が経過した際、スロットル弁開度THのサンプリングデータがある領域に偏った状態になったとしても、前述した均等重み付けの手法により、その領域での演算結果がモデルパラメータAの演算結果に反映される度合が他の領域と同じようになることで、誤差KTHERRCORの演算結果を示す曲線が変化していないことが判る。すなわち、モデルパラメータAの同定演算結果において、高い算出精度を確保できていることが判る。次に、図12を参照しながら、大気圧推定処理について説明する。この処理は、以下に述べるように、推定大気圧HPAおよびその初期値である初期推定大気圧HPAiniを算出するものであり、ECU2によって、クランキング開始以降、CRK信号が所定回数発生するごとに実行される。なお、本実施形態では、推定大気圧HPAが推定上流側圧力に、初期推定大気圧HPAiniが内燃機関を始動する際に検出された下流側圧力に相当する。同図に示すように、まず、ステップ60で、初期圧算出済みフラグF_FINHPAINIが1であるか否かを判別する。この判別結果がNOのときには、ステップ61に進み、初期推定大気圧HPAiniを算出する。このステップ61では、初期推定大気圧の前回値HPAinizと吸気圧PBAを比較し、両者のうちの大きい方が、初期推定大気圧HPAiniとして設定される。次いで、ステップ62に進み、始動モードフラグF_STMODが1であるか否かを判別する。この始動モードフラグF_STMODは、エンジン3のクランキングが終了するまでの間は1に設定され、クランキングが終了したときに0に設定される。このステップ62の判別結果がYESで、クランキング中のときには、そのまま本処理を終了する。一方、ステップ62の判別結果がNOで、クランキングが終了したときには、初期推定大気圧HPAiniの算出を終了すべきであると判定して、ステップ63に進み、それを表すために、初期圧算出済みフラグF_FINHPAINIを1に設定した後、本処理を終了する。このように、ステップ63で、初期圧算出済みフラグF_FINHPAINIが1に設定されると、前述したステップ60の判別結果がYESとなり、その場合には、ステップ64に進み、推定大気圧HPAの算出処理を以下に述べるように実行する。その後、本処理を終了する。次に、図13を参照しながら、推定大気圧HPAの算出処理について説明する。同図に示すように、まず、ステップ70で、初期設定フラグF_FINHPAINIRが1であるか否かを判別する。今回のループが本処理の1回目の制御タイミングであるときには、ステップ70の判別結果がNOとなり、その場合には、ステップ71に進み、推定大気圧HPAおよび遅延大気圧HPADをいずれも初期推定大気圧HPAiniに設定する。次いで、ステップ72に進み、初期設定フラグF_FINHPAINIRを1に設定した後、本処理を終了する。このように、ステップ72で、初期設定フラグF_FINHPAINIRが1に設定されると、前述したステップ70の判別結果がYESとなり、その場合には、ステップ73に進み、前述した図7のステップ11と同様に、スロットル弁開度THに応じて、前述した図2のマップを検索することにより、開度関数値のマップ値KTHを算出する。次いで、ステップ74に進み、流量関数値FPBAPAを算出する。このステップ74では、圧力比R_PHDを、吸気圧PBAと前述した遅延推定大気圧HPADとの比PBA/HPDAとして算出するとともに、前述した図6のマップにおいて、横軸の圧力比R_PHを圧力比R_PHDに置き換えたマップを用い、このマップを圧力比R_PHDに応じて検索することにより、流量関数値FPBAPAが算出される。次に、ステップ75に進み、前述した図7のステップ13と同様に、吸気温TAに応じて、図示しないマップを検索することにより、吸気温TAの平方根RTTAを算出する。その後、ステップ76で、エンジン回転数NEに応じて、図14に示すマップを検索することにより、回転数補正係数KTHNEを算出する。この回転補正係数KTHNEは、エアフローセンサ20よりも上流側に配置されたエアクリーナ図示せずの圧力損失を補正するための値である。ステップ76に続くステップ77で、下式19により、推定通過空気量HGAIRTHを算出する。この式19は、前述した式6において、右辺の通過空気量GAIRTHを推定通過空気量HGAIRTHに、左辺の大気圧PAを遅延推定大気圧HPADに置き換えたものに相当する。次いで、ステップ78に進み、下式20により、空気量偏差DGAIRを算出する。次に、ステップ79で、車速VPが所定車速VPLよりも大きいか否かを判別する。この判別結果がYESのときには、車両が走行中であると判定して、ステップ80に進み、空気量偏差DGAIRに応じて、図15に示すマップを検索することにより、補正項CORHPAを算出する。一方、ステップ79の判別結果がNOのときには、停車中であると判定して、ステップ81に進み、補正項CORHPAを値0に設定する。以上のステップ80または81に続くステップ82で、下式21により、更新推定大気圧HPACALを算出する。次いで、ステップ83に進み、下式22に示す加重平均演算なまし演算により、推定大気圧HPAを算出する。上式22のCA1は、重み係数であり、0<CA1<1が成立するような所定値に設定されている。また、HAPzは、推定大気圧HPAの前回値である。次に、ステップ84で、上記のように算出した推定大気圧の今回値HPAを遅延推定大気圧HPADに設定した後、本処理を終了する。以上のように、この推定大気圧HPAの算出処理では、更新推定大気圧HPACALが、推定大気圧の前回値HPAzを補正項CORHPAで補正することにより算出されるとともに、補正項CORHPAが空気量偏差DGAIRに応じて算出されるので、更新推定大気圧HPACALは、空気量偏差DGAIRが値0になるように算出される。言い換えれば、更新推定大気圧HPACALは、推定通過空気量HGAIRTHが吸入空気量GAIRに一致するように算出され、そのような更新推定大気圧HPACALと推定大気圧の前回値HPAzとの加重平均演算によって、推定大気圧HPAが算出されるので、推定大気圧HPAを実際の大気圧PAに精度よく追従するように算出することができる。以下、図16を参照しながら、吸入空気量GAIRの算出処理について説明する。この処理は、以下に述べるように、エアフローセンサ20の検出信号の出力電圧値VAFMを用いて、吸入空気量GAIRを算出するものであり、ECU2によって、CRK信号が所定回数発生するごとに実行される。まず、ステップ90で、出力電圧値VAFMに応じて、図17に示すマップを検索することにより、基本吸入空気量GAIRNを算出する。次いで、ステップ91に進み、基本吸入空気量GAIRNを補正係数KAFMERRで割った値GAIRN/KAFMERRを、吸入空気量GAIRとして設定する。この補正係数KAFMERRの算出手法については後述する。以上のようにステップ91で、吸入空気量GAIRを算出した後、本処理を終了する。なお、本実施形態では、基本吸入空気量GAIRNが基本吸気パラメータに相当し、補正係数KAFMERRが補正値および第2補正値に相当する。上記ステップ91において、補正係数KAFMERRを基本吸入空気量GAIRNを除算するための値として用いた理由は、前述した誤差補正係数KTHCORを、関数値誤差KTHERRの逆数として算出した理由と同じである。なお、図16の算出処理で算出された吸入空気量GAIRは、ECU2によって実行される各種の制御処理において用いられる。例えば、燃料噴射制御処理および点火時期制御処理において、燃料噴射量や点火時期を算出する際に用いられる。次に、図18を参照しながら、上記の補正係数KAFMERRの算出処理について説明する。この処理は、前述した関数値誤差KTHERRを用いる手法によって、補正係数KAFMERRを算出するものであり、ECU2によって、CRK信号が所定回数発生するごとに実行される。同図に示すように、まず、ステップ100で、前述した図7のステップ11と同じ手法により、開度関数値のマップ値KTHを算出する。次いで、ステップ101に進み、流量関数値FPBAPAを算出する。この場合、流量関数値FPBAPAは、前述した図6の横軸の圧力比R_PHを圧力比R_P=PBA/PAに置き換えたマップを用い、このマップを圧力比R_Pに応じて検索することにより算出される。次に、ステップ102で、前述した図7のステップ13と同じ手法により、吸気温TAの平方根RTTAを算出する。その後、ステップ103で、下式23により、開度関数値のモデル式値KTHCALを算出する。この式23は、前述した式7において、左辺のKTHをKTHCALに置き換えたものに相当する。ステップ103に続くステップ104で、関数値誤差KTHERRを、モデル式値とマップ値との比KTHCAL/KTHに設定する。次いで、ステップ105に進み、スロットル弁開度THが前述した所定開度THB未満であるか否かを判別する。この判別結果がYESのときには、そのまま本処理を終了する。一方、ステップ105の判別結果がNOのときには、ステップ106に進み、圧力比R_Pが臨界圧力比R_CRIT未満であるか否かを判別する。この判別結果がNOのときには、そのまま本処理を終了する。一方、ステップ106の判別結果がYESのとき、すなわちTHB≦THかつR_CRIT<R_Pのときには、前述した第2誤差に起因して、モデル化誤差が発生する領域にあり、補正係数KAFMERRを更新すべきであると判定して、ステップ107に進み、下式24により、補正係数KAFMERRを算出する。その後、本処理を終了する。上式24のCA2は、重み係数であり、0<CA2<1が成立するような所定値に設定されている。なお、この重み係数CA2を、エンジン回転数NEに応じてマップを検索する手法により算出してもよい。また、式24のKAFMERRzは、補正係数KAFMERRの前回値である。以上のように、図18の算出処理では、補正係数KAFMERRは、その前回値KAFMERRzと関数値誤差KTHERRの加重平均演算により算出されるので、補正係数KAFMERRを、関数値誤差KTHERRが表すモデル化誤差を反映した値として算出することができる。同じ理由により、何らかの原因で、関数値誤差KTHERRの演算結果において演算誤差が一時的に急増した場合でも、その影響を回避しながら、補正係数KAFMERRを精度よく算出することができる。次に、図19を参照しながら、吸気圧PBAの算出処理について説明する。この処理は、以下に述べるように、吸気圧センサ22の検出信号の出力電圧値VPBAを用いて、吸気圧PBAを算出するものであり、ECU2によって、CRK信号が所定回数発生するごとに実行される。まず、ステップ110で、出力電圧値VPBAに基づき、下式25により、基本吸気圧PBANを算出する。なお、下式25のα、β、γは所定値である。次いで、ステップ111に進み、吸気圧PBAを、基本吸気圧PBANと補正項PBAERRCORとの和PBANPBAERRCORに設定する。この補正項PBAERRCORの算出手法については後述する。以上のようにステップ111で、吸気圧PBAを算出した後、本処理を終了する。なお、本実施形態では、基本吸気圧PBANが基本吸気パラメータに相当し、補正項PBAERRCORが補正値および第3補正値に相当する。この場合、図19の算出処理で算出された吸気圧PBAは、ECU2によって実行される各種の制御処理において用いられる。例えば、燃料噴射制御処理および点火時期制御処理において、燃料噴射量や点火時期を算出する際に用いられる。次に、図20を参照しながら、上記の補正項PBAERRCORの算出処理について説明する。この処理は、前述した関数値誤差KTHERRを用いる手法によって、補正項PBAERRCORを算出するものであり、ECU2によって、CRK信号が所定回数発生するごとに実行される。同図を参照すると明らかなように、この処理のステップ120~125は、図18のステップ100~105と同様に構成されているので、以下の説明では、ステップ126以降を中心に説明する。この処理の場合、ステップ125の判別結果がYESで、THB≦THのときには、ステップ126に進み、圧力比R_Pが臨界圧力比R_CRIT以上であるか否かを判別する。この判別結果がNOのときには、そのまま本処理を終了する。一方、ステップ126の判別結果がYESで、R_CRIT≦R_Pのときには、前述した第3誤差に起因して、モデル化誤差が発生する領域にあり、補正項PBAERRCORを更新すべきであると判定して、ステップ127に進み、流量関数値の補正演算用値FPBAPACORを、流量関数値と関数値誤差との比FPBAPA/KTHERRに設定する。この流量関数値の補正演算用値FPBAPACORは、モデル化誤差によって生じている流量関数値FPBAPAの誤差を表す値である。次いで、ステップ128に進み、流量関数値の補正演算用値FPBAPACORに応じて、図21に示すマップを検索することにより、圧力比の補正演算用値R_PCORを算出する。このマップは、前述した図6のマップにおいて、横軸と縦軸の関係を入れ換えるとともに、縦軸を圧力比の補正演算用値R_PCORに、横軸を流量関数値の補正演算用値FPBAPACORにそれぞれ置き換えたものに相当する。以上のように算出される圧力比の補正演算用値R_PCORは、モデル化誤差によって生じている圧力比R_Pの誤差を表す値である。次に、ステップ129で、吸気圧の補正演算用値PBACORを、大気圧PAと圧力比の補正演算用値R_PCORとの積PAR_PCORに設定する。この吸気圧の補正演算用値PBACORは、モデル化誤差によって生じている吸気圧PBAの誤差を表す値である。さらに、ステップ129に続くステップ130で、吸気圧誤差PBAERRを、吸気圧の補正演算用値PBACORと基本吸気圧PBANとの偏差PBACOR-PBANに設定する。次いで、ステップ131に進み、下式26により、補正項PBAERRCORを算出した後、本処理を終了する。上式26のCA3は、重み係数であり、0<CA3<1が成立するような所定値に設定されている。なお、この重み係数CA3を、エンジン回転数NEに応じてマップを検索する手法により算出してもよい。また、式26のPBAERRCORzは、補正項PBAERRCORの前回値である。以上のように、図20の算出処理では、補正項PBAERRCORは、その前回値PBAERRCORzと吸気圧誤差PBAERRの加重平均演算により算出されるので、補正項PBAERRCORを、関数値誤差KTHERRが表すモデル化誤差を反映した値として算出することができる。同じ理由により、何らかの原因で、関数値誤差KTHERRの演算結果において演算誤差が一時的に急増した場合でも、その影響を回避しながら、補正項PBAERRCORを精度よく算出することができる。以上のように、第1実施形態の吸気パラメータ算出装置1によれば、通過空気量GAIRTHを算出する際、関数値誤差KTHERRが、開度関数値のマップ値KTHとモデル式値KTHCALとの比として算出されることにより、モデル式6のモデル化誤差を表すものとして算出される。そして、TH<THBの領域では、そのような関数値誤差KTHERRを用いて、モデルパラメータAがオンボード同定され、オンボード同定されたモデルパラメータAを用いて、補正係数KTHCORが算出されるとともに、その補正係数KTHCORで基本通過空気量GAIRTHNを補正することにより、通過空気量GAIRTHが算出されるので、通過空気量GAIRTHを、モデル式6のモデル化誤差が補正された値として算出することができる。それにより、通過空気量GAIRTHを精度よく算出することができる。また、通過空気量GAIRTHの算出式11,12において、大気圧PAに代えて、推定大気圧HPAを用いているので、大気圧PAが変動しやすい条件下でも、大気圧変動の影響を回避しながら、通過空気量GAIRTHを算出することができ、その算出精度をさらに向上させることができる。さらに、モデルパラメータAがオンボード同定されるので、スロットル弁7aにおける経年変化および個体間のばらつきなどに起因して、誤差モデル式9が誤差KTHERRCORとスロットル弁開度THとの間の実際の関係から乖離した場合、すなわちモデル化誤差を生じた場合でも、オンボード同定されたモデルパラメータAを用いることによって、そのようなモデル化誤差を迅速に補償でき、誤差モデル式9を誤差KTHERRCORとスロットル弁開度THとの実際の関係に迅速に適合させることができる。それにより、補正係数KTHCORによる補正精度を向上させることができ、通過空気量GAIRTHの算出精度をより一層、向上させることができる。また、モデルパラメータAの同定演算では、スロットル弁開度THの領域ごとに算出した値KTHERRCORTH-THB2の4つの相加平均値に対して、均等重み付けでの加重平均演算を施すことにより、加重平均値XXYTTLを算出し、領域ごとに算出した値TH-THB4の4つの相加平均値に対して、均等重み付けでの加重平均演算を施すことにより、加重平均値XXXXTTLを算出し、前者を後者で除算することにより、モデルパラメータAが算出される。それにより、モデルパラメータAの同定演算中、スロットル弁開度THのサンプリングデータがある領域に偏った状態になったとしても、その領域での演算結果がモデルパラメータAの演算結果に反映される度合を他の領域と同じようにすることができ、モデルパラメータAの演算結果において、高い演算精度を確保することができる。さらに、エンジン始動後における車両の総走行距離DISTが所定値DLEARN未満であるときに、大気圧PAに代えて初期推定大気圧HPAiniを用いながら、モデルパラメータAの同定演算が実行され、DIST≧DLEARNのときに、モデルパラメータAの同定演算が禁止される。この場合、エンジン始動後での車両の総走行距離が小さいときには、スロットル弁開度THが低開度域で推移し、TH<THBの状態になる頻度が高いことで、モデルパラメータAの同定演算頻度を高めることができ、その演算精度をさらに向上させることができる。これに加えて、開度関数値のモデル式値KTHCALの算出において、大気圧PAに代えて、初期推定大気圧HPAiniが用いられるので、大気圧PAが変化しやすいような条件下でも、関数値誤差KTHERRを精度よく算出することができ、それにより、モデルパラメータAの演算精度をより一層、向上させることができる。以上のようにモデルパラメータAの演算において、高い演算精度を確保できることで、補正係数KTHCORの算出精度すなわち通過空気量GAIRTHの算出精度をさらに向上させることができる。一方、THB≦THの領域では、モデルパラメータAの同定演算が停止されるので、モデル化誤差が発生しない領域での無用な同定演算を回避することができ、演算負荷を低減することができる。また、吸入空気量GAIRを算出する際、補正係数KAFMERRで基本吸入空気量GAIRNを補正することにより、吸入空気量GAIRが算出され、この補正係数KAFMERRは、THB≦THかつR_P<R_CRITのとき、すなわち前述した第2誤差に起因して、モデル式6のモデル化誤差が発生していると推定されるときに、そのモデル化誤差を表す関数値誤差KTHERRと、補正係数の前回値KAFMERRzとの加重平均演算によって算出される。それにより、補正係数KAFMERRを、関数値誤差KTHERRが表すモデル化誤差を反映した値として算出することができ、吸入空気量GAIRを精度よく算出することができる。さらに、吸気圧PBAを算出する際、補正項PBAERRCORで基本吸気圧PBANを補正することにより、吸気圧PBAが算出され、この補正項PBAERRCORは、R_CRIT≦R_Pのとき、すなわち前述した第3誤差に起因して、モデル式6のモデル化誤差が発生していると推定されるときに、そのモデル化誤差を表す関数値誤差KTHERRを用いて算出されるので、この補正項PBAERRCORを、関数値誤差KTHERRが表すモデル化誤差を反映した値として算出することができ、吸気圧PBAを精度よく算出することができる。また、以上のように補正された吸気圧PBAおよび吸入空気量GAIRと、前述した補正係数KTHCORとを用いて、推定通過空気量HGAIRTHが算出されるので、この推定通過空気量HGAIRTHの算出精度を向上させることができる。さらに、そのような推定通過空気量HGAIRTHから吸入空気量GAIRを減算した値である空気量偏差DGAIRを用いて、補正項CORHPAが算出され、この補正項CORHPAを用いて、推定大気圧HPAが更新されるので、その算出精度を向上させることができる。これに加えて、そのような推定大気圧HPAを用いて通過空気量GAIRTHが算出されるので、通過空気量GAIRTHの算出精度をさらに向上させることができる。なお、第1実施形態では、モデルパラメータAの同定演算手法として、図9のステップ51,52で、前述した式16,17により2つの加重平均値XXXXTTL,XXYTTLをそれぞれ算出し、ステップ53で、前述した式18により、モデルパラメータAを算出する手法を用いたが、これに代えて、以下の同定演算手法を用いてもよい。すなわち、ステップ51,52で、前述した式16,17に代えて、下式27,28を用いて、2つの加重平均値XXXXTTL,XXYTTLをそれぞれ算出した後、ステップ53で前述した式18により、モデルパラメータAを算出するように構成してもよい。上式27,28のKG1~KG4は、重み係数であり、KG1>KG2>KG3>KG4が成立し、かつKG1KG2KG3KG4=1が成立するように設定される。上式27を参照すると明らかなように、加重平均値XXXXTTLは、4つの領域における値TH-THB4の相加平均値に対して、加重平均演算を施すことにより算出されるとともに、その重み係数KG1~KG4は、スロットル弁開度THが小さい領域ほど、より大きな値に設定されている。また、加重平均値XXYTTLも、4つの領域における値KTHERRCORTH-THB2の相加平均値に対して、加重平均演算を施すことにより算出されるとともに、その重み係数KG1~KG4は、加重平均値XXXXTTLと同じ値に設定されている。この場合、関数値誤差KTHERRは、前述したように、標準品基準品に対するスロットル弁7aの開口面積のずれに起因して発生するものであり、そのずれの影響度合は、スロットル弁開度THが小さいほど、より大きくなる。したがって、4つの重み係数KG1~KG4を上記のように設定することによって、そのようなずれの影響度合を反映させながら、モデルパラメータAを同定することができ、それにより、補正係数KTHCORの算出精度を向上させることができる。なお、上記重み係数KG1~KG4は、これらのうちのいずかが互いに同じ値になるように設定してもよい。また、第1実施形態は、誤差モデル式として、Y=AX2の形式のものを用いた例であるが、本発明の誤差モデル式はこれに限らず、Y=aX2bXcの形式のものや、Y=aXbの形式のものを用いてもよい。なお、これらの形式の誤差モデル式を用いた場合、第1実施形態の誤差モデル式8を用いた場合と比較すると、演算精度の向上と演算負荷の低減とを両立させる観点からは、第1実施形態の誤差モデル式8を用いた場合の方が優れている。さらに、第1実施形態は、モデルパラメータAの同定演算を、エンジン始動後の車両の総走行距離DISTが所定値DLEARN未満であるときに実行した例であるが、本発明のモデルパラメータAの同定演算の実行条件はこれに限らず、モデルパラメータAを精度よく同定できる条件であればよい。例えば、始動後のエンジン3の運転時間が所定時間が経過していないことを、モデルパラメータAの同定演算の実行条件としてもよい。一方、第1実施形態は、相関関係モデルとして、図2のマップを用いた例であるが、本発明の相関関係モデルはこれに限らず、吸気絞り弁の開度と開度関数値との相関関係を表すものであればよい。例えば、吸気絞り弁の開度と開度関数値との相関関係を定義した数式を、オフライン同定により決定し、そのような数式を相関関係モデルとして用いてもよい。また、第1実施形態は、吸気絞り弁としてスロットル弁7aを用いた例であるが、本発明の吸気絞り弁はこれに限らず、吸気通路に設けられ、吸気絞り弁を通過する空気量を変更可能な弁であればよい。さらに、第1実施形態は、関数値比として、関数値誤差KTHERRを用いた例であるが、これに代えて、関数値比として関数値誤差KTHERRの逆数を用いてもよい。一方、第1実施形態は、各算出処理をCRK信号が所定回数発生するごとに実行した例であるが、これらの算出処理を、タイマ設定により、所定時間例えば10msecごとに実行するように構成してもよい。次に、本発明の第2実施形態に係る内燃機関の吸気パラメータ算出装置について説明する。この第2実施形態の吸気パラメータ算出装置の場合、第1実施形態の吸気パラメータ算出装置と比較すると、通過空気量GAIRTHの算出処理およびモデルパラメータ学習処理の内容のみが異なっており、これら以外の算出処理は第1実施形態の算出処理と同一であるので、以下、これらの異なる算出処理についてのみ説明する。なお、以下の説明では、第1実施形態と同じ構成については、同じ符号を付すとともに、その説明は省略する。まず、図22を参照しながら、第2実施形態における通過空気量GAIRTHの算出処理について説明する。この算出処理の場合、前述した図5の算出処理と比べて、ステップ145~147のみが図5の算出処理と異なっており、それら以外のステップは図5の算出処理のものと同一であるので、以下、これらのステップ145~147の内容を中心に説明する。同図に示すように、ステップ140の判別結果がNOで、THB≦THのときには、ステップ145で、前述した圧力比R_PHiが臨界圧力比R_CRIT未満であるか否かを判別する。この判別結果がYESのとき、すなわちTHB≦THかつR_PHi<R_CRITが成立しているときには、ステップ147に進み、補正係数KTHCORを、エア誤差KTHERRAFMの逆数1/KTHERRAFMに設定する。このエア誤差KTHERRAFMは、前述した第2誤差を表すものであり、その算出手法については後述する。以上のようにステップ146を実行した後、ステップ143,144を前述したステップ5,6と同様に実行する。その後、本処理を終了する。一方、ステップ145の判別結果がNOで、R_CRIT≦R_PHiのときには、ステップ146に進み、補正係数KTHCORを、吸気圧誤差KTHERRPBAの逆数1/KTHERRPBAに設定する。この吸気圧誤差KTHERRPBAは、前述した第3誤差を表すものであり、その算出手法については後述する。以上のようにステップ147を実行した後、前述したように、ステップ143,144を実行する。その後、本処理を終了する。次に、図23を参照しながら、第2実施形態におけるモデルパラメータ学習処理について説明する。同図のモデルパラメータ学習処理の場合、前述した図7と比べて、ステップ158~160のみが異なっており、それら以外のステップは図7の算出処理のものと同一であるので、以下、これらのステップ158~160の内容を中心に説明する。同図に示すように、ステップ156の判別結果がYESで、TH<THBのときには、ステップ157で、前述した図7のステップ17と同じ手法により、モデルパラメータAの同定演算を実行する。その後、本処理を終了する。一方、ステップ156の判別結果がNOで、THB≦THのときには、ステップ158に進み、前述した圧力比R_PHiが臨界圧力比R_CRIT未満であるか否かを判別する。この判別結果がYESのとき、すなわちTHB≦THかつR_PHi<R_CRITが成立しているときには、ステップ159に進み、下式29により、エア誤差KTHERRAFMを算出した後、本処理を終了する。上式29のCA4は、重み係数であり、0<CA4<1が成立するような所定値に設定されている。なお、この重み係数CA4を、エンジン回転数NEに応じてマップを検索する手法により算出してもよい。また、式29のKTHERRAFMzは、エア誤差KTHERRAFMの前回値である。一方、ステップ158の判別結果がNOで、R_CRIT≦R_PHiのときには、ステップ160に進み、下式30により、吸気圧誤差KTHERRPBAを算出した後、本処理を終了する。上式30のCA5は、重み係数であり、0<CA5<1が成立するような所定値に設定されている。なお、この重み係数CA5を、エンジン回転数NEに応じてマップを検索する手法により算出してもよい。また、式30のKTHERRPBAzは、吸気圧誤差KTHERRPBAの前回値である。以上のように、第2実施形態の吸気パラメータ算出装置によれば、TH<THBのときには、第1実施形態の吸気パラメータ算出装置1と同一の手法により、補正係数KTHCORが算出されるので、第1実施形態の吸気パラメータ算出装置1と同一の作用効果を得ることができる。また、THB≦THかつR_PHi<R_CRITのときには、補正係数KTHCORがエア誤差KTHERRAFMの逆数として算出され、このエア誤差KTHERRAFMが、関数値誤差KTHERRとエア誤差の前回値KTHERRAFMzとの加重平均演算[式29]によって算出されるので、補正係数KTHCORを、前述した第2誤差に起因して発生していると推定されるモデル化誤差が反映された値として算出することができる。同じ理由により、何らかの原因で、関数値誤差KTHERRの演算結果において演算誤差が一時的に増大した場合でも、その影響を回避しながら、補正係数KTHCORを精度よく算出することができる。その結果、通過空気量GAIRTHを精度よく算出することができる。さらに、R_CRIT≦R_PHiのときには、補正係数KTHCORが吸気圧誤差KTHERRPBAの逆数として算出され、この吸気圧誤差KTHERRPBAが、関数値誤差KTHERRと吸気圧誤差の前回値KTHERRPBAzとの加重平均演算[式30]によって算出されるので、補正係数KTHCORを、前述した第3誤差に起因して発生していると推定されるモデル化誤差が反映された値として算出することができる。同じ理由により、何らかの原因で、関数値誤差KTHERRの演算結果において演算誤差が一時的に増大した場合でも、その影響を回避しながら、補正係数KTHCORを精度よく算出することができる。その結果、通過空気量GAIRTHを精度よく算出することができる。次に、第3実施形態に係る内燃機関の吸気パラメータ算出装置について説明する。この第3実施形態の吸気パラメータ算出装置の場合、第1実施形態の吸気パラメータ算出装置と比較すると、通過空気量GAIRTHの算出処理の内容が異なっている点と、前述した図7のモデルパラメータ学習処理に代えて、後述するマップ更新処理を行う点とが異なっており、これらの点以外は第1実施形態の吸気パラメータ算出装置の構成と同一であるので、以下、これらの異なる点についてのみ説明する。なお、以下の説明では、第1実施形態と同じ構成については、同じ符号を付すとともに、その説明は省略する。また、本実施形態では、ECU2が、相加平均値算出手段および記憶手段に相当する。まず、図24を参照しながら、第3実施形態における通過空気量GAIRTHの算出処理について説明する。この算出処理の場合、前述した図5の算出処理と比較すると、ステップ170,172~175は、図5のステップ1,3~6と同一であり、ステップ171のみが図5の算出処理と異なっているので、以下、このステップ171の内容を中心に説明する。図24に示すように、ステップ170の判別結果がYESで、TH<THBのときには、ステップ171に進み、スロットル弁開度THに応じて、図25に示すマップを検索することにより、誤差KTHERRCORを算出する。同図に示すように、このマップでは、スロットル弁開度THの値1単位は角度から所定値THBまでの領域において、スロットル弁開度THの1 ゚ごとの値に対応して、誤差KTHERRCORのマップ値KTHERRCOR[1]~KTHERRCOR[THB-1]および値1が設定されている。これらのマップ値のうち、マップ値KTHERRCOR[1]~KTHERRCOR[THB-1]は、後述するマップ更新処理において更新されるとともに、TH=THBのときのマップ値は更新されず、常に値1に保持される。その理由は、前述したように、THB≦THの領域では、前述した第1誤差に起因するモデル化誤差が、発生しないことによる。ステップ171で、以上の手法により誤差KTHERRCORを算出した後、ステップ172で、前述した図5のステップ3と同じ手法により、補正係数KTHCORを算出する。次いで、ステップ174,175で、前述した図5のステップ5,6と同じ手法により、流量関数値FPBAPAおよび通過空気量GAIRTHをそれぞれ算出する。その後、本処理を終了する。次に、図26を参照しながら、マップ更新処理について説明する。このマップ更新処理は、前述した図25のマップ内のマップ値KTHERRCOR[1]~KTHERRCOR[THB-1]を更新するものであり、ECU2によって、CRK信号が所定回数発生するごとに実行される。同図の算出処理の場合、前述した図7のモデルパラメータ学習処理と比較すると、ステップ180,182~185は、図7のステップ11,13~16と同一であり、ステップ181186の内容のみが図7の処理と異なっているので、以下、これらのステップ181,186の内容を中心に説明する。同図に示すように、ステップ180で、前述した図7のステップ11と同じ手法によりマップ値KTHを算出した後、ステップ181で、圧力比R_PHに応じて、前述した図6のマップを検索することにより、流量関数値FPBAPAを算出する。次いで、ステップ182~185を、図7のステップ13~16と同一の手法で実行する。そして、ステップ185の判別結果がYESで、TH<THBのときには、ステップ186に進み、誤差KTHERRCORのマップ値の更新処理を実行する。この更新処理は、前述したマップ値KTHERRCOR[1]~KTHERRCOR[THB-1]のいずれかを更新するものであり、以下に述べるように実行される。まず、スロットル弁開度THの今回の検出値が、前述したマップ値KTHERRCOR[1]~KTHERRCOR[THB-1]の設定開度1 ゚~THB-1 ゚のうちのいずれに最も近いかを判定し、その判定結果に基づいて更新するマップ値を決定する。例えば、今回の検出値が0.3 ゚の場合には、マップ値KTHERRCOR[1]の更新が決定され、今回の検出値が1.8 ゚の場合には、マップ値KTHERRCOR[2]の更新が決定される。以下、マップ値KTHERRCOR[1]の更新が決定された場合を例にとって説明する。すなわち、マップ値KTHERRCOR[1]の更新が決定されると、上記ステップ184で算出した関数値誤差KTHERRから値1を減算した値KTHERR-1を、マップ値KTHERRCOR[1]の今回値KTHERRCOR[1]_TMPとし、これにRAM内に記憶されているマップ値KTHERRCOR[1]の積分値の前回値ΣKTHERRCOR[1]zを加算することにより、積分値の今回値ΣKTHERRCOR[1]を算出する。そして、算出された積分値の今回値ΣKTHERRCOR[1]をマップ値KTHERRCOR[1]のサンプリング個数で除算した値を、マップ値KTHERRCOR[1]としてRAM内のマップに記憶する。すなわち、マップ値KTHERRCOR[1]は、今回の制御タイミングまでにサンプリング算出されたマップ値KTHERRCOR[1]の相加平均値に更新される。ステップ186で、以上のようにマップ値の更新処理を実行した後、本処理を終了する。以上のように、第3実施形態の吸気パラメータ算出装置によれば、通過空気量GAIRTHの算出処理では、誤差KTHERRCORが、スロットル弁開度THに応じて、図25のマップを検索することにより算出され、この誤差KTHERRCORを用いて補正係数KTHCORが算出される。このマップ内のマップ値KTHERRCOR[n]は、その更新処理において、誤差KTHERRCORの相加平均値が、領域値nに対応して算出され、算出された相加平均値が第n領域のマップ値KTHERRCOR[n]として記憶されるので、マップ値の更新処理が進行するほど、その算出精度を高めることができ、通過空気量GAIRTHの算出精度を向上させることができる。なお、第3実施形態は、誤差KTHERRCOR検索用のマップを用いた例であるが、これに代えて、補正係数KTHCOR検索用のマップを用いるとともに、ステップ186のマップ値更新処理において、補正係数KTHCORのマップ値を更新するように構成してもよい。このように構成した場合でも、第3実施形態の吸気パラメータ算出装置と同じ作用効果を得ることができる。また、第3実施形態は、前述したマップ値KTHERRCOR[n]の設定間隔を1 ゚に設定した例であるが、マップ値の設定間隔はこれに限らず、1 ゚より大きい値や、1 ゚未満の値としてもよい。"