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