资源描述：
《实验一 常微分方程》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、年级、专业2012级数学与应用数学姓名张旭学号12160011063名单序号16实验时间2014年3月19日使用设备、软件PC,MATLAB注：实验报告的最后一部分是实验小结与收获实验一常微分方程1.分别用Euler法和ode45解下列常微分方程并与解析解比较：(1)解析法：y=dsolve('Dy=x+y','y(0)=1','x')y=2*exp(x)-x-1Euler:function[x,y]=euler(odefun,xspan,y0,h)x=xspan(1):h:xspan(2);y(1)=y0;fori=1:length(x)-1y(
2、k+1)=y(i)+h*feval(odefun,x(i),y(i));endx=x';y=y';endode45:odefun=inline('x+y','x','y');xspan=[0,3];y0=1;h=0.1;[x1,y1]=euler(odefun,xspan,y0,h);[x2,y2]=ode45(odefun,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('eu
3、ler','ode45','dsolve');ode45求得的结果与用解析法求得的结果更接近，故ode45的精度较高，Euler法求得的结果精度较低。9/92014春数学实验实验一常微分方程年级、专业2012级数学与应用数学姓名张旭学号12160011063名单序号16实验时间2014年3月19日使用设备、软件PC,MATLAB注：实验报告的最后一部分是实验小结与收获实验一常微分方程(2)令则原方程等价于方程组：，，不能解析，只能用数值法求解。Euler：function[t,y]=euler2(odefun1,odefun2,tspan,y0,h
4、)t=tspan(1):h:tspan(2);y(1,1)=y0(1);y(2,1)=y0(2);fori=1:length(t)-1k1=odefun1(t(i),y(1,i),y(2,i));k2=odefun2(t(i),y(1,i),y(2,i));y(1,i+1)=y(1,i)+h*d1;y(2,i+1)=y(2,i)+h*d2;endt=t';y=y';endode45：odefun1=inline('0*t1+0*y1+y2');odefun2=inline('-2*y1+0.01*y2^2+sin(t1)');[t1,y1]=eul
5、er2(odefun1,odefun2,[0,5],[0,1],0.1);[t2,y2]=ode45('eu',[0,5],[0,1]);plot(t1,y1(:,1),'o',t2,y2(:,1),'LineWidth',2);xlabel('t轴');ylabel('y轴');legend('euler','ode45');ode45中eu：9/92014春数学实验实验一常微分方程年级、专业2012级数学与应用数学姓名张旭学号12160011063名单序号16实验时间2014年3月19日使用设备、软件PC,MATLAB注：实验报告的最后一部分是
6、实验小结与收获实验一常微分方程functiondy=eu(t,y)dy=zeros(2,1);dy(1)=y(2);dy(2)=-2*y(1)+0.01*y(2)^2+sin(t);ode45求得的结果精度较高，euler法求得的结果在准确值上下波动。2.一通过原点的曲线，它在处的切线斜率等于若上限增为1.58，1.60会发生什么？等价于求解，且的初值问题。解析法：y=dsolve('Dy=2*x+y^2','y(0)=0','x')y=(2^(1/3)*airy(3,-2^(1/3)*x)+2^(1/3)*3^(1/2)*airy(1,-2^(1
7、/3)*x))/(3^(1/2)*airy(0,-2^(1/3)*x)+airy(2,-2^(1/3)*x))ode45：odefun=inline('2*x+y^2');subplot(1,3,1);ode45(odefun,[0,1.57],0);title('01.5之后的
8、斜率增长速度很快，若上限增为1.58，1.609/92014春数学实验实验一常微分方程年级、专业2012级数学与应用数学姓
