刘 倩，李 梅

(上海船舶电子设备研究所，上海 201108)

0 引 言

随着近年来声呐探测技术的发展，如何进行水下防御成为了国内外的研究热点。主动声呐利用发射脉冲信号在目标上的回波来探测目标，为了防止敌方探测到我方目标，就要削弱或者完全抵消目标回波。通过估计目标回波的参数，包括中心频率、初相位等，利用估计得到的参数发射反相信号可对目标回波进行抵消，但当信噪比降低时，参数估计会产生较大的误差，从而降低对目标回波的抵消性能。如何在低信噪比情况下提高目标回波参数估计性能是急需解决的问题。

传统的降噪方法有相干平均和滤波，相干平均通过将多次观测信号相加取平均以抵消噪声，该方法虽然简单、运算量小，但对非平稳随机噪声的降噪效果不佳，并且该方法需要大量的观测数据，处理进程较长。滤波方法是假定信号和噪声处在不同的频带上，通过滤波器将通带外的噪声进行滤除，但在未知信号频率范围时，无法设计滤波器参数。近年来，国内外学者提出许多通过将信号分量与噪声分量进行分离从而降噪的方法，如奇异值分解[1]（SVD）、独立分量分析[2-3]（ICA）、小波变换[4-7]、经验模态分解[8]（EMD）等。奇异值分解利用含噪信号构造Hankel 矩阵，对此矩阵进行奇异值分解，所得到奇异值中较大的对应信号成分，较小的对应噪声成分，但当信噪比较低时，信号成分与噪声成分奇异值的分界不明显，导致信号与噪声难以分离。ICA 通过构造解混系统将信号分量与噪声分量进行盲分离达到降噪目的，但由于ICA 的不确定性，分解出的信号分量幅度和相位可能产生变化。小波变换利用不同尺度的小波基将信号分解成许多子带信号，将信号能量集中到某些频带的少数系数上，通过将其他频带上的小波系数置零或是给予小的权重，即可达到抑制噪声的目的。然而，小波基的选取及分解层数的确定存在一定的盲目性。EMD 是根据信号自身特性，将一个复杂的信号自适应地分解为一系列固有模态函数（I M F），通过将噪声所在的I M F 置零达到降噪的目的，但当信噪比较低时，EMD 分解容易产生模态混叠问题，导致信号分量与噪声分量难以分离。为解决模态混叠问题，Wu[9]和Yeh[10]分别提出了集总经验模态分解（EEMD）和互补集总经验模态分解（CEEMD）算法。EEMD 虽然可以通过向信号中不断添加白噪声改变信号极值分布来避免模态混叠问题，但其集总平均次数一般要在几百次以上，非常耗时，并且会在每个IMF 分量中残留一些白噪声。而CEEMD通过向信号中添加许多对符号相反的白噪声，可以在集总次数较少的情况下解决模态混叠问题，大大节约计算时间，并且可以基本上去除残留在IMF 分量中的白噪声。

1 基于CEEMD-MI 的目标回波参数估计

1.1 算法流程

利用CEEMD-MI 对目标回波进行参数估计可分为3 个步骤。首先利用CEEMD 对目标回波进行自适应分解，再通过MI 重构目标回波，最后对重构后的信号进行参数估计。

1.2 基于CEEMD-MI 的目标回波去噪

1.2.1 CEEMD 算法

互补集总经验模态分解（CEEMD）是依据信号自身特点来对信号进行自适应分解的方法，其原理是通过向信号中施加N组正负噪声，不断地减去由信号局部极大值点和局部极小值点构造的包络均值，当满足IMF 条件，即极值点数目和过零点数目相等或最多差一个并且在该时间区间任一点包络均值为0，即可分离出添加正负噪声后信号的各个IMF，对其取平均。最后将N组噪声情况下分离出的IMF 进行集总平均，即可得到最终的IMF 分量，具体步骤如下：

步骤1对信号x(t)分别施加N对正负噪声ui(t)，即

可得到加入正负噪声后的信号为：

步骤2对信号和，分别找出其局部极大值点和局部极小值点，利用三次样条函数分别对这些局部极大值点和局部极小值点进行插值得到的上、下包络和以及的上、下包络和，可得的包络均值为的包络均值为分别减去它们的包络均值可得：

行走式压电驱动器工作原理如图1所示[2]，其位移输出机构(动子)一般由两侧箝位机构和中间驱动机构组成，当向右运动时：左箝位压紧，右箝位放松；驱动机构伸长，动子向前移动一个步距；右箝位压紧，左箝位放松；驱动机构恢复原长度，左箝位压紧，右箝位放松，恢复初始状态。

步骤3同理，找出和的局部极大和极小值点，利用三次样条函数对其插值得到上下包络后分别减去包络均值和，得到：

步骤4令

