[理学]差分方法稳定性介绍.pptVIP

偏微分方程求精确解 偏微分方程基础理论 三类典型方程差分方法稳定性 作业与上机练习 实验报告形式要求 练习1 练习2 练习3 练习4 取时间步长为0.1、空间步长为0.05,任选一稳定格式计算，并作解的图形。 实验名称、 学生姓名学号 实验目的、内容、 算法、程序、 计算结果及分析 （可提交word, pdf电子版） Coruant-Friedrichs-Lewy条件 Lax等价定理 第四节 一阶线性方程（组）的几个差分格式 迎风格式(upwind scheme) Lax-Friedrichs格式 Lax-Wendroff格式 蛙跳格式(leap-frog scheme) Fourier积分 Fourier方法判断稳定性 定理1 差分方程（迎风格式）按2-模稳定的充分必要条件是：存在M0, 对 有 （Von Neumann条件） 证明：（1）和（2）的等价性 特别取 ，有 注：此定理的证明结论可用于其他类似的差分格式。 方程组情形 定义：若A的特征值全是实数，且有p个线性无关的特征向量，则称方程组为双曲型方程组。 如： 方程组情形 Von Neumann条件 Lax-Friedrichs格式 Lax-Wendroff格式 蛙跳格式(leap-frog scheme) 问：迎风格式如何？ 迎风格式（常系数） 一阶线性变系数方程 Lax-Friedrichs格式 例 一阶线性变系数方程组 Lax-Friedrichs格式 Lax-Wendroff格式 守恒律方程 Lax-Friedrichs格式 例 守恒律方程组 Lax-Friedrichs格式 Richtmyer二步格式 参考书: 1. 数学物理方程与特殊函数, 杨奇林编, 清华大学出版社 2.数学物理方程,欧维义编,吉林大学出版社 3. 谷超豪等，数学物理方程（第二版），高等教育出版社，2002 4. Practical Partial Differential Equation(实用偏微分方程,第一版)), R. Haberman著, 2005.5 参考书: 1.偏微分方程, 周蜀林编, 北京大学出版社， 2005 2.偏微分方程, 南京大学数学系, 1979 偏微分方程的差分方法 参考书: 1.李瑞霞、何志庆,微分方程数值方法，华东理工大学出版社，2005.12 2.陆金甫、关治，偏微分方程数值解法（第二版),清华大学出版社，2004. 3. 徐翠薇，计算方法引论，高等教育出版社，1985 4.陆君安等，偏微分方程的MATLAB解法，武汉大学出版社，2001 参考书: 1.李瑞霞、何志庆,微分方程数值方法，华东理工大学出版社，2005.12 2. 徐翠薇，计算方法引论，高等教育出版社，1985 一、矩形区域上的Laplace方程 网格划分 五点差分格式 一般区域上的调和方程求解 内点（内部节点）：D内的网格节点。记为 非正则内点：相邻四个网格点不全是内点。 边界点（边界节点）：网格线与边界的交点。 正则内点：有相邻四个内点 可用五点公式建立方程 可用临近边界点的值近似 一般区域上的调和方程求解 对区域内的节点按顺序编号，并用U(i)表示与第i个节点相邻的所有节点的集合（不包括第i个节点）。于是差分格式可写成 极值原理：若对所有i, 都有 ,且不全相等，则 不会在内点取得这组数的正的最大值（负的最小值）。 证明： 反证法。设在节点i处取正最大值M。由区域的连通性，必存在某个节点k处取M，且 中至少有一个 矛盾！！ 定理1 差分方程组（1）有唯一解。 定理1 差分方程组（1）有唯一解。 定理2 设 则 推论 设 则 定理3 设 ,则存在常数C，使得 二、一维热传导方程求解 差商格式 算法的收敛性和稳定性 收敛性证明 稳定性证明 热传导方程隐式差分格式 无须对 加条件, 隐格式是无条件稳定的! 第三节 波动方程求解 差分方程 初始速度的离散 差分格式