一种双阵列声学立靶斜入射弹着点定位方法

1.本发明涉及弹着点自动检测声学立靶系统技术，具体涉及一种双阵列声学立靶斜入射弹着点定位方法。


背景技术：

2.在进行直瞄弹道武器射击过程中，能够及时准确的掌握弹丸在靶面上的弹着点位置，不仅对评价武器性能优劣的靶场试验具有重要意义，而且也是检验武器使用和训练效果的重要依据。目前，直瞄弹道武器射击在靶上的弹着点位置主要通过人工目测的方式来判断，这种方式存在时效性差、精确度不够等诸多缺点，不能满足现代化测试和训练的迫切需求。随着电子信息技术发展，弹着点定位方式越来越趋于自动化、智能化，出现了一系列无需实物靶面的弹着点自动检测电子靶系统。声学立靶就是一种非常重要且很常用的电子靶系统。声学立靶系统是基于被动声学定位原理，主要利用弹丸在超声速飞行时产生的弹道激波信号，通过多个麦克风组成阵列对弹道激波信号进行检测和定时、定向等处理，然后根据相应的数学模型计算得到弹着点在虚拟靶面上的位置坐标。声学立靶主要有以下优点：第一，构造简单，不论靶面大小都很容易组建，可以针对不同口径的直瞄弹道武器进行弹着点自动检测与定位；第二，虚拟靶面及其附近的电子设备体积小，不易被弹丸击中；第三，可全天候工作，不受能见度影响；第四，全套系统体积小、重量轻，运输方便，架设容易，机动性强。
3.公开号为cn106247862a的专利公开了一种基于弹道激波的自动声学报靶装置、测量方法及数据处理算法，该方法与文献《弹丸任意角度入射弹着点声学检测模型》均针对小型靶面，采用声学传感器检测超声速弹丸飞行过程中激发的“n”型弹道激波信号计算弹丸飞行弹道和击中靶面的位置。公开号为cn106595394a的专利公开了一种利用音爆(即弹道激波)测量超声速弹丸弹着点的方法，该方法通过采集超声速弹丸飞跃压力传感器时的音爆特征，根据cfd数值模拟方法估计超声速弹丸流场，进而确定弹丸的速度矢量和空间位置。文献《基于激波传播路径的弹头斜入射双三角阵定位模型》提出了一种基于两个三角形阵列的激波信号和到达时间差的弹着点定位方法，仅对于小靶面弹着点的定位较为精确。文献《基于开放式靶标的数学建模及传感器阵列分析》基于相同原理，改进布阵形式，提出了一种基于开放靶标的弹着点定位方法。该方法通用性较强，但是对地形要求较高、布阵复杂。
4.在大靶面声学立靶系统的设计中，近靶超声速弹丸的弹道激波到达部分麦克风的传播距离更长，易受到各种因素的干扰，使系统产生较大的弹着点定位误差。例如，声波传播速度的变化将无法忽略，这会导致弹道激波沿靶面平面向不同麦克风传播的视在速度出现明显差异。如何优化相关数学模型及解算方法来增强弹着点定位性能的抗干扰能力与解算稳定性，是大靶面声学立靶系统的设计难点。此外，为贴近靶场试验和日常训练的应用需求，也要求声学立靶系统应能在任意射击位置和运动射击条件下仍能实现精确的弹着点定位任务，如何保证任意角度斜入射弹丸的弹着点定位精度，是大靶面声学立靶系统的另一
个设计难题。


技术实现要素：

