摄动理论, 此条目的主題是一般的数学, 关于量子力学的, 請見, 微擾理論, 量子力學, 使用一些特別的数学方法來對於很多不具精确解的问题給出近似解, 这些方法从相关的較簡單问题的精确解开始入手, 將原本問題分為具有精確解的較簡單部分與不具精確解的微扰部分, 适用的问题通常具有以下性質, 通过加入一个微扰项於較簡單部分的數學表述, 可以計算出整個問題的近似解, 计算出来的解答通常会表达为一个微小参数的冪級數, 解答与精确解之间的差别, 可以用这微小参数来做数量比较, 冪級數的第一个项目是精确解的解答, 后面的项目描. 此条目的主題是一般的数学摄动理论 关于量子力学的摄动理论 請見 微擾理論 量子力學 摄动理论使用一些特別的数学方法來對於很多不具精确解的问题給出近似解 这些方法从相关的較簡單问题的精确解开始入手 摄动理论將原本問題分為具有精確解的較簡單部分與不具精確解的微扰部分 1 摄动理论适用的问题通常具有以下性質 通过加入一个微扰项於較簡單部分的數學表述 可以計算出整個問題的近似解 摄动理论计算出来的解答通常会表达为一个微小参数的冪級數 摄动理论解答与精确解之间的差别 可以用这微小参数来做数量比较 冪級數的第一个项目是精确解的解答 后面的项目描述解答的修正 这修正是因为精确解与原本问题的 完全解 之间的误差而产生的 更正式地 完全解A displaystyle A 的近似可以表達为一个級數 A ϵ 0 A 0 ϵ 1 A 1 ϵ 2 A 2 displaystyle A epsilon 0 A 0 epsilon 1 A 1 epsilon 2 A 2 cdots 在這例子裏 A 0 displaystyle A 0 是簡單又有 精確解 的問題的精確解 A 1 A 2 displaystyle A 1 A 2 代表由某种系统程序反覆地找到的高阶项目修正 因为ϵ displaystyle epsilon 的值很微小 这些高阶项目修正应该会越来越不重要 目录 1 微扰阶数 2 一阶无简并摄动理论 2 1 一阶本征值修正 2 2 一阶本征函数修正 3 参阅 4 參考資料 5 外部链接微扰阶数 编辑摄动理论的标准阐述主要是以微扰的阶数来分辨 一阶摄动理论或二阶摄动理论 再来就是以微扰的简并度来分辨 无简并或有简并 有简并的摄动 又称为奇异摄动 singular perturbation 比较难解 必须用到更进阶的理论 一阶无简并摄动理论 编辑本段落讲述微分方程的一阶微扰理论 为了简单易解 假设零微扰系统的解答是不简并的 一阶本征值修正 编辑 许多常微分方程或偏微分方程可以表达为 D g x l g x displaystyle Dg x lambda g x 1 其中 D displaystyle D 是某特定微分算子 l displaystyle lambda 是其本征值 假设微分算子可以写为 D D 0 ϵ D 1 displaystyle D D 0 epsilon D 1 其中 ϵ 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 的本征值为l i 0 displaystyle lambda i 0 的本征函数 用方程表达 D 0 f i 0 x l i 0 f i 0 x displaystyle D 0 f i 0 x lambda i 0 f i 0 x 还有 这一集合的解答 f i 0 x displaystyle f i 0 x 形成一个正交归一集 f i 0 x f j 0 x d x d i j displaystyle int f i 0 x f j 0 x dx delta ij 其中 d i j displaystyle delta ij 是克羅內克函數 取至零阶 完全解g x displaystyle g x 应该相当接近集合里一个零微扰解 设定这零微扰解为f n 0 x displaystyle f n 0 x 用方程表达 g x f n 0 x O ϵ displaystyle g x f n 0 x mathcal O epsilon 其中 O displaystyle mathcal O 采用大O符號来描述函数的渐近行为 完全解的本征值也可近似为 l l n 0 O ϵ displaystyle lambda lambda n 0 mathcal O epsilon 将完全解g x displaystyle g x 写为零微扰解的线性组合 g x m c m f m 0 x displaystyle g x sum m c m f m 0 x 2 其中 除了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 利用正交归一性 可以得到 l n 0 c n ϵ m c m f n 0 x D 1 f m 0 x d x l c n displaystyle lambda n 0 c n epsilon sum m c m int f n 0 x D 1 f m 0 x dx lambda c n 这可以很容易地改变为一个简单的线性代数问题 一个寻找矩阵的本征值的问题 给予 m A n m c m l c n displaystyle sum m A nm c m lambda c n 求l displaystyle lambda 其中 A n m displaystyle A nm 是矩阵元素 A n m l n 0 d n m ϵ f n 0 x D 1 f m 0 x d x displaystyle A nm lambda n 0 delta nm epsilon int f n 0 x D 1 f m 0 x dx 我们并不需要解析整个矩阵 注意到线性方程裡的每一个c m displaystyle c m 都是O ϵ displaystyle mathcal O epsilon 只有c n displaystyle c n 的值是O 1 displaystyle mathcal O 1 所以 取至ϵ displaystyle epsilon 一阶 线性方程可以很容易地解析为 l l n 0 ϵ f n 0 x D 1 f n 0 x d x displaystyle lambda lambda n 0 epsilon int f n 0 x D 1 f n 0 x dx 3 这就是一阶摄动理论的本征值解答 一阶本征值数修正是 l n 1 f n 0 x D 1 f n 0 x d x displaystyle lambda n 1 int f n 0 x D 1 f n 0 x dx 一阶本征函数修正 编辑 取至一阶 函数g x displaystyle g x 可以用类似的推理求得 设定 g x f n 0 x ϵ f n 1 x displaystyle g x f n 0 x epsilon f n 1 x 4 那麼 公式 1 变为 D 0 ϵ D 1 f n 0 x ϵ f n 1 x l n 0 ϵ l n 1 f n 0 x ϵ f n 1 x displaystyle left D 0 epsilon D 1 right left f n 0 x epsilon f n 1 x right lambda n 0 epsilon lambda n 1 left f n 0 x epsilon f n 1 x right 取至一阶 展开这方程 经过一番运算 可以得到 D 1 f n 0 x D 0 f n 1 x l n 0 f n 1 x l n 1 f n 0 x displaystyle D 1 f n 0 x D 0 f n 1 x lambda n 0 f n 1 x lambda n 1 f n 0 x 5 由于 f i 0 x displaystyle f i 0 x 是一个完备集 f n 1 x displaystyle f n 1 x 可以写为 f n 1 x i n C i f i 0 x displaystyle f n 1 x sum i neq n C i f i 0 x 6 请注意 这方程右手边的总和表达式 并不含有f n 0 x displaystyle f n 0 x 项目 任何f n 0 x displaystyle f n 0 x 的贡献 可以与公式 4 的零階項目相合并 将公式 6 代入公式 5 可以得到 D 1 l n 1 f n 0 x l n 0 i n C i f i 0 x D 0 i n C i f i 0 x i n l n 0 l i 0 C i f i 0 x displaystyle D 1 lambda n 1 f n 0 x lambda n 0 sum i neq n C i f i 0 x D 0 sum i neq n C i f i 0 x sum i neq n lambda n 0 lambda i 0 C i f i 0 x 将这方乘式两边都乘以f j 0 x displaystyle f j 0 x 再随著x displaystyle x 积分 利用正交归一性 可以得到 f j 0 x D 1 f n 0 x d x i n l n 0 l i 0 C i f j 0 x f i 0 x l n 0 l j 0 C j displaystyle int f j 0 x D 1 f n 0 x dx sum i neq n lambda n 0 lambda i 0 C i int f j 0 x f i 0 x lambda n 0 lambda j 0 C j 稍加编排 改变下标j displaystyle j 为m displaystyle m 那麼 一阶本征函数修正f n 1 x displaystyle f n 1 x 可以写为 f n 1 x m n f m 0 x l n 0 l m 0 f m 0 y D 1 f n 0 y d y displaystyle f n 1 x sum m neq n frac f m 0 x lambda n 0 lambda m 0 int f m 0 y D 1 f n 0 y dy 参阅 编辑多體微擾理論參考資料 编辑 William E Wiesel Modern Astrodynamics Ohio Aphelion Press 2010 107 ISBN 978 145378 1470 外部链接 编辑摄动方法简介 页面存档备份 存于互联网档案馆 作者Mark H Holmes 第二章 摄动方法简介作者Johan Bystrom Lars Erik Persson 及Fredrik Stromberg 取自 https zh wikipedia org w index php title 摄动理论 amp oldid 69162964, 维基百科,wiki,书籍,书籍,图书馆,