2025-11-24T10:10:18.436123

Computing the Hard Scaled Relative Graph of LTI Systems

Krebbekx, Baron-Prada, Tóth et al.
Scaled Relative Graphs (SRGs) provide a novel graphical frequency-domain method for the analysis of nonlinear systems, where Linear Time-Invariant (LTI) systems are the fundamental building block. To analyze feedback loops with unstable LTI components, the hard SRG is required, since it aptly captures the input/output behavior on the extended $L_2$ space. In this paper, we develop a systematic computational method to exactly compute the hard SRG of LTI systems, which may be unstable and contain integrators. We also study its connection to the Nyquist criterion, including the multivariable case, and demonstrate our method on several examples.
academic

LTI システムのハード スケール相対グラフの計算

基本情報

  • 論文ID: 2511.17297
  • タイトル: Computing the Hard Scaled Relative Graph of LTI Systems
  • 著者: Julius P. J. Krebbekx, Eder Baron-Prada, Roland Tóth, Amritam Das
  • 所属機関: アイントホーフェン工科大学、オーストリア科学技術研究所、ETH チューリッヒ、HUN-REN コンピュータ科学・制御研究所
  • 分類: eess.SY（システム・制御）、cs.SY、math.OC
  • 提出日: 2025年11月21日
  • 論文リンク: https://arxiv.org/abs/2511.17297

概要

スケール相対グラフ（SRG）は、非線形システム分析のための新しい周波数領域グラフィカル手法を提供し、線形時不変（LTI）システムが基本的な構成要素となります。不安定なLTI成分を含むフィードバックループを分析するには、拡張L₂空間上の入出力動作を適切に捉えることができるハードSRGを使用する必要があります。本論文は、不安定である可能性があり、積分器を含むLTIシステムのハードSRGを正確に計算するための体系的な計算方法を開発しています。また、多変数の場合を含むナイキスト判定法との関連性を調査し、複数の例を通じて本方法を実証しています。

研究背景と動機

研究課題

本論文は、非線形システム分析における中核的な問題を解決することを目指しています。すなわち、不安定なLTI成分を含むフィードバックシステムのスケール相対グラフ（SRG）を体系的に計算する方法です。

問題の重要性

  1. 古典的手法の限界: ナイキスト図はLTIシステム分析の基礎ですが、非線形システム分析への体系的な一般化は困難です
  2. 実務的需要: 産業応用では、不安定成分を含むフィードバックシステムの分析が頻繁に必要とされます
  3. 理論的完全性: 既存のSRGフレームワークには、不安定システムを処理するための体系的な方法が不足しています

既存手法の限界

  1. ソフトSRG: 最初に3で提案され、安定なSISO システムのみに適用可能で、不安定システムには問題があります
  2. 拡張SRG: 11でナイキスト判定法情報を追加することでSISO不安定システムの問題を解決しましたが、単変量の場合に限定されます
  3. ハードSRG: 12で定義され、理論的には不安定システムを処理できますが、体系的な計算方法が不足しています

研究動機

ハードSRGの主な利点はフィードバックループ内の不安定システムを処理することにあるため、不安定なLTIシステムのハードSRGを計算する方法を開発することは、非線形システムのSRG分析を大幅に進展させます。

中核的貢献

  1. 体系的計算方法: 不安定システムと積分器を含むシステムを含む、LTIシステムのハードSRGを正確に計算できる非近似アルゴリズム（定理2）を提案
  2. 理論的等価性証明: SISO システムについて、ハードSRGと拡張SRGが完全に等価であることを証明（定理3）
  3. MIMO拡張: 多変数（MIMO）の場合に方法を拡張し、スクエアシステムに対する完全な計算フレームワークを提供
  4. ナイキスト判定法との関連性: ハードSRGがMIMO ナイキスト判定法と一般化ナイキスト判定法（GNC）の代替方法としてどのように機能するかを実証
  5. 実用的アルゴリズム: 伝達関数表現に基づき、複素平面上の交差と実軸を中心とした円盤の除去により境界を構築

方法の詳細

タスク定義

LTI システム G: L²ₑᵖ → L²ₑᵍ が与えられ、その伝達関数が G(s) ∈ Rᵍˣᵖ(s) である場合、そのハードSRGを計算します：

