Euler方法及其改进方法.doc

《常微分方程》课内实验任务书 学生姓名: 张学阳 1009300132 及 学 号： 学 院: 理学院 班 级: 数学101 课程名称： 常微分方程 题 目: Euler方法及其改进方法 指导教师 姓名及职称： 李鹏松 教 授 朱秀丽 讲 师 方 向 实验师 2012年10月17日 《常微分方程》课内实验 实验一 Euler方法及其改进方法 一、实验目的 1．通过用Matlab编程运用Euler方法及其改进方法求解常微分方程初值问题，更进一步掌握常微分方程及其数值解法课程的理论内容，加深对数值解法的理解。 2．熟悉Matlab编程环境。 二、实验学时和类型 本次实验为2学时、设计性实验。 三、实验内容 1.实验题目 运用Euler方法及其改进方法求解常微分方程初值问题，具体问题在课内实验习题中学生任意抽取。 2.实验原理 Euler方法： Euler改进方法： 3.设计思想 在能够获得精确解的题目抽取题目，分别应用Euler方法、Euler改进方法求解数值结果，将所得结果列表或者画图，参照精确结果，对Euler方法、Euler改进方法的计算误差进行分析。 4.参考程序源代码 %fun为目标函数字符串 %x0为自变量初始值。 %y0为fun(x0); %bou=[a,b]自变量区间 %h为步长 fun=1/x^2-y^2; bou=[1,6]; a=bou(1); b=bou(2); x0=1; y0=0; h=0.1; n=ceil((b-a)/h); xx=linspace(a,b,n+1); yy=zeros(1,n+1); lengthx=length(xx); xx(1)=x0;yy(1)=y0; for i=2:n+1 x=xx(i-1);y=yy(i-1); k=eval(fun); yy(i)=yy(i-1)+h*k; end Ys=dsolve(Dy=1/x^2-y^2,y(1)=0,x); for i=1:lengthx x=xx(i); exacty(i)=eval(Ys); end YY=exacty; yend=[xx,yy,YY] p=plot(xx,yend(:,2),k-o,LineWidth,1,... MarkerEdgeColor,k,... MarkerFaceColor,g,... MarkerSize,4); hold on; pp=plot(xx,yend(:,3),r-.+,LineWidth,0.8,... MarkerEdgeColor,r,... MarkerFaceColor,m,... MarkerSize,6); legend([p,pp],Eula,jiequejie); %fun为目标函数字符串 %x0为自变量初始值。 %y0为fun(x0); %bou=[a,b]自变量区间 %h为步长 fun=1/x^2-y^2; bou=[1,6]; a=bou(1); b=bou(2); x0=1; y0=0; h=0.1; n=ceil((b-a)/h); xx=linspace(a,b,n+1); yy=zeros(1,n+1); lengthx=length(xx); xx(1)=x0;yy(1)=y0; for i =2:n+1 x=xx(i-1);y=yy(i-1); k1=eval(fun); x=xx(i);y=yy(i-1)+h*k1; k2=eval(fun); yy(i)=yy(i-1)+h/2*(k1+k2); end Ys=dsolve(Dy=1/x^2-y^2,y(1)=0,x); for i=1:lengthx x=xx(i); exacty(i)=eval(Ys); end YY=exacty; yend=[xx,yy,YY] p=plot(xx,yend(:,2),k-o,LineWidth,1,... MarkerEdgeColor,k,... MarkerFaceColor,g,...