有限差分法基础课件.pptVIP

有限差分法基础;数值离散概述;有限差分法求解示例;离散网格点;令 得到差分格式;其中;差分和逼近误差; 差分的三种形式（一阶）：; 由导数（微商）和差商的定义可知，当自变量的差分（增量）趋近于零时，就可以由差商得到导数。因此在数值计算中常用差商近似代替导数。;差分和逼近误差;差分和逼近误差;差分和逼近误差;差分和逼近误差;差分和逼近误差;差分和逼近误差;差分和逼近误差（二阶导数）; 二阶中心差分：;差分方程的建立过程;模型方程; Poisson方程：;时间一维空间一??的例子; 1.划分网格;; 时间导数用一阶向前差商近似代替：; 差分方程和其定解条件一起，称为相应微分方程 问题的差分格式。上述初值问题的差分格式可改写为：; 显式有限差分模板：; 时间推进：;第一个例子的求解结果;精确解：; 例 考虑长度为1的均匀直杆，其表面是绝热的，而且杆截面足够细，可; 时间导数用一阶向前差商近似代替：; 显式有限差分模板：;; 如仍取 而为缩短计算时间，时间步长 取 ，则最终的差分方程：;差分法的基本理论; 上式中的 就是差分方程与微分方程的差别，称之为截断误 差。显然 与 、 成正比，一般情况下，当步长趋向零时,有限差分方程的截断误差是趋向于零的，则称有限差分方程与相应的偏微分方程是相容的。 一个可用的偏微分方程的差分表达式必须是相容的。否则在 、 趋近零时，差分方程不能趋于原微分方程，差分方程的解就不能代表微分方程的解，差分求解就失去了意义! ; 2.收敛性 收敛性研究的是差分方程的解与微分方程的解之间的差别问 题。如果在求解区域中的任一离散点 上，当网格步长 、 趋于零时，有限差分方程的解趋近于所近似的微分方程解，则称有 限差分方程的解是收敛的。 ; 3.稳定性 稳定性讨论的是???分解的误差在计算过程中的发展问题。在 数值解中，引进误差是不可避免的，电子计算机也有舍入误差， 因此实际算得的有限差分方程的解是近似解。这种误差是要向其 他方向传播的，如果计算中引入的误差在以后逐层计算过程中影 响逐渐消失或者保持有界，则称差分方程是稳定的。否则就是不 稳定的。; 上式中 为差分方程的精确解,如果令 为差分方程的近似数值 解，之间的误差为 ???同样，近似数值解也满足同样的方程：; 4.Lax等价定理 　 对于一个适定的线性初值问题，如果有限差分近似是相容的，则稳定性是收敛性的充分和必要条件。这是有限差分方法最基本的定律。 　　 适用条件： 1）偏微分方程的解存在、唯一且连续地依赖于初值； 2）该定理只适用于线性问题，对非线性此定理至今未得到证明。 重要的实际意义：一般情况下，证明有限差分方程的解收敛于它所近似的偏微分方程的解比较困难。而证明有限差分方程的稳定性和相容性相对来说比较容易。根据该定理只要证明有限差分方程是相容的、稳定的，就保证了收敛性。;几种差分格式介绍;几种差分格式介绍;几种差分格式介绍;;几种差分格式介绍;几种差分格式介绍