matlab模拟平抛实验,谁有平抛物体运动模拟的MATLAB程序，简单的模拟……

匿名用户

1级

2018-11-21 回答

m=10;%质量

h=10;%高度

v0=0.5;%初速度

g=9.8;%重力加速度

u=0.15;%空气平均阻力系数

tt=0.1;%物体与地面的接触时间

a1=g-u*g;

a2=g+u*g;

N=0;s=0;

h1=figure('units','normalized');

axis([0 14 -0.5 10]);

axis off;

ht=uicontrol('parent',h1,'units','normalized','style','text','position',[.7 .9 .3 .1],'string','0','fontsize',16);

b=line(0,10,'color','r','marker','.','markersize',25,'erasemode','xor');

ball=line(0,10,'color','r','linestyle','-','linewidth',2,'erasemode','none');

line([0,14],[-0.02,-0.02],'color','b','linestyle','-','linewidth',2);

line([0 0],[0 10],'color','b','linestyle','-','linewidth',2);

while N~=m*g

N=0;

v=sqrt(2*a1*h);

t1=sqrt(2*h/a1);

for t=0:0.001:t1

sx=s+v0*t;

y=h-a1*t*t/2;

set(b,'xdata',sx,'ydata',y);

set(ball,'xdata',sx,'ydata',y);

drawnow

end

N=2*m*v*tt+m*g;

set(ht,'string',num2str(N));

t2=v/a2;

s=sx;

h0=v*v/(2*a2);

for t=0:0.001:t2

sx=s+v0*t;

y=v0*t-a2*t*t/2;

set(b,'xdata',sx,'ydata',y);

set(ball,'xdata',sx,'ydata',y);

drawnow

end

h=y;s=sx;

end