摄动理论

✍ dations ◷ 2025-04-26 13:54:22 #泛函分析,数学物理,微扰理论,微分方程

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

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

相关

  • 反比例函数一般地,如果两个变量x,y之间的关系可以表示为 y = k x
  • 联邦区域发展部(已裁撤)俄罗斯联邦区域发展部(俄语:Министерство регионального развития Российской Федерации – Минрегион Ро
  • 火焰谷国家休闲区火焰谷国家休闲区(英语:Flaming Gorge National Recreation Area),也作弗莱明峡谷国家休闲区,是一个美国的国家休闲区(英语:National Recreation Area),位于怀俄明州与犹他州。位于休
  • .укр.укр (DNS名称为".xn--j1amh")是乌克兰的乌克兰语国际化国家及地区顶级域。而乌克兰的传统国家及地区顶级域是.ua。A .ac .ad .ae .af .ag .ai .al .am .ao .aq .ar .as
  • 特雷弗·派克特雷弗·派克(英语:Trevor Pyke,1932年-1997年),英国音频工程师。他因电影从海底出击提名奥斯卡最佳音响效果奖。
  • 陈加元陈加元(1953年6月-),男,浙江永康人,中华人民共和国政治人物,曾任浙江省人民政府副省长、党组副书记,浙江省政协副主席。
  • 大仙术士 李白《大仙术士 李白》第一集封面《大仙术士 李白》,是台湾漫画家叶明轩的漫画作品,为根据唐朝诗人李白生平改编的奇幻漫画,在台湾角川《FORCE原力志》连载。剧情中还出现了杜甫、
  • 池殷成池殷成(韩语:지은성,1991年8月3日-),韩国男演员,另一译名为池恩诚。2013年出道,2016年在KBS2电视剧《我心中的花雨》中饰演满腔正义热血的法律系大学生朴善浩一角。因此角为该剧的主
  • 无人知晓的夏日清晨《无人知晓》(日语:誰も知らない),港译《谁知赤子心》,是一部2004年日本剧情片,本片的导演是是枝裕和,故事改编自1988年的日本巢鸭儿童遗弃事件。主要演员有柳乐优弥,北浦爱,木村飞影
  • 雅罗斯瓦夫·卡钦斯基雅罗斯瓦夫·亚历山大·卡钦斯基(波兰语:Jarosław Aleksander Kaczyński,1949年6月18日-)波兰政治家,法律与公正党党魁,波兰总理(2006年7月14日—2007年11月16日),曾任参议员和众议