技術要約：一般化リャプノフ方程式を用いたスイッチング系における到達可能集合および可観測集合

1. 問題設定

本論文は、スイッチング微分代数方程式（DAE）のモデル次数低減（MOR）において、一般化リャプノフ方程式（GLE）の解を用いることの理論的正当性を扱う。具体的には、以下の形式のスイッチング系を対象とする：
$$\begin{cases} E_{q(t)}\dot{x}(t) = A_{q(t)}x(t) + B_{q(t)}u(t), \\ y(t) = C_{q(t)}x(t), \end{cases}$$
ここで、$q(t)$ はスイッチング信号であり、行列 $E_{q(t)}$ は特異である可能性がある。

核心となる課題は、このようなシステムの到達可能集合および可観測集合を厳密に特徴付けることにある。近年の研究（Manucci & Unger, 2024）では、これらのシステムに対して均衡に基づくMOR手法を導出するためにGLEを用いることが提案されているが、これらのGLEの解と、スイッチングDAEの実際の到達可能・可観測部分空間との間の理論的な結びつきは、完全には確立されていなかった。本論文の目的は、GLEの解の像（image set）が、常にシステムの到達可能集合および可観測集合を包含することを証明し、それによってMORのためのグラム行列としての使用を正当化することである。

2. 手法

2.1 解空間と再定式化

著者らは、モード間の切り替え時に発生するジャンプやディラックのインパルスを扱うために、区分的に滑らかな分布（$\mathcal{D}^{\infty}_{pw}$）の空間内で作業を行う。

**準ワイエルシュトラス形式（QWF）**を用いることで、スイッチングDAEを、状態依存のジャンプとインパルスを伴う等価なスイッチング常微分方程式（ODE）のシステムへと再定式化する。再定式化されたシステム（8）は以下のように与えられる：
$$\begin{cases} \dot{x}(t) = A^{\text{diff}}_{q_k}x(t) + B^{\text{diff}}_{q_k}u(t), & t \in (t_k, t_{k+1}), \\ x(t_k^+) = \Pi_{q_k}x(t_k^-) + \Pi_{q_k}B_{q_{k-1}}U_{q_{k-1}}(t_k^-), \\ y(t) = C^{\text{diff}}_{q_k}x(t) + D_{q_k}U_{q_k}(t), \\ y[t_k] = x(t_k^-)に依存するインパルス項. \end{cases}$$
ここで、$\Pi_j$ は整合性射影（consistency projector）を表し、$U$ はシステムのべき零指数までの制御信号の微分を含む。

2.2 到達可能性と可観測性の特徴付け

著者らは、この再定式化されたシステムに対して到達可能集合および可観測集合を厳密に定義する：

• 到達可能性： 著者らは、入力依存のジャンプを考慮した「修正到達可能集合」を定義する。彼らは（補題3.3）により、到達可能集合が再帰的に計算可能であることを証明する。極めて重要なことに、定理3.4は、入力・状態依存のジャンプを持つシステムの到達可能集合が、状態依存のジャンプと拡張された入力行列を持つ簡略化されたシステムの到達可能集合に含まれることを示している。
• 可観測性： 同様に、インパルス出力項を考慮して、不可観測集合が再帰的に特徴付けられる（補題4.3）。定理4.5は、元のシステムの不可観測集合が、拡張された出力行列を持つが、出力インパルスを持たないシステムの不可観測集合と等価であることを示している。

2.3 一般化リャプノフ方程式への接続

著者らは、以下の条件を満たす、さらに簡略化されたスイッチング線形システム（35）を導入する：

1. ジャンプは純粋に状態依存である（入力依存のジャンプはない）。
2. 入力と出力は、元のシステムのインパルス効果を捉えるために拡張されている。
3. 射影行列 $\Pi_j$ は一般的に扱われる。

定理5.1は、決定的な包含関係を確立する：

• ジャンプを持つシステムの到達可能集合 $R$ は、ジャンプを持たないシステム（ただし拡張された行列を持つ）の到達可能集合 $\tilde{R}$ の部分集合である。
• 可観測集合 $O$ は、ジャンプのないシステムの可観測集合 $\tilde{O}$ の部分集合である。

拡張された行列を持つジャンプのないシステムについては、到達可能集合および可観測集合が、以下の一般化リャプノフ方程式の解 $P$ および $Q$ の像空間と一致するという既知の結果（[13]より引用）がある：
$$\begin{aligned} AP + PA^T + \sum_{j=1}^M (F_j P F_j^T + \tilde{B}_j \tilde{B}_j^T) &= 0, \\ A^T Q + QA + \sum_{j=1}^M (F_j^T Q F_j + \tilde{C}_j^T \tilde{C}_j) &= 0. \end{aligned}$$

3. 主な貢献と結果

1. 厳密な特徴付け： 本論文は、区分的に滑らかな分布の空間内における、スイッチングDAEの到達可能集合および可観測集合の形式的な定義と再帰的な特徴付けを提供する。
2. 包含定理： 主要な理論的結果は、複雑なスイッチングDAE（ジャンプとインパルスを伴う）の到達可能集合および可観測集合が、特定の拡張されたスイッチングODEシステムの到達可能集合および可観測集合の部分集合であるという証明である。
3. GLEの妥当性確認： 包含定理とスイッチングODEに関する既存の結果（[12]）を組み合わせることで、著者らは、GLEの解 $P$ および $Q$ が以下を満たすことを証明する：
   $$\text{到達可能集合} \subseteq \text{span}(P) \quad \text{および} \quad \text{可観測集合} \subseteq \text{span}(Q).$$
   これにより、GLEの解の像が真のシステム部分空間を「包囲（enclose）」していることが確認される。

4. 意義と主張

本論文は、スイッチングDAEにおけるモデル次数低減のための、GLEに基づく均衡（balancing）の使用を正当化するために必要な理論的基礎を提供すると主張している。

• MORの正当化： 著者らは、GLEの解が到達可能集合および可観測集合を包囲しているため、これらは均衡に基づくMORのためのグラム行列として適していると述べている。これにより、低次モデルが元のスイッチングDAEシステムの入出力挙動を保持することが保証される。
• 主張の範囲： 著者らの主張は控えめである。彼らは、この特定の著作において新しい数値アルゴリズムや実験的な応用を提案しているのではない。代わりに、彼らの先行研究（[12]）で提案された手法が「なぜ機能するのか」を明らかにしている。彼らは、この再定式化が入出力写像を変更しないことを明記しており、これがMORの主要な目的である。
• 限界： 論文では、再定式化が隣接するモードに依存する拡張行列を含むため、$M(M-1)^2$ 個のモードを持つシステムになる可能性があることを認めている。しかし、GLEの解の像の集合は、この特定の定式化に関わらず一貫していると述べている。

要約すると、本論文は、分布とインパルスを含むスイッチングDAEの抽象的な理論と、GLEを用いたモデル低減の実践的な適用との間の溝を埋め、後者が前者の本質的なダイナミクスを正しく捉えていることを証明している。