(1),GNはPM及び(PA−PM)の関数となる。KTP中の負荷のパラメータは(PA−PM)のみである。EGR)の影響を受けずに吸入空気量、ひいては基本噴射時間TPを算出することができる。その結果、単に吸気圧PMから吸入空気量(噴射時間TAU)を求めた場合に比べ、CPU92による噴射時間TAUの算出、及びインジェクタ25の駆動制御の精度が向上する。2実施例を図13に基づき説明する。第2実施例は、CPU92による噴射時間算出ルーチンの処理内容、特に充填効率(補正係数KTP)を求めるための処理内容が上記第1実施例と異なっている。吸気圧PMと大気圧PAとの比(PM/PA)が補正係数KTPの演算に用いられている。これは以下の理由による。strokeとして作用する。吸気バルブ11が開弁した瞬間(バルブオーバラップ開始時)の燃焼室8内の燃焼ガスの圧力は、大気圧にほぼ等しい。この燃焼ガスの圧力は、燃焼室8で再吸入される際に大気圧PAから吸気圧PMに転ずる。従って、内部EGRの量は、この圧力変化から大きく影響を受けるものと考えられる。change.【0101】ここで、吸気バルブ11が開弁した瞬間のシリンダ4内の容積をVとすると、吸気行程で燃焼ガスの圧力が大気圧PAから吸気圧PMに変化する。この式(9)を上記式(7),(8)にそれぞれ代入すると、下記式(7a),(8a)が得られる。(PA−PM)の関数ではなく、比(PA/PM)の関数となる。PM).【0103】このような観点から、第2実施例では比(PM/PA)が用いられて補正係数KTPが求められる。Throttle、スロットルバルブの全開状態)時を「1.0」とする負荷率と定義でき、値の意味が理解しやすくなるからである。噴射時間算出ルーチン以外の箇所であるエンジン1やその周辺機器等の構成は第1実施例と同様であるので、ここではその説明を省略する。9のステップ101と同様に、エンジン1の始動時であるか否かを判定する。CPU92はエンジン始動時であり吸入空気量の正確な検出が困難であると判断し、ステップ202へ移行する。この算出が終了すると、CPU92はこのルーチンを一旦終了する。routine.【0105】ステップ201の判定条件が成立していないと、CPU92はエンジン1の運転状態が始動後の状態であると判断し、ステップ203においてエンジン回転速度NEが予め設定した判定値α以上であるか否かを判定する。ここで、NE≧αの領域は、バルブオーバラップ時に圧力変化をともないながら燃焼室8から吸気通路23へ戻された燃焼ガスが次回燃焼予定の混合気に混入される第2の領域に相当する。NE<αの領域は、バルブオーバラップ時に排気通路32へ一旦排出された燃焼ガスが背圧PEXと吸気圧PMとの差圧により吸気通路23へ戻され、次回燃焼予定の混合気に混入される第1の領域に相当する。9E及び変位角θを読み込むとともに、大気圧PAと吸気圧PMとの差圧(ブースト圧Pb)を求める。NE,Pb,θに対応する補正係数KTPを求める。theta.【0107】一方、ステップ203の判定条件が成立している(NE≧α)と、CPU92はエンジン回転速度NE及び変位角θを読み込むとともに、吸気圧PMと大気圧PAとの比を求める。CPU92は図示しないマップを用い、前記の各値NE,PM/PA,θに対応する補正係数KTPを求める。9〜106と同様の処理を実行する。CPU92はステップ206で上記式(1)に従って、1回の吸入行程でシリンダ4に充填される空気の質量GNを演算する。208において上記式(4)に従って噴射時間TAUを算出し、このルーチンを一旦終了する。routine.【0109】このようにCPU92は、充填効率(補正係数KTP)を演算する(ステップ204,205)ための充填効率演算手段として機能する。for.【0110】CPU92は噴射時間TAUを算出すると、別途用意されたルーチンにおいてその値TAUに応じた駆動信号を外部出力回路97を介してインジェクタ25に出力する。この信号に応じインジェクタ25の開弁時間が制御され、その時間に応じた量の燃料が同インジェクタ25から噴射される。/PA)がパラメータとして加えられる。/PA)が加味されることにより、吸気バルブ11開弁時に燃焼ガスが大気圧PAから吸気圧PMに変化し、内部EGRの量が変化しても、充填効率が精度良く求められる。embodiment,吸気バルブ11の任意の開閉タイミングでの充填効率に応じた吸入空気量を演算することができる。11.【0112】さらに、第2実施例ではエンジン1の運転状態に応じて補正係数KTPの演算方法が切換えられる。CPU92によって実行される処理中には、エンジン回転速度NEが判定値α未満である場合に差圧(PA−PM)に基づき補正係数KTPを演算する処理(ステップ205)と、エンジン回転速度NEが判定値α以上である場合に比(PM/PA)に基づき補正係数KTPを演算する処理(ステップ204)とが含まれる。NEにより変化するが、第2実施例ではその内部EGRの量に応じた吸入空気量を全運転領域にわたり精度良く演算することができる。3実施例を図14,15に従って説明する。第3実施例では、エンジン1の運転状態に応じた目標変位角θVTAが最遅角時の値(θa)に設定されたときのその値θaと実際の変位角θとの偏差を求め、この偏差を学習値GVTとする。θVTAとが一致するようにVVT46のリングギヤ58の位置を調整する。補正後の変位角θVTを用いて充填効率及び空気質量GNを求める。この空気質量GNを用いて点火時期を算出している。これらが前記第1及び第2実施例と大きく異なっている点である。エンジン1、その周辺機器等の構成は第1実施例と同様であるので、ここではその説明を省略する。θVTを算出するとともに、その変位角θVTを用いて吸気バルブ11の開閉タイミングを調整するためのバルブタイミング制御ルーチンを示している。図15のフローチャートは、前記バルブタイミング制御ルーチンで算出した補正後の変位角θVTを利用して点火時期を算出するための点火時期算出ルーチンを示している。routine.【0114】まず、図14のバルブタイミング制御ルーチンついて説明する。CPU92はステップ301において、カム角センサ42によるカム角信号SG2の発生時の回転速度センサ40による回転速度信号SG1から、予め定められた基準の回転速度信号までに要する時間を計測し、その時間に基づきクランクシャフト7に対するカムシャフト13の実際の回転位相(変位角θ)を算出する。7.【0115】続いて、CPU92はステップ302において、目標変位角θVTAが最遅角時の値θa(例えば、零)であるか否かを判定する。θVTAは、スロットル開度TA、エンジン回転速度NE等に基づき、予め規定されたマップを参照して求められたものである。map.【0116】最遅角時の目標変位角θaは、吸気バルブ11の開閉タイミングを最も遅くするときのクランクシャフト7に対するカムシャフト13の回転位相である。11.この際、センサの取付け誤差、経時変化等がなく、実際の変位角θがこの目標変位角θaに一致すれば、リングギヤ58が第1の位置(リングギヤ58がカバー55に当接する位置)に保持されるはずである。θと目標変位角θaとが一致しなくなる場合が起こる。θが目標変位角θVTA(=θa)に収束してから所定時間が経過したか否かを判定する。CPU92はこの判定条件が成立していると、変位角θの変動が少なくほぼ一定の値になっているものと判断して学習モードを設定し、ステップ304において、前記変位角θと最遅角時の目標変位角θaとの偏差を求める。この偏差は、センサの取付け誤差、経時変化等に起因する回転位相のばらつきに相当する。ここで、目標変位角θaが零であれば変位角θが学習値GVTとなる。この学習値GVTは、目標変位角θVTAが最遅角時の値θa(零)になる毎に更新される。mode,ステップ304の処理は行わずに、ステップ305へ移行する。305.【0119】そして、CPU92はステップ305において、前記変位角θから学習値GVTを減算して同変位角θを補正する。従って、ステップ305で得られる補正後の変位角θVTは、センサの取付け誤差等の回転位相のばらつきを排除した精度の高い値となる。CPU92は、ステップ306において補正後の変位角θVTと目標変位角θVTAとが一致するようにVVT46のリングギヤ58の位置を調整する。目標変位角θVTAと補正後の変位角θVTとの偏差を求め、予め規定されたマップを参照し、あるいは所定の演算式に従いその偏差に対応した、LSV78を制御するための目標制御量DVを算出する。CPU92は、この目標制御量DVに基づいてLSV78の開度をフィードバック制御する。すると、VVT46の両油圧室59,61へ供給される油圧が調整され、リングギヤ58が移動する。moves.吸気バルブ11の開閉タイミングがエンジン1の運転状態に応じて連続的に変更され、もってバルブオーバラップが連続的に調整される。306の処理を実行した後にこのルーチンを一旦終了する。θVTAと実際の変位角θとの間にずれが存在しても、そのずれを、学習値GVTを用いた補正によって吸収し、精度の高い変位角θVTを求めることができる。And,この補正後の変位角θVTをLSV78の開度のフィードバック制御に用いるので、例えば過渡運転時において変位角θが目標変位角θVTAに収束するのに要する時間が短くなる等、目標変位角θVTAへの制御精度が向上する。(GVT)を求めるようにしている。このため、ほぼ一定値となった変位角θを用い、補正後の変位角θVTを正確に求めることができる。accurately.【0122】次に、図15の点火時期算出ルーチンについて説明する。CPU92はステップ401〜404において、前述した図13のステップ203〜206と同様の処理を実行する。CPU92はステップ401において、エンジン回転速度NEが予め設定した判定値α以上であるか否かを判定する。advance.【0123】ステップ401の判定条件が成立していない(NE<α)と、CPU92はステップ403において、エンジン回転速度NE及び補正後の変位角θVTを読み込むとともに、大気圧PAと吸気圧PMとの差圧(ブースト圧Pb)を求める。変位角θVTは前述したバルブタイミング制御ルーチンのステップ305で求めた値である。correctedに置き換えられる。補正係数KTPは上述したように充填効率に相当するものである。θVTを読み込むとともに、吸気圧PMと大気圧PAとの比を求める。403で補正係数KTPを求めると、ステップ404で上記式(1)に従って、1回の吸入行程でシリンダ4に充填される空気の質量GNを演算する。GN,ROM93内に格納されている記憶データから基本点火進角ABSEを選び出す。続いて、ステップ406において前記基本点火進角ABSEに各種補正値AKを加算し、その加算結果を点火時期AOPとして設定する。各種補正値AKとしては、例えば、冷却水温THWが低いアイドル時に暖機性能を向上させるために点火時期を進めるための補正値がある。その他にも、暖機後のアイドル時に、エンジン回転速度NEが目標値よりも低くなったときに、両者の偏差に応じて点火時期を進めるための補正値がある。点火時期AOPの算出後、CPU92はこのルーチンを一旦終了する。CPU92は充填効率(補正係数KTP)を演算する(ステップ402,403)ための充填効率演算手段として機能する。従って、第3実施例によっても第2実施例と同様に、充填効率の変化に起因する吸入空気量のばらつきを吸収し、吸気バルブ11の任意の開閉タイミングでの充填効率に応じた吸入空気量を演算することができる。calculated.【0128】また、点火時期算出ルーチンでは、バルブタイミング制御ルーチンで求められた精度の高い補正後の変位角θVTを、充填効率(補正係数KTP)の演算に用いている。このため、補正前の変位角θを用いた場合に比べ、精度の高い充填効率及び吸入空気量が得られ、適正な点火時期を算出することができる。28に一次電流の遮断信号を出力して点火時期を制御する。本発明は次のような別の実施例に具体化することができる。(PA-PM)を1つのセンサで検出し、その値をブースト圧Pbとしてもよい。Pb.【0131】(2)大気圧PAを間接的に求めてもよい。具体的な方法としては、スロットル開度TAが全開近傍にあるときの吸気圧センサ37の出力を大気圧とみなしてもよい。KTP0〜図12のマップに代えて、エンジン1回転当たりの空気質量GNや基本噴射時間TPを直接求めるためのマップを作成してもよい。14以外のバルブ駆動装置、例えば電磁駆動バルブ等、を有するエンジンにも適用可能である。by,バルブオーバラップを調整するようにしてもよい。VVT,バルブオーバラップを調整するようにしてもよい。overlap.【0134】(6)図13のステップ203,205の処理を省略してもよい。401,403の処理を省略してもよい。エンジン始動後であればエンジン回転速度NEの高低によらず、比(PM/PA)を用いて補正係数KTPを求めるようにしてもよい。used.【0135】(7)吸入空気量算出に当たって、充填効率を算出することなく、エンジン回転速度NE、吸気圧PM及びバルブタイミングから直接マップによって吸入空気量を求めるようにしてもよい。timing.【0136】(8)第3実施例において、最進角時の目標変位角と実際の変位角θとの偏差を求め、この偏差を学習値GVTとしてもよい。以上、本発明の各実施例について説明したが、各実施例から把握できる請求項以外の技術的思想について、以下にそれらの効果とともに記載する。このような構成とすることにより、補正後の作動タイミングを正確に求めることができる。respectively,その検出値に基づき充填効率又はその相関値を求める。value.この充填効率又は相関値と吸気圧とに基づき気筒への吸入空気量を求めるようにしている。pressure.【0139】このため、吸気バルブの閉じタイミングの変化や、バルブオーバラップの変化等にともない充填効率が変動しても、その充填効率又は相関値を精度良く求め、バルブの任意の開閉タイミングでの充填効率に応じた吸入空気量を演算することができる。timing,実際の作動タイミングとの偏差を求める。同偏差に基づきタイミング検出手段による作動タイミングを補正して精度を高め、その補正後の作動タイミングと目標作動タイミングとが一致するようにタイミング調整部材の位置を調整する。前記補正後の精度の高い作動タイミングを用いて充填効率を演算するようにしている。このため、充填効率及び吸入空気量をさらに精度良く算出することができる。