飞行方案大作业

88;Рq(1)=1/2*p(1)*V(1)^2; %大气密度计算公式РCx(1)=0.2+0.005*a(1)^2;РCy(1)=0.25*a(1)+0.05*det(1)*180/pi; %升力系数РY(1)=Cy(1)*q(1)*S;РX(1)=Cx(1)*q(1)*S;РSIT(1)=(P*sind(a(1))+(Y(1)-m(1)*g*cos(sit(1))))/m(1)/V(1);РQ(1)=atan(-H(1)/(30000-x(1)))+pi;Рr(1)=6708.2039;РR(1)=-V(1)*cos(Q(1));Рn(1)=Q(1)+pi;РSIT1(1)=k/r(1)*(V(1)*sin(n(1)));Рmza=-0.1; %俯仰力矩系数对攻角的偏导数Рmzdet=0.024; %俯仰力矩系数对舵偏角的偏导数Рt=0;Рi=0;Рdt=0.01;Рms=0.46; %质量秒消耗量Рwhile H>0 & H1>0 %运用迭代法求解Р i=i+1;Р t=t+dt;Р det(i+1)=k1*(sit(i)-sit1(i))+k2*(SIT(i)-SIT1(i));Р a(i+1)=-mzdet/mza*det(i)*180/pi;Р Cy(i+1)=0.25*a(i)+0.05*det(i)*180/pi;Р Cx(i+1)=0.2+0.005*a(i)^2;Р Y(i+1)=Cy(i)*q(i)*S;Р X(i+1)=Cx(i)*q(i)*S;Р m(i+1)=m(i)-ms*dt;Р sit(i+1)=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;