2025-11-12T23:07:10.120463

Symplectic model order reduction of port-Hamiltonian systems

Glas, Mamunuzzaman, Mu et al.
This work proposes a novel structure-preserving model order reduction (MOR) method for linear, time-invariant port-Hamiltonian (pH) systems. Our goal is to construct a reduced order pH system, which can still be interpreted in the physical domain of the full order model. By this we mean, that if an electrical circuit is the initial high-dimensional pH system, we want the reduced order model to be still interpretable as an electronic circuit. In the case of the well-known mass spring damper (MSD) system, there are MOR methods available, which already guarantee the preservation of this particular structure. Moreover, we show that our new structure-preserving MOR method, which is based on symplectic MOR methods, will recover the known second-order Arnoldi method in the case of MSD systems. However, for the example of an electrical circuit pH model (and more models of similar block structure), our method yields a novel model reduction method. We present numerical results on the aforementioned electronic circuit model, highlighting the advantages of the proposed method.
academic

Symplectic model order reduction of port-Hamiltonian systems

基本信息

  • 论文ID: 2203.07751
  • 标题: Symplectic model order reduction of port-Hamiltonian systems
  • 作者: Silke Glas, Mir Mamunuzzaman, Hongliang Mu, Hans Zwart
  • 分类: math.OC (Optimization and Control), math.CA (Classical Analysis)
  • 发表时间: 2022年3月 (arXiv预印本)
  • 论文链接: https://arxiv.org/abs/2203.07751

摘要

本文提出了一种新颖的结构保持模型降阶(MOR)方法，用于线性时不变端口哈密顿(pH)系统。研究目标是构建一个降阶pH系统，使其仍能在原全阶模型的物理域中进行解释。具体而言，如果初始高维pH系统是一个电路，我们希望降阶模型仍可解释为电子电路。对于质量-弹簧-阻尼(MSD)系统，已有MOR方法能保证这种特定结构的保持。本文展示了基于辛MOR方法的新结构保持MOR方法在MSD系统情况下能恢复已知的二阶Arnoldi方法。然而，对于电路pH模型等具有相似块结构的模型，该方法产生了新颖的模型降阶技术。

研究背景与动机

问题描述

端口哈密顿系统是经典哈密顿框架的推广，在工程问题中极其有用，因为它将经典哈密顿方法与开放物理系统建模的强大框架相结合。然而，在控制设计中，控制器通常与原系统具有相同的阶数，这促使了对能捕获原系统基本行为的降阶模型的需求。

现有方法的局限性

  1. 结构丢失问题: 直接应用传统MOR方法会丢失pH结构
  2. 物理解释性缺失: 现有方法无法保证降阶模型在物理域中的可解释性
  3. 块结构破坏: 对于具有特殊块结构的pH系统，现有方法无法保持这种结构

研究动机

当将物理模型重新表述为pH形式时，矩阵J、R和H通常具有特殊的块结构。本文的目标是在降阶过程中保持这种块结构，确保降阶模型仍可解释为与全阶模型相同类型的物理系统。

核心贡献

  1. 引入辛-pH方法: 提出了基于辛MOR的新颖pH系统MOR方法，适用于块结构形式(2)，使得降阶模型保持相同形式的pH系统
  2. 理论连接: 证明在MSD系统的特殊情况下，所提方法恢复了二阶Arnoldi方法(SO-Arnoldi)
  3. 数值验证: 在电子电路(LRCR电路)上展示数值结果，与文献中的其他已知方法进行比较
  4. 稳定性优势: 发现SO-Arnoldi方法在某些情况下产生不稳定的降阶模型，而本文方法保持稳定性

方法详解

任务定义

给定线性时不变pH系统：

ẋ(t) = (J - R)Hx(t) + Bu(t)
y(t) = B*Hx(t)

其中x(t) ∈ C²ⁿ是状态，u(t) ∈ Cᵖ是输入，y(t) ∈ Cᵖ是输出。目标是找到相同形式但状态空间维数更小的系统，使得两个传递函数在预定频率点相等。

核心理论框架

辛几何基础

定义1 (辛形式): 设V是C上的偶数维向量空间，辛形式Ω: V × V → C满足：

  1. 反对称性：Ω(v₁, v₂) = -Ω(v₂, v₁)
  2. 非退化性：对所有v₁ ∈ V有Ω(v₁, v₂) = 0当且仅当v₂ = 0

定义2 (辛映射): 矩阵Q称为(J₂ₙ, J₂ₖ)-辛的，如果：

Q*J₂ₙQ = J₂ₖ

无耗散情况的辛MOR

定理4: 对于无耗散pH系统(R = 0)，如果Q ∈ C²ⁿˣ²ᵏ是(J₂ₙ⁻¹, J₂ₖ⁻¹)-辛矩阵，且满足：

  • ran(B) ⊆ ran(Q)
  • (sₘI₂ₙ - J₂ₙH)⁻¹Bμₘ ∈ ran(Q)

则降阶pH系统在给定频率处与原系统相等。

有耗散情况的辛MOR

定理5: 对于有耗散的pH系统，构造扩展矩阵：

QE := [Q₁  0  ]
      [0   Q₂]

其中Q₂是(J₂ₙ, J₂ₖ)-辛的，Q₁*是Q₂的左逆。在适当条件下，降阶系统满足插值条件。

