1、 指挥控制与仿真 ()引用格式:郭子杰,王军,张贤椿,等基于遗传算法的自导水雷雷位优化指挥控制与仿真,():,():基于遗传算法的自导水雷雷位优化郭子杰,王 军,张贤椿,庞景壮(南京理工大学自动化学院,江苏 南京)摘 要:自导水雷作为一种新式智能水雷,在现代海上战争中有举足轻重的地位。为提高由自导水雷组成的雷障对通过目标的拦截概率,首先对自导水雷的目标探测系统进行建模,获得自导水雷在战场环境下的控制水域范围,然后基于遗传算法,对自导水雷的布雷位置进行优化。通过对优化前后多个目标通过雷障时的情况进行仿真,获得目标在雷位优化前后的平均触雷概率,验证了采用遗传算法优化后自导水雷布设雷位的优越性,对自
2、导水雷的雷位布设具有指导意义。关键词:雷位优化;自导水雷;遗传算法;触雷概率中图分类号:文献标志码:,(,):,:;收稿日期:修回日期:作者简介:郭子杰(),男,硕士研究生,研究方向为作战任务筹划。王 军(),男,副研究员。自导水雷又称水鱼雷,由自导鱼雷和锚雷组合而成,由于其具有智能识别目标、控制水域面积大、机动能力强等优点,在海战场上有着举足轻重的地位。因为自导水雷单雷控制的水域面积大,所以封锁同样面积的水域所需自导水雷数量大大降低。而水雷对目标的毁伤能力很大程度上取决于目标的触雷概率。因此,如何在预定水域合理布设自导水雷,会直接影响自导水雷的作战效能。自导水雷的雷位设置主要受到两个因素的影
3、响。一是自导水雷目标探测系统的探测距离,探测距离取决于探测系统的性能参数、水域的水文条件、目标参数,在制定布雷计划之前需要对目标探测系统的实际探测距离进行预测。二是预定拦截的目标数量,在拦截单目标时,通常只需要布设 道雷线即可;但在拦截多个目标组成的舰艇编队时,由于要考虑先通过雷障的目标对雷障整体性的破坏,需要在一片水域内布设多枚自导水雷才能做到对多个目标的有效拦截。由于自导水雷探测面积大、探测到目标后不能立刻对目标进行毁伤,一个目标可能会触发多个水雷,对雷障完整性造成较大的破坏。因此本文主要研究如何优化自导水雷的雷位,减小通过雷障的目标对雷障完整性的破坏,从而提高雷障对多个目标的拦截能力。本
4、文首先基于声呐方程,对自导水雷目标探测系统的探测距离进行建模。然后依据建模结果,采用遗传算法对拦截多个目标时自导水雷的雷位进行优化,将随机布设自导水雷与使用遗传算法对水雷位置进行优化后,多目标通过雷障的触雷概率进行对比,验证优化的合理性。自导水雷探测范围建模与仿真自导水雷在水下潜伏时首先打开被动声呐,当被动声呐探测到目标后,会开启主动声呐对目标进行目标识别和定位,若判断此目标为预先设定的攻击目标后,自导水雷发射战斗部对目标发起攻击。由于主动、被动两种声呐的工作方式略有区别,需要对这两种声呐分别进行仿真。第 期指挥控制与仿真 声呐工作原理被动声呐仅通过接收机接收水域中传播过来的能量信号,当接收到
5、的能量信号强度超过设定阈值即判定此时探测到目标。被动声呐方程如下:()()而主动声呐会周期地向海水中发射某种形式的声波,如果声波信号在传播中遇到目标,一部分能量被反射回来,接收机接收并处理该能量信号,识别目标并获得目标位置。主动声呐主要受到背景干扰,有环境噪声和海洋混响两种,主动声呐方程有如下两种形式:()()式中,为发射源极;为传播损失;为目标强度;为背景噪声级;为检测阈值;为海洋混响;为接收指向性指数,接收指向性指数表示的是声呐接收器抑制非目标方向的干扰能力的参数,即指向性指数越高,抑制噪声的能力越强。传播损失水声信号在水中传播的损耗随着传播距离的增加而增加。若声呐的作用距离为,海水的声吸
6、收系数为,那么所能允许的单程最大传播损失为 ()式中,为衰减系数,计算过程中一般选择 。水域噪声和混响海洋环境噪声主要可以分为风浪噪声和雨水噪声,风浪噪声的产生机理近似于由于破碎波以及海面活动所产生的气泡自然脉动;降雨落在海面上的噪声来源于贴在海面上所形成的微小气泡。式()、()分别对风浪噪声和雨水噪声建模:()()()()式中,为频率,单位为;为海面与海水的温差;为风速,单位为 ;为降雨量,单位为 。将风浪噪声和雨噪声合成为该水域的噪声干扰:()()海洋混响可分为三种,分别为体积混响、海面混响和海底混响,式()、()、()分别对其进行建模:()()()()()式中,、为体积混响、海面混响和海
7、底混响;、为体积混响强度、海面散射强度和海底散射强度,单位为;为形成混响的体积,单位为;为声速,单位为 ;为发射脉冲宽度;为声呐探测距离,单位为;为发射和接收器组合的等效合成波束宽度,单位为;为产生混响的散射面积,单位为;为等效合成束宽,为掠射角,单位为;为风速,单位为;为频率,单位为。目标强度及辐射噪声声呐主动探测时,不同的舰艇以及舰艇的不同部位对声呐的反射强度会有差异,表 为常见声呐目标强度值,可作为声呐目标探测工作中的目标强度的一般估计值。表 常见舰艇的目标强度值 目标方位 潜艇小型大型(有涂层)大型(无涂层)正横艇首或艇尾其他位置水面舰船正横其他位置被动声呐探测距离除受水文条件的影响外
8、,还受目标参数的影响。舰艇的辐射噪声主要有三个来源,分别为舰艇机械振动产生的机械噪声、水动力噪声和螺旋桨噪声,舰艇排水量的不同和航速的变化是辐射噪声的主要影响因素。舰艇对被动声呐产生的辐射噪声用经验公式计算:()式中,是敌方舰船的航速,单位为;是敌方舰船的排水量,单位为;是敌方舰船的辐射噪声频率,单位为。声呐探测范围仿真设自导水雷布设水域内的干扰主要为各向同性的环境噪声,自导水雷预定的作战时间内,预定海域可能达到的最大风速 ,最大降水量 郭子杰,等:基于遗传算法的自导水雷雷位优化第 卷,海面与海水的最大温差 。首先计算水域的环境噪声干扰,水域环境噪声 随频率变化如图 所示。图 海洋环境噪声随频
9、率变化图 设攻击的目标排水量为 ,速度为 ;自导水雷被动声呐的接收指向性指数 ,检测阈值 ,自导水雷被动声呐的工作频率为。计算得到此时目标能被自导水雷的被动声呐探测到的最大距离 。自导水雷开启主动声呐探测目标时,舰艇的声反射强度影响因素主要为舰艇的舷角,按照舰艇对声反射的强度变化规律,可将目标的弦角分为 、四个变化规律大致相同的区域。以 舷角范围为例,根据声反射强度的不同,将其划分为 个区域:)在 舷角范围内,因舰壳和尾流的遮蔽效应,其声反射强度最小,为 ;)在 舷角范围内的声反射强度为 ;)在 舷角范围内的声反射强度为 ;)在 舷角范围内的声反射强度为 ;)在 舷角范围内的声反射强度最大,为
10、。设自导水雷的主动声呐开启时的主要背景干扰为各向同性的环境噪声,主动声呐工作频率与被动声呐相同,计算得到目标在各个区域上能被自导水雷主动声呐探测到的最大距离,见表。该目标在自导水雷作战时间内,能被水雷的主动声呐和被动声呐同时探测到的最小范围如极坐标图图 所示,图中原点 为目标所在位置。表 各区域被探测到的最大距离 区域()探测距离 图 目标可被探测到的最小范围 雷障拦截能力仿真 多目标通过雷障的航迹分类多个目标通过雷障时一般会存在两种情况:)目标航迹相互独立通过若在布雷任务开始之前,未能获取目标的航迹信息,出于封锁海域的作战目标,可以将目标航迹视作在 范围内随机等概率分布,目标航迹分布的概率密
11、度函数可以表示为,()目标以纵队依次通过若在布雷任务之前已经获悉目标的航迹信息,针对此航迹信息进行布雷任务,则可以将首个目标的航迹视作在 范围内服从均值为 的正态分布,一旦首个目标触发水雷后,后续目标按照首个目标的航迹通过雷障,此时后续目标的航迹服从均值为首个目标航迹 的正态分布。目标的航迹分布的概率密度函数可以表示为第 期指挥控制与仿真(),(),()目标通过雷障仿真方法为研究方便,设目标的航向与 轴正方向相同,航迹在 轴上的 范围内,自导水雷布设在某一片长为,宽为 的海域,如图 所示。图 水雷布设区域 雷障对多个目标的有效拦截能力的仿真方法如下:)假设雷障预先设定的拦截目标个数为,前 个目
12、标通过雷障时若被水雷拦截,则会对雷障的完整性造成破坏;)每个目标通过雷障时,一旦自导水雷发射战斗部拦截目标,即视为目标触发了该水雷,该水雷对后续通过的目标不具备拦截能力;)一个目标可以触发多个水雷,目标通过雷障时,只要未被战斗部命中都可以继续航行,可触发其他水雷;由于水下爆炸威力较大,同时被攻击位置为船体较为脆弱的底部和侧面,设定目标被水雷命中后无法继续航行。)通过多次仿真,将 个目标通过被破坏雷障时的平均触雷概率,作为雷障对多个目标的拦截能力。多目标通过雷障的仿真流程图如图 所示。基于遗传算法的雷位优化 遗传算法基本原理遗传算法的算法步骤与自然界的物种进化过程类似:首先对于某个问题的可行解进
13、行编码,形成初始化的种群;接下来,用适应性函数对种群中的个体进行适应度评估;之后用选择函数对种群中的个体按照某种规则进行择优;然后让被选择的个体进行交叉,产生优秀的子代;同时为了满足“基因”的多样性,按照一定的图 多目标通过雷障的仿真流程图 规则让个体进行“变异”。重复以上过程直到达到终止条件。过程如图 所示。雷位优化步骤:生成初始组群。自导水雷布设在长为,宽为 的海域内,共布设 个自导水雷,采用随机均匀布设的方式产生初始种群。水雷的坐标为(,),()()()若考虑布置 个水雷的部署方案,初始种群中个体的个数为,将初始种群记为,每个布雷方案记为,迭代次数记为,则种群 可以表示为 郭子杰,等:基
14、于遗传算法的自导水雷雷位优化第 卷图 遗传算法流程图 ,(,),(,),(,),(,)():计算适应度。设共要拦截 个目标,目标的平均触雷概率为,为个目标通过雷障时触发水雷的目标个数。对 个目标通过雷障的情况仿真 次,求出 个目标通过雷障的平均触雷概率,即为个体适应度。:终止条件判断。终止条件一般为算法运行次数达到预设的迭代次数或者种群中存在个体达到对可行解的要求。若已经达到了终止条件,则在种群中找出最优个体作为遗传算法迭代得到的最优解;否则,继续执行。:选择个体。选择操作一般分为三个部分:首先在原先种群中按照优选概率 选出最优秀的一部分个体,不经过交叉直接放入下一代的种群中;然后在原先种群中
15、按照轮盘赌的方式选择出 ()个个体,将选出来的个体放入交叉池中;最后再次按照轮盘赌的方式选出()()个个体,用来进行变异操作。:交叉。在交叉池中 ()个父代中,按照交叉步长,选出两个个体 和(),把这两个个体中的基因,即各自导水雷的雷位进行交叉,得到下一代的个体,代个体的雷位坐标表示如下:()()()()()()()式中,为交叉系数。:变异。将选出的()()个个体进行高斯变异,代变异个体的坐标表示如下:()(),()()()()式中,为标准值为、方差为 的高斯随机变量。水雷雷位优化仿真结果设预定的拦截目标数量为 个,自导水雷布设区域为长为 、宽为 的长方形水域,自导水雷的初始雷位在该水域内随机
16、产生,自导水雷的布设深度为 ,水雷数量为 枚,初始雷位如图 所示。用遗传算法对自导水雷的雷位进行优化时,将群体大小设为 个,最大迭代次数设为 次,优选概率 ,交叉概率 ,终止条件为多个目标的平均触雷概率达到 或者达到预设的迭代次数。图 自导水雷的初始雷位 目标随机通过雷障的雷位优化结果多个目标以航迹在 内随机分布的方式通过雷障。按照上文设置的参数对雷位进行优化,算法收敛过程如图 所示,图中横坐标为迭代次数,纵坐标为多目标平均触雷概率。优化之前多目标平均触雷概率为,算法稳定后,计算最后 次迭代中目标触雷概率的平均值为 ,多个目标通过雷障时触雷概率的平均值提高了 。将在该水域布设的自导水雷数量增加
17、到 枚后,多目标的平均触雷概率明显提升,优化前触雷概率为,算法收敛过程如图 所示,算法稳定后,计算最第 期指挥控制与仿真 后 次迭代中目标触雷概率的平均值为 ,多个目标通过雷障时触雷概率的平均值提高了 。图 布设 枚水雷目标随机通过触雷概率变化图 图 布设 枚水雷目标随机通过触雷概率变化图 将布设水雷的数量从 枚增加到 枚,多个目标随机通过雷障的平均触雷概率明显提高,优化前后的平均触雷概率如表 所示。表 布设 枚水雷优化前后多目标随机通过时的平均触雷概率 水雷数量 枚优化前触雷概率 优化后触雷概率 目标纵队通过雷障的雷位优化结果多个目标以纵队通过雷障时,设首个目标的航线分布情况如图 所示,后续
18、目标在首目标触雷后按照首目标的航线通过雷障。图 首目标航线分布图 按照上文设置的参数对雷位进行优化,算法的收敛过程如图 所示。优化之前多目标平均触雷概率为 ,算法稳定后,计算最后 次迭代中目标触雷概率的平均值为 ,多个目标通过雷障时触雷概率的平均值提高了 。图 布设 枚水雷目标纵队通过平均触雷概率变化图 为进一步提高雷障的封锁能力,可以增加布设的自导水雷数量。将自导水雷数量增加为 枚后,算法的收敛过程如图 所示,优化之前多目标平均触雷概率为 ,算法稳定后,计算最后 次迭代中目标触雷概率的平均值为 ,多个目标通过雷障时触雷概率的平均值提高了 。将自导水雷数量从 枚递增到 枚,优化前后的多目标平均
19、触雷概率如表 所示。郭子杰,等:基于遗传算法的自导水雷雷位优化第 卷图 布设 枚水雷目标纵队通过平均触雷概率变化图 表 多目标纵队通过时布设 枚水雷优化前后的平均触雷概率 水雷数量 枚优化前触雷概率 优化后触雷概率 由表 和表 可知,目标采用纵队依次通过雷障的平均触雷概率较低。在经过优化后,多目标通过雷障时的平均触雷概率有了明显提升,其中目标纵队通过雷障时的平均触雷概率提升较大,优化之后雷障的拦截能力的提升近似于在雷障中多布设一枚水雷。结束语本文首先对自导水雷的声呐探测系统进行建模,对声呐探测系统的探测范围进行了仿真,然后对目标采用随机航线和编成纵队两种方式,对通过雷障时的平均触雷概率进行仿真
20、。采用遗传算法对自导水雷的雷位进行优化,并对优化前后目标的平均触雷概率进行了对比分析。据此,可以在预定拦截目标数量确定后,选择较为合适的布雷数量与雷位坐标,提高自导水雷组成的雷障对目标的拦截能力。参考文献:颜冰,董理,张晓兵 水雷作战效能分析 北京:国防工业出版社,:,袁茂钱 水雷反潜中探测技术的研究与分析 成都:西南交通大学,:,():刘伯胜,黄益旺,陈文剑 水声学原理 版 北京:科学出版社,:,吴静 海洋环境噪声建模研究 哈尔滨:哈尔滨工程大学,:,王宏,韩明连,陆达人 舰船声尾流自导鱼雷及其防御技术 声学技术,():,():朱伟,冯金鹿,李宵宇 主动声呐探测性能预报分析及软件设计 舰船电子工程,():,():,():,:,():秦长江,黄亭飞,黄金才 改进遗传算法的反无人机作战火力分配优化 国防科技,():,():(责任编辑:许韦韦)