再对加入N组噪声下得到的IMFik做集总平均，可得：

式中，k=1,2,...,m，IMFik是对x(t)经CEEMD 分解后得到的m个固有模态函数。

1.2.2 互信息

互信息（MI）是2 个随机变量间相互依赖性的度量，可以看成是一个随机变量中包含的关于另一个随机变量的信息量。假设存在一个随机变量X和另一个随机变量Y，其互信息为：

式中：p（x,y）为X和Y的联合概率密度分布；p（x）和p（y）分别为X和Y的边缘概率密度分布。x(t)经CEEMD分解后得到的各个IMF 分量中，噪声分量所在的IMF 与x(t)之间的互信息量小，信号分量所在的IMF 与x(t)之间的互信息量大，计算各个IMF 与x(t)之间的互信息，选取互信息量大的IMF 分量进行信号重构，可去除x(t)中的噪声成分，达到去噪目的。

1.3 目标回波信号参数估计

令重构后的信号xR(t)=Asin(2πf0t+φ0)，其Hilbert 变换为：

由式（19）可得xR(t)的解析信号为：

2 仿真分析

2.1 降噪性能分析

假设信号x(t)的频率为10 kHz，初相位为0 rad，脉宽为3 ms，采样率为200 kHz，信噪比为10 dB，图1为去噪前信号的时域波形。

图1 去噪前信号Fig. 1 Signal before denoising

对x(t)进行CEEMD 分解，得到的各个固有模态函数的时域波形及频谱如图2 所示。

图2 固有模态函数Fig. 2 Intrinsic mode function

图2 中IMF1～IMF6 分量与x(t)的互信息量如图3 所示。可知，IMF3、IMF4 分量与x(t)的互信息量较大，并且从图2 频谱分析可以看出，信号成分主要分布在IMF3与IMF4 分量中。故选取IMF3 分量与IMF4 分量对信号进行重构，重构后的信号xR(t)时域波形如图4 所示。

图3 各个固有模态函数分量与原信号互信息Fig. 3 Mutual information of IMFs and original signal

图4 重构信号（CEEMD-MI 去噪处理）Fig. 4 Reconstructed signal（CEEMD denoising）

将图4 与图1 对比可得，通过CEEMD-MI 算法对信号进行自适应分解，根据各个IMF 分量与原信号互信息量的大小对信号进行重构可在很大程度上去除原信号中的噪声分量，得到较好的降噪效果。

2.2 不同信噪比下参数估计性能分析

假设信号x(t)的频率为10 kHz，初相位为0 rad，脉宽为3 ms，采样率为200 kHz。当信噪比分别为10 dB、8 dB、5 dB 和3 dB 时，分别对未进行去噪处理的信号和经CEEMD-MI 去噪处理后的信号进行参数估计，估计所得的中心频率、初相位如表1 所示。可知，当信噪比较低时，未进行去噪处理估计所得的信号中心频率及初相位与真实值偏差较大，而利用CEEMD-MI 进行去噪处理后，估计所得的信号中心频率及初相位与真实值偏差很小, 中心频率估计误差小于0.2%，初相位估计误差小于2.5%。

表1 估计信号参数Tab. 1 Estimated signal parameters and echo cancellation

如图5 所示，当信噪比是10 dB 时，实线和虚线分别代表利用表1 中未进行去噪处理和利用CEEMDMI 去噪处理后估计所得参数生成的信号与x(t)的差值，前者与x(t)的差值幅度很大，而后者与x(t)的差值幅度较小。由此可得，经CEEMD-MI 算法对低信噪比信号进行去噪处理后，参数估计性能得到了很大程度上的提高。

3 试验数据验证

图6 为水听器接收信号图，图7 为对目标回波未进行去噪处理及经CEEMD-MI 去噪处理后，利用估计所得参数生成的信号对与目标回波信号的差值。从图中可以得到，在低信噪比情况下，对目标回波未进行去噪处理估计所得参数生成的信号与目标回波差值很大，由于中心频率估计的不准确使得估计信号与回波信号的差值产生了其他频率的信号成分，而利用CEEMD-MI 对目标回波进行去噪处理后估计所得参数生成的信号与目标回波差值很小，几乎只剩噪声成分。

图6 水听器接收信号Fig. 6 Hydrophone receive signal

图7 估计信号与目标回波差值Fig. 7 Difference between estimated signal and target echo

4 结 语

本文提出一种基于CEEMD-MI 的目标回波参数估计算法，该算法首先利用CEEMD 对低信噪比情况下的目标回波进行自适应分解，再利用MI 对分解得到的各个IMF 进行筛选并进行信号重构，有效地去除了信号中的噪声成分。最后对重构后的信号进行参数估计。仿真结果和湖试试验结果表明，本文所提算法可在低信噪比情况下对目标回波参数进行有效估计。