0.4 0.5 0.6 0.7 0.8 0.9 -0.4797 -0.4939 -0.5090 -0.5237 -0.5378 -0.5514 -0.4882 -0.5022 -0.5173 -0.5318 -0.5458 -0.5593 -0.5089 -0.5227 -0.5376 -0.5520 -0.5658 -0.5791 -0.5366 -0.5502 -0.5649 -0.5790 -0.5927 -0.6058 -0.5738 -0.5871 -0.6014 -0.6153 -0.6287 -0.6415 -0.6272 -0.6407 -0.6553 -0.6694 -0.6830 -0.6960

当Xg=.8896时 马赫数 Ma 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9

3.2.12 插值算法

攻角ɑ(?) 0 2 4 6 8 10 -0.6179 -0.6384 -0.6600 -0.6805 -0.6999 -0.7182 -0.6207 -0.6410 -0.6626 -0.6830 -0.7024 -0.7207 -0.6253 -0.6455 -0.6670 -0.6874 -0.7067 -0.7249 -0.6319 -0.6521 -0.6734 -0.6937 -0.7129 -0.7310 -0.6424 -0.6624 -0.6835 -0.7036 -0.7226 -0.7406 -0.6669 -0.6866 -0.7074 -0.7272 -0.7459 -0.7636 -0.6997 -0.7190 -0.7395 -0.7589 -0.7774 -0.7948 -0.7435 -0.7624 -0.7824 -0.8014 -0.8194 -0.8365 -0.8069 -0.8266 -0.8474 -0.8672 -0.8859 -0.9035 气动数据插值——等距双变元样条插值；

推力、重心、转动惯量等——不等距一元线性插值。

4.弹道计算的积分算法（即龙格库塔法）

龙格库塔四阶法的本质是利用微分方程右端子函数在若干点上的函数值的线性组

合来构造近似公式，构造时要求近似公式在任一点处的Taylor展开式与解处的Taylor展开式的前几项重合，从而使近似公式达到所需的阶数，和更高的精度。

四阶龙格库塔法的计算公式如下：

?dy?f?x,y??对于形如?dx??y(x0)?y0计算公式为：

x??a,b?y?Rn 的微分方程组

K1?f?xi,yi?

hh??

K2?f?xi?,yi?K1?22??hh??K3?f?xi?,yi?K2?22 ??

K4?f?xi?h,yi?hK3?yi?1?yi?h?K1?2K2?2K3?K4? 65.数据仿真

5.1程序编制

程序编制的基本思路如下： main函数 实参 Rungkuta子函数 插值数据 实参 主函数 结果 Shujuchuli子函数，主要对已有数据进行插值，并返回所需点的插值结果数据。

注：（1）程序中所用变量如X，Y和方程式中导弹位置对应； （2）程序中Wz为俯仰角速度；

（3）程序中Big_theta为导弹俯仰角；

（4）程序中Theta指弹道倾角；

（5）程序中alpha1指每一个迭代时刻的导弹攻角； （6）程序中M指导弹质量。

（7）程序中其它参数如cx,cy,L,S等完全与方程（1），（2）一一对应。 （8）剩余为程序中的中间变量，可在程序中轻松辨识。 5.2 仿真结果及图线

运行程序得到如下导弹飞行过程中的各变化参数的变化曲线：

导弹飞行弹道曲线150100Y/m500050010001500X/m200025003000

导弹飞行速度曲线300250200V/(m/s)1501005000246t/s81012

导弹俯仰角变化曲线2015105?/:o0-5-10-150246t/s81012