(57)【特許請求の範囲】
【請求項１】
吸気を圧縮して燃焼室に供給するコンプレッサを備える過給エンジンに適用されて、前記コンプレッサによる圧縮後の吸気の温度を推定する過給エンジンの吸気温推定装置において、
前記コンプレッサによる圧縮前の吸気の温度である圧縮前吸気温、同コンプレッサによる圧縮前の吸気の圧力である圧縮前吸気圧、及び同コンプレッサによる圧縮後の吸気の圧力である圧縮後吸気圧に基づき、前記コンプレッサによる吸気の圧縮が断熱された状態で行われたとしたときの同コンプレッサによる圧縮後の吸気温を断熱圧縮吸気温として演算する断熱圧縮吸気温演算部と、
前記断熱圧縮吸気温に対して一次遅れの値を返す関数を用いて、前記コンプレッサによる圧縮後の吸気温の推定値である推定吸気温を演算する推定吸気温演算部と、を備え、
前記断熱圧縮吸気温演算部による前記断熱圧縮吸気温の演算値をＴ１ｍ、前記推定吸気温の更新前の値をＴ１［ｏｌｄ］、前記推定吸気温の更新後の値をＴ１［ｎｅｗ］とし、Ｋを１以上の値を取る係数としたとき、
前記推定吸気温演算部は、下式に示される関係となるように前記推定吸気温の値を更新することで同推定吸気温の演算を行うとともに、前記圧縮後吸気圧の単位時間あたりの変化量が大きいときには、同変化量が小さいときよりも、前記係数の値を同変化量に応じて小さく可変設定する
ことを特徴とする過給エンジンの吸気温推定装置。
【数１】
0006376152.tif 000005
【請求項２】
前記推定吸気温演算部は、前記圧縮後吸気圧の単位時間あたりの変化量の減少に応じて前記係数の値が増加側に変化するときに、同係数の単位時間あたりの変化量を規定のガード値以下に制限するガード処理を行って前記推定吸気温の演算を行う
請求項１に記載の過給エンジンの吸気温推定装置。
【請求項３】
前記推定吸気温演算部は、前記圧縮後吸気圧の単位時間あたりの変化量が大きいときには、同変化量が小さいときよりも、前記ガード値を小さい値に設定して前記ガード処理を行う
請求項２に記載の過給エンジンの吸気温推定装置。
【請求項４】
前記推定吸気温演算部は、前記ガード処理による前記係数の単位時間あたりの変化量の制限を、前記断熱圧縮吸気温の上昇中には行わず、同断熱圧縮吸気温の低下中にのみ行う
請求項２又は３に記載の過給エンジンの吸気温推定装置。
【請求項５】
前記推定吸気温演算部は、前記圧縮後吸気圧の単位時間あたりの変化量が規定の値以上のときに、同変化量に応じた前記係数の可変設定を行う
請求項１〜４のいずれか１項に記載の過給エンジンの吸気温推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、吸気を圧縮して燃焼室に供給するコンプレッサを備える過給エンジンに適用されて、コンプレッサによる圧縮後の吸気の温度を推定する過給エンジンの吸気温推定装置に関する。
【背景技術】
【０００２】
車載等のエンジンの多くでは、吸気通路内の吸気の体積流量を検出するエアフローメータと、吸気通路に取り込まれた外気の温度を検出する外気温センサとが設けられている。そして、エアフローメータ及び外気温センサの検出結果から、燃焼室に導入された吸気の質量、すなわちシリンダ流入空気量を求めている。
【０００３】
一方、吸気を圧縮して燃焼室に供給するコンプレッサを備えた過給エンジンでは、コンプレッサの圧縮により、吸気温が上昇する。また、吸気温が上がれば、吸気の密度は低下する。そのため、過給エンジンでは、シリンダ流入空気量を正確に見積もるため、コンプレッサによる圧縮後の吸気温を確認する必要がある。
【０００４】
従来、過給エンジンでの、コンプレッサによる圧縮後の吸気温を推定する吸気温推定装置として、特許文献１に記載の装置が知られている。同文献の装置では、コンプレッサによる吸気の圧縮が断熱された状態で行われたとするときの圧縮後の吸気温の理論値である断熱圧縮吸気温を、コンプレッサによる圧縮前の吸気温及び吸気圧、コンプレッサによる圧縮後の吸気圧から求めている。そして、その断熱圧縮吸気温に対して、圧縮後吸気圧及びエンジン回転数に基づき設定された補正係数により、エンジン負荷及びエンジン回転数の影響による吸気温変化に対応した補正を行うことで、圧縮後の吸気温の推定値（推定吸気温）を求めている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平０３−２２９９５２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
ところで、実際のコンプレッサでの吸気の圧縮は、断熱的でなく、圧縮後吸気圧の変化に対する吸気温の変化には、外部（吸気管など）との熱交換等による応答遅れがある。そのため、圧縮後吸気圧が変化する過渡時には、上記断熱圧縮吸気温は、実際の吸気温から乖離する。しかしながら、上記従来の過給エンジンの吸気温推定装置では、そうした外部との熱交換の影響は考慮されておらず、過渡時における吸気温の推定結果に誤差が生じる虞がある。
【０００７】
尤も、過給エンジンの吸気系における吸気挙動の精密な物理モデルを構築し、その物理モデルを用いて吸気挙動を再現すれば、コンプレッサによる圧縮後の吸気の温度を高い精度で推定することが可能ではある。しかしながら、そうした物理モデルの構築には、非常に多くの手間が必要であり、また、そうした物理モデルを用いた吸気温の推定には、膨大な演算が必要となる。
【０００８】
本発明は、こうした実情に鑑みてなされたものであり、その解決しようとする課題は、コンプレッサによる圧縮後の吸気温を、簡易な演算で精度良く推定することのできる過給エンジンの吸気温推定装置を提供することにある。
【課題を解決するための手段】
【０００９】
上記課題を解決する過給エンジンの吸気温推定装置は、吸気を圧縮して燃焼室に供給するコンプレッサを備える過給エンジンに適用されて、コンプレッサによる圧縮後の吸気の温度を推定する。同吸気温推定装置は、コンプレッサによる圧縮前の吸気の温度である圧縮前吸気温、同コンプレッサによる圧縮前の吸気の圧力である圧縮前吸気圧、及び同コンプレッサによる圧縮後の吸気の圧力である圧縮後吸気圧に基づき、コンプレッサによる吸気の圧縮が断熱された状態で行われたとしたときの同コンプレッサによる圧縮後の吸気温を断熱圧縮吸気温として演算する断熱圧縮吸気温演算部を備える。
【００１０】
コンプレッサによる圧縮後の吸気圧の変化に対して、実際の圧縮後の吸気温は、吸気管などの外部との熱交換の影響のため、断熱圧縮吸気温に遅れて変化する。外部と吸気との温度差が一定とすると、外部と吸気との間の熱移動の量は、時間に比例して多くなる。そのため、圧縮後吸気圧の変化の量が同じであっても、その変化が短時間で行われた場合には、同変化が長い時間を掛けて行われた場合よりも、断熱圧縮吸気温の変化に対する実際の圧縮後吸気温の遅れは小さくなる。すなわち、断熱圧縮吸気温に対する圧縮後吸気温の遅れは、圧縮後吸気圧の変化が急激なときには小さくなり、圧縮後吸気圧の変化が緩やかなときには大きくなり、一次遅れ系として簡潔に表すことはできないものとなる。
【００１１】
一方、上記過給エンジンの吸気温推定装置は、断熱圧縮吸気温に対して一次遅れの値を返す関数を用いて、コンプレッサによる圧縮後の吸気温の推定値である推定吸気温を演算する推定吸気温演算部を備える。上記関数をそのまま演算に用いれば、演算された推定吸気温は、断熱圧縮吸気温に対する一次遅れ系となるため、過渡時には、実際の圧縮後吸気温から乖離した値となってしまう。
ちなみに、上記過給エンジンの吸気温推定装置における推定吸気温演算部による推定吸気量の演算は、下式に示される関係となるように推定吸気温の値を更新することで行うことができる。下式におけるＴ１ｍは、断熱圧縮吸気温演算部による断熱圧縮吸気温の演算値を示し、Ｔ１［ｏｌｄ］及びＴ１［ｎｅｗ］は、推定吸気温の更新前の値、更新後の値をそれぞれ示している。さらに、Ｋは、１以上の値を取る係数となっている。
【数１】
0006376152.tif 000002
【００１２】
その点、上記過給エンジンの吸気温推定装置における推定吸気温演算部は、圧縮後吸気圧の単位時間あたりの変化量が大きいときには、同変化量が小さいときよりも、上記関数の係数の値を同変化量に応じて小さく可変設定するようにしている。こうして演算された推定吸気温は、圧縮後吸気圧の単位時間あたりの変化量が大きいとき、すなわち同圧縮後吸気圧の変化が急なときには、断熱圧縮吸気温に対する遅れが小さくなる。また、圧縮後吸気圧の単位時間あたりの変化量が小さいとき、すなわち同圧縮後吸気圧の変化が緩やかなときには、断熱圧縮吸気温に対する遅れが大きくなる。したがって、こうした推定吸気温演算部により演算される推定吸気温は、上述した実際の圧縮後吸気温と同じように、圧縮後吸気圧の変化の緩急によって断熱圧縮吸気温に対する遅れ方が変化するものとなる。しかも、推定吸気温演算部が演算に用いる関数そのものは、一次遅れ系を表す単純なものであり、その関数の係数の値を、圧縮後吸気圧の変化の緩急に合わせて変更しているだけである。したがって、上記過給エンジンの吸気温推定装置によれば、コンプレッサによる圧縮後の吸気温を、簡易な演算で精度良く推定することができる。
【００１３】
ところで、上述のように圧縮後吸気温は、圧縮後吸気圧に遅れて変化するため、圧縮後吸気圧の変化の収束後も、その圧縮後吸気圧の変化に応じた圧縮後吸気温の変化は続いている。そのため、圧縮後吸気圧の単位時間あたりの変化量が減少した時点で、断熱圧縮吸気温に対する推定吸気温の追従性を直ちに低下させるようにすると、実際の圧縮後吸気温の変化よりも推定吸気温の変化が遅れてしまうことがある。これに対しては、上記過給エンジンの吸気温推定装置における推定吸気温演算部が、圧縮後吸気圧の単位時間あたりの変化量の減少に応じて係数の値が増加側に変化するときに、同係数の値の単位時間あたりの変化量を規定のガード値以下に制限するガード処理を行って推定吸気温の演算を行うようにするとよい。こうしたガード処理を行えば、圧縮後吸気圧の急激な変化が収束した後もしばらくは、断熱圧縮吸気温に対する推定吸気温の追従性が高い状態が続くため、上記のような推定吸気温の遅れを抑えることが可能となる。
【００１４】
なお、上記のようなガード処理の結果、圧縮後吸気圧の変化が収束した後にも、係数の値の変化が制限された状態が続けられると、推定吸気温の演算値が、実際の圧縮後吸気温よりも早く、圧縮後吸気圧の変化後の値に収束してしまうことがある。これに対しては、上記推定吸気温演算部が、圧縮後吸気圧の単位時間あたりの変化量が大きいときには、同変化量が小さいときよりも、ガード値を小さい値に設定して同ガード処理を行うようにすればよい。このようにすれば、圧縮後吸気圧の変化が収束に向かい、その単位時間あたりの変化量が減少するに従って、ガード処理による係数の変化量の制限が緩和される。そのため、圧縮後吸気圧の変化の収束に応じてより速やかに、ガード処理による係数の変化量の制限を解除することができる。
【００１５】
一方、一般に、過給エンジンの加速時における圧縮後吸気圧の上昇には、同過給エンジンの減速時における圧縮後吸気圧の低下よりも時間がかかるため、圧縮後吸気温が上昇する圧縮後吸気圧の上昇時には、同圧縮後吸気圧の急激な変化は生じにくい。そのため、上記ガード処理による上記係数の単位時間あたりの変化量の制限は、断熱圧縮吸気温の上昇中には行わず、同断熱圧縮吸気温の低下中にのみ行うようにしてもよい。
【００１６】
なお、過給エンジンの加速時における圧縮後吸気圧の上昇や、減速時における同圧縮後吸気圧の低下は、単調な上昇、減少とはならず、上昇の過程において一時的な圧縮後吸気圧の低下が生じたり、低下の過程において一時的な圧縮後吸気圧の上昇が生じたりすることがある。こうした変化過程での圧縮後吸気圧の一時的な変動に対してまで、上記係数の値を変更すると、その一時的な変動に対して鋭敏に反応して、推定吸気温の演算値に不安定な変動が現れることがある。これに対しては、上記過給エンジンの吸気温推定装置において推定吸気温演算部が、圧縮後吸気圧の単位時間あたりの変化量に応じた上記係数の可変設定を、同変化量が規定の値以上のときに行うようにするとよい。このようにすれば、圧縮後吸気圧の単位時間あたりの変化量が規定の値を超えるまでは、同変化量に対して係数が変化しない不感帯が設定されるため、上記のような一時的な圧縮後吸気圧の変動の影響が、推定吸気温の演算値に現れにくくなる。
【図面の簡単な説明】
【００１９】
【図１】過給エンジンの吸気温推定装置の一実施形態の構成と、過給エンジンの構成とを模式的に示す図。
【図２】同実施形態での推定吸気温の演算に使用する係数のマップ値と、圧縮後吸気圧の変化率との関係を示すグラフ。
【図３】上記ガード処理において設定されるガード値と、圧縮後吸気圧の変化率との関係を示すグラフ。
【図４】上記実施形態における推定吸気温の演算に係る処理の手順を示すフローチャート。
【図５】同実施形態における圧縮後吸気圧の低下時の、変化率、係数、マップ値、断熱圧縮吸気温、実際の圧縮後吸気温等の推移を模式的に示すタイムチャート。
【図６】過給エンジンが加速後に減速したときの圧縮後吸気圧の推移を示すグラフ。
【図７】圧縮後吸気圧の上昇時における上記実施形態の推定吸気温の演算結果を、実際の吸気温等と比較して示すグラフ。
【図８】圧縮後吸気圧の低下時における上記実施形態の推定吸気温の演算結果を、実際の吸気温等と比較して示すグラフ。
【発明を実施するための形態】
【００２０】
以下、過給エンジンの吸気温推定装置の一実施形態を、図１〜図８を参照して詳細に説明する。
図１に示すように、本実施形態の吸気温推定装置２０が適用される過給エンジン１０の吸気通路１１には、吸気通路１１に取り込まれた吸気に含まれる塵などの不純物を濾過するエアクリーナ１２が設けられている。吸気通路１１におけるエアクリーナ１２よりも下流側の部分には、コンプレッサ１３が設けられている。コンプレッサ１３は、過給エンジン１０の排気の流勢によるタービンホイールの回転、或いは過給エンジン１０の出力軸であるクランクシャフトの回転により駆動されて、吸気を圧縮して、吸気通路１１の下流側に吐出する。さらに、吸気通路１１におけるコンプレッサ１３よりも下流側の部分には、スロットルバルブ１４が設けられている。スロットルバルブ１４は、吸気通路１１の流路断面積を変化させて、同吸気通路１１を流れる吸気の体積流量（以下、吸入空気流量と記載する）を調整する。こうした吸気通路１１の下流端は、過給エンジン１０の燃焼室１５に接続されている。
【００２１】
また、吸気通路１１におけるエアクリーナ１２とコンプレッサ１３との間の部分には、エアフローメータ１６、温度センサ１７、及び第１圧力センサ１８の３つのセンサが設けられている。また、吸気通路１１におけるコンプレッサ１３とスロットルバルブ１４との間の部分には、第２圧力センサ１９が設けられている。エアフローメータ１６は、上述の吸入空気流量を検出し、温度センサ１７は、コンプレッサ１３による圧縮前の吸気の温度（以下、圧縮前吸気温Ｔ０と記載する）を検出する。また、第１圧力センサ１８は、コンプレッサ１３による圧縮前の吸気の圧力（以下、圧縮前吸気圧Ｐ０と記載する）を、第２圧力センサ１９は、コンプレッサ１３による圧縮後の吸気の圧力（以下、圧縮後吸気圧Ｐ１と記載する）を、それぞれ検出する。
【００２２】
本実施形態の吸気温推定装置２０は、過給エンジン１０の制御を司る電子制御ユニット２１に設けられている。吸気温推定装置２０は、上記温度センサ１７、第１圧力センサ１８、及び第２圧力センサ１９の検出結果から、コンプレッサ１３による圧縮後の吸気の温度（以下、圧縮後吸気温と記載する）を推定する。そして、電子制御ユニット２１は、吸気温推定装置２０が推定した圧縮後吸気温を、過給エンジン１０の各種制御に使用する。例えば、電子制御ユニット２１は、エアフローメータ１６の検出した吸入空気流量と、吸気温推定装置２０による圧縮後吸気温の推定値（以下、推定吸気温Ｔ１と記載する）とに基づいて、燃焼室１５に導入された吸気の質量（シリンダ流入空気質量）を求めている。そして、電子制御ユニット２１は、シリンダ流入空気質量に基づき、同シリンダ流入空気質量に対して適切な重量比率となるように、燃焼室１５に供給する燃料の質量を決定して、過給エンジン１０の燃料供給量を制御している。
【００２３】
吸気温推定装置２０は、断熱圧縮吸気温演算部２２と、推定吸気温演算部２３と、の２つの演算部を備える。吸気温推定装置２０は、これら２つの演算部（２２，２３）が行う処理を通じて、規定の演算周期ごとに、推定吸気温Ｔ１を演算する。
【００２４】
断熱圧縮吸気温演算部２２は、断熱圧縮吸気温Ｔ１ｍの演算を行う。断熱圧縮吸気温Ｔ１ｍは、コンプレッサ１３による吸気の圧縮が断熱された状態で準静的に行われたとしたときの、すなわち圧縮後吸気圧Ｐ１の変化に対する圧縮後吸気温の変化が完全な断熱変化であるとしたときの、圧縮後吸気温の理論的な値である。こうした断熱圧縮吸気温Ｔ１ｍの演算は、圧縮前吸気温Ｔ０、圧縮前吸気圧Ｐ０、及び圧縮後吸気圧Ｐ１の各検出値に基づき、下式（１）に示される関係となるように行われる。なお、下式（１）におけるＡ、Ｂは、空気の物理特性に応じて値が定められた定数である。
【００２５】
【数２】
0006376152.tif 000003
推定吸気温演算部２３は、こうして演算された断熱圧縮吸気温Ｔ１ｍに基づき、推定吸気温Ｔ１を演算する。推定吸気温Ｔ１の演算は、推定吸気温Ｔ１の値を、下式（２）に示される関係となるように、上記演算周期ごとに更新することで行われる。下式（２）におけるＴ１［ｏｌｄ］は、推定吸気温Ｔ１の更新前の値であり、Ｔ１［ｎｅｗ］は、前記推定吸気温の更新後の値である。また、下式（２）におけるＫは、１以上の値を取る係数である。
【００２６】
【数３】
0006376152.tif 000004
上式（２）の関係によれば、更新ごとに推定吸気温Ｔ１の値は、断熱圧縮吸気温Ｔ１ｍに対する推定吸気温Ｔ１の更新前の値（Ｔ１［ｏｌｄ］）の差を係数Ｋで除算した値ずつ加増される。なお、過給エンジン１０の始動後における最初の推定吸気温Ｔ１の演算に際しては、断熱圧縮吸気温Ｔ１ｍの値がそのまま推定吸気温Ｔ１の値に設定される。
【００２７】
上式（２）において、係数Ｋの値を「１」とすると、推定吸気温Ｔ１の更新後の値は、更新前の値に関わらず、断熱圧縮吸気温Ｔ１ｍと常に同じ値となる。すなわち、このときの断熱圧縮吸気温Ｔ１ｍに対する推定吸気温Ｔ１の演算値の遅れ（時定数）はゼロとなる。一方、係数Ｋの値を１から大きくしていくと、断熱圧縮吸気温Ｔ１ｍに対する推定吸気温Ｔ１の演算値の遅れ（時定数）は次第に大きくなる。このように上式（２）において係数Ｋは、断熱圧縮吸気温Ｔ１ｍに対する推定吸気温Ｔ１の追従性を定める係数となっている。
【００２８】
こうした係数Ｋを定数とすると、上式（２）は、断熱圧縮吸気温Ｔ１ｍに対して一次遅れの値を推定吸気温Ｔ１として返す関数Ｆとなる。したがって、推定吸気温演算部２３による推定吸気温Ｔ１の演算は、断熱圧縮吸気温Ｔ１ｍに対して一次遅れの値を返す関数Ｆを用いて行われることになる。
【００２９】
ただし、本実施形態では、係数Ｋは、定数ではなく、圧縮後吸気圧Ｐ１の単位時間あたりの変化量に応じて値が変化する可変値とされている。係数Ｋの値の設定は、下記の（ａ）〜（ｃ）の各処理を通じて行われる。すなわち、推定吸気温演算部２３は、係数Ｋの値の設定に際してまず、（ａ）圧縮後吸気圧Ｐ１の変化率ΔＰ１を算出する。続いて、推定吸気温演算部２３は、（ｂ）その変化率ΔＰ１に基づき、係数Ｋのマップ値Ｋｍａｐを算出する。そして、推定吸気温演算部２３は、（ｃ）マップ値Ｋｍａｐに基づき、係数Ｋの値を設定する。なお、このときの係数Ｋの値の設定では、（ｄ）同係数Ｋの値の単位時間あたりの変化量を制限するガード処理が必要に応じて実施される。
【００３０】
上記（ａ）の処理では、今回の演算周期に取得した圧縮後吸気圧Ｐ１の検出値から前回の演算周期に取得した圧縮後吸気圧Ｐ１の検出値を引いた値が圧縮後吸気圧Ｐ１の変化率ΔＰ１として算出される。この変化率ΔＰ１の値は、圧縮後吸気圧Ｐ１が減少しているときには負の値となる。そして、その変化率ΔＰ１の絶対値は、演算周期を単位時間としたときの、圧縮後吸気圧Ｐ１の単位時間あたりの変化量を表すことになる。なお、これ以外にも、演算周期ごとの圧縮後吸気圧Ｐ１の取得値の移動平均線の傾きを求め、その傾きを変化率ΔＰ１とするなど、他の方法で同変化率ΔＰ１を求めることも可能である。
【００３１】
また、上記（ｂ）の処理でのマップ値Ｋｍａｐの算出は、予め記憶されたマップを参照して行われる。このマップには、変化率ΔＰ１とマップ値Ｋｍａｐとの関係が記憶されている。
【００３２】
図２に示すように、マップ値Ｋｍａｐは、変化率ΔＰ１の絶対値が規定値α１未満の範囲では、一定の値βとなっている。変化率ΔＰ１の絶対値をα１から大きくしていくと、マップ値Ｋｍａｐの値は、βから次第に減少するようになる。なお、変化率ΔＰ１の絶対値がα１から規定値α２（＞α１）に達するまでの区間では、変化率ΔＰ１の絶対値の増加に対してマップ値Ｋｍａｐが急激に減少しているが、α２を超える区間では、変化率ΔＰ１の絶対値の増加に対してのマップ値Ｋｍａｐの減少はより緩やかとなっている。
【００３３】
さらに上記（ｃ）における係数Ｋの値の設定では、下記条件（イ）〜（ハ）の１つ以上が満たされていない場合、マップ値Ｋｍａｐの値がそのまま係数Ｋの値に設定される。
（イ）圧縮後吸気圧Ｐ１の単位時間あたりの変化量（変化率ΔＰ１の絶対値）の減少に応じて係数Ｋの値が変化されていること。すなわち、断熱圧縮吸気温Ｔ１ｍに対する推定吸気温Ｔ１の追従性が高くなる側に、係数Ｋの値が変化している過程にあること。なお、本実施形態では、図２に示すように圧縮後吸気圧Ｐ１の単位時間あたりの変化量（変化率ΔＰ１の絶対値）の減少に対するマップ値Ｋｍａｐの変化は、増加側への変化となっている。そして、このときのマップ値Ｋｍａｐの増加に応じて係数Ｋの値も増加するようになっている。
【００３４】
（ロ）断熱圧縮吸気温Ｔ１ｍが低下中であること。
（ハ）マップ値Ｋｍａｐの値をそのまま係数Ｋの値に設定した場合に、前回の演算周期からの係数Ｋの値の変化量ΔＫがガード値ΔＫｍａｘを超える値となること。
【００３５】
これに対して、上記条件（イ）〜（ハ）のすべてが満たされている場合、上記（ｄ）のガード処理が実施され、係数Ｋの値は、前回の演算周期において設定された係数Ｋの値（以下、前回値Ｋｏｌｄと記載する）からの変化量ΔＫがガード値ΔＫｍａｘとなる値に設定される。こうしたガード処理では、前回の演算周期からの係数Ｋの変化量ΔＫが、すなわち演算周期を単位時間としたときの、係数Ｋの単位時間あたりの変化量ΔＫが、ガード値ΔＫｍａｘ以下に制限されることになる。
【００３６】
なお、図３に示すように、ガード値ΔＫｍａｘの値は、圧縮後吸気圧Ｐ１の単位時間あたりの変化量（変化率ΔＰ１の絶対値）に応じた可変値とされている。すなわち、ガード値ΔＫｍａｘは、変化率ΔＰ１の絶対値が規定値γ１以下のときには、最大設定値ε１に設定され、変化率ΔＰ１の絶対値が規定値γ２（＞γ１）以上のときには、最小設定値ε２（＞ε１）に設定される。そして、変化率ΔＰ１の絶対値がγ１からγ２までの範囲にあるときには、ガード値ΔＫｍａｘの値は、同絶対値に増加に応じてε１からε２へと一定の比率で減少するように設定される。
【００３７】
以上のように、推定吸気温演算部２３は、推定吸気温Ｔ１の演算に用いる関数Ｆにおいて、断熱圧縮吸気温Ｔ１ｍに対する推定吸気温Ｔ１の追従性を定める係数Ｋの値を、圧縮後吸気圧Ｐ１の単位時間あたりの変化量（変化率ΔＰ１の絶対値）に応じて可変設定している。すなわち、推定吸気温演算部２３は、断熱圧縮吸気温Ｔ１ｍに対して一時遅れの値を返す上記関数Ｆにあって、本来は定数である係数Ｋを変数としたものを、推定吸気温Ｔ１の演算に用いている。
【００３８】
図４に、こうした吸気温推定装置２０における推定吸気温Ｔ１の演算にかかる処理のフローチャートを示す。吸気温推定装置２０では、同ルーチンの処理が、規定の演算周期ごとに繰り返し実行される。なお、同フローチャートは、各演算周期において、断熱圧縮吸気温演算部２２、及び推定吸気温演算部２３が行う処理を、一連の処理として表したものである。具体的には、同フローチャートにおけるステップＳ１００〜Ｓ１０３の処理は、断熱圧縮吸気温演算部２２により行われ、ステップＳ２００〜Ｓ２２０の処理は、推定吸気温演算部２３により行われる。また、ステップＳ２１０〜Ｓ２１４の処理は、上記（ｃ）の係数Ｋの値の設定に係る処理となっている。
【００３９】
本ルーチンの処理が開始されると、まずステップＳ１００〜Ｓ１０２において、圧縮前吸気温Ｔ０、圧縮前吸気圧Ｐ０、及び圧縮後吸気圧Ｐ１の各検出値が読み込まれる。そして、ステップＳ１０３において、読み込まれた各検出値に基づき、上式（１）に従って断熱圧縮吸気温Ｔ１ｍの演算が行われる。
【００４０】
続いて、ステップＳ２００において、圧縮後吸気圧Ｐ１の変化率ΔＰ１が算出される。そして、次のステップＳ２０１において、同変化率ΔＰ１に基づきマップ値Ｋｍａｐが算出される。
【００４１】
また、次のステップＳ２１０では、同変化率ΔＰ１に基づきガード値ΔＫｍａｘが算出され、続くステップＳ２１１では、断熱圧縮吸気温Ｔ１ｍが低下中であるか否かが、すなわち上記条件（ロ）を満たしているか否かが判定される。そして、断熱圧縮吸気温Ｔ１ｍが低下中であれば（ＹＥＳ）、ステップＳ２１２に処理が進められ、低下中でなければ（ＮＯ）、ステップＳ２１４に処理が進められる。本実施形態では、前回の演算周期における断熱圧縮吸気温Ｔ１ｍの演算値よりも今回の演算周期における断熱圧縮吸気温Ｔ１ｍの演算値の方が小さい（温度が低い）ことをもって、断熱圧縮吸気温Ｔ１ｍが低下中であると判定している。これ以外にも、演算周期ごとの断熱圧縮吸気温Ｔ１ｍの演算値の移動平均線の傾きを求め、その傾きから断熱圧縮吸気温Ｔ１ｍが低下中であるか否かを判定するなど、他の方法で同判定を行うようにしてもよい。
【００４２】
ステップＳ２１２に処理が進められると、同ステップＳ２１２において、上記ステップＳ２０１で算出したマップ値Ｋｍａｐが、現状の係数Ｋの値と上記ステップＳ２１０で算出したガード値ΔＫｍａｘとの和を超えているか否かが判定される。そして、マップ値Ｋｍａｐが上記和を超えていれば（ＹＥＳ）、ステップＳ２１３に処理が進められ、超えていなければ（ＮＯ）、ステップＳ２１４に処理が進められる。
【００４３】
なお、本ステップＳ２１２では、上記条件（イ）及び（ハ）を共に満たしているか否かを判定している。本ルーチンにおいて本ステップＳ２１２に処理が進められた時点では、係数Ｋの値の操作は行われておらず、このときの係数Ｋの値は、前回の演算周期に本ルーチンが実行されたときに設定された値、すなわち上述の前回値Ｋｏｌｄとなっている。
【００４４】
一方、図２に示すように、圧縮後吸気圧Ｐ１の単位時間あたりの変化量（変化率ΔＰ１の絶対値）の減少に対してマップ値Ｋｍａｐは増加側にのみ変化する。また、マップ値Ｋｍａｐが増加しているときには、係数Ｋの値もマップ値Ｋｍａｐの値を超えない範囲で増加する。よって、圧縮後吸気圧Ｐ１の単位時間あたりの変化量（変化率ΔＰ１の絶対値）の減少に応じて係数Ｋの値が変化（増加）しているときには、すなわち上記条件（イ）を満しているときには、（い）Ｋｍａｐ＞Ｋｏｌｄの関係となる。
【００４５】
これに対して、マップ値Ｋｍａｐの値をそのまま係数Ｋの値に設定した場合に、前回の演算周期からの係数Ｋの値の変化量ΔＫがガード値ΔＫｍａｘを超える値となる場合、すなわち上記条件（ハ）を満たす場合には、（ろ）Ｋｍａｐ＞Ｋｏｌｄ＋ΔＫｍａｘの関係となる場合と、（は）Ｋｍａｐ＜Ｋｏｌｄ−ΔＫｍａｘの場合とがある。このうち、（ろ）の関係が成立する場合には（い）の関係が自ずと成立するのに対して、（は）の関係が成立する場合には（い）の関係は必然的に不成立となる。よって、（ろ）の関係の成立、すなわちステップＳ２１２での肯定判定をもって、上記条件（イ）（ハ）を共に満たしていると判定することができる。
【００４６】
以上の処理の結果、ステップＳ２１３に処理が進められた場合、同ステップＳ２１３において、現状の係数Ｋの値とガード値ΔＫｍａｘとの和が、係数Ｋの新たな値として設定された後、ステップＳ２２０に処理が進められる。すなわち、このステップＳ２１３において、上記（ｄ）のガード処理が実施される。
【００４７】
これに対して、ステップＳ２１４に処理が進められた場合、同ステップＳ２１４において、上記ステップＳ２０１で算出されたマップ値Ｋｍａｐの値が係数Ｋの値として設定された後、ステップＳ２２０に処理が進められる。
【００４８】
ステップＳ２２０に処理が進められると、同ステップＳ２２０において、上記ステップＳ１０３で演算された断熱圧縮吸気温Ｔ１ｍと、上記ステップＳ２１３又は上記ステップＳ２１４で設定された係数Ｋの値とに基づき、上式（２）に従って、推定吸気温Ｔ１の演算が行われる。そして、推定吸気温Ｔ１の演算後、今回の演算周期における本ルーチンの処理が終了される。
【００４９】
（作用）
上図１のような過給エンジン１０では、吸気通路１１に取り込んだ吸気をコンプレッサ１３により圧縮して、燃焼室１５に送り込むようにしている。このときのコンプレッサ１３での圧縮により、吸気の温度は上昇する。したがって、燃焼室１５に流入する吸気の温度は、コンプレッサ１３での圧縮により上昇した後の温度、すなわち上述の圧縮後吸気温となる。
【００５０】
過給エンジン１０の加速や減速により、圧縮後吸気圧Ｐ１が変化すると、上記圧縮後吸気温も変化する。ただし、圧縮後吸気温が変化すると、吸気と外部（コンプレッサ１３の筐体や吸気管など）との間に温度差が生じて、圧縮後の吸気と外部との間に熱移動が発生するようになる。よって、このときの圧縮後吸気温の変化は、完全な断熱変化とはならず、圧縮後吸気圧Ｐ１の変化に対する圧縮後吸気温の変化には遅れが生じるようになる。すなわち、圧縮後吸気圧Ｐ１の変化に対する圧縮後吸気温の変化が完全な断熱変化であるとして演算された断熱圧縮吸気温Ｔ１ｍに対しても、実際の圧縮後吸気温は遅れて変化することになる。このときの断熱圧縮吸気温Ｔ１ｍに対する圧縮後吸気温の遅れは、圧縮後吸気圧Ｐ１の変化が急激なときには小さくなり、圧縮後吸気圧Ｐ１の変化が緩やかなときには大きくなる。そのため、一次遅れ系として簡潔に表すことはできないものとなる。
【００５１】
ここで、上式（２）における係数Ｋを定数として推定吸気温Ｔ１の演算を行う場合を考える。上述したように、上式（２）は、係数Ｋを定数とすると、断熱圧縮吸気温Ｔ１ｍに対して、すなわち圧縮後吸気圧Ｐ１に対して１次遅れの値を推定吸気温Ｔ１として返す関数Ｆとなる。この場合、圧縮後吸気圧Ｐ１の変化が急激なときに適合するよう係数Ｋの値を設定すれば、同変化が緩やかなときには、推定吸気温Ｔ１の演算値は、実際の圧縮後吸気温よりも早く変化するようになる。とはいえ、同変化が緩やかなときに適合するように係数Ｋの値を設定すれば、同変化が急激なときには、推定吸気温Ｔ１の演算値は、実際の圧縮後吸気温よりも遅れて変化するようになる。
【００５２】
これに対して、本実施形態の吸気温推定装置２０では、推定吸気温演算部２３は、係数Ｋの値を、圧縮後吸気圧Ｐ１の単位時間あたりの変化量（変化率ΔＰの絶対値）に応じて可変設定した上で、推定吸気温Ｔ１の演算を行うようにしている。図２のように、圧縮後吸気圧Ｐ１の単位時間あたり変化量（変化率ΔＰ１の絶対値）が大きいときには、同変化量小さいときよりも、小さい値が係数Ｋに設定される。上述のように、係数Ｋの値を１から大きくしていくと、断熱圧縮吸気温Ｔ１ｍに対する推定吸気温Ｔ１の演算値の遅れ（時定数）は次第に大きくなる。よって、係数Ｋの値を小さくするほど、断熱圧縮吸気温Ｔ１ｍに対する推定吸気温Ｔ１の追従性が高くなる。すなわち、断熱圧縮吸気温Ｔ１ｍに対する推定吸気温Ｔ１の遅れは小さくなる。
【００５３】
こうした推定吸気温演算部２３により演算された推定吸気温Ｔ１は、圧縮後吸気圧Ｐ１の単位時間あたりの変化量が大きいとき、すなわち同圧縮後吸気圧Ｐ１の変化が急なときには、断熱圧縮吸気温Ｔ１ｍに対する遅れが小さくなる。また、圧縮後吸気圧Ｐ１の単位時間あたりの変化量が小さいとき、すなわち同圧縮後吸気圧Ｐ１の変化が緩やかなときには、断熱圧縮吸気温Ｔ１ｍに対する遅れが大きくなる。したがって、推定吸気温演算部２３により演算された推定吸気温Ｔ１は、上述した実際の圧縮後吸気温と同じように、圧縮後吸気圧Ｐ１の変化の緩急によって断熱圧縮吸気温Ｔ１ｍに対する遅れ方が変化するものとなる。しかも、推定吸気温演算部２３が演算に用いる関数Ｆそのものは、一次遅れ系を表す単純なものであり、その関数Ｆの係数Ｋの値を、圧縮後吸気圧Ｐ１の変化の緩急に合わせて変更しているだけである。したがって、本実施形態の吸気温推定装置２０によれば、コンプレッサ１３による圧縮後の吸気温を、簡易な演算で精度良く推定することができる。
【００５４】
ただし、圧縮後吸気圧Ｐ１の単位時間あたりの変化量（変化率ΔＰ１の絶対値）より算出された値（マップ値Ｋｍａｐ）をそのまま係数Ｋの値に設定すると、次の問題が生じることがある。
【００５５】
図５に、圧縮後吸気圧Ｐ１の低下時における変化率ΔＰ１、係数Ｋ、マップ値Ｋｍａｐ、断熱圧縮吸気温Ｔ１ｍ、実際の圧縮後吸気温（以下、実吸気温Ｔ１ｔと記載する）等の推移を模式的に示す。同図では、時刻ｔ１から時刻ｔ４にかけて圧縮後吸気圧Ｐ１が低下しており、そのうち、ｔ２〜ｔ３の期間には、その前後のｔ１〜ｔ２の期間、及びｔ３〜ｔ４の期間よりも、圧縮後吸気圧Ｐ１の低下が急激なものとなっている。すなわち、ｔ２〜ｔ３の期間における圧縮後吸気圧Ｐ１の単位時間あたりの変化量（変化率ΔＰ１の絶対値）は、ｔ１〜ｔ２の期間、及びｔ３〜ｔ４の期間よりも大きくなっている。
【００５６】
上述のように、このときの実吸気温Ｔ１ｔの遅れは、圧縮後吸気圧Ｐ１の急激な変化に対しては小さく、緩やかな変化に対しては大きくなる。ただし、そうした実吸気温Ｔ１ｔの遅れの変化が現れるのは、圧縮後吸気圧Ｐ１の変化率ΔＰ１の変化から若干の時間が経過してからとなる。
【００５７】
一方、マップ値Ｋｍａｐは、圧縮後吸気圧Ｐ１の単位時間あたりの変化量（変化率ΔＰ１の絶対値）が減少すると、値が即座に大きくされる。そのため、上述のガード処理を行わず、マップ値Ｋｍａｐの値をそのまま係数Ｋの値に設定して、推定吸気温Ｔ１を演算すれば、同図に一点鎖線で示すように、そうして演算された推定吸気温Ｔ１の低下は、実吸気温Ｔ１ｔよりも遅くなってしまう。
【００５８】
これに対して、本実施形態では、変化率ΔＰ１の絶対値の減少に応じて上記係数Ｋの値が増加側に変化するときに、同係数Ｋの値の単位時間あたりの変化量をガード値ΔＫｍａｘ以下に制限するガード処理を行っている。そのため、圧縮後吸気圧Ｐ１の急激な変化が収束した時刻ｔ３以降もしばらくは、係数Ｋに小さい値が設定された状態が、ひいては断熱圧縮吸気温Ｔ１ｍに対する推定吸気温Ｔ１の追従性が高い状態が続く。そのため、推定吸気温Ｔ１の演算結果が、実吸気温Ｔ１ｔの挙動により近似するようになる。
【００５９】
なお、上記のようなガード処理の結果、圧縮後吸気圧Ｐ１の変化が収束した後も、係数Ｋの値の変更が制限された状態が続けられると、推定吸気温Ｔ１の演算値が、実吸気温Ｔ１ｔよりも早く、圧縮後吸気圧Ｐ１の変化後の値に収束してしまうことがある。すなわち、実吸気温Ｔ１ｔが圧縮後吸気圧Ｐ１の変化後の値に収束する時刻ｔ５よりも早い時期に、推定吸気温Ｔ１の演算値が収束してしまうようになる。その点、本実施形態では、図３に示すように、圧縮後吸気圧Ｐ１の単位時間あたりの変化量（変化率ΔＰ１の絶対値）が大きいときには、同変化量が小さいときよりも、ガード値ΔＫｍａｘを大きい値に設定している。こうした場合、圧縮後吸気圧Ｐ１の変化が収束に向かい、その単位時間あたりの変化量が減少するに従って、ガード処理による係数Ｋの変化量の制限が緩和される。そのため、圧縮後吸気圧Ｐ１の変化の収束に応じてより速やかに、ガード処理による係数Ｋの変化量の制限を解除することができる。
【００６０】
なお、図６に示すように、過給エンジン１０の加速時における圧縮後吸気圧Ｐ１の上昇には、同過給エンジン１０の減速時における圧縮後吸気圧Ｐ１の低下よりも時間がかかる。そのため、圧縮後吸気温が上昇する圧縮後吸気圧Ｐ１の上昇時には、上記ガード処理でのガード値ΔＫｍａｘによる係数Ｋの値の単位時間あたりの変化量の制限を行わずとも、推定吸気温Ｔ１の演算結果を実吸気温Ｔ１ｔに近似させることができる。そこで、本実施形態では、推定吸気温Ｔ１の演算にかかる負荷をより少なくするため、ガード処理における係数Ｋの単位時間あたりの変化量の制限を、断熱圧縮吸気温Ｔ１ｍの上昇中には行わず、同断熱圧縮吸気温Ｔ１ｍの低下中にのみ行うようにしている。
【００６１】
また、過給エンジン１０の加速時における圧縮後吸気圧Ｐ１の上昇や、減速時における同圧縮後吸気圧Ｐ１の低下は、単調な上昇、低下とはならず、上昇の過程において一時的な圧縮後吸気圧Ｐ１の低下が生じたり、低下の過程において一時的な圧縮後吸気圧Ｐ１の上昇が生じたりすることがある。こうした一時的な圧縮後吸気圧Ｐ１の変動に対して、上記係数Ｋの値を変更すると、その一時的な変動に対して推定吸気温Ｔ１の演算値が鋭敏に反応して、同演算値に不安定な変動が現れることがある。その点、本実施形態では、図２に示すように、圧縮後吸気圧Ｐ１の単位時間あたりの変化量（変化率ΔＰ１の絶対値）がα１以下の範囲では、マップ値Ｋｍａｐは一定の値βとされており、上記変化量に対して係数Ｋが変化しない不感帯が設定されている。そのため、推定吸気温Ｔ１の演算値には、上記のような圧縮後吸気圧Ｐ１の一時的な変動の影響が現れにくくなる。
【００６２】
図７は、過給エンジン１０の加速による圧縮後吸気圧Ｐ１の上昇時における、断熱圧縮吸気温Ｔ１ｍ、推定吸気温Ｔ１の演算値、及び実吸気温Ｔ１ｔの推移を示している。また、図８は、過給エンジン１０の減速による圧縮後吸気圧Ｐ１の低下時における、断熱圧縮吸気温Ｔ１ｍ、推定吸気温Ｔ１の演算値、及び実吸気温Ｔ１ｔの推移を示している。さらに両図には、比較例として、圧縮後吸気圧Ｐ１が急激に変化しているときに適合する値に係数Ｋを固定したときの推定吸気温の演算値Ｔ１ａの推移と、圧縮後吸気圧Ｐ１が緩やかに変化しているときに適合する値に係数Ｋを固定したときの推定吸気温の演算値Ｔ１ｂの推移とが併せ示されている。
【００６３】
上述のように、断熱圧縮吸気温Ｔ１ｍに対する実吸気温Ｔ１ｔの遅れは、同断熱圧縮吸気温Ｔ１ｍの、すなわち圧縮後吸気圧Ｐ１の変化の緩急により変化するため、係数Ｋを一定の値に固定すると、そうした遅れの変化に対応できない場面ができてしまう。したがって、両図に示されるように、係数Ｋの値を固定したときの推定吸気温の演算値Ｔ１ａ，Ｔ１ｂの軌跡はいずれも、実吸気温Ｔ１ｔの軌跡から大きく乖離したものとなっている。これに対して、圧縮後吸気圧Ｐ１の単位時間あたりの変化量に応じて係数Ｋの値を可変設定するようにした本実施形態の吸気温推定装置２０での推定吸気温Ｔ１の演算値の軌跡は、実吸気温Ｔ１ｔの軌跡に近似したものとなっている。
【００６４】
なお、上記実施形態は、以下のように変更して実施することもできる。
・上記実施形態では、圧縮後吸気圧Ｐ１の単位時間あたりの変化量（変化率ΔＰ１の絶対値）がα１未満のときには、係数Ｋの値をβに固定するとともに、同変化量がα１を超えて大きくなるに従って、係数Ｋの値をβから次第に小さくするようにしていた。すなわち、上記実施形態では、上記変化量が一定の大きさに達するまでは、同変化量に対して係数Ｋの値が変化しない不感帯を設定するようにしていた。そして、これにより、圧縮後吸気圧Ｐ１の変化過程での一時的な変動の影響が、推定吸気温Ｔ１の演算値に現れにくくなるようにしていた。こうした不感帯を設定せず、上記変化量のすべての区間において、同変化量に応じた係数Ｋの値の可変設定を行うようにしてもよい。
【００６５】
・上記のように不感帯を設定しないことを選択した場合において、変化過程の圧縮後吸気圧Ｐ１の一時的な変動の影響が推定吸気温Ｔ１の演算値に与える影響を無視できない場合には、次のようにするとよい。すなわち、圧縮後吸気圧Ｐ１の検出値や、その単位時間あたりの変化量の算出値にフィルタをかけて、上記一時的な変動を予め均しておくようにする。このようにすれば、上記のような不感帯が設定されていなくても、推定吸気温Ｔ１の演算結果に、圧縮後吸気圧Ｐ１の一時的な変動の影響が現れにくくなる。
【００６６】
・上記実施形態では、ガード処理を、断熱圧縮吸気温Ｔ１ｍの上昇中には行わず、同断熱圧縮吸気温Ｔ１ｍの低下中にのみ行うようにしていたが、断熱圧縮吸気温Ｔ１ｍの上昇中にもガード処理を実施するようにしてもよい。例えば、図４のフローチャートにおけるステップＳ２１１の処理を割愛し、ステップＳ２１０の処理の後、ステップＳ２１２に処理が進められるように、同フローチャートの処理手順を変更すれば、断熱圧縮吸気温Ｔ１ｍの上昇中にもガード処理が実施されるようになる。
【００６７】
・上記実施形態では、圧縮後吸気圧Ｐ１の単位時間あたりの変化量（変化率ΔＰ１の絶対値）が大きいほど、ガード処理のガード値ΔＫｍａｘを大きい値に設定していた。そして、これにより、圧縮後吸気圧Ｐ１の変化が収束に向うにつれて、ガード処理による係数Ｋの変化量の制限を緩和するようにしていた。圧縮後吸気圧Ｐ１の単位時間あたりの変化量の減少に応じてガード値ΔＫｍａｘの値を小さくする代わりに、或いはそれに加え、圧縮後吸気圧Ｐ１の単位時間あたりの変化量が一定の値まで低下した時点で、ガード処理による係数Ｋの単位時間あたりの変化量の制限を解除するようにしてもよい。このようにすれば、圧縮後吸気圧Ｐ１の変化がある程度まで収束した時点で、ガード処理による係数Ｋ１の変化量の制限が確実に解除されるようになる。
【００６８】
・上記実施形態では、圧縮後吸気圧Ｐ１の単位時間あたりの変化量（変化率ΔＰ１の絶対値）に応じてガード値ΔＫｍａｘの値を可変設定するようにしていたが、同ガード値ΔＫｍａｘを定数としてもよい。そうした場合にも、圧縮後吸気圧Ｐ１の単位時間あたりの変化量の減少に応じて係数Ｋの値が即座に減少しないようになる。そのため、圧縮後吸気圧Ｐ１の変化の収束後も変化を続ける実吸気温Ｔ１ｔの挙動に近似した、推定吸気温Ｔ１の演算結果が得られる。
【００６９】
・上記実施形態では、圧縮後吸気圧Ｐ１の単位時間あたりの変化量の減少に応じて係数Ｋの値が変化するときの、同係数Ｋの値の単位時間あたりの変化量をガード値ΔＫｍａｘ以下に制限するガード処理を行っていたが、そうしたガード処理を行わないようにしてもよい。ガード処理がその効果を発揮するのは、圧縮後吸気圧Ｐ１の変化がある程度よりも急となる場合である。よって、同変化がそこまで急とならない場合には、ガード処理を行わずとも、コンプレッサ１３による圧縮後の吸気温を精度良く推定できる。
【００７０】
・上記実施形態では、上式（２）を、断熱圧縮吸気温Ｔ１ｍに対して一次遅れの値を返す関数として用いて推定吸気温Ｔ１を演算していた。そうした関数として、上記以外の関数を用いるようにしてもよい。例えば、断熱圧縮吸気温Ｔ１ｍの加重移動平均値を返す関数なども、断熱圧縮吸気温Ｔ１ｍに対して一次遅れの値を返す関数となる。そうした場合、上記実施形態における係数Ｋの代わりに、その関数において断熱圧縮吸気温Ｔ１ｍに対する推定吸気温Ｔ１の追従性を定める係数の値を、圧縮後吸気圧Ｐ１の単位時間あたりの変化量に応じて可変設定することになる。例えば、上記加重移動平均値を返す関数の場合、複数の断熱圧縮吸気温Ｔ１ｍの算出値の各々に付される重み付け係数が、上記実施形態における係数Ｋの代わりとなる。このように使用する関数によっては、断熱圧縮吸気温Ｔ１ｍに対する推定吸気温Ｔ１の追従性を定める係数が複数存在することがある。また、使用する関数によっては、推定吸気温Ｔ１の追従性を高くする際に、係数の値を大きくするようになることもある。いずれにせよ、圧縮後吸気圧Ｐ１の単位時間あたりの変化量が大きいときには、同変化量が小さいときよりも、断熱圧縮吸気温Ｔ１ｍに対する推定吸気温Ｔ１の追従性が高くなるように、そうした係数の値を操作すれば、実吸気温Ｔ１ｔにより近似した推定吸気温Ｔ１の演算結果を得ることができる。
【符号の説明】
【００７１】
１０…過給エンジン、１１…吸気通路、１３…コンプレッサ、１５…燃焼室、１７…温度センサ、１８…第１圧力センサ１８、１９…第２圧力センサ、２０…吸気温推定装置、２２…断熱圧縮吸気温演算部、２３…推定吸気温演算部。
【図１】
0006376152.tif 000006
【図２】
0006376152.tif 000007
【図３】
0006376152.tif 000008
【図４】
0006376152.tif 000009
【図５】
0006376152.tif 000010
【図６】
0006376152.tif 000011
【図７】
0006376152.tif 000012
【図８】
0006376152.tif 000013