第5章偏微分方程.pptVIP

绕过上节显式算法的条件稳定性的一个方法是采用相应的隐式格式。 隐式差分法 获取隐式格式的一个直接的方法是，只需将空间的二阶导数用新时刻上的三点差分公式来近似代替，其它不变,这样我们得到一个修正的迭代格式 引入算符 L 将隐式公式写为 因为未知量 фn+1 出现在方程两边，所以我们可以通过解出 фn+1 来得到迭代格式，即求解 注意， L 的矩阵形式 这个格式在 Δt 的最低阶上和显式格式等价，但是这个格式比显式格式要好得多，好在它可以用大的时间步长，这是因为该隐式差分格式是无条件稳定的。 对隐式格式而言，每个时间步的计算都需要将矩阵 (1+ Δt L) 的逆矩阵作用在 фn 上。由于逆矩阵本身与时间无关，我们可以只在计算开始时求它一次，然后在所有时刻都用。 将显式公式 利用算符 L 改进的算法 写为 将隐式公式和显式公式相加 解出得 最后我们需要说明的是，用两个时刻上的三点差分公式的不同组合来近似代替空间二阶导数，还可以得到其它改进的隐式迭代格式，例如， 这个差分格式被称为Crank–Nicolson方法，是求解含时薛定谔方程的首选。它的精度比前面的直接隐式差分格式好得多（直接隐式差分格式对时间的精度是 O(Δt )，而Crank–Nicolson差分格式则是 O(Δt 2)）。 含时薛定谔方程 其中 V 是位势函数。 为了方便起见，令 ?=2m=1，将上式写为如下抛物线方程形式 其中 H 算子为 可以证明这个递归关系是幺正的，即保证波函数的模方在全空间的积分不随时间变化，这正是量子力学所要求的。 应用下面的差分格式 5.3 双曲型方程 研究两端固定的弦的振动 假设初始位移 为 ф1(x) 和初始速度 ψ1(x) 将二阶导数写为差分形式 当 c1 时，解是稳定的，当 c 1 时，解是不稳定的。 得到波动方程的差分格式 其中 利用初始条件可以先求出两行数值作为公式计算的启动值，设初始位移 为 ф1(x) 和初始速度 ψ1(x) ，则 显式公式表明，需要两行的已知的数据才能求出下一行的数值 Matlab偏微分方程工具箱简介 MATLAB提供了两种方法解决偏微分方程问题： 一个是 pdfpe 函数，可以直接求解一般的偏微分方程 另一个是matlab 的偏微分工具箱（PDE toolbox）可以比较规范的求解各种常见的二阶偏微分方程 以椭圆型方程为例，演示偏微分工具箱的用法 平面有界区域 D上的椭圆型偏微分方程边值问题 a、c、f 均为D上的关于 (x, y) 的已知函数 未知函数为 u(x, y), 它的边界条件为三类 （1）Direchlet条件 （2）Neumann条件 （3）混合边条件在边界上部分为Direchlet条件，另外部分为Neumann条件 例子 * 单击此处编辑母版标题样式 单击此处编辑母版标题样式 第五章 偏微分方程 本章内容 椭圆型方程 1 2 抛物型方程 3 双曲型方程 5.0 概论 任何一种随空间变化或随时空变化变化的物理现象都需要用偏微分方程描述 大部分物理上重要的偏微分方程是二阶的，它们可以分为 因变量用它们在自变量的许多离散的格点上之值来描述， 通过适当的离散化，偏微分方程就化为一大组差分方程。 双曲型方程 抛物型方程 椭圆型方程 泊松方程和定态薛定谔方程 波动方程 扩散方程和含时薛定谔方程 5.1 椭圆型方程 我们将考虑二维空间 (x， y) 内的场量 ф 的椭圆型边值问题，方程为 取固定边界条件，即在 (x， y) 平面内某一根很大的闭合曲线上（为方便起见取为单位正方形）规定了 ф 值，边值问题就是要用方程求出单位正方形内每处的 ф。 空间离散化 首先定义一个网格，覆盖 (x， y) 平面内的单 位正方形。为方便起见，我们取格子间隔 h (步长)是均匀的，并且在两个方向上相等，使得单位正方形被 (N+1)ⅹ (N+1) 个格点覆盖。这些格点 用指标 i，j 编号，其中 i， j = 0，1，2，3，...，N。 对每个方向上的二阶微商应用三点差分近似，方程近似为 其中 A 是出现在线性方程组中的矩阵， ф 表示内点上的 ф 值（不包含边界点），非齐次项 B 包含内点上的 S 以及边界上的 ф 值。 上式等价于关于区域的内点上的未知 ф 值的一个线性方程组，写为矩阵形式 ф 可以通过矩阵求逆来获得 一维问题的差分格式 一维椭圆型方程 差分格式为 当点数不太多时，我们可以利用矩阵求逆的方法来求解 具体的为 例子 精确解为 二维问题的差分格式 其中 取特殊的网格 h=1/3, 如右图，求差分方程的矩阵形式 差分方程的矩阵形式为 其中 解得 松弛法 首先将上式改写为“解出”фi 的形式 以一维问题为例，其差分方程为 当格点取得非常密时，选择矩阵求逆的方法来求解计算量是非常大的，需要寻求