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

Computing the Hard Scaled Relative Graph of LTI Systems

基本信息

  • 论文ID: 2511.17297
  • 标题: Computing the Hard Scaled Relative Graph of LTI Systems
  • 作者: Julius P. J. Krebbekx, Eder Baron-Prada, Roland Tóth, Amritam Das
  • 机构: Eindhoven University of Technology, Austrian Institute of Technology, ETH Zurich, HUN-REN Institute for Computer Science and Control
  • 分类: eess.SY (Systems and Control), cs.SY, math.OC
  • 提交时间: 2025年11月21日
  • 论文链接: https://arxiv.org/abs/2511.17297

摘要

Scaled Relative Graphs (SRGs) 为非线性系统分析提供了一种新颖的图形化频域方法，其中线性时不变（LTI）系统是基本构建模块。为了分析包含不稳定LTI组件的反馈回路，需要使用hard SRG，因为它能够恰当地捕获扩展L₂空间上的输入/输出行为。本文开发了一种系统化的计算方法来精确计算LTI系统的hard SRG，这些系统可能是不稳定的并包含积分器。文章还研究了其与Nyquist准则的联系，包括多变量情况，并通过多个实例演示了该方法。

研究背景与动机

研究问题

本文旨在解决非线性系统分析中的核心问题：如何系统地计算包含不稳定LTI组件的反馈系统的Scaled Relative Graph (SRG)。

问题重要性

  1. 经典方法的局限: Nyquist图是LTI系统分析的基石，但难以系统地推广到非线性系统分析
  2. 实际需求: 工业应用中经常需要分析包含不稳定组件的反馈系统
  3. 理论完整性: 现有SRG框架缺乏处理不稳定系统的系统化方法

现有方法局限

  1. Soft SRG: 最初在3中提出，仅适用于稳定的SISO系统，对不稳定系统存在问题
  2. Extended SRG: 在11中通过添加Nyquist准则信息解决了SISO不稳定系统问题，但仅限于单变量情况
  3. Hard SRG: 在12中定义，理论上可以处理不稳定系统，但缺乏系统的计算方法

研究动机

由于hard SRG的主要优势在于处理反馈回路中的不稳定系统，开发计算不稳定LTI系统hard SRG的方法将显著推进非线性系统的SRG分析。

核心贡献

  1. 系统化计算方法: 提出了一种非近似的算法（Theorem 2），可以精确计算LTI系统的hard SRG，包括不稳定系统和含积分器的系统
  2. 理论等价性证明: 证明了对于SISO系统，hard SRG与extended SRG完全等价（Theorem 3）
  3. MIMO扩展: 将方法扩展到多变量（MIMO）情况，为square系统提供了完整的计算框架
  4. 与Nyquist准则的联系: 展示了hard SRG如何作为MIMO Nyquist准则和广义Nyquist准则（GNC）的替代方法
  5. 实用算法: 基于传递函数表示，通过在复平面上相交和移除以实轴为中心的圆盘来构造边界

方法详解

任务定义

给定一个LTI系统 G: L²ₑᵖ → L²ₑᵍ，其传递函数为 G(s) ∈ Rᵍˣᵖ(s)，计算其hard SRG：

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

其中 G|_T 是G在时间区间0,T上的截断算子。

核心理论框架

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}。

算法核心组件

Proposition 2: 最大增益半径Rα的计算

