LMS自适应滤波在航空瞬变电磁信号去噪中应用.docVIP

LMS自适应滤波在航空瞬变电磁信号去噪中应用 　　摘 要：以 LMS自适应滤波理论为基础，借助Matlab编程软件成功滤除了航空瞬变电磁信号中的噪声，为航空瞬变电磁信号的后期正反演研究提供了可靠的数据。通过实验证明，LMS自适应滤波器能够合理、有效地去除天电噪声和地质噪声。 　　关键词：航空瞬变电磁信号；数据去噪；滤波方法；LMS自适应滤波 　　中图分类号：TP39；P631.3+26 文献标识码：A 文章编号：2095-1302（2018）01-00-03 　　0 引 言 　　航空瞬变电磁法（Airborne Transient Electromagnetic Method，ATEM）是20世纪中期问世的一种快速普查良导电金属矿的航空物探方法，其具有速度快，勘探成本低，探测范围广，可有效抑制复杂地形影响等优点，已成为国内外广泛使用的一种地质勘探方法[1]。但航空瞬变电磁勘探法主要观测的是二次场信号，其有效信号幅值弱，频带宽[2]，因而测得的二次场信号的真实性或准确性难以保证。此外，航空瞬变电磁法在实际应用中受噪声干扰的影响更为严重，甚至得到的观测信号是已被噪声信号掩盖的无用信号[3]，导致后期处理进入错误的方向。为保证实测信号的可靠性，对测得的数据必须先做滤波、去噪等预处理。 　　航空瞬变电磁信号的噪声类型分为天电噪声、仪器振动噪声、地质噪声以及人文噪声等[4]。由于平均滤波、中值滤波等滤波方法对天电噪声的滤除不够理想，本文提出运用LMS自适应滤波方式滤除其高频范围内的天电干扰和地质噪声。 　　1 LMS自适应滤波原理 　　自适应滤波（adaptive filtering）是信号处理领域一个非常重要的分支。自1959年Widrow提出自适应的概念以来，自适应滤波理论一直受到普遍关注，并得到了不断发展与完善。信号处理理论和应用的发展为自适应滤波理论提供了必要的理论基础，其已在通信、雷达、自动控制、图像与语音处理等领域得到了广泛应用[5]。 　　顾名思义，自适应滤波器是一种能够根据输入信号自动调整自身性能并?M行数字信号处理的数字滤波器，其最本质的特点在于具有自学习和自调整的能力，即自适应能力。与固定滤波器相比，自适应滤波器能够根据当前自身的状态和环境自动调整和校正当前滤波器的参数。由于信号与噪声的时变特性未知，因而其统计特性不确定，故可认为其是最优滤波方法。自适应滤波器具有可调整系数的滤波结构与可调整和校正滤波器系数的自适应算法[6]。 　　Widrow等人提出的最小均方算法（Least Mean Square，LMS）是一种以期望响应和滤波器输出信号之间误差的均方值最小为准则，依据输入信号在迭代过程中估计梯度质量，并更新权系数以达到最优的自适应迭代算法。LMS算法是一种梯度最速下降算法，其显著特点在于其简单性。这种算法无需计算相关矩阵及矩阵运算[5]。 　　LMS算法是一种线性自适应滤波算法，包括两个基本过程，即滤波过程与自适应过程。在滤波过程中，自适应滤波器计算其对输入的响应，并通过与期望响应比较，得到估计的误差信号。在自适应过程中，系统估计误差自动调整滤波器自身的参数。这两个过程共同组成一个反馈环，LMS自适应滤波原理如图1所示。 　　其中，s （k）是输入信号，d（k）是期望信号，y（k）是滤波后的信号，e（k）是误差信号。 　　该滤波器根据e（k）和s（k），通过自适应算法找到E[e2（k）]最小时滤波器的权值，从而找到滤波系数权值，实现自适应滤波。 　　LMS算法系数的更新规则是通过迭代寻找极值，即最速下降法，该方法后一时刻的系数由前一时刻的系数通过一定的迭代运算得出： 　　式中，W为滤波系数，k为当前迭代值，Wk和Wk+1分别代表k时刻和k+1时刻的自适应系数，X为待滤波数据，μ为放大系数，称为收敛因子，用于调整自适应迭代的步长，ε是误差信号（期望值与预测输出值之差），为误差信号平方的梯度。假设当前为k时刻，其后一时刻为k+1时刻，则k+1时刻的系数为k时刻系数加上归一化函数的负梯度。系数函数也称为性能函数，等同于期望信号的均方误差。从（3）式可得出，LMS算法与平方、平均、矩阵逆变换无关，但其收敛判断标准与最速下降方法确定的值相等。对μ的标准由下式得出： 　　n为滤波系数的阶数，R为输入的相关矩阵E[XkXkT]，tr为矩阵迹。滤波系数的阶数n和滤波增益μ是失调M和自适应处理时间常数τmse的折中。 　　2 LMS自适应滤波实现流程 　　图2是运用LMS自适应滤波算法对航空瞬变电磁去噪的实现流程[7]。由图可知，读入数据后，首先计算信号长度和功率，然后计算出步长值并将该值作为初始步长存入算法中，设置好前两个抽头系数值后进入算法的自适应计算环节，首先计算前两个值，之后再计算误