相対ゲイン行列(RGA)シミュレーター 戻る
制御工学

相対ゲイン行列(RGA)シミュレーター

複数の入力と出力を持つ多変数プロセスを、分散型の単一ループ制御で扱うとき、どの入力をどの出力に組み合わせればよいかを決めるツールです。2×2のゲイン行列を入力すると、ブリストルの相対ゲイン行列、推奨ペアリング、ループ間の相互干渉の強さがリアルタイムで分かります。

パラメータ設定
2×2 定常プロセスゲイン行列 K の各要素。入力 u_j が出力 y_i に与えるゲインが g_ij です。
ゲイン g11(u1 → y1)
ゲイン g12(u2 → y1)
入力2が出力1に与える相互干渉ゲイン
ゲイン g21(u1 → y2)
入力1が出力2に与える相互干渉ゲイン
ゲイン g22(u2 → y2)
計算結果
RGA要素 λ11
RGA要素 λ12 (=1−λ11)
推奨ペアリング
相互作用の強さ
ニーダリンスキ指数 NI
制御性の判定
プロセス構造図とRGAマップ

左の2入力(u1, u2)と右の2出力(y1, y2)を結ぶ4本のゲイン経路。線の太さは |g_ij| に比例します。右側はRGA行列のマップ(緑=1付近の良いペア、赤=負、黄=0.5付近の強干渉)。推奨ペアリングの経路をパルスで強調します。

λ11 と相互干渉ゲイン g12 の関係
RGA要素の比較(λ11・λ12・λ21・λ22)
理論・主要公式

$$\lambda_{11}=\frac{g_{11}g_{22}}{g_{11}g_{22}-g_{12}g_{21}},\qquad \text{RGA}=\begin{bmatrix}\lambda_{11}&1-\lambda_{11}\\1-\lambda_{11}&\lambda_{11}\end{bmatrix}$$

2×2プロセスの相対ゲイン行列。g_ij はゲイン行列 K の要素。RGAの各行・各列の和は常に1になり、負の要素を持つペアリングは避ける。

$$\det K = g_{11}g_{22}-g_{12}g_{21},\qquad \text{NI}=\frac{\det K}{g_{11}g_{22}}$$

ゲイン行列の行列式 detK と、対角ペアリングのニーダリンスキ指数 NI。|detK| が0に近いと行列は特異(制御不能)。NI が負なら対角ペアリングは構造的に不安定。

相対ゲイン行列(RGA)とは

🙋
「相対ゲイン行列」って、なんだか難しそうな名前ですけど、これは何を計算する道具なんですか?
🎓
ざっくり言うと「どのつまみが、どのメーターを担当すべきか」を教えてくれる行列だよ。たとえば蒸留塔みたいに、バルブが2つ(u1, u2)あって、測りたい温度や濃度も2つ(y1, y2)あるとする。普通はu1でy1を、u2でy2を、というふうに別々のPID制御ループで回したい。でも困ったことに、u1を動かすとy1だけじゃなくてy2まで動いてしまう。この「ループ同士の干渉」を数値で見える化するのが、エドガー・ブリストルが考えた相対ゲイン行列、略してRGAなんだ。
🙋
なるほど。でも、ゲイン行列 K を見れば干渉してるかどうか分かりそうですけど、わざわざRGAにするのはなぜですか?
🎓
いい疑問だね。生のゲイン行列 K には落とし穴があって、入力や出力の「単位」を変えるだけで数字がガラッと変わってしまうんだ。温度を℃で測るかKで測るか、流量をL/minにするかm³/hにするかで K の見た目が変わる。これだと「干渉が強い」と言っても客観性がない。RGAの素晴らしいところは、各要素を自分自身のゲインと、他のループを閉じたときのゲインの「比」で定義しているから、単位やスケーリングをどう取っても値が変わらないんだ。だから設計の早い段階で、機器の選定すら決まる前でも使える。
🙋
RGAの数字を見たら、具体的にどう読めばいいんですか?左で g12 を動かすと λ11 が変わっていきますね。
🎓
読み方はシンプルだよ。RGA要素が1に近ければ「そのペアはほとんど干渉しない、きれいなペア」。0.5付近なら「2つのループが互いに半分ずつ引っ張り合う、激しい干渉」。そして要素が負になったら危険信号だ。1を超えて大きい場合も「悪条件」で、わずかなモデル誤差に弱くなる。原則は「正で、いちばん1に近い要素のところで入力と出力を組む」。デフォルトの行列だと λ11≈1.09 だから、対角ペアリング(y1-u1, y2-u2)が推奨される、というわけだね。
🙋
RGA要素が「負」だと、なぜそんなにダメなんですか?
🎓
これがRGAでいちばん大事な警告なんだ。RGA要素が負ということは、「他のループを開いているとき」と「閉じているとき」で、そのペアの利得の符号が逆になることを意味する。たとえば制御器を設計するとき、他のループが手動(開)の状態を見て「u1を増やせばy1が上がる」と思って正のフィードバックを組んだのに、実際に他のループも自動(閉)にした途端、u1を増やすとy1が下がるようになってしまう。これでは制御系が暴走する。だから負のRGA要素のペアは、どんなに調整を頑張っても絶対に避けるんだ。
🙋
「ニーダリンスキ指数」というのも出てきますけど、これはRGAとは別物ですか?
🎓
補い合う関係だね。RGAは「どのペアを選ぶか」を教えてくれるけど、対角ペアリングを選んだとして、それを全ループ閉じたときに本当に安定できるかは別途チェックがいる。そこで使うのがニーダリンスキ指数 NI = detK/(g11·g22) だ。NIが負だと、積分動作を含むPID制御ではどう調整しても不安定になる、という構造的なNG判定が出る。実務では「RGAで正の要素のペアを選ぶ」「NIが正であることを確認する」、この2つをセットで通して、初めてそのループ構成で先に進む、という流れになる。

