(57)【要約】      （修正有）【課題】将来の予測値を計算するためのモデルをＰＷＡＲＸモデルで表現する場合において、演算負荷を軽減して演算資源の節約を図る。【解決手段】予測値ｙ＾（ｗｃａｎｄ＋Δ），ｙ＾（ｗｃａｎｄ−Δ）の計算（ステップＳ１０）と集積点判断（ステップＳ１２）とが１ステップごとに繰り返される。そして、集積点に到達したと判断された場合（ステップＳ１２においてＹｅｓの場合）、それまでに計算した予測値ｙ＾（ｗｃａｎｄ＋Δ），ｙ＾（ｗｃａｎｄ−Δ）に基づいて勾配∇の計算が行われる（ステップＳ１４）。【選択図】図８
2016130877.tif 000002
【特許請求の範囲】
【請求項１】
プラントの制御量の出力値をその目標値に近づけるようにフィードバック制御によって前記プラントの制御入力を決定するフィードバックコントローラと、
前記プラントと前記フィードバックコントローラとを含む閉ループシステムのモデルを用いて前記プラントの特定状態量の将来の予測値を計算すると共に、計算した予測値と、前記特定状態量に課せられた制約と、オリジナルの目標値と、修正目標値の候補とを用いて表される評価関数の最適値を最急降下法による繰り返し計算によって探索し、前記最適値を用いて前記フィードバックコントローラに与えられる前記目標値を修正するリファレンスガバナと、を備え、
前記モデルは、前記特定状態量の拡大状態ベクトルに基づいて特定される区分領域に応じて、複数のサブモデルを切り替え可能なＰＷＡＲＸモデルによって表現され、
前記リファレンスガバナは、前記オリジナルの目標値に基づいて前記修正目標値の候補を用意し、用意した修正目標値の候補から正方向および負方向に所定値だけ離れた２つの近傍値のそれぞれを前記ＰＷＡＲＸモデルに適用して前記特定状態量の予測値を２つ算出し、算出した２つの予測値に基づいて前記修正目標値の候補に対する前記評価関数の勾配を算出し、予め設定した回数に亘って算出した前記評価関数の勾配を最小化する前記修正目標値の候補を前記最適値として決定するように構成され、
前記リファレンスガバナは、算出した２つの予測値の何れもが前記サブモデルの内の少なくとも３つが接する点であって、尚且つ、隣り合う区分領域のベクトル場が互いに向かい合う集積点に到達した場合、前記特定状態量の予測値の算出を終了し、それまでに算出した前記評価関数の勾配に基づいて前記最適値を決定するように構成されていることを特徴とするプラント制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、プラント制御装置に関する。
【背景技術】
【０００２】
従来、例えば特許文献１には、プラントの制御量の出力値を目標値に近づけるようにフィードバック制御によって当該プラントの制御入力を決定するフィードバックコントローラと、当該プラントと当該フィードバックコントローラとを含む閉ループシステムのモデルを用いて当該プラントの特定状態量の将来の予測値を計算し、当該予測値と当該特定状態量に課せられた制約とに基づいて当該フィードバックコントローラに与えられる目標値を修正するリファレンスガバナと、を備えるプラント制御装置が開示されている。
【０００３】
このプラント制御装置において、リファレンスガバナは、具体的に、最急降下法を用いた繰り返し計算によって、特定状態量ｙの将来の予測値ｙ＾および制約ｙ−と、修正目標値候補ｗｃａｎｄと、オリジナルの目標値ｒとを用いて表される評価関数Ｊ（ｗｃａｎｄ）の最適値を探索する。また、リファレンスガバナは、この評価関数Ｊ（ｗｃａｎｄ）を最適化する修正目標値候補ｗｃａｎｄを、制御量の最終的な修正目標値ｗとして決定する。このようにしてオリジナルの目標値ｒが修正され、最終的な修正目標値ｗがフィードバックコントローラに入力される。
【０００４】
リファレンスガバナによる最終的な修正目標値ｗの決定手法は、次の通りである。先ず、オリジナルの目標値ｒに基づいて、修正目標値候補ｗｃａｎｄを複数用意する。続いて、修正目標値候補ｗｃａｎｄから微小摂動Δだけ離れた２つの近傍値ｗｃａｎｄ＋Δ，ｗｃａｎｄ−Δに対する予測値ｙ＾（ｗｃａｎｄ＋Δ），ｙ＾（ｗｃａｎｄ−Δ）を上述した予測モデルによって予測し、勾配∇（ｗｃａｎｄ）（＝｛Ｊ（ｗｃａｎｄ＋Δ）−Ｊ（ｗｃａｎｄ−Δ）｝／２Δ）を計算することを有限回（予め設定した回数）繰り返す。最終的な修正目標値ｗは、計算した勾配∇（ｗｃａｎｄ）を最小化する修正目標値候補ｗｃａｎｄとして決定される。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１４−１２７０８３号公報
【非特許文献】
【０００６】
【非特許文献１】秋田敏和ら、「ＰＷＡＲＸモデルによる前方車追従行動のモデル化と解析」、第５０回自動制御連合講演会、２００７年１１月２４日〜２５日
【発明の概要】
【発明が解決しようとする課題】
【０００７】
ところで、例えば非特許文献１に開示されている区分的アフィンＡＲＸモデル（ＰＷＡＲＸモデル）は、離散的な状態遷移に応じて複数のサブモデル（ＡＲＸモデル）を切り替えるシステム表現手法であり、非線形性が強くモデル化の難しいエンジンのようなプラントを扱うのに有用である。そこで、本発明者は、上述したリファレンスガバナによる修正目標値ｗの決定手法において、予測値ｙ＾を計算するための予測モデルをこのＰＷＡＲＸモデルで表現することを検討している。
【０００８】
しかし、ＰＷＡＲＸモデルはサブモデルで記述される軌道の流れ（ベクトル場）によって動きが規定されるため、各サブモデルの担当領域の交わるところで集積点となる可能性がある。また、このような集積点に予測値ｙ＾（ｗｃａｎｄ＋Δ），ｙ＾（ｗｃａｎｄ−Δ）の何れもが到達した場合はどこにも動くことができず、その集積点に留まり続けることになる。そうすると、同一集積点で勾配計算を続けることになるので、演算負荷が過大な状態が続き、演算資源を無駄に使用し続けてしまう。
【０００９】
本発明は、上述のような課題に鑑みてなされたものである。即ち、将来の予測値ｙ＾を計算するためのモデルをＰＷＡＲＸモデルで表現する場合において、演算負荷を軽減して演算資源の節約を図ることを目的とする。
【課題を解決するための手段】
【００１０】
本発明は、上記の目的を達成するため、プラント制御装置であって、
プラントの制御量の出力値をその目標値に近づけるようにフィードバック制御によって前記プラントの制御入力を決定するフィードバックコントローラと、
前記プラントと前記フィードバックコントローラとを含む閉ループシステムのモデルを用いて前記プラントの特定状態量の将来の予測値を計算すると共に、計算した予測値と、前記特定状態量に課せられた制約と、オリジナルの目標値と、修正目標値の候補とを用いて表される評価関数の最適値を最急降下法による繰り返し計算によって探索し、前記最適値を用いて前記フィードバックコントローラに与えられる前記目標値を修正するリファレンスガバナと、を備え、
前記モデルは、前記特定状態量の拡大状態ベクトルに基づいて特定される区分領域に応じて、複数のサブモデルを切り替え可能なＰＷＡＲＸモデルによって表現され、
前記リファレンスガバナは、前記オリジナルの目標値に基づいて前記修正目標値の候補を用意し、用意した修正目標値の候補から正方向および負方向に所定値だけ離れた２つの近傍値のそれぞれを前記ＰＷＡＲＸモデルに適用して前記特定状態量の予測値を２つ算出し、算出した２つの予測値に基づいて前記修正目標値の候補に対する前記評価関数の勾配を算出し、予め設定した回数に亘って算出した前記評価関数の勾配を最小化する前記修正目標値の候補を前記最適値として決定するように構成され、
前記リファレンスガバナは、算出した２つの予測値の何れもが前記サブモデルの内の少なくとも３つが接する点であって、尚且つ、隣り合う区分領域のベクトル場が互いに向かい合う集積点に到達した場合、前記特定状態量の予測値の算出を終了し、それまでに算出した前記評価関数の勾配に基づいて前記最適値を決定するように構成されていることを特徴とする。
【発明の効果】
【００１１】
本発明によれば、修正目標値の候補から正方向および負方向に所定値だけ離れた２つの近傍値のそれぞれをＰＷＡＲＸモデルに適用して求めた特定状態量の２つの予測値の何れもが集積点に到達した場合、特定状態量の予測値の算出を終了し、それまでに算出した勾配に基づいて評価関数の最適値を決定するので、集積点到達後の無駄な将来予測や勾配計算を省くことができ、演算資源を節約できる。
【図面の簡単な説明】
【００１２】
【図１】実施の形態に係る制御装置が適用されるシステムの構成を示す図である。
【図２】実施の形態に係る制御装置が有する目標値追従制御構造を示す図である。
【図３】図２に示す目標値追従制御構造を等価変形して得られたフィードフォワード構造である。
【図４】評価関数Ｊ（ｗ）の勾配∇を説明するための図である。
【図５】実施の形態に係るリファレンスガバナアルゴリズムを示す図である。
【図６】集積点の問題点を説明するための図である。
【図７】集積点の問題点を説明するための図である。
【図８】実施の形態における勾配演算アルゴリズムを示すフローチャートである。
【図９】図２に示す目標値追従制御構造を適用可能なディーゼルエンジンの入出力の例を示す図である。
【発明を実施するための形態】
【００１３】
本発明の実施の形態のプラント制御装置は、車両動力プラントであるディーゼルエンジンの制御装置である。図１は、本実施の形態に係る制御装置が適用されるシステムの構成を示す図である。図１に示すように、ディーゼルエンジンの本体２には４つの気筒が直列に備えられ、気筒ごとにインジェクタ８が設けられている。エンジン本体２には吸気マニホールド４と排気マニホールド６が取り付けられている。
【００１４】
吸気マニホールド４にはエアクリーナ２０から取り込まれた新気が流れる吸気通路１０が接続されている。吸気通路１０にはターボ過給機１４のコンプレッサが取り付けられている。このコンプレッサの下流にはインタークーラ２２が備えられ、インタークーラ２２の下流にはディーゼルスロットル２４が設けられている。排気マニホールド６にはエンジン本体２からの排気を大気中に放出するための排気通路１２が接続されている。排気通路１２にはターボ過給機１４のタービンが取り付けられている。ターボ過給機１４は可変容量型であって、タービンには可変ノズル１６が備えられている。
【００１５】
図１に示すシステムは、排気系から吸気系へ排気を再循環させるＥＧＲ装置を備えている。ＥＧＲ装置は、吸気通路１０におけるディーゼルスロットル２４の下流と排気マニホールド６とをＥＧＲ通路３０によって接続する高圧ループＥＧＲ装置である。ＥＧＲ通路３０にはＥＧＲ弁３２が設けられている。但し、ＥＧＲ装置は、吸気通路１０におけるコンプレッサの上流と、排気通路１２におけるタービンの下流とを、ＥＧＲ通路３０とは別のＥＧＲ通路によって接続する低圧ループＥＧＲ装置であってもよい。
【００１６】
図１に示すＥＣＵ（Electronic Control Unit）４０が本実施の形態に係る制御装置に相当する。ＥＣＵ４０は、ＲＡＭ（ランダムアクセスメモリ）、ＲＯＭ（リードオンリメモリ）、ＣＰＵ（マイクロプロセッサ）等を備えている。ＥＣＵ４０は、各種センサの信号を取り込み処理する。各種センサには、エンジン回転速度を検出する回転速度センサ４２や、アクセルペダルの開度に応じた信号を出力するアクセルペダル開度センサ４４などが含まれている。ＥＣＵ４０は、取り込んだ各センサの信号を処理して所定の制御プログラムに従ってアクチュエータを操作する。ＥＣＵ４０によって操作されるアクチュエータには、可変ノズル１６、ディーゼルスロットル２４、ＥＧＲ弁３２などが含まれている。
【００１７】
本実施の形態において、ＥＣＵ４０は、ディーゼルエンジンの過給圧・ＥＧＲ率制御を実行する。過給圧・ＥＧＲ率制御における制御入力（操作量）は可変ノズル開度、ＥＧＲ弁開度およびディーゼルスロットル開度であり、制御出力（特定状態量）は過給圧とＥＧＲ率である。ここで、過給圧とＥＧＲ率にはハード上或いは制御上の制約が課せられている。ＥＣＵ４０は、過給圧とＥＧＲ率がそれぞれの制約を満たし、尚且つ、それぞれの目標値に追従させるように制御入力を決定する。
【００１８】
図２は本実施の形態に係る制御装置が有する目標値追従制御構造を示す図である。なお、図２に示す目的値追従制御構造は、ＥＣＵ４０のＲＯＭに格納された制御プログラムに従いＣＰＵが動作することで仮想的に実現される構成である。この目標値追従制御構造は、目標値マップ（ＭＡＰ）、リファレンスガバナ（ＲＧ）およびフィードバックコントローラ（ＦＢＣ）を備えている。
【００１９】
目標値マップは、ディーゼルエンジン（ＤＥ）の運転条件を示す外生入力ｄ＝［エンジン回転速度；燃料噴射量］が与えられると、ディーゼルエンジンの制御量の目標値ｒ＝［ＥＧＲ率目標値；過給圧目標値］を出力する。
【００２０】
リファレンスガバナは、目標値ｒが与えられると、制御出力（特定状態量）ｙ＝［ＥＧＲ率；過給圧］に関する制約が満たされるように目標値ｒを修正し、修正目標値ｗ＝［ＥＧＲ率修正目標値；過給圧修正目標値］を出力する。リファレンスガバナの詳細については後述する。
【００２１】
フィードバックコントローラは、リファレンスガバナから修正目標値ｗが与えられると、ディーゼルエンジンの状態量ｘ＝［ＥＧＲ率；過給圧］を修正目標値ｗに近づけるように、フィードバック制御によってディーゼルエンジンの制御入力ｕ＝［ディーゼルスロットル開度；ＥＧＲ弁開度；可変ノズル開度］を決定する。フィードバックコントローラの仕様に限定はなく、公知のフィードバックコントローラを用いることができる。
【００２２】
図３は図２に示す目標値追従制御構造を等価変形して得られたフィードフォワード構造である。図２において破線で囲まれた閉ループシステムは既に設計済みであるとして、図３に示すフィードフォワード構造では１つのモデル（Ｐ）とされている。閉ループシステムのモデルは次のモデル式（１）で表される。式（１）において、ｆ，ｈはモデル式の関数である。また、ｋは離散時間ステップを表している。
【数１】
2016130877.tif 000003
【００２３】
ところで、上述したＰＷＡＲＸモデルは、非線形システムの表現能力の高さと、データクラスタリング理論やシステム同定理論との親和性が高く、非常によく用いられるシステムモデル形式である。上記閉ループシステムに対するＰＷＡＲＸモデルは、次式（２）で記述される。式（２）において、ｓはサブモデルの切り替え総数であり、θ１，・・・θｓは各サブモデルの係数であり、Ｂｗ，Ｂｄは既知の係数である。また、ｘ−（ｋ）は拡大状態ベクトル（ｘ−（ｋ）＝［ｙ（ｋ−１），ｙ（ｋ−２），…，ｙ（ｋ−Ｎｚ），１］Ｔ）である（Ｎｚは整数、Ｔは転置記号）。
【数２】
2016130877.tif 000004
【００２４】
リファレンスガバナは、上記の式（１），（２）で表される予測モデルを用いてディーゼルエンジンの制御出力ｙの将来の予測値ｙ^を計算する。上述したように、本実施の形態における制御出力ｙはＥＧＲ率および過給圧であり、制御出力ｙには制約が課せられている。制御出力ｙがその上限値ｙ−以下であることが制御出力ｙに課せられた制約である。制御出力予測値ｙ^の計算には、状態量ｘおよび外生入力ｄに加えて修正目標値ｗが用いられる。リファレンスガバナは、制御出力予測値ｙ^と制御出力上限値ｙ−とに基づき、次式（３）で表される評価関数Ｊ（ｗ）を用いて修正目標値ｗを計算する。
【数３】
2016130877.tif 000005
【００２５】
式（３）に示す評価関数Ｊ（ｗ）の右辺第１項は修正目標値候補ｗｃａｎｄを変数とする目的関数である。この目標関数はオリジナルの目標値ｒと修正目標値候補ｗｃａｎｄとの距離が小さいほど小さな値を取るように構成されている。評価関数Ｊ（ｗ）の右辺第２項はペナルティ関数である。ペナルティ関数は制御出力予測値ｙ^が制約ｙ−に抵触する場合に目的関数にペナルティを加えるように構成されている。ペナルティ関数には、ペナルティに重みを付けるための重み定数ρが設定されている。なお、右辺第２項のｙ^（ｋ＋ｉ|ｋ）は時刻ｋの時点での情報に基づく時刻ｋ＋ｉの時点の制御出力予測値を表し、Ｎｈは予測ホライズン（予測ステップ数）を表している。
【００２６】
また、本実施の形態においては、勾配法（具体的には最急降下法）を用いて式（３）の評価関数Ｊ（ｗ）を最小化する修正目標値ｗを探索している。この勾配演算では、修正目標値候補ｗｃａｎｄに対して微小摂動Δ（正の実数）を加えたときのＪ（ｗｃａｎｄ＋Δ）と、微小摂動Δを減じたときのＪ（ｗｃａｎｄ−Δ）との差を２Δで割ることで勾配∇（ｗｃａｎｄ）を算出している。つまり、勾配∇（ｗｃａｎｄ）＝｛Ｊ（ｗｃａｎｄ＋Δ）−Ｊ（ｗｃａｎｄ−Δ）｝／２Δである。図４は、評価関数Ｊ（ｗｃａｎｄ）の勾配∇（ｗｃａｎｄ）を説明するための図である。図４に示すように、勾配∇（ｗｃａｎｄ）は、修正目標値候補ｗｃａｎｄの近傍のｗｃａｎｄ＋Δ，ｗｃａｎｄ−Δの評価関数Ｊ（ｗｃａｎｄ＋Δ），Ｊ（ｗｃａｎｄ−Δ）を通る直線の傾きとして表される。なお、図４においては、評価関数Ｊ（ｗｃａｎｄ）の構造が一次元の場合を示しているが、本実施の形態では目標値ｒがＥＧＲ率と過給圧の２つあるので、実際の評価関数Ｊ（ｗｃａｎｄ）の構造は二次元となる。
【００２７】
図５は、本実施の形態に係るリファレンスガバナアルゴリズムを示す図である。図５に示すように、本実施の形態では、修正目標値候補ｗｃａｎｄに対して、現在のディーゼルエンジンの運転条件に基づく閉ループシステムの将来予測と、この予測結果を用いた評価関数Ｊ（ｗｃａｎｄ）の計算（つまり、上述した勾配∇（ｗｃａｎｄ）の計算）が有限回（予め設定した回数）反復される。そして、計算した評価関数Ｊ（ｗｃａｎｄ）を最小化する修正目標値候補ｗｃａｎｄが選択され、選択された修正目標値候補ｗｃａｎｄが最終的な修正目標値ｗとして決定される。
【００２８】
ところで、式（３）に示す評価関数の右辺第２項の制御出力予測値ｙ^（ｋ＋ｉ|ｋ）は、修正目標値候補ｗｃａｎｄの近傍値ｗｃａｎｄ＋Δ，ｗｃａｎｄ−Δのそれぞれに対して計算されるので、制御出力予測値ｙ^は、予測値ｙ＾（ｗｃａｎｄ＋Δ），ｙ＾（ｗｃａｎｄ−Δ）として表されることになる。ここで、式（２）に示したように、ＰＷＡＲＸシステムはｘ-（ｋ）が属する区分領域Ｘ１〜Ｘｓに応じて使用するサブモデルの切り替えが行われる。また、ＰＷＡＲＸモデルはサブモデルで記述される軌道の流れ（ベクトル場）によって動きが規定される。そのため、各サブモデルの担当領域の交わるところが集積点となる可能性がある。この「集積点」とは、サブモデルの内の少なくとも３つが接する点であって、尚且つ、隣り合う区分領域のベクトル場が互いに向かい合う点を意味する。
【００２９】
図６乃至図７は、集積点の問題点を説明するための図である。なお、図６乃至図７においては、５ステップ先まで将来予測を行うとしている。図６に示すように、区分領域Ｘ１，Ｘ２，Ｘ３の間には交点が存在し、これらの区分領域のベクトル場が互いに向かい合っている場合、上述した予測値ｙ＾（ｗｃａｎｄ＋Δ），ｙ＾（ｗｃａｎｄ−Δ）の計算値がこの交点に収束してしまう。このように状態軌道が集積点に到達した場合、予測値ｙ＾（ｗｃａｎｄ＋Δ），ｙ＾（ｗｃａｎｄ−Δ）の計算値がその集積点に留まり続けることになる。そうすると、同一集積点で勾配∇（ｗｃａｎｄ）の計算が続くことになる。図７に示すように、集積点到達以降は、上述した勾配∇（ｗｃａｎｄ）の計算式（＝｛Ｊ（ｗｃａｎｄ＋Δ）−Ｊ（ｗｃａｎｄ−Δ）｝／２Δ）の分子がゼロであるにも関わらず、予測値ｙ＾（ｗｃａｎｄ＋Δ），ｙ＾（ｗｃａｎｄ−Δ）や勾配∇（ｗｃａｎｄ）が計算されてしまう。このように、将来予測や勾配∇（ｗｃａｎｄ）計算を続けることは、ＣＰＵの演算負荷が過大な状態が続くことに他ならず、演算資源を無駄にしてしまう。
【００３０】
そこで、本実施の形態では、予測値ｙ＾（ｗｃａｎｄ＋Δ），ｙ＾（ｗｃａｎｄ−Δ）の計算値が集積点に到達した場合には、到達以降の将来予測をストップして、勾配∇（ｗｃａｎｄ）の評価に入ることとしている。集積点到達後の将来予測をストップすれば、無駄な将来予測や勾配計算を省くことができるので演算資源を節約できる。なお、集積点がどの点で起きるかは、式（２）に示したサブモデルの設定時において集積点解析を行うことで事前に把握できる。本実施の形態では、このようにして把握した集積点に関する情報がリストとしてＥＣＵ４０のメモリに記憶されているものとする。
【００３１】
図８は、本実施の形態における勾配演算アルゴリズムを示すフローチャートである。図８に示すように、本実施の形態では、将来予測（ステップＳ１０）と集積点判断（ステップＳ１２）とが１ステップごとに繰り返される。なお、ステップＳ１０の将来予測では、予測値ｙ＾（ｗｃａｎｄ＋Δ），ｙ＾（ｗｃａｎｄ−Δ）の計算だけでなく、評価関数Ｊ（ｗｃａｎｄ）の勾配∇（ｗｃａｎｄ）の計算も行われている。また、ステップＳ１２の集積点判断においては、上述した集積点に関するリストが参照されている。そして、集積点に到達したと判断された場合（ステップＳ１２においてＹｅｓの場合）、それまでに計算した勾配∇の評価が行われる（ステップＳ１４）。具体的には、集積点到達前に計算された勾配∇のうち最小のものが選択される。なお、選択された勾配∇（ｗｃａｎｄ）に対応する修正目標値候補ｗｃａｎｄが、最終的な修正目標値ｗとして決定される。
【００３２】
以上、本実施の形態によれば、集積点到達後の将来予測をストップして余分な演算を続けることがなくなるので、制御性能は同等でも、演算負荷が過大になることを抑制できる。従って、上述したリファレンスガバナを他のリファレンスガバナと併用する場合に、当該他のリファレンスガバナでの計算に演算資源を多く分配できるといった、演算資源の有効活用を図ることができる。
【００３３】
なお、本発明は上述の実施の形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲で種々変形して実施することができる。例えば、図２に示した目標値追従制御構造は、ディーゼルエンジンが低圧ループＥＧＲ装置と高圧ループＥＧＲ装置とを備える場合には、図９の（ａ）〜（ｄ）に示すような制御入力と制御出力との組み合わせにも適用することができる。図９の（ａ）および（ｂ）では、可変ノズル開度（ＶＮ開度）やディーゼルスロットル開度（Ｄスロ開度）の他に、低圧ループＥＧＲ装置のＥＧＲ弁開度（ＬＰＬ−ＥＧＲ弁開度）と高圧ループＥＧＲ装置のＥＧＲ弁開度（ＨＰＬ−ＥＧＲ弁開度）とが制御入力に含まれている。図９の（ｃ）および（ｄ）では、ＥＧＲ率の代わりに、低圧ループＥＧＲ装置のＥＧＲ量（ＬＰ−ＥＧＲ量）と高圧ループＥＧＲ装置のＥＧＲ量（ＨＰ−ＥＧＲ量）とが制御出力に含まれている。
【符号の説明】
【００３４】
２  エンジン本体
１６  可変ノズル
２４  ディーゼルスロットル
３２  ＥＧＲ弁
４０  ＥＣＵ
【図１】
2016130877.tif 000006
【図２】
2016130877.tif 000007
【図３】
2016130877.tif 000008
【図４】
2016130877.tif 000009
【図５】
2016130877.tif 000010
【図６】
2016130877.tif 000011
【図７】
2016130877.tif 000012
【図８】
2016130877.tif 000013
【図９】
2016130877.tif 000014