计算方法--常微分方程求解实验.docVIP

实验五 常微分方程求解实验 实验目的 通过本实验学会对给定初值我呢他，用欧拉法、改进欧拉法、四阶龙格-库塔法求数值解和误差，并比较优缺点.对给定刚性微分方程，求其数值解，并与精确解比较，分析计算结果. 实验题目 解初值问题各种方法比较 实验题目：给定初值问题 精确解为，按 欧拉法，步长 改进欧拉法，步长 四阶标准龙格-库塔法，步长 求在节点处的数值解及误差，比较各方法的优缺点. 刚性方程计算 实验题目：给定刚性微分方程 其精确解为.任选取一显示方法，取不同的步长求解，并分析计算结果. 实验原理 1.欧拉格式 由数值微分的向前差商公式可以解决初值问题（6.1）中的导数的数值计算问题： 由此可得 （6.1）实际上给出 于是有 再由得 （6.2） 递推公式（6.2）称为欧拉格式。 2.改进欧拉格式 先对欧拉格式（6.2）对进行计算，并将结果记为，再代入（6.7）可得“预报-校正”形式的差分格式： 公式（6.8）称为改进欧拉格式。 3.四阶经典龙格-库塔格式： 实验内容 1. 解初值问题各种方法比较 实验程序： function shiyan51 h=0.1; dyfun=inline(y./x+x*exp(x)); [x,y1]=maeuler1(dyfun,[1,2],0,h); [x,y2]=maeuler(dyfun,[1,2],0,h); [x,y3]=marunge4(dyfun,[1,2],0,h); y=x.*(exp(x)-exp(1)); err1=abs(y-y1); err2=abs(y-y2); err3=abs(y-y3); [x,err1,err2,err3] 3.刚性方程计算 实验程序： function shiyan53 clc %format long fun=inline(-600*y+1199.8*exp(-0.1*x)-600); tic; [x,y1]=marunge4(fun,[0,5],2,0.001); yb1=exp(-600*x)+2*exp(-0.1*x)-1; t1=toc err1=abs(yb1-y1); s1=sum(err1.^2) % tic; [x,y2]=marunge4(fun,[0,5],2,0.0001); Yb2=exp(-600*x)+2*exp(-0.1*x)-1; t2=toc err2=abs(yb2-y2); s2=sum(err2.^2) 实验结果 1. 解初值问题各种方法比较的实验结果 1.00000000000000 0 0 0 1.10000000000000 0.04264443219022 0.00097356580312 0.00000092300897 1.20000000000000 0.09520401466700 0.00199790937570 0.00000178250161 1.30000000000000 0.15865732017504 0.00306315484526 0.00000259633614 1.40000000000000 0.23411778154437 0.00415909973975 0.00000337650650 1.50000000000000 0.32284614010794 0.00527488906625 0.00000413124347 1.60000000000000 0.42626522205571 0.00639872308444 0.00000486627011 1.70000000000000 0.54597705499551 0.00751757761384 0.00000558557848 1.80000000000000 0.68378242658891 0.00861692278559 0.00000629192402 1.90000000000000 0.84170305352887 0.00968043001901 0.00000698714712 2.00000000000000 1.02200658351062 0.01068965912814 0.00000767238581 3.刚性方程计算的实验结果 t1 =8.9953 s1 =1.3221e-006 t2 =100.9231 s2 =5.3715e-014 实验结