摄动理论
摄动理论使用一些特别的数学方法来对于很多不具精确解的问题给出近似解,这些方法从相关的较简单问题的精确解开始入手。摄动理论将原本问题分为具有精确解的较简单部分与不具精确解的微扰部分。摄动理论适用的问题通常具有以下性质:通过加入一个微扰项于较简单部分的数学表述,可以计算出整个问题的近似解。
摄动理论计算出来的解答通常会表达为一个微小参数的幂级数。摄动理论解答与精确解之间的差别,可以用这微小参数来做数量比较。幂级数的第一个项目是精确解的解答。后面的项目描述解答的修正。这修正是因为精确解与原本问题的「完全解」之间的误差而产生的。更正式地,完全解formula_1的近似可以表达为一个级数:
formula_2;
在这例子里,formula_3是简单又有「精确解」的问题的精确解,formula_4代表由某种系统程序反复地找到的高阶项目修正。因为formula_5的值很微小,这些高阶项目修正应该会越来越不重要。
微扰阶数.
摄动理论的标准阐述主要是以微扰的阶数来分辨:一阶摄动理论或二阶摄动理论。再来就是以微扰的简并度来分辨:无简并或有简并。有简并的摄动,又称为奇异摄动(--
),比较难解,必须用到更进阶的理论。
一阶无简并摄动理论.
本段落讲述微分方程的一阶微扰理论。为了简单易解,假设零微扰系统的解答是不简并的。
一阶本征值修正.
许多常微分方程或偏微分方程可以表达为
formula_6;(1)
其中,formula_7是某特定微分算子,formula_8是其本征值。
假设微分算子可以写为
formula_9;
其中,formula_5是微小的度量。
又假设我们已知道formula_11的解答的完备集formula_12;其中,解答formula_13是formula_11的本征值为formula_15的本征函数。用方程表达,
formula_16。
还有,这一集合的解答formula_12形成一个正交归一集:
formula_18;
其中,formula_19是克罗内克函数。
取至零阶,完全解formula_20应该相当接近集合里一个零微扰解。设定这零微扰解为formula_21。用方程表达,
formula_22;
其中,formula_23采用大O符号来描述函数的渐近行为。
完全解的本征值也可近似为
formula_24。
将完全解formula_20写为零微扰解的线性组合,
formula_26;(2)
其中,除了formula_27以外,所有的常数formula_28的值是formula_29;只有formula_27的值是formula_31。
将公式 (2)代入公式 (1),乘以formula_32,利用正交归一性,可以得到
formula_33。
这可以很容易地改变为一个简单的线性代数问题,一个寻找矩阵的本征值的问题:给予
formula_34,求formula_8;其中,formula_36是矩阵元素:
formula_37。
我们并不需要解析整个矩阵。注意到线性方程里的每一个formula_38都是formula_29;只有formula_27的值是formula_31。所以,取至formula_5一阶,线性方程可以很容易地解析为
formula_43。(3)
这就是一阶摄动理论的本征值解答。一阶本征值数修正是
formula_44。
一阶本征函数修正.
取至一阶,函数formula_20可以用类似的推理求得。设定
formula_46。(4)
那么,公式 (1)变为
formula_47。
取至一阶,展开这方程。经过一番运算,可以得到
formula_48。(5)
由于formula_12是一个完备集,formula_50可以写为
formula_51。(6)
请注意,这方程右手边的总和表达式,并不含有formula_52项目。任何formula_52的贡献,可以与公式 (4)的零阶项目相合并。
将公式 (6)代入公式 (5),可以得到
formula_54。
将这方乘式两边都乘以formula_55,再随著formula_56积分,利用正交归一性,可以得到
formula_57。
稍加编排,改变下标formula_58为formula_59。那么,一阶本征函数修正formula_50可以写为
formula_61。
生成维基百科快照图片,大概需要3-30秒!