飞行方案大作业

)=sit(i)+(P*sind(a(i))+(Y(i)-m(i)*g*cos(sit(i))))/m(i)/V(i)*dt;Р V(i+1)=V(i)+(P*cosd(a(i))-(X(i)+m(i)*g*sin(sit(i))))/m(i)*dt;Р x(i+1)=x(i)+V(i)*cos(sit(i))*dt;Р H(i+1)=H(i)+V(i)*sin(sit(i))*dt;Р Q(i+1)=atan(-H(i)/(30000-x(i)))+pi;Р sit1(i+1)=k*(Q(i)-Q(1));Р H1(i+1)=H(i)+V(i)*sin(sit1(i));Р SIT(i+1)=(sit(i+1)-sit(i))/dt;Р r(i+1)=(H(i)^2+(30000-x(i))^2)^(1/2);Р R(i+1)=(r(i+1)-r(i))/dt;Р n(i+1)=acos(-R(i)/V(i))+pi;Р SIT1(i+1)=k/r(i)*(V(i)*sin(n(i)));Р T(i+1)=T0-0.0065*H(i+1);Р p(i+1)=p0*(T(i+1)/T0)^4.25588;Р q(i+1)=1/2*p(i+1)*V(i+1)^2;РendРplot(x,H); Рhold onР[t,y]=ode45('jieduan1',[0 39.0564],[250 0 0 7000]);Рplot(y(:,3),y(:,4));Рhold onР[t,y]=ode45('jieduan2',[39.0564 115],[192.768 -0.009 9100 2998.71]);Рplot(y(:,3),y(:,4));Р其中每一段的初始值,均为上阶段的结束值Р所以每一阶段计算结束后,需要再给出所有数据的结果,找到每一段距离相对应的数据,即为初始值。