原创力文档ode45的精度较高,
ode45的精度较高,
实验时间 2014年3月19日 使用设备、软件 PC, MATLAB
注：实验报告的最后一部分是实验小结与收获
1.分别用Euler法和ode45解下列常微分方程并与解析解比较:
(1) y 二 x y, y(0) =1,0 ::: x :: 3
解析法:
y=dsolve(Dy=x+y,y(0)=1,x) y = 2*exp(x) - x - 1
Euler: function [x,y]=euler(od efun ,xspa n,y 0,h) x=xspa n(1):h:xspa n( 2);
y(1)=y0;
for i=1:le ngth(x)-1
y(k+1)=y(i)+h*feval(od efun ,x(i),y(i));
end
x=x;
y=y;
end
ode45: odefu n=in li ne(x+y,x,y);
xspa n=[0,3];
y0=1;
h=0.1;
[x1,y1]=euler(od efun,xspan,y0,h); [x2,y2]=od e45(od efun,xspan,y0); x3=0:0.1:3;
y3=2*exp(x3)-x3-1;
Plot(x1,y1,k,x2,y2,ko,x3,y3,k*); xlabel(x 轴);
ylabel(y 轴);
legend(euler,o
文档评论（0）