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

ポート・ハミルトン系の辛シンプレクティック型モデル低次元化

基本情報

  • 論文ID: 2203.07751
  • タイトル: Symplectic model order reduction of port-Hamiltonian systems
  • 著者: Silke Glas, Mir Mamunuzzaman, Hongliang Mu, Hans Zwart
  • 分類: math.OC（最適化と制御）、math.CA（古典解析）
  • 発表時期: 2022年3月（arXiv プレプリント）
  • 論文リンク: https://arxiv.org/abs/2203.07751

要旨

本論文は、線形時不変ポート・ハミルトン（pH）系に対する新規な構造保持型モデル低次元化（MOR）手法を提案している。研究目標は、低次元化されたpH系が元の全次数モデルの物理領域において解釈可能なままであるような低次元化pH系を構築することである。具体的には、初期の高次元pH系が電子回路である場合、低次元化モデルも依然として電子回路として解釈可能であることが望ましい。質量-ばね-ダンパ（MSD）系に対しては、この特定の構造保持を保証するMOR手法が既に存在する。本論文は、辛型MOR手法に基づく新規な構造保持型MOR手法がMSD系の場合に既知の二次Arnoldi法を復元できることを示す。しかし、電子回路pHモデルなどの類似したブロック構造を有するモデルに対しては、本手法は新規なモデル低次元化技術をもたらす。

研究背景と動機

問題記述

ポート・ハミルトン系は古典的ハミルトン枠組みの一般化であり、工学問題において極めて有用である。なぜなら、古典的ハミルトン手法と開放物理系建模の強力な枠組みを結合するからである。しかし、制御設計において、制御器は通常、元の系と同じ次数を有するため、元の系の基本的な挙動を捕捉できる低次元化モデルの必要性が生じている。

既存手法の限界

  1. 構造喪失問題: 従来のMOR手法を直接適用するとpH構造が失われる
  2. 物理解釈性の欠如: 既存手法は低次元化モデルの物理領域における解釈可能性を保証できない
  3. ブロック構造の破壊: 特殊なブロック構造を有するpH系に対して、既存手法はこの構造を保持できない

研究動機

物理モデルをpH形式に再定式化する際、行列J、R、Hは通常、特殊なブロック構造を有する。本論文の目標は、低次元化プロセス中にこのブロック構造を保持し、低次元化モデルが全次数モデルと同じタイプの物理系として解釈可能であることを確保することである。

核心的貢献

  1. 辛型pH手法の導入: ブロック構造形式（2）に適用可能な辛型MORに基づくpH系の新規なMOR手法を提案し、低次元化モデルが同じ形式の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ᵖは出力である。目標は、同じ形式であるが状態空間次元がより小さい系を見つけることであり、2つの伝達関数が予定された周波数点で等しくなることを要求する。

核心的理論枠組み

シンプレクティック幾何学の基礎

定義1（シンプレクティック形式）: Vを複素数体上の偶数次元ベクトル空間とし、シンプレクティック形式Ω: 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対pH-Arnoldi: 辛型pHはほとんどの場合、より低いH₂およびH∞誤差を提供
  2. 辛型pH対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系建模、シンプレクティック幾何学、制御理論など複数の領域における重要な研究をカバーし、研究に堅実な理論的基礎を提供している。