SRGe(G)=⋃T>0SRG(G∣T)\text{SRG}_e(G) = \bigcup_{T>0} \text{SRG}(G|_T)

ここで G|_T は時間区間0,T上のGの切断演算子です。

中核的理論フレームワーク

1. 双曲凸性表現

任意のh-凸集合Cの閉包は以下のように表現できます：

cl C=⋂α∈R(DRα(α)∖Drα(α))=:G({rα},{Rα})\text{cl}\,C = \bigcap_{\alpha \in \mathbb{R}} \left(D_{R_\alpha}(\alpha) \setminus D_{r_\alpha}(\alpha)\right) =: \mathcal{G}(\{r_\alpha\}, \{R_\alpha\})

ここで Dr(x)D_r(x) は中心がx、半径がrの閉円盤を表します。

2. シフトシステムの定義

α ∈ ℝ に対して、シフトシステムを定義します：

Gα=(G0(n−q)×p)−(αI0(n−p)×p)G_\alpha = \begin{pmatrix} G \\ 0_{(n-q)\times p} \end{pmatrix} - \begin{pmatrix} \alpha I \\ 0_{(n-p)\times p} \end{pmatrix}

ここで n = max{p,q}。

アルゴリズムの中核成分

命題2: 最大ゲイン半径Rαの計算

