基于LabVIEW的弹道解算软件设计

摘要：针对铅垂平面内的外弹道模型，对基于LabVIEW的弹道解算方法进行研究，提出了以四阶龙格库塔算法为核心的解算方法。使用LabVIEW编写了弹道解算程序，在初始点状态已知的情况下，快速高效的递推出外弹道轨迹与落点坐标。程序能够对数据进行处理与显示。人机界面友好，易于操作。对于外弹道解算具有一定的应用价值。

现代导弹从出现已经过去了50多年，这50年当中导弹技术不断发展。其中包括了伴随着导弹而发展起来的新型学科-外弹道学，它不仅研究导弹飞行过程中质心的运动过程还研究了姿态控制学与制导理论等。外弹道学涉及到的技术范围广且需要各个学科相互联系。本文研究的是铅垂平面内，弹丸质心的运动轨迹，建立质心运动方程组，已知弹丸初始状态利用四阶龙格库塔方法求解发射轨迹。据此编写了LABVIEW程序实现解算过程并且将轨迹显示出来。

1 质心运动方程

研究弹丸质心的运动过程，第一步要做基本假设;

1)整个发射至击中目标的运动过程中，攻角为0。

2)弹丸是轴对称的。

3)地面为平面。

4)飞行过程中重力加速度大小不变且垂直向下。

5)科氏加速度为0。

6)大气条件为标准大气条件，无风。

在此基础之上，作用于弹丸的力只有重力与空气阻力。依此弹丸的质心运动方程为式(1)。

由于ax=cHτ(y)G(vτ)v，所以得到了弹丸质心的运动方程组为式(2)。

式中：ax为空气阻力加速度;c为具有弹丸特征的弹道系数;H(y)为空气密度函数;G(v)为阻力函数。

2 四阶龙格库塔

四阶龙格库塔法广泛用于模拟仿真应用中的弹道解算，在工程中被称为“高精度单步算法”。由欧拉公式导出，k1为一阶精度的欧拉公式，k2用xi点处的k1与xi+1的k2的

平均值作为平均斜率的近似值。就会得到二阶的欧拉公式。以此类推得到四阶龙格库塔方程如式(3)所示。龙格库塔法因单步四次计算，计算量大，因此广泛使用于计算机仿真计算。

3 程序设计

根据弹丸质心运动方程组，使用四阶龙格库塔法，在LABVIEW软件平台编写软件解算弹道轨迹。在LABVIEW中依据Runge-Kutta算法VI建立解算方程。Runge-KuttaVI如图2所示。

此VI需要初始状态的数据。根据弹丸质心运动方程组，初始点已知。当t=0时，x=0;y=0;u=V0COSθ0;w=vosinθ0。θ0为发射角已知。X(name of variables)为变量名称将质心运动方程组所涉及的变量名设置为字符串数组。time start为初始时间为0，time end为结束时间。h(steprate)为步长设置为1。X0为对应于X变量名的数据信息，这里设置u=700;w=700;y=0;x=0。time设为t。F(X，t)为对应变量名的函数名。程序设置如图3所示。

4 结论

以固定的弹体，在标准大气压条件下，c为常数;x，y为0;纵向与横向初速度分别设为100、200、300、400、500。记录解算数据：时间、最大高度、最远距离如表1所示。

初速度为707 m/s解算出的轨迹曲线如图4所示。

经过数据仿真实验，得出结论。此弹道解算软件具有可靠性，解算出的弹道数据对弹道研究具有一定的现实意义。