,也能使第二种情况水平速度达到极限值,这与前面解析分析的结果一致。注意由于计算的时间范围较大,抛体上升的高度较小,所以上升段的曲线不明显,如果缩短计算的时间,就可以更明显地表现上升段的曲线。4.思考题⑴将解微分方程的时间范围缩短为0~1s,适当改变坐标轴的设定范围,重新画图,结果有什么不同?⑵画出抛体的垂直速度随时间变化的图形,以观察它达到极限值的情况⑶et改为plot,结果有什么不同?5.参考程序主程序的文件名是znxpyd.mm=1;b=[0,0.2,0.2];p=[0,0,1];px=[4.6;4.5;4.5];py=[3.5;1.8;0.4];strdd{1}='无阻尼';strdd{2}='阻力正比于v';strdd{3}='阻力正比于v^2';figurefori=1:3[t,y]=ode45('znxpydfun',[0:0.01:10],[0,3,0,5],[],b(i),p(i),m);H{i}=max(y(:,3))T{i}=t(find(y(:,3)==H{i}))vx0{i}=y(end,2)subplot(2,1,1)axis([06-702]);holdonxlabel('x')ylabel('y')plot(y(:,1),y(:,3));subplot(2,1,2)axis([01004])holdonxlabel('t')ylabel('dx/dt')text(px(i),py(i),strdd{i});plot(t,y(:,2))end函数文件是一个独立的文件,文件名为znxpydfun.mfunctionydot=znxpydfun(t,y,flag,b,p,m)ydot=[y(2);-b/m*y(2)*(y(2).^2+y(4).^2)^(p/2);y(4);-9.8-b/m*y(4)*(y(2).^2+y(4).^2)^(p/2)];