海上舰队的防空编队建模与防空能力评估
刘 昕
(中北大学 大数据学院，山西 太原 030051)
摘 要： 海上舰艇编队巡逻是维护我国领海主权的有效方法，其编队队形决定着自身的防空能力，也影响着各种突发事件的应变能力. 在信息化战争技术迅速发展的今天，国防安全对海上舰队防空编队方案的制定提出更高的要求. 本文通过设定合理的海上巡逻任务，以及舰艇编队和来袭导弹的参数和指标，充分考虑各舰之间的数据共享和协同作战特点，在分析传统舰队编队队形特点的基础上，以使在舰队护卫舰防御范围最薄弱处的位置与指挥舰的距离最长为目标，建立并求解舰队巡逻的最佳队形编队，然后通过MATLAB软件对设定场景进行仿真实验，并对所建立的最佳舰队队形的抗饱和攻击能力进行评估，该项工作可为海上舰队防空编队的建模方法提供参考依据.
关键词： 海上舰艇编队； 多目标非线性规划； 机动汇算模型； 抗饱和攻击能力
0 引 言
在当前南海、 东海等海域面临严峻挑战的复杂形势下，海上舰艇是维护我国领海主权的最有效工具之一，其自身的对空防御能力是执行巡航任务的最基本保障[1]. 如果没有良好的防空反导能力，海上舰队的自身安全就无法保障，更不用说完成作战任务. 尽管现代舰载防空武器性能已得到较大提升，但单舰的防空能力仍然十分有限，而舰艇编队可利用其共享数据链实现全方位协同侦查和拦截，因此研究海上舰艇编队模型及其防空能力具有重要的现实意义[2-4]. 目前对于海上舰队的建模主要集中在研究防空作战中目标分配[5-7]、 通信过程[8]等方面，对舰队防空队形的建模研究较少. 虽然有一些文献对传统编队队形(如梯形、 菱形等)进行定性分析[4,9]，但缺少对其防空能力的定量评估，难以客观表示其作战效果. 李浩等人[3]虽然基于模糊层次分析法建立了舰艇防空队形的综合评价模型，但并没有给出该队形抗饱和攻击能力的评估. 根据以上海上舰队防空编队研究发展现状，本文通过分析传统编队队形特点，提出使来袭导弹在被发现时刻的位置与指挥舰的距离最长的目标，建立最佳海上舰队编队模型，同时对该模型下的最佳舰队队形进行抗饱和攻击能力评估，旨在发掘该舰队的最佳防空防御能力.
1 问题设定
本文根据现代海战和武器装备的特点，设定合理的海上巡逻任务[10]，据此任务建立海上舰队巡逻时的最佳防空编队模型，并对最佳防空队形的抗饱和攻击能力开展仿真实验评估. 巡逻任务设计如下：由1艘导弹驱逐舰(指挥舰)和4艘导弹护卫舰组成的我海上舰队正在南海海域执行巡逻任务，航向200°(以正北为0°，顺时针方向)，航速16节. 情报显示：在以驱逐舰为原点的20°～220°的扇形面内等可能地有导弹来袭，飞行速度0.9 Ma，射程230 km，导弹探测距离为30 km，搜索扇面为30°，且每批来袭导弹的数量不超过4枚. 我海上舰队各舰发现来袭导弹的概率服从均匀分布，均匀分布的范围是导弹与该舰之间距离在20～30 km，且各舰之间信息共享； 舰载防御导弹水平射程为10～80 km，速度2.4 Ma，从发现来袭导弹到发射拦截导弹的准备时间为7 s，各舰单次仅执行对单批来袭导弹的拦截任务[10].
围绕人才培养的根本任务，作为大学教学管理职能实施与运行的主体——大学教学管理组织，对结构、人员、职权及运作等一系列内容进行组织、协调、控制，使之在教学活动中得到最充分的发挥，促使教学目标的最终实现。因此，从大学教学管理组织的机构来看，笔者较为赞同，大学教学管理是一个多层次的系统，在大学内部，涉及学校、学院、系、教研室等不同层次机构，不同层次的组织机构具有不同的管理职责与权利分工[3]。大学教学管理组织既是静态的专门性组织机构，又发挥动态的管理效能，是保证大学教学活动得以正常开展的基本条件。
2 海上舰队最佳防空队形模型的建立与求解
2.1 问题分析与建模
2.1.1 传统舰队编队阵型特点分析
定义10 理想制造能力(MCI)指在企业一般的生产资料基础上，没有生产性约束条件下，所能达到的最大制造能力水平。
如图 1 所示，结合本文设定的问题情境：① 单边阵型在舰队所在直线上有较强的攻防能力，但在两侧防守能力较为薄弱； ② 双边阵型中若存在指挥舰，会使指挥舰暴露在阵型顶端或边缘，对于区域防守能力较差； ③ 多边阵型属于四周全面防护的阵型，但根据本文所设定问题，该编队阵型不能发挥在20°～220°范围内的最佳防御能力.
图 1 传统舰队编队阵型
Fig.1 Traditional fleet formation
传统小规模舰队编队多会采用单边、 双边、 多边等阵型进行日常巡航任务，如图 1 所示.
李站长说，我这里好说，顶多搞点罚款过去，要是被市林业公安的发现了，可能我也脱不了干系的。何泽说，李站长尽管放心，市里的工作，由买主去做，负责没有障碍。
2.1.2 本文提出的最佳海上舰队编队模型
综合图 3, 图 4 以及上述分析过程，可知该最佳海上舰队编队模型的抗饱和攻击能力为
图 2 舰队大致阵型方位图
Fig.2 The approximate formation azimuth map of the fleet
由于来袭导弹是超低空来袭，高度很低，远小于导弹射程，可认为水面防空是在同一个平面内，故以指挥舰为原点，正东方向为x 轴，正北方向为y 轴建立坐标系. 为方便描述，将指挥舰编号为0，其他四艘护卫舰顺时针方向分别编号为1，2，3和4. 根据上述分析可知l =d =20 km，受防空导弹射程限制，则10 km≤R d ≤80 km，则有目标函数
f (θ 1,θ 2,θ 3,θ 4,R 1,R 2,R 3,R 4)=maxR d .
为了更清晰地描述计算机仿真过程，程序中详细的逻辑关系与参数表达如下：
图 4 中O 点表示指挥舰所处位置，M 0表示敌军导弹初次被发现时的位置，D 表示OM 0的距离，即敌军导弹最初被发现时距离指挥舰的长度，也是最危险方向上第一次可侦查到来袭导弹时距指挥舰的距离.R i 表示各护卫舰距指挥舰的距离.P i 表示第i 号舰所处位置，N n 表示我方拦截敌方导弹第n 个拦截批次. 设δ i 表示OM 0与M 0P i 间的夹角，β i 表示OM 0与OP i 间的夹角，T in 表示第i 号舰拦截第n 次来袭导弹所需时间； 另设防空导弹准备时间为t ，舰载防空导弹发射速度为V 防，来袭导弹的飞行速度为V 攻. 结合图1和图2，以防空导弹与来袭导弹相遇碰撞为成功拦截一批，可继续进入下一次拦截. 抗饱和攻击能力f (n )，指在以上最佳海上舰队编队模型情况下，在220°攻击方向上按最佳防御方案拦截最多批次的导弹数量的防御能力，当导弹飞行距离大于等于OM 0即表示击中指挥舰，则防御结束. 此外，当导弹与舰艇间距离R i ≤10 km或R i ≥80 km，该舰艇无法进行防御，每批次拦截只有一艘舰艇进行防御.
2.2 模型求解
根据上述问题情景设置以及最佳海上舰队编队模型与约束条件，通过MATLAB对上述模型进行求解，得到各舰的部署位置如表 1 所示，各舰的空间分布如图 3 所示.
1384年，世界上第一份具有现代意义的保险单在意大利诞生，这份保险单承保的是，一批货物从法国南部的阿尔兹运到意大利的比萨。在这张保单中，有明确的保险标的、明确的保险责任等内容。世界上最早成立的证券交易所是荷兰阿姆斯特丹证券交易所，创立于1602年。
表 1 编队最佳队形各护卫舰的部署位置(以指挥舰为原点)
Tab.1 Deployment location of frigates in the best formation
(starting from command ships)
图 3 海上舰队最佳防空编队队形空间图
Fig.3 Spatial chart of the best air defense formation of marine fleet
3 编队防空能力评估分析与仿真实验
3.1 防空能力评估分析
由于舰载防空导弹拦截敌方来袭导弹问题属于最短路径相遇问题，可采用最短时间接近的机动绘算模型[13]进行计算，护卫舰的防空反导绘算原理如图 4 所示.
对于海上舰队的防空能力，主要从其抗饱和攻击能力方面来衡量，即我指挥舰在短时间内遭受高密度连续多批次的导弹袭击时舰队所能成功拦截来袭导弹的最大批次数[2,12]. 由于我编队各舰防空导弹射程为10～80 km，且来袭导弹可能从20°～220°的任何方向攻击，采用保守估计原则，以前述最佳编队防御能力最薄弱方向上 10～80 km 范围内的最大拦截来袭导弹批数作为编队的防空能力评价指标. 由图1中舰艇编队中各舰空间位置分布可知，当来袭导弹从20°或220°方向突防时，编队中参与防御的舰艇数最少，所以这两个方向为编队最危险的防御方向，可按其中一个方向来计算敌方导弹袭击指挥舰时的各舰对来袭导弹的防御能力.
图 4 护卫舰的防空反导绘算原理示意图
Fig.4 The schematic diagram of the frigate’s anti-aircraft and anti-missile mapping principle
回到汉武帝时代，元朔元年（前128年）、二年（前127年），匈奴人在东线攻势凶猛，多个郡县告急。但刘彻认为匈奴人活跃的河套平原始终威胁着帝都的安全，决定不去管东线的匈奴，先全力解决河套问题。
对于海上舰队来说，指挥舰相当于整个舰队的“大脑”，要尽可能最大限度地保护其安全，因此一般将指挥舰部署在安全防御范围的核心位置，将4艘护卫舰部署在指挥舰周围承担掩护任务[11]. 考虑到敌方导弹在20°～220°的扇形面内等可能来袭，为最大限度地保护指挥舰的安全，需要保证来袭导弹在20°～220°的任意方向上袭击都先被外围护卫舰发现. 对于舰队的侦测能力方面，应该尽早发现来袭导弹，为下一步行动提供足够的预警时间，即在舰艇编队的侦测最薄弱环节位置处据距离指挥舰的距离最远，且保证护卫舰一定能侦察到来袭导弹，即护卫舰之间最小监测范围应相切. 设护卫舰与驱逐舰间的距离为R i (i =1,2,3,4 ),护卫舰相对驱逐舰的方位角为θ i (i =1,2,3,4 )，舰艇编队的侦测最薄弱环节位置处据距离指挥舰的距离为R d ，此外，由于导弹在扇形区域内为等可能来袭，且各个护卫舰以及主舰的防御能力均相同，则根据以上分析可知舰队大致阵型方位如图 2 所示.
在实习期间学生的安全问题是校企双方最为担心的问题，应制定可行的规章制度来防范安全事故并承担相应的风险。学生在俱乐部既不是正式员工也不是纯粹的实习学生，应该是具有双重身份。例如学生在机械操作中意外受伤，学生在场地服务时碰坏球车等问题如何解决？经过校企双方共同协商，通过给学生全员购买意外伤害险和加强学生安全培训使这个问题得以解决，使“校企一体”教学模式顺利实施。
f (n )=maxn
约束条件为
3.2 饱和攻击仿真实验设计
根据3.1中的分析，假设某次敌方来袭导弹群从编队防御最薄弱的方向(220°)向指挥舰突防，每批来袭导弹4枚，各批次间隔时间极短可忽略，导弹批数足够将指挥舰摧毁； 根据设计的最佳编队队形，只要敌方导弹来袭方向确定，来袭导弹被发现的位置就可确定，各舰就会根据来袭导弹的位置进行自动拦截. 由于单批次内来袭导弹间距较小，我方舰艇对来袭导弹的拦截以批次为单位，且一旦相遇即视为成功拦截. 根据以上分析并通过MATLAB编程进行仿真实验，程序流程图如图 5 所示.
约束条件为
1) 舰队模型参数初始化是指初始化舰队队形参数(各护卫舰相对于指挥舰的空间位置，见表1)和舰载防空导弹参数(飞行速度、 攻击范围20～80 km 等).
2) 输入来袭导弹群的信息是指来袭导弹的飞行速度、 批数和方位(与正北方向的夹角)等.
3) 根据最佳编队模型计算在来袭导弹方向上发现来袭导弹的预警距离，即来袭导弹群与指挥舰之间距离.
图 5 仿真实验程序流程图
Fig.5 Flow chart of simulation experiment program
4) 根据各舰艇坐标计算当前来袭导弹群与各舰艇之间的空间位置关系，即图 4 所示反导绘算模型.
5) 依次判断各舰是否可以成功拦截当前位置的一批来袭导弹； 如果有舰艇能拦截当前位置来袭导弹，选取可最快拦截导弹的舰艇进行拦截，更新各舰成功拦截的来袭导弹批数和拦截位置，并返回第4)步.
6) 如果没有舰艇能拦截当前位置来袭导弹，则意味着到达最大抗饱和攻击能力，指挥舰将被来袭导弹击中摧毁，此时输出各舰成功拦截的来袭导弹批次和拦截位置.
3.3 实验结果
输入来袭导弹群的飞行速度、 导弹批数(100，远大于舰载防空导弹的防空能力)、 袭击方向(220°)和各批次拦截时间. 程序输出各舰成功拦截的来袭导弹批次和拦截位置如表 2 所示.
二是尽快全面摸清震损水利设施情况。会同当地水利部门组织技术力量全面开展险情排查，逐库、逐站、逐处、逐段、逐项摸清水库、水电站、河道堤防、供水设施、灌溉工程、水文测报等水利设施的震损情况，全面掌握震损工程的来蓄水量、地质构造、工程结构以及可能的影响范围，深入灾区了解农村人畜饮水设施损毁情况，逐村逐户核查饮水困难人口，组织专家分析研判是否存在重大地质灾害、堰塞湖等次生灾害风险。
表 2 各舰在最危险方向上拦截导弹的信息
Tab.2 Information of the warships intercepting missiles in the most dangerous direction
综上可得，对于220°方向(防御最薄弱的方向)上的来袭导弹，1～3号护卫舰分别只可拦截1批； 4号护卫舰可拦截2批来袭导弹； 指挥舰(0号)可拦截3批，编队整体舰载防空导弹总共可拦截8批. 因此，建立的最佳编队的抗饱和攻击能力为8批.
4 结 论
针对基于共享数据链的海上多舰艇编队问题，本文通过设定合理的海上巡逻任务，以编队的最大防御面积为目标，采用多目标非线性规划方法对最佳巡逻编队队形开展了建模研究，求解的最佳编队队形为四艘护卫舰分别均匀部署在指挥舰周围，以指挥舰为原点的方向，与指挥舰的舰间距为47.32 km. 然后通过MATLAB软件进行仿真实验，对所建立的最佳舰队队形的抗饱和攻击能力进行评估，仅依赖舰载防空导弹在最危险的方向上可拦截8批次来袭导弹，该项工作可为海上舰队防空编队的建模方法提供参考依据.
参考文献：
[1] 程少. 从木船到航母：世界瞩目的中国海军水面舰艇部队[J]. 百科探秘(海底世界)，2015(3)：23-27.
[2] 韩庆，李振冲. 谈谈舰船编队的队形[J]. 现代舰船，1995(7)：13-16.
[3] 李浩，王公宝. 关于水面舰艇编队防空队形的优化研究[J]. 舰船科学技术，2009，5(31)：109-112.
[4] 高善民. 水面舰艇编队队形[J]. 海军大连舰艇学院学报，2000，17(2)：12-17.
Gao Shanmin. Formation of surface ships[J]. Journal of Dalian Naval Ship College，2000，17(2)：12-17. (in Chinese)
[5] 赵国钢，孙永侃，王立强. 驱护舰编队防空火力分配决策模型研究[J]. 现代防御技术，2007，35(1)：45-48.
Zhao Guogang，Sun Yongkan，Wang Liqiang. Study on antiaircraft firepower allocation decision-making model of destroyer and frigate formation[J]. Modern Defence Technology, 2007, 35(1): 45-48. (in Chinese)
[6] 姚跃亭，赵建军，尹波波，等. 舰艇编队防空目标分配优化算法研究[J]. 计算机与数字工程，2011，39(1)：31-34.
Yao Yueting，Zhao Jianjun, Yin Bobo, et al. Algorithms of target allocation in naval fleet air defense[J] Computer and Digital Engineering，2011，39(1)：31-34. (in Chinese)
[7] 姜鲁东，余家祥，斗计华. 水面舰艇编队区域防空导弹目标分配模型研究[J]. 现代防御技术，2012，40(3)：29-33.
Jiang Ludong，Yu Jiaxiang，Dou Jihua. Air defense missile target allocation models for surface warship formation[J]. Modern Defence Technology，2012，40(3)：29-33. (in Chinese)
[8] 刘开第，庞彦军，郝继梅. 基于信息融合技术的防空通信系统效能评估[J]. 火力与指挥控制，2010，35(5)：5-9.
Liu Kaidi，Pang Yanjun，Hao Jimei. Effectiveness evaluation of the air defense communication systems based on information fusion technology[J]. Fire Control & Command Control，2010，35(5)：5-9. (in Chinese)
[9] 赵建印，刘芳. 海上编队双舰防空队形部署优化研究[J]. 兵工学报，2010，31(6)：865-869.
Zhao Jianyin，Liu Fang. Optimization of two warships formation for air defense[J]. Acta Armamentarii，2010，31(6)：865-869. (in Chinese)
[10] 2015年全国研究生数学建模竞赛A题[EB/OL]. (2016-01-04)[2019-03-07]. http:∥wenku.baidu.com/view/9d51cec9de80d4d8d15a4fc9.html
[11] 李永刚. 舰艇编队对空防御队形优化配置模型[J]. 指挥控制与仿真，2004，26(4)：65-67.
Li Yonggang. A model of optimal disposition of warship formation for air defense[J]. Information Command Control System & Simulation Technology，2004，26(4)：65-67. (in Chinese)
[12] 张训立，杨承军. 弹道导弹防御系统效能分析与措施研究[J]. 现代防御技术，2002，30(5)：9-12.
Zhang Xunli，Yang Chengjun. Efficiency analysis and measure study of ballistic missile defen ce system[J]. Modern Defence Technology，2002，30(5)：9-12. (in Chinese)
[13] 黄文雅，童雪娟. 海警舰艇接近机动自动绘算软件的设计与实现[J]. 公安海警学院学报，2015, 14(3)：76-78
Huang Wenya，Tong Xuejuan. Design and implementation of coast guard ship close maneuvering automatic plotting software[J]. Journal of China Maritime Police Academy，2015, 14(3)：76-78. (in Chinese)
Air Defense Formation Modeling and Air Defense Capability Evaluation of Sea Fleet
LIU Xin
(School of Data Science and Technology, North University of China, Taiyuan 030051, China)
Abstract ：Maritime patrol of naval vessels formation is an effective way to maintain the national territorial sovereignty. Its formation determines the air defense capability as well as affects the resilience of various emergencies. Today, with the rapid development of information warfare technology, national defense security puts forward higher requirements for the formulation of the air defense air defense formation plan for the sea fleet. By setting reasonable maritime patrol missions, as well as the parameters and indicators of ship formations and incoming missiles, this paper fully considered the characteristics of data sharing and coordinated operations between ships, based on the analysis of the formation characteristics of traditional fleet formations, so that at the weakest position of the fleet frigate defense area, the distance from the command ship was the longest, and the best formation of the fleet patrol was established and solved. Then the simulation experiment of the set scene was carried out by MATLAB software, and the most established the anti-saturation attack capability of the formation of the fleet was evaluated, and this work can provide a reference for the modeling method of the air defense formation of the naval fleet.
Key words ：naval vessel formation; multi-objective nonlinear programming; maneuvering model; anti-saturation attack capability
中图分类号： E911
文献标识码： A
doi： 10.3969/j.issn.1673-3193.2019.06.008
文章编号： 1673-3193(2019)06-0531-06
收稿日期： 2019-03-07
作者简介： 刘 昕(1989-)，男，博士生，主要从事复杂系统建模与仿真研究.
标签：海上舰艇编队论文; 多目标非线性规划论文; 机动汇算模型论文; 抗饱和攻击能力论文; 中北大学大数据学院论文;