改进间接滤波噪声抵消算法在VLF 中的应用

贾书阳，蒋宇中，牛 政，刘 刚

(海军工程大学 电子工程学院，湖北 武汉 430033)

0 引 言

虽然甚低频通信不如卫星通信等信息传输率高，但由于海上降雨频繁，卫星通信的高频段信号波长与雨滴直径较为接近，雨水会对信号造成严重衰减[1]。而甚低频具有传播距离远，在一定程度上能穿透海水及地表，传输信息稳定的能力等特点，是恶劣天气条件下卫星通信失效时，一种陆地对远洋航行中船舶传递信息的备用手段。因此，开展甚低频信号高速传输的研究，对用于对潜通信的甚低频通信技术未来发展具有一定指导意义[2]。

甚低频的通信频率为3～30 kHz，其实现传输过程中的不足为噪声干扰。超低频远洋通信成功的关键是进行有效的干扰抵消。远洋通信在接收端受到的干扰主要包括：通信信道内的非合作人工干扰（广播信号、工业干扰、大气噪声等）[3]。另外，船舶周围也存在强烈的低频干扰噪声，主要与船舶轴频电场、防腐系统、大电流设备等电流脉动有关。由于噪声形态在频域、时域和空域呈现随机变化的非平稳状态，要实现船舶电磁噪声抵消，特别是要实现很大的噪声抑制比，技术难度大，需要在多个学术领域开展多基础理论与相关技术创新研究[4]。

关于电磁噪声特征及其抵消技术，美国和俄罗斯都曾开展过较为深入的研究。然而仅从国外文献获悉有关电磁噪声较多的描述，无具体数据和技术可借鉴[5]。本文首先将信号通过带通滤波器有效滤除带外噪声。其次，对于带内噪声，采用改进的间接滤波算法进行降噪技术抑制带内噪声，传统的自适应滤波算法。

本文搭建接收信号所用到的天线模型，即一种垂直双通道的天线接收模型，对信道接收的信号进行了简单的建模；针对传统的滤波算法无法直接恢复期望信号的问题，介绍间接滤波算法的基本原理，以及结合小波变换进一步改进了间接滤波算法；针对本文提出的2 种算法进行仿真和实验验证。

1 信号与天线模型

1.1 天线模型

本文提出一种垂直双通道的天线接收模型布设方法。如图1 所示，主接收天线m与参考天线n在空间中正交放置。采用两互相垂直的天线接收信号，主天线正对信号源，参考天线平行于信号源。由于磁性天线增益最大的地方在磁芯方向，而正交于磁芯方向增益最小。主天线m直接接收到信号，参考天线n中期望信号的成分最少，理论上能接收到除了主天线方向以外空间中的所有干扰信号，但信号受到教研室有限环境影响发生反射，所以参考天线n也会收到期望信号，但十分微弱可忽略。

图1 数据传输硬件系统结构图Fig.1 Structure diagram of data transmission hardware system

1.2 信号建模

由图1得出，主天线在i时刻接收到信号x(i)可表示为：

其中：e(i)表示期望信号；n(i)表示背景噪声。船舶在航行中在甚低频频段通信的背景噪声的主要成分是大气和海浪噪声、电机工作噪声以及其他通信用户干扰。通信时这些部分干扰必须被有效抑制，主天线在甚低频频段才能保证很大容量的通信。

参考天线与主天线互相垂直，接收到的信号几乎不包含期望信号成分。且由于磁性天线间相距较近，远小于甚低频电磁波信号的波长，在此忽略两天线之间接收信号的延时带来的影响。则参考天线在i时刻接收到信号y(i)即为噪声信号n(i)。

2 干扰抵消算法

2.1 间接滤波算法

一般地，自适应滤波算法以期望信号作为滤波器输出，当信噪比极小时，往往因为期望信号微弱而导致滤波器系数偏差较大，很难恢复出理想的期望信号。与以往滤波算法不同的是间接滤波算法以噪声作为自适应滤波器输出信号。假设信号与噪声之间互不相关，由于接收机接收到的序列中信号功率小于噪声的功率，传统的自适应滤波算法很难将信号从噪声中恢复出来。本文分别在传统的最小均方自适应滤波算法（MLSA）的基础上，提出一种间接最小均方自适应滤波算法，并对处理后的结果运用小波变换分解，进一步抑制噪声提高通信信噪比（SNR），解决双通道下的低信噪比的滤波问题。间接自适应滤波器的原理框图如图2 所示。

图2 间接自适应滤波器原理图Fig.2 Schematic diagram of indirect adaptive filter

2.2 改进的间接滤波算法

改进后的间接滤波算法如图3所示。算法改进主要有两方面：为减少求解滤波系数向量计算量，将厄米特矩阵二次型应用到间接滤波算法中；同时，引入小波变换分别对自适应滤波器后的期望信号(i)及副天线接收的噪声信号n(i)进行分解重组，以进一步提升期望信号的信噪比。

