数学におけるスペクトル半径(スペクトルはんけい、英: spectral radius)とは、複素正方行列や線形位相空間上の有界線形作用素の固有値の絶対値の最小上界のことである。ギリシャ文字 ρ によって表記されることが多い。
行列のスペクトル半径および諸性質
複素正方行列
について、その固有値を
(実数または複素数)とする。このときの
のスペクトル半径
は以下のように定義される。
![{\displaystyle \rho ({\boldsymbol {A}}):=\max _{i}(|\lambda _{i}|)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4a5d277af6097c8cb1c5ee673c0d45edfb64cb83)
より一般に、単位的バナッハ環の元
について、そのスペクトル σ(A) = {λ ∈ C | λI - A は可逆でない } に含まれる数の絶対値の上限
が
のスペクトル半径と呼ばれる(ここで I はバナッハ環の単位元とする)。 有界線形作用素
と作用素ノルム ||·|| に対し、次式がなりたつ(#ノルムによる評価節を参照のこと)。
![{\displaystyle \rho ({\boldsymbol {A}})=\lim _{k\to \infty }\|{\boldsymbol {A}}^{k}\|^{1/k}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3ae2187002c33a29d6cccdb0b5167dd4f974d5b2)
複素ヒルベルト空間上の有界作用素は、そのスペクトル半径が数域半径と一致する場合、spectraloid operator と呼ばれる。このような作用素の例としては、正規作用素がある。
等比列の収束
スペクトル半径は行列の等比列の収束性と次のようにして密接に関係している。
を複素行列、
をそのスペクトル半径とすると、
のとき、およびそのときに限り
である。
これは特に、任意の行列ノルム ||・|| について
- ρ(A) < 1 ならば ||A|| → 0(ノルムの連続性により)
- ρ(A) > 1 ならば ||A|| → ∞(ノルムの同値性により)
ということを導く。
が ρ(A) < 1 を導くことは以下のようにしてわかる。 (v, λ) を行列 A の固有ベクトル-固有値の組とすると、
![{\displaystyle {\boldsymbol {A}}^{k}{\boldsymbol {v}}=\lambda ^{k}{\boldsymbol {v}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1749427e0c4236c65d604e0111c617714f6cb9b1)
であるから、
![{\displaystyle 0=(\lim _{k\to \infty }{\boldsymbol {A}}^{k}){\boldsymbol {v}}=\lim _{k\to \infty }{\boldsymbol {A}}^{k}{\boldsymbol {v}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2d00bba351c66a6ccc24e82b57f14d41fcc9f578)
ここで、v ≠ 0 であることより、
![{\displaystyle \lim _{k\to \infty }\lambda ^{k}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/45f4737182faa5a052d3e5e094f2542c7d9e8308)
でなければならないが、これは、|λ| < 1 であることを意味する。これがすべての固有値 λ に対して成立しなければならないから、ρ(A) < 1 と結論づけることができる。
一方、ρ(A) < 1 が
を導くことは以下のようにしてわかる。ジョルダン標準形の理論から、任意の複素行列 A ∈ MnC について、互いに可換な半単純行列 S とベキ零行列 N があって、A = S + N、ρ(A) = ρ(S) が成立している。 K をNK = 0 であるような自然数とすれば、任意の自然数 k について
![{\displaystyle {\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}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/afb258fe212fb9c0e5bea5ccc8a9fa35a654f8a8)
が成り立っている。ρ(S)が1より小さいため任意の j について
![{\displaystyle {\binom {k}{j}}{\boldsymbol {S}}^{k-j}\rightarrow 0\quad (k\rightarrow \infty )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c7bc228e0b15fe5a39fee9aea90af302b8ba7f9b)
であり、したがって上式右辺の有限和の各項は 0 に収束している。
ノルムによる評価
複素行列のスペクトル半径と任意の行列ノルム ||·|| に関して、次式が成立する (Gelfand, 1941)。
![{\displaystyle \rho ({\boldsymbol {A}})=\lim _{k\to \infty }\|{\boldsymbol {A}}^{k}\|^{1/k}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3ae2187002c33a29d6cccdb0b5167dd4f974d5b2)
この定理は以下のようにして示される。 ε > 0 を任意の正の実数とする。このとき、
![{\displaystyle {\tilde {\boldsymbol {A}}}=(\rho ({\boldsymbol {A}})+\epsilon )^{-1}{\boldsymbol {A}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/581a652e558ccf70ec9249eac499021d9c238588)
について
![{\displaystyle \rho ({\tilde {\boldsymbol {A}}})={\frac {\rho ({\boldsymbol {A}})}{\rho ({\boldsymbol {A}})+\epsilon }}<1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/26cb59c85d0615bcb89e639d1001e397153b068e)
だから、#等比列の収束により、
![{\displaystyle \lim _{k\to \infty }{\tilde {\boldsymbol {A}}}^{k}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3b5809642a46b706e8da21e888d76286c5b79b55)
が成り立っている。したがって、ある自然数 N1 ∈ N が存在して、
![{\displaystyle \forall k\geq N_{1}\Rightarrow \|{\tilde {\boldsymbol {A}}}^{k}\|<1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/71e8e0be6b213a93520a71e7741d25dba58d0bef)
が成り立つ。これは
![{\displaystyle \forall k\geq N_{1}\Rightarrow \|{\boldsymbol {A}}^{k}\|^{1/k}<\rho ({\boldsymbol {A}})+\epsilon .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cc31c40ba2e6b8de9697ca8cdb14c6600f985389)
ということを示している。同様にして
![{\displaystyle {\check {\boldsymbol {A}}}={\frac {\boldsymbol {A}}{\rho ({\boldsymbol {A}})-\epsilon }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/29a84187a384e08307e72c07733752d064f1b46a)
を考えることにより、ある自然数 N1 ∈ N が存在して、
![{\displaystyle \forall k\geq N_{1}\Rightarrow \|{\check {\boldsymbol {A}}}^{k}\|>1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4d6259b1fbb8c041e2c2847ae9ba2be50aa13aaf)
がわかる。以上のことから
![{\displaystyle \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 }](https://wikimedia.org/api/rest_v1/media/math/render/svg/10d7cde88311ccb10c57bb85b595db8c8aa3162a)
が言えるが、これは
![{\displaystyle \lim _{k\to \infty }\|{\boldsymbol {A}}^{k}\|^{1/k}=\rho ({\boldsymbol {A}})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3b990fb85a37454a1fb65768a2a6eea73c3bc516)
ということを表している。
さらに||·|| が 一貫性を持つ場合には、任意の複素行列 A ∈ MnC と k ∈ N に対し
![{\displaystyle \rho (A)\leq \|{\boldsymbol {A}}^{k}\|^{1/k}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/50f0c5c70fa80319f043326dee3159b5189fb575)
が成立している。これは以下のようにして示すことができる。A の固有ベクトル v と対応する固有値 λ について、行列ノルムの一貫性から次式を得る。
![{\displaystyle |\lambda |^{k}\|{\boldsymbol {v}}\|=\|\lambda ^{k}{\boldsymbol {v}}\|=\|{\boldsymbol {A}}^{k}{\boldsymbol {v}}\|\leq \|{\boldsymbol {A}}^{k}\|\cdot \|{\boldsymbol {v}}\|}](https://wikimedia.org/api/rest_v1/media/math/render/svg/48c2fbbebeee55e022c073e5520e3d9558be5213)
ここで、v ≠ 0 であるので、任意の固有値 λ に対して次式を得る。
![{\displaystyle |\lambda |^{k}\leq \|{\boldsymbol {A}}^{k}\|}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5bed8402145c388a623ee44c65360b57435c514a)
したがって、
![{\displaystyle \rho ({\boldsymbol {A}})\leq \|{\boldsymbol {A}}^{k}\|^{1/k}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7784d7173d1f11a892703d72921efc08dfd77c0b)
が成立する。また、ヒルベルト空間上の作用素ノルムについては
![{\displaystyle \rho ({\boldsymbol {A}}^{*}{\boldsymbol {A}})=\|{\boldsymbol {A}}^{*}{\boldsymbol {A}}\|=\|{\boldsymbol {A}}\|^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c22ce71a414c602cce31ffd1aa948ffd8d9422e2)
が成り立つ。
Gelfand の公式は、有限個の行列の積のスペクトル半径に対しても考えることができる。すべての行列が可換であると仮定すると、次式を得る。
![{\displaystyle \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}).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/02e68b313e3de396122570cd2d12c62a78f3bc71)
例
例: 次の行列を考える。
![{\displaystyle {\boldsymbol {A}}={\begin{bmatrix}9&-1&2\\-2&8&4\\1&1&8\end{bmatrix}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2f49ebc7809c3a0a86aa0404f32d9db22890171f)
この固有値は 5, 10, 10 であるから、定義より、スペクトル半径は ρ(A)=10 である。以下の表には、ベクトルの p-ノルムから誘導された行列の作用素ノルムおよびヒルベルト-シュミットノルム(フロベニウスノルム)に関する
の、k の増加に対する値が列挙されている(この行列の場合には
となっている)。
k | ![{\displaystyle \|.\|_{1}=\|.\|_{\infty }}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8646cb1a450de70416bf488892a8a046c28eeae5) | ![{\displaystyle \|.\|_{F}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d791c86c0c00587d1109f649a049defff37ee304) | |
1 | 14 | 15.362291496 | 10.681145748 |
2 | 12.649110641 | 12.328294348 | 10.595665162 |
3 | 11.934831919 | 11.532450664 | 10.500980846 |
4 | 11.501633169 | 11.151002986 | 10.418165779 |
5 | 11.216043151 | 10.921242235 | 10.351918183 |
| | | |
10 | 10.604944422 | 10.455910430 | 10.183690042 |
11 | 10.548677680 | 10.413702213 | 10.166990229 |
12 | 10.501921835 | 10.378620930 | 10.153031596 |
| | | |
20 | 10.298254399 | 10.225504447 | 10.091577411 |
30 | 10.197860892 | 10.149776921 | 10.060958900 |
40 | 10.148031640 | 10.112123681 | 10.045684426 |
50 | 10.118251035 | 10.089598820 | 10.036530875 |
| | | |
100 | 10.058951752 | 10.044699508 | 10.018248786 |
200 | 10.029432562 | 10.022324834 | 10.009120234 |
300 | 10.019612095 | 10.014877690 | 10.006079232 |
400 | 10.014705469 | 10.011156194 | 10.004559078 |
| | | |
1000 | 10.005879594 | 10.004460985 | 10.001823382 |
2000 | 10.002939365 | 10.002230244 | 10.000911649 |
3000 | 10.001959481 | 10.001486774 | 10.000607757 |
| | | |
10000 | 10.000587804 | 10.000446009 | 10.000182323 |
20000 | 10.000293898 | 10.000223002 | 10.000091161 |
30000 | 10.000195931 | 10.000148667 | 10.000060774 |
| | | |
100000 | 10.000058779 | 10.000044600 | 10.000018232 |
グラフのスペクトル半径
有限グラフのスペクトル半径は、その隣接行列のスペクトル半径として定義される。
この定義は、頂点の次数が有界な無限グラフ(すなわち、ある実数 C が存在して、グラフ中のすべての頂点の次数が C より小さくなる)の場合に拡張される。この場合、グラフ G に対して、その頂点集合を基底にするようなヒルベルト空間 l2(V(G)) 上に
![{\displaystyle (\gamma f)(v)=\sum _{(u,v)\in E(G)}f(u)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/286933b8882c4fa4a084ef90932c99dae8b13ff9)
によって G の隣接作用素とよばれる l2(V(G)) 上の有界作用素 γ を考えることができる。このとき、 γ のスペクトル半径のことを G のスペクトル半径という。
関連記事
参考文献
- Gert K. Pedersen (2001). Analysis Now. Graduate Texts in Mathematics (Corrected ed. ed.). Springer. ISBN 978-0387967882
- Ronald G. Douglas (1998). Banach Algebra Techniques in Operator Theory. Graduate Texts in Mathematics. Springer. ISBN 978-0387983776