直升机声信号增强方法及装置与流程

本发明涉及一种信号增强方法，尤其涉及一种基于谱稀疏约束非负矩阵分解的直升机声信号增强方法及装置。

背景技术：

直升机是战场上一种威胁性很大的武器，有“空中坦克”之称；它能够突然出现并超低空飞行，以躲避传统雷达的搜索。由于直升机在飞行过程中辐射高强度且有明显特征的噪声，因而研究者尝试采用被动声探测的方法，完成直升机目标的检测、定位和识别等一系列任务。

复杂背景噪声干扰是制约上述任务处理效果的主要因素，因而需要探索有效的噪声抑制和信号增强方法。噪声抑制算法已有近几十年的发展，常见方法包括谱减法、维纳滤波、自适应滤波噪声抵消方法和波束形成算法等。谱减法和维纳滤波均需对一定时间段的噪声谱进行估计，其基本假设是噪声在该时间段具有统计平稳性；然而，时间环境噪声是短时非平稳的，因而制约了该类算法的性能。自适应滤波算法需要选取合适的参考噪声通道，多通道差分信号是较好的选择，但在目标方位未知时无法获取；波束形成信号增强方法同样需要事先知道目标方位，同时其性能也受到阵元个数和阵列孔径的限制。

有鉴于此，需要一种能够使声探测设备在复杂背景噪声中提取质量较高的目标直升机声信号，降低环境噪声影响的方案。

技术实现要素：

本发明的目的在于提供一种基于谱稀疏约束非负矩阵分解的直升机声信号增强方法及装置，使声探测设备在复杂背景噪声中提取质量较高的目标直升机声信号，降低环境噪声的影响，为增加目标探测距离、方位参数精确估计和目标准确识别奠定基础。

本发明提供了一种直升机声信号增强方法，包括：

步骤一：采集典型应用环境下的背景噪声信号作为噪声训练数据，对噪声时频幅度谱进行离线字典学习，得到用信号字典矩阵和噪声字典矩阵组合表示的带噪信号时频幅度谱；

步骤二：采集带噪直升机信号，固定所述噪声字典矩阵，通过所述信号字典矩阵对带噪直升机信号时频幅度谱进行在线字典学习，并对学习到的直升机信号字典施加谱稀疏约束，通过估计得到直升机信号时频幅度谱；

步骤三：结合带噪直升机信号的相位及估计得到的直升机信号时频幅度谱，通过逆变换得到增强的目标直升机信号。

进一步地，在所述步骤一中，采用非负矩阵分解方法对所述噪声时频幅度谱进行字典学习。

进一步地，所述步骤一包括：

在典型应用环境录制一定时长的背景噪声波形，并采用按间隔抽点的方式将噪声信号降采样，以与直升机声信号的采样频率相匹配；

对采集到的噪声信号分帧、加窗，对每个短时窗信号进行傅里叶变换，获得噪声时频幅度谱。

进一步地，所述步骤二包括：

对采集到的带噪直升机声信号分帧、加窗，对每个短时窗信号进行傅里叶变换，获得带噪直升机信号时频幅度谱。

本发明还提供了一种直升机声信号增强装置，包括：

离线字典学习模块，用于采集典型应用环境下的背景噪声信号作为噪声训练数据，对噪声时频幅度谱进行离线字典学习，得到用信号字典矩阵和噪声字典矩阵组合表示的带噪信号时频幅度谱；

信号增强模块，用于采集带噪直升机信号，固定所述噪声字典矩阵，通过所述信号字典矩阵对带噪直升机信号时频幅度谱进行在线字典学习，并对学习到的直升机信号字典施加谱稀疏约束，通过估计得到直升机信号时频幅度谱；

目标信号获取模块，用于结合带噪直升机信号的相位及估计得到的直升机信号时频幅度谱，通过逆变换得到增强的目标直升机信号。

进一步地，所述离线字典学习模块采用非负矩阵分解方法对所述噪声时频幅度谱进行字典学习。

进一步地，所述离线字典学习模块在典型应用环境录制一定时长的背景噪声波形，并采用按间隔抽点的方式将噪声信号降采样，以与直升机声信号的采样频率相匹配；

对采集到的噪声信号分帧、加窗，对每个短时窗信号进行傅里叶变换，获得噪声时频幅度谱。

进一步地，所述信号增强模块对采集到的带噪直升机声信号分帧、加窗，对每个短时窗信号进行傅里叶变换，获得带噪直升机信号时频幅度谱。

与现有技术相比本发明的有益效果是：

能够从复杂背景噪声中抽取出质量较高的目标信号，可作为目标检测、定向与识别的前端处理，通过对实际数据的处理，对复杂背景噪声的抑制能力可达到15db以上，从而大大提升了直升机目标声信号的信噪比，进一步可提高系统目标探测距离、方位估计精度和目标识别准确率。

附图说明

图1是本发明一种直升机声信号增强方法的流程图；

图2是本发明一种直升机声信号增强装置的结构框图。

具体实施方式

下面结合附图所示的各实施方式对本发明进行详细说明，但应当说明的是，这些实施方式并非对本发明的限制，本领域普通技术人员根据这些实施方式所作的功能、方法、或者结构上的等效变换或替代，均属于本发明的保护范围之内。

参图1所示，本实施例提供了一种基于谱稀疏约束非负矩阵分解的直升机声信号增强方法，包括：

步骤s1：采集典型应用环境下的背景噪声信号作为噪声训练数据，对噪声时频幅度谱进行离线字典学习，得到用信号字典矩阵和噪声字典矩阵组合表示的带噪信号时频幅度谱。

步骤s2：采集带噪直升机信号，固定所述噪声字典矩阵，通过所述信号字典矩阵对带噪直升机信号时频幅度谱进行在线字典学习，并对学习到的直升机信号字典施加谱稀疏约束，通过估计得到直升机信号时频幅度谱。

步骤s3：结合带噪直升机信号的相位及估计得到的直升机信号时频幅度谱，通过逆变换得到增强的目标直升机信号。

通过该基于谱稀疏约束非负矩阵分解的直升机声信号增强方法，能够从复杂背景噪声中抽取出质量较高的目标信号，可作为目标检测、定向与识别的前端处理，通过对实际数据的处理，对复杂背景噪声的抑制能力可达到15db以上，从而大大提升了直升机目标声信号的信噪比，进一步可提高系统目标探测距离、方位估计精度和目标识别准确率。

在本实施例中，在步骤s1中，采用非负矩阵分解方法(nonnegativematrixfactorization，nmf)对所述噪声时频幅度谱进行字典学习。

在本实施例中，步骤s1包括：

在典型应用环境录制一定时长的背景噪声波形，并采用按间隔抽点的方式将噪声信号降采样，以与直升机声信号的采样频率相匹配；

对采集到的噪声信号分帧、加窗，对每个短时窗信号进行傅里叶变换，获得噪声时频幅度谱。

在本实施例中，步骤s2包括：

对采集到的带噪直升机声信号分帧、加窗，对每个短时窗信号进行傅里叶变换，获得带噪直升机信号时频幅度谱。

本实施例还提供了一种基于谱稀疏约束非负矩阵分解的直升机声信号增强装置，包括：

离线字典学习模块10，用于采集典型应用环境下的背景噪声信号作为噪声训练数据，对噪声时频幅度谱进行离线字典学习，得到用信号字典矩阵和噪声字典矩阵组合表示的带噪信号时频幅度谱；

信号增强模块20，用于采集带噪直升机信号，固定所述噪声字典矩阵，通过所述信号字典矩阵对带噪直升机信号时频幅度谱进行在线字典学习，并对学习到的直升机信号字典施加谱稀疏约束，通过估计得到直升机信号时频幅度谱；

目标信号获取模块30，用于结合带噪直升机信号的相位及估计得到的直升机信号时频幅度谱，通过逆变换得到增强的目标直升机信号。

通过该基于谱稀疏约束非负矩阵分解的直升机声信号增强装置，能够从复杂背景噪声中抽取出质量较高的目标信号，可作为目标检测、定向与识别的前端处理，通过对实际数据的处理，对复杂背景噪声的抑制能力可达到15db以上，从而大大提升了直升机目标声信号的信噪比，进一步可提高系统目标探测距离、方位估计精度和目标识别准确率。

在本实施例中，离线字典学习模块10采用非负矩阵分解方法(nonnegativematrixfactorization，nmf)对所述噪声时频幅度谱进行字典学习。

在本实施例中，离线字典学习模块10在典型应用环境录制一定时长的背景噪声波形，并采用按间隔抽点的方式将噪声信号降采样，以与直升机声信号的采样频率相匹配；

对采集到的噪声信号分帧、加窗，对每个短时窗信号进行傅里叶变换，获得噪声时频幅度谱。

在本实施例中，信号增强模块20对采集到的带噪直升机声信号分帧、加窗，对每个短时窗信号进行傅里叶变换，获得带噪直升机信号时频幅度谱。

下面对本发明作进一步详细说明。

本发明利用nmf方法在时频域分别对噪声和信号的特征进行字典学习，带噪信号的时频幅度谱可由信号字典和噪声字典的组合表示，保留信号部分进行重构可获得增强的目标信号。首先，在典型应用环境采集一定时长的纯背景噪声，利用普通nmf方法对噪声时频幅度谱进行离线字典学习；然后，固定噪声字典矩阵，信号字典矩阵则对带噪信号时频幅度谱在线学习，并对学习到的字典施加谱稀疏约束以突出直升机目标信号中的离散谐波部分(直升机信号谱为离散谐波谱和连续宽带谱的检测，其中前者是目标检测、定向与识别所利用的信号分量，因此也是信号增强所关心的部分。因此，可在字典学习过程中施加谱稀疏约束，从而可增强直升机信号中的谐波分量)，该方法可称为谱稀疏约束非负矩阵分解(spectralsparsenessconstrainednonnegativematrixfactorization，简记为sscnmf)；最后，保留信号部分的时频幅度谱，并结合带噪信号的相位，通过短时傅里叶逆变换重建目标声信号。具体流程如下：

1、噪声字典学习

1)在典型应用环境利用声音采集设备录制一定时长(如20min)的背景噪声波形，并采用按一定间隔抽点的方式将噪声信号降采样为1024hz，以与直升机声信号的采样频率相匹配。

2)对采集到的噪声信号n(t)分帧、加hamming窗，窗长取256点(对应0.25s)，帧间重叠50％，然后对每个短时窗信号进行256点傅里叶变换，获得时频幅度谱|n(f,t)|。窗长的选择同时考虑了直升机声信号的线谱分辨率和环境噪声的短时非平稳特性。

3)采用非负矩阵分解方法对噪声时频幅度谱进行字典学习。nmf将已知非负矩阵x(此处为时频幅度谱|n(f,t)|)，寻找一种分解方法，使其近似等于两个非负矩阵d和c的乘积，即

x≈λ＝d·c(1)

式中，和矩阵d为字典矩阵(dictionarymatrix)，其列向量为基向量，描述不同谱模式，矩阵x中的列向量即是由矩阵d里面的基向量线性组合构成；矩阵c称之为编码矩阵(codematrix)，包含了由基向量构建矩阵x的组合方式。由于具有幅度不变性，音频应用中经常将kl散度准则作为矩阵分解的目标函数：

字典矩阵d和编码矩阵c依据乘法更新规则(multiplicativeupdatingrule)循环迭代获得：对于需要估计的参数θ，其迭代公式为其中即代表表达式中正的部分，代表表达式中负的部分。依据上述方法很容易推导得到矩阵c和d的迭代公式：

上式中.*和./分别代表两矩阵元素间的相乘和相乘，i是与x同大小、元素全部为1的矩阵。

对目标函数收敛性的判断一般可以通过计算目标函数的相对变化比值来完成，如果目标函数相对变化比值ε小于阈值，则可认为目标函数收敛[13]，如下式：

根据经验，取ε＝10^-3。

利用nmf算法可学习得到噪声字典矩阵，其算法流程如下：

a)取字典个数r为个，n0值根据实验效果确定，这里取n0＝5；通过[0,1]均值分布噪声初始化和进一步对矩阵d的每个基关于l1范数归一化；对cost(0)取较大的正实数，并根据公式(2)计算cost(1)；

b)n＝1；

c)

d)n＝n+1；

e)根据式(3)和(4)更新矩阵c和d；

f)根据公式(2)计算cost(n)；

g)endwhile；

h)输出矩阵c和d。

以背景噪声时频幅度谱|n(f,t)|为输入，可学习得到噪声字典矩阵，以表示dn。

2、直升机信号增强

1)对采集到的带噪直升机声信号x(t)分帧、加窗，加hamming窗，窗长取256点(对应0.25s)，帧间重叠50％，然后对每个短时窗信号进行256点傅里叶变换，获得时频幅度谱|x(f,t)|(下面简记为x)和相位谱∠x(f,t)。

2)利用sscnmf算法估计直升机声信号时频谱。假设环境噪声n(t)是加性噪声，且与目标信号s(t)不相关，则带噪目标信号x(t)可表示为：

x(t)＝s(t)+n(t)(6)

对带噪信号进行stft，并可将其分成纯净目标信号和噪声两部分：

x＝|stft(x(t))|(7)

式中，x为带噪信号短时幅度谱；s、w分别为纯净目标信号和噪声短时幅度谱；ds、cs表示纯净信号的基矩阵和编码矩阵，而dn、cn则表示噪声基矩阵和目标信号。如果可以分离出属于目标信号部分，通过他们的乘积：

使其还原出纯净目标信号，达到环境噪声抑制的目的。

sscnmf算法对学习的直升机信号字典的每个谱模式施加稀疏约束，以突出线谱分量，此时代价函数包含kl散度和谱稀疏约束之和：

cost＝costkl+λ·costss(10)

式中，m为矩阵x的行数(即频点数)，rs为信号字典ds基个数。λ为谱稀疏约束系数，根据实验数据处理，可取λ＝0.1。costss表示为字典ds每个元素的值与其方差的值指标，由于ds关于l1范数进行归一化，因而其方差越大其值越小，对于其离散程度越高。

保持dn固定，其他参数cn、cs和ds的迭代公式同样依据乘法更新规则推导得到：

利用sscnmf方法估计信号时频幅度谱，算法步骤如下:

a)确定信号字典基个数rs，可取与噪声信号相同的个数；通过[0,1]均值分布噪声初始化cn、cs和ds，ds进一步关于其l1范数归一化；对cost(0)取较大的正实数，并根据公式(10)计算cost(1)；

b)n＝1；

c)

d)n＝n+1；

e)根据式(13)～(19)更新矩阵cn、cs和ds；

f)根据公式(10)计算cost(n)；

g)endwhile；

h)输出矩阵

3)结合估计得到的信号幅度谱以及原始信号的相位∠x(f,t)，通过短时傅里叶逆变换得到增强的目标直升机信号

通过该直升机升信号增强方法，能够从复杂背景噪声中抽取出质量较高的目标信号，可作为目标检测、定向与识别的前端处理。

上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明，它们并非用以限制本发明的保护范围，凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。

对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。