摄动理论

摄动理论使用一些特別的数学方法來對於很多不具精确解的问题給出近似解，这些方法从相关的較簡單问题的精确解开始入手。摄动理论將原本問題分為具有精確解的較簡單部分與不具精確解的微扰部分。摄动理论适用的问题通常具有以下性質：通过加入一个微扰项於較簡單部分的數學表述，可以計算出整個問題的近似解。 摄动理论计算出来的解答通常会表达为一个微小参数的冪級數。摄动理论解答与精确解之间的差别，可以用这微小参数来做数量比较。冪級數的第一个项目是精确解的解答。后面的项目描述解答的修正。这修正是因为精确解与原本问题的「完全解」之间的误差而产生的。更正式地，完全解 A {\displaystyle A\,\!} 的近似可以表達为一个級數： A = ϵ 0 A 0 + ϵ 1 A 1 + ϵ 2 A 2 + ⋯ {\displaystyle A=\epsilon ^{0}A_{0}+\epsilon ^{1}A_{1}+\epsilon ^{2}A_{2}+\cdots \,\!} ； 在這例子裏， A 0 {\displaystyle A_{0}\,\!} 是簡單又有「精確解」的問題的精確解， A 1 , A 2 , {\displaystyle A_{1},\,A_{2},\,\!} 代表由某种系统程序反覆地找到的高阶项目修正。因为 ϵ {\displaystyle \epsilon \,\!} 的值很微小，这些高阶项目修正应该会越来越不重要。 摄动理论的标准阐述主要是以微扰的阶数来分辨：一阶摄动理论或二阶摄动理论。再来就是以微扰的简并度来分辨：无简并或有简并。有简并的摄动，又称为奇异摄动（singular perturbation），比较难解，必须用到更进阶的理论。 本段落讲述微分方程的一阶微扰理论。为了简单易解，假设零微扰系统的解答是不简并的。 许多常微分方程或偏微分方程可以表达为 D g ( x ) = λ g ( x ) {\displaystyle Dg(x)=\lambda g(x)\,\!} ；(1) 其中， D {\displaystyle D\,\!} 是某特定微分算子， λ {\displaystyle \lambda \,\!} 是其本征值。 假设微分算子可以写为 D = D ( 0 ) + ϵ D ( 1 ) {\displaystyle D=D^{(0)}+\epsilon D^{(1)}\,\!} ； 其中， ϵ {\displaystyle \epsilon \,\!} 是微小的度量。 又假设我们已知道 D ( 0 ) {\displaystyle D^{(0)}\,\!} 的解答的完备集 { f i ( 0 ) ( x ) } {\displaystyle \{f_{i}^{(0)}(x)\}\,\!} ；其中，解答 f i ( 0 ) ( x ) {\displaystyle f_{i}^{(0)}(x)\,\!} 是 D ( 0 ) {\displaystyle D^{(0)}\,\!} 的本征值为 λ i ( 0 ) {\displaystyle \lambda _{i}^{(0)}\,\!} 的本征函数。用方程表达， D ( 0 ) f i ( 0 ) ( x ) = λ i ( 0 ) f i ( 0 ) ( x ) {\displaystyle D^{(0)}f_{i}^{(0)}(x)=\lambda _{i}^{(0)}f_{i}^{(0)}(x)\,\!} 。 还有，这一集合的解答 { f i ( 0 ) ( x ) } {\displaystyle \{f_{i}^{(0)}(x)\}\,\!} 形成一个正交归一集： ∫ f i ( 0 ) ( x ) f j ( 0 ) ( x ) d x = δ i j {\displaystyle \int f_{i}^{(0)}(x)f_{j}^{(0)}(x)\,dx=\delta _{ij}\,\!} ； 其中， δ i j {\displaystyle \delta _{ij}\,\!} 是克羅內克函數。 取至零阶，完全解 g ( x ) {\displaystyle g(x)\,\!} 应该相当接近集合里一个零微扰解。设定这零微扰解为 f n ( 0 ) ( x ) {\displaystyle f_{n}^{(0)}(x)\,\!} 。用方程表达， g ( x ) = f n ( 0 ) ( x ) + O ( ϵ ) {\displaystyle g(x)=f_{n}^{(0)}(x)+{\mathcal {O}}(\epsilon )\,\!} ； 其中， O {\displaystyle {\mathcal {O}}\,\!} 采用大O符號来描述函数的渐近行为。 完全解的本征值也可近似为 λ = λ n ( 0 ) + O ( ϵ ) {\displaystyle \lambda =\lambda _{n}^{(0)}+{\mathcal {O}}(\epsilon )\,\!} 。 将完全解 g ( x ) {\displaystyle g(x)\,\!} 写为零微扰解的线性组合， g ( x ) = ∑ m c m f m ( 0 ) ( x ) {\displaystyle g(x)=\sum _{m}c_{m}f_{m}^{(0)}(x)\,\!} ；(2) 其中，除了 c n {\displaystyle c_{n}\,\!} 以外，所有的常数 c m ,   m ≠ n {\displaystyle c_{m},\ m\neq n\,\!} 的值是 O ( ϵ ) {\displaystyle {\mathcal {O}}(\epsilon )\,\!} ；只有 c n {\displaystyle c_{n}\,\!} 的值是 O ( 1 ) {\displaystyle {\mathcal {O}}(1)\,\!} 。 将公式 (2)代入公式 (1)，乘以 f n ( 0 ) ( x ) {\displaystyle f_{n}^{(0)}(x)\,\!} ，利用正交归一性，可以得到 λ n ( 0 ) c n + ϵ ∑ m c m ∫ f n ( 0 ) ( x ) D ( 1 ) f m ( 0 ) ( x ) d x = λ c n {\displaystyle \lambda _{n}^{(0)}c_{n}+\epsilon \sum _{m}c_{m}\int f_{n}^{(0)}(x)D^{(1)}f_{m}^{(0)}(x)\,dx=\lambda c_{n}\,\!} 。 这可以很容易地改变为一个简单的线性代数问题，一个寻找矩阵的本征值的问题：给予 ∑ m A n m c m = λ c n {\displaystyle \sum _{m}A_{nm}c_{m}=\lambda c_{n}\!\,\!} ，求 λ {\displaystyle \lambda \,\!} ；其中， A n m {\displaystyle A_{nm}\,\!} 是矩阵元素： A n m = λ n ( 0 ) δ n m + ϵ ∫ f n ( 0 ) ( x ) D ( 1 ) f m ( 0 ) ( x ) d x {\displaystyle A_{nm}=\lambda _{n}^{(0)}\delta _{nm}+\epsilon \int f_{n}^{(0)}(x)D^{(1)}f_{m}^{(0)}(x)\,dx\,\!} 。 我们并不需要解析整个矩阵。注意到线性方程裡的每一个 c m {\displaystyle c_{m}\,\!} 都是 O ( ϵ ) {\displaystyle {\mathcal {O}}(\epsilon )\,\!} ；只有 c n {\displaystyle c_{n}\,\!} 的值是 O ( 1 ) {\displaystyle {\mathcal {O}}(1)\,\!} 。所以，取至 ϵ {\displaystyle \epsilon \,\!} 一阶，线性方程可以很容易地解析为 λ = λ n ( 0 ) + ϵ ∫ f n ( 0 ) ( x ) D ( 1 ) f n ( 0 ) ( x ) d x {\displaystyle \lambda =\lambda _{n}^{(0)}+\epsilon \int f_{n}^{(0)}(x)D^{(1)}f_{n}^{(0)}(x)\,dx\,\!} 。(3) 这就是一阶摄动理论的本征值解答。一阶本征值数修正是 λ n ( 1 ) = ∫ f n ( 0 ) ( x ) D ( 1 ) f n ( 0 ) ( x ) d x {\displaystyle \lambda _{n}^{(1)}=\int f_{n}^{(0)}(x)D^{(1)}f_{n}^{(0)}(x)\,dx\,\!} 。 取至一阶，函数 g ( x ) {\displaystyle g(x)\,\!} 可以用类似的推理求得。设定 g ( x ) = f n ( 0 ) ( x ) + ϵ f n ( 1 ) ( x ) {\displaystyle g(x)=f_{n}^{(0)}(x)+\epsilon f_{n}^{(1)}(x)\,\!} 。(4) 那麼，公式 (1)变为 ( D ( 0 ) + ϵ D ( 1 ) ) ( f n ( 0 ) ( x ) + ϵ f n ( 1 ) ( x ) ) = ( λ n ( 0 ) + ϵ λ n ( 1 ) ) ( f n ( 0 ) ( x ) + ϵ f n ( 1 ) ( x ) ) {\displaystyle \left(D^{(0)}+\epsilon D^{(1)}\right)\left(f_{n}^{(0)}(x)+\epsilon f_{n}^{(1)}(x)\right)=(\lambda _{n}^{(0)}+\epsilon \lambda _{n}^{(1)})\left(f_{n}^{(0)}(x)+\epsilon f_{n}^{(1)}(x)\right)\,\!} 。 取至一阶，展开这方程。经过一番运算，可以得到 D ( 1 ) f n ( 0 ) ( x ) + D ( 0 ) f n ( 1 ) ( x ) = λ n ( 0 ) f n ( 1 ) ( x ) + λ n ( 1 ) f n ( 0 ) ( x ) {\displaystyle D^{(1)}f_{n}^{(0)}(x)+D^{(0)}f_{n}^{(1)}(x)=\lambda _{n}^{(0)}f_{n}^{(1)}(x)+\lambda _{n}^{(1)}f_{n}^{(0)}(x)\,\!} 。(5) 由于 { f i ( 0 ) ( x ) } {\displaystyle \{f_{i}^{(0)}(x)\}\,\!} 是一个完备集， f n ( 1 ) ( x ) {\displaystyle f_{n}^{(1)}(x)\,\!} 可以写为 f n ( 1 ) ( x ) = ∑ i ≠ n C i f i ( 0 ) ( x ) {\displaystyle f_{n}^{(1)}(x)=\sum _{i\neq n}C_{i}f_{i}^{(0)}(x)\,\!} 。(6) 请注意，这方程右手边的总和表达式，并不含有 f n ( 0 ) ( x ) {\displaystyle f_{n}^{(0)}(x)\,\!} 项目。任何 f n ( 0 ) ( x ) {\displaystyle f_{n}^{(0)}(x)\,\!} 的贡献，可以与公式 (4)的零階項目相合并。 将公式 (6)代入公式 (5)，可以得到 ( D ( 1 ) − λ n ( 1 ) ) f n ( 0 ) ( x ) = λ n ( 0 ) ∑ i ≠ n C i f i ( 0 ) ( x ) − D ( 0 ) ∑ i ≠ n C i f i ( 0 ) ( x ) = ∑ i ≠ n ( λ n ( 0 ) − λ i ( 0 ) ) C i f i ( 0 ) ( x ) {\displaystyle (D^{(1)}-\lambda _{n}^{(1)})f_{n}^{(0)}(x)=\lambda _{n}^{(0)}\sum _{i\neq n}C_{i}f_{i}^{(0)}(x)-D^{(0)}\sum _{i\neq n}C_{i}f_{i}^{(0)}(x)=\sum _{i\neq n}(\lambda _{n}^{(0)}-\lambda _{i}^{(0)})C_{i}f_{i}^{(0)}(x)\,\!} 。 将这方乘式两边都乘以 f j ( 0 ) ( x ) {\displaystyle f_{j}^{(0)}(x)\,\!} ，再随著 x {\displaystyle x\,\!} 积分，利用正交归一性，可以得到 ∫ f j ( 0 ) ( x ) D ( 1 ) f n ( 0 ) ( x ) d x = ∑ i ≠ n ( λ n ( 0 ) − λ i ( 0 ) ) C i ∫ f j ( 0 ) ( x ) f i ( 0 ) ( x ) = ( λ n ( 0 ) − λ j ( 0 ) ) C j {\displaystyle \int \,f_{j}^{(0)}(x)D^{(1)}f_{n}^{(0)}(x)\,dx=\sum _{i\neq n}(\lambda _{n}^{(0)}-\lambda _{i}^{(0)})C_{i}\int \,f_{j}^{(0)}(x)f_{i}^{(0)}(x)=(\lambda _{n}^{(0)}-\lambda _{j}^{(0)})C_{j}\,\!} 。 稍加编排，改变下标 j {\displaystyle j\,\!} 为 m {\displaystyle m\,\!} 。那麼，一阶本征函数修正 f n ( 1 ) ( x ) {\displaystyle f_{n}^{(1)}(x)\,\!} 可以写为 f n ( 1 ) ( x ) = ∑ m ≠ n f m ( 0 ) ( x ) λ n ( 0 ) − λ m ( 0 ) ∫ f m ( 0 ) ( y ) D ( 1 ) f n ( 0 ) ( y ) d y {\displaystyle f_{n}^{(1)}(x)=\sum _{m\neq n}{\frac {f_{m}^{(0)}(x)}{\lambda _{n}^{(0)}-\lambda _{m}^{(0)}}}\int \,f_{m}^{(0)}(y)D^{(1)}f_{n}^{(0)}(y)\,dy\,\!} 。 多體微擾理論 William E. Wiesel. Modern Astrodynamics. Ohio: Aphelion Press. 2010: 107. ISBN 978-145378-1470.  摄动方法简介 （页面存档备份，存于互联网档案馆）作者Mark H. Holmes 第二章：摄动方法简介作者Johan Byström,Lars-Erik Persson，及Fredrik Strömberg