飞行力学纵向配平的计算需要解非线性方程组,解非线性方程组要将方程写成F(X1,X2,...)=。。。。的形式,需要注意的是有可能求得的值为局部最优,这与初始值相关,
function F=trim(X)
s=75.12 ;rho=0.4135 ;v=70 ;m=11984 ;g=9.8 ;
T=X(1);alpha=X(2);delta_e=X(3);
F(1)=T*cos(alpha)-0.5*rho*v^2*s*(1.5684*alpha^2-1.17e-6*alpha+0.0102+2*(-3.0489*alpha^2-0.1173*alpha-0.0039)*delta_e^3+2*(-5.0376*alpha^2+0.1937*alpha-0.0065)*delta_e^2+2*(-3.6877*alpha^2+0.1418*alpha-0.0048)*delta_e);
F(2)=T*sin(alpha)+0.5*rho*v^2*s*(4.4955*alpha-0.0008-2*0.0940*delta_e^3-2*0.0005*delta_e^2+2*0.1827*delta_e)-m*g;
F(3)=-0.0969*alpha+0.0003*2*delta_e^2+0.0112*2*delta_e^3-0.0218*2*delta_e;
%注意函数名和文件名要一致,在命令窗口输入,是否成功会有提示
<<X0=[10000,0.1,-0.3]
<<X=fsolve('trim',X0),F=trim(X)