よくある質問

相対ゲイン行列(RGA)は、エドガー・ブリストルが提案した、多入力多変数プロセスのループ間の相互干渉を評価する行列です。2×2系では λ11 = (g11·g22)/(g11·g22−g12·g21) で計算し、RGA全体は [[λ11, 1−λ11],[1−λ11, λ11]] となります。RGAの大きな特徴は、入力や出力の単位・スケーリングに依存しないことです。本ツールはゲイン行列からRGAを計算し、どの入力をどの出力に組み合わせるべきかを判定します。
原則は「RGA要素が正で、かつ1に最も近い対」を組み合わせることです。2×2系では λ11 ≥ 0.5 なら対角ペアリング(y1-u1, y2-u2)、λ11 < 0.5 なら非対角ペアリング(y1-u2, y2-u1)を選びます。RGA要素が負のペアは絶対に避けます。負の要素のループを閉じると、他のループを開いた状態と閉じた状態で利得の符号が反転し、不安定化を招くためです。
λ11 が0.5付近の値は、2つのループが互いに強く干渉していることを意味します。ある入力を動かすと両方の出力が同程度に変化するため、片方のループを単独で調整しても、もう片方が乱れて元に戻ってしまいます。この状態では分散型の単一ループ制御は非常に難しく、整定に時間がかかったり振動が止まらなかったりします。対策としては、デカップリング(非干渉化補償器)の導入や、多変数制御(モデル予測制御など)への切り替えを検討します。
ニーダリンスキ指数 NI = detK/(g11·g22) は、対角ペアリングで全ループを閉じたときの構造的な安定性を判定する指標です。NI が負であれば、その対角ペアリングは積分動作を含む限りどう調整しても不安定になります(ニーダリンスキの不安定条件)。NI が正であることは安定の必要条件であり、十分条件ではありません。RGAで負の要素がないことと、NIが正であることの両方を満たして初めて、そのペアリングが構造的に成立します。

実世界での応用

蒸留塔の制御:RGAが最も古典的に使われてきたのが蒸留塔です。還流量と再沸器熱量という2つの操作変数で、塔頂と塔底の組成という2つの制御変数を保ちます。両者は本質的に強く干渉し、典型的なLV構成ではRGA要素が大きな正の値や負の値になることが知られています。RGAの符号と大きさから、組成制御をどのバルブに割り当てるか、デカップリングが必要かを設計初期に判断します。

空調・HVACシステム:1台の空調機で温度と湿度を同時に制御する場合、冷却コイルの能力と再熱・加湿の操作が互いに干渉します。冷却は温度を下げると同時に除湿も進めてしまうため、温度ループと湿度ループのペアリングを誤ると、両者が永遠に追いかけ合います。RGAでどちらの操作をどちらの測定に割り当てるべきかを定量的に決められます。

化学反応器・混合プロセス:反応器の温度と圧力、あるいはブレンドタンクの流量比と総流量など、操作変数同士が幾何学的・熱力学的に結合したプロセスは多くあります。RGAは反応器のジャケット温度・原料流量の割り当てや、pH中和プロセスの酸・アルカリ投入ループの選定に使われ、相互干渉が制御性を損なう前に構造を見直す判断材料になります。

