微分方程数值.docVIP

微分方程数值解 指导教师：李晓峰 学 院： 应用科学学院 班 级： 信息与计算科学131801 姓 名： 王慧兰 学 号： 201318030120 目录 一、实验名称 错误!未定义书签。 二、实验目的 错误!未定义书签。 三、实验内容 错误!未定义书签。 四、实验要求 错误!未定义书签。 五、实验步骤、源程序及运行截图 错误!未定义书签。 1. Euler法（向前） 错误!未定义书签。 1.1源程序 错误!未定义书签。 1.2实验步骤 错误!未定义书签。 1.3运行截图 错误!未定义书签。 2. Euler法（向后） 错误!未定义书签。 2.1源程序 错误!未定义书签。 2.2实验步骤 错误!未定义书签。 2.3运行截图 错误!未定义书签。 3.梯形迭代法公式 错误!未定义书签。 3.1源程序 错误!未定义书签。 3.2实验步骤 错误!未定义书签。 3.3运行截图 错误!未定义书签。 4.改进的Euler公式 错误!未定义书签。 4.1源程序 10 4.2实验步骤 错误!未定义书签。 4.3运行截图 11 一、实验名称 Euler法、梯形法求解初值问题。 二、实验目的 掌握求解常微分方程初值问题的单步法；比较不同算法所得的数值结果，体会步长缩小对问题解得影响。 三、实验内容 求解初值问题 的数值解。 四、实验要求 编写程序，步长取h=0.1，分别用Euler法、梯形法迭代格式、4阶RK方法求解该问题。整理比较各节点处的数值解、误差、相对误差、体会算法对数值解梯度的改进。 （1）改变步长h=0.2和h=0.5，通过比较各节点处数值解，在不同的步长下，算法优劣结论一致。 （2）按同一算法不同步长整理数据，比较在同一节点按不同步长计算所得的数值解的误差变化趋势，体会缩小步长对数值解的影响。 （3）改进Euler法程序使其适用任意初值问题。 五、源程序、实验核心步骤运行截图 1、Euler法（向前） 1.1源程序 QEuler.m function [h,k,X,Y,P,REn]=Qeuler1(funfcn,x0,y0,b,n,tol) x=x0; h=(b-x)/n; X=zeros(n,1); y=y0; Y=zeros(n,1); k=1; X(k)=x; Y(k)=y; for k=2:n+1 fxy=feval(funfcn,x,y); delta=norm(h*fxy,inf); wucha=tol*max(norm(y,inf),1.0); if delta=wucha x=x+h; y=y+h*fxy; X(k)=x;Y(k)=y; end plot(X,Y,rp) grid label(自变量 X), ylabel(因变量 Y) title(用向前欧拉（Euler）公式计算dy/dx=f(x,y)，y(x0)=y0在[x0,b]上的数值解) end P=[X,Y]; syms dy2, REn=0.5*dy2*h^2; 1.2实验步骤 （1）建立并保存以funfcn.m文件命名的M文件函数 function f = funfcn(x,y) f = 8*x-3*y-7; （2）建立并保存以QEuler.m文件命名的M文件函数。 （3）输入程序： subplot(2,1,1) x0=0;y0=1;b=1-1.e-4; n=100;tol=1.e-4; [h1,k1,x1,Y1,P1,Ren1]=QEuler1(@funfcn,x0,y0,b,n,tol) hold on S1= 8/3*x1-29/9+38/9*exp(-3*x1), plot(x1,S1,b-) title(用向前欧拉公式计算dy/dx=8x-3y-7，y(0)=1在[0,1]上的数值解) legend(n=100时，dy/dx=8x-3y-7，y(0)=1在[0,1]上的数值解, dy/dx=8x-3y-7，y(0)=1在[0,1]上的精确解) hold off jdwuc1=S1-Y1; jwY1=S1-Y1; xwY1=jwY1./S1;k1=1:n;k=[0,k1]; P1=[k,x1,Y1,S1,jwY1,xwY1] subplot(2,1,2) n1=10; [h2,k2,x2,Y2,P2,Ren2]=QEuler1(@funfcn,x0,y0,b,n1,tol) hold on S1 = 8/3*x2-29/9+38/9*exp(-3*x2), plot(x2,S1,b-)