基于匈牙利算法的轰炸机目标分配汇编.docxVIP

基于匈牙利算法的轰炸机目标分配姓名：钱 磊学号：1212471339专业：软 件 工 程引言2015年5月21日，我所服役的中国人民解放军空军某部出动中国最新轰六K型战略轰炸机编队穿越宫古海峡，前往西太平洋执行远海训练任务。这是中国空军首次穿越宫古海峡。飞越过程中，日本航空自卫队从那霸基地紧急起飞F-15战斗机对我机进行伴随监视飞行，当天晚上还公布了近距离拍摄的中国轰炸机照片。这还仅仅只是一次训练，真正的战场上危险程度只会更甚。虽然战略轰炸机航程远，载弹量也很大，但它毕竟是亚音速飞机，基本不具备与歼击机对垒的能力，因此如何以最小的风险完成多目标轰炸任务对于一个轰炸机编队来说尤为关键。模型所谓多目标攻击，就是一架飞机同时攻击多个目标或多架飞机协同攻击多个目标。多目标分配是有关多机空战中目标分配的一个最优化问题。对于一架轰炸机，其优势函数可通过Harold E.Bullock创造的函数利用相对距离、相对速度、相对离轴角度进行刻画。在实际中经常会遇到这样的问题，有n 项不同的任务，需要n 个人分别完成其中的一项，但由于任务的性质和各人的专长不同，因此各人去完成不同的任务的效率（或花费的时间或费用）也就不同。于是产生了一个问题，应指派哪个人去完成哪项任务，使完成n 项任务的总效率最高（或所需时间最少），这类问题称为分配问题或指派问题。多目标攻击的分配问题可描述为：设有n架轰炸机要轰炸n个目标，每架轰炸机只轰炸其中一个目标。用Cij表示第i架轰炸机轰炸第j个目标的优势函数值（通过Harold函数得到，优势数值越小，优势越高）。如何给每架轰炸机指派目标，是总优势数值和最小（总优势数值越小，总优势越高）为每一个轰炸机/目标对对应一个二值函数Xij,使得：1, 第架轰炸机对第个目标实施轰炸X（i,j）= 0，第架轰炸机对第个目标不实施轰炸则此目标分配问题的数学模型为：Z=min算法根据上述模型，此问题是特殊的0-1规划，也是特殊的运输问题，可以用这两种问题的求解方法求解，但根据分配问题的特殊结构，以及战场上对求解速度的要求，可使用匈牙利算法。以下结合实例来说明匈牙利算法的具体步骤：我部有四架轰六-K型战略轰炸机飞赴西太平洋对敌四个重要战略目标实施打击。根据某时刻敌我相对距离、相对速度、相对离轴角度经过Harold函数计算得到以下优势数值：目标飞机12341671122459833110445982第一步：缩减矩阵在价值矩阵中，每行元素减去该行最小元素完成行缩减，每列元素减去该列最小元素完成列缩减，显然最后得到的缩减矩阵每行每列均有0元素；-5第二步：初始分配在缩减矩阵中按“先少后多，先上后下，先左后右”的规则圈出0元素，每行每列至多圈一个，同时划去同行列其他0元素。当圈出数满n个时，则已得最优解，否则到第三步。第三步：行列标号无圈行标s；已标行i中，划0所在列j处标i;已标列j中，圈0所在行i处标j;交替2、3，直至无法标号。此时，若有一无圈0列标上号，则称该行突破，到第五步，否则到第四步。第四步：继续缩减用直线覆盖未标行和已标列；取未覆盖元素最小元素δ为调整量；已标行各元素减去δ，已标列各元素加上δ，新出现0元素划去，抹去直线，到上一步。第五步：调整分配1.若无圈0列j1标号为i1，则（i1，j1）处的划0改为圈0；2.若得圈0行i2标号为i2,则（i2,j2）处的圈0改为划0；交替进行1、2，直至标号s的行中得圈0，显然圈0数增加一个。此时当圈0数满n个时，则已得最优分配，否则转第三步。根据匈牙利算法，解得最优解矩阵为上图。对应最优解为：一号飞机轰炸四号目标，二号飞机轰炸一号目标，三号飞机轰炸二号目标，四号飞机轰炸三号目标，此时优势值和最小，即总优势最高。论证匈牙利算法是由匈牙利数学家Edmonds于1965年提出，因而得名。匈牙利算法是基于Hall定理中充分性证明的思想，它是部图匹配最常见的算法，该算法的核心就是寻找增广路径，它是一种用增广路径求二分图最大匹配的算法。首先，给出分配问题最优解性质：如果将分配问题价值矩阵每一行（列）的各元素减去该行（列）的最小元素，得到一个新的矩阵，则以该矩阵为价值矩阵的分配问题，其最优解与原问题最优解相同。利用该性质，可使原价值矩阵变换为含有多个0元素的新价值矩阵，而最优解不变。若能在新的价值矩阵中找到n个不同行且不同列的0元素，则可令其对应的Xij等于1，其他Xij为0，显然该解一定是最优解。匈牙利算法的基本思想正式指出一条选满n个异行异列0元素的途径。结论匈牙利算法在小规模分配问题中求解简便快捷，环境要求低，十分适应于现代战场对多目标分配问题快速求解的需要。结合Harold函数对空军轰炸机多目标打击分配问题具有较高的实用意义。结束语通过这学期对《运筹学概论》这门课的学习，我对运筹学这门学