movement.m

有关matlab r2007的实用教程

function rdot=movement(t,r)
c=0.01;
g=9.81;
m=1;                                          	%空气阻力系数及重力加速度（m/s^2）
vm=sqrt(r(3)^2+r(4)^2);                        	%速度大小
rdot=[r(3);r(4);-c*vm*r(3)/m;-c*vm*r(4)/m-g;]       	%运动方程