太理水利工程计算及设计第七章-常微分方程初值问题的数值方法.pptVIP

第一节 有限差分解法 一、有限差分的基本概念 二、欧拉法及其简单改进 三、欧拉法的误差分析 * 第七章 常微分方程初值问题的数值方法 本章内容 许多实际问题的数学模型是微分方程的定解问题, 而能用解析方法求出精确解的微分方程为数不多,有的方程即使有解析解,也可能由于解的表达式非常复杂而不易计算,因此有必要研究微分方程的数值解法。 问题的描述和基本概念 一阶常微分方程初值问题的形式 初值问题：给定起始点（t=t0）上未知量y的值而确定其未来（tt0）值的问题叫初值问题。 边值问题：未知量y的值给定在多个点上的问题叫边值问题。 解的存在性与唯一性 设有初值问题 若f(y,t)在区域G：0≤t≤T，︱y︱∞内连续，且存在常数L，使得对于0≤t≤T及所有的y、y*均有 成立，初值问题（8-4）的解存在且唯一。 （8-4） 里普希兹条件 （8-5） 初值问题数值解法的提法 所谓初值问题的数值解法，就是寻求问题的解y(t)在一系列节点t0t1t2…tn=T上的近似值y0,y1,y2, …,yn的近似算法，相邻节点的距离hj=tj-tj-1称为步长。 数值方法给出有限个近似值yj（j=1，2，…，n），利用这n个值，通过插值，可以得到整个区间上的y(t) 的近似值。 单步法与多步法 对初值问题（8-4）的数值解，若y(tj+1)的近似值yj+1是用其前面的K个已知值yj,yj-1, …,y j-(K-1)计算的，就称该方法为K步法。如K=1就是单步法，称K1的方法为多步法。 如果计算yj+1要用到yj+1本身，则称之为隐式法，否则为显式法。 隐式法与显式法 第一节 有限差分解法 有限差分法是求解微分方程的基本数值方法之一，是以小而有限的间距上变量的变化（差分）的比值（差商）代替某点的导数（或偏导数），即令 从而使连续的初值问题离散化为节点上的代数方程。由此建立的代数方程称为差分方程。 差分近似式可由泰勒级数得到。 第一节 有限差分解法 1. 向前差分近似式 解 得 舍去O(△t)项，得到导数的向前差分近似式 （8-7） 截断误差为O(△t） （8-6） 第一节 有限差分解法 2. 向后差分近似式 解 ，并舍去O(△t)项，得到导数的向后差分近似式 （8-9） 截断误差为O(△t） （8-8） 第一节 有限差分解法 （8-11） 截断误差为O(△t2） （8-10） 两式相减得 3. 中心差分近似式 （8-6） （8-8） 第一节 有限差分解法 4. 二阶导数的差分近似式 （8-12） 截断误差为O((△t)2） 两式相加得 （8-6） （8-8） 第一节 有限差分解法 如记变量t的等距离散值tj=t0+jh，其中h=tj-tj-1, j=1,2, …,n， 则上述各差分近似可简记为 向前差分 向后差分 中心差分 二阶差分 这里yj=y(tj)。 （8-13） （8-14） （8-15） （8-16） 5. 变量t为等距离散值的差分近似式 第一节 有限差分解法 第一节 有限差分解法 （1） 相容性 是指当h=△t→0时，差分方程应逼近微分方程。相容性保证了差分方程的截断误差与步长h一同趋于零。截断误差趋于零愈快，差分方程的准确度愈高，因而数值解的准确性愈好。 6. 差分近似式的选择以及由此构造的差分解法应满足相容性、收敛性及稳定性的要求。 （2）收敛性 是指当h=△t→0时，差分方程的数值解应趋于微分方程的准确解。 （3）稳定性 是指对于固定的步长h，当t充分大时，数值解的误差能限制在一定范围内，不致使误差淹没解。 第一节 有限差分解法 欧拉法是基于差分法的最简单的显式单步法，尽管其精度较低，但算法简单，起源最早，也是其它较复杂数值解法的一个很好的引导。 一阶常微分方程初值问题的形式 用向前差分式 1. 欧拉法 第一节 有限差分解法 代替导数，整理得 （8-18） 计算式（8-18）就是欧拉法公式。 欧拉法的几何意义：欧拉法常称为折线法。 第一节 有限差分解法 Y=y(t) t0 t1 t2 t3 p0 p1 p2 p3 2.向后欧拉公式 以t=tj+1处的向后差分式 取代初值问题（8-4）中的导数，则有 从而得到向后欧拉公式 （8-20） 第一节 有限差分解法 3.中点欧拉公式 以t=tj处的中心差分式 取代初值问题（8-4）中的导数，就可得到中点欧拉公式 （8-21） 第一节 有限差分解法 4. 梯形欧拉公式