pingpao.m

matlab仿真的平抛运动

%  平抛运动
%  copyright : yinbo
%  datetime   2006-5-12
h0=100;
g=9.8;
tf1=0:0.1:sqrt(2*h0/g);
x1=2*tf1;
h1=100-0.5*g*tf1.^2;
v1=g*tf1;
set(gcf,'DoubleBuffer','on');
axis([0 50 0 100])
grid on
text(3,100,0.8,'平抛运动(非弹性碰撞)','fontsize',20,'color','r');
handles=line('color','b','linestyle','.',...
    'markersize',30,'erasemode','xor');
   Tz=title('t=0');    
   tt=0;
   dtt=0.01;
                                                 %使小球运动
for i=2:length(h1)
    set(handles,'xdata',x1(i),'ydata',h1(i));    %小球位置
    drawnow;   
      tt=tt+dtt;
    set(Tz,'string',['t=',num2str(tt)]);
    pause(0.05)                                   %刷新屏幕
end
hold on
v2=0.6*v1(length(h1));
tf=0:0.1:2*v2/g;
x2=2*tf+x1(length(h1)); 


while v2>0.01
  h2=v2*tf-0.5*g*tf.^2;
set(gcf,'DoubleBuffer','on');

for j=2:length(h2)
     set(handles,'xdata',x2(j),'ydata',h2(j));
     drawnow
      tt=tt+dtt;
     set(Tz,'string',['t=',num2str(tt)]);
     pause(0.05)
 end
 hold on
 x2=x2(length(tf));
 v2=0.6*v2;
 tf=0:0.1:2*v2/g;
x2=x2+0.9*2*2*tf;
end