资源描述：
《导弹总体设计实验报告》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、导弹总体设计《防空导弹主级燃料相对质量因数计算》——实验报告02010603班杜泽弘060500指导老师：龚春林实验结果：调用函数代码：dv/du的表达程序：function[V]=f1(u,Vo,y,alpha,theta)Is=2156;P=2.2;p0=5880;g=9.801;ys=340;Cxx=[0.00430.05110.06510.08470.11200.03600.04360.05580.07360.09730.03080.03720.04810.06410.08490.02650.03230.04190.05600.
2、07460.02220.02720.03560.04780.0644];Max=[1.52.12.73.34.0];alphax=[246810]*3.14/180;Cx=interp2(alphax,Max,Cxx,alpha,Vo/ys,'spline');V=Is/(1-u)-rho(y)*Vo^2*Cx*Is/2/P/p0/(1-u)-Is*sin(theta)/P;dθ/du的表达程序：function[theta]=f2(u,V,y,alpha,thetao)Is=2156;P=2.2;p0=5880;g=9.801;ys=
3、340;VT=420;yT=22000;xT=45070;q0=atan(yT/xT);theta=VT/yT*(Is/P/g+1/V^2/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的表达程序：function[y]=f3(V,theta)Is=2156;P=2.2;g=9.801;y=Is/P/g*sin(theta)*V;dx/du的表达程序：function[x]=f4(V
4、,theta)Is=2156;P=2.2;g=9.801;x=Is/P/g*cos(theta)*V;攻角的计算程序：function[alpha]=f5(u,V,y,alphao,theta)Cyy=57.3*[0.03020.03040.03060.03090.03110.03130.02790.02800.02840.02860.02880.02900.02610.02640.02670.02690.02720.02740.02470.02480.02510.02540.02570.02590.02260.02270.02310.
5、02330.02360.02380.02090.02100.02130.02160.02190.0221];ys=340;Is=2156;P=2.2;p0=5880;g=9.801;May=[1.52.02.53.03.54.0];alphay=[1246810]*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)*V^2*Cyalpha
6、*Is/2/P/p0/(1-u)+Is/(1-u));密度的计算程序：function[rhorho]=rho(y)T0=288.15;rho0=0.12495;ify<11000T=288.15-0.0065*y;rhorho=rho0*(T/T0)^4.25588;elseify<20000T=216.65;rhorho=0.0371/exp((y-11000)/6341.62);elseT=216.65+0.01*(y-20000);rhorho=0.088035*(216.65/T)^35.1632;endend相对质量因数迭代
7、计算程序：function[u]=f(x,y)Is=2156;P=2.2;p0=5880;g=9.801;VT=420;yT=22000;xT=45070;q0=atan(yT/xT);u=(P*g*yT/Is/VT*(cot(q0)-x/y));主程序代码：V0=500;theta0=28*3.14/180;alpha0=1.5*3.14/180;x0=674;y0=329;Is=2156;g=9.801;P=2.2;p0=5880;VT=420;q0=atan(22000/45070);yT=22000;h=0.009;%h为迭代步
8、长u0=0;%假设燃料相对质量因数的初始值while(1)V1=V0+h*f1(u0,V0,y0,alpha0,theta0);theta1=theta0+h*f2(u0,V0,y0,alpha0,thet
