摄动理论

✍ dations ◷ 2025-08-06 10:24:31 #泛函分析,数学物理,微扰理论,微分方程

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

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

相关

  • 卢·蒙特利卢·蒙特利(英语:Louis J. Montulli II)是一位以开发网页浏览器而闻名的程序员。在1991年和1992年,他与堪萨斯大学的迈克尔·格罗贝兹和查尔斯·雷扎克合作开发了一个名为Lynx的
  • 首尔市立交响乐团首尔市立交响乐团成立于1945年,是韩国最古老和有名的交响乐团。它的主要表演场地是世宗文化会馆。现任音乐总监是郑明勋先生。
  • 皇家交易所伦敦皇家交易所(Royal Exchange)是位于伦敦市的商业中心的一处交易所。由商人托马斯·格雷沙姆始建于16世纪,梯形平面,两侧的康希尔(Cornhill)和针线街聚集在伦敦市的心脏银行交叉
  • 海德拉巴海得拉巴(Hyderabad)可以指:
  • 瑞美关系瑞典-美国关系,是指瑞典王国与美利坚合众国之间的双边关系(英语:Bilateral relations)。两者的联系起始于17世纪的瑞典在美殖民活动,美国独立战争之后,瑞典是首批承认美国的国家之
  • 波尔图 (消歧义)波尔图(葡萄牙语:Porto)可以指:
  • 采尔湖 (魏尔海姆-雄高县)坐标:47°51′45″N 11°3′50″E / 47.86250°N 11.06389°E / 47.86250; 11.06389采尔湖(德语:Zellsee),是德国的湖泊,位于该国东南部,由巴伐利亚州负责管辖,处于魏尔海姆-雄高县,
  • 罗亨利罗亨利,GCB,GCMG(Henry Brougham Loch, 1st Baron Loch,1827年5月23日-1900年6月20日),英国军人、殖民地总督。第二次鸦片战争时,任额尔金的私人秘书。罗亨利、巴夏礼等39人在通州
  • 甘友信甘友信,广东保昌人,明朝政治人物、进士出身。洪武十八年,登乙丑科进士,授县学训导。
  • 系统尺度展开系统尺度展开,又称van Kampen's展开或者Ω-展开,是由Nico van Kampen开创运用于随机过程分析的数学方法。它能够对一个具有非线性变化率系统的主方程的解进行估计。这种展开的