以射弹中心(瞄准点)为原点建立二维坐标系，射弹正面为X轴，纵深为Y轴，而射击误差在二维平面上的概率密度函数为

p ( x , y ) = ρ 2 π E x E y e ρ 2 ( x 2 E x 2 + y 2 E y 2 ) (1)

其中， ρ 是数值为0.477的常数。 E x 和 E y 分别是X轴和Y轴上的射击误差。通常使用射弹散布偏差来表示射击误差，用子母E表示，单位是m。

同理，跑道的观测误差也可用散布概率偏差来表示，以跑道中心为原点建立二维坐标系，也可得到与射击误差类似的结果。

首先需要计算弹药的等效封控面积，结合跑道的毁伤封控原理，等效封控面积可等效为每发弹药对跑道内部最小起降窗口的封控或影响面积，所以等效封控面积要比实际毁伤面积大得多，并且每发弹对机场跑道的封控区域基本可以认为是均匀的，假设二维平面中某点在跑道内部的概率为 P m ( x , y ) ，用来表征跑道的观测误差，同时弹药的射击散布误差用该点被弹药的封控区域覆盖的概率 P n ( x , y ) 表示，所以弹药的等效封控面积为：

S n = ∬ P m ( x , y ) × P n ( x , y ) d x d y ( x , y ) ∈ R 2 (2)

后面将对 P m ( x , y ) 和 P n ( x , y ) 分别进行求解。

目标为一正面长度为B、纵深长度为A的矩形跑道，在其二维平面上随机取一点 P ( x , y ) ，以该点为中心，作平行于X轴和Y轴长度为A和B的矩形，把该区域记为Ω，即该点落在跑道范围内的概率为 [ 7 ] ：

P m ( x , y ) = ∯ p ( x , y ) d x d y (3)

已知式(1)中 2 ρ E x x 和 2 ρ E y y 服从相互独立的标准正态分布，即

( 2 ρ E x x ， 2 ρ E y y ) ~ N ( 0 , 0 , 1 , 1 , 0 ) (4)

代入式(3)得：

P m ( x , y ) = { Φ [ 2 ρ E x ( x + B 2 ) ] − Φ [ 2 ρ E x ( x − B 2 ) ] } × { Φ [ 2 ρ E y ( y + A 2 ) ] − Φ [ 2 ρ E y ( y − A 2 ) ] } (5)

假设子弹药封控区域为一个 a × b 的类矩形区域，该区域通常由目标机场的最小起降窗口尺寸和弹坑大小决定的。设弹着点坐标为 Q ( x i , y i ) ，任取点 P ( x , y ) 相对于弹着点的坐标为 ( x − x i , y − y i ) ，根据式(5)可以计算得点 P ( x , y ) 被单发弹药封控的概率为：

P i ( x , y ) = { Φ [ 2 ρ E x ( x − x i + b 2 ) ] − Φ [ 2 ρ E x ( x − x i − b 2 ) ] } × { Φ [ 2 ρ E y ( x − y i + a 2 ) ] − Φ [ 2 ρ E y ( x − y i − a 2 ) ] } (6)

可知点 P ( x , y ) 被所有弹药封控的概率为：

P n ( x , y ) = 1 − ∏ i = 1 n [ 1 − P i ( x , y ) ] (7)

需要计算的近似最优值可以抽象为：如何确定多发弹药的位置使得弹药对机场跑道的封控面积S更大，或封控的最小起降窗口更多，即

