Characterization of stability radii for robustly asymptotically stable dissipative Hamiltonian differential-algebraic systems

技術的概要：ロバスト漸近安定な散逸ハミルトニアン微分代数方程式系の安定半径の特性評価

問題定義
本論文は、$E \dot{x} = (J - R)Qx$ の形式を持つ線形時不変の散逸ハミルトニアン微分代数方程式（dHDAE）系のロバスト漸近安定性に取り組む。このような系の有限固有値が閉左半平面に存在し、純虚数固有値（ゼロを除く）が半単純であることはよく知られているが、これらの系が必ずしもロバストに漸近安定であるわけではない。係数行列の微小な摂動により、正則性が失われたり、系の指数が増加したり、あるいは固有値が虚軸上に移動して不安定化したりする可能性がある。

中心的な問題は、与えられたロバスト漸近安定な dHDAE 系から、そのような系の集合の境界までの「距離」を定量化することである。具体的には、著者らは以下の事態を引き起こすために必要な最小の構造化摂動を特性評価することを目的としている：

1. 特異性（正則性の喪失）。
2. 微分代数指数の 1 超への増加。
3. 純虚数固有値の出現（漸近安定性の喪失）。

本解析では、2 種類の構造化摂動を考慮する：

• $S_d$（散逸的）： $E$ と $R$ の符号決定性（すなわち $\Delta E \le 0, \Delta R \le 0$）および $J$ の斜エルミート構造を保持する摂動。
• $S_i$（不定）： エルミート/斜エルミート構造を保持するが、摂動自体の符号決定性の制約は必ずしも保持せず、結果として得られる行列が半正定値であることのみを要求する摂動。

手法
著者らは、構造行列論、スペクトル解析、および最適化手法の組み合わせを採用する：

1. 構造的特性評価： 本論文は、dHDAE ペンシルの正則性と指数を解析するために、クラネッカー標準形およびユニタリ変換を通じて導出された特定の「階段形」（補題 8）を利用する。ロバスト漸近安定性は、正則性、指数 $\le 1$、および負のスペクトル半幅という 3 つの条件によって特性評価される（定理 11）。
2. 標準形への還元： 解析において、$Q$ が可逆であると仮定し、系を $E \dot{x} = (J - R)x$ の形に変換する。ここで $E=E^* \ge 0$、$R=R^* \ge 0$、$J=-J^*$ である。
3. 写像問題： 不安定性までの距離を計算するために、著者らは特定の代数制約（例えば、行列を特異にすること、または虚軸上に固有値を生成すること）を満たす最小ノルム摂動 $(\Delta E, \Delta J, \Delta R)$ を求める問題として定式化する。エルミートおよび斜エルミート写像問題に対する最小ノルム解を活用する（定理 15 および 16）。
4. 最適化定式化： 距離は、一般化されたレイリー商を含む最適化問題として表現される。$E$ が固定されている場合（部分的摂動）については、厳密な特性評価が導出される。完全摂動（$E$ が変化する）については、著者らは境界値を導出する。
   • 特異性および高指数への距離は、ランク条件に制約された摂動のノルム二乗和の最小化を含む。
   • 虚数固有値への距離は、周波数パラメータ $\omega$ およびベクトル $x$ に関する最小化を含む。
5. 数値実装： 導出された最適化問題は、2 段階の戦略を用いて解かれる。内側ループでは、レベルシフト付き自己無撞着場（SCF）反復法を通じて非線形固有値問題（NEPv）を解き、外側ループでは標準的な最適化ルーチン（例：fminsearch、fmincon）を用いて周波数パラメータ $\omega$ について最小化する。

主要な貢献と結果

• 部分的摂動に対する厳密な特性評価： 行列 $E$ が固定され、$J$ と $R$ のみが摂動される場合、本論文は特異性、高指数、および虚数固有値への構造化距離の厳密な数式を提供する（定理 24–26）。
  • $R > 0$ である $S_i$ 摂動の場合、特異性への距離は $N$ が $E$ の核空間を張るものとして、$\sqrt{\lambda_{\min}(N^*(R^2 - J^2)N)}$ に等しい。
  • $S_i$ に対する虚数固有値への距離は、$\inf_{\omega} \sigma_{\min}([R, i\omega E - J])$ で与えられる。
• 完全摂動に対する境界値： すべての行列（$E, J, R$）で摂動が許容される場合、著者らは明示的な下限および上限を導出する：
  • 特異性への距離： $S_d$ については厳密な数式が提供され（定理 20）、$S_i$ については境界値が与えられる（定理 21）。
  • 高指数への距離： 2 段階の摂動アプローチ（まず $E$ を摂動して核空間の次元を増加させ、次に $J$ と $R$ を摂動する）を採用することで上限が導出される（定理 22）。著者らは、2 次摂動項の存在により、これらは現時点では厳密な特性評価ではなく上限であることを指摘する。
  • 虚数固有値への距離： $S_d$ および $S_i$ 両方に対して下限が導出される（定理 18）。これらには、特定の行列ペンシルの一般化レイリー商および固有値に関する複雑な最適化が関与する。
• アルゴリズム的枠組み： 標準的な線形代数分解（SVD、スペクトル分解）を用いてロバスト漸近安定性をチェックし、関連する距離の境界値を計算するためのステップバイステップのアルゴリズム（セクション 2.7）が提示される。

意義と主張
本論文は、構造保存摂動下における dHDAE 系の安定半径の包括的な特性評価を提供する最初の試みであると主張する。一般的な非構造化 DAE に対する安定半径の計算は困難であることが知られているが、著者らは、散逸ハミルトニアン系の特定の幾何学的および代数的性質により、明示的な数式および計算可能な境界値の導出が可能であることを実証する。

その意義は、構造力学、流体力学、多孔質媒質流れなどの物理モデルが、モデリングおよび離散化誤差に対してどの程度の耐性を持つかを定量化できる点にある。最も近い特異系、高指数系、または不安定系までの距離を決定することにより、エンジニアおよび科学者は dHDAE モデルの信頼性を評価できる。

著者らは、完全摂動（特に高指数への距離）に対する境界値の緊密性については謙虚であり、2 次摂動項の複雑さにより現在の結果は上限であることを認めている。これらの場合に対する厳密な数式の決定およびより効率的な数値最適化手法の開発を、今後の研究に必要な方向性として特定している。