微分方程数值解部分讲义.docxVIP

2月27日 第一章 常微分方程初值问题的数值解法 引言 本章研究常微分方程初值问题的主要数值解法，包括基本方法和基本理论问 题。 目标在于给出解在一些离散点上的近似值。 考虑常微分方程初值问题 字=/（兀,刃,珀）v兀v X （1.1） ax y仁=〉心0） （L2） 【注】在后面的讨论中，我们总认为这个初值问题的解存在、唯一且连续依赖于 初值条件，即初值问题（1.1）, （1.2）是适定的。 欧拉法（Euler法） 1欧拉法 将解的存在区间N等分，得到N个小区间。任取一个小区间［耳,兀+J,由原 方程? = /（*,刃，结合牛顿-莱布尼兹公式，得到 dx y（￡+i）- 心）=「y\x）dx =「/（x, y）djc 在区间［xn，￡+1 ］上，用/（兀〉）在左端点暫的值来代替/（x?, y（￡））来代替/（%, ,y）, 得到 y（陥）-y（￡）=广 f（xn，y （兀）皿=/（暫，曲“ ））（￡+】-暫）=〃（捡，y （%）） 其中，力二兰二乞为步长。这样， N 》（￡+】）=》‘（￡）+僻（￡*（￡）） （*）可见几何图形 图1.1 这样， ＞＜兀12歹（兀0）+妙（兀0*（兀0）） = %+妙（兀0，%），其中，北=）心（））。记 H 二 yo+E（Xo，）b），即丿（兀1）=必。 y（x2）? y（￡） + hf（X\, y （西））=x + hf（X\, x）,记 y2 = y{ + hf{xx,必）,即 y（x2）^y2. y（x3） ? y（x2）4- hf{x2, y（x2）） ? + hf（x2, y2）,记 y3 = y2^hf（x2,y2）,即 y（x3） ? y3 o 依次类推,我们便在（*）中用几代替y（xn）,并由77的任意性，得到 f ＞；+1 =儿 + 妙K，儿）E = 0,1,2,… bo = ySo） 这就是欧拉公式。 欧拉公式也可以直接由泰勒展开式得到，事实上 y（￡+】）=y（xn）+$ （耳）（兀曲 一 ￡）+*〉‘（￡）（￡+】—暫尸+… =y （兀）+?/（￡）〃+* _/（兀）力2+… 忽略高阶项，就得到 y （捡+】）=丁（￡）+ ）/（￡）〃 =心）+ f（xn,y（xn））h 再在上式中用几替代y（xn）,就得到 =儿+/（￡，几） 这是步长相等的情形。如果不要求步长相等，那么一般情况下，有 = y（%）+/（￡』（￡））（￡+】-暫） 用儿替代畑,就得到 儿+i =儿 + /（￡，儿）（￡+】一?） 图1?2 例1 ?以h = OA为步长，用欧拉法求解初值问题 dx ,y（o）= i 的数值解，并与精确解y⑴二丄（F+2）厂相比较。 【解答】 所建立的欧拉格式为 I ）；+1 =儿 + 0.1（百八—儿）=片 + 0.1（0.1/7^ -几）,2 0丄 2,… 10 = y（o）= 1 利用这个格式来该常微分方程初值问题的数值解，并与精确解相比较，计算结果 如下表: % 0 1 1 0」000 0.9000 0.9094 0.2000 0.8190 0.8351 0.3000 0.7535 0.7742 0.4000 0.7004 0.7239 0.5000 0.6572 0.6823 0.6000 0.6218 0.6476 0.7000 0.5925 0.6182 0.8000 0.5680 0.5931 0.9000 0.5472 0.5712 1.000 0.5291 0.5518 图像: 4 力=0」 /2二0.01 MATLAB 程序： N 二 10; h=l/N; x=0; y=i; y 1 =(xA2+2)*exp(-x)/2; X 二[x]; Y=[y]; Yl=[yl]; disp([x,y,yl]); for k=l:N y=y+(x*exp(?x)?y)*h; x=x+h; y 1 =(xA2+2)*exp(-x)/2; disp([x,y,yl]); X 二[X,x]; Y 二[Y,刃; Yl=[Yl,yl]; end plot（X,Y；r\X,Yl；b,）; legend（EyOp!6a7xl/iE- 14a1）; 以下是欧拉格式的误差估计： 误差的产生的原因： （1） 计算格式本身不能准确描述原来的方程。 （2） 计算机本身引入的误差。（计算机在进行数据运算时总是有截断误差的问 题。） 因此，计算机输岀的是欧拉方程的近似解工，而不是微分方程的精确解 雳一歹（￡）=（氏一儿）+（儿一心）） 可见，为了使计算得到的解工是丿（益）很好的好的近似，我们要求： （D欧拉方法的精确解儿是微分方程精确解y（兀）的很好近似。 ⑵儿是yn好的近似。 问题（1）称为格式的收敛性问题。 问题（2）称为格式的稳定性问题。 2.收敛性研究 欧拉格式收敛的定