攝動理論

攝動理論使用一些特別的數學方法來對於很多不具精確解的問題給出近似解，這些方法從相關的較簡單問題的精確解開始入手。攝動理論將原本問題分為具有精確解的較簡單部分與不具精確解的微擾部分。攝動理論適用的問題通常具有以下性質：通過加入一個微擾項於較簡單部分的數學表述，可以計算出整個問題的近似解。

攝動理論計算出來的解答通常會表達為一個微小參數的冪級數。攝動理論解答與精確解之間的差別，可以用這微小參數來做數量比較。冪級數的第一個項目是精確解的解答。後面的項目描述解答的修正。這修正是因為精確解與原本問題的「完全解」之間的誤差而產生的。更正式地，完全解$A\,\!$的近似可以表達為一個級數：

$A=\epsilon^0 A_0 + \epsilon^1 A_1 + \epsilon^2 A_2 + \cdots\,\!$。

在這例子裏，$A_0\,\!$是簡單又有「精確解」的問題的精確解，$A_1,\, A_2, \,\!$代表由某種系統程序反覆地找到的高階項目修正。因為$\epsilon\,\!$的值很微小，這些高階項目修正應該會越來越不重要。

## 微擾階數

攝動理論的標準闡述主要是以微擾的階數來分辨：一階攝動理論或二階攝動理論。再來就是以微擾的簡併度來分辨：無簡併或有簡併。有簡併的攝動，又稱為奇異攝動（singular perturbation），比較難解，必須用到更進階的理論。

## 一階無簡併攝動理論

本段落講述微分方程的一階微擾理論。為了簡單易解，假設零微擾系統的解答是不簡併的。

### 一階本徵值修正

許多常微分方程或偏微分方程可以表達為

$Dg(x)=\lambda g(x) \,\!$；(1)

其中，$D\,\!$是某特定微分算子，$\lambda\,\!$是其本徵值。

假設微分算子可以寫為

$D=D^{(0)}+\epsilon D^{(1)} \,\!$。

其中，$\epsilon\,\!$是微小的度量。

又假設我們已知道$D^{(0)}\,\!$的解答的完備集$\{f^{(0)}_i (x)\}\,\!$；其中，解答$f^{(0)}_i(x)\,\!$是$D^{(0)}\,\!$的本徵值為$\lambda^{(0)}_i\,\!$的本徵函數。用方程表達。

$D^{(0)} f^{(0)}_i (x)=\lambda^{(0)}_i f^{(0)}_i (x) \,\!$。

還有，這一集合的解答$\{f^{(0)}_i (x)\}\,\!$形成一個正交歸一集：

$\int f^{(0)}_i (x) f^{(0)}_j (x) \,dx = \delta_{ij}\,\!$。

其中，$\delta_{ij}\,\!$是克羅內克函數。

取至零階，完全解$g(x)\,\!$應該相當接近集合里一個零微擾解。設定這零微擾解為$f^{(0)}_n (x) \,\!$。用方程表達。

$g(x)=f^{(0)}_n (x) + \mathcal{O}(\epsilon)\,\!$。

其中，$\mathcal{O}\,\!$採用大O符號來描述函數的漸近行為。

完全解的本徵值也可近似為

$\lambda=\lambda^{(0)}_n + \mathcal{O}(\epsilon)\,\!$。

將完全解$g(x)\,\!$寫為零微擾解的線性組合。

$g(x)=\sum_m c_m f^{(0)}_m (x)\,\!$；(2)

其中，除了$c_n\,\!$以外，所有的常數$c_m,\ m\ne n\,\!$的值是$\mathcal{O}(\epsilon)\,\!$；只有$c_n\,\!$的值是$\mathcal{O}(1)\,\!$。

將公式 (2)代入公式 (1)，乘以$f^{(0)}_n (x) \,\!$，利用正交歸一性，可以得到

$\lambda^{(0)}_n c_n+ \epsilon \sum_m c_m \int f^{(0)}_n(x) D^{(1)} f^{(0)}_m(x)\,dx =\lambda c_n\,\!$。

