质点弹道.doc

基于matlab质点弹道数值计算1. 无控弹质点弹道模型(不考虑姿态运动影响)向量形式 其中：、分别为地面坐标系、轴上的单位向量。地面坐标系是与地球表面固连的坐标系。坐标原点定义在炮口断面中心，轴沿水平线指向射击方向，轴垂直向上，轴按右手法则确定。为速度向量，为风速向量，、分别为速度向量在地面坐标系三轴上投影向量的幅值，风速向量的分量定义类同。为空气密度为特征面积，对于轴对称弹取，为弹径。为弹箭质量为重力加速度向量，一般计算时取。 为阻力系数(作业中忽略诱导阻力) 分量形式表示 质点位置(射程方向) (飞行高度) (侧偏) 2. 仿真模型参数空气密度：按照美国标准大气推广委员会(COESA)提出的模型进行计算(随高度变化)，建议采用SIMULINK中的COESA Atmosphere Model模块。当地音速： 按照美国标准大气推广委员会(COESA)提出的模型进行计算(随高度变化)。弹径105mm质量14.97kg阻力系数随飞行马赫数的变化如下表：马赫数00.8750.9250.9650.991.0251.0851.191.351.822.5阻力系数0.1240.1240.150.20.350.3750.4150.415.3850.3350.3180.276 仿真时采用线性插值的方式得到对应飞行状态的阻力系数，建议采用SIMULINK中的Lookup Table。模块。仿真条件：炮口初速，发射点海拔高度为0，落点海拔为0。3. 作业内容分别在射角、的情况下绘制弹丸飞行速度、飞行高度随时间变化曲线，并求出全飞行时间、弹道顶点高度、落点侧偏、射程(发射点到落点的水平距离)、落速(落点飞行速度)、落角(落点弹道倾角)。采用ODE45算法求解有关微分方程。4. 作业提交方式交作业时以WORD文档形式打印提交，必须包含图形数据(以图形形式)和计算结果(以表格形式)。源程序以文档附件形式提交：应包含Simulink仿真程序截图及有关的.m脚本文件，保证作业结果可以复现。空军标准气象条件子函数:当地声速 C.mfunction C=C(y);T0=288.34;Ty=T0-0.00586.*y;Cy=20.046.*sqrt(Ty);C=Cy;End重力加速度 g.mfunction g=g(y);g0=9.806;R=6371000;gy=g0.*(1-2.*y/R);g=gy;end密度 rou.mfunction rou=rou(y);rou0=1.225;rouy=rou0.*(1-2.0323e-5.*y)4.83;rou=rouy;end解微分方程组主函数 maguo.m function d=maguo(t,x);m=14.97;d=0.105;S=pi/4.*d2;Ma=0 0.875 0.925 0.965 0.99 1.025 1.085 1.19 1.35 1.8 2 2.5;Cx0=0.124 0.124 0.15 0.2 0.35 0.375 0.415 0.415 0.385 0.335 0.318 0.276;p1=polyfit(Ma,Cx0,5);d= -(polyval(p1,x(7).*S.*rou(x(5).*x(8).*(x(1)-15)./m./2; -(polyval(p1,x(7).*S.*rou(x(5).*x(8).*(x(2)-0)./m./2-g(x(5); -(polyval(p1,x(7).*S.*rou(x(5).*x(8).*(x(3)-10)./m./2; x(1); x(2); x(3); x(7)-x(8)./C(x(5); x(8)-sqrt(x(1)-15).2+(x(2)-0).2+(x(3)-10).2); % x(9)-sqrt(x(1).2+(x(2).2+(x(3).2); %end调用：clearclcfor k=1:3sita0=30,45,60;t=37.848,54.465,70.81;v0=493;vx0=v0*cos(sita0(k)/180*pi);vy0=v0*sin(sita0(k)/180*pi);vz0=0;x0=0;y0=0;z0=0;Ma0=v0/C(0); %Cvm0=sqrt(vx0-15).2+(vy0-0).2+(vz0-10).2); %tspan=0,t(k);Xx0= vx0 vy0 vz0 x0 y0 z0 Ma0 vm0 v0; %M=eye(7);M(7,7)=0;M(8,8)=0;M(9,9)=0;options=odeset(Mass,M);t,x=ode15s(maguo,tspan,Xx0,options); %maguo.mfiguresubplot(2,3,1);plot(t,x(:,9);grid on subplot(2,3,2)plot(t,x(:,5);grid on%subplot(2,3,3); plot(t,x(:,4);grid on subplot(2,3,4);plot(t,x(:,5);grid on%subplot(2,3,5);plot(t,x(:,6);grid onsubplot(2,3,5);plot(t,x(:,8);grid onshegao=max(x(:,5)luosu=x(end,9)shec=sqrt(x(end,4).2+x(end,5).2+x(end,6).2)luojiao=acos(x(end,1)./(sqrt(x(end,1).2)+x(end,2).2)./pi.*180cepian=x(end,6)figureplot3(x(:,6),x(:,4),x(:,5);grid onclearend射角全飞行时间（S）弹道顶点高度落角射程(KM)落速落点侧偏