求汽车空气阻力的MATLAB程序,有空气阻力抛射体的运动(matlab)

有空气阻力抛射体的运动

微固学院 2012030010023 曾鹏灏

由于空气阻力大小和方向始终与抛体的速度有关，会不停地发生改变，对求解抛射体的运动轨迹增加极大的难度。因此，我们须借助matlab等数学工具软件，从数量上的累积来解决此问题。

Matlab程序如下：

g=9.8; m=0.145; C=0.5; r=0.0366; A=pi*r^2; rou=1.2; D=(rou*C*A)/2; detat=0.01;

theta=(35/180*pi); x(1)=0； y(1)=0; V(1)=50;

Vx(1)=V(1)*cos(theta); Vy(1)=V(1)*sin(theta); n=0 ; y(1)=0.0001; while(y(n+1)>0) n=n+1;

V(n)=sqrt(Vx(n)^2+Vy(n)^2); ax(n)=-(D/m)*V(n)*Vx(n); ay(n)=-g-(D/m)*V(n)*Vy(n);

Vx(n+1)=Vx(n)+ax(n)*detat; Vy(n+1)=Vy(n)+ay(n)*detat; if(Vy(n+1)<=0.1&&Vy(n+1)>=0) M=n+1 end

x(n+1)=x(n)+Vx(n)*detat+0.5*ax(n)*detat^2; y(n+1)=y(n)+Vy(n)*detat+0.5*ay(n)*detat^2; if(y(n+1)<=0.1&&y(n+1)>