譜半徑

數學上，矩陣或有界線性算子的譜半徑（spectral radius）是其特徵值絕對值中的最大值（也就是矩陣的譜中元素絕對值中的最小上界），會表示為ρ(·)。

## 矩陣

令λ1, ..., λn是矩陣A ∈ Cn×n中的特徵值，則其譜半徑 ρ(A) is 定義為：

$$\rho(A) = \max \left \{ |\lambda_1|, \dotsc, |\lambda_n| \right \}.$$

$A$的條件數可以用譜半徑表示，公式為$\rho(A) \rho(A^{-1})$。

譜半徑是矩陣所有範數的一種下確界（infimum）。另一方面，$\rho(A) \leqslant \|A\|$對每一個矩陣範數 $\|\cdot\|$都成立，Gelfand公式指出$\rho(A) = \lim_{k\to\infty} \|A^k\|^{1/k}$。不過，針對任意向量$\mathbf{v} \in \mathbb{C}^n$，譜半徑不一定會滿足$\|A\mathbf{v}\| \leqslant \rho(A) \|\mathbf{v}\|$。若要說明原因，可以令$r>1$為任意數，考慮矩陣$C_r = \begin{pmatrix} 0 & r^{-1} \\ r & 0 \end{pmatrix}$。$C_r$的特徵多項式是$\lambda^2 - 1$，因此其特徵值為$\pm 1$，且$\rho(C_r) = 1$。不過$C_r \mathbf{e}_1 = r \mathbf{e}_2$，因此$\| C_r \mathbf{e}_1 \| = r > 1 = \rho(C_r) \|\mathbf{e}_1\|$，其中$\|\cdot\|$是$\mathbb{C}^n$上的任何$\ell^p$範數。至於可以當$k \to \infty$時，讓$\|C_r^k\|^{1/k} \to 1$的原因是$C_r^2 = I$，因此當$k \to \infty$時，使$\|C_r^k\|^{1/k} \leqslant \|C_r\|^{1/k} = r^{1/k} \to 1$。

$\|A\mathbf{v}\| \leqslant \rho(A) \|\mathbf{v}\|$針對所有$\mathbf{v} \in \mathbb{C}^n$

成立的條件是$A$為埃爾米特矩陣及$\|\cdot\|$為歐幾里得範數。

## 圖

有限圖的譜半徑定義為其鄰接矩陣的譜半徑。

此一定義可以擴散到無限圖，但是其每個頂點都只連接有限個頂點（存在一實數C使得每一個頂點的度都小於C）。此情形下，針對圖G可定義：

$$\ell^2(G) = \left \{ f : V(G) \to \mathbf{R} \ : \ \sum\nolimits_{v \in V(G)} \left \|f(v)^2 \right \| < \infty \right \}.$$

令γ是 G的鄰接算子：

$$\begin{cases} \gamma : \ell^2(G) \to \ell^2(G) \\ (\gamma f)(v) = \sum_{(u,v) \in E(G)} f(u) \end{cases}$$

G的譜半徑定義為有界線性算子γ的譜半徑。

## 上界

### 矩陣譜半徑的上界

以下的命題指出了一個簡單但是有用的矩陣譜半徑上界：

命題：令A ∈ Cn×n，其譜半徑為ρ(A)，以及相容（Consistent）矩陣範數。則針對每一個整數$k \geqslant 1$:

$$\rho(A)\leq \|A^k\|^{\frac{1}{k}}.$$

證明

令(v, λ)為矩陣A的特徵值-特徵向量對。利用矩陣範數的次可乘性（sub-multiplicative property），可得：

$$|\lambda|^k\|\mathbf{v}\| = \|\lambda^k \mathbf{v}\| = \|A^k \mathbf{v}\| \leq \|A^k\|\cdot\|\mathbf{v}\|$$

因為v ≠ 0，可得

$$|\lambda|^k \leq \|A^k\|$$

因此

$$\rho(A)\leq \|A^k\|^{\frac{1}{k}}.$$

### 圖譜半徑的上界

有關n個頂點，m個邊的圖，有許多的譜半徑的上界公式。例如，若

$$\frac{(k-2)(k-3)}{2} \leq m-n \leq \frac{k(k-3)}{2}$$

其中$3 \le k \le n$為整數，則：

