摄动理论

✍ dations ◷ 2025-04-04 11:24:32 #泛函分析,数学物理,微扰理论,微分方程

摄动理论使用一些特别的数学方法来对于很多不具精确解的问题给出近似解,这些方法从相关的较简单问题的精确解开始入手。摄动理论将原本问题分为具有精确解的较简单部分与不具精确解的微扰部分。摄动理论适用的问题通常具有以下性质:通过加入一个微扰项于较简单部分的数学表述,可以计算出整个问题的近似解。

摄动理论计算出来的解答通常会表达为一个微小参数的幂级数。摄动理论解答与精确解之间的差别,可以用这微小参数来做数量比较。幂级数的第一个项目是精确解的解答。后面的项目描述解答的修正。这修正是因为精确解与原本问题的“完全解”之间的误差而产生的。更正式地,完全解 A {\displaystyle A\,\!} 的近似可以表达为一个级数:

在这例子里, A 0 {\displaystyle A_{0}\,\!} 是简单又有“精确解”的问题的精确解, A 1 , A 2 , {\displaystyle A_{1},\,A_{2},\,\!} 代表由某种系统程序反复地找到的高阶项目修正。因为 ϵ {\displaystyle \epsilon \,\!} 的值很微小,这些高阶项目修正应该会越来越不重要。

摄动理论的标准阐述主要是以微扰的阶数来分辨:一阶摄动理论或二阶摄动理论。再来就是以微扰的简并度来分辨:无简并或有简并。有简并的摄动,又称为奇异摄动(singular perturbation),比较难解,必须用到更进阶的理论。

本段落讲述微分方程的一阶微扰理论。为了简单易解,假设零微扰系统的解答是不简并的。

许多常微分方程或偏微分方程可以表达为

其中, D {\displaystyle D\,\!} 是某特定微分算子, λ {\displaystyle \lambda \,\!} 是其本征值。

假设微分算子可以写为

其中, ϵ {\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)}\,\!} 的本征函数。用方程表达,

还有,这一集合的解答 { f i ( 0 ) ( x ) } {\displaystyle \{f_{i}^{(0)}(x)\}\,\!} 形成一个正交归一集:

其中, δ i j {\displaystyle \delta _{ij}\,\!} 是克罗内克函数。

取至零阶,完全解 g ( x ) {\displaystyle g(x)\,\!} 应该相当接近集合里一个零微扰解。设定这零微扰解为 f n ( 0 ) ( x ) {\displaystyle f_{n}^{(0)}(x)\,\!} 。用方程表达,

其中, O {\displaystyle {\mathcal {O}}\,\!} 采用大O符号来描述函数的渐近行为。

完全解的本征值也可近似为

将完全解 g ( x ) {\displaystyle g(x)\,\!} 写为零微扰解的线性组合,

其中,除了 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)\,\!} ,利用正交归一性,可以得到

这可以很容易地改变为一个简单的线性代数问题,一个寻找矩阵的本征值的问题:给予 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}\,\!} 是矩阵元素:

我们并不需要解析整个矩阵。注意到线性方程里的每一个 c m {\displaystyle c_{m}\,\!} 都是 O ( ϵ ) {\displaystyle {\mathcal {O}}(\epsilon )\,\!} ;只有 c n {\displaystyle c_{n}\,\!} 的值是 O ( 1 ) {\displaystyle {\mathcal {O}}(1)\,\!} 。所以,取至 ϵ {\displaystyle \epsilon \,\!} 一阶,线性方程可以很容易地解析为

这就是一阶摄动理论的本征值解答。一阶本征值数修正是

取至一阶,函数 g ( x ) {\displaystyle g(x)\,\!} 可以用类似的推理求得。设定

那么,公式 (1)变为

取至一阶,展开这方程。经过一番运算,可以得到

由于 { f i ( 0 ) ( x ) } {\displaystyle \{f_{i}^{(0)}(x)\}\,\!} 是一个完备集, f n ( 1 ) ( x ) {\displaystyle f_{n}^{(1)}(x)\,\!} 可以写为

请注意,这方程右手边的总和表达式,并不含有 f n ( 0 ) ( x ) {\displaystyle f_{n}^{(0)}(x)\,\!} 项目。任何 f n ( 0 ) ( x ) {\displaystyle f_{n}^{(0)}(x)\,\!} 的贡献,可以与公式 (4)的零阶项目相合并。

将公式 (6)代入公式 (5),可以得到

将这方乘式两边都乘以 f j ( 0 ) ( x ) {\displaystyle f_{j}^{(0)}(x)\,\!} ,再随著 x {\displaystyle x\,\!} 积分,利用正交归一性,可以得到

稍加编排,改变下标 j {\displaystyle j\,\!} m {\displaystyle m\,\!} 。那么,一阶本征函数修正 f n ( 1 ) ( x ) {\displaystyle f_{n}^{(1)}(x)\,\!} 可以写为

相关

  • PAS染色法PAS染色法(Periodic Acid-Schiff stain)在组织学上,主要用来检测组织中的糖类,过碘酸把糖类相邻两个碳上的羟基氧化成醛基,再用Schiff试剂和醛基反应使呈现紫红色。通常用于霉菌
  • 昼夜节律紊乱昼夜节律睡眠障碍(CRSD)是影响睡眠时间及其他方面的一种睡眠障碍 。患有昼夜节律睡眠障碍者无法在一般工作、学校和社会需求所需要的入睡时间睡着,或起床时间起床。如果允许他
  • 捷克君主列表这是一份捷克君主的列表。波希米亚,摩拉维亚,西里西亚和卢萨蒂亚在历史上都是以捷克人为主的波希米亚王国(前身是波希米亚公国)的领地。 请参见主条目:
  • The Big O《THE Big O》(日语:THE ビッグオー,另译:魅影巨神、魔神比古奥)是一部由日本日升公司制作的机器人动画,1999年10月13日到2000年1月19日在WOWOW电视台首播。另有由有贺等所作漫画,
  • 塞尔维亚社会主义共和国塞尔维亚社会主义共和国( 错误:{{lang-xx}}:文本有斜体标记(帮助))是南斯拉夫社会主义联邦共和国的一个加盟共和国,也是塞尔维亚共和国的前身。
  • 欧伊根·埃利希欧伊根·埃利希(Eugen Ehrlich,1862年9月14日-1922年5月2日),是一名奥地利法学家,也是公认的社会法学派创始人之一。他出生于法律世家,并在维也纳大学学习法律。1899年起,他在切尔诺
  • 施塔德尔贝格山坐标:46°50′10″N 16°01′50″E / 46.836077°N 16.030532°E / 46.836077; 16.030532施塔德尔贝格山(德语:Stadelberg),是中欧的山峰,位于斯洛文尼亚和奥地利接壤的边境,由罗加
  • 戴维斯·古根海姆戴维斯·古根海姆(Philip Davis Guggenheim,1964年-)是一位美国电影导演与制片人,亦是一名奥斯卡奖得主。古根海姆是奥斯卡奖得主查尔斯·古根海姆与马里恩·古根海姆的儿子。他
  • 殷姓殷姓是一个中国人的姓氏,在《百家姓》中排名第74位。
  • 阿尔卡季·帕夫洛维奇·罗森戈尔茨阿尔卡季·帕夫洛维奇·罗森戈尔茨(俄语:Аркадий Павлович Розенгольц,1889年11月4日-1938年3月15日),苏联政治家、布尔什维克。出生于维捷布斯克,是一个