数值分析第9章-1.pptVIP

* * 结束 第九章 常微分方程初值问题的数值解法 科学研究和工程技术中的问题往往归结为求某个常微分方程的定解问题. 常微分方程的理论指出，很多方程的定解虽然存在，但可能十分复杂难于计算，也可能不能用简单的初等函数表示,因此常求其能满足精度要求的近似解. 常微分方程的数值解法常用来求近似解，由于它提供的算法能通过计算机便捷地实现，因此近年来得到迅速的发展和广泛的应用常微分方程数值解法的特点是：对求解区间进行剖分，然后把微分方程离散成在节点上的近似公式或近似方程，最后结合定解条件求出近似解.因此数值解法得到的近似解是一个离散的函数表. §9.1 引 言 * 结束 9.1.1 基本知识复习 一阶常微方程的初值问题： 其中f (x,y)是已知函数，(2)是定解条件，常微分方程的理论指出：当f (x,y)定义在区域G=(a≤x≤b,｜y｜＜∞)，若存在正的常数L，使： ｜f (x,y1) – f (x,y2) ｜≤L｜y1- y2| 对任意的x∈[a,b]和任意的y1, y2成立，则称f (x,y)对y满足李普希兹(Lipschitz)条件,L称李普希兹常数. 若f (x,y)在区域G连续，关于y满足李普希兹条件，则上述一阶常微分方程的初值问题的解存在且唯一.我们以下的讨论，都在满足上述条件下进行. * 结束 9.1.2 一阶常微分方程组和高阶常微分方程 一阶常微分方程组常表述为： 写成向量形式为 * 结束 还有高阶常微分方程定解问题，如二阶定解问题： 可化作一阶常微方程组的定解问题： 本章着重讨论一阶常微方程初值问题的各种数值解法.这些解法都可以写成向量形式而用于一阶常微分方程组的初值问题.这样，也就解决了高阶方程的定解问题. * 结束 §9.2 欧拉(Euler)方法 对问题 最简单而直观的方法是欧拉方法.欧拉方法在精度要求不高时，仍不失为一实用方法.同时通过欧拉方法的讨论，容易弄清常微方程初值问题数值解法的一些基本概念和构造方法的思路. 9.2.1 欧拉方法的导出 把区间［a,b］分为n个小区间，取步长h=(b-a)/ n ，节点xi=x0+ih,i=0,1,2,…，n，其中x0=a，又设y (x)为上述问题的解. * 结束 把(9.1)式两端积分，积分区间取［xi,xi+1］ 记g (x)=f (x, y (x))对(9.3) 用左矩形公式(8.3) 得到： 舍去Ri，并令yi=y(xi),得到： 这就是欧拉公式. * 结束 欧拉公式还可由泰勒展开得到.假定y (x)二阶连续可导，把 y(xi+1)在xi点展开： 舍去h2项，并令yi= y(xi) ，注意到y’(xi)=f (xi, yi),也得到 欧拉公式有明显的几何意义，如图9.1，过点(x0, y0)的曲线是解y(x).欧拉方法是在(x0, y0)作y(x)的切线，它与直线x=x1交于 (x1, y1)，过(x1, y1)作过此点的积分曲线的切线，又与x=x2交于(x2, y2)，…如此下去，得到一条折线，欧拉方法就是用这条折线近似地代替曲线y(x)，故欧拉方法有时也称欧拉折线法. 9.2.2 欧拉隐式公式和欧拉中点公式 对(9.3)式的右端用右矩形公式代替积分，可得到 * 结束 推出： 不难看出，(9.5)式中公式两端都含有yi+1,一般情况不能由yi的值计算yi+1,而需要解方程，故称为欧拉隐式公式. 推出： 称为欧拉中点公式.容易看出,中点公式计算yi+1时,不仅需要yi的值，还需要yi-1的值. 把(9.3)式积分区间改为[xi-1,xi+1]，并对右端用中矩形公式代替积分，可得到： * 结束 凡是从已知(或已算出)的y0,y1,…,yi能直接从公式算出yi+1的公式称显式，否则称隐式公式(9.4),(9.6)都是显式，(9.5)是隐式在计算yi+1时，只需要yi的值，则称公式为单步法；若除yi之外还需要以前的yi-1等多个值，则称多步法公式.(9.4)、(9.5)是单步法，(9.6)是多步法，确切说是二步法. 9.2.3 局部截断误差和方法的阶 以上三个公式都有误差.假定yi,yi-1的值没有误差，用公式计算yi+1时产生的误差称为局部截断误差，也称为单步截断误差，假定yi+1以前的y值都没有误差，局部截断误差表示为： 对欧拉公式： * 结束 对欧拉隐式公式： 对欧拉中点公式的局部截断误差： 凡是局部截断误差为O(hp+1)的方法称为p阶方法，这样欧拉方法和欧拉隐式方法称为一阶方法，而欧拉中点方法称为二阶方法.一般讲高阶方法比低阶方法要精确. * 结束 9.2.4 梯形公式及其预估校正法 对积分式(9.3)的右端用梯形积分公式近似代替，得 把y(xi+1)记为yi+1，并略去余项-h3y(ξ)/12得到