制御工学・システム制御の分野を学ぶ際に必ず登場する概念の一つが「利得行列」です。
フィードバック制御・状態フィードバック・最適制御など、現代制御理論の核心部分において利得行列は中心的な役割を果たしています。
ロボットの動作制御・航空機の飛行制御・自動車の安定制御・産業用機械の精密制御など、私たちの身の回りの高度な制御システムの多くが利得行列を活用した制御則に基づいて動いています。
この記事では、利得行列とは何かという基本的な定義から、その計算方法・状態フィードバックとの関係・LQR制御など最適制御への応用・具体的な計算例まで、体系的かつ丁寧に解説していきます。
制御工学を学ぶ学生の方から、制御システム設計の実務に携わるエンジニアの方まで、幅広くお役に立てる内容となっているでしょう。
目次
利得行列とは状態フィードバック制御における制御入力を決める行列のこと
それではまず、利得行列の基本的な定義と概念について解説していきます。
利得行列（Gain Matrix）とは、制御理論においてシステムの状態変数から制御入力を計算するための行列です。
状態フィードバック制御では、制御入力u（t）をシステムの状態ベクトルx（t）の線形結合として設計します。
状態フィードバック制御の基本式
u（t） ＝ −K × x（t）
u：制御入力ベクトル（m次元）
K：利得行列（m×n行列、mは入力の数、nは状態変数の数）
x：状態変数ベクトル（n次元）
この式において、Kが利得行列であり、状態変数をどの程度制御入力に反映させるかを決定します。
利得行列Kの各要素（ゲイン値）をどのように設定するかが制御設計の核心であり、閉ループ系の安定性・応答速度・定常偏差・外乱抑制性能などに直結します。
適切な利得行列を設計することで、不安定なシステムを安定化したり、所望の応答特性を実現したりすることができます。
状態空間表現と利得行列の関係
利得行列を理解するためには、状態空間表現（State-Space Representation）の基礎を押さえておく必要があります。
線形時不変システム（LTIシステム）の状態空間表現は以下の連立微分方程式で表されます。
状態方程式：ẋ（t） ＝ A × x（t） ＋ B × u（t）
出力方程式：y（t） ＝ C × x（t） ＋ D × u（t）
A：システム行列（n×n）
B：入力行列（n×m）
C：出力行列（p×n）
D：直達行列（p×m）
x：状態ベクトル（n次元）、u：入力ベクトル（m次元）、y：出力ベクトル（p次元）
状態フィードバック制御 u ＝ −Kx を状態方程式に代入すると、閉ループ系の状態方程式が得られます。
ẋ ＝ Ax ＋ B(−Kx) ＝ (A − BK)x
閉ループ系のシステム行列：Acl ＝ A − BK
閉ループ系の安定性は Acl の固有値（極）によって決まります。
すべての固有値の実部が負であれば、閉ループ系は漸近安定です。
つまり、利得行列Kを適切に設計することで閉ループ極を所望の位置に配置し、所望の安定性と応答特性を実現することが状態フィードバック設計の目標です。
利得行列の次元と各要素の意味
利得行列Kはm×n行列（mは制御入力の数、nは状態変数の数）であり、その各要素が重要な意味を持ちます。
たとえば1入力・2状態変数のシステムでは、Kは1×2行列（行ベクトル）K ＝ [k₁ k₂] となります。
この場合、制御入力は u ＝ −k₁x₁ − k₂x₂ と表され、k₁は状態変数x₁に対するゲイン、k₂は状態変数x₂に対するゲインを意味します。
ゲイン値が大きいほどその状態変数が制御入力に強く反映され、応答が速くなる一方でオーバーシュートや振動が生じやすくなります。
ゲイン値が小さすぎると応答が遅くなり、外乱に弱くなります。
適切なゲインのバランスを見つけることが制御設計の技術的な核心です。
利得行列の計算方法
続いては、利得行列を実際に計算・設計する方法について確認していきます。
利得行列の設計には複数のアプローチがあり、システムの要求仕様に応じた適切な方法を選択することが重要です。
極配置法による利得行列の設計
極配置法（Pole Placement）は、閉ループ系の極（固有値）を所望の位置に配置するように利得行列を設計する方法です。
極配置法の手順：
手順1：制御仕様から所望の閉ループ極 {p₁, p₂, …, pn} を決定する
（例：2次系で減衰係数ζ＝0.7、固有角周波数ωn＝10 rad/sを目標とする場合）
（所望極：p₁,₂ ＝ −ζωn ± jωn√(1−ζ²) ＝ −7 ± 7.14j）
手順2：所望の特性多項式を計算する
（A−BK）の特性多項式 ＝ (s−p₁)(s−p₂)…(s−pn)
手順3：特性多項式の係数を比較してK（利得行列）を求める
Ackermannの公式などを使って系統的に計算できます。
極配置法は直感的でわかりやすい設計手法ですが、多入力システムでは利得行列の解が一意に定まらないという問題があります。
また、所望極を複素平面の左半面に正しく配置することが安定化の必要条件であり、どの程度左に配置するかが応答速度・ゲインの大きさ・制御入力の飽和可能性に影響します。
LQR（線形二次レギュレータ）による最適利得行列の設計
現代制御理論で最も広く使われる利得行列設計手法の一つがLQR（Linear Quadratic Regulator：線形二次レギュレータ）です。
LQRは、以下の評価関数（コスト関数）を最小化する最適な利得行列Kを求める手法です。
LQRの評価関数：
J ＝ ∫₀^∞ [xᵀQx ＋ uᵀRu] dt
Q：状態変数に対する重み行列（n×n、半正定値対称行列）
R：制御入力に対する重み行列（m×m、正定値対称行列）
xᵀQx：状態変数の偏差を小さくする項（制御精度に関わる）
uᵀRu：制御入力を小さくする項（省エネ・アクチュエーター保護に関わる）
LQRの最適利得行列Kは、リカッチ方程式（Riccati Equation）を解くことで求められます。
代数リカッチ方程式（ARE）：
AᵀP ＋ PA − PBR⁻¹BᵀP ＋ Q ＝ 0
Pを解いた後、最適利得行列は：
K ＝ R⁻¹BᵀP
Qの対角要素を大きくすると対応する状態変数の偏差を小さくしようとする制御になり応答が速くなりますが制御入力が大きくなります。
Rを大きくすると制御入力を小さく抑えた省エネルギー的な制御になりますが、応答が遅くなります。
このQとRの重みのバランスを調整することで、応答性と省エネのトレードオフを設計者が直感的に調整できる点がLQRの大きな特徴です。
具体的な利得行列の計算例
倒立振子の制御を例に、利得行列の計算を具体的に示します。
簡略化した倒立振子の状態空間モデル（2状態変数の場合）：
A ＝ [[0, 1], [g/l, 0]]（g：重力加速度≒9.8 m/s²、l：振子長さ＝1mとする）
B ＝ [[0], 
] ＝ [[0], ]
A ＝ [[0, 1], [9.8, 0]]
所望閉ループ極：p₁,₂ ＝ −3 ± 3j（減衰係数0.707、固有角周波数約4.24 rad/s）
所望特性多項式：(s＋3＋3j)(s＋3−3j) ＝ s² ＋ 6s ＋ 18
閉ループ行列 A−BK ＝ [[0,1],[9.8−k₁,−k₂]] の特性多項式：
s² ＋ k₂s ＋ (k₁−9.8) ＝ s² ＋ 6s ＋ 18
係数比較より：k₂ ＝ 6、k₁−9.8 ＝ 18 → k₁ ＝ 27.8
利得行列 K ＝ [27.8, 6]
この例からわかるように、极配置法では所望の閉ループ極から逆算して利得行列の各要素を求めることができます。
利得行列の応用分野と実用例
続いては、利得行列が実際の制御システムにどのように応用されているかを確認していきます。
利得行列を活用した制御系は、現代のさまざまな高度な自動化・制御システムの中核を担っています。
ロボット制御における利得行列の活用
産業用ロボットや自律型ロボットの関節制御において、利得行列は非常に重要な役割を果たします。
多関節ロボットアームの制御では、各関節の角度・角速度・トルクなどを状態変数として状態空間モデルを構築し、LQRなどの手法で最適な利得行列を設計します。
関節ごとの利得値（Kp：比例ゲイン、Kd：微分ゲイン、Ki：積分ゲイン）を適切に設定することで、高精度・高応答性・外乱抑制性能を兼ね備えたロボット制御が実現されます。
近年の協働ロボット（人間と共に作業するロボット）では、接触力のフィードバックを含む利得行列設計が安全性の確保にも重要な役割を果たしています。
航空機・宇宙機の飛行制御への応用
航空機の自動操縦装置（オートパイロット）や宇宙機の姿勢制御システムにも、利得行列を用いた状態フィードバック制御が広く使われています。
航空機の縦方向の運動（ピッチ制御）を例にとると、機首上下の角度（ピッチ角）・角速度・高度・速度などを状態変数として含む高次元の状態空間モデルが使われます。
LQRやH∞制御などの最適制御手法によって利得行列を設計することで、乱気流などの外乱に対しても安定した飛行を維持できる制御系が実現されます。
ロケットの誘導制御・人工衛星の姿勢制御・ドローンの自律飛行なども、利得行列を活用した現代制御理論の代表的な応用例です。
電力システムと自動車制御への応用
電力システムの安定化制御においても、利得行列は重要な役割を担います。
発電機の励磁制御（AVR）・電力系統安定化装置（PSS）・FACTS（柔軟交流送電システム）などの制御系設計において、多変数制御のための利得行列設計が不可欠です。
自動車のアクティブサスペンション制御・電動パワーステアリング・横すべり防止装置（ESC）なども、状態フィードバック制御と利得行列に基づいた制御システムの代表例です。
| 応用分野 | 状態変数の例 | 制御入力の例 | 利得行列設計手法 |
|---|---|---|---|
| ロボットアーム制御 | 関節角度・角速度 | 各関節モータートルク | PD制御・LQR |
| 航空機姿勢制御 | ピッチ・ロール・ヨー角、角速度 | 舵面偏向角 | LQR・H∞制御 |
| 倒立振子安定化 | 振子角度・角速度・台車位置 | 台車駆動力 | 極配置法・LQR |
| 電力系統安定化 | 発電機回転数・電圧 | 励磁電圧・PSS出力 | LQR・ロバスト制御 |
| 自動車横すべり制御 | ヨーレート・横加速度・舵角 | 各輪ブレーキ力 | 状態フィードバック・MPC |
このように、利得行列を活用した制御系は現代の自動化・メカトロニクスシステムの至るところに応用されています。
まとめ
この記事では、利得行列の基本定義から状態空間表現との関係・極配置法・LQRによる最適設計・具体的な計算例・多様な応用分野まで幅広く解説しました。
利得行列とは状態フィードバック制御において制御入力と状態変数を結びつける行列であり、その設計がシステムの安定性・応答性・ロバスト性を左右する制御工学の核心概念です。
極配置法・LQR・リカッチ方程式などの設計手法を理解することで、多様な制御システムに対して体系的に利得行列を設計できる力が身につくでしょう。
ロボット・航空宇宙・自動車・電力システムなど広範な応用分野での利得行列の活用を知ることで、制御理論の実用的な重要性をより深く実感できたのではないでしょうか。