5.本发明的目的在于提供一种双阵列声学立靶斜入射弹着点定位方法。
6.实现本发明目的的技术解决方案为：第一方面。本发明提供一种双阵列声学立靶斜入射弹着点定位方法，步骤如下：
7.步骤1、两个三维立体麦克风阵列分站m1和m2构成双阵列声学立靶系统，每个分站的阵元数l不小于4；两个分站同步采集入射超声速弹丸的弹道激波信号，获得弹道激波到达各麦克风的波达时刻以及相对于各分站的波达方向；
8.步骤2、对步骤1得到的弹道激波波达时刻信息和波达方向信息，解算弹道激波信号分别传往两个分站的弹道脱离点p1和p2的三维点坐标；
9.步骤3、对步骤2得到的弹道脱离点p1和p2的坐标，两点连线构成近靶弹道线段，取近靶弹道线段中点在靶面的投影点作为声学立靶的弹着点坐标估计。
10.第二方面，本发明提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现第一方面所述的方法的步骤。
11.第三方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现第一方面所述的方法的步骤。
12.第四方面，本发明提供一种计算机程序产品，包括计算机程序，该计算机程序被处理器执行时实现第一方面所述的方法的步骤。
13.与现有技术相比，本发明的有益效果为：1)两个l元三维立体麦克风阵列分站同步采集超声速弹丸弹道激波波达时刻信息和波达方向信息，可在任意射击位置和运动射击条件下实现精确的弹着点定位；2)基于弹丸弹道激波脱离点模型，可构建近靶弹道线段，得到线段方程及其单位方向矢量；3)采用对较大测量误差范围低敏感的小孔径阵列分站组网探测技术及其信号处理算法，可有效消除声速变化、风力方向等环境影响，具有更强的环境适应性；4)本发明的方法实现过程便捷，易于实施。
附图说明
14.图1是一种双阵列声学立靶斜入射弹着点定位方法的流程图。
15.图2是分站系统采用的一种4元三维立体麦克风阵列结构示意图。
16.图3是超声速弹丸弹道几何关系示意图。
17.图4是靶面弹着点分布图。
18.图5(a)、图5(b)是弹丸弹着点靶面横向和纵向坐标平均误差等高线图。
19.图6(a)、图6(b)是弹丸弹着点靶面横向和纵向坐标均方根误差等高线图。
具体实施方式
20.本发明提出一种基于双阵列声学立靶斜入射弹着点定位方法，以严格同步的波达时刻信息和高度精确的波达方向信息作为超声速弹丸的近靶弹道激波参数，通过建立近靶
弹道线段的数学模型，解算出准确的弹着点坐标，从而在很大程度上抵消温湿度等环境因素的干扰。
21.结合图1，本发明的一种双阵列声学立靶斜入射弹着点定位方法，步骤如下：
22.步骤1、如图2所示，两个三维立体麦克风阵列分站均以4元正三棱锥阵型为例，水平孔径和俯仰孔径均为0.5m。两个分站同步采集入射超声速弹丸的弹道激波信号，获得弹道激波到达各麦克风的波达时刻以及相对于各分站的波达方向。两个三维立体麦克风阵列分站模型和本发明办法具有普适性，可适用于任意数量麦克风组成的任意形状和尺寸的立体阵列方案。本发明以分站m1为例，4个麦克风m1～m4组成正三棱锥阵型。水平方向的3个阵元m2～m4均匀分布在直径不大于0.5m的水平圆周上，俯仰方向的阵元m1与水平圆心之间的连线与z轴方向平行。阵元m1定义为分站m1的本地参考阵元，即分站m1的本地坐标系原点和时标参考点。定义阵元m1同时也是整个声学立靶系统的全局参考阵元和时标参考点。分站m2的阵列结构与m1完全一致，4个麦克风编号为m5～m8，相互之间的位置关系分别对应分站m1中的m1～m4。其中，中心俯仰方向的阵元m5定义为分站m2的本地坐标系原点和时标参考点。
23.步骤1-1、结合图2和图3，测量弹道激波到达两个分站内各麦克风的波达时刻为tj,j＝1～8。其中，麦克风m1～m4属于分站m1，m1是其本地参考阵元；麦克风m5～m8属于分站m2，m5是其本地参考阵元。m1同时也是整个双站立靶系统的全局参考阵元。计算得到两个分站m1和m2内各麦克风与系统全局参考阵元之间的波达时间差(tdoa)为δt
j1
＝t
j-t1,j＝1～8；
24.步骤1-2、根据步骤1-1得到的tdoa信息，分别解算弹道激波相对于分站m1和m2的波达方向，提取对应的单位方向矢量和
25.步骤1-2-1、以分站m1为例，根据三维立体麦克风阵列分站结构相关预设常数和tdoa方程，估计弹道激波相对于分站m1的波达方向：
[0026][0027]
式中常数a1～a7为与三维立体麦克风阵列结构相关的预设常数，δt
21
、δt
31
、δt
41
为正三棱锥分布的阵列分站m1内麦克风m2、m3、m4相对中间的参考麦克风m1的tdoa测量值；
[0028]
步骤1-2-2、将步骤1-2-1得到相对于分站m1的波达方向矢量a1β(x
s1
,y
s1
,z
s1
)
t
归一化处理，得到波达方向单位矢量
[0029]
步骤2、对步骤1得到的弹道激波波达时刻信息和波达方向信息，解算弹道激波信号分别传往两个分站的弹道脱离点p1和p2的三维点坐标；
[0030]
步骤2-1、定义两个分站m1和m2的位置坐标矢量分别为和两个弹道脱离点p1和p2到各自接收分站m1和m2之间的传播距离分别为和根据步骤1-2得到的弹道激波相对于两个分站的单位方向矢量和获得p1和p2两点之间的连线即近靶弹道线段矢量：
[0031][0032]
步骤2-2、由于两个分站相距仅有数米，近靶弹道线段一般很短，因而可认为超声速弹丸在p1和p2两点之间匀速直线飞行，p1和p2两点处的弹道激波锥角即马赫角θ1和θ2近似相等。定义θ1＝θ2＝θ，得到关系式：
[0033][0034]
式中ma和v分别表示近靶弹道线段的超声速弹丸马赫数和飞行速度，c为大气声速，《x,y》表示矢量x和y的内积。将代入步骤2-1中的近靶弹道线段表达式，得到两个弹道脱离点p1和p2到达各自接收分站m1和m2的传播距离之和满足：
[0035][0036]
步骤2-3、由于两个分站m1和m2的本地参考阵元分别为1和5，弹道激波信号在两个本地参考阵元之间的tdoa满足以下关系：
[0037][0038]
式中c1和c2分别为弹道激波脱离弹道后传向两个分站的大气声速。在无风理想条件下，c1βc2＝c。进而得到两个弹道脱离点p1和p2到各自接收分站m1和m2的传播距离之差满足：
[0039][0040]
步骤2-4、结合步骤2-2和步骤2-3求解得到和进而得到两个弹道激波脱离点p1和p2的位置坐标矢量：
[0041][0042]
步骤3、对步骤2得到的弹道脱离点p1和p2坐标，两点连线构成近靶弹道线段，取近靶弹道线段中点在靶面的投影点作为声学立靶的弹着点坐标估计，具体包括以下步骤：
[0043]
步骤3-1、定义步骤2-4得到的两个弹道脱离点坐标分别为p1(x
p1
,y
p1
,z
p1
)
t
和p2(x
p2
,y
p2
,z
p2
)
t
，建立p1和p2两点连线得到的斜入射超声速弹丸近靶弹道线段，线段中点坐标为p3(x
p3
,y
p3
,z
p3
)
t
：
[0044][0045]
步骤3-2、定义靶平面为空间平面s，靶平面方程为：x＝0。据步骤3-1可知近靶弹道线段中点p3(x
p3
,y
p3
,z
p3
)
t
不在靶平面上，将其投影到靶平面的点坐标定义为ts(x
t
,y
t
,z
t
)
t
：
[0046][0047]
即为声学立靶的弹着点坐标估计。
[0048]
结合附图，下面给出具体实施例：
[0049]
本实施例提供了一种基于两个三维立体麦克风阵列分站和超声速弹丸弹道激波信号的斜入射弹着点定位方法，具体实施通过图1所示步骤实现，参数设置和效果通过以下实验进一步说明。
[0050]
参照图3，对基于双阵列声学立靶系统的超声速弹丸定位实验参数设置如下：采用右手系三维直角坐标系，坐标系原点距离地面高度为5m，方位角定义为待测矢量在oxy平面上的投影线与x轴之间的夹角，俯视顺时针方向为正；俯仰角定义为待测矢量与其在oxy平面上的投影线之间的夹角，向上方向为正。平面oyz(即x＝0)为默认靶面，靶心t0(0,0,0)为坐标系原点，靶面区域为以靶心为中心的10m
×
10m区域。分站m1和m2架设高度均为1.5m，分站m1的中心阵元坐标作为分站m1的坐标m1(0,-5,3.5)，分站m2的中心阵元坐标作为分站m2的坐标m2(0,-5,-3.5)，分站坐标单位均为m。定义各麦克风三维点坐标的自定位误差均服从ε～n(0,0.1)的高斯分布，单位为m；定义各麦克风波达时刻的测量误差均服从ε
t
～n(0,2.5)的高斯分布，单位为ms；定义各分站激波方位角测量误差和俯仰角测量误差分别服从ε
α
～n(0,0.25)和的高斯分布，单位为度；定义超声速弹丸的近靶飞行速度为v＝800m/s；定义大气声波传播速度c＝342m/s。
[0051]
参照图4，在10m
×
10m的靶面上以0.5m为间距等距取点，作为真实弹着点位置。对靶面共计441个预设弹着点分别单独进行100次蒙特卡洛试验，得到弹着点坐标估计值与真实值之间的误差。图5(a)和5(b)分别给出靶面弹着点在y维度和z维度的定位平均误差分布，可以看出针对10m
×
10m的大靶面，两个维度的定位平均误差均能达到5cm以下，靶心处的定位平均误差在两个维度均可以达到3cm。靶面弹着点在y维度和z维度的定位均方根误差分布图分别如图6(a)和6(b)所示，可以看出针对10m
×
10m的大靶面，两个维度的定位均方根误差均能达到6cm以下，靶心处的定位均方根误差在两个维度均可以达到3cm。
[0052]
本发明借助弹道模型和解析几何理论，同步采集入射超声速弹丸的弹道激波，根
据弹道激波到达各麦克风的波达时刻以及相对于各个分站的波达方向信息，解算分别对应两个分站的弹道激波脱离点的三维点坐标，再通过两个弹道脱离点连线构成近靶弹道线段，进而得到近靶弹道线段的直线方程及其单位方向矢量，最后将其与靶面方程的交点坐标作为声学立靶的弹着点坐标估计。本发明能够在无需实物靶面的前提下，针对大靶面斜入射超声速弹丸的弹着点进行被动式的高精度自动定位，同时具备机动性能好、系统结构简单、架收便利、不易损坏、可全天候工作、成本低廉等优点，对于满足靶场试验和日常训练的现代化技术保障需求具有重要意义。