摄动理论由用于寻找无法精确求解的问题的近似解的数学方法组成,这些方法从相关问题的精确解开始入手。摄动理论适用于可以通过加入一个微扰项到一个可以精确求解的问题上而表述的问题。 摄动理论计算出来的解答通常会表达为一个微小参数的幂级数。摄动理论解答与精确解之间的差别,可以用这微小参数来做数量比较。幂级数的第一个项目是精确解的解答。后面的项目描述解答的修正。这修正是因为精确解与原本问题的“完全解”之间的误差而产生的。更正式地,完全解
的近似可以表达为一个级数:
;
在这例子里,
是简单又有“精确解”的问题的精确解,
代表由某种系统程序反复地找到的高阶项目修正。因为
的值很微小,这些高阶项目修正应该会越来越不重要。
微扰阶数
摄动理论的标准阐述主要是以微扰的阶数来分辨:一阶摄动理论或二阶摄动理论。再来就是以微扰的简并度来分辨:无简并或有简并。有简并的摄动,又称为奇异摄动(singular perturbation),比较难解,必须用到更进阶的理论。
一阶无简并摄动理论
本段落讲述微分方程的一阶微扰理论。为了简单易解,假设零微扰系统的解答是不简并的。
一阶本征值修正
许多常微分方程或偏微分方程可以表达为
;(1)
其中,
是某特定微分算子,
是其本征值。 假设微分算子可以写为
;
其中,
是微小的度量。 又假设我们已知道
的解答的完备集
;其中,解答
是
的本征值为
的本征函数。用方程表达,
。
还有,这一集合的解答
形成一个正交归一集:
;
其中,
是克罗内克函数。 取至零阶,完全解
应该相当接近集合里一个零微扰解。设定这零微扰解为
。用方程表达,
;
其中,
采用大O符号来描述函数的渐近形为。 完全解的本征值也可近似为
。
将完全解
写为零微扰解的线性组合,
;(2)
其中,除了
以外,所有的常数
的值是
;只有
的值是
。 将公式 (2)代入公式 (1),乘以
,利用正交归一性,可以得到
。
这可以很容易地改变为一个简单的线性代数问题,一个寻找矩阵的本征值的问题:给予
,求
;其中,
是矩阵元素:
。
我们并不需要解析整个矩阵。注意到线性方程里的每一个
都是
;只有
的值是
。所以,取至
一阶,线性方程可以很容易地解析为
。(3)
这就是一阶摄动理论的本征值解答。一阶本征值数修正是
。
一阶本征函数修正
取至一阶,函数
可以用类似的推理求得。设定
。(4)
那么,公式 (1)变为
。
取至一阶,展开这方程。经过一番运算,可以得到
。(5)
由于
是一个完备集,
可以写为
。(6)
请注意,这方程右手边的总和表达式,并不含有
项目。任何
的贡献,可以与公式 (4)的零阶项目相合并。 将公式 (6)代入公式 (5),可以得到
。
将这方乘式两边都乘以
,再随著
积分,利用正交归一性,可以得到
。
稍加编排,改变下标
为
。那么,一阶本征函数修正
可以写为
。
参阅
- 多体摄动理论
外部连接
- 摄动方法简介作者Mark H. Holmes
- 第二章:摄动方法简介作者Johan Byström,Lars-Erik Persson,及Fredrik Strömberg
![[程序员段子] 1927年12月31日午夜一秒谜案](https://img.cnbeta.com/newsimg/110812/14571201438079743.png)
(1928年1月1日0点0分,当地时间由
我叫IE6,都用过Windows XP吧?用过,就一定知道我长什么样儿。别问我是干嘛的,问就说明你OUT至少十年了。告诉你,我是让你能上网的,能看新闻玩网游的,在浏览器界辛苦耕耘十年的小老儿。
十年啊!你想过十年的变化没?十年前,我来自微软,和搭档Windows XP初露锋芒,跟它一起红遍全球。
十年前,网页就是摆几张图,弄几坨字。瞧瞧现在,看个小片儿,偷棵小菜儿,听首小曲儿,团点小货儿,织条小博儿,好玩的网站,真是太多太多啦……
十年前,我走在互联网的大道上,如履平地!现在呢,网页里要搭配更多的插件,何况盗版商们还不忘在我身上多加几个按钮!我走起来如履薄冰啊!
我刚出道的时候,谁知道什么互联网的标准呢?后来,标准普及了,越来越健全,但我已步入晚年。
现在那些安全的浏览器,套在我身上一点都不给力。因为我出道的时候,病毒没这么凶,杀些儿小怪还可以,想我对付BOSS?哪有那么好的装备啊……
能为大家服务十年,已付出了我毕生的心血,尤其感谢中国,中国人民最可爱,因为直到今日,依然对我不离不弃。可我得说,我还是要退休了。十年啊,拉登都顶不住了,何况我一个IE6呢?
我的接班人,拥有强劲性能的老八和老九,都是值得你们信任的浏览器。
是时候说再见了,谢谢大家颁给我这个奖项,我爱你们! 附全球IE6倒计时网站中国版:
