znxpyd.m

来自「研究阻尼斜抛运动的MATLAB程序」· M 代码 · 共 32 行

m=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';

figure
for i=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,3,i)
axis([0 6 -70 2]);
hold on
xlabel('x')
ylabel('y')
comet(y(:,1),y(:,3));
plot(y(:,1),y(:,3));
subplot(2,3,i+3)
axis([0 10 0 4])
hold on
xlabel('t')
ylabel('dx/dt')
text(px(i),py(i),strdd{i});

comet(t,y(:,2))
plot(t,y(:,2))
end