\bar{\sigma}(G_\alpha) & \text{if } G \in RH_\infty^{q\times p} \\ \infty & \text{otherwise} \end{cases}$$ **主要な考え方**: - 安定システムの場合、H∞ノルム（最大特異値）を使用 - 不安定または非正則システムの場合、Rα = ∞ #### 命題3: 最小ゲイン半径rαの計算 $$r_\alpha = \begin{cases} \underline{\sigma}(G_\alpha) & \text{if } G_\alpha \text{ is minimum-phase} \\ 0 & \text{otherwise} \end{cases}$$ **主要な考え方**: - 逆システム $G_\alpha^{-1}$ を求めて命題2を適用 - Gαが最小位相（すべての伝達零点が開左半平面にある）の場合のみ、rαは非ゼロ - Smith-McMillan形式を利用：Gαの極点は $G_\alpha^{-1}$ の零点であり、その逆も成立 ### メインアルゴリズム（定理2） **スクエアLTI システム** G: L²ₑᵖ → L²ₑᵖ、伝達関数 G(s) ∈ Rᵖˣᵖ(s) に対して： $$\text{SRG}_e(G) = \mathcal{G}(\{r_\alpha\}, \{R_\alpha\})$$ **アルゴリズムステップ**: 1. 各α ∈ ℝ に対して、Gα = G - αI を構築 2. 命題2を使用してRαを計算 3. 命題3を使用してrαを計算 4. すべての環形領域の交差を計算 ### 技術的革新点 #### 1. 拡張L₂空間の処理 - 従来の方法はL₂空間で機能し、発散信号を処理できません - ハードSRGは拡張空間L²ₑで定義され、周期的および発散信号を含みます - すべてのT > 0に対して切断信号を考慮することで不安定性を処理 #### 2. 最小位相検出 - 主要な革新：rαの計算をGαが最小位相であるかどうかの検出に帰着させる - 最小位相 ⟺ すべての伝達零点が開左半平面にある - 非最小位相システムのrα = 0は、フィードバック内で引き起こす可能性のある問題を反映 #### 3. ナイキスト判定法との統一 SISO システムについて、以下を証明しました： $$\text{SRG}_e(G) = \text{SRG}'(G)$$ ここで SRG'(G) は拡張SRGであり、ソフトSRGとナイキスト判定法情報を含みます。 #### 4. 特殊ケースの処理 - **積分器**: 自動的にRα = ∞をもたらす - **非正則システム**: Smith-McMillan形式によるランク不足の検出 - **時間遅延**: rα = 0をもたらし、Rαは影響を受けない ## 実験設定 ### サンプルシステム #### SISO システム 1. **G₁(s) = 1/(s² + s + 1)**: 安定な2次システム 2. **G₂(s) = 1/(s(s+1))**: 積分器を含むシステム 3. **G₃(s) = 3/((s-2)(s/10+1))**: 不安定システム（右半平面極点） #### MIMO システム 1. **G₄(s)**: 2×2安定システム $$G_4(s) = \begin{pmatrix} \frac{1}{s+1} & \frac{1}{s+2} \\ \frac{1}{s+4} & \frac{1}{s+3} \end{pmatrix}$$ 2. **G₅(s)**: 2×2不安定システム $$G_5(s) = \begin{pmatrix} \frac{s+7}{s-1} & \frac{s-5}{(s+2)^2} \\ \frac{1}{(s+4)^3} & \frac{s}{(s+3)^2} \end{pmatrix}$$ 3. **G₆(s)**: 積分器を含む3×3システム ### 可視化方法 - 複素平面上にハードSRG境界（灰色領域）をプロット - ソフトSRG境界（濃い灰色領域）を重ねる - -1点などの重要な点にマーク（安定性分析用） ## 実験結果 ### 主要な結果 #### SISO システムの結果 1. **G₁**: ハードSRGとソフトSRGが一致（システムが安定のため）、有界領域を形成 2. **G₂**: 積分器のため、ハードSRGは無限に延伸し、半実軸 [0,∞) を含む 3. **G₃**: 不安定システムのハードSRGは追加領域を含み、ナイキスト囲み情報を反映 **拡張SRG等価性の検証**: - G₂について、$\mathcal{N}_{G_2} \cap \mathbb{R} = [0,\infty)$ - これは図2bのハードSRGと完全に一致 - 定理3の正確性を検証 #### MIMO システムの結果 **G₄（安定な2×2システム）**: - ハードSRGとソフトSRGの境界は非常に接近 - 両者とも-1点から分離しており、フィードバック安定性を示唆 **G₅（不安定な2×2システム）**: - ハードSRGは**2つの穴**を表示（図2e） - これらの穴は安定な閉ループ構成領域に対応 - より小さい穴はゲイン設計に利用可能 **G₆（積分器を含む3×3システム）**: - ハードSRGは部分的に無界（積分器のため） - ソフトSRG境界はIm方向で±j∞に向かう傾向 ### 安定性分析の応用 #### ゲイン設計例（G₅に基づく） 論文は3つのゲイン選択戦略を実証しています： 1. **純粋スケーリング**: k₁ = 0, k₂ = 1/5 - 区間[-6, 0) ⊈ SRGₑ(G₅) - スケーリング後、-1点はSRGから分離 2. **反転とスケーリング**: k₁ = 0, k₂ = -2 - 右側の小さい穴を-1点上に反転 - ハードSRGの幾何学的特性を利用 3. **平行移動**: k₁ = -1.5, k₂ = 1 - 小さい穴を左に移動して-1点を含める - 閉ループ安定性を保証 各ケースは安定性条件を満たします： $$\text{dist}(-1, \text{SRG}_e(k_1I + k_2G_5)) \geq r > 0$$ ### 主要な発見 1. **精度**: SISO および正規行列MIMO システムについて、境界は正確（非近似） 2. **数値安定性**: 有限αセットに起因する数値境界は常に外部近似であり、SRGを過小評価しません 3. **幾何学的直感性**: ハードSRGは直感的な幾何学的解釈を提供し、安定性マージンを可視化 4. **モジュール性**: GのハードSRGを計算すると、異なるゲイン設計に再利用でき、再計算は不要 ## 関連研究 ### SRGフレームワークの発展 1. **オリジナルSRG** [2,3]: Ryu等およびChaffey等により提案、L₂空間に基づくソフトSRG 2. **拡張SRG** [4,11]: ナイキスト判定法情報を追加することで不安定なSISO システムに拡張 3. **ハードSRG** [12]: Chenにより提案、L₂[0,T]上で定義、理論的には不安定システムを処理可能 ### SRG応用分野 1. **非線形分析**: 一般化円判定法[4]、小ゲイン定理[3] 2. **リセット制御**: 分析[8]および設計[9] 3. **分散システム**: 安定性条件[10] 4. **周波数領域特性**: 非線形ボード線図と帯域幅定義[5] 5. **位相分析**: 位相進み/遅れ概念[6] ### 古典的手法との関係 1. **ナイキスト判定法**: ハードSRGはSISO の場合、拡張SRGと等価であり、後者はナイキスト情報を含む 2. **MIMO ナイキスト判定法** [13]: ハードSRGは代替方法を提供し、巻き数の計算が不要 3. **一般化ナイキスト判定法（GNC）** [14]: ハードSRGは特性軌跡の手動接合の問題を回避 ### 本論文の利点 1. **体系性**: 発見的方法ではなく、明確なアルゴリズムを提供 2. **汎用性**: 安定、不安定、積分器を含むシステムを処理 3. **精度**: 非近似方法、正確な境界を提供 4. **拡張性**: MIMO の場合に自然に拡張 ## 結論と考察 ### 主要な結論 1. **アルゴリズムの完全性**: LTI システム（不安定システムを含む）のハードSRGを正確に計算するための初の体系的方法を開発 2. **理論的統一**: SISO の場合、ハードSRGと拡張SRGの等価性を証明し、2つのフレームワークを統一 3. **実用的価値**: ハードSRGはMIMO ナイキスト判定法およびGNCの実用的代替案を提供し、以下の利点があります： - 巻き数の計算が不要 - ゲイン変化時にSRGを再利用可能 - L₂ゲイン界と安定性マージンを直接提供 4. **幾何学的直感性**: 円盤交差に基づく幾何学的構築は、明確な可視化と直感的理解を提供 ### 限界 #### 1. モデル依存性 - 現在の方法は伝達関数表現に基づき、システムのパラメトリックモデルが必要 - ナイキスト判定法の利点は周波数応答測定データを直接使用できることにあります - **将来の需要**: 入出力データからrαおよびRαを直接計算 #### 2. 最小位相検出の課題 - 命題2は不安定極点の存在を検出するだけで（計数不要） - 命題3は非最小位相零点を検出する必要があります - **既知の困難**: データから非最小位相零点を検出することは課題です[22] #### 3. 非スクエアシステム - 現在の完全な理論はスクエアシステム（p = q）のみに適用 - ワイド場合（p > q）: Gαは常に非自明な代数核を持ち、rα = 0をもたらす - トール場合（p < q）: 左疑似逆を使用可能ですが、多くの追加技術詳細が必要 #### 4. 適切性の仮定 - ハードSRG安定性分析は、フィードバックシステムの適切性（well-posedness）の仮定が必要 - ソフトおよび拡張SRGはゲイン界の他に適切性を保証 - これはハードSRGの拡張SRGに対する相対的な欠点 ### 将来の方向 #### 1. データ駆動方法 周波数応答または入出力データからハードSRGを直接計算し、特に： - データから不安定極点を検出 - データから非最小位相零点を識別 #### 2. 非スクエアシステムの拡張 - トール システム（p < q）の理論を完成 - ワイド システムに適用可能な代替方法を開発 #### 3. 非線形システムの応用 - 非線形MIMO システム分析に結果を適用[7] - [12]のハード SRG非線形フィードバック分析フレームワークと結合 #### 4. 計算最適化 - 効率的な数値アルゴリズムを開発 - 高次元MIMO システムの計算複雑性に対処 - 精度と計算コストのバランスを取るための適応的なα集合選択 #### 5. ロバスト性分析 - 不確実性のあるシステムに拡張 - ロバスト制御理論と結合 ## 深い評価 ### 利点 #### 1. 理論的厳密性 - **完全な数学的証明**: すべての主要結果（命題1-3、定理2-3）に厳密な証明がある - **幾何学的基礎**: 双曲凸性に基づく堅実な理論的基礎 - **等価性証明**: 定理3は古典的手法との明確な関連性を確立 #### 2. 方法の革新性 - **統一フレームワーク**: 安定および不安定なLTI システムを処理するための統一計算方法を初めて提供 - **巧妙な逆システムトリック**: 逆Gαを求めてrαを計算することは優雅な解決策 - **幾何学的構築**: 円盤交差方法は直感的かつ実装しやすい #### 3. 実用的価値 - **アルゴリズムの実装可能性**: 方法は標準的な線形システムツール（特異値、最小位相検出）に基づく - **モジュール設計**: 計算されたSRGは異なるゲイン設計に再利用可能 - **直感的可視化**: 複素平面上の幾何学的表現は理解と応用を容易にする #### 4. 例の充分性 - SISO および MIMO、安定および不安定、積分器を含むなど複数のケースをカバー - ゲイン設計応用を実際に実証 - 可視化は明確で効果的 ### 不足 #### 1. 理論的限界 - **スクエア制限**: 主要結果はスクエアシステムのみに適用可能で、応用範囲を制限 - **適切性**: 拡張SRGのように適切性を保証できない - **正規行列仮定**: 命題1について、一般的なMIMO システムでは境界が不正確である可能性 #### 2. 実践的課題 - **モデル依存**: 正確な伝達関数モデルが必要 - **非最小位相検出**: データからの検出は課題 - **計算複雑性**: 高次元MIMO システムについて、連続的なαに対する計算が必要 #### 3. 実験の不足 - **計算時間分析の欠如**: アルゴリズムの計算効率が報告されていない - **数値精度の議論**: 有限αセットの数値誤差が定量化されていない - **GNCとの実際の比較**: 理論的比較のみで、実際のケース比較が不足 #### 4. 提示の問題 - **記号が密集**: 多くの数学記号は可読性に影響する可能性 - **アルゴリズム疑似コードの欠如**: 方法は明確ですが、直接実装可能なアルゴリズム説明が不足 - **図の説明**: 一部の図（図2bの数値アーティファクトなど）はより多くの説明が必要 ### 影響力評価 #### 分野への貢献 1. **理論的完全性**: SRGフレームワークにおける不安定システム処理の空白を埋める 2. **方法論的進歩**: 非線形システム分析のための新しいツールを提供 3. **古典理論の現代化**: ナイキスト判定法に対する現代的な幾何学的解釈を提供 #### 実用的価値 1. **制御設計**: 不安定成分を含むフィードバックシステムの設計に直接適用可能 2. **安定性分析**: 従来の方法より直感的な安定性マージン可視化を提供 3. **教育的価値**: 幾何学的方法はフィードバックシステムの動作理解を支援 #### 再現性 - **理論の再現性**: すべての証明は付録に完全に記載 - **実装の可行性**: 標準ツールに基づき、実装が容易 - **例の明確性**: 具体的な伝達関数を提供し、検証が容易 - **コードの欠如**: オープンソース実装が提供されていない（残念） ### 適用シーン #### 最適な応用 1. **複雑なフィードバックシステム設計**: 特に不安定成分を含むシステム 2. **ロバスト性分析**: 安定性マージンの可視化が必要な場合 3. **多変数システム**: MIMO システムのゲイン調整と安定性分析 4. **非線形システム**: 非線形SRG分析の基本的な構成要素として #### あまり適さないシーン 1. **純粋なデータ駆動**: 正確なモデルが得られない場合 2. **非スクエアシステム**: p ≠ q の場合、方法が不完全 3. **リアルタイム応用**: 計算複雑性が重要な制約の場合 4. **極めて高次元システム**: 次元が高い場合、計算が困難になる可能性 ### 総合評価 これは**高品質な理論制御論文**であり、スケール相対グラフフレームワークに重要な貢献をしています。主な利点は： - 不安定なLTI システムのハードSRGを処理するための初の体系的方法を提供 - 理論が厳密で証明が完全 - 古典的なナイキスト判定法との明確な関連性を確立 - 非線形システム分析の基礎を構築 主な改善の余地： - 非スクエアシステムへの拡張 - データ駆動版の開発 - オープンソース実装とより多くの実際のケースの提供 **推奨指数**: ⭐⭐⭐⭐☆ (4/5) この研究は、非線形制御、ロバスト制御、システム分析に従事する研究者にとって重要な参考価値があり、特にSRG方法に関心がある、または不安定システムを処理する必要がある学者にとって有用です。 ## 参考文献（精選） 1. **[2] Ryu et al., 2022**: "Scaled relative graphs: Nonexpansive operators via 2D Euclidean geometry" - SRGの原始的な数学的定義 2. **[3] Chaffey et al., 2023**: "Graphical Nonlinear System Analysis" - SRGを制御理論に導入 3. **[11] Krebbekx et al., 2025**: "Scaled Relative Graph Analysis of General Interconnections" - 拡張SRGの定義 4. **[12] Chen, 2025**: "Soft and Hard Scaled Relative Graphs" - ハードSRGの初の定義 5. **[13] MacFarlane & Postlethwaite, 1977**: MIMO ナイキスト判定法の古典的文献 6. **[14] Desoer & Wang, 1980**: 一般化ナイキスト判定法（GNC）