导弹总体设计导弹总体设计防空导弹主级燃料相对质量因数计算防空导弹主级燃料相对质量因数计算实验报告实验报告02010603 班班杜泽弘杜泽弘060500指导老师：龚春林指导老师：龚春林实验结果：调用函数代码调用函数代码：dv/du 的表达程序： functionV=f1(u,Vo,y,alpha,theta) Is=2156;P=2.2;p0=5880;g=9.801;ys=340; Cxx =0.0043 0.0511 0.0651 0.0847 0.11200.0360 0.0436 0.0558 0.0736 0.09730.0308 0.0372 0.0481 0.0641 0.08490.0265 0.0323 0.0419 0.0560 0.07460.0222 0.0272 0.0356 0.0478 0.0644; Max=1.5 2.1 2.7 3.3 4.0; alphax=2 4 6 8 10*3.14/180; Cx=interp2(alphax,Max,Cxx,alpha,Vo/ys,spline); V=Is/(1-u)-rho(y)*Vo2*Cx*Is/2/P/p0/(1-u)-Is*sin(theta)/P;d/du 的表达程序： functiontheta=f2(u,V,y,alpha,thetao) Is=2156;P=2.2;p0=5880;g=9.801;ys=340; VT=420;yT=22000;xT=45070;q0=atan(yT/xT); theta=VT/yT*(Is/P/g+1/V2/sin(thetao)*(V*f3(V,thetao)-y*f1(u,V,y,alpha,thetao)/(1+(cot(q0)- u*Is*VT/P/g/yT)*cot(thetao);dy/du 的表达程序： functiony=f3(V,theta) Is=2156;P=2.2;g=9.801; y=Is/P/g*sin(theta)*V;dx/du 的表达程序： functionx=f4(V,theta) Is=2156;P=2.2;g=9.801; x=Is/P/g*cos(theta)*V;攻角的计算程序： functionalpha=f5(u,V,y,alphao,theta) Cyy=57.3*0.0302 0.0304 0.0306 0.0309 0.0311 0.03130.0279 0.0280 0.0284 0.0286 0.0288 0.02900.0261 0.0264 0.0267 0.0269 0.0272 0.02740.0247 0.0248 0.0251 0.0254 0.0257 0.02590.0226 0.0227 0.0231 0.0233 0.0236 0.02380.0209 0.0210 0.0213 0.0216 0.0219 0.0221; ys=340;Is=2156;P=2.2;p0=5880;g=9.801; May=1.5 2.0 2.5 3.0 3.5 4.0; alphay=1 2 4 6 8 10*3.14/180; Ma0=V/ys; Cyalpha=interp2(alphay,May,Cyy,alphao,Ma0,spline); alpha=(V*f2(u,V,y,alphao,theta)+Is*cos(theta)/P)/(rho(y)*V2*Cyalpha*Is/2/P/p0/(1-u)+Is/(1-u);密度的计算程序： functionrhorho=rho(y) T0=288.15;rho0=0.12495; if y0.009u0=u1;V0=V1;theta0=theta1;y0=y1;x0=x1;alpha0=alpha1; elsebreak end endsprintf(燃料相对质量因数为:u=%4.4f,u1)中间计算的迭代所得的参数： 1.攻角的迭代结果： 10.9057 10.1890 9.9535 9.7525 9.5843 9.4410 9.3191 9.2169 9.1332 9.0676 9.0198 8.9897 8.9774 8.9832 9.0076 9.0511 9.1146 9.1990 9.3057 9.4361 9.5920 9.7757 9.9896 10.2370 10.5212 10.8465 11.2176 11.7528（单位为） 2.阻力系数的迭代结果 Cx： -0.0234 0.1246 0.1131 0.1086 0.1046 0.1011 0.0980 0.0953 0.0929 0.0908 0.0889 0.0873 0.0859 0.0848 0.0838 0.0831 0.0826 0.0824 0.0823 0.0825 0.0830 0.0837 0.0848 0.0861 0.0879 0.0901 0.0928 0.09603.升力线斜率的迭代结果：yC1.7470 1.7967 1.7781 1.7644 1.7513 1.7388 1.7267 1.7150 1.7037 1.6927 1.6821 1.6719 1.6619 1.6523 1.6429 1.6338 1.6250 1.6164 1.6081 1.6000 1.5920 1.5842 1.5765 1.5689 1.5612 1.5534 1.5453 1.5369 1.5281 1.5186 1.5083 1.49624导弹主级燃料相对质量因数u=0.3467。