摄动理论

✍ dations ◷ 2025-04-02 19:14:29 #泛函分析,数学物理,微扰理论,微分方程

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

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

相关

  • 胸腔医学人体解剖学 - 人体生理学 组织学 - 胚胎学 人体寄生虫学 - 免疫学 病理学 - 病理生理学 细胞学 - 营养学 流行病学 - 药理学 - 毒理学呼吸医学为探讨呼吸器官疾病之一门学
  • 甲烷生成产甲烷作用,又称甲烷生成,指合成甲烷是微生物代谢的重要的和广泛的形式。可以生成甲烷的微生物称作产甲烷菌(英语:Methanogen)。这些微生物都属于原核生物中的古菌域,这是在系统发
  • 硫化钴硫化钴是硫钴二元化合物组成的混合物,化学式可以看作CoS。它由硫化钠和Co2+盐溶液反应,沉淀得到的CoS组成为Co(SH,OH)2,放置后该非晶态沉淀转变为晶态的CoxS(x<1)和Co9S8的混合物
  • 坦干伊喀湖耀鲶坦干伊喀湖耀鲶,为辐鳍鱼纲鲶形目塘虱鱼科的其中一种,被IUCN列为濒危保育类动物,分布于非洲坦干伊喀湖流域,栖息在沙泥底质底层水域,体长可达11.8公分,生活习性不明。 维基物种中
  • 鲍勃·卡尔鲍勃·卡尔,或译卜卡(英语:Bob Carr, 1947年9月28日-),全名罗伯特·约翰·卡尔(英语:Robert John Carr)是澳大利亚工党籍政治家,曾担任新南威尔士州州长(1995年4月-2005年8月),及澳大利亚外
  • 李代芳李代芳(1894年9月18日-1975年7月5日),字芸轩,山东青岛人,政治人物、企业家。曾任青岛市参议会参议长、青岛市商会理事长、制宪国大代表等职。1894年9月18日,李代芳生于山东省莱州府
  • 德·狄维士德·狄维士(Douglas L. de Vos,1964年-),出生于密歇根州的大急流市,1986年毕业于普杜大学管理系。现任安达高及安利公司总裁,安利创始人理查·狄维士第三子。狄维士曾担任安利公司
  • 史梦蛟史梦蛟,生卒年待考,约干嘉年间在世,清中期地方官、刊刻家、历史学家。浙江余姚人,祖为东钱湖史氏,南宋时迁居余姚,成为两城史氏。以辽州知州调任山西,清道光六年54岁时任忻州知州。
  • 乌姆朗格索乌姆朗格索(Umrangso),是印度阿萨姆邦North Cachar Hills县的一个城镇。总人口9024(2001年)。该地2001年总人口9024人,其中男性4946人,女性4078人;0—6岁人口1273人,其中男622人,女651
  • 博彩交易所博彩交易所(Bet Exchange),是一个可以让参与者自由买卖、处分彩券的交易所。他能撮合不同赔率的彩券,具交易所的功能,为彩券提供可作二手交易的次级市场,而不是彩券发行的初级市场