基于3阶小波张量的船舶辐射噪声识别

郭德鑫，康春玉，夏志军，张 忆

(1. 海军大连舰艇学院 研究生队，辽宁 大连 116018；2. 海军大连舰艇学院 水武与防化系，辽宁 大连 116018)

0 引 言

利用船舶在航行过程中产生的辐射噪声进行特征提取来实现目标识别是水下目标分类识别的一种重要手段，但是由于船舶型号、航速、螺旋桨工况的不同，目标的辐射噪声特征变化较大，特别是在恶劣海况下，识别效果更不令人满意，提取更加具有鲁棒性的特征及研究与特征相匹配的分类识别方法是急需解决的问题。

近年来，随着张量研究的不断深入与发展，张量分解在图像识别，语音信号处理和医学等领域[1]中得到广泛应用。在图像识别领域，周春光等[2]在线性主成分分析的基础上提出了稀疏张量主成分分析算法，有效降低了遮挡对人脸特征提取的影响，并且在不同表情的人脸识别中优于原算法；计雨含利用奇异值分解对人物、表情、向量长度构建的3阶张量进行识别，提高了人脸在光照、角度发生改变时的识别精确度[3]和分类识别时的速度[4]；同时在循环次数、特征维数以及训练样本和测试样本数目差别等不同方面对非负张量分解和非负矩阵分解进行的对比也表明了非负张量分解具有的有效性和高效性[5]；任卿龙[6]利用不同波段上的遥感图像构建高维张量，对飞机目标和舰船目标进行分级识别也取得了较好的实验结果。在声学领域，王珊[7]在混合矩阵估计中引入张量分解，利用分解后的各个矩阵实现了对多通道混合音乐的盲分离；杨立东[8-9]用语音命令张量经过Tucker分解得到的特征进行识别，平均识别率达到95%，同样对声矢量传感器输出的多维阵列信号进行张量奇异值分解，相比于传统的矢量奇异值分解能够更好地抑制噪声，并对信号子空间的估计更加准确[10]。医学方面，利用张量方法表征脑电信号的时域、频率、空域特征[11]，处理肌电信号[12]，提取心电信号[13]的操作，均比传统方法具有更高的准确率和相对的优越性。机器学习方面也有关于张量应用的研究[14-15]。

张量分解在图像识别，语音信号处理和医学等领域取得了较好的成果，但在与其具有类似高维结构的船舶辐射噪声信号特征提取和分类识别领域中尚未见到相关报道。本文在船舶辐射噪声特征提取中引入张量，利用张量的空间坐标系具有旋转不变性，以及直接对高维空间进行整体分析和操作的特点进行船舶辐射噪声特征提取和张量构造，并经过张量分解实现对船舶辐射噪声的分类识别。

1 张量分解

张量可以直观理解为一个多维数组，是多个矢量的空间乘积，张量的阶是张量不同维上元素个数的数量表达，也被称为way或mode。通常认为向量是张量的1阶表达,矩阵是张量的2阶表达，具有3阶或更高阶数的张量则称为高阶张量。本文提及的张量指的都是高阶张量，用粗体大写字母表示，如；向量用粗体小写字母表示，如；矩阵用斜体大写字母表示，如。一个N阶张量可以表示为其中表示第i阶上的维数。

张量分解的思想最早由Hitchcock在1927年提出。1963年，Tucker首次提出了Tucker分解，并在随后进行了完善。1970年，Carroll，Chang和Harshman提出了CANDECOMP/PARAFAC分解，即CP分解[1]。近年来，尹凤等[16]在张量CP分解和Tucker分解的基础上，提出一种利用交替最小二乘法计算高阶张量秩的分解方法。Tucker分解将张量分解为核张量和因子矩阵的形式，类似于矩阵的奇异值分解，因此Tucker分解又被称为高阶奇异值分解（HOSVD）[17]，并可以将其看做是奇异值分解（SVD）的高阶形式。本文所采用的张量分解方法就是Tucker分解。

图1 三阶张量Tucker分解模型Fig. 1 Tucker decomposition model of a third-order tensor

同奇异值分解一样，高阶奇异值分解中核张量代表原始张量的某些特性，因子矩阵代表着这些特性的重要程度，核张量和因子矩阵也是本文进行舰船辐射噪声分类识别的关键。