$$\rho(G) \leq \sqrt{2 m-n-k+\frac{5}{2}+\sqrt{2 m-2 n+\frac{9}{4}}}$$

## 乘冪數列

### 定理

譜半徑和矩陣乘冪數列是否收斂有緊密的關係。以下的定理會成立：

定理：令A ∈ Cn×n，其譜半徑ρ(A)。則ρ(A) < 1若且唯若

$$\lim_{k \to \infty} A^k = 0.$$

另一方面，若ρ(A) > 1, $\lim_{k \to \infty} \|A^k\| = \infty$。上述敘述針對Cn×n上的任何矩陣範數都有效。

### 定理證明

假設問題中的極限值為零，可以證明ρ(A) < 1。令(v, λ)為A的特徵值和特徵向量對。因為Akv = λkv可得：

$$\begin{align} <pre> 0 &= \left(\lim_{k \to \infty} A^k \right) \mathbf{v} \\ &= \lim_{k \to \infty} \left(A^k\mathbf{v} \right ) \\ &= \lim_{k \to \infty} \lambda^k\mathbf{v} \\ &= \mathbf{v} \lim_{k \to \infty} \lambda^k </pre> \end{align}$$

因為假設v ≠ 0，會得到

$$\lim_{k \to \infty}\lambda^k = 0$$

表示|λ| < 1。因為這對任何一個特徵值都會成立，因此可知ρ(A) < 1。

接下來假設A的譜半徑小於1。根據若爾當標準型定理，可以知道針對所有的A ∈ Cn×n，存在V, J ∈ Cn×n以及非奇異的V和J分塊對角矩陣使得：

$$A = VJV^{-1}$$

而

$$J=\begin{bmatrix} J_{m_1}(\lambda_1) & 0 & 0 & \cdots & 0 \\ 0 & J_{m_2}(\lambda_2) & 0 & \cdots & 0 \\ \vdots & \cdots & \ddots & \cdots & \vdots \\ 0 & \cdots & 0 & J_{m_{s-1}}(\lambda_{s-1}) & 0 \\ 0 & \cdots & \cdots & 0 & J_{m_s}(\lambda_s) \end{bmatrix}$$

其中

$$J_{m_i}(\lambda_i)=\begin{bmatrix} \lambda_i & 1 & 0 & \cdots & 0 \\ 0 & \lambda_i & 1 & \cdots & 0 \\ \vdots & \vdots & \ddots & \ddots & \vdots \\ 0 & 0 & \cdots & \lambda_i & 1 \\ 0 & 0 & \cdots & 0 & \lambda_i \end{bmatrix}\in \mathbf{C}^{m_i \times m_i}, 1\leq i\leq s.$$

因此可得

$$A^k=VJ^kV^{-1}$$

因為J是分塊對角矩陣

$$J^k=\begin{bmatrix} J_{m_1}^k(\lambda_1) & 0 & 0 & \cdots & 0 \\ 0 & J_{m_2}^k(\lambda_2) & 0 & \cdots & 0 \\ \vdots & \cdots & \ddots & \cdots & \vdots \\ 0 & \cdots & 0 & J_{m_{s-1}}^k(\lambda_{s-1}) & 0 \\ 0 & \cdots & \cdots & 0 & J_{m_s}^k(\lambda_s) \end{bmatrix}$$

而$m_i \times m_i$若爾當方塊矩陣k次方可以得到，針對$k \geq m_i-1$：

$$J_{m_i}^k(\lambda_i)=\begin{bmatrix} \lambda_i^k & {k \choose 1}\lambda_i^{k-1} & {k \choose 2}\lambda_i^{k-2} & \cdots & {k \choose m_i-1}\lambda_i^{k-m_i+1} \\ 0 & \lambda_i^k & {k \choose 1}\lambda_i^{k-1} & \cdots & {k \choose m_i-2}\lambda_i^{k-m_i+2} \\ \vdots & \vdots & \ddots & \ddots & \vdots \\ 0 & 0 & \cdots & \lambda_i^k & {k \choose 1}\lambda_i^{k-1} \\ 0 & 0 & \cdots & 0 & \lambda_i^k \end{bmatrix}$$

因此，若$\rho(A) < 1$，則針對所有的i，$|\lambda_i| < 1$都會成立。因此針對所有的i，可得：

$$\lim_{k \to \infty}J_{m_i}^k=0$$

這也表示

