On the Lyapunov equation with the state matrix in companion form

技術要約：コンパニオン形式の状態行列を持つリアプノフ方程式について

問題提起
本論文は、以下の連続時間リアプノフ方程式の唯一の解 $X$ の成分ごとの非負性について研究している：
$$XA + A^\top X = -Q$$
ここで、$A \in \mathbb{R}^{n \times n}$ は多項式 $a(s)$ に付随するフルリッツ（Hurwitz）コンパニオン行列であり、$Q = Q^\top \succeq 0$ は正定値行列である。

標準的なリアプノフ理論は、$X$ が対称かつ正定値（ローエナー順序において）であることを保証するが、この性質は $X$ が成分ごとに非負であること（すなわち、すべての $i, j$ に対して $X_{ij} \geq 0$ であること）を意味するわけではない。本論文は、ルイージ・フォルトゥナ教授によって提唱された予想 1.1、すなわち「$A$ がフルリッツ・コンパニオン行列であり、$Q \succeq 0$ ならば、解 $X$ は成分ごとに非負である」という命題に取り組んでいる。この性質は正定性とは異なるものであり、正のシステム（positive systems）の理論において特に重要である。なぜなら、成分ごとの非負性は、質量、濃度、あるいは個体数といった物理量の保存を保証するためである。

手法
証明戦略は、問題を初等対称多項式とコーシー型行列を含む特定の代数的構造へと還元することに基づいている。手法は以下のステップで進行する：

1. 問題の簡略化:

   • ランク1への還元: リアプノフ方程式の線形性と、任意の正定値行列 $Q$ がランク1行列の和（$Q = \sum q_i q_i^\top$）に分解できることを利用し、著者らは一般性を失うことなく $Q = qq^\top$ と仮定する。
   • 対角化可能性: 著者らは、$X$ は $A$ の連続関数であると論じている。したがって、結果を証明するには、対角化可能なコンパニオン行列（固有値が相異なる場合）について証明すれば十分である。非対角化可能なケースは、係数が極めて近い対角化可能なケースによって近似できる。

2. スペクトル分解:

   • $A$ が相異なる実固有値 $-\lambda_i$ （$0 < \lambda_1 < \dots < \lambda_n$ と順序付けられたもの）を持つ対角化可能な行列であると仮定し、著者らは $A^\top$ の固有ベクトルを利用する。これらの固有ベクトル $v_j$ は、固有値の初等対称多項式 $\sigma_m^{(j)}$ を用いて構成される。
   • 強制ベクトル $q$ は、この固有基底において $q = \sum a_i v_i$ と展開される。

3. 積分表示と二次形式:

   • リアプノフ解の積分公式 $X = \int_0^\infty e^{A^\top t} Q e^{At} dt$ を用いて、著者らは成分 $X_{kl}$ の明示的な表現を導出する。
   • 成分 $X_{kl}$ は、以下の要素を持つ行列 $M$ による二次形式 $a^\top M a$ であることが示される：
     $$M_{ij} = \frac{\sigma_{n-k}^{(i)}\sigma_{n-l}^{(j)} + \sigma_{n-k}^{(j)}\sigma_{n-l}^{(i)}}{\lambda_i + \lambda_j}$$
   • この行列 $M$ は、$M = \Delta C \Delta$（$\Delta$ は対角行列、$C$ は $C_{ij} = \frac{r_i + r_j}{\lambda_i + \lambda_j}$ を要素とするコーシー型行列）と因数分解される。

4. 正値性の検証:

   • $X_{kl}$ の非負性は、$M$ が正定値であることから従う。$\Delta$ が非特異であるため、これは $C \succeq 0$ であることを証明することと同値である。
   • 著者らは、このようなコーシー型行列が正定値であることを確認した既知の結果 [8] を引用している。

主要な貢献および結果

• 定理 1 (主結果): 本論文は、コンパニオン行列 $A$ が実固有値のみを持つという特定の条件下において、予想 1.1 を証明している。この仮定の下では、一意の解 $X$ は成分ごとに非負である。
• 全正値性解析: 著者らはより強い性質である全正値性（すべての小行列式が非負であること）についても調査している。
  • 彼らは、$A$ の固有値が実数であっても、$X$ が一般に全正値ではないことを示す反例を提示している。
  • 定理 2: 著者らは、全正値性の十分条件を確立している。すなわち、強制ベクトル $q$ が $A^\top$ の固有基底において非負の展開を持つ場合（具体的には $q = \sum \alpha_j v_j$ かつ $\alpha_j \geq 0$）、$X$ は全正である。これは、$X$ を全正行列の積（ヴァンデルモンド行列とコーシー行列を含む）へと因数分解することによって証明される。

意義および主張
本論文は、コンパニオン行列に関するリアプノフ解の構造的性質に関する長年の予想を解決したと主張している。

• 動機: $X$ の成分ごとの非負性は、正のシステムの最適制御における応用から動機づけられている。具体的には、非負の入力行列を持つ単一入力システムにおいて、関連する代数リカッチ方程式（リアプノフ方程式として書き換えられる）の成分ごとに非負な解が存在することは、最適レギュレータが正の電圧供給のみを用いたハードウェアで実装可能であることを意味する。
• 範囲と限界: 著者らは、自身の証明の一般性について謙虚な姿勢を示している。彼らは、実固有値という仮定が制限であることを明示的に述べている。本論文は、複素固有値を持つフルリッツ・コンパニオン行列に対しても予想が成立するかどうかという問いを未解決のまま残している。
• 理論的価値: 本研究は、リアプノフ理論をコーシー型行列および初等対称多項式の理論と結びつけ、$X$ の成分ごとの非負性が、Kwong（またはアンチ・ローエナー）行列に双対な特定のクラスの行列の正定性に由来することを明らかにしている。

要約すると、本論文は、実スペクトルを持つフルリッツ・コンパニオン行列に対するリアプノフ解が成分ごとに非負であることを厳密に証明し、一方で全正値性には強制ベクトルに対する追加の制約が必要であることを明確にしている。