2 小波变换张量特征

研究表明，船舶信息主要集中在1 000 Hz以下的船舶辐射噪声中[18]，而1 000 Hz以上主要是螺旋桨空化噪声，因此对低频噪声进行分解可以得到更多的船舶信息，更便于船舶特征的提取。离散小波变换可以对信号低频进行分解，因此本文采用的方法是将船舶辐射噪声经过分帧、小波变换、美尔频率倒谱系数（MFCC）提取处理后构建一个帧数×小波分解分量×MFCC特征维数的3阶辐射噪声张量。

2.1 分帧

2.2 小波变换

小波变换（WT）是在傅里叶变换（FT）的基础上发展来的，将FT中无限长的三角函数基换成有限长的会衰减的小波基，能够在获得频率信息的同时定位到时间，其与信号做卷积的形式表示为下式：

对信号进行小波变换是改变小波的尺度，实现在频域内的伸缩；改变平移量实现在时域内的移动，通过分析与在时刻的相似程度得到该时刻的频率信息。在信号处理中，一般采用离散小波变换，把式(4)中的参数，都进行离散化，令其中，通常取分别对应小波的伸缩变换与平移，从而把连续小波变为离散小波(DW)[18]。

则离散小波变换(DWT)为：

图中，L为低通滤波器，H为高通滤波器，↓2为系数为2的降采样滤波器。由于一直在进行降采样，所以虽然滤波器L，H不变，但其滤波带宽一直在减半。如此经过次分解后，信号就可以表示为：

在小波变换中，Morlet小波，Haar小波，Daubechies小波，Meyer小波等是常用的小波，本文进行小波变换时采用的是Daubechies小波，常记为是其阶数。小波不具有解析式，由满足一定条件的滤波器系数迭代逼近小波，一般先给出低通滤波器系数再由式（8）推导得到高通滤波器系数。

图2 离散小波变换分解流程Fig. 2 Decomposition process of discrete wavelet transformation

2.3 MFCC特征提取

MFCC特征是一种根据人耳听觉对不同频率的声波灵敏度不同的特性提出的一种听觉感知特性，在信号处理中可以理解为信号的能量在不同频率范围的分布系数占比。计算过程如下[19]：

将频谱幅度平方后得到的能量谱通过美尔（Mel）频率滤波器组得到各分量的Mel频谱,并对Mel频谱取自然对数得到对数频谱：

2.4 张量特征构造

通过上述处理，就可以构建一个帧数×小波分解分量×MFCC特征维数的3阶船舶辐射噪声张量，其结构如图3所示。

图3 三阶张量Fig. 3 A third-order tensor

3 张量特征分类识别

4 船舶辐射噪声分类识别实验

分类识别实验中，共有Ⅰ，Ⅱ，Ⅲ类船舶4 506个辐射噪声样本，采样频率为25 kHz，每个样本长为6.553 6 s。从整个辐射噪声样本中每间隔10个样本取1个作为训练样本集，余下的作为测试样本集，得到451个训练样本，4 055个测试样本。对每个样本采用汉明窗进行分帧，每帧长度为81 92个采样点，帧移为1 024个采样点，分帧后共得到153帧。在小波分解过程中，采用db3小波进行3层小波分解，小波分解后得到4个分量。对每个分量提取12维美尔频率倒谱系数（MFCC）。然后使用上述步骤得到的153帧、4个小波分解分量、1 2维M F C C特征系数构建一个的3阶小波张量，最后按分类识别方法完成分类识别实验，测试样本集总的正确分类识别概率为详细的分类识别概率如表1所示。

表1 分类识别概率（%）Tab. 1 Probability of classification and recognition (%)

5 结 语

辐射噪声特征提取是进行目标识别的重要组成部分，提取能够充分表示噪声信息的特征是提高识别率的关键。本文提出的舰船辐射噪声张量特征提取方法，利用小波变换及美尔频率倒谱系数（MFCC）构建具有较高鲁棒性的3阶张量。通过HOSVD张量分解得到训练样本的核张量和各维度上的因子矩阵，将因子矩阵的转置与测试样本构建的张量进行矩阵乘运算得到投影张量，最后通过范数比较核张量和投影张量之间的差异性，确定测试样本所属类别。实验结果表明该方法识别率高，具有一定的优越性。