fun.m

利用四阶runge-kutta法

%M文件fun.m
function f=fun(t,y)
v=sqrt(y(2)^2+y(4)^2)
Cy=[(v/340)^2,v/340,1]*[-0.0173;0.0434;0.3275]*(pi/90);
Cx=[(v/340)^2,v /340,1]*[0.0002,0.0038,0.1582;-0.0022,-0.0132,-0.8520;-0.0173,0.0434,1.9712]*[(pi/90)^2;pi/90;1];
p=1.225*exp(-0.00015*y(3));
s=1.2;
f=[y(2),-(p*y(2)^2*s*Cx)/460,y(4),(p*y(2)^2*s*Cy)/460-9.8]'