$$\lim_{k \to \infty} J^k = 0.$$

因此

$$\lim_{k \to \infty}A^k=\lim_{k \to \infty}VJ^kV^{-1}=V \left (\lim_{k \to \infty}J^k \right )V^{-1}=0$$

另一方面，若$\rho(A)>1$，當k增加時，在J中至少有一個元素無法維持有界，因此證明了定理的第二部份。

## Gelfand公式

### 定理

以下的定理可以用矩陣範數的極限來計算T譜半徑

定理（Gelfand公式，1941年）：令任何矩陣範數，可得

$\rho(A)=\lim_{k \to \infty} \left \|A^k \right \|^{\frac{1}{k}}.$.

### 證明

令任意ε > 0，先建構以下二個矩陣：

$$A_{\pm}= \frac{1}{\rho(A) \pm\varepsilon}A.$$

則:

$$\rho \left (A_{\pm} \right ) = \frac{\rho(A)}{\rho(A) \pm \varepsilon}, \qquad \rho (A_+) < 1 < \rho (A_-).$$

先將之前的定理應用到A+:

$$\lim_{k \to \infty} A_+^k=0.$$

這表示，根據級數極限定理，一定存在N+ ∈ N使得針對所有的k ≥ N+，下式都成立

$$\begin{align} \left \|A_+^k \right \| < 1 \end{align}$$

因此

$$\begin{align} \left \|A^k \right \|^{\frac{1}{k}} < \rho(A)+\varepsilon. \end{align}$$

將之前的定理用在A−，表示$\|A_-^k\|$無界，一定存在N− ∈ N使得針對所有的k ≥ N−，下式都成立

$$\begin{align} \left\|A_-^k \right \| > 1 \end{align}$$

因此

$$\begin{align} \left\|A^k \right\|^{\frac{1}{k}} > \rho(A)-\varepsilon. \end{align}$$

令，可得：

$$\forall \varepsilon>0\quad \exists N\in\mathbf{N} \quad \forall k\geq N \quad \rho(A)-\varepsilon < \left \|A^k \right \|^{\frac{1}{k}} < \rho(A)+\varepsilon$$

因此，依定義，可得下式

$$\lim_{k \to \infty} \left \|A^k \right \|^{\frac{1}{k}} = \rho(A).$$

## 舉例

考慮以下矩陣

$$A=\begin{bmatrix} 9 & -1 & 2\\ -2 & 8 & 4\\ 1 & 1 & 8 \end{bmatrix}$$

其中的特徵值為5, 10, 10。依照定義，ρ(A) = 10。在以下的表中，會以四個最常用的矩陣範式，在k增加時，計算$\|A^k\|^{\frac{1}{k}}$（注意，因為此矩陣特殊的形式，$\|.\|_1=\|.\|_\infty$）:

k｜$\|.\|_1=\|.\|_\infty$｜$\|.\|_F$｜$\|.\|_2$

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

$\vdots$｜$\vdots$｜$\vdots$｜$\vdots$

10｜10.604944422｜10.455910430｜10.183690042

11｜10.548677680｜10.413702213｜10.166990229

12｜10.501921835｜10.378620930｜10.153031596

$\vdots$｜$\vdots$｜$\vdots$｜$\vdots$

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

$\vdots$｜$\vdots$｜$\vdots$｜$\vdots$

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

$\vdots$｜$\vdots$｜$\vdots$｜$\vdots$

1000｜10.005879594｜10.004460985｜10.001823382

2000｜10.002939365｜10.002230244｜10.000911649

3000｜10.001959481｜10.001486774｜10.000607757

$\vdots$｜$\vdots$｜$\vdots$｜$\vdots$

10000｜10.000587804｜10.000446009｜10.000182323

20000｜10.000293898｜10.000223002｜10.000091161

30000｜10.000195931｜10.000148667｜10.000060774

$\vdots$｜$\vdots$｜$\vdots$｜$\vdots$

100000｜10.000058779｜10.000044600｜10.000018232

## 有界線性算子

針對有界線性算子 A 及算子範數

·

可以得到

$$\rho(A) = \lim_{k \to \infty}\|A^k\|^{\frac{1}{k}}.$$

（複數希爾伯特空間上的）有界算子若其譜半徑等於數值半徑，可以稱為「譜算子」（spectraloid operator）。其中一個例子是正規算子。