第八章常微分方程的数值解_第1页
第八章常微分方程的数值解_第2页
第八章常微分方程的数值解_第3页
第八章常微分方程的数值解_第4页
第八章常微分方程的数值解_第5页
已阅读5页,还剩17页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

第八章常微分方程的数值解第八章 常微分方程的数值解引言简单的数值方法欧拉方法梯形方法在高等数学中我们见过以下常微分方程:8.1引言(1),(2)式称为初值问题,(3)式称为边值问题。(其中L为Lipschitz常数)则初值问题(1)存在唯一的连续解。

考虑一阶常微分方程初值问题其中,y=y(x)是未知函数,y(x0)=y0

是初值条件,而f(x,y)是给定的二元函数.由常微分方程理论知,若f(x)在x[a,b]连续且f满足对y的Lipschitz条件:常微分方程的数值解法有单步法和多步法之分:单步法:在计算yn+1时只用到前一点yn

的值;多步法:计算yn+1时不仅利用yn,还要利用yn-1,yn-2,…….,一般k步法要用到yn,yn-1,yn-2,…….,yn-k+1。求问题(1)的数值解,就是要寻找解函数在一系列离散节点x1<x2<……<xn<xn+1

上的近似值y1,y

2,…,yn

。为了计算方便,可取xn=x0+nh,(n=0,1,2,…),h称为步长。8.2简单的数值方法一、欧拉(Euler)方法在x=x0处,用差商代替导数:由得同理,在x=xn

处,用差商代替导数:由得若记则上式可记为此即为求解初值问题的Euler方法,又称显式Euler方法。Pn+1yOxx0x1x2xnP0P1P2Pny=y(x)xn+1Euler方法的几何意义:(Euler折线法)例:用Euler方法求解常微分方程初值问题并将数值解和该问题的解析解比较。解:Euler方法的具体格式:

xn

y(xn) yn

yn-y(xn) 0.0 0 0 00.2 0.1923 0.2000 0.00770.4 0.3448 0.3840 0.03920.6 0.4412 0.5170 0.07580.8 0.4878 0.5824 0.0946 1.0 0.5000 0.5924 0.09241.2 0.4918 0.5705 0.07871.4 0.4730 0.5354 0.0624取h=0.2,xn=nh,(n=0,1,2…,15),f(x,y)=y/x–2y2

计算中取f(0,0)=1.计算结果如下:xn

y(xn) yn

yn-y(xn)1.6 0.4494 0.4972 0.0478 1.8 0.4245 0.4605 0.03592.0 0.4000 0.4268 0.02682.2 0.3767 0.3966 0.0199 2.4 0.3550 0.3698 0.01472.6 0.3351 0.3459 0.01082.8 0.3167 0.3246 0.0079 3.0 0.3000 0.3057 0.0057由表中数据可以看到,微分方程初值问题的数值解和解析解的误差一般在小数点后第二位或第三位小数上,这说明Euler方法的精度是比较差的。O:数值解;—

:准确解数值解和解析解的图示比较如下:若直接对y`=f(x,y)在[xn,xn+1]积分,利用数值积分中的左矩形公式:此即为Euler公式。设y(xn)=

yn,则得若用右矩形公式:得上式称后退的Euler方法,又称隐式Euler方法。可用迭代法求解:初值:迭代:k=0,1,……因故当hL<1时,迭代法收敛。二、梯形方法由利用梯形求积公式:得上式称梯形方法,是一种隐式方法。用迭代法求解:初值:迭代:k=0,1,……因故当hL/2<1时,迭代法收敛。由以上分析可以看出,隐式方法的计算比显式方法复杂,需要用迭代法求解非线性方程才能得出计算结果。可采用将显式Euler格式与梯形格式结合使用的方法来避免求解非线性方程。记再用梯形格式计算:--预测--校正上面两式统称预测-校正法,又称改进的Euler方法。三、单步法的局部截断误差和精度单步法的一般形式为:(

与f有关)显式单步法形式为:整体截断误差:从x0开始,考虑每一步产生的误差,直到xn,则有误差称为数值方法在节点xn处的整体截断误差。但en不易分析和计算,故只考虑从xn到xn+1的局部情况。定义:设y(x)是初值问题(1)的精确解,则称为显式单步法在节点xn+1处的局部截断误差。若存在最大整数p使局部截断误差满足则称显式单步法具有p阶精度或称p阶方法。注:将Tn+1表达式各项在xn处作Taylor展开,可得具体表达式。Euler方法的局部截断误差:故Tn+1=O(h2),p=1,(设yn=y(xn))其中称局部截断误差主项。即Eul

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论