マトリックスS-手順を用いたデータ駆動型構造化コントローラ設計

基本情報

• 論文ID: 2503.14949
• タイトル: Data-Driven Structured Controller Design Using the Matrix S-Procedure
• 著者: Zhaohua Yang, Yuxing Zhong, Nachuan Yang, Xiaoxu Lyu, Ling Shi
• 分類: math.OC（数学最適化と制御）
• 発表日時: 2025年10月14日（Automaticaに投稿されたプレプリント）
• 論文リンク: https://arxiv.org/abs/2503.14949

要約

本論文は、$H_2$性能と$H_\infty$性能を考慮した離散時間線形時不変（LTI）システムのデータ駆動型最適構造化コントローラ設計に焦点を当てている。具体的には3つのシナリオを研究する：(i)モデルベースの構造化制御、(ii)データ駆動型の非構造化制御、(iii)データ駆動型の構造化制御。$H_2$性能については、シナリオ(ii)と(iii)を主に研究する。これはシナリオ(i)が既に文献で広く研究されているためである。$H_\infty$性能については、3つのシナリオすべてを考察する。構造化制御に対しては、線形化技術を導入し、元の非凸問題を半定値計画法（SDP）問題に変換し、これに基づいて反復線形行列不等式（ILMI）アルゴリズムを開発した。

研究背景と動機

問題背景

1. ネットワーク化制御システムの課題：大規模ネットワーク化制御システムにおいて、集中型制御戦略は膨大な通信負荷をもたらすため、疎な構造を持つコントローラの設計が必要である。
2. 構造化コントローラ設計の困難性：構造制約を伴うコントローラ設計はNP困難な複雑性を持ち、既存の方法はLMIベースと勾配ベースの方法に大別される。
3. データ駆動型制御の必要性：実際の応用では、多くのシステムの動力学モデルが得られないことが多く、これがデータ駆動型制御の研究を促進している。

既存方法の限界

1. 従来的方法は既知モデルに依存：既存の構造化コントローラ設計方法はシステムモデルが既知であることを仮定している。
2. データ駆動型方法の保守性：Miller等15の既存データ駆動型方法には以下の問題がある：
   • 瞬時有界ノイズをエネルギー有界ノイズ列に過度に近似し、追加の保守性を導入する
   • 非凸実行可能領域の凸部分集合内でのみ構造化コントローラを設計する
   • 性能単調性の保証が欠ける

核心的貢献

1. 改善されたデータ駆動型フレームワーク：
   • すべての可能なシステム行列の集合を各データ点集合の交集合として表現し、15の集合より小さい集合を得る
   • マトリックスS-手順を利用してすべての可能なシステム行列を処理し、実行可能領域が15の実行可能領域を包含することを証明する
   • データ長に対する性能の単調非減少性を証明する
2. 新規な線形化技術：
   • 元の非凸問題をSDP問題に変換する線形化技術を開発する
   • 構造化コントローラを計算するためのILMIアルゴリズムを提案する
   • 15より大きな実行可能領域を探索できる
3. 包括的な$H_\infty$制御研究：
   • モデルベースの構造化制御、データ駆動型非構造化制御、データ駆動型構造化制御を初めて包括的に研究する

方法の詳細

タスク定義

離散時間LTIシステムを考える：

x_{k+1} = Ax_k + Bu_k + Gd_k
y_k = Cx_k + Du_k + Hd_k

ここで$x_k \in \mathbb{R}^{n_x}$はシステム状態、$u_k \in \mathbb{R}^{n_u}$は制御入力、$y_k \in \mathbb{R}^{n_y}$は制御出力、$d_k \in \mathbb{R}^{n_d}$は外部外乱である。

状態フィードバックコントローラ$u_k = Kx_k$を使用する。ここで$K \in \mathbb{R}^{n_u \times n_x}$は構造制約$K \in S$を満たす必要がある。

データ収集メカニズム

行列$A, B$は未知だが確定的であり、行列$C, D, G, H$は既知であると仮定する。収集されたデータは以下を含む：

X = [x_0, x_1, ..., x_T] ∈ ℝ^{n×(T+1)}
U = [u_0, u_1, ..., u_{T-1}] ∈ ℝ^{m×T}

重要な仮定：プロセスノイズは有界である、すなわち$\|w_k\|_2 \leq \epsilon$。

核心的技術革新

1. 改善されたシステム行列集合表現

第$i$番目の時間ステップのデータ$(x_i, u_i, x_{i+1})$に基づいて、以下を定義する：

Σ_i = {(A,B) | x_{i+1} = Ax_i + Bu_i + w_i, w_i w_i^T ⪯ ε²I}

データと一致するすべてのシステム行列の集合は：

Σ = ∩_{i=0}^{T-1} Σ_i

これは15の集合より正確であり、保守性を低減する。

2. マトリックスS-手順の適用

$H_2$制御に対して、マトリックスS-手順を利用すると、条件(5a)がすべての$(A,B) \in Σ$に対して成立するための十分条件は、非負の$α_0, α_1, ..., α_{T-1}$と$β > 0$が存在して以下が成立することである：