スペクトル半径

数学における**スペクトル半径**（スペクトルはんけい、）とは、複素正方行列や線形位相空間上の有界線形作用素の固有値の絶対値の最小上界のことである。ギリシャ文字 *[[ρ]]* によって表記されることが多い。

## 行列のスペクトル半径および諸性質
複素正方行列\boldsymbol{A}\in{\mathbb C}^{n\times n}について、その固有値を\lambda_1, \lambda_2, \ldots, \lambda_n
（実数または複素数）とする。このときの\boldsymbol{A}
の**スペクトル半径 <math>\rho(\boldsymbol{A})
</math>**は以下のように定義される。

:\rho(\boldsymbol{A}) := \max_i(|\lambda_i|)

より一般に、単位的バナッハ環の元 *<math>\boldsymbol{A}
</math>*について、その**スペクトル** *&sigma;*(***A***) = {*&lambda;* &isin; **C** | *&lambda;*'I*** - ***A*** は可逆でない } に含まれる数の絶対値の上限 '''<math>\rho(\boldsymbol{A})
</math>'''が ''<math>\boldsymbol{A}
</math>'' のスペクトル半径と呼ばれる(ここで ***I*** はバナッハ環の[[単位元]]とする)。
有界線形作用素 ''<math>\boldsymbol{A}
</math>'' と[[作用素ノルム]] ||·|| に対し、次式がなりたつ([[#ノルムによる評価]]節を参照のこと)。

:<math>\rho(\boldsymbol{A}) = \lim_{k \to \infty}\|\boldsymbol{A}^k\|^{1/k}.</math>

複素ヒルベルト空間上の有界作用素は、そのスペクトル半径が[[数域半径]]と一致する場合、'''spectraloid operator''' と呼ばれる。このような作用素の例としては、[[正規作用素]]がある。

== 等比列の収束 ==
スペクトル半径は行列の等比列の収束性と次のようにして密接に関係している。 <math>\boldsymbol{A}\in{\mathbb C}^{n\times n}</math>を複素行列、'''<math>\rho(\boldsymbol{A})
</math>'''をそのスペクトル半径とすると、

:<math>\lim_{k \to \infty} \boldsymbol{A}^k=0</math> のとき、およびそのときに限り <math>\rho(\boldsymbol{A})<1</math>である。

これは特に、任意の[[行列ノルム]] ||・|| について
*''&rho;''(***A***) < 1 ならば ||***A***|| &rarr; 0（ノルムの連続性により）
*''&rho;''(***A***) > 1 ならば ||***A***|| &rarr; &infin;（ノルムの同値性により）
ということを導く。

<math>\lim_{k \to \infty} \boldsymbol{A}^k = 0</math>が ''&rho;''(***A***) < 1 を導くことは以下のようにしてわかる。 (***v***, ''λ'') を行列 ***A*** の[[固有ベクトル]]-[[固有値]]の組とすると、

:<math>\boldsymbol{A}^k\boldsymbol{v} = \lambda^k\boldsymbol{v}</math>

であるから、
:<math>0 = (\lim_{k \to \infty} \boldsymbol{A}^k)\boldsymbol{v}= \lim_{k \to \infty} \boldsymbol{A}^k\boldsymbol{v}</math>

ここで、***v*** ≠ 0 であることより、

:<math>\lim_{k \to \infty}\lambda^k = 0</math>

でなければならないが、これは、|''λ''| < 1 であることを意味する。これがすべての固有値 ''λ'' に対して成立しなければならないから、''ρ''(***A***) < 1 と結論づけることができる。

一方、''&rho;''(***A***) < 1 が<math>\lim_{k \to \infty} \boldsymbol{A}^k = 0</math>を導くことは以下のようにしてわかる。[[ジョルダン標準形]]の理論から、任意の複素行列 ***A*** ∈ M<sub>''n''</sub>'''C''' について、互いに可換な半単純行列 ***S*** とベキ零行列 ***N*** があって、***A*** = ***S*** + ***N***、''&rho;''(***A***) = ''&rho;''(***S***) が成立している。 ''K'' を***N***<sup>''K''</sup> = 0 であるような自然数とすれば、任意の自然数 ''k'' について
:<math>\boldsymbol{A}^k = \boldsymbol{S}^k + k \boldsymbol{S}^{k-1} \boldsymbol{N} + \binom{k}{2} \boldsymbol{S}^{k-2} \boldsymbol{N}^2 + \dotsb + \binom{k}{K-1} \boldsymbol{S}^{k-K+1} \boldsymbol{N}^{K-1}</math>
が成り立っている。''&rho;''(***S***)が1より小さいため任意の ''j'' について
:<math> \binom{k}{j} \boldsymbol{S}^{k-j} \rightarrow 0 \quad (k \rightarrow \infty) </math>
であり、したがって上式右辺の有限和の各項は 0 に収束している。

=== ノルムによる評価 ===
複素行列のスペクトル半径と任意の行列ノルム ||·|| に関して、次式が成立する (Gelfand, 1941)。

:<math>\rho(\boldsymbol{A})=\lim_{k \to \infty}\|\boldsymbol{A}^k\|^{1/k}.</math>

この定理は以下のようにして示される。 ''ε'' > 0 を任意の正の実数とする。このとき、

:<math>\tilde{\boldsymbol{A}}=(\rho(\boldsymbol{A}) + \epsilon)^{-1} \boldsymbol{A}.</math>

について

:<math>\rho(\tilde{\boldsymbol{A}}) = \frac{\rho(\boldsymbol{A})}{\rho(\boldsymbol{A})+\epsilon} < 1</math>

だから、[[#等比列の収束]]により、

:<math>\lim_{k \to \infty}\tilde{\boldsymbol{A}}^k=0</math>

が成り立っている。したがって、ある自然数 ''N''<sub>1</sub> ∈ '''N''' が存在して、

:<math>\forall k\geq N_1 \Rightarrow \|\tilde{\boldsymbol{A}}^k\| < 1</math>

が成り立つ。これは

:<math>\forall k\geq N_1 \Rightarrow \|\boldsymbol{A}^k\|^{1/k} < \rho(\boldsymbol{A})+\epsilon.</math>

ということを示している。同様にして

:<math> \check{\boldsymbol{A}} = \frac{\boldsymbol{A}}{\rho(\boldsymbol{A}) - \epsilon}</math>

を考えることにより、ある自然数 ''N''<sub>1</sub> ∈ '''N''' が存在して、

:<math>\forall k\geq N_1 \Rightarrow \|\check{\boldsymbol{A}}^k\| > 1</math>

がわかる。以上のことから

:<math>\forall \epsilon>0, \exists N\in\mathbb{N}: \forall k\geq N \Rightarrow \rho(\boldsymbol{A}) - \epsilon < \|\boldsymbol{A}^k\|^{1/k} < \rho(\boldsymbol{A})+\epsilon</math>

が言えるが、これは

:<math>\lim_{k \to \infty}\|\boldsymbol{A}^k\|^{1/k} = \rho(\boldsymbol{A})</math>
ということを表している。

さらに||·|| が [[行列ノルム|一貫性]]を持つ場合には、任意の複素行列 ***A*** ∈ M<sub>''n''</sub>'''C''' と ''k'' ∈ '''N''' に対し

:<math>\rho(A)\leq \|\boldsymbol{A}^k\|^{1/k}</math>

が成立している。これは以下のようにして示すことができる。***A*** の[[固有ベクトル]] ***v*** と対応する[[固有値]] ''λ'' について、行列ノルムの一貫性から次式を得る。

:<math>|\lambda|^k\|\boldsymbol{v}\| = \|\lambda^k \boldsymbol{v}\| = \|\boldsymbol{A}^k \boldsymbol{v}\| \leq \|\boldsymbol{A}^k\|\cdot\|\boldsymbol{v}\|</math>

ここで、'''v''' ≠ 0 であるので、任意の固有値 λ に対して次式を得る。

:<math>|\lambda|^k\leq \|\boldsymbol{A}^k\|</math>

したがって、

:<math>\rho(\boldsymbol{A})\leq \|\boldsymbol{A}^k\|^{1/k}</math>

が成立する。また、ヒルベルト空間上の作用素ノルムについては

:<math>\rho(\boldsymbol{A}^* \boldsymbol{A}) = \| \boldsymbol{A}^* \boldsymbol{A} \| = \| \boldsymbol{A} \|^2</math>
が成り立つ。

Gelfand の公式は、有限個の行列の積のスペクトル半径に対しても考えることができる。すべての行列が可換であると仮定すると、次式を得る。

:<math>
\rho(\boldsymbol{A}_1 \boldsymbol{A}_2 \dotsb \boldsymbol{A}_n) \leq \rho(\boldsymbol{A}_1) \rho(\boldsymbol{A}_2) \dotsb \rho(\boldsymbol{A}_n).
</math>

== 例 ==
'''例''': 次の行列を考える。
:<math>\boldsymbol{A}=\begin{bmatrix}
9 & -1 & 2\\
-2 & 8 & 4\\
1 & 1 & 8
\end{bmatrix}</math>

この固有値は 5, 10, 10 であるから、定義より、スペクトル半径は ''ρ''(***A**'')=10 である。以下の表には、ベクトルの ''p''-ノルムから誘導された行列の作用素ノルムおよびヒルベルト-シュミットノルム(フロベニウスノルム)に関する <math>\|\boldsymbol{A}^k\|^{1/k}</math> の、''k'' の増加に対する値が列挙されている（この行列の場合には<math>\|\boldsymbol{A}^k\|_1=\|\boldsymbol{A}^k\|_\infty</math> となっている）。



==グラフのスペクトル半径==
有限[[グラフ理論|グラフ]]の**スペクトル半径*'は、その[[隣接行列]]のスペクトル半径として定義される。

この定義は、頂点の次数が有界な無限グラフ（すなわち、ある実数 *C* が存在して、グラフ中のすべての頂点の次数が *C* より小さくなる）の場合に拡張される。この場合、グラフ *G* に対して、その頂点集合を基底にするようなヒルベルト空間 l<sup>2</sup>(*V*(*G*)) 上に
:<math> (\gamma f)(v) = \sum_{(u,v) \in E(G)} f(u) </math>
によって *G* の[[隣接作用素]]とよばれる l<sup>2</sup>(*V*(*G*)) 上の有界作用素 *&gamma;* を考えることができる。このとき、 *γ* のスペクトル半径のことを *G'' のスペクトル半径という。

## 関連記事
* 作用素環論
* バナッハ空間

## 参考文献
* 
*