# 行列の固有値評価と不等式：ゲルシュゴリンからペロン・フロベニウス、カントロヴィッチまで > ゲルシュゴリンの円板定理、クーラン・フィッシャーとワイル・交互定理、シュール補行列とウッドベリー公式、アダマールとコーシー・ビネ、ペロン・フロベニウス、リャプノフ方程式、カントロヴィッチの不等式を証明する。 > https://rikai.mugen-giken.com/mathematics/linear-algebra/matrix-inequalities ## 0. この記事の要点 行列の固有値を実際に求めるのは大変ですが、「固有値は正か」「どのくらいの範囲に入るか」だけ分かれば十分なことがよくあります。この記事では、行列の成分や部分行列、行列式といった手に入りやすい量から、固有値を挟み込む道具を順に紹介します。読み終えると、固有値を計算しなくても、小さな行列なら正定値かどうかや固有値の範囲を見積もれるようになります。 - **ゲルシュゴリンの円板定理**: $A$ の固有値は、中心 $a_{ii}$、半径 $R_i=\sum_{j\ne i}|a_{ij}|$ の円板の和集合に入ります。固有値を計算せずに、対角優位な行列の正則性や[正定値](/mathematics/linear-algebra/spectral-theorem#def-positive-definite)性が判定できます。 - エルミート行列の固有値 $\lambda_1\ge\cdots\ge\lambda_n$ は、**クーラン・フィッシャーの最大最小原理**でレイリー商の言葉で書けます。そこから**ワイルの不等式**（摂動で固有値は $\|E\|$ 以上動かない）と**コーシーの交互定理**（主部分行列の固有値は元の固有値の間に挟まる）が従います。 - ブロック行列の**シュール補行列** $A-BD^{-1}C$ は、行列式・逆行列・正定値性・慣性を支配します。**ウッドベリーの公式**と**シルヴェスターの行列式恒等式** $\det(I+XY)=\det(I+YX)$ はその系で、低ランク更新の計算を小さな行列に帰着します。 - **アダマールの不等式** $|\det A|\le\prod_j\|\boldsymbol{a}_j\|$ と**コーシー・ビネの公式**は、行列式の大きさと小行列式の和の評価です。 - 非負行列の**ペロン・フロベニウスの定理**は、既約なら最大固有値が単純で正の固有ベクトルをもつことを述べ、**コリャッツ・ヴィーラントの公式**で上下から評価できます。**シルヴェスター方程式**の可解条件と**リャプノフ方程式**による安定性判定、**カントロヴィッチの不等式**まで、一貫して「固有値を求めずに評価する」道具です。 ## 1. 動機 — 固有値を計算せずに、固有値について知りたい 実用では、行列の固有値そのものより「固有値が正か」「最大固有値の上限はいくらか」「摂動で固有値がどれだけ動くか」が重要になる場面が多くあります。たとえば、数値計算の結果が数学的に正しいことを保証する**認証計算**では、丸め誤差を含んだ近似固有値に対して、真の固有値が存在する範囲を厳密に示す必要があります。多体量子系のフェルミオン行列式の評価では、行列式の絶対値を主要部分行列の大きさで抑えるアダマールの不等式が使われます。連続時間の線形系 $\dot{\boldsymbol{x}}=A\boldsymbol{x}$ の安定性は、$A$ のすべての固有値の実部が負であること、すなわち特性多項式の根の位置の問題です。 この記事の定理は、固有値を直接求める代わりに、行列の成分、部分行列、行列式といった量から固有値を挟み込む道具です。扱う行列は、断らない限り複素正方行列で、エルミート行列の固有値は $\lambda_1(A)\ge\cdots\ge\lambda_n(A)$ と大きい順に並べます。ベクトルの内積は [スペクトル定理](/mathematics/linear-algebra/spectral-theorem) と同じ $\langle\boldsymbol{x},\boldsymbol{y}\rangle=\boldsymbol{y}^{*}\boldsymbol{x}$、ノルムは $\|\boldsymbol{x}\|$、行列の作用素ノルムは $\|A\|=\max_{\|\boldsymbol{x}\|=1}\|A\boldsymbol{x}\|$（）です。 ## 2. ゲルシュゴリンの円板定理 $A=(a_{ij})\in M_n(\mathbb{C})$ に対し、$R_i=\sum_{j\ne i}|a_{ij}|$（第 $i$ 行の非対角成分の絶対値の和）とおき、 $$ D_i=\{z\in\mathbb{C}\mid|z-a_{ii}|\le R_i\}\qquad(i=1,\ldots,n) $$ を**ゲルシュゴリン円板**という。 $A\in M_n(\mathbb{C})$ のすべての固有値は $D_1\cup\cdots\cup D_n$ に含まれる。 $\lambda$ を固有値、$\boldsymbol{x}=(x_1,\ldots,x_n)^{\mathsf{T}}\ne0$ を固有ベクトルとし、絶対値が最大の成分を $x_i$ とします（$|x_i|>0$）。第 $i$ 成分の式 $\sum_ja_{ij}x_j=\lambda x_i$ を $(\lambda-a_{ii})x_i=\sum_{j\ne i}a_{ij}x_j$ と書き直すと、 $$ |\lambda-a_{ii}|\,|x_i|\le\sum_{j\ne i}|a_{ij}||x_j|\le\sum_{j\ne i}|a_{ij}||x_i|=R_i|x_i| $$ です。$|x_i|>0$ で割って $|\lambda-a_{ii}|\le R_i$、すなわち $\lambda\in D_i$ です。
a11 = 1 a22 = 3 + 2i a33 = 2 - 2i D1 D2 D3
$A$ と $A^{\mathsf{T}}$ は同じ固有値をもつので、列について和をとった円板でも同じことが言えます。円板定理の最も使いやすい系は次の 2 つです。 1. （レヴィ・ドゥスプランクの定理）すべての $i$ について $|a_{ii}|>R_i$（**狭義対角優位**）ならば、$A$ は正則である。 2. $A$ がエルミート行列で、すべての $i$ について $a_{ii}>R_i$ ならば $A$ は正定値である。 3. すべての固有値の絶対値は $\max_i\sum_j|a_{ij}|$ 以下である。 1. $0\in D_i$ ならば $|a_{ii}|\le R_i$ となって仮定に反するので、 により $0$ は固有値でなく、$A$ は正則です（）。2. エルミート行列の固有値は実数で、$D_i\cap\mathbb{R}\subseteq[a_{ii}-R_i,a_{ii}+R_i]$ の左端 $a_{ii}-R_i$ が正ですから、固有値はすべて正です（）。3. $z\in D_i$ なら $|z|\le|a_{ii}|+R_i=\sum_j|a_{ij}|$ です。 $A=\begin{pmatrix}10&1&0\\1&8&2\\0&2&5\end{pmatrix}$ は対称で、$R_1=1,R_2=3,R_3=2$ で、円板は $D_1=[9,11]$、$D_2=[5,11]$、$D_3=[3,7]$ です（実対称行列なので固有値は実数で、実軸上の区間として考えれば足ります）。したがって固有値はすべて $[3,11]$ に入り、$A$ は正定値、最大固有値は $11$ 以下で、条件数 $\lambda_1/\lambda_3\le 11/3$ です。固有値を求めずに、これだけの情報が得られます。 円板の和集合が互いに交わらない $k$ 個の円板からなる連結成分をもつとき、その成分にはちょうど $k$ 個の固有値が含まれることも知られています。これは固有値が成分の連続関数であることを使う議論（対角成分だけの行列から非対角成分を連続的に入れていく）で示され、この記事では扱いません。 ## 3. 固有値の変分的特徴づけ [エルミート行列](/mathematics/linear-algebra/spectral-theorem#def-hermitian-unitary) $A$ の固有値は、正規直交固有ベクトル $\boldsymbol{u}_k$ で $A\boldsymbol{u}_k=\lambda_k\boldsymbol{u}_k$ となる列でした。レイリー商の最大・最小（）は $\lambda_1$ と $\lambda_n$ を与えます。途中の $\lambda_k$ を、固有ベクトルを使わずに部分空間の言葉で書くのが次の定理です。 $A$ を $n$ 次エルミート行列、$\lambda_1\ge\cdots\ge\lambda_n$ をその固有値とする。$1\le k\le n$ に対し $$ \lambda_k=\max_{\dim S=k}\ \min_{\boldsymbol{x}\in S,\ \|\boldsymbol{x}\|=1}\langle A\boldsymbol{x},\boldsymbol{x}\rangle=\min_{\dim T=n-k+1}\ \max_{\boldsymbol{x}\in T,\ \|\boldsymbol{x}\|=1}\langle A\boldsymbol{x},\boldsymbol{x}\rangle . $$ ここで $S,T$ は $\mathbb{C}^n$ の部分空間を動く。 証明の方針は次のとおりです。固有ベクトルで展開すると $\langle A\boldsymbol{x},\boldsymbol{x}\rangle=\sum\lambda_i|c_i|^2$ と書けるので、$\lambda_k$ を与える部分空間の例はすぐ見つかります。残りは、ほかのどんな $k$ 次元部分空間 $S$ をとっても最小値が $\lambda_k$ を超えないことです。これは、$S$ と「$\lambda_k$ 以下の固有値に対応する固有ベクトルが張る空間」の次元の和が $n$ を超えるため、共通部分に少なくとも 1 本の単位ベクトルが入ることから分かります。 正規直交固有ベクトル $\boldsymbol{u}_1,\ldots,\boldsymbol{u}_n$ をとります（）。$\boldsymbol{x}=\sum c_i\boldsymbol{u}_i$ なら $\langle A\boldsymbol{x},\boldsymbol{x}\rangle=\sum\lambda_i|c_i|^2$ です。単位球面は有界閉集合でコンパクトであり（）、連続関数は最大値・最小値をもつので（）、最大・最小が存在します。 **最大最小。** $S_0=\operatorname{span}(\boldsymbol{u}_1,\ldots,\boldsymbol{u}_k)$ では、単位ベクトル $\boldsymbol{x}=\sum_{i\le k}c_i\boldsymbol{u}_i$ に対し $\langle A\boldsymbol{x},\boldsymbol{x}\rangle=\sum_{i\le k}\lambda_i|c_i|^2\ge\lambda_k$ で、$\boldsymbol{x}=\boldsymbol{u}_k$ で等号なので、$\min=\lambda_k$ です。したがって右辺の最大は $\lambda_k$ 以上です。一方、$\dim S=k$ の任意の $S$ と、$N=\operatorname{span}(\boldsymbol{u}_k,\ldots,\boldsymbol{u}_n)$（次元 $n-k+1$）に対し、次元公式（ の 1）から $\dim(S\cap N)\ge k+(n-k+1)-n=1$ です。単位ベクトル $\boldsymbol{x}\in S\cap N$ をとると $\langle A\boldsymbol{x},\boldsymbol{x}\rangle=\sum_{i\ge k}\lambda_i|c_i|^2\le\lambda_k$ ですから、$S$ 上の最小値は $\lambda_k$ 以下です。ゆえに最大最小は $\lambda_k$ です。 **最小最大。** $T_0=\operatorname{span}(\boldsymbol{u}_k,\ldots,\boldsymbol{u}_n)$ では最大値が $\lambda_k$ です。次元 $n-k+1$ の任意の $T$ は $\operatorname{span}(\boldsymbol{u}_1,\ldots,\boldsymbol{u}_k)$ と自明でない共通部分をもち、その単位ベクトル $\boldsymbol{x}$ で $\langle A\boldsymbol{x},\boldsymbol{x}\rangle\ge\lambda_k$ なので、$T$ 上の最大値は $\lambda_k$ 以上です。 これを使うと、固有値の摂動に関する基本的な不等式が短く示せます。 $A,B$ を $n$ 次エルミート行列とする。 1. $i+j-1\le n$ ならば $\lambda_{i+j-1}(A+B)\le\lambda_i(A)+\lambda_j(B)$。 2. $i+j-n\ge1$ ならば $\lambda_{i+j-n}(A+B)\ge\lambda_i(A)+\lambda_j(B)$。 3. 任意の $k$ について $|\lambda_k(A+B)-\lambda_k(A)|\le\|B\|$。 1. $A$ の正規直交固有ベクトルを $\boldsymbol{u}_i$、$B$ のものを $\boldsymbol{w}_j$ とし、$X=\operatorname{span}(\boldsymbol{u}_i,\ldots,\boldsymbol{u}_n)$（次元 $n-i+1$）、$Y=\operatorname{span}(\boldsymbol{w}_j,\ldots,\boldsymbol{w}_n)$（次元 $n-j+1$）とします。$\boldsymbol{x}\in X$ で $\langle A\boldsymbol{x},\boldsymbol{x}\rangle\le\lambda_i(A)\|\boldsymbol{x}\|^2$、$\boldsymbol{x}\in Y$ で $\langle B\boldsymbol{x},\boldsymbol{x}\rangle\le\lambda_j(B)\|\boldsymbol{x}\|^2$ です。次元公式から $\dim(X\cap Y)\ge(n-i+1)+(n-j+1)-n=n-(i+j-1)+1$ ですから、$X\cap Y$ の中に次元がちょうど $n-(i+j-1)+1$ の部分空間 $T$ をとれて、$T$ 上で $\langle(A+B)\boldsymbol{x},\boldsymbol{x}\rangle\le(\lambda_i(A)+\lambda_j(B))\|\boldsymbol{x}\|^2$ です。 の最小最大（$k=i+j-1$）から $\lambda_{i+j-1}(A+B)\le\max_{\boldsymbol{x}\in T,\|\boldsymbol{x}\|=1}\langle(A+B)\boldsymbol{x},\boldsymbol{x}\rangle\le\lambda_i(A)+\lambda_j(B)$ を得ます。 2. 1 を $-A$、$-B$ に使います。$\lambda_p(-C)=-\lambda_{n+1-p}(C)$ なので、添字 $i'=n+1-i$、$j'=n+1-j$ に対する 1 の不等式 $\lambda_{i'+j'-1}(-A-B)\le\lambda_{i'}(-A)+\lambda_{j'}(-B)$ は、$-\lambda_{i+j-n}(A+B)\le-\lambda_i(A)-\lambda_j(B)$ に等しく（$n+1-(i'+j'-1)=i+j-n$）、符号を戻して 2 を得ます。 3. 1 で $j=1$、2 で $j=n$ とすると $\lambda_k(A)+\lambda_n(B)\le\lambda_k(A+B)\le\lambda_k(A)+\lambda_1(B)$ です。エルミート行列 $B$ では $\|B\|^2=\max_{\|\boldsymbol{x}\|=1}\sum\lambda_i(B)^2|c_i|^2=\max_i\lambda_i(B)^2$ なので、$\|B\|=\max(|\lambda_1(B)|,|\lambda_n(B)|)$ で、3 が従います。 3 は「摂動 $B$ の大きさ $\|B\|$ 以上には、どの固有値も動かない」ことを意味し、近似固有値の誤差評価（事後誤差評価）の基礎になります。 $A$ を $n$ 次エルミート行列とし、$P$ を $n\times m$ 行列で列が正規直交するもの（$P^{*}P=I_m$）とする。$B=P^{*}AP$ とすると、$1\le k\le m$ について $$ \lambda_k(A)\ge\lambda_k(B)\ge\lambda_{k+n-m}(A). $$ とくに $B$ が $A$ の $m$ 次主部分行列（同じ行と列を選んで残した行列）ならば、これが成り立つ。 $B$ は $m$ 次エルミート行列で、$\langle B\boldsymbol{x},\boldsymbol{x}\rangle=\langle AP\boldsymbol{x},P\boldsymbol{x}\rangle$、$\|P\boldsymbol{x}\|=\|\boldsymbol{x}\|$ です。$P$ は単射なので、$\mathbb{C}^m$ の $k$ 次元部分空間 $S$ は $\mathbb{C}^n$ の $k$ 次元部分空間 $PS$ に写り、$S$ 上の $B$ のレイリー商は $PS$ 上の $A$ のレイリー商と同じ値の集合をとります。 の最大最小より $\lambda_k(B)=\max_{S}\min_{\boldsymbol{x}\in S}\langle B\boldsymbol{x},\boldsymbol{x}\rangle=\max_S\min_{\boldsymbol{y}\in PS}\langle A\boldsymbol{y},\boldsymbol{y}\rangle\le\lambda_k(A)$（$PS$ は $\mathbb{C}^n$ の $k$ 次元部分空間の一部にすぎないので、全部にわたる最大値 $\lambda_k(A)$ を超えない）。最小最大では、$\mathbb{C}^m$ の次元 $m-k+1$ の部分空間 $T$ が $PT$（次元 $m-k+1=n-(k+n-m)+1$）に写るので $\lambda_k(B)=\min_T\max_{\boldsymbol{y}\in PT}\langle A\boldsymbol{y},\boldsymbol{y}\rangle\ge\lambda_{k+n-m}(A)$ です。主部分行列は、選んだ添字に対応する標準基底ベクトルを列にもつ $P$ に対する $P^{*}AP$ です。 $A=\begin{pmatrix}2&1&0\\1&2&1\\0&1&2\end{pmatrix}$ の固有値は $2+\sqrt2,\ 2,\ 2-\sqrt2$ です（特性多項式 $(t-2)^3-2(t-2)=(t-2)((t-2)^2-2)$）。最後の行と列を除いた $B=\begin{pmatrix}2&1\\1&2\end{pmatrix}$ の固有値は $3,1$ で、 $$ 2+\sqrt2\ (\approx3.41)\ \ge\ 3\ \ge\ 2\ \ge\ 1\ \ge\ 2-\sqrt2\ (\approx0.59) $$ と交互に並びます。主部分行列の固有値を調べることで、元の行列の固有値を上下から押さえられるわけです。 ## 4. シュール補行列・ウッドベリーの公式・シルヴェスターの行列式恒等式 ブロック行列 $M=\begin{pmatrix}A&B\\C&D\end{pmatrix}$（$A$ が $p$ 次、$D$ が $q$ 次の正方行列）を、ガウスの消去法をブロックで行うことで分解します。 $P$ を $p$ 次、$R$ を $q$ 次の正方行列、$Q$ を $p\times q$ 行列とすると $$ \det\begin{pmatrix}P&Q\\0&R\end{pmatrix}=\det P\cdot\det R . $$ 転置を考えれば、下三角のブロック行列 $\begin{pmatrix}P&0\\Q&R\end{pmatrix}$ についても同じである。 ライプニッツの公式（）$\sum_\sigma\operatorname{sgn}(\sigma)\prod_jm_{\sigma(j)j}$ で、$j\le p$ のとき $\sigma(j)>p$ なら $m_{\sigma(j)j}=0$（左下ブロックが零）なので、非零の項では $\sigma$ は $\{1,\ldots,p\}$ をそれ自身に写します。このとき $\sigma$ は $\{p+1,\ldots,n\}$ もそれ自身に写し、$\sigma=\sigma_1\times\sigma_2$ と分かれます。逆順対は各部分の中にしかない（$i\le p を使います。 $M=\begin{pmatrix}A&B\\C&D\end{pmatrix}$ とする。 1. $D$ が正則のとき、$S_D=A-BD^{-1}C$（$D$ に関する**シュール補行列**）とおくと $$ M=\begin{pmatrix}I&BD^{-1}\\0&I\end{pmatrix}\begin{pmatrix}S_D&0\\0&D\end{pmatrix}\begin{pmatrix}I&0\\D^{-1}C&I\end{pmatrix} $$ であり、$\det M=\det D\cdot\det S_D$。$M$ が正則であることと $S_D$ が正則であることは同値で、このとき $M^{-1}$ の左上ブロックは $S_D^{-1}$ である。 2. $A$ が正則のとき、$S_A=D-CA^{-1}B$ とおくと $$ M=\begin{pmatrix}I&0\\CA^{-1}&I\end{pmatrix}\begin{pmatrix}A&0\\0&S_A\end{pmatrix}\begin{pmatrix}I&A^{-1}B\\0&I\end{pmatrix},\qquad\det M=\det A\cdot\det S_A, $$ $M$ が正則ならば $M^{-1}$ の左上ブロックは $A^{-1}+A^{-1}BS_A^{-1}CA^{-1}$ である。 3. $M=\begin{pmatrix}A&B\\B^{*}&D\end{pmatrix}$ がエルミート行列で $D$ が正則なら、$M$ の正の固有値・負の固有値の個数は、$D$ と $S_D$ のそれらの和である。とくに $D$ が正定値のとき、$M$ が正定値（半正定値）であることと $S_D$ が正定値（半正定値）であることは同値である。 証明の方針は、1 と 2 では右辺の 3 つの行列を実際に掛けて $M$ に戻ることを確かめ、そこから行列式と逆行列を読み取ることです。3 では、1 の分解が $M=L\,\operatorname{diag}(S_D,D)\,L^{*}$ の形（両側から $L$ と $L^{*}$ で挟んだ形）になっていることを使い、シルヴェスターの慣性法則（この形の変形では正・負の固有値の個数が変わらない、という法則）に結びつけます。 1. 右辺の 3 つの積を計算すると、$\begin{pmatrix}S_D&0\\0&D\end{pmatrix}\begin{pmatrix}I&0\\D^{-1}C&I\end{pmatrix}=\begin{pmatrix}S_D&0\\C&D\end{pmatrix}$、これに左から $\begin{pmatrix}I&BD^{-1}\\0&I\end{pmatrix}$ を掛けると $\begin{pmatrix}S_D+BD^{-1}C&B\\C&D\end{pmatrix}=M$ です。 により、両端の行列の行列式は $1$、真ん中は $\det S_D\det D$ で、積の定理（）から $\det M=\det D\det S_D$ です。$\det D\ne0$ なので $\det M\ne0\iff\det S_D\ne0$（）。逆行列は $M^{-1}=\begin{pmatrix}I&0\\-D^{-1}C&I\end{pmatrix}\begin{pmatrix}S_D^{-1}&0\\0&D^{-1}\end{pmatrix}\begin{pmatrix}I&-BD^{-1}\\0&I\end{pmatrix}$ で、左上ブロックは $S_D^{-1}$ です。 2. 同様に積を計算して確かめられます（$\begin{pmatrix}A&0\\0&S_A\end{pmatrix}\begin{pmatrix}I&A^{-1}B\\0&I\end{pmatrix}=\begin{pmatrix}A&B\\0&S_A\end{pmatrix}$、これに左から $\begin{pmatrix}I&0\\CA^{-1}&I\end{pmatrix}$ を掛けると $\begin{pmatrix}A&B\\C&CA^{-1}B+S_A\end{pmatrix}=M$）。逆行列は $M^{-1}=\begin{pmatrix}I&-A^{-1}B\\0&I\end{pmatrix}\begin{pmatrix}A^{-1}&0\\0&S_A^{-1}\end{pmatrix}\begin{pmatrix}I&0\\-CA^{-1}&I\end{pmatrix}$ の左上ブロック $A^{-1}+A^{-1}BS_A^{-1}CA^{-1}$ です。 3. エルミートなら $C=B^{*}$、$D=D^{*}$ で、$L=\begin{pmatrix}I&BD^{-1}\\0&I\end{pmatrix}$ とおくと 1 の分解は $M=L\begin{pmatrix}S_D&0\\0&D\end{pmatrix}L^{*}$ です（$L^{*}=\begin{pmatrix}I&0\\D^{-1}B^{*}&I\end{pmatrix}$）。$M$ は $\operatorname{diag}(S_D,D)$ と合同なので、シルヴェスターの慣性法則（）により、符号の個数が等しく、ブロック対角行列の符号の個数は各ブロックの和です。$D$ が正定値のときは、$D$ の寄与が $n_+(D)=q$ のみなので、$M$ が正定値（$n_+(M)=p+q$）であることは $n_+(S_D)=p$ と同値です。半正定値の場合は $\langle M\boldsymbol{x},\boldsymbol{x}\rangle=\langle\operatorname{diag}(S_D,D)L^{*}\boldsymbol{x},L^{*}\boldsymbol{x}\rangle$ で、$L^{*}$ が全単射であることから、$M\ge0\iff\operatorname{diag}(S_D,D)\ge0\iff S_D\ge0$ です。 シュール補行列の意味は、$M\begin{pmatrix}\boldsymbol{x}\\\boldsymbol{y}\end{pmatrix}=\begin{pmatrix}\boldsymbol{b}\\\boldsymbol{c}\end{pmatrix}$ という連立方程式から $\boldsymbol{y}$ を消去して残る $\boldsymbol{x}$ だけの方程式 $S_D\boldsymbol{x}=\boldsymbol{b}-BD^{-1}\boldsymbol{c}$ の係数行列、ということです。 $A$（$n$ 次）、$C$（$k$ 次）を正則行列、$U$ を $n\times k$、$V$ を $k\times n$ とする。このとき $A+UCV$ が正則であることと $C^{-1}+VA^{-1}U$ が正則であることは同値で、そのとき $$ (A+UCV)^{-1}=A^{-1}-A^{-1}U\bigl(C^{-1}+VA^{-1}U\bigr)^{-1}VA^{-1}, $$ $$ \det(A+UCV)=\det A\cdot\det C\cdot\det\bigl(C^{-1}+VA^{-1}U\bigr) $$ が成り立つ。 証明のねらいは、大きな行列 $N$ を 1 つ用意して、シュール補行列を「$D$ を消す側」と「$A$ を消す側」の 2 通りで計算し、結果を等しいとおくことです。$N$ は、一方の消去で $A+UCV$、もう一方で $C^{-1}+VA^{-1}U$ が現れるように選んであります。 $N=\begin{pmatrix}A&-U\\V&C^{-1}\end{pmatrix}$ とおきます。 の 1 を $D=C^{-1}$ として使うと、$S_D=A-(-U)C\,V=A+UCV$ で、$\det N=\det C^{-1}\det(A+UCV)$、$N$ の正則性は $A+UCV$ のそれと同値、$N^{-1}$ の左上ブロックは $(A+UCV)^{-1}$ です。2 を使うと、$S_A=C^{-1}-VA^{-1}(-U)=C^{-1}+VA^{-1}U$ で、$\det N=\det A\det(C^{-1}+VA^{-1}U)$、$N^{-1}$ の左上ブロックは $A^{-1}+A^{-1}(-U)S_A^{-1}VA^{-1}$ です。$N$ の正則性の同値と、行列式・左上ブロックの二通りの表示を等置すれば、主張が得られます（$\det C^{-1}=1/\det C$）。 $k=1$、$C=1$、$U=\boldsymbol{u}$、$V=\boldsymbol{v}^{*}$ のとき（シャーマン・モリソンの公式）、$1+\boldsymbol{v}^{*}A^{-1}\boldsymbol{u}\ne0$ ならば $$ (A+\boldsymbol{u}\boldsymbol{v}^{*})^{-1}=A^{-1}-\frac{A^{-1}\boldsymbol{u}\boldsymbol{v}^{*}A^{-1}}{1+\boldsymbol{v}^{*}A^{-1}\boldsymbol{u}},\qquad\det(A+\boldsymbol{u}\boldsymbol{v}^{*})=\det A\,(1+\boldsymbol{v}^{*}A^{-1}\boldsymbol{u}). $$ $A=I_2$、$\boldsymbol{u}=(1,2)^{\mathsf{T}}$、$\boldsymbol{v}=(1,1)^{\mathsf{T}}$ なら $I+\boldsymbol{u}\boldsymbol{v}^{\mathsf{T}}=\begin{pmatrix}2&1\\2&3\end{pmatrix}$ で、$1+\boldsymbol{v}^{\mathsf{T}}\boldsymbol{u}=4=\det$（$2\cdot3-1\cdot2=4$）。逆行列は $I-\tfrac14\begin{pmatrix}1&1\\2&2\end{pmatrix}=\begin{pmatrix}3/4&-1/4\\-1/2&1/2\end{pmatrix}$ で、元の行列との積が $I$ になります。 $A$ の逆行列が既知で $UCV$ が低ランクなら、$n$ 次の逆行列を $k$ 次の逆行列に帰着できるので、逐次更新（カルマンフィルタ、最小二乗法の逐次版）で使われます。 $X$ を $m\times n$、$Y$ を $n\times m$ の行列とする。 $$ \det(I_m+XY)=\det(I_n+YX). $$ さらに $t^n\det(tI_m-XY)=t^m\det(tI_n-YX)$。すなわち $XY$ と $YX$ の $0$ でない固有値は、重複度も込めて一致する。 $N=\begin{pmatrix}I_m&-X\\Y&I_n\end{pmatrix}$ とします。 の 1（$D=I_n$）で $S_D=I_m+XY$、2（$A=I_m$）で $S_A=I_n+YX$ なので、$\det N=\det(I_m+XY)=\det(I_n+YX)$ です。$t\ne0$ のとき、$X$ と $-Y/t$ に適用して $\det(I_m-XY/t)=\det(I_n-YX/t)$ で、$\det(tI_m-XY)=t^m\det(I_m-XY/t)$、$\det(tI_n-YX)=t^n\det(I_n-YX/t)$ だから $t^{n}\det(tI_m-XY)=t^{m}\det(tI_n-YX)$ が $t\ne0$ で成り立ちます。両辺は $t$ の多項式で、無限個の $t$ で一致するので多項式として等しいです。 ## 5. アダマールの不等式とコーシー・ビネの公式 1. $A=(a_{ij})$ が $n$ 次の半正定値エルミート行列ならば $\det A\le a_{11}a_{22}\cdots a_{nn}$。$A$ が正定値のとき等号は $A$ が対角行列のときに限る。 2. 一般の $A\in M_n(\mathbb{C})$ の列ベクトルを $\boldsymbol{a}_1,\ldots,\boldsymbol{a}_n$ とすると $|\det A|\le\prod_{j=1}^{n}\|\boldsymbol{a}_j\|$。 1. まず $A$ が正定値の場合を $n$ についての帰納法で示します。$n=1$ は自明です。$A=\begin{pmatrix}A'&\boldsymbol{b}\\\boldsymbol{b}^{*}&a_{nn}\end{pmatrix}$ と分けると、$A'$（主部分行列）は正定値です（$\langle A'\boldsymbol{x},\boldsymbol{x}\rangle=\langle A(\boldsymbol{x},0)^{\mathsf{T}},(\boldsymbol{x},0)^{\mathsf{T}}\rangle>0$）。 の 2 より $\det A=\det A'\cdot(a_{nn}-\boldsymbol{b}^{*}A'^{-1}\boldsymbol{b})$。$A'^{-1}$ も正定値（固有値が $A'$ の固有値の逆数）なので $\boldsymbol{b}^{*}A'^{-1}\boldsymbol{b}\ge0$ で、$\det A'>0$ とあわせて $\det A\le a_{nn}\det A'\le a_{nn}\prod_{i0$）が正定値で、$\det(A+\varepsilon I)\le\prod(a_{ii}+\varepsilon)$ の $\varepsilon\to0$ の極限（多項式は連続）で得られます。 2. $A^{*}A$ は半正定値エルミートで、対角成分は $\|\boldsymbol{a}_j\|^2$、行列式は $|\det A|^2$ です（、$\det A^{*}=\overline{\det A}$）。1 を $A^{*}A$ に適用します。 幾何的には、「各辺の長さが決まっている平行体の中で、体積が最大になるのは辺が直交するときである」という主張です。同値な言い方で、$|\det A|$ は列の長さの積以下です。 $A$ を $m\times n$、$B$ を $n\times m$ 行列とする。$m>n$ ならば $\det(AB)=0$。$m\le n$ ならば $$ \det(AB)=\sum_{S}\det A_{[m],S}\cdot\det B_{S,[m]} $$ が成り立つ。和は $\{1,\ldots,n\}$ の $m$ 個の元からなる部分集合 $S$ 全体にわたり、$A_{[m],S}$ は $A$ から列 $S$ を、$B_{S,[m]}$ は $B$ から行 $S$ を選んだ $m$ 次正方行列である。 $m>n$ のとき、$AB$ の像は $A$ の像に含まれ、$A$ の像は $n$ 本の列で張られるので次元は $n$ 以下です。$AB$ が正則なら像が $\mathbb{C}^m$ 全体になり矛盾するので、$AB$ は正則でなく、行列式は $0$ です（）。 $m\le n$ とします。$A$ の第 $k$ 列を $\boldsymbol{A}_k$ と書くと、$AB$ の第 $j$ 列は $\sum_{k=1}^{n}b_{kj}\boldsymbol{A}_k$ です。行列式は列について多重線形なので（）、 $$ \det(AB)=\sum_{k_1,\ldots,k_m=1}^{n}b_{k_11}b_{k_22}\cdots b_{k_mm}\,\det(\boldsymbol{A}_{k_1},\ldots,\boldsymbol{A}_{k_m}) $$ です。$k_1,\ldots,k_m$ に重複があれば等しい列があるので行列式は $0$ です。重複がなければ、$S=\{k_1,\ldots,k_m\}=\{s_1<\cdots の 4）、$\det(\boldsymbol{A}_{k_1},\ldots,\boldsymbol{A}_{k_m})=\operatorname{sgn}(\tau)\det A_{[m],S}$ です。よって $$ \det(AB)=\sum_S\det A_{[m],S}\sum_{\tau\in S_m}\operatorname{sgn}(\tau)\prod_{j=1}^{m}b_{s_{\tau(j)}j}=\sum_S\det A_{[m],S}\det B_{S,[m]} $$ です（最後の等号は $B_{S,[m]}$ の $(i,j)$ 成分が $b_{s_ij}$ であることとライプニッツの公式）。 $A$ を $n\times m$ 行列（$m\le n$）とすると $$ \det(A^{*}A)=\sum_{S}|\det A_{S,[m]}|^2 $$ （和は行の $m$ 個の部分集合）。とくに $A$ の階数が $m$ であることと、$0$ でない $m$ 次小行列式が存在することは同値である。 を $A^{*}$ と $A$ に適用すると、$(A^{*})_{[m],S}=(A_{S,[m]})^{*}$ なので $\det(A^{*}A)=\sum_S\overline{\det A_{S,[m]}}\det A_{S,[m]}$ です。$\det(A^{*}A)\ne0$ は $A^{*}A$ が正則であること、すなわち $A$ の列が一次独立であることと同値です（$A^{*}A\boldsymbol{x}=0\Rightarrow\|A\boldsymbol{x}\|^2=0$）。 $A=\begin{pmatrix}1&4\\2&5\\3&6\end{pmatrix}$（列 $\boldsymbol{v}=(1,2,3)$、$\boldsymbol{w}=(4,5,6)$）とすると、$\det A^{\mathsf{T}}A=\|\boldsymbol{v}\|^2\|\boldsymbol{w}\|^2-(\boldsymbol{v}\cdot\boldsymbol{w})^2=14\cdot77-32^2=54$ です。一方、$2\times2$ 小行列式（行 $\{1,2\},\{1,3\},\{2,3\}$）は $-3,-6,-3$ で、平方の和は $9+36+9=54$ と一致します。これは $\|\boldsymbol{v}\times\boldsymbol{w}\|^2$ にほかなりません（ の係数 $-3,-6,-3$ と同じ小行列式）。また、$A^{\mathsf{T}}A$ にアダマールの不等式（ の 1）を使うと $54\le14\cdot77$ です。 ## 6. ペロン・フロベニウスの定理 成分がすべて非負の実行列を**非負行列**といい、$A\ge0$ と書きます。実ベクトルについても $\boldsymbol{x}\ge0$（すべての成分が $\ge0$）、$\boldsymbol{x}>0$（すべて $>0$）、$\boldsymbol{x}\ge\boldsymbol{y}$（$\boldsymbol{x}-\boldsymbol{y}\ge0$）と書き、行列の正値性 $A>0$ は全成分が正であることを表します。非負行列は、マルコフ連鎖の遷移確率、ネットワークの隣接行列、人口モデルなどに現れます。非負性だけから固有値の強い性質が出ます。 $A=(a_{ij})\ge0$ が**既約**であるとは、任意の $i,j$ に対し、$i=i_0,i_1,\ldots,i_k=j$ で $a_{i_{s-1}i_s}>0$ がすべての $s$ で成り立つものが存在する（$k=0$ を許す）ことをいう。 既約性は「[有向グラフ](/mathematics/graph-theory/connectivity-and-flows#def-network)（$a_{ij}>0$ のとき $i\to j$ に辺を引く）が強連結」という条件です。
123 123 既約（強連結） 既約でない
$A\ge0$ を $n\ge2$ 次の既約行列とすると、$B=(I+A)^{n-1}$ は $B>0$ である。さらに、$\boldsymbol{x}\ge0$、$\boldsymbol{x}\ne0$ ならば $B\boldsymbol{x}>0$ であり、$A$ の各行は零でない。 $(A^k)_{ij}=\sum a_{ii_1}a_{i_1i_2}\cdots a_{i_{k-1}j}$ は非負の項の和なので、$(A^k)_{ij}>0$ であることは、$i$ から $j$ への長さ $k$ の道（辺 $a_{\cdot\cdot}>0$ をたどる列）が存在することと同値です。$i\ne j$ のとき最短の道は同じ頂点を 2 度通らないので長さは $n-1$ 以下で、その長さ $k$ に対し $(A^k)_{ij}>0$ です。$B=\sum_{k=0}^{n-1}\binom{n-1}{k}A^k$ は非負行列の和なので $B_{ij}\ge\binom{n-1}{k}(A^k)_{ij}>0$。$i=j$ では $k=0$ の項の $1$ があります。$\boldsymbol{x}\ge0$ が正の成分をもてば $(B\boldsymbol{x})_i=\sum_jB_{ij}x_j>0$ です。既約なら各 $i$ から別の頂点への道があるので、$a_{ij}>0$ となる $j$ が存在し、行は零ではありません。 $\boldsymbol{x}\ge0$、$\boldsymbol{x}\ne0$ に対し、 $$ r_A(\boldsymbol{x})=\min_{i:\,x_i>0}\frac{(A\boldsymbol{x})_i}{x_i}=\max\{t\ge0\mid A\boldsymbol{x}\ge t\boldsymbol{x}\} $$ と定めます（二つの表示が等しいのは、$x_i=0$ の行では $(A\boldsymbol{x})_i\ge0=tx_i$ が自動的に成り立つからです）。$r_A(c\boldsymbol{x})=r_A(\boldsymbol{x})$（$c>0$）です。これを**コリャッツ・ヴィーラントの関数**といいます。 $A\ge0$ を $n\ge2$ 次の既約行列とし、$\rho=\rho(A)=\max\{|\lambda|:\lambda\text{ は }A\text{ の固有値}\}$ を**スペクトル半径**とする。 1. $\rho>0$ は $A$ の固有値であり、成分がすべて正の固有ベクトル $\boldsymbol{y}>0$ が存在する。 2. （コリャッツ・ヴィーラント） $$ \rho=\max_{\boldsymbol{x}\ge0,\ \boldsymbol{x}\ne0}r_A(\boldsymbol{x})=\min_{\boldsymbol{x}>0}\ \max_i\frac{(A\boldsymbol{x})_i}{x_i}. $$ 3. $\rho$ は $A$ の特性多項式の単純根（代数的重複度 $1$）で、固有空間は $\boldsymbol{y}$ が張る 1 次元空間である。非負の固有ベクトル（固有値は任意）は $\boldsymbol{y}$ の正の定数倍に限る。 この証明は 5 つの段に分かれています。第 1 段で、$r_A$ を最大にするベクトル $\boldsymbol{y}_0$ の存在を、コンパクト性を使って示します（$B$ を掛けて正のベクトルだけの集合に移すのは、そこで $r_A$ が連続になり、最大値が存在するようにするためです）。第 2 段では、最大値 $L$ の「これ以上大きくできない」性質から、$A\boldsymbol{x}\ge L\boldsymbol{x}$ なら等号が成り立つことを示します。第 3 段で $L=\rho$ と 1 を、第 4 段で 2 の第 2 の等号を、第 5 段で 3 を示します。 $B=(I+A)^{n-1}>0$（）とし、$\Delta=\{\boldsymbol{x}\ge0\mid\sum x_i=1\}$（有界閉集合なのでによりコンパクト）とします。 **第 1 段（$r_A$ の最大の存在）。** $\boldsymbol{x}\in\Delta$ について $A\boldsymbol{x}\ge r_A(\boldsymbol{x})\boldsymbol{x}$ の両辺に $B\ge0$ を掛けると、$AB=BA$ より $A(B\boldsymbol{x})\ge r_A(\boldsymbol{x})B\boldsymbol{x}$、すなわち $r_A(\boldsymbol{x})\le r_A(B\boldsymbol{x})$ です。$N=B(\Delta)$ はコンパクト集合の連続像でコンパクト（）、その元はすべて正ベクトル（）で、正ベクトルの集合上では $r_A(\boldsymbol{y})=\min_i(A\boldsymbol{y})_i/y_i$ は連続です。ゆえに $r_A$ は $N$ 上で最大値 $L=r_A(\boldsymbol{y}_0)$（$\boldsymbol{y}_0\in N$）をとります。$r_A$ はスケール不変なので、$L=\sup_{\boldsymbol{x}\ge0,\boldsymbol{x}\ne0}r_A(\boldsymbol{x})$ で、$\boldsymbol{x}\in\Delta$ では $r_A(\boldsymbol{x})\le r_A(B\boldsymbol{x})\le L$、$\boldsymbol{y}_0$ で等号です。 **第 2 段（$\star$）$\boldsymbol{x}\ge0$、$\boldsymbol{x}\ne0$、$A\boldsymbol{x}\ge L\boldsymbol{x}$ ならば $A\boldsymbol{x}=L\boldsymbol{x}$ かつ $\boldsymbol{x}>0$。** $\boldsymbol{z}=A\boldsymbol{x}-L\boldsymbol{x}\ge0$ が $\boldsymbol{0}$ でないと仮定します。$\boldsymbol{w}=B\boldsymbol{x}>0$ について $A\boldsymbol{w}=BA\boldsymbol{x}=L\boldsymbol{w}+B\boldsymbol{z}$ で、$B\boldsymbol{z}>0$ です。$\varepsilon=\min_i(B\boldsymbol{z})_i/w_i>0$ とおけば $A\boldsymbol{w}\ge(L+\varepsilon)\boldsymbol{w}$ で $r_A(\boldsymbol{w})\ge L+\varepsilon>L$ となり、$L$ の最大性に反します。ゆえに $\boldsymbol{z}=0$。このとき $B\boldsymbol{x}=(1+L)^{n-1}\boldsymbol{x}$ で、左辺は正なので $\boldsymbol{x}>0$ です。 **第 3 段（1）。** $r_A(\boldsymbol{y}_0)=L$ は $A\boldsymbol{y}_0\ge L\boldsymbol{y}_0$ を意味するので、（$\star$）から $A\boldsymbol{y}_0=L\boldsymbol{y}_0$、$\boldsymbol{y}_0>0$。$A$ の行は零でなく $\boldsymbol{y}_0>0$ なので $(A\boldsymbol{y}_0)_i>0$、したがって $L>0$ です。次に $L=\rho$ を示します。$L$ は固有値なので $L\le\rho$。逆に $A\boldsymbol{x}=\lambda\boldsymbol{x}$（$\boldsymbol{x}\ne0$）とすると、$|\lambda||x_i|=|\sum_ja_{ij}x_j|\le\sum_ja_{ij}|x_j|$ より $A|\boldsymbol{x}|\ge|\lambda||\boldsymbol{x}|$（$|\boldsymbol{x}|$ は成分ごとの絶対値）、すなわち $r_A(|\boldsymbol{x}|)\ge|\lambda|$ で $|\lambda|\le L$ です。よって $\rho=L$、$\boldsymbol{y}=\boldsymbol{y}_0$ とすれば 1 が示され、2 の第 1 の等号も従います。 **第 4 段（2 の第 2 の等号）。** $A^{\mathsf{T}}$ も非負の既約行列（道の向きを逆にする）で、$A$ と同じ特性多項式をもつので $\rho(A^{\mathsf{T}})=\rho$、第 3 段より $A^{\mathsf{T}}\boldsymbol{w}=\rho\boldsymbol{w}$ となる $\boldsymbol{w}>0$ があります。$\boldsymbol{x}>0$ で $\max_i(A\boldsymbol{x})_i/x_i=s$ とすると $A\boldsymbol{x}\le s\boldsymbol{x}$ で、$\rho\,\boldsymbol{w}^{\mathsf{T}}\boldsymbol{x}=\boldsymbol{w}^{\mathsf{T}}A\boldsymbol{x}\le s\,\boldsymbol{w}^{\mathsf{T}}\boldsymbol{x}$、$\boldsymbol{w}^{\mathsf{T}}\boldsymbol{x}>0$ から $\rho\le s$。$\boldsymbol{x}=\boldsymbol{y}$ では $s=\rho$ なので、最小は $\rho$ です。 **第 5 段（3）。** 実の固有ベクトル $\boldsymbol{x}$（固有値 $\rho$）をとり、$t=\min_ix_i/y_i$、$\boldsymbol{v}=\boldsymbol{x}-t\boldsymbol{y}\ge0$ とします。$A\boldsymbol{v}=\rho\boldsymbol{v}$ で、$\boldsymbol{v}$ にはある成分が $0$ です。$\boldsymbol{v}\ne0$ なら $\boldsymbol{v}\ge0$ に（$\star$）（あるいは $B\boldsymbol{v}=(1+\rho)^{n-1}\boldsymbol{v}>0$）を使って全成分が正となり矛盾するので、$\boldsymbol{v}=0$、$\boldsymbol{x}=t\boldsymbol{y}$ です。複素ベクトルは実部・虚部（どちらも固有ベクトルまたは $0$）に分けて、固有空間は $\boldsymbol{y}$ の張る空間です（幾何的重複度 1）。非負の固有ベクトル $\boldsymbol{x}\ge0$（固有値 $\lambda$）は、$\rho\,\boldsymbol{w}^{\mathsf{T}}\boldsymbol{x}=\boldsymbol{w}^{\mathsf{T}}A\boldsymbol{x}=\lambda\boldsymbol{w}^{\mathsf{T}}\boldsymbol{x}$、$\boldsymbol{w}^{\mathsf{T}}\boldsymbol{x}>0$ から $\lambda=\rho$ で、したがって $\boldsymbol{y}$ の倍数です。 代数的重複度: ジョルダン標準形（）で、固有値 $\rho$ の細胞の個数は幾何的重複度 $1$ に等しいので、$\rho$ の細胞は 1 つ、その大きさが代数的重複度です。大きさが 2 以上なら $(A-\rho I)^2\boldsymbol{v}=0$ かつ $(A-\rho I)\boldsymbol{v}\ne0$ となる $\boldsymbol{v}$ があり、$(A-\rho I)\boldsymbol{v}$ は固有空間の元なので $c\boldsymbol{y}$（$c\ne0$）です。左から $\boldsymbol{w}^{\mathsf{T}}$ を掛けると $0=\boldsymbol{w}^{\mathsf{T}}(A-\rho I)\boldsymbol{v}=c\,\boldsymbol{w}^{\mathsf{T}}\boldsymbol{y}\ne0$ となり矛盾します。よって大きさは 1、$\rho$ は単純根です。 $A>0$（全成分が正）ならば、$A$ は既約で、$\rho$ 以外の固有値 $\lambda$ は $|\lambda|<\rho$ を満たす。 $A>0$ なら $a_{ij}>0$ なので道の長さ 1 で既約です。$A\boldsymbol{x}=\lambda\boldsymbol{x}$、$|\lambda|=\rho$ とします。$A|\boldsymbol{x}|\ge|\lambda||\boldsymbol{x}|=\rho|\boldsymbol{x}|$ に の証明の（$\star$）（$L=\rho$）を使うと $A|\boldsymbol{x}|=\rho|\boldsymbol{x}|$、$|\boldsymbol{x}|>0$ です。したがって各 $i$ で $|\sum_ja_{ij}x_j|=\rho|x_i|=\sum_ja_{ij}|x_j|$、すなわち複素数 $c_j=a_{ij}x_j$ について $|\sum c_j|=\sum|c_j|$ で、すべて $c_j\ne0$ です。$c_j=r_je^{i\theta_j}$ とおくと $|\sum c_j|^2=\sum_{j,k}r_jr_k\cos(\theta_j-\theta_k)\le(\sum r_j)^2$ で、等号は全ての $\cos(\theta_j-\theta_k)=1$ のとき、つまり偏角がすべて等しいときだけです。よって $x_j=e^{i\theta}|x_j|$（$\theta$ は共通）、$\boldsymbol{x}=e^{i\theta}|\boldsymbol{x}|$ で $A\boldsymbol{x}=e^{i\theta}A|\boldsymbol{x}|=\rho\boldsymbol{x}$、すなわち $\lambda=\rho$ です。 $A=\begin{pmatrix}1&1\\1&0\end{pmatrix}$ は既約です。$\boldsymbol{x}=(3,2)^{\mathsf{T}}>0$ に対し $A\boldsymbol{x}=(5,3)^{\mathsf{T}}$、比 $5/3\approx1.667$、$3/2=1.5$ なので $$ 1.5\le\rho(A)\le1.667 $$ です（ の 2）。真の値は黄金比 $\frac{1+\sqrt5}{2}\approx1.618$ で、対応する固有ベクトルは $(\varphi,1)$ です。固有値を計算せず、ベクトル 1 本の比だけで上下から挟めるのがこの公式の利点で、$\boldsymbol{x}$ を $A\boldsymbol{x}$ に更新する（冪乗法）たびに挟み込みが狭まります。 マルコフ連鎖の遷移行列 $P\ge0$ が既約で各行の和が $1$ なら、$P\boldsymbol{1}=\boldsymbol{1}$ で正の固有ベクトル $\boldsymbol{1}$ があるので、 の 3 により固有値 $1$ が $\rho(P)$ で、単純です。定常分布は $P^{\mathsf{T}}$ の対応する正の固有ベクトルとして一意に存在します。 ## 7. シルヴェスター方程式とリャプノフ方程式 $A\in M_m(\mathbb{C})$、$B\in M_n(\mathbb{C})$ とし、$\mathcal L:M_{m,n}(\mathbb{C})\to M_{m,n}(\mathbb{C})$ を $\mathcal L(X)=AX+XB$ とする。次は同値である。 1. $\mathcal L$ は全単射、すなわちすべての $C$ に対し $AX+XB=C$ はただ一つの解 $X$ をもつ。 2. $A$ の固有値と $B$ の固有値の符号を変えたもの $-\mu$ に共通するものがない。すなわち、$A$ の固有値 $\lambda$ と $B$ の固有値 $\mu$ で $\lambda+\mu=0$ となるものが存在しない。 $\mathcal L$ は有限次元空間 $M_{m,n}$ の自分自身への線形写像なので、全単射であることと単射であることは同値です（次元定理）。 (2)$\Rightarrow$(1): $AX+XB=0$ とします。帰納的に $A^kX=X(-B)^k$ ですから、多項式 $p$ に対し $p(A)X=Xp(-B)$ です。$p=\varphi_A$ とすると、ケイリー・ハミルトンの定理（）で $\varphi_A(A)=0$ なので $X\varphi_A(-B)=0$ です。$\varphi_A(t)=\prod_i(t-\lambda_i)$（$\lambda_i$ は $A$ の固有値）と因数分解して $\varphi_A(-B)=\prod_i(-B-\lambda_iI)=(-1)^m\prod_i(B+\lambda_iI)$ です。各因子 $B+\lambda_iI$ は、$-\lambda_i$ が $B$ の固有値でないので正則です。積は正則なので $X=0$ です。 (1)$\Rightarrow$(2)（対偶）: $\lambda$ が $A$ の固有値、$-\lambda$ が $B$ の固有値とします。$A\boldsymbol{x}=\lambda\boldsymbol{x}$（$\boldsymbol{x}\ne0$）で、$B^{\mathsf{T}}$ も $B$ と同じ固有値をもつので $B^{\mathsf{T}}\boldsymbol{\eta}=-\lambda\boldsymbol{\eta}$（$\boldsymbol{\eta}\ne0$）、すなわち $\boldsymbol{\eta}^{\mathsf{T}}B=-\lambda\boldsymbol{\eta}^{\mathsf{T}}$ となる $\boldsymbol{\eta}$ があります。$X=\boldsymbol{x}\boldsymbol{\eta}^{\mathsf{T}}\ne0$ に対し $AX+XB=\lambda X-\lambda X=0$ で、$\mathcal L$ は単射でありません。 $B=-A$ なら固有値の組 $\lambda$、$-\mu=\lambda$ が必ず共通なので $AX-XA=C$ は一般に可解でなく、$AX=XA$ の解（$A$ と可換な行列）が非自明に存在します。$B=A^{*}$（リャプノフ方程式の場合）なら条件は $\lambda_i+\overline{\lambda_j}\ne0$ です。 連続時間の線形系 $\dot{\boldsymbol{x}}=A\boldsymbol{x}$ が安定（すべての解が $t\to\infty$ で $0$ に収束）であるための条件は、$A$ の固有値の実部がすべて負であることです。これをリャプノフ方程式 $A^{*}X+XA=-Q$ で判定します。まず、離散版（スタイン方程式）を片づけます。 $C\in M_n(\mathbb{C})$ のスペクトル半径が $\rho(C)<1$ ならば、$c>0$ と $\theta<1$ が存在して、すべての $k$ で $\|C^k\|\le c\,\theta^k$。 シュール分解 $C=UTU^{*}$（）をとります。$T$ の対角成分は固有値で、絶対値は $\rho(C)=:\rho$ 以下です。$\tau=\max_i\sum_{j>i}|t_{ij}|$ とし、$0<\varepsilon\le1$ を $\theta:=\rho+\varepsilon\tau<1$ となるように選びます。$D_\varepsilon=\operatorname{diag}(1,\varepsilon,\ldots,\varepsilon^{n-1})$ に対し $T_\varepsilon=D_\varepsilon^{-1}TD_\varepsilon$ の $(i,j)$ 成分は $\varepsilon^{j-i}t_{ij}$ で、行ごとの絶対値の和は $|t_{ii}|+\sum_{j>i}\varepsilon^{j-i}|t_{ij}|\le\rho+\varepsilon\tau=\theta$ です。行和ノルム $\|X\|_\infty=\max_i\sum_j|x_{ij}|$ は積について劣乗法的（$\sum_j|\sum_kx_{ik}y_{kj}|\le\sum_k|x_{ik}|\|Y\|_\infty\le\|X\|_\infty\|Y\|_\infty$）なので、$\|T_\varepsilon^k\|_\infty\le\theta^k$ です。$T^k=D_\varepsilon T_\varepsilon^kD_\varepsilon^{-1}$ の成分は $|(T^k)_{ij}|=\varepsilon^{i-j}|(T_\varepsilon^k)_{ij}|\le\varepsilon^{-(n-1)}\theta^k$ で、作用素ノルムはフロベニウスノルム以下（コーシー・シュワルツ）、ユニタリ変換で不変なので $\|C^k\|=\|T^k\|\le n\,\varepsilon^{-(n-1)}\theta^k$ です。 $C\in M_n(\mathbb{C})$ が $\rho(C)<1$ を満たすとする。任意の正定値エルミート行列 $Q'$ に対し、方程式 $$ X-C^{*}XC=Q' $$ はただ一つの解 $X$ をもち、$X=\sum_{k=0}^{\infty}(C^{*})^kQ'C^k$ は正定値エルミート行列である。 より $\|(C^{*})^kQ'C^k\|\le c^2\theta^{2k}\|Q'\|$ で級数は収束し、各項は半正定値エルミートです。$k=0$ の項が $Q'>0$ なので和 $X$ は正定値エルミートです（$\langle X\boldsymbol{x},\boldsymbol{x}\rangle\ge\langle Q'\boldsymbol{x},\boldsymbol{x}\rangle>0$）。項別に $C^{*}XC=\sum_{k\ge1}(C^{*})^kQ'C^k=X-Q'$ です。一意性: $X_1,X_2$ が解なら差 $Y$ は $Y=C^{*}YC$ を満たし、繰り返して $Y=(C^{*})^kYC^k\to0$、$Y=0$ です。 $A\in M_n(\mathbb{C})$ について次は同値である。 1. $A$ のすべての固有値の実部が負である。 2. 任意の正定値エルミート行列 $Q$ に対し、$A^{*}X+XA=-Q$ はただ一つの解 $X$ をもち、それは正定値エルミート行列である。 3. ある正定値エルミート行列 $Q,X$ について $A^{*}X+XA=-Q$。 証明の要点は、(1)$\Rightarrow$(2) でケイリー変換 $C=(A-I)^{-1}(A+I)$ を使うことです。これは「固有値の実部が負」という条件を「固有値の絶対値が $1$ 未満」に写す変換で、リャプノフ方程式を直前のスタイン方程式に書き換えられます。そうすれば解の存在・一意性・正定値性は から従います。 (2)$\Rightarrow$(3) は $Q=I$ とすれば明らか。 (3)$\Rightarrow$(1): $A\boldsymbol{x}=\lambda\boldsymbol{x}$（$\boldsymbol{x}\ne0$）とすると $\boldsymbol{x}^{*}A^{*}=\bar\lambda\boldsymbol{x}^{*}$ で、 $$ -\boldsymbol{x}^{*}Q\boldsymbol{x}=\boldsymbol{x}^{*}(A^{*}X+XA)\boldsymbol{x}=(\bar\lambda+\lambda)\boldsymbol{x}^{*}X\boldsymbol{x}=2\operatorname{Re}\lambda\cdot\boldsymbol{x}^{*}X\boldsymbol{x} $$ です。左辺は負、$\boldsymbol{x}^{*}X\boldsymbol{x}>0$ なので $\operatorname{Re}\lambda<0$ です。 (1)$\Rightarrow$(2): 実部が負なら $1$ は固有値でなく $A-I$ は正則です。ケイリー変換 $C=(A-I)^{-1}(A+I)$ をとります。$A=UTU^{*}$（シュール分解）とすると $C=U(T-I)^{-1}(T+I)U^{*}$ で、上三角行列の逆と積は上三角で対角成分が $\frac{t_{ii}+1}{t_{ii}-1}$ ですから、$C$ の固有値は $\frac{\lambda+1}{\lambda-1}$ です。$|\lambda+1|^2-|\lambda-1|^2=4\operatorname{Re}\lambda<0$ なので、これらの絶対値は $1$ 未満、$\rho(C)<1$ です。また $C-I=2(A-I)^{-1}$ は正則で、$A=(C+I)(C-I)^{-1}$ すなわち $A(C-I)=C+I$ が成り立ちます。この式と随伴 $(C-I)^{*}A^{*}=(C+I)^{*}$ を使うと、任意の行列 $X$ に対し $$ (C-I)^{*}(A^{*}X+XA)(C-I)=(C^{*}+I)X(C-I)+(C^{*}-I)X(C+I)=2\bigl(C^{*}XC-X\bigr) $$ です。したがって $A^{*}X+XA=-Q$ は、$C-I$ が正則であることから、$X-C^{*}XC=Q'$、$Q'=\frac12(C-I)^{*}Q(C-I)$ と同値で、$Q>0$ なら $\langle Q'\boldsymbol{x},\boldsymbol{x}\rangle=\frac12\langle Q(C-I)\boldsymbol{x},(C-I)\boldsymbol{x}\rangle>0$ より $Q'>0$ です。 が $Q'$ に対し、一意の解 $X$ とその正定値性を与えます。 $A=\begin{pmatrix}-1&2\\0&-3\end{pmatrix}$（固有値 $-1,-3$）、$Q=I$ とし、実対称行列 $X=\begin{pmatrix}a&b\\b&c\end{pmatrix}$ で $A^{\mathsf{T}}X+XA=-I$ を解きます。$XA=\begin{pmatrix}-a&2a-3b\\-b&2b-3c\end{pmatrix}$、$A^{\mathsf{T}}X=(XA)^{\mathsf{T}}$ なので、和は $\begin{pmatrix}-2a&2a-4b\\2a-4b&4b-6c\end{pmatrix}=-I$ です。$a=\tfrac12$、$b=\tfrac14$、$c=\tfrac13$。$X=\begin{pmatrix}1/2&1/4\\1/4&1/3\end{pmatrix}$ の行列式は $\tfrac16-\tfrac1{16}=\tfrac5{48}>0$ で $a>0$ ですから正定値で、 の結論どおりです。 ## 8. カントロヴィッチの不等式 $A$ を正定値エルミート行列とし、その固有値がすべて閉区間 $[m,M]$（$0 $\|\boldsymbol{x}\|=1$ として示せば足ります。$A$ の正規直交固有ベクトル $\boldsymbol{u}_k$（固有値 $\lambda_k\in[m,M]$）で $\boldsymbol{x}=\sum c_k\boldsymbol{u}_k$、$p_k=|c_k|^2$（$\sum p_k=1$）とすると、$\langle A\boldsymbol{x},\boldsymbol{x}\rangle=\sum p_k\lambda_k=:a$、$\langle A^{-1}\boldsymbol{x},\boldsymbol{x}\rangle=\sum p_k/\lambda_k=:b$ です。$\lambda\in[m,M]$ なら $(\lambda-m)(\lambda-M)\le0$ から $\lambda^2-(M+m)\lambda+Mm\le0$、$\lambda>0$ で割って $\lambda+\dfrac{Mm}{\lambda}\le M+m$ です。$\lambda=\lambda_k$ について $p_k$ 倍して足すと $a+Mm\,b\le M+m$ です。相加平均・相乗平均の不等式（$a,b>0$）で $2\sqrt{Mm\,ab}\le a+Mm\,b\le M+m$ なので、平方して $ab\le\dfrac{(M+m)^2}{4Mm}$ です。 [条件数](/mathematics/numerical-linear-algebra/norms-and-conditioning#def-condition-number) $\kappa=M/m$ を使うと右辺は $\frac{(\kappa+1)^2}{4\kappa}$ です。この不等式は、正定値対称行列の 2 次関数の最急降下法（厳密な直線探索つき）の収束率の評価、$f(\boldsymbol{x}_{k+1})-f^{\ast}\le\bigl(\frac{\kappa-1}{\kappa+1}\bigr)^2\bigl(f(\boldsymbol{x}_k)-f^{\ast}\bigr)$ を導く鍵になります（この記事では導出しません）。$A=I$（$\kappa=1$）では等号です。 ## 9. 演習 $A=\begin{pmatrix}6&2&1\\2&5&1\\1&1&4\end{pmatrix}$ が正定値であることを、ゲルシュゴリンの円板定理で示せ。さらに最大固有値の上界を与えよ。 $A$ は対称で、$R_1=3$、$R_2=3$、$R_3=2$ です。円板は $D_1=[3,9]$、$D_2=[2,8]$、$D_3=[2,6]$ で、固有値（実数）はすべて $[2,9]$ に入ります（）。最小固有値は $2$ 以上で正なので $A$ は正定値、最大固有値は $9$ 以下です（実際、$\lambda_1\le\max_i\sum_j|a_{ij}|=\max(9,8,6)=9$ とも一致します）。 $A=\begin{pmatrix}2&1&0\\1&2&1\\0&1&2\end{pmatrix}$ と $E=\begin{pmatrix}0&0&0\\0&0&0\\0&0&\varepsilon\end{pmatrix}$（$\varepsilon>0$）について、ワイルの不等式 の 1 と 2 により、$\lambda_2(A+E)$ の範囲を求めよ（$A$ の固有値は $2+\sqrt2,2,2-\sqrt2$）。 $E$ の固有値は $\lambda_1(E)=\varepsilon$、$\lambda_2(E)=\lambda_3(E)=0$ です。 の 1 で $i=2$、$j=1$: $\lambda_2(A+E)\le\lambda_2(A)+\lambda_1(E)=2+\varepsilon$。2 で $i=2$、$j=3$（$i+j-n=2$）: $\lambda_2(A+E)\ge\lambda_2(A)+\lambda_3(E)=2$。よって $2\le\lambda_2(A+E)\le2+\varepsilon$ です。$\|E\|=\varepsilon$ なので、 の 3 から得られる評価 $|\lambda_2(A+E)-2|\le\varepsilon$ と同じ幅ですが、$E$ が半正定値であることを使ったので下側が $2$ で押さえられている点が精密です。 $A=\begin{pmatrix}1&0\\1&1\\0&1\end{pmatrix}$ について、$\det(A^{\mathsf{T}}A)$ を直接計算し、コーシー・ビネの公式（）の小行列式の和と比較せよ。またアダマールの不等式との関係を述べよ。 $A^{\mathsf{T}}A=\begin{pmatrix}2&1\\1&2\end{pmatrix}$ で $\det=3$。$A$ の $2\times2$ 小行列式（行 $\{1,2\},\{1,3\},\{2,3\}$）は $1\cdot1-0\cdot1=1$、$1\cdot1-0\cdot0=1$、$1\cdot1-1\cdot0=1$ で、平方の和は $3$ と一致します。アダマールの不等式（ の 1）では $\det(A^{\mathsf{T}}A)\le2\cdot2=4$ で、$3\le4$ です（列 $(1,1,0)$ と $(0,1,1)$ は直交しておらず、等号になりません）。 $A=\begin{pmatrix}0&2\\1&1\end{pmatrix}$ は既約である。$\boldsymbol{x}=(2,3)^{\mathsf{T}}$ でコリャッツ・ヴィーラントの評価を行い、$\rho(A)$ を挟め。また $\boldsymbol{x}=(1,1)^{\mathsf{T}}$ ではどうか。 $\boldsymbol{x}=(2,3)^{\mathsf{T}}$ で $A\boldsymbol{x}=(6,5)^{\mathsf{T}}$、比は $3$ と $5/3$ なので $5/3\le\rho\le3$ です（ の 2）。$\boldsymbol{x}=(1,1)^{\mathsf{T}}$ では $A\boldsymbol{x}=(2,2)^{\mathsf{T}}$、比はどちらも $2$ なので、上下の評価が一致して $\rho(A)=2$ です（特性多項式 $t^2-t-2=(t-2)(t+1)$ と整合し、$(1,1)$ は固有ベクトル）。 ## 参考文献 - R. A. Horn and C. R. Johnson, *Matrix Analysis*, 2nd ed., Cambridge University Press, 2013 — ゲルシュゴリンの定理、非負行列（ペロン・フロベニウス）、行列方程式を扱う各章。 - R. Bhatia, *Matrix Analysis*, Graduate Texts in Mathematics 169, Springer, 1997 — 第 III 章（固有値の変分原理、ワイルの不等式、交互定理）。 - A. Berman and R. J. Plemmons, *Nonnegative Matrices in the Mathematical Sciences*, SIAM, 1994 — 非負行列の理論。 - F. Zhang (ed.), *The Schur Complement and Its Applications*, Springer, 2005 — シュール補行列、ウッドベリーの公式。 - S. Gerschgorin, "Über die Abgrenzung der Eigenwerte einer Matrix", *Izvestiya Akademii Nauk SSSR, Seriya Matematicheskaya* 6 (1931), 749–754.