JP3277915B2 - 内燃機関の吸気管圧力予測方法および装置 - Google Patents

Description

【発明の詳細な説明】 【０００１】 【発明の属する技術分野】本発明は、内燃機関の吸気管
圧力予測方法および装置に係り、より詳しくは、スロッ
トル開度と機関回転速度とに基づいて吸気管圧力を検出
する内燃機関の吸気管圧力予測方法および装置に関す
る。 【０００２】 【従来の技術】従来より、内燃機関の吸気管予測方法
は、内燃機関の燃料噴射量制御方法などに用いられてい
る。従来の内燃機関の燃料噴射量制御方法は、吸気管圧
力と機関回転速度とを検出し、この検出された吸気管圧
力と機関回転速度とに基づいて基本燃料噴射時間を演算
し、この基本燃料噴射時間を吸気温や機関冷却水温等に
応じて補正して燃料噴射時間を求め、この燃料噴射時間
に相当する時間燃料噴射弁を開弁することにより燃料噴
射量を制御するものである。この燃料噴射量制御方法に
おいては、ダイヤフラム式の圧力センサを吸気管に取付
け、機関脈動成分を除去するために時定数が３〜５msec
のフイルタを介して圧力センサ出力を処理することによ
り吸気管圧力を検出し、この検出された吸気管圧力と回
転速度センサで検出された機関回転速度とに基づいて基
本燃料噴射時間を演算するようにしている。 【０００３】しかしながら、圧力センサのダイヤフラム
による応答遅れおよびフイルタの時定数による応答遅れ
が存在するため、加減速時等の過渡運転時には、実際の
吸気管圧力の変動に対して検出された吸気管圧力の変動
に時間遅れが生ずる。このため、加速時にはスロツトル
弁が急開されて実際の吸気管圧力が急激に上昇するのに
対して検出された吸気管圧力に時間遅れが生じ、実際の
吸気管圧力より小さい値の吸気管圧力によって基本燃料
噴射時間が演算されることになるため、空燃比がオーバ
リーンになり加速応答性が悪化すると共に排気エミツシ
ヨンが悪化する。逆に、減速時にはスロツトル弁が急閉
されることから吸気管圧力が急激に低下するため実際の
吸気管圧力より大きな値の吸気管圧力によって基本燃料
噴射時間が演算されることになり、空燃比がオーバリツ
チになってドライバビリテイが悪化すると共に排気エミ
ツシヨンが悪化する。この空燃比のオーバリツチおよび
オーバリーンを防止するために、加速増量や減速減量等
の各種の増減量補正を行なっているが、過渡時には検出
された吸気管圧力に時間遅れが存在するため全運転領域
で完全に目標空燃比に制御することが不可能であった。 【０００４】一方、実際の値に対して時間遅れのない物
理量としてスロツトル開度を用い、スロツトル開度と機
関回転速度とに基づいて基本燃料噴射時間を演算する方
法（特開昭５９−２８０３１号公報）やスロツトル開度
と機関回転速度とに応じた吸気管圧力を記憶しておいて
圧力センサより得られる信号に応じて排ガス再循環時の
排ガスの分圧を考慮してこの吸気管圧力を補正して燃料
噴射量を制御することが行なわれている（特開昭５９−
３９９４８号公報）。 【０００５】 【発明が解決しようとする課題】しかしながら、通常ス
ロツトル弁は圧力センサ取付位置および機関燃焼室から
離れた上流側の位置に配置されており、スロツトル弁を
通過した空気が圧力センサ取付位置や機関燃焼室へ到達
するまでに時間遅れが生じ、また、スロツトル弁と吸気
弁との間の容積のためスロツトル開度は実吸入空気量の
変化に対して位相が進むことになる。このため、スロツ
トル開度と機関回転速度とで定められた吸気管圧力Ｐ
（ＴＡ、ＮＥ）は図３に示すように実際の吸気管圧力Ｐ
より位相が進んだ値となる。なお、ＰＭは圧力センサか
ら得られる吸気管圧力である。また、図４に示すよう
に、スロツトル開度と機関回転速度とで定められた基本
燃料噴射量ＴＰ（ＴＡ、ＮＥ）はスロツトル開度の変化
が実吸入空気量の変化に対して位相が進んでいるため要
求燃料噴射量よりも多くなる。このためスロツトル開度
と機関回転速度とに基づいて燃料噴射量を制御すると、
加速時には燃料噴射量が要求値より多くなって空燃比が
オーバリツチになり、減速時には燃料噴射量が要求値よ
り少なくなって空燃比がオーバリーンになる。また、加
速増量補正を行なった場合においても増量値は図４の斜
線で示すようになり、上記の位相進みを補正することは
できない。 【０００６】本発明は上記問題点を解決すべく成された
もので、実際の吸気管圧力の変化に対して応答遅れのな
いスロツトル開度を用いて位相進みおよび遅れのない吸
気管圧力、すなわち実際の吸気管圧力を予測することが
できるようにした内燃機関の吸気管圧力推定方法および
装置を提供することを目的とする。 【０００７】 【課題を解決するための手段】上記目的を達成するため
に請求項１の発明は、以下の式にしたがって演算された
加重平均値ＰＭＳＭ _i を次回の演算において以下の式の
過去に演算された加重平均値ＰＭＳＭ _i-1 として用いて
所定回数繰り返し演算して加重平均値ＰＭＳＭ _i を求
め、繰り返し演算後の加重平均値ＰＭＳＭ _i を吸気管圧
力予測時点での吸気管圧力として予測するものである。 ＰＭＳＭ _i ＝（（ｎ−１）・ＰＭＳＭ _i-1 ＋ＰＭＴＡ）/ｎ ただし、ＰＭＴＡは、スロットル開度と機関回転速度と
に基づいて演算されると共に、繰り返し演算中は一定と
される定常状態での吸気管圧力、ｎは、過渡時の吸気管
圧力の変化に関する時定数と前記演算周期とを用いて演
算される重みに関する係数である。 【０００８】請求項２記載の発明は、請求項１記載の発
明において、前記定常状態での吸気管圧力ＰＭＴＡを、
スロットル開度の微分値及び機関回転速度の微分値の少
なくとも一方を用いて演算される吸気管圧力予測時点に
おける定常状態での吸気管圧力としたものである。 【０００９】請求項３記載の発明は、請求項１または請
求項２記載の発明において、前記演算周期を一定とし、
スロットル開度と機関回転速度とに基いて、又は、定常
状態での吸気管圧力と機関回転速度とに基いて、前記重
みに関する係数ｎを演算するものである。 【００１０】請求項４記載の発明は、請求項１乃至３請
求項の何れか１項に記載の発明において、前記所定回数
を、現在時刻から吸気管圧力予測時点までの時間を演算
周期で除算して得られる値としたものである。 【００１１】 【００１２】請求項５記載の発明は、スロットル開度を
検出するスロットル開度検出手段と、機関回転速度を検
出する回転速度検出手段と、前記検出されたスロットル
開度と前記検出された機関回転速度とに基づいて所定の
演算周期で定常状態での吸気管圧力ＰＭＴＡを演算する
吸気管圧力演算手段と、前記演算周期を一定とし、演算
された定常状態での吸気管圧力と検出された機関回転速
度とに基いて、または検出されたスロットル開度と検出
された機関回転速度とに基づいて、過渡時の吸気管圧力
の変化に関する時定数と前記演算周期とで定まる重みに
関する係数ｎを演算する重み係数演算手段と、繰り返し
演算中は定常状態での吸気管圧力ＰＭＴＡを一定とする
と共に、以下の式にしたがって演算された加重平均値Ｐ
ＭＳＭ _i を次回の演算において以下の式の過去に演算さ
れた加重平均値ＰＭＳＭ _i-1 として用て所定回数繰り返
し演算して加重平均値ＰＭＳＭ _i を求め、繰り返し演算
後の加重平均値ＰＭＳＭ _i を吸気管圧力予測時点での吸
気管圧力として予測する吸気管圧力予測手段と、を備え
ている。 ＰＭＳＭ _i ＝（（ｎ−１）・ＰＭＳＭ _i-1 ＋ＰＭＴＡ）/ｎ 【００１３】請求項６記載の発明は、請求項５記載の発
明において、前記定常状態での吸気管圧力ＰＭＴＡを、
スロットル開度の微分値及び機関回転速度の微分値の少
なくとも一方を用いて演算される吸気管圧力予測時点に
おける定常状態での吸気管圧力としたものである。 【００１４】請求項７記載の発明は、請求項５記載の発
明において、前記所定回数を、現在時刻から吸気管圧力
予測時点までの時間を演算周期で除算して得られる値と
したものである。 【００１５】 【００１６】 【００１７】最初に、参考として、実際の吸気管圧力を
求める原理について説明する。図１に示すように、スロ
ツトル弁ＴｈからサージタンクＳを介して機関Ｅの吸気
弁までの吸気系を考え、吸気系内の空気の圧力（吸気管
絶対圧力）をＰ［mmHgabs.］、吸気系の容積をＶ
［ｌ］、吸気系内に存在する空気の重量をＱ［ｇ］、吸
気系内の空気の絶対温度をＴ［°Ｋ］、大気圧をＰｃ
［mmHgabs.］とすると共に、吸気系から機関Ｅの燃焼室
に吸入される単位時間当りの空気重量をΔＱ₁［ｇ/se
c］、スロツトル弁Ｔｈを通過して吸気系内に吸入され
る単位時間当りの空気重量をΔＱ₂［ｇ/sec］とし、微
小時間Δｔ内に吸気系の空気の重量が（ΔＱ₂−ΔＱ ₁)
・Δｔ変化し、このとき吸気系内の空気の圧力がΔＰ変
化したものとして、吸気系内の空気にボイル・シヤルル
の法則を適用すると以下の（1）式に示すようになる。 【００１８】 【数１】 【００１９】ただし、Ｒは気体定数である。 【００２０】一方、ＰＶ＝Ｑ・Ｒ・Ｔであるから上記
（1）式を変形すると、以下の（2）式が得られる。 【００２１】 【数２】 【００２２】ここで、流量係数をψ、スロツトル弁の開
口面積（スロツトル開度）をＡとするとスロツトル弁を
通過する単位時間当りの空気重量ΔＱ₂は以下の（3）式
で表わされ、行程容積をＶＳ、機関回転速度をＮＥ［ｒ
ｐｍ］、吸入効率をηとすると機関の燃焼室に吸入され
る単位時間当りの空気重量ΔＱ₁は以下の（4）式で表わ
される。 【００２３】 【数３】 【００２４】上記（3）、（4）式を（2）式に代入する
と次の（5）式が得られる。 【００２５】 【数４】 【００２６】ここで、Δｔ→０の極限をとると、 【００２７】 【数５】 【００２８】となる。 【００２９】今、圧力Ｐ₀（≠Ｐ_C）近傍での応答を考え
て圧力がＰ₀からＰ₀＋Ｐに変化したものとして、上記
（6）式のＰに代えてＰ₀＋Ｐ（ただし、Ｐは微小値）を
代入すると、以下の（7）式が得られる。 【００３０】 【数６】 【００３１】ここで、 【００３２】 【数７】【００３３】であるから、上記（7）式は以下の（9）式
のようになる。 【００３４】 【数８】 【００３５】ここで、 【００３６】 【数９】 【００３７】とすると、上記（9）式は次のようにな
る。 【００３８】 【数１０】 【００３９】上記（12）式を次の（13）式のように変形
して両辺を積分し、積分定数をＣとすると以下の（14）
式が得られる。 【００４０】 【数１１】 【００４１】ここでｔ＝０のときＰの初期値はＰ₀であ
るから上記（14）式より積分定数Ｃは次のようになる。 【００４２】 【数１２】 【００４３】上記（14）式と（15）式からＰを求めると
次のようになる。 【００４４】 【数１３】 【００４５】ただし、ｅは自然対数の底である。 【００４６】従って、スロツトル弁の開口面積Ａすなわ
ちスロツトル開度ＴＡ、機関回転速度ＮＥおよびスロツ
トル開度変化時点からの経過時間ｔを測定して上記（1
6）式に代入すれば、実際の吸気管圧力Ｐを求めること
ができる。そして、このようにして求めた実際の吸気管
圧力Ｐと機関回転速度ＮＥとに基づいて、例えば以下の
式に示す演算を行なって基本燃料噴射時間ＴＰを求め、
この基本燃料噴射時間ＴＰを吸気温や機関冷却水温等に
応じて補正して燃料噴射時間を求め、この燃料噴射時間
に相当する時間燃料噴射弁を開弁することにより機関が
要求する量の燃料を噴射することができる。 【００４７】 【数１４】 【００４８】ただし、Ｋは定数である。 【００４９】ところで、上記（16）式の吸気管圧力Ｐを
グラフで表わすと図２に示すようになり、ｔ＝０でＰ＝
Ｐ₀、ｔ→∞の極限（定常状態）ではＰ＝ｂ／ａ（定常
状態での吸気管圧力ＰＭＴＡ）となる１次遅れ要素の出
力である。従って、スロツトル開度ＴＡと機関回転速度
ＮＥとに基づいて定常状態での吸気管圧力ＰＭＴＡを演
算し、定常状態での吸気管圧力ＰＭＴＡを以下の（17）
式の伝達関数Ｇ(s) で表わされる１次遅れ要素で処理す
ることにより実際の吸気管圧力を演算するようにしても
よい。 【００５０】 【数１５】 【００５１】ただし、ｓはラプラス変換の演算子、Ｔは
時定数である。 【００５２】すなわち、所定周期毎にスロツトル開度と
機関回転速度とに基づいて定常状態での吸気管圧力を演
算し、演算された定常状態での吸気管圧力を１次遅れ要
素で処理することにより前記経過時間を変数とする吸気
管圧力を演算することもできる。 【００５３】このように、実際の吸気管圧力を予測して
この吸気管圧力と機関回転速度とに基づいて燃料噴射量
を制御すれば、実際の吸入空気量に応じた量の燃料を噴
射することができ、これによって空燃比を目標空燃比に
制御して過渡時の空燃比のオーバリツチ、オーバリーン
を防止することができる。 【００５４】次に、図６を参照して、他の参考例を説明
する。図６に示すように、スロツトル開度検出手段で検
出されたスロツトル開度ＴＡと回転速度検出手段で検出
された機関回転速度とに基づいて吸気管圧力演算手段Ａ
によって定常状態での吸気管圧力ＰＭＴＡが演算され
る。吸気管圧力演算手段Ａによって演算された定常状態
での吸気管圧力ＰＭＴＡは、補正手段Ｂによって過渡時
の吸気管圧力の応答遅れ分の補正が行なわれる。この補
正手段としては１次遅れ要素を用いることができる。補
正手段Ｂによって補正された吸気管圧力は基本燃料噴射
時間演算手段Ｃに入力され、基本燃料噴射時間演算手段
に入力されている機関回転速度ＮＥとに基づいて基本燃
料噴射時間ＴＰが演算される。そして、基本燃料噴射時
間ＴＰに基づいて燃料噴射量制御手段によって燃料噴射
量が制御される。 【００５５】このように、圧力センサやフイルタを用い
ていないので、構造簡単でかつ精度よく実際の吸気管圧
力を予測して機関が要求する量の燃料を噴射することが
できる。 【００５６】 【００５７】 【００５８】 【００５９】 【００６０】 【００６１】 【００６２】 【００６３】 【００６４】 【００６５】 【００６６】次に、本発明の原理を説明する。１次遅れ
要素をブロツク図で表わすと図５に示すようになり、入
力をｘ(t) とし、出力をｙ(t) とし、時定数をＴとする
と、図５の入出力の関係は以下の式で表わされる。 【００６７】 【数１６】 【００６８】ここで、ｔ₂を現在の演算タイミング、ｔ₁
を過去の演算タイミングとすると次の(21)式が得られ
る。 【００６９】 【数１７】 【００７０】上記(21)において、ｘ(t₂)を定常状態での
吸気管圧力ＰＭＴＡ、ｙ(t₂)を現在の実際の吸気管圧力
ＰＭＳＭ_i、ｙ(t₁)を過去の実際の吸気管圧力ＰＭＳＭ
_i-1、ｔ₂−ｔ₁（＝Δｔ）を演算周期とすれば、 【００７１】 【数１８】 【００７２】となり、Ｔ／Δｔ＝ｎとすると、以下の(2
3)式が得られる。 【００７３】 【数１９】 【００７４】すなわち、上記(23)式は、過去の実際の吸
気管圧力ＰＭＳＭ_i-1の重みをｎ−１とし、定常状態で
の吸気管圧力ＰＭＴＡの重みを１とした加重平均を求め
ることにより、現在の実際の吸気管圧力ＰＭＳＭ_iを演
算することができることを示している。また、重みに関
する係数ｎは時定数Ｔと演算周期Δｔとの比で求められ
る。 【００７５】従って、スロツトル開度と機関回転速度と
に基づいて所定周期Δｔで定常状態での吸気管圧力ＰＭ
ＴＡを演算し、過渡時の吸気管圧力の変化に関する時定
数Ｔと所定周期Δｔとで重みに関する係数ｎを演算し、
過去に演算された加重平均値ＰＭＳＭ_i-1の重みを重く
して過去に演算された加重平均値ＰＭＳＭ_i-1と定常状
態での吸気管圧力ＰＭＴＡと重みに関する係数ｎとで上
記(23)式に従って加重平均値ＰＭＳＭ_iを演算すれば、
現在の実際の吸気管圧力が求められることになる。 【００７６】そして、上記のようにして演算された加重
平均値（現在の実際の吸気管圧力）と機関回転速度とに
基づいて基本燃料噴射時間を演算し、演算された基本燃
料噴射時間に基づいて燃料噴射量を制御するようにして
もよい。 【００７７】なお、上記（10）、（16）式から理解され
るように、時定数Ｔ＝１／ａは機関回転速度ＮＥが大き
くなる程小さくなり、スロツトル開度ＴＡが大きくなる
程小さくなる。このように、時定数はスロツトル開度Ｔ
Ａと機関回転速度ＮＥを変数とする関数で表わされる。
従って演算周期Δｔを一定とすれば、重みに関する係数
ｎはスロツトル開度ＴＡと機関回転速度ＮＥとを変数と
する関数で定めることができる。なお、スロツトル開度
ＴＡと機関回転速度ＮＥとで定常状態での吸気管圧力Ｐ
ＭＴＡが一義的に定まるから、スロツトル開度ＴＡと機
関回転速度ＮＥとに代えて定常状態での吸気管圧力ＰＭ
ＴＡと機関回転速度ＮＥとに応じて重みに関する係数ｎ
を定めるようにしてもよい。 【００７８】ところで、機関燃焼室に供給される空気量
が確定するのは、吸気終了時点すなわち吸気弁閉弁時で
ある。しかしながら、燃料噴射時間を演算するために所
定時間必要であると共に、燃料噴射弁から噴射された燃
料が燃焼室に到達するまでに所定の飛行時間が必要であ
り、燃焼室に供給される空気量が確定したときに燃料噴
射量を演算すると時間遅れが生じるため、従来では、燃
焼室に供給される空気量が確定する前の吸気管圧力を用
いて基本燃料噴射時間を演算している。このため、実際
に燃焼室内に吸入された空気量に適合した量の燃料が噴
射されなくなり、加速時には吸入空気量が確定する吸気
管圧力より小さい値の吸気管圧力によって燃料噴射量が
制御されるため、空燃比がリーンとなり、減速時には吸
入空気量が確定する吸気管圧力より大きい値の吸気管圧
力によって燃料噴射量が制御されるため、空燃比がリツ
チとなる。 【００７９】一方、上記(23)式においてスロツトル開度
ＴＡと機関回転速度ＮＥとが変化しないものと仮定する
と、加重平均値演算時から吸入空気量が確定するまでの
間、すなわち加重平均値演算時から所定時間先までの間
定常状態での吸気管圧力ＰＭＴＡは一定である。従っ
て、上記(23)式の加重平均値を繰り返し演算することに
よって吸入空気量確定時の実際の吸気管圧力を予測する
ことができる。このため請求項１の発明では、上記(23)
式にしたがって演算された加重平均値ＰＭＳＭ _i を次回
の演算において上記(23)式の過去に演算された加重平均
値ＰＭＳＭ _i-1 として用いて所定回数繰り返し演算して
加重平均値ＰＭＳＭ _i を求め、繰り返し演算後の加重平
均値ＰＭＳＭ _i を吸気管圧力予測時点での吸気管圧力と
している。また、請求項５の発明は、スロットル開度を
検出するスロットル開度検出手段と、機関回転速度を検
出する回転速度検出手段と、前記検出されたスロットル
開度と前記検出された機関回転速度とに基づいて所定の
演算周期で定常状態での吸気管圧力ＰＭＴＡを演算する
吸気管圧力演算手段と、前記演算周期を一定とし、演算
された定常状態での吸気管圧力と検出された機関回転速
度とに基いて、または検出されたスロットル開度と検出
された機関回転速度とに基づいて、過渡時の吸気管圧力
の変化に関する時定数と前記演算周期とで定まる重みに
関する係数ｎを演算する重み係数演算手段と、繰り返し
演算中は定常状態での吸気管圧力ＰＭＴＡを一定とする
と共に、上記(23)式にしたがって演算された加重平均値
ＰＭＳＭ _i を次回の演算において上記(23)式の過去に演
算された加重平均値ＰＭＳＭ _i-1 として用て所定回数繰
り返し演算して加重平均値ＰＭＳＭ _i を求め、繰り返し
演算後の加重平均値ＰＭＳＭ _i を吸気管圧力予測時点で
の吸気管圧力として予測する吸気管圧力予測手段と、を
備えている。 請求項１及び請求項５の発明では、定常状
態での吸気管圧力を演算した時点から機関に吸入される
空気量が確定するまでの時間を演算周期Δｔで除算する
ことにより演算回数を求め、この演算回数だけ上記(23)
式の加重平均の演算を繰り返すことにより、機関に吸入
される空気量が確定する時点での加重平均値すなわち機
関に吸入される空気量が確定する時点での実際の吸気管
圧力を予測して燃料噴射量を制御するのが好ましい。 【００８０】なお、上記では燃料噴射時間演算時から機
関に吸入される空気量が確定するまでの間スロツトル開
度と機関回転速度とが変化しないものと仮定したが、ス
ロツトル開度や機関回転速度が変化する場合には、燃料
噴射時間演算時でのスロツトル開度の微分値および／ま
たは機関回転速度の微分値を用いて次の燃料噴射時間演
算時点でのスロツトル開度および／または機関回転速度
を予測して、吸入空気量が確定するときの定常状態での
吸気管圧力を予測し、上記のように加重平均値の演算を
繰り返して実際の吸気管圧力を予測すれば、スロツトル
開度や機関回転速度変動時の実際の吸気管圧力の予測値
の精度が更に向上する（上記請求項２記載の発明参
照）。 【００８１】また、燃料噴射弁から噴射された燃料は、
インテークマニホールド内壁面等の機関壁面に付着して
噴射された燃料の全てが燃焼室に供給されないので、こ
の燃料付着分を補正して燃料噴射量を制御するのが好ま
しい。この燃料付着量は、吸気管圧力の大きさに依存
し、吸気管圧力が小さいと燃料の蒸発量が多くなるため
燃料付着量は少なくなり、吸気管圧力が大きいと燃料の
蒸発量が少なくなるため燃料付着量は多くなる。このた
め本態様では、加重平均によって演算された実際の吸気
管圧力から機関壁面への燃料付着量の変化量を予測し、
この変化量に相当する量の燃料噴射量を補正して機関に
吸入される実際の吸入空気量に対応した量の燃料を機関
に供給するのが好ましい。なお、燃料の壁面への付着量
は機関温度や機関回転速度によっても変化する（機関温
度が高いと燃料の蒸発量が多くなるため燃料付着量は少
なくなり、機関回転速度が速くなると空気流速が速くな
って蒸発量が多くなるため燃料付着量は少なくなる）た
め、機関温度や機関回転速度の関数として燃料付着量の
変化量を定めるようにしても良く、また壁面への燃料付
着量は瞬時に安定しないため燃料噴射量の補正量を時間
減衰させて今回噴射したときの燃料付着量を次回以降の
噴射に反映させるようにしても良い。 【００８２】以上説明したように本発明においては、所
定周期で加重平均値を演算することにより実際の吸気管
圧力を予測しているため、スロツトル開度変化時点から
の経過時間を計測することなく実際の吸気管圧力を予測
することができる。 【００８３】これによって過渡時においても空燃比を目
標空燃比に制御すれば、加速応答性、ドライバビリテイ
および排気エミツシヨン等の悪化を防止することができ
る。 【００８４】 【００８５】 【発明の実施の形態】以下図面を参照して本発明の実施
の形態を詳細に説明する。図７は本実施の形態に係る内
燃機関の概略図である。 【００８６】エアクリーナ（図示せず）の下流側にはス
ロツトル弁８が配置されている。このスロツトル弁８に
は、スロツトル弁８の開度を検出するスロツトル開度セ
ンサ１０が取付けられている。スロツトル開度センサ１
０は、図８の等価回路に示すように、スロツトル弁８の
回動軸に固定された接触子１０Ｂと一端に電源が接続さ
れかつ他端が接地された可変抵抗１０Ａとで構成されて
おり、スロツトル弁８の開度が変化するに伴って、接触
子１０Ｂと可変抵抗１０Ａとの接触状態が変化し、スロ
ツトル弁８の開度に応じた電圧が接触子１０Ｂから得ら
れるように構成されている。スロツトル弁８の上流側の
吸気管壁には、吸入空気の温度を検出するサーミスタで
構成された温度センサ１４が取付けられている。スロツ
トル弁８の下流側にはサージタンク１２が配置されてお
り、サージタンク１２はインテークマニホールド１８、
吸気ポート２２および吸気弁２３を介して機関本体２０
の燃焼室２５に連通されている。このインテークマニホ
ールド１８には、各気筒に対応するように燃料噴射弁２
４が取付けられており、各気筒独立にまたは各気筒グル
ープ毎にまたは全気筒一斉に燃料を噴射できるように構
成されている。 【００８７】燃焼室２５は、排気弁２７、排気ポート２
６およびエキゾーストマニホールド２８を介して三元触
媒を充填した触媒装置（図示せず）に連通されている。
このエキゾーストマニホールド２８には、排ガス中の残
留酸素濃度を検出して理論空燃比に対応する値を境に反
転した信号を出力するＯ₂センサ３０が取付けられてい
る。 【００８８】シリンダブロツク３２には、ウオータジヤ
ケツト内に突出するように機関温度を代表する機関冷却
水温を検出するサーミスタ等で構成された冷却水温セン
サ３４が取付けられている。シリンダブロツク３６に
は、各々の燃焼室２５内に突出するように点火プラグ３
８が取付けられている。点火プラグ３８はデイストリビ
ユータ４０および点火コイルを備えたイグナイタ４２を
介してマイクロコンピユータ等で構成された制御回路４
４に接続されている。デイストリビユータ４０には、デ
イストリビユータシヤフトに固定されたシグナルロータ
とデイストリビユータハウジングに固定されたピツクア
ツプとで各々構成された気筒判別センサ４６および回転
角センサ４８が取付けられている。気筒判別センサ４６
は、例えば７２０°ＣＡ毎に気筒判別信号を出力し、回
転角センサ４８は、例えば３０°ＣＡ毎に回転角信号を
出力する。そして、この回転角信号の周期から機関回転
速度を演算することができる。 【００８９】マイクロコンピユータ等で構成された制御
回路４４は、図９に示すように、マイクロプロセツシン
グユニツト（ＭＰＵ）６０、リード・オンリ・メモリ
（ＲＯＭ）６２、ランダム・アクセス・メモリ（ＲＡ
Ｍ）６４、バツクアツプＲＡＭ（ＢＵ−ＲＡＭ）６６、
入出力ポート６８、入力ポート７０、出力ポート７２、
７４およびこれらを接続するデータバスやコントロール
バス等のバス７５を備えている。入出力ポート６８に
は、アナログ−デジタル（Ａ／Ｄ）変換器７８およびマ
ルチプレクサ８０が順に接続されており、このマルチプ
レクサ８０には、バツフア８２を介して吸気温センサ１
４が接続されると共に、バツフア８４およびバツフア８
５をそれぞれ介して水温センサ３４およびスロツトル開
度センサ１０が接続されている。また、入出力ポート６
８は、Ａ／Ｄ変換器７８およびマルチプレクサ８０に接
続されて、ＭＰＵからの制御信号に応じて吸気温センサ
１４、水温センサ３４及びスロツトル開度センサ１０出
力を順次所定周期でＡ／Ｄ変換するように制御する。 【００９０】入力ポート７０には、コンパレータ８８お
よびバツフア８６を介してＯ₂センサ３０が接続される
と共に波形整形回路９０を介して気筒判別センサ４６お
よび回転角センサ４８が接続されている。そして、出力
ポート７２は駆動回路９２を介してイグナイタ４２に接
続され、出力ポート７４は駆動回路９４を介して燃焼噴
射弁２４に接続されている。 【００９１】上記ＲＯＭ６２には、以下で説明する制御
ルーチンのプログラムや図１０に示すスロツトル開度Ｔ
Ａと機関回転速度ＮＥとで定められた定常状態での吸気
管圧力ＰＭＴＡのマツプ、図１１に示す機関回転速度Ｎ
Ｅと定常状態での吸気管圧力ＰＭＴＡ（またはスロツト
ル開度ＴＡ）とで定められた重みに関する係数ｎのマツ
プ、および実際の吸気管圧力ＰＭＳＭと機関回転速度Ｎ
Ｅとで定められた基本燃料噴射時間ＴＰのマツプが予め
記憶されている。図１０に示す定常状態での吸気管圧力
ＰＭＴＡのマツプは、スロツトル開度ＴＡと機関回転速
度ＮＥとを設定し、設定したスロツトル開度ＴＡと機関
回転速度ＮＥに対応する吸気管圧力を測定し、吸気管圧
力が安定したときの値を用いることにより作成される。
図１１に示す重みに関する係数ｎのマツプは、スロツト
ル弁をステツプ状に開いたときの吸気管圧力の応答（イ
ンデシヤル応答）時の時定数Ｔを測定し、この測定値と
図１３に示す演算ルーチンの実行周期Δｔsec とからＴ
／Δｔ（ ≒ｎ）を機関回転速度ＮＥと実際の吸気管圧
力ＰＭＴＡ（またはスロツトル開度ＴＡ）とに対応して
求めることにより作成される。そして図１２の基本燃料
噴射時間ＴＰのマツプは、機関回転速度と吸気管圧力と
を設定し目標空燃比となる基本燃料噴射時間ＴＰを測定
することにより作成される。 【００９２】次に、本実施の形態の繰り返し演算におけ
る吸気管圧力予測ルーチンを説明する前に、予測しない
吸気管圧力を用いて燃料噴射時間を演算するルーチンを
図１３に基いて説明する。このルーチンは、所定時間
（例えば、８msec）毎に実行される。ステツプ１００に
おいてＡ／Ｄ変換されたスロツトル開度ＴＡ（例えば、
８msec毎にＡ／Ｄ変換する）および機関回転速度ＮＥを
取込み、ステツプ１０２において図１０のマツプからス
ロツトル開度ＴＡと機関回転速度ＮＥに対応する定常状
態での吸気管圧力ＰＭＴＡを演算する。次のステツプ１
０４では、ステツプ１０２で演算された吸気管圧力ＰＭ
ＴＡとステツプ１００で取込んだ機関回転速度ＮＥとに
基づいて図１１に示すマツプから重みに関する係数ｎを
演算する。なお、スロツトル開度と機関回転速度で重み
に関する係数ｎのマツプを定めた場合には、ステツプ１
０４でステツプ１００で取込んだスロツトル開度ＴＡと
機関回転速度ＮＥとで重みに関する係数ｎを演算するよ
うにしてもよい。次のステツプ１０６では、ステツプ１
０２で演算された吸気管圧力ＰＭＴＡとステツプ１０４
で演算された重みに関する係数ｎと前回のこのルーチン
の実行時にステツプ１０６で演算された前回の加重平均
値ＰＭＳＭ_i-1とを用いて上記で説明した(23)式に従っ
て今回の加重平均値ＰＭＳＭ_iを演算する。次のステツ
プ１０８では今回の加重平均値ＰＭＳＭ_iと機関回転速
度ＮＥとに基づいて図１２に示すマツプから基本燃料噴
射時間ＴＰを演算する。そして、次のステツプ１１０に
おいて吸気温や機関冷却水温等で定まる補正係数ＦＫを
基本燃料噴射時間ＴＰに乗算することによって燃料噴射
時間ＴＡＵ演算する。そして、図示しない制御ルーチン
において所定クランク角になったときに燃料噴射時間Ｔ
ＡＵに相当する時間燃料噴射弁を開弁して燃料噴射を実
行する。 【００９３】図１４は所定クランク角毎に割込みによっ
て点火進角θを演算するルーチンを示すものである。な
お、図１４において図１３と同一部分については同一符
号を付して説明を省略する。ステツプ１１２では、今回
演算された加重平均値ＰＭＳＭ_iと機関回転速度ＮＥと
により基本点火進角θ _BASE を演算する。この基本点火進
角θ _BASE は、演算式によって演算してもよく、また基本
燃料噴射時間と同様にマツプを作成してこのマツプから
演算するようにしてもよい。そして、次のステツプ１１
４において基本点火進角θＢＡＳＥに吸気温や機関冷却
水温等で定まる補正係数ＩＫを乗算して点火進角θを求
める。そして図示しない点火時期制御ルーチンにおいて
基本点火進角θでイグナイタをオフすることにより点火
を実行する。 【００９４】図１５（１）、（２）に、加速時における
従来での加速増量を行なわない場合の空燃比の変化と本
実施の形態での空燃比の変化とを比較して示すと共に燃
料噴射量を求めるための本実施の形態での加重平均値Ｐ
ＭＳＭと従来の検出した吸気管圧力ＰＭとの相異を示
す。図１５から理解されるように、従来例の空燃比は加
速時にリーンスパイクが生じているが、本実施の形態の
空燃比は略フラツトになっている。 【００９５】以上説明したように、本参考例では実際の
吸気管圧力を予測して燃料噴射量と点火時期とを制御す
ることにより圧力センサやフイルタを用いることなく、
精度良い燃料噴射量制御と点火時期制御を行なうことが
できる。 【００９６】次に上記内燃機関に本発明を適用した第１
の実施の形態を説明する。この実施の形態は加重平均値
の演算を所定回繰り返すことによって吸入空気量確定時
（吸気弁全閉時）の実際の吸気管圧力を予測するもので
ある。なお、本実施の形態では更にこの予測した吸気管
圧力によって燃料噴射量を制御するようにしたものであ
る。図１６は所定時間（本実施の形態では８msec）毎に
実行して吸入空気量確定時の吸気管圧力の予測値ＰＭＳ
Ｍ２を演算するルーチンを示すものである。ステツプ２
００において機関回転速度ＮＥを取込むと共に、スロツ
トル開度ＴＡのＡ／Ｄ変換を行なってスロツトル開度Ｔ
Ａを取込む。ステツプ２０２では図１０に示すマツプか
ら機関回転速度ＮＥとスロツトル開度ＴＡとに対応する
定常状態での吸気管圧力ＰＭＴＡを演算する。次のステ
ツプ２０４では図１１に示すマツプから重み付けに関す
る係数ｎを演算する。次のステツプ２０６とステツプ２
０８では、レジスタＰＭＳＭ１に記憶されている前回演
算した加重平均値ＰＭＳＭ_i-1をＲＡＭから読出して上
記(23)式に基づいて今回の加重平均値ＰＭＳＭ_iを演算
し、ステツプ２１０においてこの加重平均値ＰＭＳＭ_i
をレジスタＰＭＳＭ１に記憶しておく。次のステツプ２
１２では、現在時刻から吸気管圧力予測時点までの時間
Ｔmsecを図１６のルーチンの演算周期Δｔ（＝８msec）
で除算することにより演算回数Ｔ／Δｔを演算する。こ
の予測時間Ｔmsecは、現在時刻から吸入空気量確定まで
の時間すなわち現在時刻から吸気弁が閉じるまでの時間
を採用することができ、各気筒独立に燃料を噴射しない
場合には燃料噴射弁から燃焼室までの燃料の飛行時間等
も考慮して決定されるが、現時点から予測先までのクラ
ンク角が同一であってもこの予測時間Ｔmsecは機関回転
速度が速くなると短くなるので機関回転速度等の運転条
件によって可変することが好ましい（例えば、機関回転
速度が速くなるに従って短くする）。次のステツプ２１
４では、演算回数Ｔ／Δｔ回上記(23)式の演算を繰り返
して実行し、ステツプ２１６においてこの演算した値を
吸気管圧力の予測値ＰＭＳＭ２とする。このように加重
平均値を繰り返して実行することにより最新の加重平均
値は定常運転状態での吸気管圧力値に近づくので、加重
平均値の演算回数を上記のように定めることにより現在
時刻からＴmsec先の吸気管圧力（現時点より定常状態に
近い状態での吸気管圧力）を予測することができる。 【００９７】図１７は所定クランク角（例えば、１２０
°ＣＡ）毎に燃料噴射時間ＴＡＵを演算するルーチンを
示すもので、機関回転速度ＮＥとステツプ２１６で演算
された吸気管圧力の予測値ＰＭＳＭ２とに基づいて図１
２に示すマツプから基本燃料噴射時間ＴＰを演算する。
そして、ステツプ２２０において上記のステツプ１１０
と同様にして燃料噴射時間ＴＡＵを演算する。 【００９８】なお、現在時刻からＴmsec経過した時点で
はスロツトル開度や機関回転速度が変化することがある
ため、スロツトル開度の微分値や機関回転速度の微分値
を用いてＴmsec先のスロツトル開度や機関回転速度を予
測してＴmsec先の定常状態での吸気管圧力を予測し、上
記の加重平均値の演算を繰り返すようにすれば、さらに
精度は向上する。 【００９９】上記のように演算したときの加重平均値お
よびＴmsec経過後の予測値ＰＭＳＭ２を図１８および図
１９に示す。図１８では、１６msec先の予測値と理論値
とを示したが、予測値は理論値と略等しくなっている。
なお、スロツトル開度のＡ／Ｄ変換タイミングは燃料噴
射時間演算タイミングと一致する場合もあるが、最大演
算周期Δｔに相当する時間ずれる。従って、このずれ時
間を平均（０＋Δｔ）／２して Ｔ± Δｔ／２時間先
の吸気管圧力を予測するようにしても良い。 【０１００】次に、第２実施の形態を説明する。本実施
の形態は燃料の機関壁面付着量を予測し燃料噴射量を補
正するようにしたものである。 【０１０１】機関燃焼室へ吸入されないで機関壁面に付
着している燃料付着量は、吸気弁閉弁時の吸気管圧力に
よって定まり、例えば、吸気管圧力がＰＭ１の状態から
ＰＭ２の状態へ加速した場合、各々の吸気管圧力での燃
料付着厚さをＴ１、Ｔ２とすると、燃料の付着厚さをＴ
１からＴ２に増加させるために必要な、壁面への燃料供
給量は、スロツトル開放速度および燃料噴射回数等に関
係なく定まる。そこで、本実施の形態では、ある基準吸
気管圧力（例えば、０mmHgabs ）から任意の吸気管圧力
まで変化させたときの壁面に供給すべき噴射量の総付着
量を図２２に示すように吸気弁全閉時での吸気管圧力に
対して予めマツプの形でＲＯＭに記憶させておく。 【０１０２】図２０は、本実施の形態の所定クランク角
（３６０°ＣＡ）毎に実行される燃料噴射量演算ルーチ
ンを示すもので、ステツプ２３０においては上記図１６
で演算された吸気管圧力の予測値ＰＭＳＭ２と機関回転
速度ＮＥとから基本燃料噴射時間ＴＰを上記と同様に演
算する。次のステツプ２３２では、吸気温や機関冷却水
温等によって定まる燃料噴射量の補正係数ＦＫを算出す
る。次のステツプ２３４では、図２２のマツプから吸気
管圧力の予測値ＰＭＳＭ２に対応する機関壁面への燃料
付着量ＦＭＷＥＴを算出する。そして次のステツプ２３
６において基本燃料噴射時間と補正係数ＦＫとを乗算す
ると共に、今回求めた燃料の付着量ＦＭＷＥＴから前回
燃料の付着量ＦＭＷＥＴ _OLD を減算した値を補正加算値
として加算することにより燃料噴射時間ＴＡＵを求め
る。この補正加算量は吸気管圧力の変化によって生ずる
燃料付着量の変化量を表わしている。そしてステツプ２
３８において今回求めた燃料の付着量ＦＭＷＥＴを前回
の付着量ＦＭＷＥＴ _OLD としてＲＡＭに記憶する。 【０１０３】以上のように燃料噴射量を制御することに
よって、図２１に示すように斜線で示す部分の量の燃料
が増量され、これによって機関の内壁面に燃料付着厚さ
だけの燃料が付着しても補正加算量によって機関に供給
する燃料は要求値になる。なお、図２４はスロツトル開
度、吸気管圧力の予測値および空燃比の変化を示すもの
であり、本実施の形態では破線で示す従来例のようにリ
ーンスパイクは生ぜず空燃比の変動が少なくなってい
る。 【０１０４】次に第３実施の形態を説明する。上記の第
２実施の形態では噴射毎の燃料付着量によって燃料噴射
量を制御するようにしてたが、機関壁面への燃料の付着
は瞬時には安定しないことを考慮して本実施の形態では
各噴射での補正加算量を時間減衰させることによって次
回以降の噴射にも反映させることにより燃焼室への燃料
の供給量を要求値と等しくするようにしている。図２５
は本実施の形態の燃料噴射演算ルーチンを示すもので、
例えば所定クランク角（３６０°ＣＡ）毎に実行され
る。なお、図２５において図２０と同一部分には同一符
号を付して説明を省略する。ステツプ２３４で燃料付着
量ＦＭＷＥＴを算出した後は、ステツプ２４０において
以下の式に従って補正加算量ＦＡＥを算出する。 ＦＡＥ＝０．２・ＦＡＥ _OLD ＋ＦＭＷＥＴ −ＦＭＷＥＴ _OLD …(24) なお、ＦＡＥ _OLD は前回演算した補正加算量、ＦＭＷＥ
Ｔ _OLD は前回演算した燃料の壁面への付着量である。 【０１０５】上記(24)式では、前回の補正加算量ＦＡＥ
_OLD に０. ２を乗算しているので、前回の補正加算量を
８０％減衰させて前回の補正加算量の２０％を今回の補
正加算量に反映させている。なお、この減衰の仕方は機
関によって最適な方法が選択され、上記のように所定ク
ランク角（上記の例では３６０°ＣＡ）毎に所定量づつ
減衰させても良く、また所定時間毎に所定量づつ減衰さ
せるようにしても良い。 【０１０６】次のステツプ２４２では、上記と同様にし
て基本燃料噴射時間と補正係数ＦＫと補正加算量ＦＡＥ
とを用いて燃料噴射時間ＴＡＵを演算する。そして、ス
テツプ２４４において補正加算量ＦＡＥを前回の補正加
算量ＦＡＥ _OLD としてＲＡＭに記憶すると共に、燃料付
着量ＦＭＷＥＴを前回の燃料付着量ＦＭＷＥＴ _OLD とし
てＲＡＭに記憶する。 【０１０７】なお、上記図２２では、吸気弁全閉状態で
の吸気管圧力に応じて燃料付着量を定める例について説
明したが、燃料付着量は機関回転速度に応じても変化す
るため図２３に示すように吸気管圧力と機関回転速度を
変数として変化するマツプとして記憶させても良い。ま
た、燃料付着量は機関温度によっても変化し、機関温度
が低い程燃料付着量が多くなるので更にこの機関温度を
変数として定めるようにしても良い。また、上記実施の
形態では、加重平均値によって吸気管圧力を予測する例
について説明したが、上記（16）式に従って吸気管圧力
を予測しても良く、定常状態の吸気管圧力を１次遅れ要
素で処理して吸気管圧力を予測しても良い。 【０１０８】 【発明の効果】以上説明したように本発明は、所定期間
先の吸気管圧力の予測値を求めるようにしたので、所定
期間として現在から吸入空気量が確定するまでの期間を
取れば、吸入空気量確定時の実際の吸気管圧力を予測で
きる、という効果を有する。