块结构保持算法

算法1 (辛-pH MOR):

  1. 根据插值数据构造V₁
  2. 计算V₂ := HV₁
  3. 按块结构划分V₁和V₂
  4. 通过公式(27)构造Ψ₁, Ψ₂, Φ₁, Φ₂
  5. 构建Q₁和Q₂
  6. 应用Petrov-Galerkin投影

构造公式：

Ψ₂ = orth(V₂₂)
Ψ₁ = V₁₂(Ψ₂*V₁₂)⁻¹
Φ₂ = V₂₁(Ψ₂*Ĵₙ*V₂₁)⁻¹Ĵₖ*
Φ₁ = V₁₁(Φ₂*V₁₁)⁻¹

技术创新点

  1. 辛结构保持: 通过辛映射保持哈密顿系统的基本几何结构
  2. 块对角结构维护: 特殊的Q₁和Q₂构造确保物理意义的块结构得以保持
  3. 插值条件: 保证降阶系统在预定频率点精确匹配原系统
  4. 统一框架: 将MSD系统的SO-Arnoldi方法作为特例包含在内

实验设置

数值算例：LRCR电路

考虑图1所示的LRCR电路，包含n=50个单元，参数设置：

  • 电容：Cᵢ = 1 pF
  • 电感：Lᵢ = 1 nH
  • 电阻：Rᵢ,₁ = 1 mΩ, Rᵢ,₂ = 1 kΩ

评价指标

  • H₂误差: 
    H₂error = √(1/2π ∫ Trace((G(iω) - Gred(iω))*(G(iω) - Gred(iω)))dω)
  • H∞误差: 
    H∞error = max‖G(iω) - Gred(iω)‖₂

对比方法

  1. pH-Arnoldi: 基于定理2的结构保持方法
  2. SO-Arnoldi: 二阶Arnoldi方法
  3. 辛-pH: 本文提出的方法

插值点选择

在区间10⁻², 10²内均匀选择对数尺度上的插值点，例如6个插值点： 10⁻²i, -10⁻²i, 10⁰i, -10⁰i, 10²i, -10²i

实验结果

主要结果

从图3可以观察到：

  1. 辛-pH vs pH-Arnoldi: 辛-pH在大多数情况下提供更低的H₂和H∞误差
  2. 辛-pH vs SO-Arnoldi:
    • 降阶维数<24时性能相近
    • 降阶维数=32时SO-Arnoldi性能急剧下降
    • 降阶维数28和36时SO-Arnoldi产生不稳定系统

稳定性分析

图4显示了SO-Arnoldi方法在r=28和r=36时产生的不稳定极点，而辛-pH方法在所有测试维数下保持稳定性。

理论验证：MSD系统等价性

定理7证明了对于MSD系统，辛-pH方法与SO-Arnoldi方法通过酉变换等价： 存在酉矩阵Θ使得zred(t) = Θẑred(t)

相关工作

pH系统MOR发展

  1. 早期工作: 基于正实传递函数的被动系统方法
  2. 矩量匹配: Polyuga和van der Schaft的工作
  3. Dirac结构近似: 结构保持的不同途径
  4. 谱分解方法: 基于Popov函数的方法
  5. 优化方法: 基于参数优化和黎曼优化的方法

辛MOR发展

  1. 自治哈密顿系统: Peng和Mohseni的开创性工作
  2. 非线性方法: 基于流形的辛降阶
  3. 数据驱动方法: 算子推断和系统辨识

结论与讨论

主要结论

  1. 成功开发了保持pH结构和块结构的MOR方法
  2. 理论上统一了MSD系统的不同MOR方法
  3. 数值实验验证了方法的有效性和稳定性优势

局限性

  1. 插值方法: 目前仅匹配传递函数值，未考虑导数匹配
  2. 块结构限制: 方法专门针对特定块结构设计
  3. 计算复杂度: 构造辛矩阵的计算成本需要进一步分析

未来方向

  1. 导数匹配: 扩展到匹配传递函数的导数以提高拟合质量
  2. 自适应选择: 开发插值点的自适应选择策略
  3. 非线性扩展: 扩展到非线性pH系统
  4. 大规模应用: 优化算法以处理更大规模的系统

深度评价

优点

  1. 理论严谨: 基于辛几何的坚实数学基础
  2. 结构保持: 同时保持pH结构和物理意义的块结构
  3. 统一性: 将现有方法作为特例统一在一个框架下
  4. 实用性: 方法具有明确的算法实现和数值验证

不足

  1. 适用范围: 方法专门针对特定块结构，通用性有限
  2. 计算效率: 辛矩阵构造的计算复杂度分析不够充分
  3. 大规模验证: 缺乏更大规模系统的数值验证

影响力

  1. 理论贡献: 为pH系统MOR提供了新的理论视角
  2. 实用价值: 在电路分析等工程应用中具有直接价值
  3. 方法学意义: 展示了几何方法在MOR中的重要作用

适用场景

  1. 电路分析: 大规模电子电路的降阶建模
  2. 机械系统: 质量-弹簧-阻尼系统的控制设计
  3. 多物理场系统: 具有能量保持特性的复杂物理系统

参考文献

论文包含31篇参考文献，涵盖了MOR理论、pH系统建模、辛几何和控制理论等多个领域的重要工作，为研究提供了坚实的理论基础。