数学物理方程第章 数学物理中的近似解法.pdf

第 9 章 数学物理中的近似解法 我们在前面几章讨论数学物理问题时，求得的都是解析解，或说是精确解，这些解法一 般只对较简单的问题和较规则的区域才有效，但实际问题往往很复杂，这时要解出精确解就 很困难，有时甚至不可能，另一方面，在建立数学模型时，我们已作了很多近似，所以求出 的精确解也知识推导出的数学问题的精确解，并非真正实际问题的精确解。因此，我们有必 要研究近似解法，只要使所求得的近似解与精确解之间的误差在规定的范围内，则仍能满足 实际的需要，电子计算机技术的飞速发展，为求得近似解提供了极大的可能。近似解一般可 分为近似解析解和近似数值解。近似解析解是指对解析解在规定误差范围内的近似，近似数 值解是指解在讨论的区域上“足够多”的点处的近似值。 9.1 解析近似解 9.1.1 正则摄动法求解非线性偏微分方程 由于大多数非线性方程无法求得其解析解，这就使我们不得不借助于解析近似法来求其近似 解析解本节将介绍一种求非线性方程近似解析解的正则摄动法。 摄动法的解题思想是：如果非线性方程的非线性项是高阶小量，则作为初步近似，可将其略 去，从而非线性问题便化为了线性问题；求解该线性问题，并将所得的解作为非线性问题的 零级近似；再把原非线性定解问题的解看作它的零级近似解与一个待求的含有小参量的解的 和，代入原定解问题，略去更高阶小量，得到关于摄动解的线性定解问题；求解该线性问题 并将求得的解作为原定解问题的高一级的近似—一级近似；仿此步骤进行下去便可得到各级 近似解。这种求近似解析解的方法被称为摄动法或小参数法，而以上这种带有小参量的定解 问题被称为摄动问题。 我们用一个具体的例子来介绍正则摄支法解题的步骤。 例 1 试在单位圆内求解定解问题 u u ⎧ ∂ ∂ u ε 0(r 1,0 ε 1) ⎪∆ + ≤ r ϕ ⎨ ∂ ∂ （9.1.1 ） ⎪ u(r ,ϕ) r 1 cosϕ ⎩ 解 这是一个含有小参量的非线性问题。当ε 0 时，显然，式（9.1.1 ）变为了线性定 解问题 ⎧∆u 0, r 1 ⎨ （9.1.2 ） u(r,ϕ) r 1 cosϕ ⎩ 由分离变量法即可求得其解，视为零级近似，并记作u0 ，即 u0 r cosϕ （9.1.3 ） 又设 u u0 +εu1 r cosϕ+εu1 （9.1.4 ） 2 将式（9.1.4 ）代入式（9.1.1 ），并略去o(ε ) 项，得 ⎧ 1 ⎪∆u1 r sin 2ϕ ⎨ 2 （9.1.5 ） ⎪ u 0 ⎩ 1 r 1 这是泊松方程的狄氏问题，我们采用下述方法来求解。由式（8.1.5 ）右边函数的形式，我们 猜测方程的解应为 u (r,ϕ) f (r) sin 2ϕ 1 式中，f (r) 为待定函数，将之代入式（9.1.5 ）,得