【図面の簡単な説明】 【図１】参考例の原理を説明するための線図である。 【図２】吸気系内の実際の吸気管圧力の時間に対する変
化を示す線図である。 【図３】従来のスロツトル開度と機関回転速度とで定ま
る吸気管圧力と実際の吸気管圧力との相異を示す線図で
ある。 【図４】従来のスロツトル開度と機関回転速度とで定ま
る燃料噴射量と要求燃料噴射量との相異を示す線図であ
る。 【図５】他の参考例を説明するためのブロツク図であ
る。 【図６】更に他の参考例を説明するためのブロツク図で
ある。 【図７】内燃機関を示す概略図である。 【図８】スロツトル開度センサの等価回路図である。 【図９】図８の制御回路の詳細を示すブロツク図であ
る。 【図１０】定常状態での吸気管圧力のマツプを示す線図
である。 【図１１】加重平均値の重み付けに関する係数のマツプ
を示す線図である。 【図１２】基本燃料噴射時間のマツプを示す線図であ
る。 【図１３】参考例の燃料噴射量演算ルーチンを示す流れ
図である。 【図１４】参考例の点火進角演算ルーチンを示す流れ図
である。 【図１５】（１）、（２）は従来例と参考例との空燃比
および吸気管圧力の変化を示す線図である。 【図１６】第１の実施の形態の吸気管圧力の予測値を演
算するルーチンを示す流れ図である。 【図１７】第１の実施の形態の燃料噴射時間演算ルーチ
ンを示す流れ図である。 【図１８】第１の実施の形態の吸気管圧力の予測値等の
変化を示す線図である。 【図１９】第１の実施の形態の吸気管圧力の予測値等の
変化を示す線図である。 【図２０】第２の実施の形態の燃料噴射時間演算ルーチ
ンを示す流れ図である。 【図２１】燃料の壁面付着厚さと吸気管圧力との関係を
示す線図である。 【図２２】補正噴射量のマツプを示す線図である。 【図２３】補正噴射量のマツプを示す線図である。 【図２４】第２の実施の形態の空燃比等の変化を従来例
と比較して示す線図である。 【図２５】第３の実施の形態の燃料噴射量演算ルーチン
を示す流れ図である。 【符号の説明】 ８ スロツトル弁 １０ スロツトル開度センサ ４８ 回転角センサ

