摄动理论

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

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

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

相关

  • 吡硫吡硫
  • 法国国立现代艺术美术馆国立现代艺术博物馆(法语:Musée National d'Art Moderne)是法国的国家现代美术的博物馆,为庞毕度国家艺术和文化中心的组成部分,主要展区位于巴黎第四区的庞毕度中心。博物馆的
  • 鲁道夫·朱利安尼鲁道夫·威廉·路易斯·“鲁迪”·朱利安尼,KBE(英语:Rudolph William Louis "Rudy" Giuliani,1944年5月28日-)是美国律师、检察官、商人及共和党的美国政治人物,出身于纽约州。朱
  • 稻属已知18种,见内文。稻属(学名:Oryza)是禾本科的一属,有18个物种,包括2个栽培种和17个原生种。分布于全球的热带与亚热带地区。稻属是由卡尔·林奈在1753年确定,此后稻属物种的归属与
  • 指定幸存者指定幸存者(英语:Designated survivor)是在美国总统和国家其他高级领导人聚集在同一处时故意安排的一名缺席内阁阁员,此举是为了确保在总统、副总统、国会领袖及内阁阁员全体丧
  • 若瑟·卡普里奥若瑟·卡普里奥(英语:Giuseppe Caprio,中文名:高理耀;1914年11月15日-2005年10月14日),是意大利的天主教司铎级枢机及耶路撒冷圣墓骑士团总团长。卡普里奥于1939年4月8日在意大利阿
  • 月光菩萨月光菩萨(梵语:चन्द्रप्रभ Candraprabha),又称月光遍照菩萨、月净菩萨或月光普照菩萨,出现于大乘佛教经典的人物,最早见于《佛说决定毘尼经》。比较为人熟知为药师佛的
  • 朱塞佩·维塔利朱塞佩·维塔利(1875年8月26日-1932年2月29日),意大利数学家,因实数不可测集的存在性上的维塔利定理而出名。他的维塔利覆盖引理也是测度理论的一条基本原理。维塔利生于拉韦纳,卒
  • 那嵩那嵩(?-1659年),明朝末年清朝初年云南元江人,字维岳,一字岳生,元江府土知府,傣族。一说那嵩原为元江土司族属,世守他郎寨,土司绝嗣,那嵩杀其族长老木龙遂为土知府;一说是“世袭土府”。那
  • 李崴李崴(1953年-),汉族,中华人民共和国政治人物、第十一届全国政协委员,中国国民党革命委员会党员,曾担任十一届全国政协常务委员、民革中央委员、广东省委副主委。2008年,当选第十一届