常微分方程的数值解法教学.pptxVIP

常微分方程的数值解法PPT教学课件会计学第一节 引言第1页/共35页●在常微分方程课程里面讨论的是一些典型方程求解解析解的基本方法。●然而在生产实践和科学研究中遇到的微分方程往往比较复杂，在很多情况下，不能给出解的解析表达式；有时候即时能用解析表达式来表示，又因为计算量太大而不实用，有时候一些是已经有了求解的基本方法的典型方程，但实际使用时也是有困难的。●以上情况说明用求解解析解的基本方法来求微分方程的解往往是不适宜的，甚至很难办到。●实际问题中，对于求解微分方程，一般只要求得到解的若干个点上的近似值或者解的便于计算的近似表达式。●本章研究微分方程的数值解法，而且着重讨论微分方程中最简单的一类问题——一阶方程的初值问题。第一节 引言第2页/共35页1、一阶方程的初值问题假定上式在区间[a, b]上存在唯一且足够 光滑的解y(x)。●所谓数值解法就是寻求解y(x)在一系列离散点，也称为节点处的值：要计算出解函数 y(x) 在一系列节点 a = x0 x1… xn= b处的近似值第一节 引言第3页/共35页●节点间距，即步长为： 通常采用等距节点，即hi = h (常数)●等间距节点●在这些节点上采用离散化方法（通常用数值积分、微分、泰勒展开等）将上述初值问题化成关于离散变量的相应问题。把这个相应问题的解yn作为y(xn)的近似值。这样求得的yn就是上述初值问题在节点xn上的数值解。一般说来，不同的离散化导致不同的方法。记为第二节 欧拉方法第4页/共35页一、欧拉法Euler1、向前差商近似导数第二节 欧拉方法第5页/共35页2、举例例1 用欧拉法求初值问题当h = 0.02时在区间[0, 0.10]上的数值解。解：根据欧拉公式可以得到：此外，可以得到方程的真解：nxnyny(xn)?n = y(xn) - yn001.00001.0000010.020.98200.98250.000520.040.96500.96600.000530.060.94890.95030.001440.080.93360.93540.001850.100.91920.9230.0021第二节 欧拉方法第6页/共35页求解过程如下：第二节 欧拉方法第7页/共35页根据已知条件：曲线y(x)上的点（x0,y0）及该点处曲线的导数f(x0,y0)，则可以得到过该点的直线：3、欧拉方法的几何意义该直线与x＝x1的交点P1，则P1的纵坐标y1为：就用y1作为y(x1)的近似值逐次进行后可以得到一条折线P0P1…Pn，该折线看作是初值问题的积分曲线的近似，因此欧拉方法也称为欧拉折线法…A第二节 欧拉方法第8页/共35页从上述几何意义上得知，由Euler法所得的折线明显偏离了积分曲线，可见此方法非常粗糙即误差太大。4、欧拉法的局部截断误差（1）截断误差定义在假设 yi = y(xi)，即第 i 步计算是精确的前提下，考虑的截断误差 Ri+1 = y(xi+1) ? yi+1 ，称为局部截断误差如图所示：APi+1即为欧拉方法在xi+1点的截断误差（2）如果某种方法的局部截断误差是则称该方法具有p阶精度第二节 欧拉方法第9页/共35页（3）则截断误差的大小？写出y(xn+1)的泰勒展开式：由欧拉方法可以得到：则上面两个公式相减得到：具有1阶精度第二节 欧拉方法第10页/共35页二、改进的欧拉法一阶方程的初值问题与如下积分方程是等价的：当x = x1时 可以借助于数值积分，求y(x1)的值 1、用矩形公式第二节 欧拉方法第11页/共35页可以推导出：用矩形法计算右端的积分与用欧拉法计出的结果完全相同 2、用梯形公式 则可以推导出：第二节 欧拉方法第12页/共35页梯形公式的截断误差：梯形公式具有二阶精度，比欧拉方法有了进步第二节 欧拉方法第13页/共35页 和欧拉公式相比较，梯形公式在计算yi+1时候也只用到前一步的值yi，但是若yi已知，将yi带入公式求解时候，一般不能直接得到yi+1，而需要通过其他方法（比如迭代法）求解，所以梯形公式被称为隐式公式。3、改进的欧拉方法 梯形公式是隐式的，一般用迭代法求解，计算量较大。实际中常将欧拉公式和梯形公式联合使用，先用欧拉公式得出一个y(xi+1)的近似值 称为预估值，然后对预估值使用梯形公式对它进行精确化，得到较为精确的近似值yi+1，称之为校正值，计算公式为：这样的预估校正系统称为改进的欧拉方法。第二节 欧拉方法第14页/共35页为了便于编写程序，常将上面的公式改写为如下式：第二节 欧拉方法第15页/共35页4、举例 P90，例题5-1在区间[0, 1.5]上，取h = 0.1，求解。解：（1）用欧拉法计算公式如下： （2）用改进欧拉法计算公式如下： 第16页/共35页xn欧拉法yn迭代一次改进欧拉法yn准确解01