多変数制御への切り替え判断:RGAは単一ループ制御を続けるか、モデル予測制御(MPC)のような多変数制御に切り替えるかの判断にも使います。RGA要素が1付近に収まっていれば分散型PIDで十分ですが、0.5付近や負の値、あるいは極端に大きい値が並ぶ「悪条件」のプロセスでは、デカップラやMPCへの投資が正当化されます。プラント全体の制御アーキテクチャを決める最初のスクリーニングツールなのです。

よくある誤解と注意点

まず多いのが、「RGAだけ見ればペアリングが完全に決まる」という思い込みです。RGAは定常状態(周波数ゼロ)のゲインだけを見た指標で、プロセスの動特性、つまり時定数やむだ時間、伝達関数の位相は一切反映していません。定常RGAでは良いペアリングに見えても、片方のループが極端に速く、もう片方が遅い場合には、動的には強く干渉することがあります。本格的な設計では、対象周波数(クロスオーバー付近)での周波数依存RGAを併用するのが定石です。本ツールの定常RGAは、あくまで構造選定の第一スクリーニングと位置づけてください。

次に、「ニーダリンスキ指数が正なら安定」という誤解です。NIが正であることは安定の必要条件にすぎず、十分条件ではありません。NIが負なら「その対角ペアリングは積分制御を含む限り必ず不安定」と断言できますが、NIが正でも実際には不安定になり得ます。NIはあくまで「構造的にあり得ないペアリングを除外する」ふるい分けの道具です。最終的な安定性は、制御器を実際に設計したうえで、ナイキスト判定やシミュレーションで確認する必要があります。

最後に、「RGA要素が大きいほど干渉が強い」と単純に考えてしまうことです。干渉が最も激しいのは、実は λ11 が0.5付近のときです。λ11 が1付近、あるいは0付近なら干渉は弱く、きれいに分離できます。一方で λ11 が1を大きく超える、あるいは負になる場合は「強い・悪条件」の干渉で、これはモデル誤差に対する感度が極端に高い状態を意味します。つまりRGA要素は「1からどれだけ離れているか」と「符号」の両方で読む必要があり、単に絶対値の大小で判断するのは誤りです。本ツールは |λ11−0.5| と符号の両方から相互作用の強さを分類しています。

使い方ガイド

  1. プロセス伝達関数の4要素(g11、g12、g21、g22)を複素数または周波数応答データで入力します。化学反応器の温度・圧力制御では、例えば温度から圧力への相互干渉ゲイン g12=0.85 のように設定します。
  2. 各ゲイン要素の周波数範囲を指定して、RGA(相対ゲイン行列)のブリストル行列を計算します。制御周波数帯0.1~10rad/sでのゲイン変動を評価することで、ループペアリングの周波数依存性を確認できます。
  3. シミュレーター出力のRGA要素λ11値と推奨ペアリングに基づき、制御ペアリング設計(u1→y1とu2→y2、または対角ペアリング)を選定します。λ11が0.5に近い場合は相互干渉が強く、デカップリング制御の導入を検討します。

具体的な計算例

蒸留塔の液面・組成同時制御の場合:g11=2.5(リボイラー蒸気流量→塔底組成)、g12=0.32、g21=0.18、g22=1.8とします。計算結果λ11=0.68、λ12=0.32となり、推奨ペアリングは対角型(リボイラー→組成、オーバーヘッド→液面)です。ニーダリンスキ指数NI=0.42で相互作用は中程度、制御ループ設計時に軽微なデカップリング補償器の追加で十分安定化できます。

実務での注意点

  1. RGA要素の周波数依存性:低周波(0.01rad/s)ではλ11=0.75でも高周波(50rad/s)で0.45に低下する場合があります。ボイラー給水制御では、低周波での遅延支配と高周波での直結特性を区別して制御則を設計してください。
  2. 非最小位相ゼロの検出:推奨ペアリング表示が不安定化警告を示す場合、逆応答ゲインの存在を示唆しています。熱交換器入口温度制御では、バイパス流量と本流の逆動作特性を確認し、積分補償器の遅延時間を調整します。
  3. ニーダリンスキ指数NI>0.5の強相互干渉系では、単一ループPID制御の応答余裕が不足します。マルチループロバスト制御(H∞またはμ合成)の導入判定に本指数を活用してください。