Claims (1)

1. (57)【特許請求の範囲】 １．以下の式にしたがって演算された加重平均値ＰＭＳ
   Ｍ _i を次回の演算において以下の式の過去に演算された
   加重平均値ＰＭＳＭ _i-1 として用いて所定回数繰り返し
   演算して加重平均値ＰＭＳＭ _i を求め、繰り返し演算後
   の加重平均値ＰＭＳＭ _i を吸気管圧力予測時点での吸気
   管圧力として予測する内燃機関の吸気管圧力予測方法。 ＰＭＳＭ _i ＝（（ｎ−１）・ＰＭＳＭ _i-1 ＋ＰＭＴＡ）/ｎ ただし、ＰＭＴＡは、スロットル開度と機関回転速度と
   に基づいて演算されると共に、繰り返し演算中は一定と
   される定常状態での吸気管圧力、 ｎは、過渡時の吸気管圧力の変化に関する時定数と前記
   演算周期とを用いて演算される重みに関する係数であ
   る。 ２．前記定常状態での吸気管圧力ＰＭＴＡは、スロット
   ル開度の微分値及び機関回転速度の微分値の少なくとも
   一方を用いて演算される吸気管圧力予測時点における定
   常状態での吸気管圧力である請求項１記載の内燃機関の
   吸気管圧力予測方法。 ３．前記演算周期を一定とし、スロットル開度と機関回
   転速度とに基いて、又は、定常状態での吸気管圧力と機
   関回転速度とに基いて、前記重みに関する係数ｎを演算
   する請求項１または２記載の内燃機関の吸気管圧力予測
   方法。 ４．前記所定回数を、現在時刻から吸気管圧力予測時点
   までの時間を演算周期で除算して得られる値とした請求
   項１〜３のいずれか１項記載の内燃機関の吸気管圧力予
   測方法。 ５．スロットル開度を検出するスロットル開度検出手段
   と、 機関回転速度を検出する回転速度検出手段と、 前記検出されたスロットル開度と前記検出された機関回
   転速度とに基づいて所定の演算周期で定常状態での吸気
   管圧力ＰＭＴＡを演算する吸気管圧力演算手段と、前記演算周期を一定とし、演算された定常状態での吸気
   管圧力と検出された機関回転速度とに基いて、または検
   出されたスロットル開度と検出された機関回転 速度とに
   基づいて、過渡時の吸気管圧力の変化に関する時定数と
   前記演算周期とで定まる重みに関する係数ｎを演算する
   重み係数演算手段と、 繰り返し演算中は定常状態での吸気管圧力ＰＭＴＡを一
   定とすると共に、以下の式にしたがって演算された加重
   平均値ＰＭＳＭ _i を次回の演算において以下の式の過去
   に演算された加重平均値ＰＭＳＭ _i-1 として用て所定回
   数繰り返し演算して加重平均値ＰＭＳＭ _i を求め、繰り
   返し演算後の加重平均値ＰＭＳＭ _i を吸気管圧力予測時
   点での吸気管圧力として予測する吸気管圧力予測手段
   と、 を備えた内燃機関の吸気管圧力予測装置。 ＰＭＳＭ _i ＝（（ｎ−１）・ＰＭＳＭ _i-1 ＋ＰＭＴＡ）/ｎ ６．前記定常状態での吸気管圧力ＰＭＴＡは、スロット
   ル開度の微分値及び機関回転速度の微分値の少なくとも
   一方を用いて演算される吸気管圧力予測時点における定
   常状態での吸気管圧力である請求項５記載の内燃機関の
   吸気管圧力予測装置。 ７．前記所定回数を、現在時刻から吸気管圧力予測時点
   までの時間を演算周期で除算して得られる値とした請求
   項５または６記載の内燃機関の吸気管圧力予測装置。