matlab实现简单物理实验模拟资料

ball.mclear vx=input('please input the speed ration of Vx and Vi:')

一. ball.m clear vx=input('pleaseinputthespeedrationofVxandVi:'); k=input('coefficientofrebound(<1):'); ifk>=1return,end tm=(1+k)/(1-k); xm=vx*tm; figure plot([0,xm],[0,0],'LineWidth',3) gridon axisequal axis([0,xm,0,1]) title('Trajectoeyofball','FontSize',16) xlabel('horizontaldistance(m)','FontSize',16) ylabel('Verticalheight(m)','FontSize',16) dt=0.001; tm=1; t=0:dt:tm; x=2*vx*t; y=1-t.^2; i=0; holdon whiletm>0.01 comet(x,y) plot(x,y,'LineWidth',2) i=i+1; tm=2*k.^i; t=0:0.001:tm; x=x(end)+2*vx*t; y=2*k^i*t-t.^2; end 二. collision.m clear m=0:0.5:2; n=length(m); theta=0:5:90; th=theta*pi/180; [M,TH]=meshgrid(m,th); PHI=2*asin(sin(TH/2)./(1+M))*180/pi; figure(1) plot(theta,PHI(:,1),'o-',theta,PHI(:,2),'s-',theta,PHI(:,3),'d-',... theta,PHI(:,4),'p-',theta,PHI(:,5),'h-') xlabel('\it\theta\rm/\circ','FontSize',16) ylabel('\it\phi\rm/\circ','FontSize',16) title('Angleofballandobject','FontSize',16) gridon legend([repmat('Ratioofmass\itM/m\rm=',n,1),num2str(m')],2) DE=M./(1+M).*(1-cos(TH)); figure(2) plot(theta,DE(:,1),'o-',theta,DE(:,2),'d-',theta,DE(:,3),'s-',... theta,DE(:,4),'p-',theta,DE(:,5),'h-') xlabel('\it\theta\rm/\circ','FontSize',16) ylabel('\Delta\itE/mgl','FontSize',16) title('MechanicalEnergy','FontSize',16) gridon legend([repmat('Rationofmass\itM/m\rm=',n,1),num2str(m')],2) 三. niudunhuan.m clear rm=5; r=-rm:0.01:rm; [X,Y]=meshgrid(r); R=sqrt(X.^2+Y.^2);