(57)【要約】【課題】内燃機関の制御装置において、複数の制御入力に基づいて内燃機関の各系の状態量を演算することができる物理モデルを搭載することである。【解決手段】内燃機関制御装置のエンジン制御ユニット４０は、制御入力４２をエンジンシステムの物理モデルに適用し、各要素における状態量を算出する状態量算出部６０を有する。状態量算出部６０は、内燃機関のサイクリックな現象に基づき、各気筒の吸気−燃焼−排気の１サイクルを複数の離散的時間に分割し、エンジンシステム全体についての離散的時間の物理モデルである統合離散的時間モデル４４を作成する。統合離散的時間モデル４４は、その下位階層物理モデルとして、吸気系離散的時間モデル７０、燃焼系離散的時間モデル７２、排気系離散的時間モデル７４、燃料系離散的時間モデル７６、並進回転系離散的時間モデル７８、エンジン性能モデル８０を有する。【選択図】図２
2008175198.tif 000002
【特許請求の範囲】
【請求項１】
複数の制御入力を取得する手段と、
取得された制御入力を内燃機関の物理モデルに適用して、内燃機関の各要素における状態量を算出する状態量算出手段と、
算出された状態量に基づき、内燃機関の性能を算出する手段と、
を備え、
状態量算出手段は、
内燃機関の気筒数に応じて吸気−燃焼−排気の１サイクルを複数の離散的時間に分割する手段と、
複数の状態量について、内燃機関全体についての離散的時間の物理モデルである統合離散的時間モデルを作成する手段と、
統合離散的時間モデルに、取得された制御入力を適用し、各離散的時間における内燃機関全体の複数の状態量を逐時的に算出する手段と、
を有することを特徴とする内燃機関の性能制御装置。
【請求項２】
複数の制御入力を取得する手段と、
取得された制御入力を内燃機関の物理モデルに適用して、内燃機関の各要素における状態量を算出する状態量算出手段と、
算出された状態量に基づき、内燃機関の性能を算出する手段と、
を備え、
状態量算出手段は、
内燃機関の気筒数に応じて吸気−燃焼−排気の１サイクルを複数の離散的時間に分割する手段と、
状態量を質量と圧力として、内燃機関全体についての離散的時間の物理モデルである統合離散的時間モデルを作成する手段と、
統合離散的時間モデルに、取得された制御入力を適用し、各離散的時間における内燃機関全体の質量と圧力を逐時的に算出する手段と、
を有することを特徴とする内燃機関の性能制御装置。
【請求項３】
請求項１または２に記載の内燃機関の性能制御装置において、
統合離散的時間モデルを作成する手段は、
内燃機関の各気筒ごとの状態量である気筒毎状態量に関する離散的時間の物理モデルを作成することを特徴とする内燃機関の性能制御装置。
【請求項４】
請求項３に記載の内燃機関の性能制御装置において、
統合離散的時間モデルを作成する手段は、
気筒数と同じ数の離散的時間の物理モデルを作成することを特徴とする内燃機関の性能制御装置。
【請求項５】
請求項１または２に記載の内燃機関の性能制御装置において、
統合離散的時間モデルを作成する手段は、
内燃機関の吸気−燃焼−排気のサイクリックな複数の遷移過程ごとの状態量である遷移過程毎状態量に関する離散的時間の物理モデルを作成することを特徴とする内燃機関の性能制御装置。
【請求項６】
請求項５に記載の内燃機関の性能制御装置において、
統合離散的時間モデルを作成する手段によって作成された離散的時間の物理モデルの数は、気筒数に依存しないことを特徴とする内燃機関の性能制御装置。
【請求項７】
請求項１または２に記載の内燃機関の性能制御装置において、
統合離散的時間モデルを作成する手段は、離散的時間の物理モデルを作成する方法の選択を行う選択手段を有し、
選択手段は、選択の指示に応じて、
内燃機関の各気筒ごとの状態量である気筒毎状態量に関する離散的時間の物理モデルを作成する方法であって、気筒数と同じ数の離散的時間の物理モデルを作成する方法と、
内燃機関の吸気−燃焼−排気のサイクリックな複数の遷移過程ごとの状態量である遷移過程毎状態量に関する離散的時間の物理モデルを作成する方法であって、気筒数に依存しない数の離散的時間の物理モデルを作成する方法との間で選択を行うことを特徴とする内燃機関の性能制御装置。
【請求項８】
請求項１から７のいずれか１に記載の内燃機関の性能制御装置において、
統合離散的時間モデルを作成する手段は、統合離散的時間モデルの下位階層の物理モデルとして、少なくとも、
吸気行程についての離散的時間の物理モデルである吸気系離散的時間モデルを作成し、
排気行程についての離散的時間の物理モデルである排気系離散的時間モデルを作成し、
燃焼行程についての離散的時間の物理モデルである燃焼系離散的時間モデルを作成することを特徴とする内燃機関の性能制御装置。
【請求項９】
請求項１から８のいずれか１に記載の内燃機関の性能制御装置において、
離散的時間に分割する手段は、
複数の気筒が同時に吸気行程を行う期間と、１つの気筒のみが吸気行程を行う期間とを区別して、吸気行程を複数の離散的時間に分割することを特徴とする内燃機関の性能制御装置。
【請求項１０】
請求項１から８のいずれか１に記載の内燃機関の性能制御装置において、
離散的時間に分割する手段は、
複数の気筒が同時に排気行程を行う期間と、１つの気筒のみが排気行程を行う期間とを区別して、排気行程を複数の離散的時間に分割することを特徴とする内燃機関の性能制御装置。
【請求項１１】
請求項８または請求項９に記載の内燃機関の性能制御装置において、
吸気系離散的時間モデルの作成は、
サージタンクを含む吸気系の全体積と、吸気バルブが開いたときの気筒の体積とに基づいて、吸気行程の開始の離散的時間における全質量と全エネルギを計算し、
サージタンクを含む吸気系の全体積と、吸気バルブが閉じたときの気筒の体積とに基づいて、吸気行程の終了の離散的時間における全質量と全エネルギを計算し、
吸気行程の開始と終了との間を複数の離散的時間に分割する場合に、サージタンクを含む吸気系の全体積と、分割された途中の離散的時間における気筒の体積とに基づいて、途中の離散的時間における全質量と全エネルギを計算し、
吸気行程の開始から終了までの間について、各離散的時間ごとに全質量の収支と全エネルギの収支とをそれぞれ計算することで、状態量についての物理モデルを作成することを特徴とする内燃機関の性能制御装置。
【請求項１２】
請求項８または請求項１０に記載の内燃機関の性能制御装置において、
排気系離散的時間モデルの作成は、
排気管を含む排気系の全体積と、排気バルブが開いたときの気筒の体積とに基づいて、排気行程の開始の離散的時間における全質量と全エネルギを計算し、
排気管を含む排気系の全体積と、排気バルブが閉じたときの気筒の体積とに基づいて、排気行程の終了の離散的時間における全質量と全エネルギを計算し、
排気行程の開始と終了との間を複数の離散的時間に分割する場合に、排気管を含む排気系の全体積と、分割された途中の離散的時間における気筒の体積とに基づいて、途中の離散的時間における全質量と全エネルギを計算し、
排気行程の開始から終了までの間について、各離散的時間ごとに全質量の収支と全エネルギの収支とをそれぞれ計算することで、状態量についての物理モデルを作成することを特徴とする内燃機関の性能制御装置。
【請求項１３】
請求項８に記載の内燃機関の性能制御装置において、
燃焼系離散的時間モデルの作成は、
吸気弁が閉じたときから排気弁が開くときまでについて、質量保存則とエネルギ保存則とに基づいて、近似計算を利用して微分方程式を解析的に解いて物理モデルを作成することを特徴とする内燃機関の性能制御装置。
【請求項１４】
請求項１から７のいずれか１に記載の内燃機関の性能制御装置において、
複数の離散的時間に分割する手段は、
複数の気筒が同時に吸気行程を行う期間と、１つの気筒のみが吸気行程を行う期間とを区別して、吸気行程を複数の離散的時間に分割し、さらに、分割された離散的時間を適当に結合して、状態量算出のサンプル間隔を内燃機関全体の１サイクルを気筒数で除した間隔に近い間隔とし、
複数の気筒が同時に排気行程を行う期間と、１つの気筒のみが排気行程を行う期間とを区別して、排気行程を複数の離散的時間に分割し、さらに、分割された離散的時間を適当に結合して、状態量算出のサンプル間隔を内燃機関全体の１サイクルを気筒数で除した間隔に近い間隔とすることで、
状態量算出のサンプル間隔が、内燃機関全体の１サイクルを気筒数で除したものより長短があることを特徴とする内燃機関の性能制御装置。
【請求項１５】
請求項４に記載の内燃機関の性能制御装置において、
統合離散的時間モデルを作成する手段は、
遷移過程毎の制御入力と、気筒毎の制御入力との間を変換する制御入力変換手段と、
遷移過程毎状態量と、気筒毎状態量との間を変換する状態量変換手段と、
を有することを特徴とする内燃機関の性能制御装置。
【請求項１６】
請求項１５に記載の内燃機関の性能制御装置において、
制御入力変換手段は、
制御入力が入力される遷移過程と、同じ制御入力が入力される気筒との関連付けに基づく演算作用素を用いて、遷移過程毎の制御入力と、気筒毎の制御入力との間を変換することを特徴とする内燃機関の性能制御装置。
【請求項１７】
請求項１５に記載の内燃機関の性能制御装置において、
状態量変換手段は、
任意の遷移過程と、その遷移過程が生じる気筒との関連付けに基づく変換行列を用いて、遷移過程毎の状態量と、気筒毎の状態量との間を変換することを特徴とする内燃機関の性能制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、内燃機関の性能制御装置に係り、特に、複数の制御入力を用いて内燃機関の複数の性能を制御する内燃機関の性能制御装置に関する。
【背景技術】
【０００２】
内燃機関の性能を制御するためには、内燃機関の物理モデルを用いて制御入力に対する性能応答を演算によって出力し、これを用いることが行われる。
【０００３】
例えば、特許文献１には、運転者の要求するトルクに対応して筒内充填空気量、燃料噴射量、点火時期等を制御するトルクディマンド制御において、スロットル開度の変化と実際の筒内充填空気量の変化に時間遅れが生じることに関し、これをフィードバック制御にフィードフォワード制御を加えて、システムの特性ばらつきや経時変化、運転条件や環境条件の変化による動特性の変化に対応する方法について述べられている。ここでは、エンジン運転条件であるスロットル通過空気量と吸気圧力に基づいて実際に筒内に充填される空気量（実筒内充填空気量）を推定し、これを、スロットル開度を入力とし、筒内充填空気量を出力として、むだ時間と遅れ要素の線形数学モデルで近似的に表現している。そして、目標筒内充填空気量の規範値に基づき、フィードフォワード的に目標スロットル開度を演算するとともに、この目標筒内充填空気量の規範値と実筒内充填空気量の推定値との偏差に基づいて、目標スロットル開度に対するフィードバック補正量を演算して補正し、最終目標スロットル開度を求めることが述べられている。
【０００４】
【特許文献１】特開２００２−３０９９９０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
特許文献１においては、スロットル開度を入力として、むだ時間と遅れ要素の線形数学モデルを用い、筒内充填空気量を演算によって出力することが行われる。内燃機関の挙動は、複数の制御入力が、吸気系、排気系、燃焼系等の各系における複数の状態量、あるいは複数の性能出力に関係するので、その挙動の全体を物理モデルで記述し、これをエンジン制御コンピュータに搭載して制御演算を実行させるのは大変である。そこで、特許文献１の１入力１出力のモデルのように、個別に解決したい性能出力に対する簡単な物理モデルあるいは数学モデルを構築し、これらをいくつか用いて、内燃機関の性能制御を実行することが行われている。
【０００６】
このように、従来技術では、内燃機関の性能制御において、複数の制御入力に応じた各系の複数の状態量の応答及びエンジン性能の応答等を統合的に演算することができる物理モデルをエンジン制御コンピュータに搭載することが困難である。その結果、全制御要求を満足する複数の制御入力の最適解を、相互関係や操作量の制約等の拘束を考慮しながら、すべてを同時に決定することが困難となっている。
【０００７】
本発明の目的は、複数の制御入力に基づいて内燃機関の各系の状態量を演算することを可能とする内燃機関の制御装置を提供することである。また、他の目的は、複数の制御入力に基づいて内燃機関の各系の状態量を演算することができる物理モデルを搭載する内燃機関の制御装置を提供することである。また、他の目的は、全制御要求を満足する複数の制御入力の最適解を、相互関係や操作量の制約等の拘束を考慮しながら、すべてを同時に決定することを可能にする内燃機関の制御装置を提供することである。また、他の目的は、演算量を低減して、複数の制御入力に基づいて内燃機関の各系の状態量を演算することを可能とする内燃機関の制御装置を提供することである。以下の手段は、上記目的の少なくとも１つに貢献する。
【課題を解決するための手段】
【０００８】
本発明に係る内燃機関の性能制御装置は、複数の制御入力を取得する手段と、取得された制御入力を内燃機関の物理モデルに適用して、内燃機関の各要素における状態量を算出する状態量算出手段と、算出された状態量に基づき、内燃機関の性能を算出する手段と、を備え、状態量算出手段は、内燃機関の気筒数に応じて吸気−燃焼−排気の１サイクルを複数の離散的時間に分割する手段と、複数の状態量について、内燃機関全体についての離散的時間の物理モデルである統合離散的時間モデルを作成する手段と、統合離散的時間モデルに、取得された制御入力を適用し、各離散的時間における内燃機関全体の複数の状態量を逐時的に算出する手段と、を有することを特徴とする。
【０００９】
本発明に係る内燃機関の性能制御装置は、複数の制御入力を取得する手段と、取得された制御入力を内燃機関の物理モデルに適用して、内燃機関の各要素における状態量を算出する状態量算出手段と、算出された状態量に基づき、内燃機関の性能を算出する手段と、を備え、状態量算出手段は、内燃機関の気筒数に応じて吸気−燃焼−排気の１サイクルを複数の離散的時間に分割する手段と、状態量を質量と圧力として、内燃機関全体についての離散的時間の物理モデルである統合離散的時間モデルを作成する手段と、統合離散的時間モデルに、取得された制御入力を適用し、各離散的時間における内燃機関全体の質量と圧力を逐時的に算出する手段と、を有することを特徴とする。
【００１０】
また、本発明に係る内燃機関の性能制御装置において、統合離散的時間モデルを作成する手段は、内燃機関の各気筒ごとの状態量である気筒毎状態量に関する離散的時間の物理モデルを作成することが好ましい。
【００１１】
また、本発明に係る内燃機関の性能制御装置において、統合離散的時間モデルを作成する手段は、気筒数と同じ数の離散的時間の物理モデルを作成することを特徴とする。
【００１２】
また、本発明に係る内燃機関の性能制御装置において、統合離散的時間モデルを作成する手段は、内燃機関の吸気−燃焼−排気のサイクリックな複数の遷移過程ごとの状態量である遷移過程毎状態量に関する離散的時間の物理モデルを作成することが好ましい。
【００１３】
また、本発明に係る内燃機関の性能制御装置において、統合離散的時間モデルを作成する手段によって作成された離散的時間の物理モデルの数は、気筒数に依存しないことを特徴とする。
【００１４】
また、本発明に係る内燃機関の性能制御装置において、統合離散的時間モデルを作成する手段は、離散的時間の物理モデルを作成する方法の選択を行う選択手段を有し、選択手段は、選択の指示に応じて、内燃機関の各気筒ごとの状態量である気筒毎状態量に関する離散的時間の物理モデルを作成する方法であって、気筒数と同じ数の離散的時間の物理モデルを作成する方法と、内燃機関の吸気−燃焼−排気のサイクリックな複数の遷移過程ごとの状態量である遷移過程毎状態量に関する離散的時間の物理モデルを作成する方法であって、気筒数に依存しない数の離散的時間の物理モデルを作成する方法との間で選択を行うことが好ましい。
【００１５】
また、本発明に係る内燃機関の性能制御装置において、統合離散的時間モデルを作成する手段は、統合離散的時間モデルの下位階層の物理モデルとして、少なくとも、吸気行程についての離散的時間の物理モデルである吸気系離散的時間モデルを作成し、排気行程についての離散的時間の物理モデルである排気系離散的時間モデルを作成し、燃焼行程についての離散的時間の物理モデルである燃焼系離散的時間モデルを作成することが好ましい。
【００１６】
また、本発明に係る内燃機関の性能制御装置において、離散的時間に分割する手段は、複数の気筒が同時に吸気行程を行う期間と、１つの気筒のみが吸気行程を行う期間とを区別して、吸気行程を複数の離散的時間に分割することが好ましい。
【００１７】
また、本発明に係る内燃機関の性能制御装置において、離散的時間に分割する手段は、複数の気筒が同時に排気行程を行う期間と、１つの気筒のみが排気行程を行う期間とを区別して、排気行程を複数の離散的時間に分割することが好ましい。
【００１８】
また、本発明に係る内燃機関の性能制御装置において、吸気系離散的時間モデルの作成は、サージタンクを含む吸気系の全体積と、吸気バルブが開いたときの気筒の体積とに基づいて、吸気行程の開始の離散的時間における全質量と全エネルギを計算し、サージタンクを含む吸気系の全体積と、吸気バルブが閉じたときの気筒の体積とに基づいて、吸気行程の終了の離散的時間における全質量と全エネルギを計算し、吸気行程の開始と終了との間を複数の離散的時間に分割する場合に、サージタンクを含む吸気系の全体積と、分割された途中の離散的時間における気筒の体積とに基づいて、途中の離散的時間における全質量と全エネルギを計算し、吸気行程の開始から終了までの間について、各離散的時間ごとに全質量の収支と全エネルギの収支とをそれぞれ計算することで、状態量についての物理モデルを作成することが好ましい。
【００１９】
また、本発明に係る内燃機関の性能制御装置において、排気系離散的時間モデルの作成は、排気管を含む排気系の全体積と、排気バルブが開いたときの気筒の体積とに基づいて、排気行程の開始の離散的時間における全質量と全エネルギを計算し、排気管を含む排気系の全体積と、排気バルブが閉じたときの気筒の体積とに基づいて、排気行程の終了の離散的時間における全質量と全エネルギを計算し、排気行程の開始と終了との間を複数の離散的時間に分割する場合に、排気管を含む排気系の全体積と、分割された途中の離散的時間における気筒の体積とに基づいて、途中の離散的時間における全質量と全エネルギを計算し、排気行程の開始から終了までの間について、各離散的時間ごとに全質量の収支と全エネルギの収支とをそれぞれ計算することで、状態量についての物理モデルを作成することが好ましい。
【００２０】
また、本発明に係る内燃機関の性能制御装置において、燃焼系離散的時間モデルの作成は、吸気弁が閉じたときから排気弁が開くときまでについて、質量保存則とエネルギ保存則とに基づいて、近似計算を利用して微分方程式を解析的に解いて物理モデルを作成することが好ましい。
【００２１】
また、本発明に係る内燃機関の性能制御装置において、複数の離散的時間に分割する手段は、複数の気筒が同時に吸気行程を行う期間と、１つの気筒のみが吸気行程を行う期間とを区別して、吸気行程を複数の離散的時間に分割し、さらに、分割された離散的時間を適当に結合して、状態量算出のサンプル間隔を内燃機関全体の１サイクル７２０度を気筒数で除した間隔に近い間隔とし、複数の気筒が同時に排気行程を行う期間と、１つの気筒のみが排気行程を行う期間とを区別して、排気行程を複数の離散的時間に分割し、さらに、分割された離散的時間を適当に結合して、状態量算出のサンプル間隔を内燃機関全体の１サイクル７２０度を気筒数で除した間隔に近い間隔とすることで、状態量算出のサンプル間隔が、内燃機関全体の１サイクル７２０度を気筒数で除したものより長短があってもよい。
【００２２】
また、本発明に係る内燃機関の性能制御装置において、統合離散的時間モデルを作成する手段は、遷移過程毎の制御入力と、気筒毎の制御入力との間を変換する制御入力変換手段と、遷移過程毎状態量と、気筒毎状態量との間を変換する状態量変換手段と、を有することが好ましい。
【００２３】
また、本発明に係る内燃機関の性能制御装置において、制御入力変換手段は、制御入力が入力される遷移過程と、同じ制御入力が入力される気筒との関連付けに基づく演算作用素を用いて、遷移過程毎の制御入力と、気筒毎の制御入力との間を変換することが好ましい。
【００２４】
また、本発明に係る内燃機関の性能制御装置において、状態量変換手段は、任意の遷移過程と、その遷移過程が生じる気筒との関連付けに基づく変換行列を用いて、遷移過程毎の状態量と、気筒毎の状態量との間を変換することが好ましい。
【発明の効果】
【００２５】
上記構成により、内燃機関の性能制御装置は、内燃機関の気筒数に応じて吸気−燃焼−排気の１サイクルを複数の離散的時間に分割し、複数の状態量について、内燃機関全体についての離散的時間の物理モデルである統合離散的時間モデルを作成する。そして、統合離散的時間モデルに、複数の制御入力を適用し、各離散的時間における内燃機関全体の複数の状態量を逐時的に算出する。ここで、質量と圧力とを、複数の状態量とすることができる。また、質量と圧力以外に、温度等を状態量の組合せに加えてもよい。
【００２６】
ここで、内燃機関は、一般的に気筒数が定まると、各気筒の吸気−燃焼−排気の１サイクルは、一定の繰り返しを行い、これの総合が内燃機関全体としてのサイクルとなる。したがって、内燃機関のサイクリックな現象に基づいて、内燃機関のサイクルを複数の離散的時間に分割して物理モデルを構成することが可能で、その離散的時間の分割は、気筒数に依存し、内燃機関の大きさ等のパラメータに依存しない。また、気筒数が異なっても、そのサイクルの分割の基本的方法は変わらない。つまり、気筒数や具体的なエンジンパラメータが変化しても、離散的時間に分割する定義自体に変更がなく、物理モデル自体の構成は変わらず、例えば、物理モデルを構成する行列の要素の変更ですむ。したがって、この内燃機関のサイクリックな現象に基づいて離散的時間に分割した物理モデルを、内燃機関全体の統合離散的時間モデルとできる。このようにして、記述量を少なくして統合離散的物理モデルを構成することができる。また、内燃機関の基本構成に他のシステムが付加されても、内燃機関のサイクリックな現象は変わらないので、統合離散的時間モデルに対する変更を要しない。この統合離散的時間モデルに、複数の制御入力を適用し、内燃機関全体としての状態量として、質量と圧力を逐時的に算出する。このようにして、複数の制御入力に基づいて内燃機関の各系の状態量を演算することが可能となる。
【００２７】
また、統合離散的時間モデルの下位階層の物理モデルとして、少なくとも、吸気行程についての離散的時間の物理モデルである吸気系離散的時間モデルを作成し、排気行程についての離散的時間の物理モデルである排気系離散的時間モデルを作成し、燃焼行程についての離散的時間の物理モデルである燃焼系離散的時間モデルを作成する。これにより、内燃機関の吸気−燃焼−排気の１サイクルを構成する各系の詳細な質量と圧力を逐次的に算出でき、これを統合して、内燃機関全体の状態量の逐時的算出を裏付けることができる。
【００２８】
また、内燃機関の性能制御装置において、内燃機関の各気筒ごとの状態量である気筒毎状態量に関する離散的時間の物理モデルを作成する。これによって、各気筒ごとに、離散的時間における状態量を把握できる。
【００２９】
また、内燃機関の性能制御装置において、内燃機関の吸気−燃焼−排気のサイクリックな複数の遷移過程ごとの状態量である遷移過程毎状態量に関する離散的時間の物理モデルを作成する。サイクリックな複数気筒の内燃機関においては、どの気筒も同じ遷移過程を異なるタイミングで実行する。したがって、気筒毎状態量を用いる場合に比べ、演算量を少なくできる。
【００３０】
また、内燃機関の性能制御装置において、気筒毎状態量に関する離散的時間の物理モデルを作成するときは、気筒数と同じ数の離散的時間の物理モデルを作成する。一方、遷移過程毎状態量に関する離散的時間の物理モデルを作成するときは、作成された離散的時間の物理モデルの数は、気筒数に依存しない。このように、遷移過程毎状態量に関する離散的時間の物理モデルを作成するときには、作成される物理モデルの数を抑制することができる。
【００３１】
また、内燃機関の性能制御装置において、統合離散的時間モデルを作成する手段は、離散的時間の物理モデルを作成する方法の選択を行う選択手段を有する。選択手段は、選択の指示に応じて、内燃機関の各気筒ごとの状態量である気筒毎状態量に関する離散的時間の物理モデルを作成する方法であって、気筒数と同じ数の離散的時間の物理モデルを作成する方法と、内燃機関の吸気−燃焼−排気のサイクリックな複数の遷移過程ごとの状態量である遷移過程毎状態量に関する離散的時間の物理モデルを作成する方法であって、気筒数に依存しない数の離散的時間の物理モデルを作成する方法との間で選択を行う。これにより、作成する物理モデル数を抑制したいときには、遷移過程毎状態量に関する離散的時間の物理モデルを作成する方法を選択すること等ができる。
【００３２】
また、内燃機関の性能制御装置において、複数の気筒が同時に吸気行程を行う期間と、１つの気筒のみが吸気行程を行う期間とを区別して、吸気行程を複数の離散的時間に分割する。これにより、吸気行程における異なる内容に応じて離散的時間の分割が行え、実際の内燃機関の動作に適合した物理モデルを作成することができる。排気行程についても同様である。
【００３３】
また、内燃機関の性能制御装置において、サージタンクを含む吸気系の全体積と、吸気バルブが開いたときの気筒の体積とに基づいて、吸気行程の開始の離散的時間における全質量と全エネルギを計算し、サージタンクを含む吸気系の全体積と、吸気バルブが閉じたときの気筒の体積とに基づいて、吸気行程の終了の離散的時間における全質量と全エネルギを計算し、吸気行程の開始と終了との間を複数の離散的時間に分割する場合に、サージタンクを含む吸気系の全体積と、分割された途中の離散的時間における気筒の体積とに基づいて、途中の離散的時間における全質量と全エネルギを計算し、吸気行程の開始から終了までの間について、各離散的時間ごとに全質量の収支と全エネルギの収支とをそれぞれ計算することで、状態量についての物理モデルを作成する。体積の変化に着目することで、質量保存則とエネルギ保存則を適用して状態量の変化を求めることができる。これにより、複雑な微分方程式を解くことを要せずに、吸気系の物理モデルを作成することができる。
【００３４】
同様に、排気管を含む排気系の全体積と、排気バルブが開いたときの気筒の体積とに基づいて、排気行程の開始の離散的時間における全質量と全エネルギを計算し、排気管を含む排気系の全体積と、排気バルブが閉じたときの気筒の体積とに基づいて、排気行程の終了の離散的時間における全質量と全エネルギを計算し、排気行程の開始と終了との間を複数の離散的時間に分割する場合に、排気管を含む排気系の全体積と、分割された途中の離散的時間における気筒の体積とに基づいて、途中の離散的時間における全質量と全エネルギを計算し、排気行程の開始から終了までの間について、各離散的時間ごとに全質量の収支と全エネルギの収支とをそれぞれ計算し、状態量についての物理モデルを作成することで、複雑な微分方程式を解くことを要せずに、排気系の物理モデルを作成することができる。
【００３５】
また、複数の気筒が同時に吸気行程を行う期間と、１つの気筒のみが吸気行程を行う期間とを区別して、吸気行程を複数の離散的時間に分割し、さらに、分割された離散的時間を適当に結合して、状態量算出のサンプル間隔を内燃機関全体の１サイクル７２０度を気筒数で除した間隔に近い間隔とし、複数の気筒が同時に排気行程を行う期間と、１つの気筒のみが排気行程を行う期間とを区別して、排気行程を複数の離散的時間に分割し、さらに、分割された離散的時間を適当に結合して、状態量算出のサンプル間隔を内燃機関全体の１サイクル７２０度を気筒数で除した間隔に近い間隔とすることで、状態量算出のサンプル間隔が、内燃機関全体の１サイクル７２０度を気筒数で除したものより長短がある。したがって、単純に１サイクル７２０度を内燃機関の気筒数で除した離散的時間の分割でなく、実際の内燃機関の現象に適合した離散的時間の分割を用いることができる。
【００３６】
また、内燃機関の性能制御装置において遷移過程毎の制御入力と、気筒毎の制御入力との間を変換する手段と、遷移過程毎状態量と、気筒毎状態量との間を変換する手段を有する。これによって、気筒毎の記述と制御過程毎の記述との間の変換を容易に行うことができる。
【００３７】
また、制御入力が入力される遷移過程と、同じ制御入力が入力される気筒との関連付けに基づく演算作用素を用いて、遷移過程毎の制御入力と、気筒毎の制御入力との間を変換する。これよって、関連付け表す演算作用素を用いた演算によって、容易に制御入力の記述の変換を行うことができる。
【００３８】
また、任意の遷移過程と、その遷移過程が生じる気筒との関連付けに基づく変換行列を用いて、遷移過程毎の状態量と、気筒毎の状態量との間を変換する。これによって、関連付けを表す行列を用いた演算によって、容易に状態量の記述の変換を行うことができる。
【００３９】
上記のように、本発明に係る内燃機関の制御装置によれば、複数の制御入力に基づいて内燃機関の各系の状態量を演算することが可能となる。また、本発明に係る内燃機関の制御装置によれば、記述量が少ない離散的時間の物理モデルを用いるので、複数の制御入力に基づいて内燃機関の各系の状態量を演算することができる物理モデルを搭載することができる。また、本発明に係る内燃機関の制御装置によれば、遷移過程毎の状態量を用いることで、演算量を低減して、複数の制御入力に基づいて内燃機関の各系の状態量を演算することが可能となる。
【発明を実施するための最良の形態】
【００４０】
以下に図面を用いて、本発明に係る実施の形態につき、詳細に説明する。以下では、６気筒の内燃機関のための性能制御装置として説明するが、これは代表的説明としての１例であって、他の気筒数であってもよい。また、離散的時間の分割のサンプル点の設定も説明のための１例であって、内燃機関のサイクリックな現象に適合する他のサンプル点の設定を行うこともできる。
【実施例１】
【００４１】
図１は、内燃機関の制御装置３０と、その制御対象の内燃機関である６気筒エンジンシステム１０とを示す図である。制御対象の６気筒エンジンシステム１０は、６つの気筒１２と、吸気側のエアフィルタ１４と、吸気量を制御するためのスロットル１６と、サージタンク１８と、サージタンク１８から各気筒１２に分岐する吸気管２０と、各気筒１２ごとに設けられる燃料噴射弁２１及び点火栓２２と、クランク回転角度を検出する回転角度センサ２４と、各気筒１２ごとに設けられ、エキゾーストマニホールド集合部を含む排気管（以下、これを単にエキゾーストマニホールド集合部２６ということにする）と、触媒２８とを含んで構成される。
【００４２】
内燃機関の制御装置３０は、目標回転数や目標トルク等の目標指令入力部３２と、エンジンシステム１０から伝送されるクランク回転角度等を受け取る観測量取得部３４と、エンジン制御ユニット（ＥＣＵ）４０と、ＥＣＵ４０からエンジンシステム１０にスロットル開度等の操作量すなわち制御入力を出力する制御出力部３６とを含んで構成される。
【００４３】
ＥＣＵ４０は、目標指令入力部３２から受け取った目標指令と、エンジンシステム１０から取得した観測量とに基づき、エンジンシステム１０に操作量すなわち制御入力を出力する機能を有する。具体的には、エンジンシステム１０に対する全体的な物理モデルである統合離散的時間モデルを有し、この統合離散的時間モデルに基づいた状態量算出部６０によってエンジン性能を算出し、あるいは、逆に、状態量算出部６０にエンジン性能４６を入力して与え、これを満たす制御入力４２を算出する機能を有する。かかるＥＣＵ４０は、適当なコンピュータで構成することができる。
【００４４】
図２は、ＥＣＵ４０の詳細な構成を説明するブロック図である。ＥＣＵ４０は、制御対象のエンジンシステムの気筒数等を含む具体的なエンジンパラメータ５０を取得するエンジンパラメータ取得部と、スロットル開度θt、燃料噴射量Ｆi、点火時期θS等の制御入力４２を取得する制御入力取得部と、制御入力４２をエンジンシステムの物理モデルに適用してエンジンシステムの各要素における状態量を算出する状態量算出部６０と、算出された状態量と制御入力４２とに基づき、エンジンシステムのトルク、燃料消費量等のエンジン性能を算出する性能算出部４６とを含んで構成される。
【００４５】
状態量算出部６０は、制御対象のエンジンシステムに対する離散的時間の物理モデルを作成し記憶する物理モデル作成記憶部６２と、質量・圧力・回転角速度算出部８２とを含んで構成される。ここで、質量・圧力・回転角速度算出部８２は、物理モデル作成記憶部６２で作成された離散的時間モデルに複数の制御入力４２を適用して、エンジンシステムの全体及び各要素における質量と圧力とを逐時的に算出し、また、算出された質量と圧力の逐時的変化をピストン機構についての並進回転系離散的時間モデルに適用してエンジンシステムの回転角速度を逐時的に算出する機能を有する。
【００４６】
物理モデル作成記憶部６２は、内燃機関のサイクリックな現象に基づき、各気筒の吸気−燃焼−排気の１サイクルを複数の離散的時間に分割する離散的時間分割部６６と、状態量を質量と圧力として、エンジンシステム全体についての離散的時間の物理モデルを作成する物理モデル作成部６８とを含んで構成される。
【００４７】
ここで、エンジンシステム全体についての離散的時間の物理モデルを、統合離散的時間モデル４４と呼ぶことにすると、統合離散的時間モデル４４は、その下位階層の物理モデルとして、吸気行程についての離散的時間の物理モデルである吸気系離散的時間モデル７０、燃焼行程についての離散的時間の物理モデルである燃焼系離散的時間モデル７２、排気行程についての離散的時間の物理モデルである排気系離散的時間モデル７４、燃料系についての離散的時間の物理モデルである燃料系離散的時間モデル７６、並進回転系についての離散的時間の物理モデルである並進回転系離散的時間モデル７８、エンジン性能の物理モデルであるエンジン性能モデル８０を有する。
【００４８】
以下に、図３から図７を用いて、離散的時間の分割の内容、及び各離散的時間モデルの内容を説明する。図３は、統合離散的時間モデルにおけるサンプル点の様子を説明する図、図４と図５は、吸気系離散的時間モデルの様子を説明する図、図６と図７は、排気系離散的時間モデルの様子を説明する図である。
【００４９】
図３は、６気筒のエンジンシステムにおいて、各気筒の吸気、燃焼、排気がどのように進行して繰り返されるかを示し、そのサイクリックな現象に基づいて離散的時間に分割して、制御のためのサンプル点の位置を定める様子を説明する図である。図３（ａ）は、６つの気筒について示され、（ｂ）は、ＴＤＣを上死点、ＢＤＣを下死点として、３番気筒のクランク角度の変化を示したものである。
【００５０】
周知のように、６気筒のエンジンシステムは、クランクの７２０度回転を１サイクルとして、６つの気筒が順次に吸気−燃焼−排気を繰り返す。そこで、図３において、横軸は、クランクの１サイクルの７２０度回転を２０度ごとに区切ってある。１から３６の番号は、その分割された２０度の期間に付された順序番号である。縦軸は、吸気、燃焼、排気の３つの行程を並べてある。そして、各気筒を区別するために、気筒番号を１から６として、各気筒について、クランク角度と行程との関係が示されている。なお、ここでは、圧縮行程と膨張行程とをまとめて燃焼行程として示している。
【００５１】
例えば、３番気筒は、吸気行程が３１番目の期間から始まり、次の７２０度回転サイクルの５番目の期間で終わる。図３（ｂ）のクランク角度を参照すると、例えば、吸気行程は、ＴＤＣのところで吸気弁が開き（ＩＶＯ＝０ＢＴＤＣ：Ｉｎｔａｋｅ  ｖａｌｖｅ  ｏｐｅｎ  ０  ｄｅｇｒｅｅ  ｂｅｆｏｒｅ  ＴＤＣ）、ＢＤＣより４０度遅れて吸気弁が閉じる（ＩＶＣ＝４０ＡＢＤＣ：Ｉｎｔａｋｅ  ｖａｌｖｅ  ｃｌｏｓｅ  ４０  ｄｅｇｒｅｅ  ａｆｔｅｒ  ＢＤＣ）。そして燃焼行程は６番目の期間から始まり、１９番目の期間で終わる。また、排気行程は、２０番目の期間から始まり、３０番目の期間で終わる。そして次の吸気行程が３１番目の期間から始まる。したがって、吸気行程は合計で１１×２０度＝２２０度、燃焼行程は合計で１４×２０度＝２８０度、排気行程は合計で１１×２０度＝２２０度である。
【００５２】
このように、各気筒は、２２０度の吸気行程、２８０度の燃焼行程、２２０度の排気行程を気筒番号の順番で繰り返す。この繰り返しは、この６気筒エンジンシステムに共通のサイクリックな現象であり、エンジンパラメータ等に関係なく行われる。したがって、このサイクリックな現象に基づいて、７２０度の１サイクルを複数の離散的時間で分割すれば、エンジンパラメータ等に依存せずに、エンジンシステム全体の物理モデルを構成することができる。統合離散的時間モデルは、このような考えの下に、エンジンシステムの１サイクルを複数の離散的時間で分割して、物理モデルを構成するものである。
【００５３】
一般的には、７２０度を気筒数で除した分割を行うことができる。６気筒の場合は、７２０度／６＝１２０度で離散的時間の分割を行うことができる。また、エンジンシステムの１サイクルは、吸気−燃焼−排気の各行程で構成されるので、これらの行程の特徴に着目した離散的時間の分割を行うことが好ましい。
【００５４】
たとえば、図３に示されるように、３番気筒の吸気行程は、その前半において２番気筒と３番気筒がともに吸気行程にあり、後半においては３番気筒と４番気筒とがともに吸気行程にある。そして、その中間においては、３番気筒のみが吸気行程にある。このように、吸気行程は、３つの種類に分けることができる。１２０度の分割を考慮すると、前半と中間とをまとめ、後半と区別して２つの種類に分けることができる。ただし、中間と後半とをまとめ、前半と区別して２つの種類に分けるものとしてもよい。
【００５５】
したがって、離散的時間の分割に際しては、吸気行程における少なくともこの２つの種類を区別できるように、分割点、すなわちエンジンの状態量算出のためのサンプル点を設定することが好ましい。
【００５６】
また、３番気筒の排気行程も、吸気行程と同様に、その前半において２番気筒と３番気筒がともに排気行程にあり、後半においては３番気筒と４番気筒とがともに排気行程にある。そして、その中間においては、３番気筒のみが排気行程にある。したがって、排気行程においても、離散的時間の分割、すなわちサンプル点の設定は、排気行程における少なくとも２つの種類を区別できるように行われることが好ましい。
【００５７】
なお、燃焼行程においては、途中の計算値を他の気筒が利用することがないため、吸気行程及び排気行程で述べたような種類分けを行う必要がない。もちろん、ノック制御等で、他の気筒の計算結果を利用する必要があるときは、その状況にあわせ、燃焼行程を複数の種類に区別し、サンプル点の適切な設定を行うことができる。
【００５８】
この吸気行程の２種類の区別、排気行程の２種類の区別を考慮して、７２０度の１サイクルについて６つのサンプル点を設定すると、図３における期間番号で、６，１２，１８，２４，３０，３６の近傍がよい。この分割は１２０度ごとであるが、その前後±２０度の範囲のサンプル点を集めることで、上記の吸気行程と排気行程の種類の異なる内容を区別できる。このように、６気筒のエンジンシステムにおける離散的時間の分割は、１２０度ごとに、その前後±２０度の範囲で、サンプル点を集めるものとすることがよい。図３では、その離散的時間の分割におけるサンプル点を集める領域を太枠で囲み、ｋ＝０，１，２，３，４，５で示してある。また、計算のためのサンプル点の位置を右上がりの斜線で示してある。
【００５９】
図３の例では、エンジンシステム全体におけるサンプル点は、２０度おきの期間の番号で、（３５，３６，１），（５，６，７），（１１，１２，１３），（１７，１８，１９），（２３，２４，２５），（２９，３０，３１）である。このように離散的時間の分割を行い、サンプル点を設定することで、吸気行程の２種類の状態、排気行程の２種類の状態を区別したサンプル点を集めることができる。
【００６０】
このようにして設定されたサンプル点の集合について、隣接するサンプル点の間隔は、必ずしも１２０度となっていない。例えば、３番気筒についてサンプル点は、２０度おきの期間の番号で、５，１９，２５，３０，３６である。つまり、吸気行程の２種類の状態、排気行程の２種類の状態を区別し、燃焼行程は１種類で扱うことから、隣接する離散的時間のサンプル間隔は、エンジンシステム全体の１サイクル７２０度を気筒数の６で除した１２０度に対し、長短があることになる。
【００６１】
このようにして、エンジンシステム全体のサイクリックな現象に基づいて、吸気−燃焼−排気の１サイクルを複数の離散的時間で分割して、複数のサンプル点を集め、離散的時間状態を定義して統合離散的時間モデルを構成する。
【００６２】
統合離散的時間モデルは、状態量として質量ｍと圧力ｐを用い、次の式（１）で示される。
【数１】
2008175198.tif 000003
ここで、例えば、ｍ（ｌ，４）とは、図３におけるｌ，ｋ＝４のサンプル点における質量ｍを示す。同様に、例えば、ｐ（ｌ＋１，０）とは、図３におけるｌ＋１，ｋ＝０のサンプル点における圧力を示す。また、例えば、ＦC（ｌ，３）とは、図３におけるｌ，ｋ＝３のサンプル点における筒内流入燃料量を示す。Ａ，ｇ，ｂは、それぞれエンジンパラメータ等で定められるパラメータ行列で、Ａは８行８列、ｇは８行１列、ｂは８行６列の行列である。これらに付される添字は、行列の相違を示すためのものである。
【００６３】
そして、ｍは、吸気系と排気系と６つの気筒の合計８個の要素における質量状態を示す８行１列の行列である。ｐも、吸気系と排気系と６つの気筒の合計８個の要素における圧力状態を示す８行１列の行列である。ＦCは、６つの気筒における合計６個の要素における筒内流入燃料量状態を示す６行１列の行列である。ｍ，ｐ，ＦCは、次の式（２）で示される。ここで、添字のａは吸気系、ｅは排気系、ｃは気筒を示し、例えばｃ３は３番気筒を示す。
【数２】
2008175198.tif 000004
【００６４】
式（１）から分かるように、任意のサンプル点における質量ｍは、離散的時間の分割における１つ前のサンプル点における質量状態と、２つ前のサンプル点における質量状態と、２つ前のサンプル点における筒内流入燃料量等に基づいて記述することができる。また、任意のサンプル点における圧力ｐは、離散的時間の分割における１つ前のサンプル点における圧力状態と、２つ前のサンプル点における圧力状態等に基づいて記述することができる。このようにして、エンジンシステム全体の離散的時間のモデルとして、式（１）、（２）によって、統合離散的時間モデルを記述することができる。
【００６５】
図４、図５は、統合離散的時間モデルの下位階層モデルである吸気系離散的時間モデルの内容を説明する図である。吸気系離散的時間モデルは、スロットルから気筒に至るまでの吸気系について、質量保存則とエネルギ保存則に基づいて構築される。例えば、ある吸気弁のタイミングに固定した場合、微分方程式を解析的に解くことに代えて、サージタンクと気筒の２つの容積に注目し、吸気弁が開いた状態と、定常状態と仮定する吸気弁が閉じた状態との間で、逐時的に質量保存則とエネルギ保存則を適用して、離散的時間モデルを構成することができる。
【００６６】
図４は、吸気行程の最初において吸気弁が開くときの様子を示す図である。ここでは、注目している気筒１３が上死点近くにあるので、気筒１３の容積を小さく示してある。図５は、吸気行程の最後において吸気弁が閉じるときの様子を示す図である。ここでは、気筒１３が下死点近くにあるので、気筒１３の容積を大きく示してある。図４において白抜き矢印は、質量の移動を示し、斜線付き矢印は、熱あるいはシリンダの移動による仕事等のエネルギの移動を示す。図４、図５において、Ｐは圧力、Ｍは質量、Ｔは温度を示し、添字のｓはスタート（開始時期）、ｆはファイナル（終了時期）、ａは吸気系、ｃは気筒をそれぞれ示す。
【００６７】
このように、吸気行程においては、逐時的に、サージタンク１８と気筒１３との間の容積変化とともに、質量とエネルギの移動が生じる。そこで、吸気行程の離散的時間のそれぞれにおいて、サージタンク１８と気筒１３の間において、質量保存則とエネルギ保存則を適用して、吸気行程における離散的時間モデルを次の式（３）、（４）のように記述できる。ここで、吸気行程の内容には、上記のように、２つの種類があるので、その２つを区別して、式（３）は１つの気筒のみが吸気行程にある場合、式（４）は２つの気筒が吸気行程にある場合の離散的時間モデルの記述式をそれぞれ示している。
【数３】
2008175198.tif 000005
ここで、Ｍtは、スロットル通過流量である。２つの気筒を区別するために、添字に、ｊ，ｋを用いた。その他の記号、添字の意味は図４、図５で説明したものと同じである。行列の各要素は、エンジンパラメータ等から定めることができる。
【００６８】
式（３）から分かるように、１つの気筒が吸気行程にある場合には、吸気行程の最後のサンプル点におけるサージタンクと気筒における圧力と質量は、吸気行程の最初のサンプル点におけるサージタンクと気筒における圧力と質量、及びスロットル通過流量とに基づいて記述することができる。
【００６９】
また、式（４）から分かるように、２つの気筒が吸気行程にある場合には、吸気行程の最後のサンプル点におけるサージタンクと２つの気筒における圧力と質量は、吸気行程の最初のサンプル点におけるサージタンクと２つの気筒における圧力と質量、及びスロットル通過流量とに基づいて記述することができる。
【００７０】
式（３）、（４）は、吸気行程における最後のサンプル点についての記述式を示しているが、途中の離散的時間においても、同様な内容で、吸気行程における各要素の圧力と質量を記述することができる。
【００７１】
このように、吸気行程において、吸気弁が閉じた状態を定常状態と仮定することで計算を容易にすることができる。吸気弁の開閉タイミングを制御したい場合は、一般的な微分方程式を解析的に解けるように工夫することで、同様な離散的時間モデルを導き出せばよく、その他のモデルの修正をする必要がない。また、インテークマニホールド、インテークポート、吸気弁等の温度を状態量として加えることで、さらに精度を向上させることもできる。
【００７２】
次に、燃焼行程においては、吸気弁が閉じたときから、排気弁が開くときまでについて、質量保存則とエネルギ保存則に基づいて、離散的時間モデルを導き出すことになる。ここで、熱発生の計算には適当な周知の燃焼モデルを用いて計算することができるが、たとえば、周知のＷｉｅｂｅ関数を用いた場合、微分方程式を解析的に解くことができないので、一部近似を利用する。式（５）は、そのような近似を用い、燃焼行程における離散的時間モデルとして、圧力と質量について記述する式である。
【数４】
2008175198.tif 000006
ここで、Ｐは圧力、Ｍは質量、ＦCは筒内流入燃料量、θSは点火時期、αは空燃比、ＩＶＣは吸気弁が閉じた状態、ＥＶＯは排気弁が開いた状態を示す。添字の意味は、式（１）等で説明した内容と同じである。
【００７３】
式（５）から分かるように、排気弁が開くときの気筒における圧力は、吸気弁が閉じたときの気筒における圧力と、点火時期θSと、点火時期θSと筒内流入燃料量ＦCと空燃比αとで定まる関数とによって記述することができる。また、排気弁が開くときの気筒における質量は、吸気弁が閉じたときの気筒における質量と、筒内流入燃料量ＦCあるいは燃料噴射量Ｆiとで記述することができる。また、シリンダ壁等の温度を状態量として加えて、冷却損失を陽に考慮することで、さらに精度を向上させることができる。
【００７４】
式（５）は、燃焼行程における最後のサンプル点についての記述式を示しているが、途中の離散的時間においても、同様な内容で、燃焼行程における気筒の圧力と質量を記述することができる。
【００７５】
次に排気行程について説明する。図６、図７は、統合離散的時間モデルの下位階層モデルである排気系離散的時間モデルの内容を説明する図である。排気系離散的時間モデルは、気筒からエキゾーストマニホールド集合部に至るまでの排気系について、質量保存則とエネルギ保存則に基づいて構築される。例えば、ある排気弁のタイミングに固定した場合、微分方程式を解析的に解くことに代えて、気筒とエキゾーストマニホールド集合部の２つの容積に注目し、排気弁が開いた状態と、定常状態と仮定する排気弁が閉じた状態との間で逐時的に質量保存則とエネルギ保存則を適用して、離散的時間モデルを構成することができる。
【００７６】
図６は、排気行程の最初において排気弁が開くときの様子を示す図である。ここでは、注目している気筒１３が下死点近くにあるので、気筒１３の容積を大きく示してある。図７は、排気行程の最後において排気弁が閉じるときの様子を示す図である。ここでは、気筒１３が上死点近くにあるので、気筒１３の容積を小さく示してある。図６において白抜き矢印は、質量の移動を示し、斜線付き矢印は、熱あるいはシリンダの移動による仕事等のエネルギの移動を示す。図６、図７において、Ｐは圧力、Ｍは質量、Ｔは温度を示し、添字のｓはスタート（開始時期）、ｆはファイナル（終了時期）、ｅは排気系、ｃは気筒をそれぞれ示す。
【００７７】
このように、排気行程においては、逐時的に、気筒１３とエキゾーストマニホールド集合部２６との間の容積変化とともに、質量とエネルギの移動が生じる。そこで、排気行程の離散的時間のそれぞれにおいて、気筒１３とエキゾーストマニホールド集合部２６の間において、質量保存則とエネルギ保存則を適用して、排気行程における離散的時間モデルを次の式（６）、（７）のように記述できる。ここで、排気行程の内容には、上記のように、２つの種類があるので、その２つを区別して、式（６）は１つの気筒のみが排気行程にある場合、式（７）は２つの気筒が排気行程にある場合の離散的時間モデルの記述式をそれぞれ示している。
【数５】
2008175198.tif 000007
ここで、各記号、添字の意味は、図６、図７で説明したものと同じである。行列の各要素は、エンジンパラメータ等から定めることができる。
【００７８】
式（６）から分かるように、１つの気筒が排気行程にある場合には、排気行程の最後のサンプル点におけるエキゾーストマニホールド集合部と気筒における圧力と質量は、排気行程の最初のサンプル点におけるエキゾーストマニホールド集合部と気筒における圧力と質量とに基づいて記述することができる。
【００７９】
また、式（７）から分かるように、２つの気筒が排気行程にある場合には、排気行程の最後のサンプル点におけるエキゾーストマニホールド集合部と２つの気筒における圧力と質量は、排気行程の最初のサンプル点におけるエキゾーストマニホールド集合部と２つの気筒における圧力と質量とに基づいて記述することができる。
【００８０】
式（６）、（７）は、排気行程における最後のサンプル点についての記述式を示しているが、途中の離散的時間においても、同様な内容で、排気行程における各要素の圧力と質量を記述することができる。
【００８１】
このように、排気行程においても、排気弁が閉じた状態を定常状態と仮定することで計算を容易にすることができる。排気弁の開閉タイミングを制御したい場合は、一般的な微分方程式を解析的に解けるように工夫することで、同様な離散的時間モデルを導き出せばよく、その他のモデルの修正をする必要がない。また、エキゾーストマニホールド集合部、エキゾーストポート、排気弁等の温度を状態量として加えることで、さらに精度を向上させることもできる。
【００８２】
統合離散的時間モデルの下位階層のエンジン性能モデルは、エンジンシステム全体の１サイクル当りの図示仕事、トルク、燃料消費量等のエンジン性能を予測する物理モデルである。このエンジン性能モデルは、エンジンシステムの各状態量の逐時的変化に基づいて、サンプル間隔及び１サイクル当たりの図示仕事を計算することができる。この図示仕事から、エンジンシステム全体のトルク、回転数、出力、熱効率、燃料消費量等を予測することができる。
【００８３】
上記構成の内燃機関の制御装置を用いて、エンジンシステムの各要素における圧力と質量の逐時的変化を計算することができる。図８は、サージタンク内の圧力と質量についての計算結果の１例であり、図９は、気筒内の圧力と質量についての計算結果の１例である。これらの図において、横軸は時間をとり、各計算の点は、離散的時間の分割におけるサンプル点に対応する。○印で示されるものが、離散的時間モデルによる計算結果で、線分で示されるものが、連続時間での保存則を用いて計算した結果である。なお、気筒内における圧力と質量についての離散的時間モデルの計算は、燃焼途中のデータが未計算であるので、×印で示してある。図８、図９から分かるように、離散的時間モデルによる計算結果は、連続時間モデルによる計算結果とよい一致を示している。
【００８４】
上記構成の内燃機関の制御装置は、図１で説明したように、エンジンシステム１０に対する全体的な物理モデルである統合離散的時間モデル４４を有し、この統合離散的時間モデル４４に制御入力４２を適用してエンジン性能４６を算出し、あるいは、逆に、統合離散的時間モデル４４にエンジン性能４６を与えて、これを満たす制御入力４２を算出する機能を有する。したがって、この機能を、内燃機関のフィードバック制御、フィードフォワード制御に組み込むことができる。図１０、図１１は、それぞれ、エンジンシステム１０に対する回転数制御、トルク制御において、統合離散的時間モデルに基づく内燃機関の制御機能をフィードフォワード部９０，９１及びフィードバック部９２，９３に組み込んで構成する様子を示す図である。これらの組み込みに際しては、回転数、スロットル通過空気量、各気筒のピストン位置等の観測量から、エンジンシステムの各要素の状態を予測する手段の他に、エンジン性能である回転数、トルク、熱効率、燃料消費量等を予測する手段を含む。さらに、エンジン性能を目的である目標回転数レギュレーション、目標トルク追従等に合わせて最適にするように、複数の制御入力を同時に最適計算する手段を含む。
【実施例２】
【００８５】
上記では、式（１）に示すように、状態量である質量ｍと圧力ｐの統合離散的時間モデルを記述した。ここで質量ｍと圧力ｐは、吸気系と排気系と６つの気筒の合計８つの要素における質量状態、圧力状態を示す行列である。すなわち、式（１）は気筒ごとに状態量を記述しているために、大きな規模の関数または行列となり、サイクリックな現象を時不変系と見る場合、状態量についての演算量が大きな規模となっている。ところで、内燃機関の各気筒が同一の構造を有し、これらの各気筒がそれぞれ吸気−燃焼−排気のサイクリックな動作をしている場合には、気筒ごとに状態量を考えるのではなく、吸気−燃焼−排気のサイクリックな過程ごとに状態量を考えることで、この記述量を少なくすることができる。例えば、吸気の過程における状態量は、吸気にある気筒が第１気筒でも第２気筒であっても、あるいは第６気筒であっても同じく定義できる。同様に燃焼の過程における状態量も、燃焼過程にある気筒がどの気筒であっても同じで、排気の過程における状態量も、排気過程にある気筒がどの気筒であっても同じである。
【００８６】
したがって、統合離散的時間モデルの各サンプル点における状態量を、気筒ごとに考え、サイクリックな現象を時不変系と見る場合、気筒の数に応じて、大きな規模の関数または行列を必要とするが、これを内燃機関の吸気−燃焼−排気のサイクリックの遷移過程ごとに考えることにすれば、気筒数に関係なく小さな規模の関数または行列で済むことになる。内燃機関のサイクルは、吸気−燃焼−排気の３つで代表されるが、これを状態量の観点から見ると、バルブが開閉するときの状態量と、その途中の状態量とは異なるものと考えられる。このことから、例えば６気筒の場合は、内燃機関の吸気−燃焼−排気のサイクリックの遷移過程は、吸気途中−吸気終了−燃焼途中−燃焼終了−排気途中−排気終了の６つの遷移過程に区分することがよい。このように、統合離散的時間モデルの各サンプル点における状態量を、気筒ごとに考えるのではなく、内燃機関のサイクリックの複数の遷移過程ごとに考えることで、状態方程式の記述に要する関数または行列を小規模とすることができる。以下では、その内容を詳細に説明する。
【００８７】
最初に、気筒ごとの状態量と、遷移過程ごとの状態量の関係の一般的説明を行い、遷移過程ごとの状態量を説明し、その後に、気筒ごとの状態量と遷移過程ごとの状態量との間の変換等について説明する。
【００８８】
いま、式（１）を状態量ｘ、入力変数ｕ、状態量関数形をｆ、観測量をｙ、観測量関数形をｇとして、状態方程式を一般形に書き直すと次の式（８）で表せる。ここで、状態量ｘ、入力変数ｕについては、気筒ごとの状態量、入力変数であることを示すために、ｘ及びｕの頭部に山形（ハット）記号を付してある。
【数６】
2008175198.tif 000008
【００８９】
式（８）を、サンプル点を示すｋを用いて表すと式（９）のように示される。
【数７】
2008175198.tif 000009
【００９０】
ここで、式（９）の３番目に示されていることは、関数ｆl、kが、サンプル点を示すｋが異なるとその関数形を変えることを意味している。これは、図３に示されるように、サンプル点を示すｋが異なると、各気筒の状態がそれぞれ異なっていることから来ている。すなわち、各気筒ごとの状態量である気筒毎状態量を考える場合には、関数ｆl、kの関数形がサンプル点と共に、つまり時間と共に変化する。このことから、各気筒ごとの状態量の関数形は、内燃機関のサイクリックの時間経過とともに変化する時変系である、といえる。
【００９１】
これに対し、内燃機関の吸気−燃焼−排気のサイクリックな遷移過程ごとの状態量は、上記のように、気筒がどれであっても同じ内容であるので、気筒に依存しない。したがって、遷移過程ごとの状態量の関数形は、内燃機関のサイクリックの時間経過によっては変化せず、そのことから、時不変系である、といえる。そのことを式（９）と対比する形で式(１０)に示す。
【数８】
2008175198.tif 000010
【００９２】
式（１０）では、式（９）と同様に、状態量ｘ、入力変数ｕ、状態量関数形をｆ、観測量をｙ、観測量関数形をｇとして、状態方程式が示されている。ここで、状態量ｘ、入力変数ｕは、頭部に山形のハット記号を付した式（９）の気筒ごとのものと区別するため、ハット記号を付していない。そして、関数形ｆ、関数形ｇは、サンプル点を示すｋに依存せず、全てのｋに対し同一の関数形を用いることができる。このことが、上記で述べた時不変系を示している。
【００９３】
このように、内燃機関のサイクリックな遷移過程ごとの状態量である遷移過程毎状態量を考えることで、状態方程式を時不変系で記述することができる。次に、図３で説明した内燃機関のサイクリックな現象を例として、実際の遷移過程毎状態量について説明する。
【００９４】
図１２は、図３のうち、ｋ＝１の周辺について抜き出し、そのサイクリックな現象について、遷移過程毎状態量を説明する図である。なお、状態量とは、上記のように、質量ｍあるいは圧力ｐ等である。ここでは、遷移過程毎状態量として、６つの遷移過程ごとの状態量を示してある。すなわち、遷移過程として、吸気途中、吸気終了、燃焼途中、燃焼終了、排気途中、排気終了について、それぞれの状態量を、ｘ（３）、ｘ（２）、ｘ（１）、ｘ（６）、ｘ（５）、ｘ（４）が示されている。それぞれの括弧の中の数字は、６つの遷移過程を区別するものであるが、吸気途中が１番でなく３番となっている。これは、サイクリックな現象を説明するために用いた図３が、サンプル点を示す番号であるｋ＝０において、吸気途中の気筒が３番気筒である等の理由である。したがって、気筒番号と、サンプル点番号との決め事によって、遷移過程の番号の定義も変わってくる。ここでは、図３を用いて遷移過程を説明するので、上記のような番号付けになっている。
【００９５】
図１２の左側に示すサイクリックな現象を示す図においては、図３で説明したように、各気筒の状態の変化が示されている。この図で１とは１番気筒であり、２とは２番気筒である。したがって、ｋ＝１のサンプル点についてみると、吸気終了の気筒は、３番気筒である。また、燃焼終了の気筒は１番気筒である。排気終了の気筒は５番気筒である。同様に、吸気途中の気筒は４番気筒、燃焼途中の気筒は２番気筒、排気途中の気筒は６番気筒である。
【００９６】
上記のことを換言すれば、例えば、吸気終了の遷移過程における状態量ｘ（２）は、ｋ＝１において、３番気筒の状態量と同じ内容であることを示している。ここで、吸気終了の遷移過程における状態量ｘ（２）は、サンプル点が変わって、別の気筒の状態量に対応するものとなっても、その状態量の内容は変わらない。つまり、内燃機関を構成する各気筒の構造が同じで、かつ、内燃機関が全体として各気筒がサイクリックに動作している場合には、どのような局面においても、吸気終了の状態量は何番気筒であっても同じである。つまり、サンプル点であるｋが変化しても、吸気終了の状態量は変化しない。このことは、他の遷移過程の状態量についても同様である。
【００９７】
図１２の右の表では、気筒毎状態量として、１番気筒の状態量をｘ（１）、２番気筒の状態量をｘ（２）、３番気筒の状態量をｘ（３）、４番気筒の状態量をｘ（４）、５番気筒の状態量をｘ（５）、６番気筒の状態量をｘ（６）として示してある。ここで括弧の中の数字は気筒番号を示している。そして、それぞれを遷移過程毎状態量と区別するため頭部に山形のハット記号を付してある。例えば、ｋ＝１の場合、吸気途中の遷移過程の状態量は４番気筒の状態量と同じで、吸気終了の遷移過程の状態量は３番気筒の状態量と同じであること等が示されている。
【００９８】
図１２の右の表の下側には、遷移過程ごとに考えた入力変数の様子が示されている。ここで、入力変数とは、例えば、点火時期、燃料噴射量、スロットル開度等であるが、ここでは、点火時期と燃料噴射量が示されている。図３では、点火時期を各気筒ごとに考えていたが、これを遷移過程について考えることにすると、点火時期は、遷移過程においては燃焼開始の時期に決定される。また、燃料噴射量は、遷移過程においては、燃焼終了の時期で決定される。これをｋ＝１の場合で見ると、燃焼開始の時期にあるのは、３番気筒であり、燃焼終了に相当するのは、１番気筒である。
【００９９】
上記のことを換言すれば、例えば、点火時期の遷移過程における入力変数をｕsとして、これはｋ＝１において、３番気筒の入力変数である点火時期と同じ内容であることを示している。ここで、点火時期がどの遷移過程であるかは、サンプル点が変わって、別の気筒の点火時期に対応するものとなっても、常に燃焼開始のときである。つまり、内燃機関を構成する各気筒の構造が同じで、かつ、内燃機関が全体として各気筒がサイクリックに動作している場合には、どのような局面においても、点火時期は何番気筒であっても、同じように燃焼開始のときである。つまり、サンプル点であるｋが変化しても、入力変数である点火時期は、遷移過程としては常に燃焼開始となり、変化しない。このことは、他の入力変数、例えば、燃料噴射量等についても同様である。
【０１００】
図１２の右の表の下側には、気筒毎の点火時期として、３番気筒の点火時期をｕs（３）、１番気筒の燃料噴射量をｕf（１）として示してある。そして、それぞれを遷移過程ごとの入力変数である点火時期ｕs、燃料噴射量ｕfと区別するため頭部に山形のハット記号を付してある。
【０１０１】
図１２ではｋ＝１について、遷移過程毎状態量とこれに対応する気筒毎状態量、遷移過程で考えた入力変数とこれに対応する気筒ごとに考えた入力変数について説明したが、同様のことは、他のｋについても同様に対応付けを行うことができる。図１３は、各ｋについてまとめた様子を示す図である。ここでは、左側に遷移過程毎状態量、遷移過程で考えた入力変数が示され、これに対応する気筒毎状態量と、気筒で考えた入力変数とが、各サンプル点について示されている。例えば、サンプル点ｋ＝１においては、図１２で説明したように、吸気途中の遷移過程毎状態量が４番気筒の状態量に対応し、点火時期は３番気筒の入力変数に対応することが示されている。
【０１０２】
そして、サンプル点を移動させると、例えば、ｋ＝０で遷移過程としては吸気途中であった３番気筒は、ｋ＝１で遷移過程が吸気終了に移り、さらに、ｋ＝２で燃焼途中の遷移過程となり、ｋ＝３で燃焼終了の遷移過程、ｋ＝４で排気途中の遷移過程、ｋ＝５で排気終了の遷移過程となる。図１３では、３番気筒の移動状態を斜線で示したが、このように、各気筒は、サイクリックに各遷移過程を移動する。そして、例えば、吸気途中の状態量は、ｋ＝０で３番気筒、ｋ＝１で４番気筒、ｋ＝２で５番気筒、ｋ＝３で６番気筒、ｋ＝４で１番気筒、ｋ＝５で２番気筒というように、サイクリックに気筒が変わるが、どの気筒においても、その状態量は、吸気途中の状態量であって、気筒番号に依存しない。同様に、吸気終了等の他の遷移過程における状態量も、気筒番号に依存しない。勿論、吸気途中、吸気終了、燃焼途中、燃焼終了、排気途中、排気終了のそれぞれにおける状態量は相互において相違するが、状態量の関数形は同じであり、遷移過程の相違によって、その値が相違するのみである。
【０１０３】
入力変数もまた、遷移過程で考えることで、気筒番号に依存しない。例えば、点火時期は、上記のように、どの気筒においても燃焼開始の遷移過程である。図１３では、点火時期にある気筒が、サンプル点ｋが変わると変化する様子が示されている。すなわち、ｋ＝０では２番気筒、ｋ＝１では３番気筒、ｋ＝２では４番気筒、ｋ＝３では５番気筒、ｋ＝４では６番気筒、ｋ＝５では１番気筒が点火時期にある気筒である。このように、点火時期にある気筒はサイクリックに変化するが、この場合でも、遷移過程で考えると、どの気筒においても、燃焼開始のときである。
【０１０４】
上記のように、遷移過程毎状態量は、気筒毎状態量とサイクリックに変化する対応関係で対応付けられている。また、遷移過程で考えた入力変数も、サイクリックに変化する気筒番号と対応付けられる。したがって、内燃機関のサイクリックな状態変化に基づく対応関係を対応付け関数として示すことで、その対応付け関数を用いて、遷移過程毎状態量と気筒毎状態量との間の変換を演算によって行うことができる。同様に、遷移過程で考えた入力変数と、対応する気筒に対する入力変数との間の変換を演算によって行うことができる。
【０１０５】
遷移過程毎状態量と気筒毎状態量との間の変換は、変換行列を用いて行うことができる。状態量として、質量ｍ、圧力ｐを選び、変換行列をＺとすると、変換関係式は次の式（１１）で示すことができる。なお、ここでは、気筒毎状態量について、遷移過程毎状態量と区別するため、頭部に山形のハット記号を付してある。
【数９】
2008175198.tif 000011
【０１０６】
ここで、状態量である質量ｍ、圧力ｐは、吸気系と排気系と複数の気筒の行列である。例えば、６気筒の場合、８行１列の行列である。ここでは、吸気系と排気系を省略し、気筒についてのみ考えると、次の式（１２）のように、各状態量は６行１列の行列で示される。
【数１０】
2008175198.tif 000012
【０１０７】
変換行列Ｚは、ｋ＝０からｋ＝５の間を順次サイクリックに変換するものであるので、直交行列の一種である置換行列を用いることができる。すなわち、置換行列をＱとして、ｋに応じて、そのべき乗数を変更してサイクリックに変換を実行できる変換行列とすることができる。その様子を、図１４に示す。ここでは、べき乗数をｋの関数という意味で、ｋの頭部に横線（バー）が付されている。べき乗数は、サンプル点が所定数であるので、その所定数で変換行列が繰り返すように設定される。例えば、図１４の場合、ｋ＝０，１，２，３，４，５の６つがサンプル点であるので、べき乗数は、最小が０で、最大が５となる。
【０１０８】
べき乗数は、上記のように、サンプル点の数で上限が決められるので、これをｋｍｏｄＡで定義できる。ｋｍｏｄＡとは、数ｋと数Ａとの間で、ｋ／Ａを計算したときの余りの数である。例えば、上記のようにＡ＝６の場合は、ｋ＝２のときのｋｍｏｄ６は、２／６の余りであるので、２である。このような変換行列Ｚは、次の式（１３）で与えることができる。
【数１１】
2008175198.tif 000013
【０１０９】
式（１３）を満たす変換行列Ｚは、図１３の対応関係を満たすような適当な置換行列Ｑを求めることで得られる。６気筒の場合、置換行列Ｑは、次の式（１４）を満たす必要がある。
【数１２】
2008175198.tif 000014
【０１１０】
図１３の場合においては、遷移過程の種類が６つ、気筒数が６つであるので、置換行列Ｑは６行６列の行列になる。また、サンプル点の数は６つであるので、置換行列Ｑのべき乗であるｋｍｏｄＡにおいて、Ａ＝６である。これらの条件の下で、図１３の対応関係を満たす置換行列Ｑの１例を図１５に示す。
【０１１１】
このようにして、置換行列Ｑと、べき乗であるｋｍｏｄＡを用いて、変換行列Ｚを定めることができる。そしてこの変換行列Ｚを用いて、遷移過程毎の状態量と、気筒毎の状態量との間を演算により変換することができる。
【０１１２】
遷移過程を考えた入力変数と、これに対応する気筒における入力変数との間の変換を表す次の式（１５）は、作用素Ψを用いて行うことができる。作用素Ψは、例えば点火時期である燃焼開始時期が、任意のｋにおいて、どの気筒の点火時期に対応するかを定める演算要素である。その関係は、式（１６）で示すことができる。
【数１３】
2008175198.tif 000015
【数１４】
2008175198.tif 000016
【０１１３】
図１３で示す例では、点火時期に対応する遷移過程は、ｋ＝０で２番気筒、ｋ＝１で３番気筒で、以下、ｋの変化と共に気筒番号が進んでゆく。したがって、作用素Ψは、上記のｋｍｏｄＡの関数を用いて、（（ｋ＋１）ｍｏｄ６）＋１と示すことができる。同様に、燃料噴射量に対応する遷移過程は、ｋ＝０で６番気筒、ｋ＝１で１番気筒、ｋ＝２で２番気筒であり、以下、ｋの変化と共に気筒番号が進んでゆく。したがって、これもｋｍｏｄＡの関数を用いて表すことができる。上記の式（１６）は、そのことをまとめたものである。
【０１１４】
これらが図１５にまとめてある。ここでは、上記のように、遷移過程毎状態量と、気筒毎状態量との間を変換する状態量変換手段として、変換行列Ｚが示されている。そして、図１３で説明した対応付けを変換する変換行列Ｑの１例が示されている。また、遷移過程を考えた制御入力と、気筒毎の制御入力との間を変換する制御入力変換手段として、作用素Ψが示されている。そして、図１３で説明した対応付けを変換する作用素Ψの１例が示されている。この変換行列、作用素を用いることで、演算により、遷移過程毎の状態変化に基づく状態方程式と、気筒毎の状態変化に基づく状態方程式との間の変換を演算により実行することができる。
【図面の簡単な説明】
【０１１５】
【図１】本発明に係る実施の形態における内燃機関の制御装置と、その制御対象の内燃機関である６気筒エンジンシステムとを示す図である。
【図２】本発明に係る実施の形態におけるＥＣＵの詳細な構成を説明するブロック図である。
【図３】本発明に係る実施の形態において、統合離散的時間モデルにおけるサンプル点の様子を説明する図である。
【図４】本発明に係る実施の形態の吸気系離散的時間モデルにおいて、吸気弁が開くときの様子を説明する図である。
【図５】本発明に係る実施の形態の吸気系離散的時間モデルにおいて、吸気弁が閉じるときの様子を説明する図である。
【図６】本発明に係る実施の形態の排気系離散的時間モデルにおいて、排気弁が開くときの様子を説明する図である。
【図７】本発明に係る実施の形態の排気系離散的時間モデルにおいて、排気弁が閉じるときの様子を説明する図である。
【図８】本発明に係る実施の形態において、サージタンク内の圧力と質量についての計算結果の１例を示す図である。
【図９】本発明に係る実施の形態において、気筒内の圧力と質量についての計算結果の１例を示す図である。
【図１０】本発明に係る実施の形態における統合離散的時間モデルに基づく内燃機関の制御機能を、エンジンシステムに対する回転数制御のフィードフォワード部及びフィードバック部に組み込む様子を示す図である。
【図１１】本発明に係る実施の形態における統合離散的時間モデルに基づく内燃機関の制御機能を、エンジンシステムに対するトルク制御のフィードフォワード部及びフィードバック部に組み込む様子を示す図である。
【図１２】他の実施形態において、遷移過程毎状態量を説明する図である。
【図１３】他の実施形態において、遷移過程毎状態量と気筒毎状態量の対応付け等を説明する図である。
【図１４】他の実施形態において、変換行列の様子を説明する図である。
【図１５】他の実施形態において、変換行列と作用素の内容をまとめた図である。
【符号の説明】
【０１１６】
１０  エンジンシステム、１２，１３  気筒、１４  エアフィルタ、１６  スロットル、１８  サージタンク、２０  吸気管、２１  燃料噴射弁、２２  点火栓、２４  回転角度センサ、２６  エキゾーストマニホールド集合部、２８  触媒、３０  内燃機関の制御装置、３２  目標指令入力部、３４  観測量取得部、３６  制御出力部、４０  ＥＣＵ（エンジン制御ユニット）、４２  制御入力、４４  統合離散的時間モデル、４６  性能算出部、５０  エンジンパラメータ、６０  状態量算出部、６２  物理モデル作成記憶部、６６  離散的時間分割部、６８  物理モデル作成部、７０  吸気系離散的時間モデル、７２  燃焼系離散的時間モデル、７４  排気系離散的時間モデル、７６  燃料系離散的時間モデル、７８  並進回転系離散的時間モデル、８０  エンジン性能モデル、８２  質量・圧力・回転角度算出部、９０，９１  フィードフォワード部、９２，９３  フィードバック部。
【図１】
2008175198.tif 000017
【図２】
2008175198.tif 000018
【図３】
2008175198.tif 000019
【図４】
2008175198.tif 000020
【図５】
2008175198.tif 000021
【図６】
2008175198.tif 000022
【図７】
2008175198.tif 000023
【図８】
2008175198.tif 000024
【図９】
2008175198.tif 000025
【図１０】
2008175198.tif 000026
【図１１】
2008175198.tif 000027
【図１２】
2008175198.tif 000028
【図１３】
2008175198.tif 000029
【図１４】
2008175198.tif 000030
【図１５】
2008175198.tif 000031