気筒と、前記気筒に設けられた吸気弁とを有し、間欠吸気を行うエンジンと、前記吸気弁が開弁状態となるタイミングでの吸入空気量と、前記吸気弁が閉弁状態となるタイミングでの吸入空気量とを検出する吸入空気量検出手段と、前記吸気弁が開弁状態となるタイミングでの吸入空気量と、前記吸気弁が閉弁状態となるタイミングでの吸入空気量との差分を算出する差分算出手段と、前記吸気弁が開弁してから閉弁するまでの時間幅を算出する時間幅算出手段と、前記差分及び前記時間幅に基づき、吸気通路の脈動波形を算出する脈動波形算出手段と、前記脈動波形に基づき、前記気筒の燃焼室に吸入されると推定される吸入空気量を算出する吸入空気量推定手段と、を備えることを特徴とする内燃機関の制御装置。前記気筒は複数存在し、前記脈動波形算出手段は、前記気筒の各々の前記時間幅と、前記気筒の各々の吸気弁が開弁状態となるタイミングと、前記気筒の各々の吸気弁が閉弁状態となるタイミングとに基づき前記脈動波形を算出する請求項1に記載の内燃機関の制御装置。前記脈動波形は、吸入空気量の時間関数であり、前記吸入空気量推定手段は、前記脈動波形が示す吸入空気量を、前記吸気弁が開弁している期間で積分することで、前記吸入空気量を算出する請求項1または2に記載の内燃機関の制御装置。前記エンジンは、360度位相がずれた2気筒エンジンであることを特徴とする請求項1乃至3のいずれか一項に記載の内燃機関の制御装置。360度位相が異なる2気筒エンジンなどの間欠吸気を行うエンジンでは、吸気脈動に起因して、吸入空気量を正確に算出することができず、充填効率の実値と算出値とが乖離してしまうという問題があった。本発明は、上記のような課題を解決するためになされたものであり、吸気脈動を考慮した吸入空気量を適切に算出することが可能な内燃機関の制御装置を提供することを主な目的とする。本発明の1つの観点では、気筒と、前記気筒に設けられた吸気弁とを有し、間欠吸気を行うエンジンと、前記吸気弁が開弁状態となるタイミングでの吸入空気量と、前記吸気弁が閉弁状態となるタイミングでの吸入空気量とを検出する吸入空気量検出手段と、前記吸気弁が開弁状態となるタイミングでの吸入空気量と、前記吸気弁が閉弁状態となるタイミングでの吸入空気量との差分を算出する差分算出手段と、前記吸気弁が開弁してから閉弁するまでの時間幅を算出する時間幅算出手段と、前記差分及び前記時間幅に基づき、吸気通路の脈動波形を算出する脈動波形算出手段と、前記脈動波形に基づき、前記気筒の燃焼室に吸入されると推定される吸入空気量を算出する吸入空気量推定手段と、を備える。上記の内燃機関の制御装置は、車両に搭載され、エンジンと、吸入空気量検出手段と、差分算出手段と、時間幅算出手段と、脈動波形算出手段と、吸入空気量推定手段とを備える。エンジンは、気筒を有し、間欠吸気を行う。吸入空気量検出手段は、例えばECUElectronicControlUnitであり、エアフロメータなどにより、吸気弁が開弁全開状態となるタイミングでの吸入空気量と、前記吸気弁が閉弁全閉状態となるタイミングでの吸入空気量とを検出する。差分算出手段は、例えばECUであり、吸気弁が開弁状態となるタイミングでの吸入空気量と、吸気弁が閉弁状態となるタイミングでの吸入空気量との差分を算出する。時間幅算出手段は、例えばECUであり、吸気弁が開弁してから閉弁するまでの時間幅を算出する。脈動波形算出手段は、吸入空気量の差分と、吸気弁が閉弁するまでの時間幅とに基づき、吸気通路の脈動波形を算出する。吸入空気量推定手段は、例えばECUであり、算出した脈動波形に基づき、気筒の燃焼室に吸入されると推定される吸入空気量を算出する。一般に、吸気脈動が生じると、吸気通路内で検出した吸入空気量が逆流の影響により誤差が生じる。これを勘案し、上記の内燃機関の制御装置は、開弁状態となるタイミングでの吸入空気量と、吸気弁が閉弁状態となるタイミングでの吸入空気量との差分を用いることで、逆流の影響を好適に除外し、かつ、吸気弁が開弁してから閉弁するまでの時間幅を用いて、脈動波形を的確に求め、吸入空気量の推定値を高精度に算出することが可能となる。上記の内燃機関の制御装置の一態様では、前記気筒は複数存在し、前記脈動波形算出手段は、前記気筒の各々の前記時間幅と、前記気筒の各々の吸気弁が開弁状態となるタイミングと、前記気筒の各々の吸気弁が閉弁状態となるタイミングとに基づき前記脈動波形を算出する。一般に、各気筒の吸気弁が開弁状態の場合には脈動波形は極大となり各気筒の吸気弁が閉弁状態の場合には脈動波形は極小となる。従って、この態様により、内燃機関の制御装置は、好適に、脈動波形を算出することができる。上記の内燃機関の制御装置の他の一態様では、前記脈動波形は、吸入空気量の時間関数であり、前記吸入空気量推定手段は、前記脈動波形が示す吸入空気量を、前記吸気弁が開弁している期間で積分することで、前記吸入空気量を算出する。これにより、内燃機関の制御装置は、好適に、気筒内に流入した吸入空気量の推定値を算出し、充填効率を算出することが可能となる。上記の内燃機関の制御装置の他の一態様では、前記エンジンは、360度位相がずれた2気筒エンジンである。この態様により、内燃機関の制御装置は、好適に、本発明を実施することができる。以下、図面を参照して本発明の好適な実施の形態について説明する。[内燃機関の概略構成]図1は、本発明に係る内燃機関の制御装置が適用された内燃機関エンジンの概略構成図を示す。図中の実線矢印はガスの流れの一例を示している。内燃機関100は、主に、吸気通路1と、スロットルバルブ2と、エアフロメータ3と、排気通路4と、第1気筒10と、第2気筒20とを備える。吸気通路1には外部から導入された吸気空気が通過し、スロットルバルブ2は吸気通路1を通過する吸気の流量を調整する。スロットルバルブ2は、ECU50から供給される制御信号によって開度が制御される。また、吸気通路1には、第1気筒10及び第2気筒20に供給される吸入空気量を検出するエアフロメータ3が設けられている。エアフロメータ3は、検出した吸入空気量以後、検出吸入空気量Qとも呼ぶ。に相当する検出信号S3を、ECU50へ送信する。吸気通路1を通過した吸気は、第1気筒10の燃焼室17及び第2気筒20の燃焼室27に供給される。また、燃焼室17には、燃料噴射弁インジェクタ15によって噴射された燃料が供給され、燃焼室27には、燃料噴射弁25によって噴射された燃料が供給される。燃料噴射弁15、25は、それぞれ、ECU50から送信される制御信号S15、S25に基づき、燃料噴射量が制御される。更に、燃焼室17には、吸気弁11と排気弁13とが設けられ、燃焼室27には、吸気弁21と排気弁23とが設けられている。吸気弁11は、開閉することによって、吸気通路1と燃焼室17との連通及び遮断を制御し、吸気弁21は、開閉することによって、吸気通路1と燃焼室27との連通及び遮断を制御する。また、吸気弁11には、吸気弁11のバルブ開度を検出するバルブ開度センサ12が設けられ、吸気弁21には、吸気弁21のバルブ開度を検出するバルブ開度センサ22が設けられている。バルブ開度センサ12、22は、検出した吸気弁11、21のバルブ開度の検出信号S12、S22をECU50へ送信する。排気弁13は、開閉することによって、燃焼室17と排気通路4との連通及び遮断を制御し、排気弁23は、開閉することによって、燃焼室27と排気通路4との連通及び遮断を制御する。燃焼室17、27では、吸気行程において上記のように供給された吸気と燃料との混合気が、圧縮行程を経た後、それぞれ点火プラグ14、24によって点火されることにより膨張行程で燃焼される。この場合、燃焼によってピストン18、28が往復運動し、当該往復運動がコンロッド19、29を介してクランク軸不図示に伝達され、クランク軸が回転する。燃焼室17、27での燃焼により発生した排気ガスは、排気行程において排気通路4へ排出される。ここで、第1気筒10と第2気筒20とは、360度位相がずれて動作し、互いの吸気行程がオーバーラップしない。従って、内燃機関100では、間欠吸気を行うことに起因して、吸気脈動が生じる。ECU50は、図示しないCPUCentralProcessingUnit、ROMReadOnlyMemory及びRAMRandomAccessMemoryなどを備え、内燃機関100の各構成要素に対して種々の制御を行う。例えば、ECU50は、上記のようにして供給された検出信号に基づいて、燃料噴射弁15、25等に対する制御を行う。そして、ECU50は、本発明における吸入空気量検出手段、差分算出手段、時間幅算出手段、脈動波形算出手段、及び吸入空気量推定手段として機能する。なお、図1の構成は、一例であり、本発明が適用可能な構成は、必ずしもこれに限定されない。例えば、内燃機関100は、吸気ポートに燃料噴射弁15、25が設置されたポート噴射であったが、これに代えて、またはこれに加えて、燃焼室17、27に直接燃料を噴射する筒内噴射であってもよい。また、内燃機関100は、点火プラグ14、24を備えない自己着火方式であってもよい。[制御方法]次に、ECU50が実行する制御について具体的に説明する。概略的には、ECU50は、吸気脈動が極大及び極小となるタイミングでの検出吸入空気量Qに基づき、吸気脈動の波形を推定し、当該吸気脈動波形に基づき燃焼室17、27内に流入した吸入空気量の推定値推定吸入空気量Qeとも呼ぶ。を算出する。これにより、ECU50は、空気の充填効率を精度よく推定し、空燃比の制御性を向上させる。これについて、図2を参照して具体的に説明する。図2は、吸入空気量、第1気筒10の行程、第2気筒20の行程を示すタイムチャートの一例である。なお、各Q1~Q10は、各時刻t1~t10でエアフロメータ3により検出した検出吸入空気量を示し、グラフB1は、推定される吸入空気量の波形即ち、脈動波形を示す。まず、ECU50は、吸入空気量が極大及び極小となるタイミングで検出吸入空気量Qを計測する。具体的には、ECU50は、吸気弁11、21の開度が最大全開となる状態開弁状態とも呼ぶ。、及び当該開度が最小全閉となる状態閉弁状態とも呼ぶ。となるタイミングで、検出吸入空気量Qを計測する。図2では、ECU50は、吸気弁11が開弁状態となる時刻t2、t6、t10でそれぞれ吸入空気量Q2、Q6、Q10を計測し、吸気弁21が開弁状態となる時刻t4、t8でそれぞれ吸入空気量Q4、Q8を計測する。また、ECU50は、吸気弁11、21が共に閉弁状態となる時刻t1、t3、t5、t7、t9でそれぞれ吸入空気量Q1、Q3、Q5、Q7、Q9を計測する。次に、ECU50は、計測した検出吸入空気量Qから吸気脈動の波形を推定する。ここでは、ECU50は、正弦波SIN波を吸気脈動の波形と推定する。具体的には、ECU50は、吸気弁11、21の開弁状態での検出吸入空気量QをQc、吸気弁11、21の閉弁状態での検出吸入空気量QをQoとすると、これらの差分dQ=Qc-Qoを算出すると共に、吸気行程の時間幅に相当する吸気弁11、21が開弁してから閉弁するまでの時間幅dTを算出する。そして、ECU50は、差分dQと時間幅dTとを用いて、吸気脈動を考慮した吸入空気量の時間tの関数吸入空気量時間関数Qthとも呼ぶ。を、以下の式1に示すように求める。Qth=dQsin2πt/4/dT式1ここで、式1に示す吸入空気量時間関数Qthは、極大値がdQとなる周期4dTの時間関数である。このように、吸気脈動を式1の吸入空気量時間関数Qthにより表すことにより、吸気脈動に起因した吸気通路1内での逆流の影響を除外し、即ち逆流の影響を含めて積算してしまうことによる誤差を抑制し、より正確な吸入空気量を推定することが可能となる。次に、図2に示す具体例を用いて、吸気脈動の波形の推定方法について説明する。図2の時刻t2では、ECU50は、検出吸入空気量Q1、Q2に基づき差分dQ=Q2-Q1を定めると共に、吸気弁11が開弁してから閉弁するまでの時間幅を時間幅dTとする。そして、ECU50は、算出した差分dQと時間幅dTとに基づき、式1により脈動波形に相当する吸入空気量時間関数Qthを求める。時刻t4、時刻t6、時刻t8、時刻t10でも同様に、ECU50は、同時刻で取得した検出吸入空気量Qと、その直前に取得した検出吸入空気量Qとに基づき、差分dQを算出すると共に、吸気弁11が開弁してから閉弁するまでの時間幅を時間幅dTに定める。そして、ECU50は、算出した差分dQ及び時間幅dTに基づき、式1により脈動波形に相当する吸入空気量時間関数Qthを求める。そして、ECU50は、式1を吸気弁11、21の各々の開弁期間で積分することにより、第1気筒10及び第2気筒20のそれぞれに流入した空気量の推定値である推定吸入空気量Qeを算出し、当該推定吸入空気量Qeに基づき第1気筒10及び第2気筒20の充填効率を求める。具体的には、ECU50は、吸気弁11の開弁期間で式1を時間積分することで、第1気筒10の燃焼室17に流入した推定吸入空気量Qeを算出する。そして、ECU50は、当該推定吸入空気量Qeに基づき、予め作成された式又はマップを参照して、充填効率を算出する。同様に、ECU50は、吸気弁21の開弁期間で式1を時間積分することで、第2気筒20の燃焼室27に流入した推定吸入空気量Qeを算出する。そして、ECU50は、当該推定吸入空気量Qeに基づき、予め作成された式又はマップを参照して、充填効率を算出する。そして、ECU50は、吸気脈動を考慮して算出した上述の充填効率と、目標空燃比とに基づき、例えば予め作成された式又はマップを参照して、燃料噴射弁15、25の燃料噴射量を決定し、制御信号S15、S25により燃料噴射弁15、25の制御を行う。このようにすることで、ECU50は、吸気脈動が生じる場合であっても、吸気脈動を考慮したより的確な充填効率を算出し、空燃比を目標空燃比に精度よく合わせることができる。従って、ECU50は、吸気脈動による空燃比の目標値と実値とのずれに起因した燃費悪化、エミッション悪化などを抑制することができる。処理フロー図3は、本実施形態における処理手順を示すフローチャートである。図3に示すフローチャートは、ECU50により所定の周期に従い繰り返し実行される。まず、ECU50は、第1気筒10の吸気弁11又は第2気筒20の吸気弁21が開弁状態又は閉弁状態であるか否か判定するステップS101。具体的には、ECU50は、各吸気弁11、21が開弁状態か閉弁状態かをバルブ開度センサ12、22の検出信号S12、S22に基づき特定する。そして、ECU50は、吸気弁11又は吸気弁21が開弁状態又は閉弁状態ではないと判断した場合ステップS101;No、フローチャートの処理を終了する。次に、ECU50は、吸気弁11又は吸気弁21が、開弁状態又は閉弁状態であると判断した場合ステップS101;Yes、検出吸入空気量Qを計測するステップS102。具体的には、ECU50は、吸気弁11又は吸気弁22が開弁状態又は閉弁状態の場合に、エアフロメータ3から受信した検出信号S3に基づき、当該検出吸入空気量Qを特定する。そして、ECU50は、吸気脈動波形に相当する吸入空気量時間関数Qthを算出するステップS103。具体的には、ECU50は、上述した差分dQ及び時間幅dTに基づき、式1を用いて、吸入空気量時間関数Qthを算出する。次に、ECU50は、吸入空気量時間関数Qthを積分し、推定吸入空気量Qeを算出し、さらに当該推定吸入空気量Qeに基づき充填効率を算出するステップS104。具体的には、ECU50は、吸入空気量時間関数Qthを、吸気弁11、22が開弁している即ち、閉弁状態ではない期間で時間積分することにより、当該期間で燃焼室17、27に流入した推定吸入空気量Qeを推定し、当該推定吸入空気量Qeから充填効率を求める。そして、ECU50は、算出した充填効率と、目標空燃比とに基づき、燃料噴射弁15、25で噴射すべき燃料噴射量を算出し、燃料噴射弁15、25の燃料噴射制御を行うステップS105。[変形例]図1等の説明では、内燃機関100は、360度位相の異なる2気筒を有していた。しかし、本発明が適用可能な構成は、これに限定されない。これに代えて、内燃機関100は、2気筒以外の気筒数を有する構成、または/及び、各気筒が360度以外の所定度だけ異なる位相を有する構成であってもよい。この場合であっても、内燃機関100は、各気筒の吸気行程がオーバーラップせず、間欠吸気を行う構成を有する。そして、ECU50は、各気筒の吸気弁が開弁状態又は閉弁状態となるタイミングでの検出吸入空気量Qに基づき、式1に相当する吸気脈動波形を推定し、当該吸気脈動波形を吸気弁の開弁期間で積分することで、推定吸入空気量Qeを算出する。これによっても、ECU50は、吸気脈動を考慮して充填効率を高精度に算出することが可能となる。