常微数值解实验报告1(模板).docVIP

实验目的： 1．Matlab编程运用Euler方法及其改进方法求解常微分方程初值问题，更进一步掌握常微分方程及其数值解法课程的理论内容，加深对数值解法的理解。 2．atlab编程环境。 实验内容： 运用Euler方法及其改进方法求解常微分方程初值问题。 实验题目：，。 实验原理： 1.Euler方法： 实验程序： %fun为目标函数字符串 %x0为自变量初始值。 %y0为fun(x0); %bou=[a,b]自变量区间 %h为步长 fun=y/x+x^2; bou=[4,100]; a=bou(1); b=bou(2); x0=1; y0=4; h=4; 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=y/x+x^2,y(1)=4,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,jingquejie); 2. Euler改进方法： 实验程序： %fun为目标函数字符串 %x0为自变量初始值。 %y0为fun(x0); %bou=[a,b]自变量区间 %h为步长 fun=y/x+x^2; bou=[4,100]; a=bou(1); b=bou(2); x0=1; y0=4; h=4; 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=y/x+x^2,y(1)=4,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],gjeula,jingquejie); 实验结果： 表1 Euler方法与精确解的比较 节点 欧拉法 精确解 1 4 4 8 24 284 12 292 906 16 965.3333333 2104 20 2230.666667 4070 24 4276.8 6996