图3 改进的间接滤波算法原理图Fig.3 Schematic diagram of improved indirect filtering algorithm

2.2.1 厄米特二次型分解的间接滤波算法

2.2.2 小波变换分解

小波变换算法可以很好地对信号的细节进行分析，其在降噪应用领域十分广泛。小波基的选取对降噪效果影响很大，但目前只是往往采用经验和不断试验的方法选取[8]。本文小波基的选取采用Daubechies 小波，简称db 小波，记为 φ(i)。之所以选用db 小波的原因是db 小波更接近于正弦信号的波形，能更好分解重构原始信号。db 小波随着阶次增大消失矩阶数也越大，有利于频域的分析，但同时也牺牲运算时间的代价，实施性较差[9]，因此这里选用db10 小波做基波，如图4 所示。将小波基进行伸缩平移后得到一个小波序列[10]:

图4 Db10 小波基函数与频谱Fig.4 Db10 wavelet basis function and spectrum

其中：a为伸缩因子，b为平移因子。任意函数的连续小波变换为该函数与某小波序列的卷积。

2）噪声信号n(i)n层小波分解；

改进的间接滤波算法输出信噪比定义为经过自适应滤波算法处理后的期望信号xr(i)平均功率与残留的噪声干扰nr(i)功率之比：

3 仿真及实测结果验证

为测试算法的效果，本文分别用模拟仿真和实验室测量2 种方法进行验证，并与原始滤波算法分别进行比较分析。

3.1 仿真

本文仿真数据实现的具体过程为：分别生成幅度为0.1 V、频率为20 kHz 的期望信号，信号的采样频率为1 MHz。噪声信号服从标准正态分布，为使信号完全淹没于噪声之中，设定信噪比S NR=−20 dB。仿真主天线收到的信号为期望信号与噪声的叠加，副天线为同样均值方差的正态分布的噪声。图5（a）为仿真主天线和副天线接收到的信号功率谱。选取滤波器的系数为20 阶，将间接滤波算法与改进的滤波算法在仿真条件下的噪声抑制效果呈现如图5（b）所示。

图5 仿真信号功率谱Fig.5 Power spectrum of simulation signa

通过观察图5（b）发现，在以上2 种仿真算法实验中期望信号的功率基本相同，而在第2 种方法中噪声功率得到了明显的抑制，可见小波变换在一定程度上抑制了信噪的高频信号，提高了SNR。

3.2 实验

为进一步验证提出算法的效果，在实验室环境中搭建发射信号装置以及接收信号摆放天线结构。

信号源有信号发生器产生，频率与仿真信号频率相同，幅度为100 mV。主副天线分别与FPGA 搭建的信号接收设备相连，其采样频率与信号模拟频率相同，能将接收到的数据存入电脑以分析处理。信号源的发射线圈距离主天线中心为5 m。

图6（a）为实测主天线和副天线接收到的信号功率谱。图6（a）中，通过与仿真信道的对比发现信号的频率为2.016 kHz，并非理想的2 kHz，原因可能是FPGA 作为接收设备的晶振频率存在一定的偏差导致分频的采样频率不精确。

图6 实验信号功率谱Fig.6 Power spectrum of experimental signal

由图6（b）可知，两类滤波算法对于实验室环境下噪声信号干扰均存在抑制效果，但改进的间接滤波算法能够更好地将低频中高于信号的部分降低，能够在信噪比较低的情况下抵消更多的干扰信号，突出信号的功率。

将2 种算法的仿真与实验信噪比计算测量各10 次后取平均结果列入表1。模拟数据中理论计算与测量结果偏差较大，可能是小波重构的信号功率不及原信号的功率大引起的；其次可能是仿真噪声采用高斯白噪声没有现实噪声复杂真实；再次是受统计次数不足，测量结果有误差以及保留有效数位影响。实验测量数据与理论计算值较为接近，验证算法可行性。虽然改进的算法性能提升效果不明显，但是从功率谱不难看出，改进的算法能更好地抑制部分噪声功率。

表1 两种算法的信噪比10 次试验结果平均值Tab.1

4 结 语

本文描述了甚低频通信时存在的问题，结合小波变换算法提出改进的间接滤波算法来解决问题。首先搭建垂直双通道的天线接收模型，简单介绍了间接滤波器算法和小波变换算法并结合对算法提出了改进。提出的改进的间接滤波算法算法结合了小波变换，仿真和实验结果表明，改进后的算法相比于传统算法，能够在信噪比较低的情况下抵消更多的干扰信号，比普通间接滤波算法信噪比更高。2 种算法计算与实验数据基本吻合，验证了算法的可行性。本文存在的主要问题是只在实验室进行测量验证，未在远距离通信、更复杂的干扰环境下对算法进行实验，下一步将对此进行改进完善。