{ max x 1 , y 1 , ⋯ , x n , y n max x 1 , y 1 , ⋯ , x n , y n S x 1 , y 1 , ⋯ , x n , y n ϵ R (8)

首先需要对跑道内部所有坐标值(变量)进行编解码，其目的是建立变量(表现型)到基因序列(基因型)的一一映射关系。在实际打击过程中，弹着点必须在跑道范围内才能起到封控的作用，即变量取值要在跑道矩形区域内部才有意义。另外跑道中心的观测误差在X轴和Y轴上均服从相互独立的标准正态分布，

所以跑道在两个方向上误差超过3σ的概率不超过0.3%，故变量的取值范围分别为 [ − 3 E x , A + 3 E x ] 、 [ − 3 E y , B + 3 E y ] 。

本文中遗传算法采用二进制对坐标值进行编解码，以 轴的变量取值为例，设 ε 为数值精度，则有 Δ x = [ A + 6 E x ε ] ，结合机场跑道尺寸和最小起降窗口封控特性，为确保封控结果的可靠性，数值精度和

变量取值要贴切实际，精度不能太低。确定变量取值后，将范围内的整数与固定长度的二进制数进行一一映射，二进制码的长度为变量最大值对应的二进制数的长度，Y轴变量的编码原理同上。

理论上种群规模由随机函数生成，但考虑到跑道尺寸以及精度要求，本文设定种群规模为200，基因串的适应度函数为：

f = S ( g x , g y ) (9)

其中 g x 和 g y 为X轴和Y轴上取值的二进制编码。

通过归一化方法将适应度转化为基因遗传给下一代的概率，即：

P i = f i ∑ f i (10)

根据遗传概率，迭代若干代之后，遗传概率低、封控效果差的基因串会被逐渐淘汰，概率高且封控效果好的会在种群中扩散，最终趋于稳定 [ 4 ] 。

在遗传算法中利用基因突变来克服过早收敛问题，若基因突变概率为p，每个基因串的长度为n，则每次基因突变的基因个数为 200 × p × n ，将所有基因串排成一排，在其上随机选取位置，该位置上的基因发生突变 [ 8 ] 。

遗传算法基本步骤如图1。

本文采用反跑道火箭弹，为精确制导弹药，可认为每发弹药的射击散布误差相同且相对较小 [ 9 ] 。在仿真实验中，每颗弹药的等效封控面积与目标机场的最小起降窗口尺寸以及弹坑大小有关，该区域可近似为(600 m × 30 m)的矩形，机场跑道设定为(2400 m × 50 m)的矩形区域，另外反跑道火箭弹的射击误差和跑道的观测误差均设定为10米。

图1. 遗传算法流程图

在考虑散布偏差的情况下，单发弹药的等效封控面积示意图如图2所示，用区域和颜色代表等效封控面积以及各点被“封控”的概率。不难看出每个点被该弹药封控的概率与点到散布中心的距离有关，越靠近散布中心“封控”概率越大。仿真所计算的是总等效封控面积达到最大时的多个反跑道火箭弹弹着点坐标。

图2. 单发反跑道火箭弹的等效封控面积示意图

图2中心黑色区域表示弹着点(弹坑)，其余颜色区域表示“封控”区域，颜色越深，“封控”概率越大，白色区域表示跑道的一部分。

该过程是在指定弹药消耗量的条件下计算弹着点的近似最优位置，且本文是针对上述第一阶段的研究，所以弹药消耗量不能大于完成封控的最低弹数。但弹药消耗量较少时结果存在多种情况，难以体现结果的准确性，所以该仿真实验以7发弹药消耗量为例进行，并展示仿真结果，弹着点分布图如表1。

表1. 弹着点近似最优解坐标

自图3可以看出，在误差的影响下，各个弹着点的等效封控区域之间几乎没有重叠区域，弹着点位置与理想情况基本一致；图4显示随着遗传算法迭代次数的增加，未被封控区域面积迅速下降，反应速度较快，并最终趋于稳定。

图3. 弹着点近似最优解及封控概率示意图

图4. 迭代次数与封控面积关系示意图

在确保准确性的条件下遗传算法相比于传统概率方法和蒙特卡洛算法效率上大大提升，缩短计算时间，而迭代次数能够较好的反应效率问题，所以本文还将在指定弹药消耗量情况下，分析跑道尺寸与迭代次数之间的关系，并给出不同跑道尺寸对应所需的迭代次数，并与相应数值精度下的蒙特卡洛算法的迭代搜索次数作对比。

跑道尺寸通常由目标机场起降的机型和飞机流量决定，不同机场跑道之间宽度相差不大，对迭代次数的影响相对较弱。该仿真实验仍然以7发弹药消耗量为例进行，跑道宽度为50 m不变，长度为1800 m~4000 m的变量，值得注意的是，遗传算法求解结果为近似最优解，所以本文所记录的迭代次数为未被封控区域下降到可接受误差范围时的数据，可接受范围为误差不超过5%，仿真结果如图5。

图5. 跑道长度与迭代次数关系示意图

自图5不难看出，随着跑道长度的增大迭代次数均逐渐增大，并且增长趋势慢慢变缓；作为求解全局近似最优解的良好算法，相对于蒙特卡洛算法，遗传算法在搜索效率上有显著提高，对于正常尺寸的跑道，大约只需要迭代2000~5000次未被封控等效面积便可满足可接受的误差范围。