一种脉冲声信号识别分类方法
阅读：745发布：2020-05-13
专利汇可以提供一种脉冲声信号识别分类方法专利检索，专利查询，专利分析的服务。并且本 发明 火炮脉冲声 信号 识别技术领域，具体涉及一种脉冲声信号识别分类方法，包括：对检测到的脉冲声信号分别进行低通和高通滤波，通过高频和低频 能量 对比初步判定脉冲声信号的类型；对低频能量高或高频能量高的脉冲声信号进行火炮声信号判定；对判定为火炮声信号的脉冲声信号进行参数化 功率谱 建模，将模型参数作为 特征向量 分别输入与初步判定的脉冲声信号类型相对应的多元高斯混合模型中，最终判定脉冲声信号的类型。本发明能够对膛口波、弹道波、爆炸波三种脉冲波进行识别分类，提升了声探测装备的识别率，降低了探测虚警率，还可解决武器口径识别分类的问题。，下面是一种脉冲声信号识别分类方法专利的具体信息内容。
1.一种脉冲声信号识别分类方法，其特征在于，包括：
对检测到的脉冲声信号分别进行低通和高通滤波，通过高频和低频能量对比初步判定脉冲声信号的类型；
对低频能量高或高频能量高的脉冲声信号进行火炮声信号判定；
对判定为火炮声信号的脉冲声信号利用AR模型进行参数化功率谱建模，将AR模型参数a1，…，ap作为特征向量分别输入与初步判定的脉冲声信号类型相对应的多元高斯混合模型中，最终判定脉冲声信号的类型。
2.根据权利要求1所述的一种脉冲声信号识别分类方法，其特征在于，所述方法具体包括：
将低频能量高的脉冲声信号初步判定为膛口波或是爆炸波信号；
基于非参数化功率谱特征提取及BP神经网络二分类对低频能量高的脉冲声信号进行火炮声信号判定；
对判定为火炮信号的脉冲声信号，利用AR模型进行参数化功率谱建模，将AR模型参数a1，…，ap作为特征向量分别输入膛口波模型及爆炸波模型，若膛口波模型输出结果大时，则判定为膛口波，若爆炸波模型输出结果大时，则判定为爆炸波。
3.根据权利要求1所述的一种脉冲声信号识别分类方法，其特征在于，所述方法具体包括：
将高频能量高的脉冲声信号初步判定为弹道波信号；
基于非参数化功率谱特征提取及BP神经网络二分类对高频能量高的脉冲声信号进行火炮声信号判定；
对判定为火炮信号的弹道波信号，利用AR模型进行参数化功率谱建模，将AR模型参数a1，…，ap作为特征向量分别输入常规火炮弹道波模型及火箭炮弹道波模型，若常规火炮弹道波模型输出结果大时，则判定为常规火炮弹道波，若火箭炮弹道波模型输出结果大时，则判定为火箭炮弹道波。
4.根据权利要求3所述的一种脉冲声信号识别分类方法，其特征在于，还包括：
对判定为火炮信号的弹道波信号，进行形状特征提取，识别弹丸口径。
5.根据权利要求2或4所述的一种脉冲声信号识别分类方法，其特征在于，采用低通和高通滤波器对检测到的脉冲声信号分别进行低通和高通滤波，所述低通和高通滤波器为椭圆滤波器，所述低通和高通滤波器的分界频率点为60Hz。
6.根据权利要求5所述的一种脉冲声信号识别分类方法，其特征在于，基于周期图法进行非参数化功率谱特征提取。
7.根据权利要求6所述的一种脉冲声信号识别分类方法，其特征在于，所述BP神经网络二分类包括：
将炮声信号作为正样本，非炮声信号作为负样本，将正、负样本特征向量作为BP神经网络的输入，进行二分类识别。
8.根据权利要求4所述的一种脉冲声信号识别分类方法，其特征在于，还包括：
根据不同的弹丸尺寸参数和测点距弹道线垂直距离建立弹丸口径识别模板。
一种脉冲声信号识别分类方法
技术领域
[0001] 本发明火炮脉冲声信号识别技术领域，具体涉及一种脉冲声信号识别分类方法。
背景技术
[0002] 火炮在发射时火药爆炸会产生强烈的冲击波，在空气中以脉冲声波的形式向外传播；弹丸在空气中超音速飞行时与空气强烈摩擦会产生强烈的弹道激波，也以脉冲声波的形式向外传播；弹丸爆炸也是产生脉冲声。另外自然界中脉冲声源多、难以区分。对于声探测装备来说，如果不能有效的识别目标脉冲声信号，将提升装备的探测虚警率，将会使火控系统打击错误目标。
[0003] 火炮产生的脉冲声的识别分类主要有这几个方面的内容：对发射时由于火药爆炸产生的热空气向外膨胀形成的冲击波进行识别，这种波称为膛口波；对弹丸超音速飞行过程中与空气剧烈摩擦产生的弹道激波进行识别，这种波称为弹道波；对弹丸爆炸产生的声波进行识别，这种波称为爆炸波。除了这几个方面，通过对不同口径火炮弹丸激波进行分析，达到对武器口径、种类进行识别也是一项重要内容。脉冲声识别分类属于模式识别的范畴，一般包括特征选择和提取、分类器设计等几个步骤。
[0004] 目前，对于脉冲声识别分类的方法主要有：求取脉冲声信号的时域自相关系数，在具有区分性的延迟时间段计算自相关系数的平均值，经过大量分析求取每种脉冲声信号自相关系数平均值的区间，不同的脉冲声信号的平均自相关系数的值处于不同的区间，从而达到识别分类的目的，这种方法必须要求区分区间不重叠才能达到很好的分类效果，在实际当中为了保证识别区间不重叠，区间范围会很小，这将会降低识别率。
[0005] 在时域对脉冲声进行识别的另一种方法是提取波形的形状特征，各种形状特征组合形成特征矢量空间，通过大量的波形分析，提取足够多的特征矢量形成特征矢量集合，采取支持向量机或是神经网络的方法对特征矢量集合进行离线训练，训练完毕后应用到实际当中对脉冲声进行识别分类，这种方法提取的多种形状特征中有的可区分性并不好，这将会影响离线训练的结果，使得应用在实际当中的识别率很低。
[0006] 在变换域的脉冲声识别分类方法主要是在频域，对时域脉冲信号进行频谱变换，在频域提取目标的能量、谱重心、极大值频率等特征，将这些特征组合形成特征向量，经过大量的数据分析获得特征矢量集合，采用神经网络的方法对特征矢量集合进行离线训练，然后应用的实际当中，这种方法和时域特征向量法类似，同样存在选择特征可分性不好，影响实际识别率的问题。
发明内容
[0007] 本发明的目的在于提供一种脉冲声信号识别分类方法，以完成对火炮脉冲声的有效识别分类，提升声探测装备的识别率，降低探测虚警率。
[0008] 本发明提供了一种脉冲声信号识别分类方法，包括：
[0009] 对检测到的脉冲声信号分别进行低通和高通滤波，通过高频和低频能量对比初步判定脉冲声信号的类型；
[0010] 对低频能量高或高频能量高的脉冲声信号进行火炮声信号判定；
[0011] 对判定为火炮声信号的脉冲声信号利用AR模型进行参数化功率谱建模，将AR模型参数a1，…，ap作为特征向量分别输入与初步判定的脉冲声信号类型相对应的多元高斯混合模型中，最终判定脉冲声信号的类型。
[0012] 进一步地，该方法具体包括：
[0013] 将低频能量高的脉冲声信号初步判定为膛口波或是爆炸波信号；
[0014] 基于非参数化功率谱特征提取及BP神经网络二分类对低频能量高的脉冲声信号进行火炮声信号判定；
[0015] 对判定为火炮信号的脉冲声信号，利用AR模型进行参数化功率谱建模，将AR模型参数a1，…，ap作为特征向量分别输入膛口波模型及爆炸波模型，若膛口波模型输出结果大时，则判定为膛口波，若爆炸波模型输出结果大时，则判定为爆炸波。
[0016] 进一步地，该方法具体包括：
[0017] 将高频能量高的脉冲声信号初步判定为弹道波信号；
[0018] 基于非参数化功率谱特征提取及BP神经网络二分类对高频能量高的脉冲声信号进行火炮声信号判定；
[0019] 对判定为火炮信号的弹道波信号，利用AR模型进行参数化功率谱建模，将AR模型参数a1，…，ap作为特征向量分别输入常规火炮弹道波模型及火箭炮弹道波模型，若常规火炮弹道波模型输出结果大时，则判定为常规火炮弹道波，若火箭炮弹道波模型输出结果大时，则判定为火箭炮弹道波。
[0020] 进一步地，该方法还包括：
[0021] 对判定为火炮信号的弹道波信号，进行形状特征提取，识别弹丸口径。
[0023] 进一步地，基于周期图法进行非参数化功率谱特征提取。
[0024] 进一步地，BP神经网络二分类包括：
[0025] 将炮声信号作为正样本，非炮声信号作为负样本，将正、负样本特征向量作为BP神经网络的输入，进行二分类识别。
[0026] 进一步地，该方法还包括：
[0027] 根据不同的弹丸尺寸参数和测点距弹道线垂直距离建立弹丸口径识别模板。
[0029] 图1是本发明一种脉冲声信号识别分类方法的流程图；
[0030] 图2是本发明低通滤波器归一化幅频响应特性示意图；
[0031] 图3是本发明高通滤波器归一化幅频响应特性示意图；
[0032] 图4是本发明GMM似然度计算过程示意图；
[0033] 图5是本发明弹道波状态图；
[0034] 图6是本发明口径识别模板示意图。
具体实施方式
[0035] 下面结合附图所示的各实施方式对本发明进行详细说明，但应当说明的是，这些实施方式并非对本发明的限制，本领域普通技术人员根据这些实施方式所作的功能、方法、或者结构上的等效变换或替代，均属于本发明的保护范围之内。
[0036] 本实施例提供了一种脉冲声信号识别分类方法，包括：
[0037] 对检测到的脉冲声信号分别进行低通和高通滤波，通过高频和低频能量对比初步判定脉冲声信号的类型；
[0038] 对低频能量高或高频能量高的脉冲声信号进行火炮声信号判定；
[0039] 对判定为火炮声信号的脉冲声信号进行参数化功率谱建模，将模型参数作为特征向量分别输入与初步判定的脉冲声信号类型相对应的多元高斯混合模型中，最终判定脉冲声信号的类型。
[0040] 通过该脉冲声信号识别分类方法，能够对膛口波、弹道波、爆炸波三种脉冲波进行识别分类，提升了声探测装备的识别率，降低了探测虚警率，还可解决武器口径识别分类的问题。
[0041] 下面对本发明进行详细说明。
[0042] 1、脉冲声识别分类总体方案。
[0043] 脉冲声信号的检测和提取过程不在本实施例中叙述，当已经检测到脉冲声信号后，对提取的脉冲声波形分别进行低通和高通滤波，一般火炮脉冲声信号中的膛口波和爆炸波主要能量集中在低频段，而弹道波的主要能量响应的集中在高频段，通过对比高频和低频能量可以初步判定该脉冲信号可能是那种类型。
[0044] 当判定该脉冲信号低频信号强时，实际上还不能确定是不是火炮膛口波或是爆炸波信号，对其进行非参数化功率谱特征提取，并作为BP神经网络的输入，进行二分类识别，根据识别结果判定是否是火炮声信号；当判定是火炮声信号时，对其进行参数化功率谱建模，将模型参数作为特征向量分别输入多元高斯混合模型GMM1和GMM2，GMM1代表膛口波模型，GMM2代表爆炸波模型，GMM1输出结果大时表示该段信号为膛口波，反之为爆炸波。
[0045] 当判定该脉冲声信号为高频信号强时，同样经过非参数化功率谱特征提取、BP神经网络识别两步骤判定该段信号是否是弹道波；当判定是弹道波时一方面提取形状特征进行口径识别，另一方面对其进行参数化功率谱建模，将模型参数作为特征向量分别输入高斯混合模型GMM3和GMM4，GMM3代表常规火炮弹道波模型，GMM4代表火箭炮弹道波模型，GMM3输出结果大时表示该段信号为常规火炮弹道波，反之为火箭炮弹道波。整体识别流程如图1所示。
[0046] 2、滤波器设计
[0047] 经过大量火炮声信号的分析，一般弹道波在60Hz～200Hz还有很强的能量，而爆炸波和膛口波主要能量集中在5Hz～60Hz，所以用60Hz作为分界频率点设计滤波器。由于椭圆滤波器在相同阶数下能获得更窄的过度带宽和较小的阻带波动，所以本实施例采用该种滤波器分别设计高通和低通滤波器，其归一化幅频响应特性如图2和图3所示。
[0048] 3、非参数化功率谱特征分析
[0049] 采用周期图法进行非参数化功率谱计算，公式如下：
[0050]
[0051]
[0052] 式中，x(n)为时域采样序列，PX(w)为求得的功率谱，c(n)为窗函数。周期图方法估计的功率谱为有偏估计，为了减少偏差，窗函数c(n)采用的是Hamming窗。
[0053] 对信号做1倍采样率长度的率谱，当判断为低频信号时取5Hz～60Hz的频段范围共56个频点幅值作为低频特征向量FBf＝(f1,f2,…,f56)。当判断为高频信号时取60Hz～200Hz频段，每10个频点对应的幅值做平均，重叠率0.5，共可得到38个数据，将其作为高频特征向量FHf＝(f1,f2,…,f38)。
[0054] 4、BP神经网络分类器设计
[0055] 使用神经网络进行识别分类时首先必须采用大量的样本对其进行充分的训练，使其能适应多种状态下的目标特性。本实施例中BP神经网络分类器的作用是对炮声信号和非炮声信号进行二分类，将炮声信号作为正样本，非炮声信号作为负样本，负样本建立方式与正样本建立方式一致。
[0056] 将正负样本特征向量作为神经网络的输入，神经网络输入层一般取常值为1的变换函数为神经元，维数与特征向量的维数一致。由于是二分类，所以输出层维数为2，取比例函数为神经元。隐含层神经元取双曲函数，其数量会影响到最终的识别性能，经过大量的分析测试，本实施例中两个神经网络隐含层个数都取25个，最终都能达到较好的性能。采用160个膛口波和爆炸波组成的低频炮声信号组成的正样本，负样本噪声包括车辆、飞机、风等共500个，弹道波信号组成的高频正样本包括100mm、122mm、130mm、155mm等口径的弹丸飞行激波共100个，高频神经网络和低频神经网络共用负样本，在对每个神经网络进行充分的训练后，采用一定数量的测试样本对其进行测试，测试结果如表1所示。
[0057] 表1 BP神经网络识别性能测试
[0058]信号类型 个数 识别正确
膛口波\爆炸波 30 27
弹道波(100迫、122、130、155) 30 28
其它噪声包括(车辆、飞机、风等) 30 26
膛口波\爆炸波 30 27
弹道波(100迫、122、130、155) 30 28
其它噪声包括(车辆、飞机、风等) 30 26
[0059] 5、参数化功率谱分析
[0060] 参数化功率谱估计主要是利用AR模型(Auto Regressive Mode)对样本序列进行建模，然后通过谱估计子计算功率谱密度。根据Wold分解定理，任何一个具有有限方差的ARMA或MA过程都可以表示成惟一的、阶数有可能无穷大的AR过程。相比ARMA模型不仅需要AR和MA阶数的确定，而且还需要AR和MA参数估计(其中MA参数的估计还需要求解非线性方程)，AR模型只涉及定阶和参数估计，因此可用AR模型来对信号进行建模。相应的，功率谱密度为：
[0061]
[0062] 这样，只要求得激励白噪声的方差σ2和AR参数a1,…,ap，就可求得Px(ω)。设采样序列为s(n)n＝0,1,…N-1，AR参数的估计可使用Yule-Walker方程，并通过Levinson递推求解该方程，步骤如下：
[0063]
[0064] 一阶AR模型的Yule-Walker矩阵方程为：
[0065]
[0066] 可解a11与 分别为
[0067] a11＝-Rx(1)/Rx(0)
[0068]
[0069] 二阶AR模型的矩阵方程：
[0070]
[0071] 解得a22，a21，σ2分别为
[0072]
[0073]
[0074]
[0075] 以此类推得递推公式：
[0076]
[0077] aki＝ak-1,i+akkak-1,k-i
[0078]
[0079] AR参数a1,…,ap可以完全代表信号的功率谱特征，阶数p越大表征的功率谱特征越详细，单会带来很多冗余的信息，越小表征的功率谱特征越粗略，但会损失一定的细节信息，适当的p值是必要的，经过多次尝试，取p＝48可以达到较好的功率谱特征表示，同时也能达到较好的分类识别效果。
[0080] 6、多元高斯混合模型建模
[0081] 多元高斯混合模型(GMM)是状态数为1的连续型隐马尔科夫模型，该状态对应的观察序列服从多元高斯概率密度函数。假设一个具有M个高斯分量D元高斯混合模型λ中，观察序列为a，则该序列在λ下似然度表示为：
[0082]
[0083] 其中第i个高斯分量定位为：
[0084]
[0085] 式中，M混合的多元高斯模型数量，αi为混合权重，满足 bi(a)为多元高斯概率密度函数，μi为多元高斯概率密度函数的均值向量，Σi为协方差矩阵。因此GMM模型可由均值向量、协方差矩阵和混合权重表示，另λ＝{αi,μi,Σi}，则GMM模型计算过程可表示如图4所示。
[0086] 本实施例建立4个GMM模型，GMM1是针对膛口波建立的模型，GMM2是针对爆炸波建立的模型，GMM3是针对常规火炮弹道波建立的模型，GMM4是针对火箭炮弹道波建立的模型。将48阶AR参数a＝[a1,…,a48]作为观察序列，得到D＝48，混合的多元高斯模型数量取M＝
10，所有模型λjj＝1,2,3,4的值通过模型训练确定，设每种模型取100个观察序列集合则每种模型的最大似然概率为：
10，所有模型λjj＝1,2,3,4的值通过模型训练确定，设每种模型取100个观察序列集合则每种模型的最大似然概率为：
[0087]
[0088] 模型训练的过程就是找到参数λj，使得似然概率最大：常采用期望最大化算法(Expectation Maximization，EM)算法来对其进行估计。当模型参数λj都训练完毕后，在实际使用时输入观察序列a＝[a1,…,a48]，看哪种模型的似然度最大，那么该脉冲信号就属于哪种类型。
[0089] 7、口径识别
[0090] 本实施例中采用弹道波对弹丸的口径进行识别分类，弹道波是弹丸在空气中超音速飞行过程中，与空气剧烈碰撞类似与撕裂空气而产生的声波，其形状类似字母“N”，一般具备连续的6个状态，如图5所示，分别为①开始升压状态，②负压最大状态，③正压最大状态，④恢复状态。
[0091] 取①和④之间的距离为弹道波持续时间T，②和③的绝对值的平均值为声压值P，则T和P理论上分别满足如下公式：
[0092]
[0093]
[0094] 式中：
[0095] K—弹形系数；
[0096] Ma—马赫数；
[0097] d—直径(按翼展计算)；
[0098] c—声速；
[0099] r—测点距弹道垂直距离；
[0100] l—长度；
[0101] p,p0－激波声压和大气环境压力。
[0102] 根据不同的弹丸尺寸参数和测点距弹道线垂直距离可建立弹丸口径识别模板，如图6所示，经过实际数据的识别测试，识别正确率85％以上。
[0103] 本实施例提出的脉冲声信号识别分类方法，可完成对火炮脉冲声中的膛口波、弹道波、爆炸波进行有效识别分类，并可判断弹丸的口径，提升了声探测装备的识别率，降低了探测虚警率。
[0104] 上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明，它们并非用以限制本发明的保护范围，凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。
[0105] 对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 一种航天器热控产品的三维快速建模方法 | 2020-05-08 | 589 |
| 一种有效的图像镜面反射移除方法 | 2020-05-08 | 893 |
| 一种参数化功能梯度立方点阵结构设计方法及材料 | 2020-05-14 | 792 |
| 一种直升机旋翼运动分析系统 | 2020-05-14 | 28 |
| 一种高速连接器的PCB及其优化方法和系统 | 2020-05-11 | 543 |
| 一种基于控制参数化的无人机路径规划方法 | 2020-05-12 | 200 |
| 一种乱序高性能核的混合分支预测装置及方法 | 2020-05-08 | 230 |
| 基于GNSS和GEO卫星的电离层投影函数建模方法 | 2020-05-13 | 279 |
| 用于通过高斯过程进行机器学习的方法 | 2020-05-08 | 445 |
| 一种基于迭代法构建相关滤波响应的目标跟踪方法 | 2020-05-08 | 317 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
分析报告
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
QQ群二维码
意见反馈
意见反馈