\bar{\sigma}(G_\alpha) & \text{if } G \in RH_\infty^{q\times p} \\ \infty & \text{otherwise} \end{cases}$$ **关键思想**: - 对于稳定系统，使用H∞范数（最大奇异值） - 对于不稳定或非正则系统，Rα = ∞ #### Proposition 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}$ 并应用Proposition 2 - 只有当Gα是最小相位（所有传输零点在开左半平面）时，rα才非零 - 利用Smith-McMillan形式：Gα的极点是 $G_\alpha^{-1}$ 的零点，反之亦然 ### 主算法（Theorem 2） **对于square LTI系统** G: L²ₑᵖ → L²ₑᵖ，传递函数 G(s) ∈ Rᵖˣᵖ(s)： $$\text{SRG}_e(G) = \mathcal{G}(\{r_\alpha\}, \{R_\alpha\})$$ **算法步骤**: 1. 对每个α ∈ ℝ，构造 Gα = G - αI 2. 使用Proposition 2计算Rα 3. 使用Proposition 3计算rα 4. 计算所有环形区域的交集 ### 技术创新点 #### 1. 扩展L₂空间的处理 - 传统方法在L₂空间上工作，不能处理发散信号 - Hard SRG在扩展空间L²ₑ上定义，包含周期和发散信号 - 通过对所有T > 0考虑截断信号来处理不稳定性 #### 2. 最小相位检测 - 关键创新：将rα的计算归结为检测Gα是否为最小相位 - 最小相位 ⟺ 所有传输零点在开左半平面 - 非最小相位系统的rα = 0，反映了其在反馈中可能引起的问题 #### 3. 与Nyquist准则的统一 对于SISO系统，证明了： $$\text{SRG}_e(G) = \text{SRG}'(G)$$ 其中SRG'(G)是extended SRG，包含soft SRG和Nyquist准则信息。 #### 4. 处理特殊情况 - **积分器**: 自动导致Rα = ∞ - **非正则系统**: 通过Smith-McMillan形式检测秩缺失 - **时延**: 导致rα = 0，Rα不受影响 ## 实验设置 ### 示例系统 #### SISO系统 1. **G₁(s) = 1/(s² + s + 1)**: 稳定的二阶系统 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包含积分器的系统 ### 可视化方法 - 在复平面上绘制hard SRG边界（灰色区域） - 叠加soft SRG边界（深灰色区域） - 标记关键点如-1点（用于稳定性分析） ## 实验结果 ### 主要结果 #### SISO系统结果 1. **G₁**: Hard SRG与soft SRG一致（因为系统稳定），形成有界区域 2. **G₂**: 由于积分器，hard SRG延伸至无穷，包含半实轴 [0,∞) 3. **G₃**: 不稳定系统的hard SRG包含额外区域，反映Nyquist封闭信息 **验证Extended SRG等价性**: - 对于G₂，$\mathcal{N}_{G_2} \cap \mathbb{R} = [0,\infty)$ - 这与图2b中的hard SRG完全一致 - 验证了Theorem 3的正确性 #### MIMO系统结果 **G₄（稳定2×2系统）**: - Hard SRG与soft SRG的边界非常接近 - 两者都与-1点分离，表明反馈稳定 **G₅（不稳定2×2系统）**: - Hard SRG显示**两个孔洞**（图2e） - 这些孔洞对应于稳定的闭环配置区域 - 较小的孔洞可用于设计增益 **G₆（含积分器的3×3系统）**: - Hard SRG部分无界（由于积分器） - Soft SRG边界在Im方向趋向±j∞ ### 稳定性分析应用 #### 增益设计示例（基于G₅） 文章展示了三种增益选择策略： 1. **纯缩放**: k₁ = 0, k₂ = 1/5 - 区间[-6, 0) ⊈ SRGₑ(G₅) - 缩放后-1点与SRG分离 2. **翻转和缩放**: k₁ = 0, k₂ = -2 - 将右侧小孔洞翻转到-1点上 - 利用hard SRG的几何特性 3. **平移**: k₁ = -1.5, k₂ = 1 - 将小孔洞左移使其包含-1点 - 保证闭环稳定性 每种情况都满足稳定性条件： $$\text{dist}(-1, \text{SRG}_e(k_1I + k_2G_5)) \geq r > 0$$ ### 关键发现 1. **精确性**: 对于SISO和normal矩阵MIMO系统，边界是精确的（非近似） 2. **数值稳定性**: 有限α集合导致的数值边界总是外逼近，不会低估SRG 3. **几何直观性**: Hard SRG提供了直观的几何解释，可视化稳定性裕度 4. **模块化**: 一旦计算了G的hard SRG，可以重复用于不同增益设计，无需重新计算 ## 相关工作 ### SRG框架发展 1. **原始SRG** [2,3]: 由Ryu等人和Chaffey等人提出，基于L₂空间的soft SRG 2. **Extended SRG** [4,11]: 通过添加Nyquist准则信息扩展到不稳定SISO系统 3. **Hard SRG** [12]: Chen提出，在L₂[0,T]上定义，理论上可处理不稳定系统 ### SRG应用领域 1. **非线性分析**: 广义圆准则[4]、小增益定理[3] 2. **复位控制**: 分析[8]和设计[9] 3. **分散系统**: 稳定性条件[10] 4. **频域特性**: 非线性Bode图和带宽定义[5] 5. **相位分析**: 相位超前/滞后概念[6] ### 与经典方法的关系 1. **Nyquist准则**: Hard SRG在SISO情况下等价于extended SRG，后者包含Nyquist信息 2. **MIMO Nyquist准则** [13]: Hard SRG提供替代方法，无需计算围绕次数 3. **广义Nyquist准则（GNC）** [14]: Hard SRG避免了特征轨迹的手工拼接问题 ### 本文优势 1. **系统化**: 提供了明确的算法，而非启发式方法 2. **通用性**: 处理稳定、不稳定、含积分器的系统 3. **精确性**: 非近似方法，提供精确边界 4. **可扩展性**: 自然扩展到MIMO情况 ## 结论与讨论 ### 主要结论 1. **算法完整性**: 开发了首个系统化方法来精确计算LTI系统（包括不稳定系统）的hard SRG 2. **理论统一**: 证明了SISO情况下hard SRG与extended SRG的等价性，统一了两种框架 3. **实用价值**: Hard SRG为MIMO Nyquist准则和GNC提供了实用替代方案，具有以下优势： - 无需计算围绕次数 - 增益变化时可重用SRG - 直接提供L₂增益界和稳定性裕度 4. **几何直观性**: 基于圆盘相交的几何构造提供了清晰的可视化和直观理解 ### 局限性 #### 1. 模型依赖性 - 当前方法基于传递函数表示，需要系统的参数模型 - Nyquist准则的优势在于可直接使用频率响应测量数据 - **未来需求**: 从输入/输出数据直接计算rα和Rα #### 2. 非最小相位检测挑战 - Proposition 2仅需检测是否存在不稳定极点（不需要计数） - Proposition 3需要检测非最小相位零点 - **已知困难**: 从数据检测非最小相位零点具有挑战性[22] #### 3. 非square系统 - 当前完整理论仅适用于square系统（p = q） - Wide情况（p > q）: Gα总有非平凡代数核，导致rα = 0 - Tall情况（p < q）: 可使用左伪逆，但需要大量额外技术细节 #### 4. 适定性假设 - Hard SRG稳定性分析需要假设反馈系统的适定性（well-posedness） - Soft和extended SRG除了增益界外还能保证适定性 - 这是hard SRG相对于extended SRG的一个劣势 ### 未来方向 #### 1. 数据驱动方法 从频率响应或输入/输出数据直接计算hard SRG，特别是： - 从数据检测不稳定极点 - 从数据识别非最小相位零点 #### 2. 非square系统扩展 - 完善tall系统（p < q）的理论 - 开发适用于wide系统的替代方法 #### 3. 非线性系统应用 - 将结果应用于非线性MIMO系统分析[7] - 结合[12]中的hard SRG非线性反馈分析框架 #### 4. 计算优化 - 开发高效的数值算法 - 处理高维MIMO系统的计算复杂性 - 自适应选择α集合以平衡精度和计算成本 #### 5. 鲁棒性分析 - 扩展到不确定系统 - 结合鲁棒控制理论 ## 深度评价 ### 优点 #### 1. 理论严谨性 - **完整的数学证明**: 所有主要结果（Propositions 1-3, Theorems 2-3）都有严格证明 - **几何基础**: 基于双曲凸性的坚实理论基础 - **等价性证明**: Theorem 3建立了与经典方法的明确联系 #### 2. 方法创新性 - **统一框架**: 首次提供了处理稳定和不稳定LTI系统的统一计算方法 - **巧妙的逆系统技巧**: 通过求逆Gα来计算rα是一个优雅的解决方案 - **几何构造**: 圆盘相交方法既直观又易于实现 #### 3. 实用价值 - **算法可实现**: 方法基于标准的线性系统工具（奇异值、最小相位检测） - **模块化设计**: 计算的SRG可重用于不同增益设计 - **直观可视化**: 复平面上的几何表示便于理解和应用 #### 4. 示例充分性 - 涵盖了SISO和MIMO、稳定和不稳定、含积分器等多种情况 - 实际演示了增益设计应用 - 可视化清晰有效 ### 不足 #### 1. 理论局限 - **Square限制**: 主要结果仅适用于square系统，限制了应用范围 - **适定性**: 不能像extended SRG那样保证适定性 - **Normal矩阵假设**: Proposition 1对于一般MIMO系统，边界可能不精确 #### 2. 实践挑战 - **模型依赖**: 需要精确的传递函数模型 - **非最小相位检测**: 从数据检测具有挑战性 - **计算复杂度**: 对于高维MIMO系统，需要对连续的α进行计算 #### 3. 实验不足 - **缺少计算时间分析**: 未报告算法的计算效率 - **数值精度讨论**: 有限α集合的数值误差未量化 - **与GNC的实际对比**: 仅理论比较，缺少实际案例对比 #### 4. 呈现问题 - **符号密集**: 大量数学符号可能影响可读性 - **缺少算法伪代码**: 虽然方法清晰，但缺少可直接实现的算法描述 - **图形说明**: 某些图（如图2b的数值伪影）需要更多解释 ### 影响力评估 #### 对领域的贡献 1. **理论完备性**: 填补了SRG框架中处理不稳定系统的空白 2. **方法学进展**: 为非线性系统分析提供了新工具 3. **经典理论现代化**: 为Nyquist准则提供了现代几何解释 #### 实用价值 1. **控制设计**: 可直接用于含不稳定组件的反馈系统设计 2. **稳定性分析**: 提供了比传统方法更直观的稳定性裕度可视化 3. **教学价值**: 几何方法有助于理解反馈系统行为 #### 可复现性 - **理论可复现**: 所有证明在附录中完整给出 - **实现可行**: 基于标准工具，易于实现 - **示例清晰**: 提供了具体的传递函数，便于验证 - **缺少代码**: 未提供开源实现（这是一个遗憾） ### 适用场景 #### 最适合的应用 1. **复杂反馈系统设计**: 特别是包含不稳定组件的系统 2. **鲁棒性分析**: 需要可视化稳定性裕度的场合 3. **多变量系统**: MIMO系统的增益调整和稳定性分析 4. **非线性系统**: 作为非线性SRG分析的基础构建块 #### 不太适合的场景 1. **纯数据驱动**: 当无法获得精确模型时 2. **非square系统**: 当p ≠ q时方法不完整 3. **实时应用**: 如果计算复杂度是关键约束 4. **极高维系统**: 当维度很高时计算可能变得困难 ### 总体评价 这是一篇**高质量的理论控制论文**，在Scaled Relative Graph框架中做出了重要贡献。主要优势在于： - 提供了首个系统化方法来处理不稳定LTI系统的hard SRG - 理论严谨，证明完整 - 与经典Nyquist准则建立了清晰联系 - 为非线性系统分析奠定了基础 主要改进空间： - 扩展到非square系统 - 开发数据驱动版本 - 提供开源实现和更多实际案例 **推荐指数**: ⭐⭐⭐⭐☆ (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" - Extended SRG的定义 4. **[12] Chen, 2025**: "Soft and Hard Scaled Relative Graphs" - Hard SRG的首次定义 5. **[13] MacFarlane & Postlethwaite, 1977**: MIMO Nyquist准则的经典文献 6. **[14] Desoer & Wang, 1980**: 广义Nyquist准则（GNC）