(57)【特許請求の範囲】
【請求項１】
第１の多変数プロセスモデルを用いてプロセス制御を実行するプロセスコントローラにおいて用いるために二つ以上の単一入力単一出力（ＳＩＳＯ）モデルからなる第１の多変数プロセスモデルを制御対象であるプロセスに適応させる方法であって、
前記第１の多変数プロセスモデルに含まれる前記ＳＩＳＯモデルの中から適応させるＳＩＳＯモデルのサブセットを選択することと、
選択された前記ＳＩＳＯモデルのサブセットの各々を別々に適応させることと、
適応された前記ＳＩＳＯモデルを含む第２の多変数プロセスモデルを生成し、前記第２の多変数プロセスモデルを前記第１の多変数プロセスモデル又はその部分に代えて用いるために前記プロセスコントローラに供与することとを含み、
前記ＳＩＳＯモデルのサブセットを選択することが、
前記プロセスへの複数のプロセス入力の各々に対して、測定されたプロセス出力値と前記第１の多変数プロセスモデルにより作成されるプロセス出力値との間の誤差評価値を求め、該誤差評価値と前記プロセス入力との間の相関評価値を決定することと、
該相関評価値に基づいて前記複数のプロセス入力のうちの一つを選択することと、
選択された前記プロセス入力のうちの一つを前記測定された前記プロセス出力に関連付けするＳＩＳＯモデルを、前記ＳＩＳＯモデルのサブセットのうちの一つとして、前記第１の多変数プロセスモデルに含まれる前記ＳＩＳＯモデルの中から選択することとを含む、
方法。
【請求項２】
前記ＳＩＳＯモデルのサブセットを選択する前に、前記プロセスのプロセス入力またはプロセス出力に基づいて、モデル適応を実行する時期を決定することをさらに含む、請求項１記載の多変数プロセスモデル適応方法。
【請求項３】
前記モデル適応を実行する時期を決定することが、
前記プロセス入力または前記プロセス出力を表すデータを収集および格納することと、
前記プロセス入力または前記プロセス出力の変化を検出して、前記変化が所定の量を越えるときに前記モデル適応を実行する時期と決定すべく収集された前記データを解析することとを含む、請求項２記載の多変数プロセスモデル適応方法。
【請求項４】
前記相関評価値に基づいて前記複数のプロセス入力のうちの一つを選択することが、
前記プロセス入力のうちの第一の入力が所定の量の変化を受けているか否かを判定することと、前記プロセス入力のうちの第一の入力が所定の量の変化を受けていると判定されたときに、選択された前記複数のプロセス入力のうちの一つとして前記プロセス入力の第一の入力を選択することとを含む、請求項１記載の多変数プロセスモデル適応方法。
【請求項５】
前記相関評価値に基づいて前記複数のプロセス入力のうちの一つを選択することが、
前記誤差評価値といずれのプロセス入力が最も強い相関関係を有するかを決定することと、
前記相関関係に関する該決定に基づいて、最も強い相関関係を有する前記プロセス入力と関連する前記ＳＩＳＯモデルのうちの一または複数を選択することとを含む、請求項１記載の多変数プロセスモデル適応方法。
【請求項６】
選択された前記ＳＩＳＯモデルのサブセットの各々を別々に適応させることが、
選択された前記ＳＩＳＯモデルのサブセットのうちの少なくとも一つに対してモデルスイッチング適応技法を実行することを含む、請求項１記載の多変数プロセスモデル適応方法。
【請求項７】
選択された前記ＳＩＳＯモデルのサブセットの各々を別々に適応させることが、
選択された前記ＳＩＳＯモデルのサブセットのうちの少なくとも一つに対して属性補間適応技法を実行することを含む、請求項１記載の多変数プロセスモデル適応方法。
【請求項８】
前記第２の多変数プロセスモデルを前記プロセスコントローラに提供するまえに、
適応された前記ＳＩＳＯモデルのサブセットを有する前記第２の多変数プロセスモデルの妥当性を検査することをさらに含む、請求項１記載の多変数プロセスモデル適応方法。
【請求項９】
適応された前記ＳＩＳＯモデルのサブセットを有する前記第２の多変数プロセスモデルの妥当性を検査することが、
適応された前記ＳＩＳＯモデルのサブセットを有する前記第２の多変数プロセスモデルが、適応された前記ＳＩＳＯモデルのサブセットを有しない前記第１の多変数プロセスモデルよりも小さなモデリング誤差を有するか否かを決定することを含む、請求項８記載の多変数プロセスモデル適応方法。
【請求項１０】
適応された前記ＳＩＳＯモデルのサブセットを有する前記第２の多変数プロセスモデルの妥当性を検査することが、
適応された前記ＳＩＳＯモデルのサブセットを有する前記第２の多変数プロセスモデルが同一のプロセス入力データおよびプロセス出力データに対して作用して、適応された前記ＳＩＳＯモデルのサブセットを有しない前記第１の多変数プロセスモデルよりもさらに優れた制御を行うか否か決定することを含む、請求項８記載の多変数プロセスモデル適応方法。
【請求項１１】
選択された前記ＳＩＳＯモデルのサブセットの各々を別々に適応させることが、
パラメータに基づくＳＩＳＯモデルを適応させることを含む、請求項１記載の多変数プロセスモデル適応方法。
【請求項１２】
前記パラメータに基づくＳＩＳＯモデルが、一次遅れ＋むだ時間モデルである、請求項１１記載の多変数プロセスモデル適応方法。
【請求項１３】
選択された前記ＳＩＳＯモデルのサブセットの各々を別々に適応させることが、
ノンパラメトリックベースのＳＩＳＯモデルを適応させることを含む、請求項１記載の多変数プロセスモデル適応方法。
【請求項１４】
前記ノンパラメトリックベースのＳＩＳＯモデルが、ステップ応答モデルおよびインパルス応答モデルのうちの一つである、請求項１３記載の多変数プロセスモデル適応方法。
【請求項１５】
選択された前記ＳＩＳＯモデルのサブセットの各々を別々に適応させることが、
前記ノンパラメトリックベースのＳＩＳＯモデルの属性を調整することを含む、請求項１３記載の多変数プロセスモデル適応方法。
【請求項１６】
前記ノンパラメトリックベースのＳＩＳＯモデルの属性を調整することが、
前記ノンパラメトリックベースＳＩＳＯモデルの、
（１）応答が最初に観察されるまでの時間、
（２）すべてのステップ応答係数における各１以上のステップ係数値、および
（３）応答曲線の傾斜
のうちの一つを調整することを含む、請求項１５記載の多変数プロセスモデル適応方法。
【請求項１７】
プロセスの制御に用いられるプロセス制御システムであって、
前記プロセスの制御に用いられる二つ以上の単一入力単一出力（ＳＩＳＯ）モデルからなる第１の多変数プロセスモデルを有するプロセスコントローラと、
前記プロセスコントローラに通信可能に接続されているモデル適応ユニットとを備えており、
該モデル適応ユニットが、
前記第１の多変数プロセスモデルに含まれる前記ＳＩＳＯモデルの中から適応させる前記ＳＩＳＯモデルのサブセットを選択するように構成された第一のユニットと、
選択された前記ＳＩＳＯモデルのサブセットの各々を変更するように構成された第二のユニットと、
変更された前記ＳＩＳＯモデルを含む第２の多変数プロセスモデルを生成するように構成された第三のユニットと、
前記第１の多変数プロセスモデル又はその部分に代えて用いるべく、前記第２の多変数プロセスモデルを前記プロセスコントローラに供与するように構成された第四のユニットとを有してなり、
前記第一のユニットが、複数のプロセス入力の各々に対して、測定されたプロセス出力と前記第１の多変数プロセスモデルにより作成されるプロセス出力との間の誤差評価値を求め、該誤差評価値と前記プロセス入力との間の相関評価値を決定するとともに、
該相関評価値を利用して前記複数のプロセス入力のうちの一つを選択し、
選択された前記プロセス入力のうちの一つを前記測定されたプロセス出力に関連付けするＳＩＳＯモデルを、前記ＳＩＳＯモデルのサブセットのうちの一つとして、前記第１の多変数プロセスモデルに含まれる前記ＳＩＳＯモデルのうちから選択するように構成されている、プロセス制御システム。
【請求項１８】
前記モデル適応ユニットが、適応サイクルを開始する時期を決定すべく、前記プロセスの一または複数のプロセス入力またはプロセス出力の変化を検出するように構成されたスーパバイザユニットをさらに備えてなる、請求項１７記載のプロセス制御システム。
【請求項１９】
前記スーパバイザユニットが、データ収集ユニットを備えているとともに、該データ収集ユニットに前記一または複数のプロセス入力またはプロセス出力のうちの一つを表すデータを収集および格納し、
前記一または複数のプロセス入力またはプロセス出力のうちの一つの変化が所定の量を越えたときに前記適応サイクルを開始する時期と決定するために、収集された前記データを解析するように構成されている、請求項１８記載のプロセス制御システム
【請求項２０】
前記第二のユニットが、選択された前記ＳＩＳＯモデルのサブセットのうちの少なくとも一つに対してモデルスイッチングモデル適応技法を実行してなる、請求項１７記載のプロセス制御システム。
【請求項２１】
前記第二のユニットが、選択された前記ＳＩＳＯモデルのサブセットのうちの少なくとも一つに対して属性補間モデル適応技法を実行してなる、請求項１７記載のプロセス制御システム。
【請求項２２】
前記モデル適応ユニットは、前記第四のユニットが前記第２の多変数プロセスモデルを前記プロセスコントローラに提供するまえに、変更された前記ＳＩＳＯモデルのサブセットを用いて前記第２の多変数プロセスモデルを検査するように構成された妥当性検査ユニットをさらに備えてなる、請求項１７記載のプロセス制御システム。
【請求項２３】
選択された前記ＳＩＳＯモデルのサブセットのうちの少なくとも一つが、パラメータに基づくＳＩＳＯモデルである、請求項１７記載のプロセス制御システム。
【請求項２４】
前記パラメータに基づくＳＩＳＯモデルが、一次遅れ＋むだ時間モデルである、請求項２３記載のプロセス制御システム。
【請求項２５】
選択された前記ＳＩＳＯモデルのサブセットのうちの少なくとも一つが、ノンパラメトリックベースのＳＩＳＯモデルである、請求項１７記載のプロセス制御システム。
【請求項２６】
前記ノンパラメトリックベースのＳＩＳＯモデルが、ステップ応答モデルおよびインパルス応答モデルのうちの一つである、請求項２５記載のプロセス制御システム。
【請求項２７】
前記第二のユニットが、前記ノンパラメトリックベースのＳＩＳＯモデルの、
（１）応答が最初に観察されるまでの時間、
（２）すべてのステップ応答係数における各１以上のステップ係数値、および
（３）応答曲線の傾斜
のうちの一つを調整することによって、前記ノンパラメトリックベースのＳＩＳＯモデルの属性を調整するように構成されている、請求項２５記載のプロセス制御システム。
【請求項２８】
前記プロセスコントローラが、モデル予測制御型コントローラである、請求項１７記載のプロセス制御システム。
【請求項２９】
二つ以上の単一入力単一出力（ＳＩＳＯ）モデルからなる第１の多変数プロセスモデルを有するプロセス制御システムにおいて用いられるモデル適応ユニットであって、
コンピュータ読み取り可能媒体と、
プロセッサ上で実行されるように構成されかつ前記コンピュータ読み取り可能媒体に格納されているプログラムとを備えており、
該プログラムが、
前記第１の多変数プロセスモデルに含まれる前記ＳＩＳＯモデルの中から適応させる前記ＳＩＳＯモデルのサブセットを選択するように構成された第一のルーチンと、
選択された前記ＳＩＳＯモデルのサブセットの各々を適応させるように構成された第二のルーチンと、
適応された前記ＳＩＳＯモデルを含む第２の多変数プロセスモデルを生成するように構成された第三のルーチンと
前記第２の多変数プロセスモデルを前記第１の多変数プロセスモデル又はその部分に代えて用いるために前記プロセス制御システムに供与するように構成された第四のルーチンとを有しており、
前記第一のルーチンが、前記プロセスへの複数のプロセス入力の各々に対して、測定されたプロセス出力と前記第１の多変数プロセスモデルにより作成されるプロセス出力との間の誤差評価値を求め、該誤差評価値と前記プロセス入力との間の相関評価値を決定し、
該相関評価値を用いて前記複数のプロセス入力のうちの一つを選択し、
選択された前記プロセス入力のうちの一つを前記測定されたプロセス出力に関連付けするＳＩＳＯモデルを、前記ＳＩＳＯモデルのサブセットのうちの一つとして、前記第１の多変数プロセスモデルに含まれる前記ＳＩＳＯモデルのうちから選択する、モデル適応ユニット。
【請求項３０】
前記プログラムが、適応サイクルを開始する時期を決定すべく、前記プロセス制御システムの制御対象であるプロセスのプロセス入力またはプロセス出力の変化を検出するように構成されたスーパバイザルーチンをさらに備えてなる、請求項２９記載のモデル適応ユニット。
【請求項３１】
前記第二のルーチンが、選択された前記ＳＩＳＯモデルのサブセットのうちの少なくとも一つに対してモデルスイッチング適応技法を実行してなる、請求項２９記載のモデル適応ユニット。
【請求項３２】
前記第二のルーチンが、選択された前記ＳＩＳＯモデルのサブセットのうちの少なくとも一つに対して属性補間適応技法を実行してなる、請求項２９記載のモデル適応ユニット。
【請求項３３】
前記プログラムは、前記第四のルーチンが、前記第２の多変数プロセスモデルを前記プロセス制御システムに提供するまえに、適応された前記ＳＩＳＯモデルのサブセットを用いて前記第２の多変数プロセスモデルを検査する妥当性検査ルーチンをさらに備えてなる、請求項２９記載のモデル適応ユニット。
【請求項３４】
選択された前記ＳＩＳＯモデルのサブセットのうちの少なくとも一つが、ノンパラメトリックベースのＳＩＳＯモデルであり、前記第二のルーチンが、前記ノンパラメトリックベースのＳＩＳＯモデルの、
（１）応答が最初に観察されるまでの時間、
（２）すべてのステップ応答係数における各１以上のステップ係数値、および
（３）応答曲線の傾斜
のうちの一つを調整することによって前記ノンパラメトリックベースのＳＩＳＯモデルの属性を調整してなる、請求項２９記載のモデル適応ユニット。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明にかかる方法および装置は、一般的にプロセス制御技術に関するものであり、さらに詳細には、モデル適応プロシージャ中にモデルスイッチングおよびパラメータ補間のうちの一方または両方を実行する適応型多変数制御技術に関するものである。
【背景技術】
【０００２】
本発明は、ここで参照することによりその全内容を明示的に援用する、２００３年４月２１日に出願された、表題が「適応フィードバック／フィードフォワードＰＩＤコントローラ（Ａｄａｐｔｉｖｅ  Ｆｅｅｄｂａｃｋ／Ｆｅｅｄｆｏｒｗａｒｄ  ＰＩＤ  Ｃｏｎｔｒｏｌｌｅｒ）」である米国特許出願番号第10/419,582号に関連している。また、この米国特許出願番号第10/419,582号は、ここで参照することによりその全内容を明示的に援用する、２０００年６月２０日に出願され２００３年６月１０日に米国特許第6,577,908B1号 として特許された、表題が「適応フィードバック／フィードフォワードＰＩＤコントローラ（Ａｄａｐｔｉｖｅ  Ｆｅｅｄｂａｃｋ／Ｆｅｅｄｆｏｒｗａｒｄ  ＰＩＤ  Ｃｏｎｔｒｏｌｌｅｒ）」である米国特許出願番号第09/597,611号の一部継続出願である。
【０００３】
大規模な製造工場および化学精製所などの如き自動化されたシステムにおいて適応プロセス制御を実現するためにロジックに基づくコントローラスイッチング戦略を用いることが当該技術分野において知られている。ロジックに基づくスイッチング戦略の説明の一例として、たとえば、モールス（Ｍｏｒｓｅ）、エフ・エム・ペイト（Ｆ．Ｍ．Ｐａｉｔ），およびエス・アール・ウェラ（Ｓ．Ｒ．Ｗｅｌｌｅｒ）、「自己調整制御のためのロジックに基づくスイッチング戦略（Ｌｏｇｉｃ−Ｂａｓｅｄ  Ｓｗｉｔｃｈｉｎｇ  Ｓｔｒａｔｅｇｉｅｓ  ｆｏｒ  Ｓｅｌｆ−Ａｄｊｕｓｔｉｎｇ  Ｃｏｎｔｒｏｌ）」、決定および制御に関する第３３回ＩＥＥＥ会議（ＩＥＥＥ  ３３rd  Ｃｏｎｆｅｒｅｎｃｅ  ｏｎ  Ｄｅｃｉｓｉｏｎ  ａｎｄ  Ｃｏｎｔｒｏｌ）（１９９４年１２月）が挙げられる。一般的にいえば、ロジックに基づくコントローラ−スイッチング戦略は、プリルーテッドコントローラによるアプローチとモデルに基づくパラメータ化コントローラによるアプローチとして一般的に区別されている二つのアプローチのうちの一つに分類されうる。
【０００４】
プリルーテッドコントローラチューニングとは、原則的に、事前定義された一組の有力なコントローラに含まれている有力なコントローラを評価することである。この評価は、満足のいく機能を有するコントローラが特定されると終了する。プリルーテッドコントローラチューニングシステムは設計が容易であり、コントローラの構造に関する要求項目がほとんどない。しかしながら、プリルーテッドコントローラチューニングシステムのこのような利点も、チューニング時間に関する性能が劣悪であるために、見劣りがするものとなっている。すなわち、事前定義された一組の有力なコントローラからの最適なコントローラの選択に必要な時間が非常に長い。
【０００５】
一般的に、モデルに基づくパラメータ化コントローラは、二つ以上のパラメータ依存サブシステムと、出力推定誤差を生成するモデルと、内部コントローラとから構成されている。動作状態では、適切に定義された一組のモデルによる推定に基づく制御信号が、被制御プロセスに伝達される。モデルに基づくパラメータ化コントローラにより、「サイクリックスイッチング」の概念に基づいてコントローラスイッチング戦略が具象化される。サイクリックスイッチングは、プロセスに対するさらなる励起信号の提供の有無に関係なく用いられうる。
【０００６】
プロセス適応制御に対するサイクリックスイッチングアプローチに関する有益な説明として、ケイ・エス・ナレンドラ（Ｋ．Ｓ．Ｎａｒｅｎｄｒａ）およびジェイ・バラクリシュナン（Ｊ．Ｂａｌａｋｒｉｓｈｎａｎ）、「複数モデルを用いた適応制御（Ａｄａｐｔｉｖｅ  Ｃｏｎｔｒｏｌ  Ｕｓｉｎｇ  Ｍｕｌｔｉｐｌｅ  Ｍｏｄｅｌｓ）」、自動制御に関するＩＥＥＥ報告（ＩＥＥＥ  Ｔｒａｎｓａｃｔｉｏｎｓ  ｏｎ  Ａｕｔｏｍａｔｉｃ  Ｃｏｎｔｒｏｌ）、第４２巻、第２号、１７７〜１８７ページ（１９９７年２月）が挙げられる。この文献は、複数のパラメータと、当該複数のパラメータに対応するモデルパラメータを有しかつ並列に動作するＮ個のモデルとによって特徴付けされるコントローラを備えたプロセス制御システムを開示している。いずれの時点においても、スイッチングルールにより、一つのモデルとそれに対応するパラメータ化コントローラとが選択され、それに対応する制御入力がプロセスの制御に用いられる。これらのモデルは、プロセス要件、オペレータの必要要件、およびその他の適切な考慮すべき事項に応じて、固定パラメトリックモデルであってもよいし、適応パラメトリックモデルであってもよい。固定パラメトリックモデル制御システムは、未知のプロセスのパラメータに十分に近いパラメータにより特徴付けられる少なくとも一つのモデルの存在を担保する簡易かつ有効な手段を提供する。
【０００７】
固定パラメトリックモデルを用いる、サイクリックスイッチングに基づくプロセス制御システムは、高速な適応速度を提供するが、プロセスコントローラ内に非常に多くのモデルを設計・格納する必要がある。なお、固定モデルは、有限の数のプロセス環境または状態のみを正確に表現できる。また固定モデルは、プロセス精度を漸近的に向上させるために、適応モデルを用いる必要がある。
【０００８】
現実的にいえば、モデルに基づくスイッチング戦略は、妥当なプロセス近似が非常に多くのモデルを必要とすることに起因する複数の問題をもたらす。たとえば、固定モデル型自己チューナを含む簡易な単一入力単一出力（ＳＩＳＯ）システムは、十分なプロセス性能を達成するには何百もの固定モデルを含むことが当然予想される。したがって、システムがさらに複雑になると、たとえば多変数システムになると、カスタマイズされる固定モデルの必要数は指数関数的に増えることになり、それにより、システム設定時間およびシステム必要メモリが増大する。より効果的に解決するには、特定のモデル構造およびコントローラタイプを検討する必要があり、簡易スイッチング戦略をさらに改良されたプロシージャに切り替えることが考えられる。
【０００９】
ダーリン（Ｄａｈｌｉｎ）コントローラに対する改良されたモデルに基づくアプローチが、ゲンドロン（Ｇｅｎｄｒｏｎ）、「原因不明のまたは変動する遅延を有するプラントのデッドタイム補償器（Ｉｍｐｒｏｖｉｎｇ  ｔｈｅ  Ｒｏｂｕｓｔｎｅｓｓ  ｏｆ  Ｄｅａｄ−Ｔｉｍｅ  Ｃｏｍｐｅｎｓａｔｉｏｎ  ｆｏｒ  Ｐｌａｎｔｓ  ｗｉｔｈ  Ｕｎｋｎｏｗｎ  ｏｒ  Ｖａｒｙｉｎｇ  Ｄｅｌａｙ）」、制御システム会議９０（Ｃｏｎｔｒｏｌ  Ｓｙｓｔｅｍ  ９０  Ｃｏｎｆｅｒｅｎｃｅ）（１９９０年、ヘルシンキ）、において提案されている：。この文献は、デッドタイム変化に基づくプロセス適応を実現する簡易な一次遅れ＋むだ時間プロセスモデルを開示している。簡易なモデルスイッチングに依存するというより、むしろ、上記のコントローラは、さまざまなデッドタイムにより特徴付けされる一組のモデルの重み付け合計に基づいたプロセスモデルを用いる。この一組のモデルの各々は、プロセス出力値の予測値を生成すべくなしてあり、その対応する重みは、予測誤差の関数の単関数として、自動的に調整される。この基本概念が、拡張されて、ダーリンコントローラ構成内のプロセスゲインおよびデッドタイムの変動を含むようになってきている。
【００１０】
一般的に、適応ＰＩＤコントローラの設計のための有力なアプローチは、直接アプローチと、間接アプローチまたはモデルに基づくアプローチとである。上述のように、モデルに基づくアプローチは、スイッチング戦略を用いる制御システムに対して有効であり、かつ、適応スイッチングＰＩＤコントローラを設計するための適切な出発点を提供する。また、モデルに基づくアプローチは、モデルパラメータの変動を追跡するための、逐次最小自乗（ＲＬＳ）推定器と結合された、モデル同定に基づく適応ＰＩＤコントローラを提供することが知られている。逐次モデル同定器に付随する典型的な問題としては、初期パラメータの選定の難しさ、励起の不十分さ、フィルタリング、パラメータワインドアップ、パラメータ追跡速度の緩慢さが挙げられる。これらの変数が複雑であるとともに正確な推定量の計算が困難であるので、公知の、モデルに基づく適応ＰＩＤコントローラは、プロセスモデルを簡略化することにより、その性能を向上しうるということが当該技術分野においてよく知られている。
【００１１】
簡略化された、モデルに基づく適応コントローラの一例として、がアストロム（Ａｓｔｒｏｍ）およびハグランド（Ｈａｇｇｌｕｎｄ）、「周波数応答技術に基づく工業用適応コントローラ（Ｉｎｄｕｓｔｒｉａｌ  Ａｄａｐｔｉｖｅ  Ｃｏｎｔｒｏｌｌｅｒｓ  Ｂａｓｅｄ  ｏｎ  Ｆｒｅｑｕｅｎｃｙ  Ｒｅｓｐｏｎｓｅ  Ｔｅｃｈｎｉｑｕｅｓ）」、オートマチカ（Ａｕｔｏｍａｔｉｃａ）、第２７巻、第４号、５９９〜６０９ページ（１９９１年）を挙げることができる。概して、この文献は、周波数ドメインにおいてプロセスモデル適応を実行し、設定値の変更および、自然現象による外乱に応答してチューニングを行うように設計されたコントローラを開示している。さらに詳細には、チューニング周波数は、プロセスの入出力に対して帯域通過フィルタを適用することにより選択され、そのフィルタの通過帯域の周波数は、オートチューナ（オンデマンド型チューナ）により定義される。このオートチューナは、適応チューナの動作に先立ってリレイ振動技術を用いて最終周期を定義し、簡易化されたＲＬＳ推定器を用いてそのチューニング周波数に対するプロセスゲインを定義する。このオートチューナは、プロセスゲインの変化を追跡する機能を有している。しかしながら、デッドタイムまたは時定数に変化が発生すると、追跡されたポイントは、もはや、−π位相を示しておらず、コントローラによるチューニングが不正確なものとなる。
【００１２】
さらに、複数のチューニング周波数を適用し、インターポレータを用いて−π位相を有する周波数を確定することにより、チューニングを向上させることが知られている。あるいは、単一のチューニング周波数を適用し、各チューニングサイクルのあと、−π位相に追従するために周波数を調整することが可能である。両方のチューナモデルは、その後の設定値変更および自然現象による外乱に対応し、コントローラ出力部または設定値入力部において外部励起を入力しうる。このようなオートチューナは、従来技術におけるようなサイズ制約およびセットアップ制約を有していないが、従来技術よりも非常に複雑である。
【００１３】
さらに、両方のチューナモデルは、最終ゲインおよび最終周期の二つのパラメータのみしか認識しない初期の適応モデルを用いている。このような簡易な２パラメータ適応モデルを組み込んでいるチューナモデルは、Ｚｉｅｇｌｅｒ−Ｎｉｃｈｏｌｓチューニングまたはその改良バージョンに対して適切であるが、内部モデル制御（ＩＭＣ）またはＬ
ａｍｂｄａチューニングが好ましい用途においては不適切である。ＲＬＳ同定器をフィードフォワード制御の静的ゲインの決定に用いてもよいが、このＲＬＳ同定器アプローチは、適切なフィードフォワード制御に必要なプロセスフィードフォワードダイナミックスに欠けている。さらに、このアプローチは、フィードフォワード信号が負荷外乱であり、また、摂動信号をフィードバック経路に入力することができないため、励起が十分でないという問題を有している。
【００１４】
フィードフォワード適応に対する他の解決策が、ブリストルおよびハンセンにより、表題が「多変数適応フィードフォワードコントローラ（Ｍｕｌｔｉｖａｒｉａｂｌｅ  Ａｄａｐｔｉｖｅ  Ｆｅｅｄｆｏｒｗａｒｄ  Ｃｏｎｔｒｏｌｌｅｒ）」である米国特許第5,043,863号に開示されている。この特許は、外乱負荷を含むように設計された微分方程式に基づくモデルを開示している。このプロセスモデルは、測定されたプロセスデータに基づいて周期的に更新される。ここで、外乱は、プロジェクション法によりもたらされるモーメント関係およびコントロール関係により特徴づけされる。一般的に、導かれる解決策は、非常に複雑であり、かなりの励振を必要とし、上述のＲＬＳ同定器によるアプローチとほとんど同じである。さらに、導かれる解決策は、フィードフォワード制御に対してのみ適切であり、フィードバック付き適応コントローラに対して適用不能である。
【００１５】
米国特許第6,577,908号には、適応スイッチングプロシージャを大幅に修正、改良したものが開示されている。この適応スイッチングプロシージャは、モデル補間に代えて、パラメータ補間・モデル位置調整技術を導入しており、この技術は、適応に用いられるモデルの数を飛躍的に減少させることを可能にする。この特許に開示されているコントローラは、適応時間の短縮、プロセスモデル評価の完全性、およびプロセス励起に関する要求項目の削減という特徴を有している。
【００１６】
今日まで、以上で説明したモデルおよびパラメータの適応およびスイッチング技術は、単純なプロセスモデルには適するが、単一入力単一出力（ＳＩＳＯ）コントローラシステムを適応するためだけの利用に留まっており、コントローラが一または複数の入力に基づいて二つ以上の変数を同時に制御するような多変数制御状況に対して適用されてこなかった。しかしながら、多変数制御、とくに多変数制御の一つの特定の形式であるモデル予測制御（ＭＰＣ）は、急速に、複雑な制御状況に対する有用な制御戦略となってきている。したがって、多変数コントローラに対して適用することができ、これらのコントローラがランタイム中に適応することを可能とし、これにより、プロセスや設定値などの変更に対する反応を向上させ、さらに優れた制御を実現させるような適応戦略を提供することが望まれている。
【発明の開示】
【課題を解決するための手段】
【００１７】
適応ＭＰＣプロセス制御システムの如き適応多変数プロセス制御システムは、二つ以上のＳＩＳＯモデルからなるセットから構成されているかまたは二つ以上のＳＩＳＯモデルからなるセットとして特徴付けされる多変数プロセスモデルを有する多変数プロセスコントローラと、該多変数プロセスモデルを適応させる適応システムとを備えている。一般的にいえば、適応システムは、スーパバイザを有しており、該スーパバイザは、適応サイクルを開始するのに十分なプロセス入力の変化を検出し、このような変化を検出した場合、適応の実行に必要なプロセス入力データおよびプロセス出力データを収集する。次いで、スーパバイザは、前記多変数プロセスモデル内のＳＩＳＯモデルのうちの適応させるべきサブセットを決定し、それによって、前記多変数プロセスモデルを適応させうる。このようなＳＩＳＯモデルの選択は、各プロセス出力に対して、どのプロセス入力が、実際の（測定された）プロセス出力と多変数プロセスモデルにより作成されたプロセス出力との間の誤差と最も強い相関関係を有しているかを判定し、これらの最も強い相関関係を有する入力をモデリング誤差を有するプロセス出力と関連づけるＳＩＳＯモデルを選択することによって行われうる。所望ならば、適応には、変化が所定の量である入力と関連するＳＩＳＯモデルのみが選択されてもよい。
【００１８】
次に、適応システムは、標準的なまたは公知のモデルスイッチング技法およびパラメータ補間技法を実行し、選択されたＳＩＳＯモデルの各々を適応させる。ＳＩＳＯモデルのうちの一または複数が適応されたあと、結果として得られた多変数プロセスモデルは、当該多変数プロセスモデル（すなわち、適応されたＳＩＳＯモデルおよび適応されていないＳＩＳＯモデルの両方を有するプロセスモデル）が、同一のプロセス入力データおよびプロセス出力データに対して動作して現在のまたは適応されていない多変数プロセスモデルよりもさらに優れた制御を行うか否か、および／または所定の量だけモデル誤差を削減するか否かを判定することによりその妥当性が検査される。適応された多変数プロセスモデルが現在の多変数プロセスモデルよりも性能が優れている場合、当該適応された多変数プロセスモデル（または、それに付随する適応されたＳＩＳＯモデル）は、多変数コントローラにより用いられる形式に変更され、次いで、プロセス制御において多変数プロセス制御システムで用いられるべく提供される。
【００１９】
もちろん、適応プロシージャは、十分な量のプロセス入力またはプロセス出力の変化が検出される毎に、プロセス入力またはプロセス出力が強制的に変更されたあとなどの如き必要または所望に応じて、繰り返されてもよい。さらに、プロセスの複雑さに応じておよび充足させるべき演算要求事項に応じて適応プロシージャの速度を変更することができるようにするために、いかなる数のＳＩＳＯモデルをいかなる特定の適応プロシージャ中に適応してもよい。
【００２０】
本明細書に記載の多変数適応プロシージャは、一次遅れ＋むだ時間モデルの如きパラメータに基づくＳＩＳＯモデルを用いて、またはステップ応答モデルもしくはインパルス応答モデルの如きノンパラメトリック型モデルを用いて適用されてうる。後の方のケースの場合、適応プロシージャは、選択されたプロセスモデルの各々を進展または適応させるために、モデルスイッチングおよびモデル属性補間（モデルパラメータ補間はモデル属性補間のサブセットである）を用いうる。このケースの場合、パラメータに基づくモデルのモデルゲイン、時定数、およびデッドタイムの如き一または複数のパラメータを変更する代わりに、ノンパラメトリックモデルの応答が最初に観察されるまでの時間、等級、勾配などの如きノンパラメトリックモデルの属性が変更され、新規のまたは適応されたノンパラメトリックモデルを決定するために、これらの属性の変更に基づいて補間が実行されうる。
【発明を実施するための最良の形態】
【００２１】
図１は、プロセス10において多変数制御を提供すべくモデル予測制御（ＭＰＣ）を用いる適応多変数プロセス制御システム8の一例を示している。この制御システム8は、ここで参照することによりその全内容を明示的に援用するところの米国特許第6,577,908号に開示されているフィードバック／フィードフォワード（ＦＢ／ＦＣ）ＰＩＤ制御システムといくぶん設計が類似している。事実、制御システム8は、以下でさらに詳細に説明するように、米国特許第6,577,908号に開示されているＰＩＤ制御システムコンポーネントを用いうる。
【００２２】
図１に示されているように、適応多変数プロセス制御システム8は、標準的なまたは典型的なＭＰＣシステム12を備えており、当該ＭＰＣシステム12は、プロセス10の制御に用いられる操作変数（ＭＶ）信号を生成するＭＰＣコントローラに接続されたプロセスモデルブロック14有している。一般的にいえば、プロセスモデルブロック14は、一組の単一入力単一出力（ＳＩＳＯ）モデルからなりうる多変数プロセスモデル14aを有しており、ＳＩＳＯモデルの各々は、一次遅れ＋むだ時間モデルの如きパラメトリックモデルまたはステップ応答曲線の如きノンパラメトリックモデルでありうる。また、ＳＩＳＯモデルの各々は、制御変数（ＣＶ）または補助変数（ＡＶ）とも呼ばれるプロセス出力のうちの一つに関して、ＭＰＣコントローラ16により生成されるＣＶおよび測定される外乱変数（ＤＶ）、すなわちプロセス入力のうちの一つに対する関係または応答を定義する。動作状態では、プロセスモデルブロック14は、ＭＶおよび測定されたＤＶを受信し、その内部の多変数プロセスモデル14aを用いて、プロセス10の予測される未来応答値（未来ＣＶ値および未来ＡＶ値）を規定する予測されるまたは未来の出力軌跡信号を生成する。この予測された出力軌跡信号はベクトル加算器18へ伝達される。
【００２３】
また、ＭＰＣ制御システム12は、設定値予測ユニットまたは設定値予測ブロック20を備えている。この設定値予測ユニット20は、プロセス出力（ＣＶおよびＡＶ）のうちの一または複数の各々に対する個別の設定値信号でありうるプロセス設定値信号ＳＰを受信すべくなしてあるとともに、予測される未来設定値を規定する設定値軌跡信号を生成すべくなしてある。また、たとえば、設定値予測ユニット20は、プロセス10に対する未来設定値の変更を示すオプチマイザ（図示せず）に接続されうる。また、設定値予測ユニット20により生成される設定値軌跡信号もベクトル加算器18に伝達され、当該ベクトル加算器は、予測されたプロセス出力軌跡信号と予測された設定値軌跡信号との間の差を定義する誤差軌跡信号を演算する。この誤差軌跡信号は、ＭＰＣコントローラ16の入力部に伝達される。ＭＰＣコントローラ16は、公知の方法でこの誤差軌跡信号を用いて操作変数（ＭＶ）を生成する。次いで、これらの操作変数が、プロセス10の入力部に伝達され、それにより、プロセス10が制御される。
【００２４】
通例のように、制御変数（ＣＶ）および所望ならば補助変数（ＡＶ）を含むプロセス10の出力値は、加算器22に戻され、そこで同一期間にわたって予測される（多変数プロセスモデル14aにより生成される）モデル出力値と組み合わされ、現在のモデル出力誤差が求められる。次いで、このモデル出力誤差は、プロセスモデルブロック14に戻され、そこで、プロセスモデルブロック14により予測された出力値を補正するために、標準的な方法で処理される。
【００２５】
いうまでもなく、ＭＰＣ制御システム12は、慣用技術であり、当該技術分野において公知のものである。したがって、プロセス10のＭＰＣ制御または他の多変数制御を実現すべく公知の方法でＭＰＣ制御システム12を作成または修正することは当業者にとって明らかなことである。さらに、本明細書記載の適応技術は、ＭＰＣコントローラとの関連において用いられているが、同様に、他のモデルに基づく多変数プロセス制御システムにおいて用いられてもかまわない。
【００２６】
さらに図１を説明すると、モデル適応ネットワーク30を利用することにより、以前までは単一入力単一出力（ＳＩＳＯ）モデルに基づく制御システムでの使用だけに制限されていたモデルスイッチング技法および／またはパラメータ（もしくはモデル属性）補間技法を用いて、ＭＰＣ制御システム12において利用される多変数プロセスモデル14aを適応させうる。
【００２７】
図１に示されているように、モデル適応ネットワーク30は、適応プロセスを制御しうるスーパバイザ32と、多変数プロセスモデル14aに使用するために評価されうる、異なる潜在的な能力を備えた一組のモデルを格納しうるモデルセットブロック34と、このモデルセットブロック34内のモデルを評価し、適応プロセスにおいて用いるべきモデルパラメータまたはモデル属性の重みづけを決定するモデル評価ブロック36とを有している。また、モデル適応ネットワーク30は、多変数プロセスモデル14aのコンポーネントにおいて用いられるべきモデルパラメータ値またはモデル属性値を決定するパラメータ／属性補間ブロック38と、適応されたときでありかつプロセスモデルブロック14に新規の多変数プロセスモデルを代入するまえに、その新規のモデルの動作の妥当性を検査するために用いられうるモデル妥当性検査ブロック40と、ＭＰＣ制御システム12において利用すべくプロセスモデルブロック14により用いられる形態に新規のまたは適応されたプロセスモデルの形態を変更するモデル形式変更ブロック42とを含んでいる。
【００２８】
適応ネットワーク30の動作は、典型的なＭＰＣアプリケーションでプロセス10をモデル化すべくブロック14内の多変数プロセスモデル14aが動作する様式をまず説明することにより最も良く理解されるであろう。このために、図２は、複数の入力ＩＮi（ｉ＝１〜ｊ）と複数の出力ＯＵＴi（ｉ＝１〜ｎ）とを有する多変数プロセスモデル14aの一般化された形態を示している。入力ＩＮiは、操作変数（ＭＶ1、ＭＶ2など）および測定外乱変数（ＤＶ1、ＤＶ2など）でありうる一方、出力ＯＵＴiは、制御変数（ＣＶ1、ＣＶ2など）および補助変数（ＡＶ1、ＡＶ2など）でありうる。
【００２９】
さまざまな入力信号ＩＮi（ｉ＝１〜ｊ）により予測出力信号ＯＵＴi（ｉ＝１〜ｎ）がもたらされるまたは生成される様式をモデル化する一つの技法としては、まず入力／出力の各対に対して異なる単一入力単一出力（ＳＩＳＯ）モデルを定義することにより、各入力信号ＩＮiに対する各出力信号ＯＵＴiの応答を関連づけすることが挙げられる。したがって、図２のプロセスモデルでは、ｎ個のＳＩＳＯモデルのｊ倍の数のモデルが存在することになり、各ＳＩＳＯモデルが、入力信号ＩＮiのうちの一つのみの入力信号に対する出力信号ＯＵＴiのうちの一つの応答を定義する。本明細書で用いられる場合、ＳＩＳＯモデルが、Ｍｏｄｅｌ（Ｉｎ、Ｏｕｔ）で示されており、変数「Ｉｎ」が１からｊ（全入力数）の範囲であり、変数「アウト」が１からｎ（全出力数）の範囲である。したがって、入力信号ＩＮ3に対する出力信号ＯＵＴ2の時間経過応答を定義するＳＩＳＯモデルはＭｏｄｅｌ（３，２）ということになる。
【００３０】
次に、図３に示されているように、ＳＩＳＯモデル、すなわちＩｎ＝１〜ｊの範囲であるＭｏｄｅｌ（Ｉｎ、ｉ）の各々の出力の合計として、特定の出力信号ＯＵＴiが決定されてもよい。換言すれば、多変数線形プロセスをモデル化する一つのアプローチでは、プロセス出力が、プロセス入力の各々とそのプロセス出力との間の関係を定義するＳＩＳＯモデルの出力を重ね合わせたもの（総和）であると仮定する。いうまでもなく、すべてのプロセス出力ＯＵＴi（ｉ＝１〜ｎ）をこの方法でモデル化しうるし、ＯＵＴｉ信号の各々とＩＮｉ信号の各々とが、定常状態に至るまでの時間の如き特定の期間にわたる関係を定義する時間依存信号であってもよい。この結果として、図１の多変数プロセスモデルブロック14は、ある時間範囲にわたるＯＵＴ信号の各々の値を定義する軌跡信号またはベクトルを定義する。
【００３１】
一般的にいえば、図３に示されているＳＩＳＯモデルＭＯＤＥＬ（Ｉｎ，Ｏｕｔ）は、一次遅れ＋むだ時間モデル（または、その他のパラメトリックモデル）の如きパラメータの基づくモデルまたはパラメトリックモデルであってもよく、または対応する入力の変化に応答するある経過時間範囲の出力信号を定義するステップ応答モデルの如きノンパラメトリックモデルであってもよい。もちろん、本明細書記載のモデルに加えてまたはそれらに代えて、他のタイプのパラメトリックモデルおよびノンパラメトリックモデル（たとえば、インパルスモデル、ランプモデルなど）を用いてもよい。
【００３２】
さらに、図１を説明すると、適応システム30は、ＭＰＣコントローラシステム12に用いられるべき新規の多変数プロセスモデルを生成するために、多変数プロセスモデル14aで用いられるＳＩＳＯモデルのうちの一または複数を適応させるように動作する。一般的にいえば、各適応サイクルまたは各適応スキャン中、適応システム30は、プロセス10の（多変数プロセスモデル14aにより作成される）予測プロセス出力値と実際の（測定される）出力値との間の誤差に対する個々のＳＩＳＯモデルの寄与に基づいて、ＳＩＳＯモデルのうちのいずれを適応させるべきかを決定する。次いで、適応システム30は、モデルセット34、モデル妥当性検査ブロック36、およびパラメータ／属性補間ブロック38を用いて、適応させるべく特定された各ＳＩＳＯモデルを中心としてプロセスモデル化を実行する。このモデルセンタリングは  特定されたＳＩＳＯモデルを修正または変更する様式を決定するモデルパラメータ補間技法またはモデル属性補間技法の形態で達成されうる。適応システム30は、このプロシージャ中またはその終了時点において、新規に適応されたＳＩＳＯモデルの妥当性検査をモデル妥当性検査ブロック40を用いて行い、当該適応されたＳＩＳＯモデル（すなわち、当該適応されたＳＩＳＯモデルを組み合わせて得られる適応多変数プロセスモデル）がＭＰＣ制御システム12により現在用いられている多変数プロセスモデル14aよりもよい性能を有することを確認する。
【００３３】
新規のもしくは適応された多変数プロセスモデルまたは当該多変数プロセスモデルを構成する一もしくは複数のＳＩＳＯモデルが、プロセスモデルブロック14内で用いられている現在の多変数プロセスモデル14aよりも小さなモデル誤差を実現するものとしてその妥当性が検証された場合、モデル形式変更ブロック42は、適応された多変数プロセスモデルの形式または当該適応された多変数プロセスモデルを構成する適応されたＳＩＳＯモデルの形式を、プロセスモデルブロック14により用いられうる形式に変更する。したがって、たとえば、モデル形式変更ブロック42は、パラメトリックタイプのＳＩＳＯモデルをステップ応答モデルに形式変更し、ＭＰＣ制御システム12において用いるためこのステップ応答モデルをプロセスモデルブロック14に供しうる。
【００３４】
プロセス10の動作中、図１のスーパバイザ32は、ＣＶ信号、ＡＶ信号，ＤＶ信号，ＭＶ信号の各々の表示データを受信し、プロセス出力（ＣＶおよびＡＶ）の変化、プロセス入力（ＭＶ）の変化、および外乱入力（ＤＶ）の変化を検出すべく動作する。プロセス入力のうちのいずれか（たとえば、ＭＶまたはＤＶ）の変化が最低レベルを超えた場合、スーパバイザ32は、モデル評価プロシージャを開始し、多変数プロセスモデル14aに関連するＳＩＳＯモデルのうちのいずれを適応させる必要があるか否かを決定する。検出された変化は、瞬間的なもの、すなわち一回のスキャン期間にわたるものである場合や、または段階的なものであり、所定の期間（または、所定の回数のコントローラスキャン）にわたり検出される場合がある。事実、多変数プロセスでは、たいていの場合、変化は複数の入力と出力とに同時に発生する。このことは、本明細書記載の適応プロシージャの動作にとって理想的である。
【００３５】
プロセス入力において、（所望ならば、監視されている信号に基づいて変えてもよい）一または複数の所定量の変化を検出したあと、スーパバイザ32は、たとえば定常状態に至るまでの時間の如き特定の期間のプロセス入力データおよびプロセス出力データを収集・格納する。適応プロシージャを開始させると考えられる変化が、段階的なものである場合、すなわち複数スキャンにわたって発生する場合、スーパバイザ32は、適応プロシージャの開始に必要とされる所定量にその変化がいまだに達していない場合であったとしても、適応プロシージャ中に適応済多変数プロセスモデルの動作の分析に必要となるデータを入手するために、その変化の開始地点に関連する入力データおよび出力データを収集する必要がある。
【００３６】
いずれの場合であっても、プロセス入力データおよびプロセス出力データが収集され、メモリ（図１に明示的には図示せず）に格納されたあと、スーパバイザ32は、プロセスモデルブロック14において用いられている現在のプロセスモデル14aとプロセス10に対する測定入力または収集入力に基づくプロセス10からの実際の測定出力との間の誤差を評価し、必要があれば、多変数プロセスモデル14a内のＳＩＳＯモデルのうちのいずれを適応させる必要があるかを決定しうる。このように、すべてのＳＩＳＯモデルをいずれかの特定の適応プロセス中に適応させうるが、このようなことが、必ずしも必要とは限らず、実際、適応サイクルが長くなってしまう場合もある。これに代えて、スーパバイザ32は、以下で説明する適応プロシージャが適用されるＳＩＳＯモデルのサブセットをまず決定しうる。ここで、このサブセットは、その数が、プロセスモデル14a内のＳＩＳＯモデルの数よりも少ない。したがって、（プロセスモデルブロック14において用いられているＳＩＳＯモデルをすべて組み合わせたものである）多変数プロセスモデルを全体的にまたは部分的に適応しうる。ある種の出力に対するモデル誤差が非常に小さい場合、特定の出力に対するモデル誤差とそのプロセス入力のうちの一または複数との間に相関関係が存在しない場合、またはある種のプロセス入力の変化が非常に小さい場合に、部分的に適応させることが有利である。このような条件は、ＭＰＣ制御プロセスにおいて一般的なものであり、したがって部分的適応、すなわち上記ＳＩＳＯモデル全部の数よりも少ない数のＳＩＳＯモデルを特定の適応サイクルまたは適応プロシージャ中に適応させることが一般的でありうる。
【００３７】
ＳＩＳＯモデルの選択を実行するために、多変数適応プロシージャは、適応モデル構成、すなわち多変数モデルのうちの適応される部分（たとえば特定のＳＩＳＯモデル）をまず定義するモデル選択プロシージャを実施する。このモデル選択プロシージャには、入力データおよび出力データが収集される期間中においてスキャン毎に各出力ＯＵＴiに対するモデル出力誤差を計算することが含まれている。事実、この誤差は通常、モデル出力シフトとして通常のＭＰＣコントローラにおいて計算されており、プロセス測定値に依存するモデル出力の調整に用いられる。いずれの場合であっても、この誤差は以下の式を用いて計算されうる：
【００３８】
【数１】
0004722461.tif 000002
【００３９】
この式中、εj（ｔ）は時間ｔにおけるモデル出力誤差である。
【００４０】
ｙi（ｔ）は時間ｔにおけるプロセスの測定出力である。
【００４１】
【数２】
0004722461.tif 000003
【００４２】
は、多変数プロセスモデルにより生成される時間ｔにおける予測プロセス出力である。
【００４３】
次に、誤差εj（ｔ）とすべてのプロセス入力ｕj（ｔ）との間の相互相関ｒN（τ）をたとえば以下の式を用いて計算しうる：
【００４４】
【数３】
0004722461.tif 000004
【００４５】
この式中、Ｎは定常状態に至るまでの時間により定義される期間中に収集されるサンプルの数にほぼ等しい。
【００４６】
【数４】
0004722461.tif 000005
【００４７】
は、τがモデルデッドタイムとモデル時定数の半分とを足したものにほぼ等しいことを示している。この値は、実質的に、入力と出力との間の時間シフトにほぼ等しく、この時間シフトで、出力が入力の変化に対して最も敏感であり、したがって相関関数が最大値を示す。
【００４８】
ここで、各出力ＯＵＴiに関連するｊ個のＳＩＳＯモデルに対して、出力誤差と最高の相互相関を有する入力に関連する一または複数のモデルが適応候補として選択されうる。所望ならば、最高のまたは最大の相互相関を有する入力と関連するＳＩＳＯモデルが適応候補として選択されてもよく、所定の値を超える相互相関値を有するＳＩＳＯモデルが適応候補として選択されてもよく、または、他の相互相関基準に基づいてＳＩＳＯモデルが適応候補として選択されてもよい。
【００４９】
いかなる特定のＳＩＳＯモデルが適応候補として選択されてもよいが、選択されたＳＩＳＯモデルのすべてが実際の適応に最適であるというわけではない。具体的にいえば、特定のＳＩＳＯモデルに関連する入力に対して、適応期間中に、多量の変化または少なくとも必要最低量の変化が生じなかった場合、このモデルは、相互相関の値が高いとしても、よい適応候補とはいえない場合がある。というのは、実質的に、その入力の変化量が小さすぎるため、多量の変化が生じたものの相互相関の値が低い他の入力と関連しているＳＩＳＯモデルと比べて当該ＳＩＳＯモデルの方が出力誤差に対しより大きな影響を与えている指標として、上記の相互相関の計算結果に対し信頼を置くことができないからである。したがって、選択されたモデルは以下の条件を充足する必要がある：
【００５０】
【数５】
0004722461.tif 000006
【００５１】
通常動作データの変化が小さい場合、ＭＰＣコントローラ12は、選択されたプロセスの入力部にステップパルスを入力してもよいし、または選択されたプロセス出力に対して設定値の変更を行い、入力の変化が所定値を越えるようにし、それにより、適応が起こりうる状態を生成する。オプチマイザが動作状態である場合、このオプチマイザは、上記のステップパルスが入力部に与えられたあと、そのプロセス出力を初期値（最適値）に自動的に戻す。それに対し、ＭＰＣコントローラがオプチマイザを備えていない場合、当該ＭＰＣコントローラは、ある期間（パルス時間）ののちに、入力を最初の位置に自動的に設定しうる。入力および出力に対する最適パルス幅は、検討中の出力に対する主要な遅延時間とほぼ一致するものとしてヒューリスティックに特定される。
【００５２】
いずれの場合であっても、スーパバイザ32は、特定されたＳＩＳＯモデルの各々に関連する入力に、適応期間中（データ収集期間中）、所定量の変化または最低量の変化が生じたか否かを判定する。上記入力に所定量の変化が生じていなければ、そのＳＩＳＯモデルは実際に適応されるモデルとして選択されない。この場合、所望ならば、検討中の出力と二番目に高い相互相関を有する上記入力に関連するＳＩＳＯモデルを選択し、そのモデルの入力をチェックし、その入力に所定量の変化または最低量の変化が生じたか否かを判定してもよい。このプロセスは、特定の出力の各入力毎に繰り返され、その特定の出力に対して適応されるべき少なくとも一つのＳＩＳＯモデルを決定しうる。しかしながら、場合によっては、特定の出力に対して、その出力の出力誤差と入力の各々とが低い相互相関を有するため、適応期間（たとえば、データ収集期間）中に十分な変化を生じなかった入力と各ＳＩＳＯモデルが関連しているため、または計算により求められた出力誤差がそもそも小さいために、最終的に、どのＳＩＳＯモデルも、適応されるべきモデルとして決定されないことがある。
【００５３】
特定の出力ＯＵＴiに対して少なくとも一つのＳＩＳＯモデルが、適応される必要があるモデルとして決定される場合、そのＯＵＴiに対する適応に選択されなかった他のＳＩＳＯモデルは、モデルプロセスの出力に対してこれらの非選択モデルがどれだけ寄与しているかを算出すべく分析される。次いで、この寄与分は、プロセス出力の測定データから減算され、それにより、プロセス出力の測定データのうちの、適応された一または複数のＳＩＳＯモデルに実際に関連する（または、それらに起因する）部分が求められる。このプロシージャにより、適応システム30は、実際のプロセスの出力に対する特定のＳＩＳＯモデルの予測寄与分を分離することができ、それにより、予測寄与分と適応されたモデルにより生成される実際の寄与分（出力）との間の誤差を特定することができる。
【００５４】
いうまでもなく、ＳＩＳＯモデル適応・選択プロセスは、すべてのプロセス出力に対して繰り返してまたは別々に実行されるべくなしてあり、したがって、（必ずしもいつもとは限らないが）一般的にいえば、少なくとも一つのＳＩＳＯモデルが、プロセス10の各々のかつすべての出力ＯＵＴiに対して選択され、適応させられる。
【００５５】
適応させるべきモデルとして一または複数のＳＩＳＯモデルが決定されたあと、適応システム30は、いかなる公知のまたは所望の適応プロシージャを、これらのモデルに対して連続的にまたは平行して適用してもよい。一つの実施形態では、適応させるべく選択された各ＳＩＳＯモデルに対して、米国特許番号第6,577,908号に記載されているモデル適応プロシージャを適用してもよい。あるいは、本明細書において参照することによりその開示内容をここで明示的に援用する米国特許出願番号第10/419,582に記載されているモデル適応プロセスを用いてＳＩＳＯモデルを適応させてもよい。
【００５６】
ここで、これらのプロシージャを繰り返すことはしないが、一般的にいえば、これらの適応プロシージャにより、適応させるべきモデルに対するモデルセットが定義され、このモデルセットは図１のモデルセットブロック34内に格納される。このモデルがパラメトリックモデルである場合、モデルセット34は、個別のサブモデルを有しうるし、各サブモデルは、上記のモデルに関連するパラメータに対して異なる値を有しうる。したがって、モデルが三つのパラメータを有しそれぞれのパラメータが三つの異なる値をとりうる場合、モデルセットは、実質的に、２７個の異なるサブモデルを有する。これらのサブモデルの各々を上記収集されたデータを用いて分析することにより、プロセス10の実際の出力と比較して上記パラメータの各々に対するパラメータ値のうちいずれが最も小さな誤差を有するかが決定される。この誤差を決定するプロシージャが、モデルセット34内の特定のモデルの出力とプロセス10の実際の出力との間の差を求める加算器48により実行されるものとして、図１に示されている。図１に示されている加算器48への接続によりモデルセット34の出力とプロセス10の出力との間の比較がリアルタイムで実行されるように示されているが、これらの比較は、実際のところ（そして通常そうであるように）、スーパバイザ32により前もって収集・格納されたプロセス入力データおよびプロセス出力データに基づいて、適応プロシージャの開始時点において実行されてもよい。
【００５７】
米国特許番号第6,577,908号により明らかなように、適応されるＳＩＳＯモデルを構成する一組のモデルは、任意の所望のかつ事前定義されたスイッチングルールにより自動的に選択されうる複数のモデルサブセットを有しているものとして構築される。個々のモデルは、複数のパラメータを備えうるし、各パラメータは、そのパラメータ対応する一組の所定の初期値から選択される値をとりうる。モデルセット34内の個々のモデルの評価には、モデル二乗誤差、すなわちノルムの計算が含まれる。このノルムは、評価されるモデルにおいて表現されるすべてのパラメータに対して割り当てられる。モデルが繰り返し評価されていくとともに、モデル評価の過程においてパラメータに割り当てられたノルムのすべての合計である累積ノルムが各パラメータに対して算出されていく。そのあと、各パラメータに割り当てられた初期値の重み付け平均である適応パラメータ値が各パラメータに対して算出されうる。
【００５８】
米国特許出願番号第10/419,582号に記載されているように、モデルセットコンポーネントはまた、複数のプロセス領域を定義する状態変数と、該複数のプロセス領域にグループ分けされる複数のプロセスモデルとを有しうる。このケースの場合、プロセスモデルの各々はまた、複数のパラメータを有しており、各パラメータは、そのパラメータに割り当てられた一組の所定の初期値から選択される値を有している。上記領域の各々は、その領域に対して定義された（異なる）一組の標準パラメータ値を有しうる。この場合、モデル評価ブロック36は、モデルセット34内のプロセスモデルの出力とプロセス10の出力との間の差を表すモデル誤差信号を解析し、プロセスモデルに対応するモデル二乗誤差を算出し、このモデル二乗誤差をプロセスモデル内に表されるパラメータ値に帰属させうる。パラメータ／属性補間器38は、モデル評価ブロック36に通信可能に接続され、プロセスモデルにおいて表されているパラメータ（または、属性）に対して適応プロセスパラメータ値（または、モデル属性値）を算出する。
【００５９】
いずれの場合であっても、いうまでもなく、モデル評価ブロック36は、誤差を最も小さくさせるように、モデルセット内のモデルの各々のパラメータ（または、他の属性）に対する重みを決定する。米国特許番号第6,577,908号に記載されているように、これらの重みは、異なるパラメータ値の各々により生成される誤差に対する寄与に基づいて算出されてもよい。さらに、米国特許番号第6,577,908号に記載されているように、算出されたパラメータは、その設計において想定される、中央パラメータ値ｐk（ａ）（ｋ＝１、…ｍ）とパラメータ変化の範囲とを含んだ新規のモデルセットを規定しうる。変更範囲は、±Δ％で規定されると都合がよいばあいがある。その範囲内において、少なくとも二つのパラメータを定義する必要がある。実際的には、さらに二つのパラメータ、ｐk（ａ）＋Δ％ｐk（ａ）／１００、ｐk（ａ）−Δ％ｐk（ａ）／１００が、ｐk（ａ）の前後に定義されうる。すべてのパラメータは、適応下限と適応上限とを定義しており、パラメータｐk（ａ）がその限界値を超えるような場合には、その限界値に戻される。もちろん、各モデルが最小限の数のパラメータにより表されることが好ましい。この技法を用いて数多くのプロセスをモデル化するのに、一次遅れ＋むだ時間プロセスモデルが適していると見なされている。
【００６０】
パラメータに対する重み付けがモデル評価ブロック36により決定されたあと、パラメータ／属性補間ブロック38は、このモデル評価ブロック36により決定されたパラメータ重み付け値を用いて生成されるパラメータ値を有する新規のまたは適応されたプロセスモデルを定義すべく動作する。
【００６１】
特定のプロセス出力に対し選択されたすべてのＳＩＳＯモデルが適応されたあとは、次のプロセス出力が選択され、上記のプロシージャが、適応される必要があるとして選択されたＳＩＳＯモデルのすべてに対しパラメータ補間が実行されるまで、繰り返される。そのあと、適応されたＳＩＳＯモデルおよび適応されていないＳＩＳＯモデルからなる多変数モデルは、すべての出力に対して別々に妥当性検査が行われうる。モデル妥当性検査ブロック40は、上記の妥当性検査プロシージャを実行し、その適応結果を受け入れるかまたは拒絶する。所望ならば、適応された多変数プロセスモデルが受理されるべきか否かを決定するために、二段階の判断基準を設定してもよい。この場合、適応されたモデルの出力誤差は、事前に定義された値よりも小さいことが必要であり、適応されたモデルの出力誤差は、現在のモデルの出力誤差よりも小さいことが必要である。いうまでもなく、現在のモデルと適応されたモデルとの分析には同一のデータセットを用いることが必要である。
【００６２】
多変数モデルが妥当性検査をされ受理されたあと、必要ならば、モデル形式変更ブロック42は、適応されたモデルを、ＭＰＣコントローラにより実際に用いられている形式に変更してもよい。したがって、パラメトリックモデルは、（ＭＰＣコントローラにおいて通常そうであるように）ＭＰＣコントローラがステップ応答モデルを用いている場合には、ステップ応答モデルに形式変更されることが必要である。これに代えて、ステップ応答モデルは、ＭＰＣコントローラがパラメトリックモデルを用いている場合には、パラメトリックモデルに形式変更されることが必要である。もちろん、形式変更されたモデルは、そのあと、プロセスモデルブロック14で使用されている現在の多変数プロセスモデル14aと交換すべく用いられる。この場合、多変数プロセスモデル14a全体を交換してもよいし、または適応プロシージャ中に実際に適応されたＳＩＳＯモデルのみを交換してもよい。
【００６３】
以上では、多変数プロセスモデル適応プロセスが記載されており、またそこで用いられるとともに米国特許番号第6,577,908号に記載されているＳＩＳＯモデル適応プロセスが、パラメータ補間が実行されるパラメトリックモデルを用いるものとして記載されているが、本明細書に記載されている技法を用いると、ステップ応答モデルの如きノンパラメトリックモデルの形式のＳＩＳＯモデルを使用して適応することも可能である。図４は、ＳＩＳＯモデル適応がノンパラメトリックモデルに対して実行されうる様式の一例となるノンパラメトリックステップ応答モデル50を示している。具体的にいえば、一または複数のモデルパラメータの値を所定量だけ変更することによりパラメトリックモデルを分析することに代えて、ノンパラメトリックモデルに関連するモデル属性は、適応プロセスにおいて検討されるさまざまなノンパラメトリックモデルを定義すべく、所定の方法で変更されうる。
【００６４】
したがって、図４のステップ応答モデル50に関して、（図１の）モデルセット34に検討されるべきさまざまなモデルは、図４のステップ応答曲線の一部または全部にわたるその特定の性状または属性を変更することにより定義されうる。一例としては、図４の線52により示されているように、すべてのステップ応答係数において、同一の割合だけステップ係数値を増加または減少させることにより、さらに二つのステップ応答を定義しうる。この変更は、ステップ応答曲線のスケーリングの変更であるが、一次遅れ＋むだ時間パラメトリックモデルのゲインを変更することと同等である。同様に、図４の線54により示されているように、図４のステップ応答内で出力応答を最初に見られる時間を時間的に左側または右側に移動することができる。このモデル属性の変更は、一次遅れ＋むだ時間パラメトリックモデルのデッドタイムを変更（増加または減少）させることと同じである。さらに、ステップ応答の起きるポイントにおいて揺れ動く線により定義される値の右側または左側にステップ応答値を移動させることにより、見かけステップ遅れを変更することができる。すなわち、図４の線56により示されているように、応答曲線の傾斜を変更（増加または減少）することにより、さまざまな遅れ時間または応答時間を有するステップ応答を定義することができる。この変更は、一次遅れ＋むだ時間パラメトリックモデルの時定数を増加または減少させることと同じである。
【００６５】
このアプローチの利点は、当該アプローチにより、パラメータによって表すことができないような場合がある非常に複雑なＳＩＳＯモデルの変更や適応が可能となり、それにより、ノンパラメトリックモデルまたは他の複雑なプロセスモデルを用いた多変数プロセスモデルの適応が可能となることである。もちろん、本明細書においてステップ応答曲線に対する三つの特定の変更が示されているが、本明細書に記載の適応プロセスに検討するさまざまなモデルを定義するために、ステップ応答曲線の他の属性に対する他の変更および、他のノンパラメトリックモデルの属性に対する同一の変更または異なる変更が適用されてもよい。さらに、いうまでもなく、ノンパラメトリックＳＩＳＯモデル適応が本明細書において説明されたような属性変更を用いて実行される場合、図１の適応システム30は、モデル適応の実行のために、パラメータ設定値の代わりに属性設定値を内挿する。その結果、モデルセット34は、パラメトリックＳＩＳＯモデルに対して前もって設定されたパラメータ値を格納・切り換えするのと同一の方法で、ノンパラメトリックモデルに対して適応を実行するときに二つ以上の異なる属性値を格納・切り換えうる。
【００６６】
本発明は、当該発明の教示および例示を意図する実施形態の一例を参照して記載されているが、開示されている適応多変数コントローラはこれらの示威し形態に限定されるものではない。当業者により、本発明の範疇から逸脱することなく、さまざまな変更、改良、および追加が加えられうる。
【００６７】
たとえば、上述のように、コントローラ適応は、制御されるプロセスの数学的モデルの構築に用いられたパラメータまたは属性の統計的内挿を基づいている。このプロセスは三つのパラメータまたは属性により特徴付けされるし、これらのパラメータの各々は三つの値が割り当てられうるが、明らかに、開示されている適応多変数コントローラは、他のパラメータおよび／または異なる数のパラメータ（または、属性）を有すべく拡張されるように構成されており、各々のパラメータは異なる数の割り当て値を潜在的に有しうる。さらに、所望ならば、パラメータ／属性補間38は、モデル評価ブロック36により定義されるパラメータ重みまたは属性重みに基づいてパラメータまたは属性すべての間において内挿することに代えて、最も小さな関連付け誤差を有する二つのパラメータ間の如きパラメータ（または、属性）からなるサブセット間において内挿を行ってもよい。このようにすることにより、各ＳＩＳＯモデル適応に関連する内挿時間および演算回数が減少されうる。
【００６８】
同様に、適応システム30は、パラメトリックモデルの全部のもしくは一部のパラメータに基づいて、またはノンパラメトリックモデルの全部のまたは一部の属性に基づいて適応させるように動作しうる。具体的にいえば、適応システム30は、適応プロシージャ中に、（モデルセット34内の）任意の特定のモデルに存在する一または複数の「重要な」パラメータまたは属性に対して焦点を当て、一または複数のその他のパラメータまたは属性には適応または変更を加えないことが可能である。このようなパラメータまたは属性を選択する技術を用いることにより、適応システム30が特定のＳＩＳＯモデルを適応させるのに必要な時間を著しく削減することが可能となる。このことは、適応させる必要があるＳＩＳＯモデルの数が大きい多変数プロセスにおいては重要なことである。
【００６９】
さらに、モデル評価およびパラメータ補間が、モデルセット24、モデル評価ブロック36、スーパバイザ32、パラメータ／属性補間器38、モデル妥当性検査ブロック40、およびモデル形式変更ブロック42として特定される個別のコンポーネントとして示されている。当業者にとって明らかなように、個別のコンポーネントへの分割はコントローラの実装および運転に関する責務を有する人により自由裁量されるものであり、また、これらの機能のすべてはいかなる所望の方法で実現されてもよい。さらに、本明細書に記載される適応多変数プロセス制御システムはソフトウェアにより実現されるのが好ましいが、その全部または一部が、ハードウェア、ファームウェアなどにより実現されてもよく、プロセス制御システムに関連するその他のプロセッサにより実行されてもよい。したがって、本明細書記載のエレメントは、標準型多目的ＣＰＵにより、またはたとえば特定用途向け集積回路（ＡＳＩＣ）もしくは所望な他のハードワイヤードデバイスの如き特別に設計されたハードウェアもしくはファームウェアにより実現されてもよい。ソフトウェアで実現される場合、該ソフトウェアルーチンは、磁気ディスク、レーザーディスク（たとえば、ＣＤ、ＤＶＤなど）、または他の記憶媒体などのコンピュータ読取り可能メモリに、コンピュータあるいはプロセッサのＲＡＭまたはＲＯＭに、任意のデータベースなどに格納されてもよい。同様に、このソフトウェアは、たとえばコンピュータ読取り可能ディスク、スマートカードメモリ、または他の伝達可能なコンピュータ記憶メカニズム、または、電話回線、インターネットなどの通信チャネルを用いることを含む既知のまたは所望の伝達方法を介して、（上記の電話回線、インターネットなどの通信チャネルを用いて伝達することは、移動可能記憶媒体を介してかかるソフトウェアを伝達することと同様または相互交換可能であるとみなされる）ユーザまたはプロセスプラントへ伝達されてもよい。
【００７０】
上記の特定のアプローチは本発明の上記実施例と実質的に変わらないと認められる。したがって、特許請求の範囲は、本発明の真の精神および範疇に該当するすべての変形、変更、および改良ならびにそれらと略均等なものを内包すると解釈されるのが正しい。したがって、本明細書において詳細に記載されていないが、本発明の他の実施例であっても、特許請求の範囲に定義されている本発明の範疇内であると考えられる。
【図面の簡単な説明】
【００７１】
【図１】プロセスモデルのパラメータまたは属性の補間に基づいて動作する適応多変数制御システムを示す機能ブロック線図である。
【図２】多入力多出力プロセスコントローラを示す概念図である。
【図３】一組の単一入力単一出力モデルが多入力多出力プロセスモデルの出力を生成すべく用いられうる一つの様態を示す概念図である。
【図４】ステップ応答モデルの如きノンパラメトリックモデルがモデル属性補間技法を用いて適応されうる様態を示すグラフである。
【符号の説明】
【００７２】
8  制御システム
10  プロセス
12  ＭＰＣ制御システム
14  プロセスモデルブロック
14a  多変数プロセスモデル
16  ＭＰＣコントローラ
18  ベクトル加算器
20  設定値予測ブロック、設定値予測ユニット
24  モデルセット
30  モデル適応ネットワーク、適応システム
32  スーパバイザ
34  モデルセットブロック
36  モデル評価ブロック
38  パラメータ／属性補間、パラメータ／属性補間器
40  モデル妥当性検査ブロック
42  モデル形式変更ブロック
48  加算器
50  ノンパラメトリックステップ応答モデル
52、54、56  線
【図１】
0004722461.tif 000007
【図２】
0004722461.tif 000008
【図３】
0004722461.tif 000009
【図４】
0004722461.tif 000010