CN114036783A - 基于截断机理的线目标功能毁伤判断方法及毁伤计算模型 - Google Patents
基于截断机理的线目标功能毁伤判断方法及毁伤计算模型 Download PDFInfo
- Publication number
- CN114036783A CN114036783A CN202210009447.XA CN202210009447A CN114036783A CN 114036783 A CN114036783 A CN 114036783A CN 202210009447 A CN202210009447 A CN 202210009447A CN 114036783 A CN114036783 A CN 114036783A
- Authority
- CN
- China
- Prior art keywords
- damage
- line
- circle
- damaged
- line target
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于军事运筹领域,特别涉及目标毁伤计算方法。其技术方案提出一种基于截断机理的线目标功能毁伤判断方法,从功能毁伤层面上对线目标进行毁伤判断,并且巧妙通过构造线目标的两侧长边与各个毁伤圆的相互相交关系矩阵,寻找两侧长边之间的毁伤区域通路,判断线目标是否在长方向功能流向上被截断,从而为线目标毁伤运筹提供一种新的思路方法。同时提出一种基于截断机理的线目标功能毁伤计算模型,利用蒙特卡洛仿真计算方法,并按提出的基于截断机理的线目标功能毁伤判断方法,计算各次仿真线目标是否功能丧失,得到线目标功能毁伤概率。本方法相对于传统方法大辐减少了线目标的实际物理毁伤长度,减少了武器使用数量。
Description
技术领域
本发明属于军事运筹领域,特别涉及目标毁伤计算方法。
背景技术
目标毁伤计算是指对目标遭受武器打击后的效果进行预测,通常用目标达到一定毁伤要求的概率来描述。通常对长宽比较大的矩形目标或功能流向为线性的道路、桥梁等交通类目标作为线目标。线目标在现实中,不仅仅体现为简单的直线段,同时存在大量由多段直线段目标相连的折线目标。对于这类目标,典型的计算方法是采用平均相对毁伤长度来表示对其的毁伤效果(参考文献《交通重点目标毁伤模拟计算》,国防交通工程与技术2010年6期),即武器毁伤圆半径对线目标长度上的覆盖长度比。该方法主要是从物理长度上对线目标毁伤情况进行描述,适合于毁伤圆半径大于线目标宽度的大威力武器,对于小威力武器,即武器毁伤圆直径相比线目标宽度较小的情况,导致毁伤圆难以完全覆盖线目标宽方向,若仍采用平均相对毁伤长度进行毁伤计算,在武器的命中误差的影响下,小威力武器难以实现对线目标长方向上的毁伤意图,所需武器数量也将非常多。
发明内容
本发明的目的是针对小威力武器打击以功能流向为主的线目标,提供一种基于截断机理的线目标功能毁伤判断方法及毁伤计算模型,实现以较少弹量达到线目标功能丧失的毁伤意图。
为实现上述发明目的,本发明的发明构思是:本发明从功能毁伤层面上对线目标进行毁伤判断;并且巧妙通过构造线目标的两侧长边与各个毁伤圆的相互相交关系矩阵,寻找两侧长边之间的毁伤区域通路,判断线目标是否在功能流向上被截断,从而为线目标毁伤运筹提供一种新的思路方法。
本发明的技术方案是:一种基于截断机理的线目标功能毁伤判断方法,其特征是:
定义:所述线目标的功能流向边界为长边;武器落点位置为圆心,与毁伤半径构成毁伤圆;
根据明确的武器随机落点位置及毁伤圆大小,通过所述毁伤圆与所述线目标功能流向两侧长边相交情况与所述毁伤圆之间相互相交情况,构造线目标的两侧长边与各个毁伤圆的相互相交关系矩阵,寻找两侧长边之间是否存在由毁伤圆连接的通路;若存在所述通路则所述线目标功能流向被截断,所述线目标功能丧失;反之,所述线目标功能未丧失。
进一步的,对于所述线目标两条长边为折线时,将位于所述功能流向同侧折线长边的各直线段与所述毁伤圆相交情况的并集,作为该侧折线长边与所述毁伤圆的相交情况。
进一步的,上述具体步骤为:
第一步,计算长边与毁伤圆相交情况
第二步,计算各毁伤圆相互相交情况
第三步,构建相交情况矩阵H,H为M×(M+2)矩阵,其中第1列、第(M+2)分别表示M个毁伤圆与a、b长边的相交情况、第2至第(M+1)列分别各毁伤圆之间的相交情况;
第四步,根据相交关系矩阵,寻找两侧长边之间是否存在由毁伤圆连接的通路,判断线目标是否被截断;若存在所述通路则所述线目标功能流向被截断,所述线目标功能丧失;反之,所述线目标功能未丧失。
进一步的,所述第四步中寻找是否存在通路采用欧拉通路算法。
进一步的,判断所述毁伤圆与所述线目标功能流向两侧长边相交情况与所述毁伤圆之间相互相交情况时,考虑一定允许的容差。
进一步的,一种基于截断机理的线目标功能毁伤计算模型,其特征是:利用蒙特卡洛仿真计算方法,通过多次仿真,按照武器二维正态分布命中规律,生成各次仿真M枚随机落点的位置及毁伤圆大小,并按上述基于截断机理的线目标功能毁伤判断方法计算各次仿真线目标是否功能丧失;进而统计线目标功能丧失的次数,得到线目标功能毁伤概率。
本发明的有益效果:本发明提出的一种基于截断机理的线目标功能毁伤判断方法,从功能毁伤层面上对线目标进行毁伤判断,并且巧妙通过构造线目标的两侧长边与各个毁伤圆的相互相交关系矩阵,寻找两侧长边之间的毁伤区域通路,判断线目标是否在功能流向上被截断,从而为线目标毁伤运筹提供一种新的思路方法;而且该思路方法可适于直线目标、折线目标,甚至可扩展到任何线目标。在该思路方法下,对线目标的毁伤转化为只需要对线目标中的一段或多个子段的截断毁伤来实现,相对于传统平均相对毁伤长度方法而言,大辐减少了线目标的实际毁伤长度,从而减少了武器使用数量,特别是征对于小威力武器打击线目标问题解决优越性特别明显。
附图说明
图3为本发明实施例1中相交情况矩阵H以图的方式展现节点间相连关系的示意图;
图4为本发明实施例1中复杂折线目标两侧边示意图;
图5为本发明实施例2中线目标与已知的武器落点毁伤圆位置示意图;
图6为本发明实施例2中相交情况矩阵H以图的方式展现节点间相连关系的示意图。
具体实施方式
下面结合附图和具体实施例对本发明的技术方案作进一步具体的说明。
以下结合附图1,进一步说明本发明的具体实施方式。本发明不限于以下实施例的描述。
实施例1:本实施例给出一种具体的示例。
第一步,按照武器命中落点随机误差规律,生成每枚武器的落点位置,结合其毁伤半径,生成对应毁伤圆。
第二步,计算每个毁伤圆是否与两条长边相交。考虑到线目标可以为非直线的折线目标,在计算时,可将折线目标的每一段直线段分别进行处理。具体为:
(2)分别针对a方向长边各个直线段i,计算线段所在直线方程。直线段i的a方向长边线段的端点分别为()和(),构建线段所在直线方程Ax+By+C=0,其中,A=,B=,C=。同理计算b方向长边的直线方程。因为计算得到的各个直线方程在整个计算中不会变化,因此只需计算一次即可,后续计算均可使用第一次计算得到的直线方程。
(3)循环计算每个毁伤圆j,计算其与a方向长边直线段i所在直线Ax+By+C=0的相交情况。本发明使用毁伤圆心(到直线Ax+By+C=0的距离d是否小于等于毁伤圆半径与允许的容差距离之和作为是否相交的判断依据。
若直线与毁伤圆相交,进一步判断计算a方向长边直线段i是否被毁伤圆j(圆心为(,半径由变更为)覆盖。具体方法为:利用解析几何方法直接计算直线Ax+By+C=0与毁伤圆j解析公式 的交点()和(),交点坐标的具体解算方法为公知常识。其中,若直线与毁伤圆相切,则两个交点相同。
因为直线段()->()与上述交点形成的直线段()->()在同一直线上,若两个直线段存在交集,说明毁伤圆j与直线段i真实相交,反之说明未相交。两个直线段是否存在交集的判断有很多,如只要存在某线段任一端点与另一线段两个端点的矢量方向是否相反来判断,若相反,说明相交;反之,说明不相交。
在真实相交情况下,记录毁伤圆j与a方向长边的交差情况,可不再继续计算毁伤圆j与其他直线段相交情况,即只要有相交即可;若不相交,再继续计算毁伤圆j与直线段(i+1)的相交情况。也就是对于折线目标的a、b方向各长边的任一折线段与毁伤圆j相交,则认为该方向长边与毁伤圆j相交。
(4)若不相交,则继续计算毁伤圆j与a方向长边其他直线段的相交情况。
(5)同理,计算毁伤圆j与b方向长边相交情况。
第四步,构建相交情况矩阵H。H维数为M×(M+2)(行×列)。其形式为:
该矩阵也可表示为图,如图3所示相交情况矩阵H以图的方式展现节点间相连关系的示意图,图3中各个节点表示各个毁伤圆,若毁伤圆相交,则两个节点间存在连线。三角形节点表示与a方向长边相交的节点,星形节点表示与b方向长边相交的节点,圆形节点表示未与a、b方向长边相交的节点,方形节点表示同时与a、b方向长边相交的节点。
第五步,求解任一三角形节点与任一星形节点之间是否存在通路。求解方法可使用图论中的欧拉通路相关算法求解。若存在通路,表示线目标被截断,无法通行。若不存在通路,表示线目标未被(彻底)截断。若存在方形节点,即同时与a、b方向长边相交的节点,表示线目标被截断。
第六步,利用蒙特卡洛仿真计算方法,对步骤一到步骤五进行多次仿真计算,得到线目标被截断的毁伤概率。
本发明并不要求a、b两个方向的长边线段数量完全相同,同样可用于多边形描述的线目标,包含a、b方向长边线段数量不相同,线目标各段宽度不相同的情况。如图4所示复杂折线目标两侧边示意图,其中红色实线是一个方向的长边,由三条直线段构成,红色虚线是另一个方向的长边,由四条直线段构成。在处理时,只需针对各个方向长边各自的线段数量分别进行处理即可。该思路方法可适于直线目标、折线目标,甚至可扩展到任何线目标。
实施例2:
为了进一步说明本发明的技术方案,用一个具体目标来呈现截断判断方法。
如图5所示线目标与已知的武器落点毁伤圆位置示意图,折线目标中心线为(0,10)->(50,10)->(70,20),宽为10米。
有10个毁伤圆,落点坐标分别为(31,15)、(37,14)、(37,10)、(39,7)、(42,11)、(45,11)、(47,12)、(48,7)、(50,5)、(56,13)。
a边可分为两个线段,端点分别为(0,5)->(50,5)和(50,5)->(70,15),对应直线方程分别为0x-50y+250=0和10x-20y-400=0。
b边可分为两个线段,端点分别为(0,15)->(50,15)和(50,15)->(70,25),对应直线方程分别为0x-50y+750=0和10x-20y-200=0。
10个毁伤圆与4条线段的相交情况为:
10个毁伤圆相互两两相交情况为:
构建的相交情况矩阵H如下:
对应的相交关系示意图如图6所示。
通过欧拉通路解算可知,存在两条通路,一条是2->3->4,如图6中红色实线路径。一条是2->3->4->5->6->8->9,如图6所示相交情况矩阵H以图的方式展现节点间相连关系的示意图,图中红色实线+虚线路径。
由此可知,当前实例中线目标被截断。
Claims (6)
1.一种基于截断机理的线目标功能毁伤判断方法,其特征是:
定义:所述线目标的功能流向边界为长边;武器落点位置为圆心,与毁伤半径构成毁伤圆;
根据明确的武器随机落点位置及毁伤圆大小,通过所述毁伤圆与所述线目标功能流向两侧长边相交情况与所述毁伤圆之间相互相交情况,构造线目标的两侧长边与各个毁伤圆的相互相交关系矩阵,寻找两侧长边之间是否存在由毁伤圆连接的通路;若存在所述通路则所述线目标功能流向被截断,所述线目标功能丧失;反之,所述线目标功能未丧失。
2.如权利要求1所述的一种基于截断机理的线目标功能毁伤判断方法,其特征是:对于所述线目标两侧长边为折线时,将位于所述功能流向同侧折线长边的各直线段与所述毁伤圆相交情况的并集,作为该侧折线长边与所述毁伤圆的相交情况。
3.如权利要求1所述的一种基于截断机理的线目标功能毁伤判断方法,其特征是:具体步骤为:
第一步,计算长边与毁伤圆相交情况
第二步,计算各毁伤圆相互相交情况
第三步,构建相交情况矩阵H,H为M×(M+2)矩阵,其中第1列、第(M+2)分别表示M个毁伤圆与a、b长边的相交情况、第2至第(M+1)列分别各毁伤圆之间的相交情况;
第四步,根据相交关系矩阵,寻找两侧长边之间是否存在由毁伤圆连接的通路,判断线目标是否被截断;若存在所述通路则所述线目标功能流向被截断,所述线目标功能丧失;反之,所述线目标功能未丧失。
4.如权利要求3所述的一种基于截断机理的线目标功能毁伤判断方法,其特征是:所述第四步中寻找是否存在通路采用欧拉通路算法。
5.如权利要求1所述的一种基于截断机理的线目标功能毁伤判断方法,其特征是:判断所述毁伤圆与所述线目标功能流向两侧长边相交情况与所述毁伤圆之间相互相交情况时,考虑一定允许的容差。
6.一种基于截断机理的线目标功能毁伤计算模型,其特征是:利用蒙特卡洛仿真计算方法,通过多次仿真,按照武器二维正态分布命中规律,生成各次仿真M枚武器随机落点的位置及毁伤圆大小,并按权利要求3所述的一种基于截断机理的线目标功能毁伤判断方法计算各次仿真线目标是否功能丧失;进而统计线目标功能丧失的次数,得到线目标功能毁伤概率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210009447.XA CN114036783B (zh) | 2022-01-06 | 2022-01-06 | 基于截断机理的线目标功能毁伤判断方法及毁伤计算模型 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210009447.XA CN114036783B (zh) | 2022-01-06 | 2022-01-06 | 基于截断机理的线目标功能毁伤判断方法及毁伤计算模型 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114036783A true CN114036783A (zh) | 2022-02-11 |
CN114036783B CN114036783B (zh) | 2022-03-22 |
Family
ID=80147331
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210009447.XA Active CN114036783B (zh) | 2022-01-06 | 2022-01-06 | 基于截断机理的线目标功能毁伤判断方法及毁伤计算模型 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114036783B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114419450A (zh) * | 2022-03-29 | 2022-04-29 | 中国人民解放军96901部队 | 一种基于影像特征分析的线性目标毁伤效能快速评估方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101220185B1 (ko) * | 2012-04-30 | 2013-01-21 | 한국과학기술원 | 무선통신 중계를 위한 주행 로봇용 중계기 발사장치 |
CN104143028A (zh) * | 2014-08-12 | 2014-11-12 | 中国人民解放军装甲兵工程学院 | 一种基于仿真实验的装甲装备毁伤规律分析方法 |
DE102013219056A1 (de) * | 2013-09-23 | 2015-03-26 | Fico Cables Lda | Verfahren zur herstellung einer stützmatte |
CN107588685A (zh) * | 2017-09-21 | 2018-01-16 | 中国舰船研究设计中心 | 一种确定悬浮式拦截弹最佳布阵间距的方法 |
CN107957319A (zh) * | 2017-11-17 | 2018-04-24 | 湘潭大学 | 均匀荷载面曲率的简支梁裂纹损伤识别方法 |
CN108666374A (zh) * | 2018-05-18 | 2018-10-16 | 通威太阳能(安徽)有限公司 | 一种背面钝化矩阵点式激光开槽导电结构 |
CN110119580A (zh) * | 2019-05-17 | 2019-08-13 | 中国人民解放军军事科学院国防工程研究院 | 一种地面工程目标毁伤评估系统及方法 |
CN110808863A (zh) * | 2019-11-08 | 2020-02-18 | 中国人民解放军陆军工程大学 | 基于区域毁伤的基础设施网络抗毁性评估方法和模型 |
CN111090931A (zh) * | 2019-12-05 | 2020-05-01 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种面目标区域打击规划方法 |
CN112685944A (zh) * | 2021-01-07 | 2021-04-20 | 北京理工大学 | 一种同类多枚杀爆弹对面目标瞄准点规划方法 |
-
2022
- 2022-01-06 CN CN202210009447.XA patent/CN114036783B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101220185B1 (ko) * | 2012-04-30 | 2013-01-21 | 한국과학기술원 | 무선통신 중계를 위한 주행 로봇용 중계기 발사장치 |
DE102013219056A1 (de) * | 2013-09-23 | 2015-03-26 | Fico Cables Lda | Verfahren zur herstellung einer stützmatte |
CN104143028A (zh) * | 2014-08-12 | 2014-11-12 | 中国人民解放军装甲兵工程学院 | 一种基于仿真实验的装甲装备毁伤规律分析方法 |
CN107588685A (zh) * | 2017-09-21 | 2018-01-16 | 中国舰船研究设计中心 | 一种确定悬浮式拦截弹最佳布阵间距的方法 |
CN107957319A (zh) * | 2017-11-17 | 2018-04-24 | 湘潭大学 | 均匀荷载面曲率的简支梁裂纹损伤识别方法 |
CN108666374A (zh) * | 2018-05-18 | 2018-10-16 | 通威太阳能(安徽)有限公司 | 一种背面钝化矩阵点式激光开槽导电结构 |
CN110119580A (zh) * | 2019-05-17 | 2019-08-13 | 中国人民解放军军事科学院国防工程研究院 | 一种地面工程目标毁伤评估系统及方法 |
CN110808863A (zh) * | 2019-11-08 | 2020-02-18 | 中国人民解放军陆军工程大学 | 基于区域毁伤的基础设施网络抗毁性评估方法和模型 |
CN111090931A (zh) * | 2019-12-05 | 2020-05-01 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种面目标区域打击规划方法 |
CN112685944A (zh) * | 2021-01-07 | 2021-04-20 | 北京理工大学 | 一种同类多枚杀爆弹对面目标瞄准点规划方法 |
Non-Patent Citations (3)
Title |
---|
周勇: "机载布撒器子弹药抛撒参数设计与散布规律研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 * |
赵国宏: "美军核打击作战任务规划研究", 《指挥与控制学报》 * |
高润芳 等: "目标毁伤效果评估仿真系统研究", 《弹箭与制导学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114419450A (zh) * | 2022-03-29 | 2022-04-29 | 中国人民解放军96901部队 | 一种基于影像特征分析的线性目标毁伤效能快速评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114036783B (zh) | 2022-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11694558B2 (en) | Systems and methods for continuous replanning of vehicle trajectories | |
CN103247041B (zh) | 一种基于局部采样的多几何特征点云数据的分割方法 | |
CN114036783B (zh) | 基于截断机理的线目标功能毁伤判断方法及毁伤计算模型 | |
CN105068550B (zh) | 一种基于拍卖模式的水下机器人多目标选择方法 | |
CN110955953B (zh) | 基于结构化网格的多种类侵爆弹对建筑目标毁伤评估方法 | |
US9460540B2 (en) | Method for animating characters, with collision avoidance based on tracing information | |
CN107063255A (zh) | 一种基于改进果蝇优化算法的三维航路规划方法 | |
CN103528585A (zh) | 一种不等距分割可通行区域的路径规划方法 | |
CN103676944A (zh) | 基于Dubins路径和稀疏A*搜索的无人机航迹规划方法 | |
CN103940420B (zh) | 等航程大圆航线设计方法 | |
CN104317305A (zh) | 一种面向复杂战场威胁的射前航迹确定方法 | |
CN105069245A (zh) | 基于多重搜索技术的三维乘波体快速设计方法 | |
CN110531782A (zh) | 用于社区配送的无人机航迹路径规划方法 | |
CN111737871A (zh) | 一种结合岩土材料特性的落石轨迹三维预测分析方法 | |
CN114139023A (zh) | 一种基于Louvain算法的海上态势生成多目标层次化分群方法 | |
Bernard | The hairpin vortex illusion | |
Cook et al. | Intelligent cooperative control for urban tracking | |
CN112002012B (zh) | 城市区域的可视性分析方法 | |
Benacka | Spreadsheet application showing the proper elevation angle, points of shot and impact of a projectile | |
CN107239559A (zh) | 一种基于矢量网格进行空间动目标位置计算的方法 | |
CN114862152A (zh) | 基于复杂网络的目标重要度评估方法 | |
CN113504793A (zh) | 一种基于Voronoi图的无人机集群路径规划方法 | |
CN106557602B (zh) | 一种基于gis的空间距离校核三维仿真方法 | |
CN118115134A (zh) | 一种基于栅格化的陆战场态势推演时间片优选的方法 | |
CN111815762B (zh) | 一种页岩气压裂裂缝延展过程的三维仿真可视化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |