声音特征提取--Demon谱分析（舰船水声）

DEMON(噪声包络调制检测)

1.调制

**调制：**调制是将能量低的消息信号与能量高的载波信号进行混合，产生一个新的高能量信号的过程，该信号可以将信息传输到很远的距离。或者说，调制是根据消息信号的幅度去改变载波信号的特性（幅度、频率或者相位）的过程。

一个信号包括了幅度、频率和相位，我们可以根据消息信号的幅度来改变载波信号的幅度、频率和相位，即我们所熟知的调幅、调频和调相。消息信号（基带信号）与载波信号进行混合产生的信号称为调制信号。

参考什么是调制、为什么要调制、调制有哪些？

2.DEMON 谱

当舰船航速达到某个极限时，螺旋桨的高速旋转将产生空化噪声，由螺旋桨产生的一系列周期信号被空化噪声调制，在频域上表现为将低频的离散线谱线性搬移到较高频率噪声所处的频率周围，这些周期信号在频域是离散的线谱系列，对应舰船的轴频、叶频以及谐波频率，由于舰船辐射噪声低频部分还有其它噪声源，很难从低频段直接提取出螺旋桨产生的线谱系列，而是要先对高频段的调制信号进行解调，通过低通滤波得到信号的低频包络谱，即DEMON谱，又称噪声包络信号识别。【1】

即由于螺旋桨叶片速率低频线谱成分是调制再螺旋桨空化噪声的高频连续谱上的，从调制谱的高频成分中能提取低频分量，即DEMON谱分析技术，获得螺旋桨空化噪声的调制线谱。【6】

目标辐射噪声信号模型可以表示为下式【1】：
$$x(t)=[1+f(t)]\ast n(t)+g(t)$$

• 周期调制信号f(t)

  f(t)在频域表现为螺旋桨旋转产生的叶片速率线谱。舰船航行，由于螺旋桨高速旋转，产生一些了周期信号，在频域上表现为离散的线谱系列，对应了舰船螺旋桨的轴频、叶频以及它们谐波的频率，一般以调制的形式存在高频段舰船噪声中。

  周期信号f(t)的线谱频率f_m=mns,f_m是线谱系列的第m次谐波频率，n为叶片数，s为螺旋桨轴频，单位Hz，使用DEMON谱分析可以将低频线谱从高频调制连续谱中分离出来，从线谱中可以提取螺旋桨轴频、叶片数、叶频等特征。

• 宽带信号n(t)

  n（t）为连续谱对应的宽带时域信号，主要对应螺旋桨空化噪声。舰船常规速度航行时螺旋桨都是空化的，且空化噪声为舰船辐射噪声连续谱的主要成分。空化噪声连续谱在100-1000hz存在一个峰值，峰值右侧，谱级随频率增大而下降，左侧随频率增大而上升，噪声谱峰值位置受航速和吃水深度影响。

• 线谱信号g(t)

  线谱对应的时域信号g(t),为周期叠加成分，包括往复运动产生的机械噪声、螺旋桨叶片共振线谱和叶片速率线谱、水动力引起的共振线谱。这些线谱频率大多在1000 Hz 以下

3.DEMON谱分析流程

对于螺旋桨噪声的分析，首先要对螺旋桨噪声目标频段进行带通滤波，通常用平方检波解调获得包络谱，低通滤波后经过快速傅里叶变换得到目标频段解调谱。DEMON谱具体分析步骤如下【2】：

（1）带通滤波，原始信号中携带特征信息的频带有限，截取具有调制分量的频带信号，作为 DEMON 谱分析的宽带信号。

（2）线性检波，利用正交解调的方式获取调制信号的包络信息，常用的有平方检波、绝对值检波。

（3）低通滤波，舰船辐射噪声的轴频、叶频及其谐波分类引起的调制主要分布在0~100Hz 范围内，故采用低通滤波截取有效信息。

（4）傅里叶变换，调制信息多表现为单频线谱及其谐波分量的调制，故而对其调制信息的频率分布的研究具有实际意义。简言之，通过傅里叶变换将已调制的低频解调谱转换至频谱，方便后续的调制信息提取。

想在DEMON谱上获得清晰、效果好的线谱，上述方法并不能满足需求，原因是噪声的不同频带能量不稳定性以及信道对辐射噪声信号的畸变影响。水下辐射噪声信号在不同频率范围内调制强度不同；不同频带的DEMON谱线相互之间也不一样，有的频带内会发生线谱遗失，有的频带内的线谱并不是真正的谱线。同时，得到的DEMON谱信号中有谱成分不明显的问题，这是由于噪声中存在连续谱，从而导致部分距离比较近的线谱无法有效区分，给整个DEMON图谱带来模糊性。通常用分频段综合的DEMON提取方法解决上述问题。【4】【5】

分频段综合：

将目标的辐射噪声信号分成几个频段，再将各个频段的DEMON处理结果结合，由于这样得到的DEMON谱中存在连续谱，使得线谱容易被覆盖，需要在噪声的DEMON谱中减去连续谱，并对得到的结果进行线谱净化处理，最后可得到待分析目标辐射噪声较为清晰的线谱。

线谱净化：

可使用α双向滤波器，其作用是隔开两种不同类型的谱，流程如下：用该滤波器得到待分析目标辐射噪声DEMON谱的趋势项，为了使线谱和连续谱分离要设置一个门限，门限的值设置为目标 DEMON 谱趋势项与一个特定数的乘积，过门限的谱线保留，没有过门限的谱线去除，便可达到线谱、连续谱分离目的，得到理想的线谱。

α双向滤波器：

是一个一阶递归滤波器，如待平滑的目标辐射噪声信号为s(n),其长度为N，当n=1，2，…,N-1，m=1,2,…,N时，α双向滤波器原理如下式：

对信号序列进行DEMON分析的具体实现流程如下[7]【4】【5】：

1. 选择目标噪声可能具有较强调幅的频带，对于船舶水下辐射噪声来说，一般选择大于1kHz某频带范围信号，将该频带内的信号按照平均或比例或频率范围等规则分成 N个频带 ，在时域中进行滤波后 ，再基于平方检波算法对各频带信号进行处理 ，从而得到包络信号 ，最后对包络信号进行快速傅里叶变换 FFT，获得该频带的ＤＥＭＯＮ谱。假设为P_i,i=1,2,…N

目标的轴频、叶频、叶片数等信息，关系如下：【3】

参考资料

【1】舰船辐射噪声特征提取与识别研究-牛芳
【2】基于潜标平台的舰船目标探测技术研究-张亚斌
【3】基于特征自动提取的水声目标识别方法研究-丰玉霖
【4】水下小型UUV辐射噪声特征提取技术研究-张洪岩
【5】船舶水下噪声特征获取方法研究-刑军华
【6】舰船目标识别技术研究-沈广楠
【7】目标通过特性的LOFAR及DEMON 分析-朱世才