這可以很容易地改變為一個簡單的線性代數問題，一個尋找矩陣的本徵值的問題：給予$\sum_m A_{nm}c_m = \lambda c_n\!\,\!$，求$\lambda\,\!$；其中，$A_{nm}\,\!$是矩陣元素：

$A_{nm} = \lambda^{(0)}_n\delta_{nm} + \epsilon \int f^{(0)}_n(x) D^{(1)} f^{(0)}_m(x)\,dx \,\!$。

我們並不需要解析整個矩陣。注意到線性方程裡的每一個$c_m\,\!$都是$\mathcal{O}(\epsilon)\,\!$；只有$c_n\,\!$的值是$\mathcal{O}(1)\,\!$。所以，取至$\epsilon\,\!$一階，線性方程可以很容易地解析為

$\lambda=\lambda^{(0)}_n + \epsilon \int f^{(0)}_n(x) D^{(1)} f^{(0)}_n(x)\,dx \,\!$。(3)

這就是一階攝動理論的本徵值解答。一階本徵值數修正是

$\lambda^{(1)}_n=\int f^{(0)}_n(x) D^{(1)} f^{(0)}_n(x)\,dx \,\!$。

### 一階本徵函數修正

取至一階，函數$g(x)\,\!$可以用類似的推理求得。設定

$g(x)=f^{(0)}_n(x) + \epsilon f^{(1)}_n(x)\,\!$。(4)

那麼，公式 (1)變為

$\left(D^{(0)} +\epsilon D^{(1)}\right) \left( f^{(0)}_n(x) + \epsilon f^{(1)}_n(x) \right) = ( \lambda^{(0)}_n+ \epsilon \lambda^{(1)}_n) \left( f^{(0)}_n(x) + \epsilon f^{(1)}_n(x) \right) \,\!$。

取至一階，展開這方程。經過一番運算，可以得到

$D^{(1)} f^{(0)}_n (x) +D^{(0)}f^{(1)}_n (x)=\lambda^{(0)}_n f^{(1)}_n (x) +\lambda^{(1)}_n f^{(0)}_n (x)\,\!$。(5)

由於$\{f^{(0)}_i (x)\}\,\!$是一個完備集，$f^{(1)}_n (x)\,\!$可以寫為

$f^{(1)}_n (x)=\sum_{i\ne n} C_i f^{(0)}_i (x)\,\!$。(6)

請注意，這方程右手邊的總和表達式，並不含有$f^{(0)}_n (x)\,\!$項目。任何$f^{(0)}_n (x)\,\!$的貢獻，可以與公式 (4)的零階項目相合併。

將公式 (6)代入公式 (5)，可以得到

$(D^{(1)} - \lambda^{(1)}_n) f^{(0)}_n (x)=\lambda^{(0)}_n \sum_{i\ne n} C_i f^{(0)}_i (x) - D^{(0)}\sum_{i\ne n} C_i f^{(0)}_i (x)=\sum_{i\ne n} (\lambda^{(0)}_n - \lambda^{(0)}_i) C_i f^{(0)}_i (x)\,\!$。

將這方乘式兩邊都乘以$f^{(0)}_j (x)\,\!$，再隨著$x\,\!$積分，利用正交歸一性，可以得到

$\int\, f^{(0)}_j (x) D^{(1)}f^{(0)}_n (x)\, dx=\sum_{i\ne n} (\lambda^{(0)}_n - \lambda^{(0)}_i) C_i \int\, f^{(0)}_j (x) f^{(0)}_i (x)=(\lambda^{(0)}_n - \lambda^{(0)}_j) C_j \,\!$。

稍加編排，改變下標$j\,\!$為$m\,\!$。那麼，一階本徵函數修正$f^{(1)}_n (x)\,\!$可以寫為

$f^{(1)}_n(x) = \sum_{m \ne n} \frac {f^{(0)}_m (x)} {\lambda^{(0)}_n- \lambda^{(0)}_m} \int\, f^{(0)}_m(y) D^{(1)} f^{(0)}_n(y) \,dy\,\!$。

## 參閱

- 多體微擾理論

## 外部連結

- 攝動方法簡介 作者Mark H. Holmes
- 第二章：攝動方法簡介作者Johan Byström,Lars-Erik Persson，及Fredrik Strömberg