3第三章_自适应数字滤波器.ppt

LMS算法的matlab实现 clc clear all close all length = 1024 * 8; N = 200; x = zeros (1, length + N); for i = 0:0.005:0.5 if i==0 x = 0.5* cos(2* pi * i* (0:length + N - 1)); else x = x + cos(2 * pi * i * (0:length + N - 1)); end end x = x;%输入信号 d = zeros( 1, length + N) ;%期望信号 for i = 0:0.005:0.2 if i==0 d = d + 0.5*cos(2*pi* i*(0:length + N-1)); else d = d + cos(2 * pi * i * (0:length + N- 1)); end end d = d; w = zeros( N, 1) ; E = zeros(1, length) ; mu = 0.00001; for i = 1: length y(i)=x(i:i + N-1)*w;%滤波器的实际输出 E(i)=d(i)-y(i);%估计误差 w = w + 2*mu*E(i) *x(i:i + N- 1);%迭代计算权矢量 end figure (1); plot( (1: length) , E) ; title(误差) ; 图 3.4.4 心电图电源干扰的消除 (a) 原始输入； (b) 参考输入； (c) 噪声对消输出 2) 胎儿心电监护 监测胎儿心音，需要将母亲的心音及背景干扰去除。母亲的心音强度通常是胎儿心音的2倍到10倍，肌肉活动及胎儿运动产生的背景干扰也常常大于胎儿心音的强度。可以采用自适应对消的方法来监护胎儿心电图。将母亲胸部的信号作为参考输入，主要包括母亲心音和背景噪声，也有胎儿的心音信号，如图3.4.5(a)所示；从母亲腹部取出的信号作为原始输入，包括母亲心音、胎儿心音和背景干扰，如图3.4.5(b)所示，可以看出，这个系统是一个有信号分量泄漏到参考输入端的自适应对消系统。对消后的输出如图3.4.5(c)所示，可以看到，母体的心电信号的强度小于胎儿的心电信号，并且胎儿的心电信号比较清楚。 Ultrasound Doppler Device Fetal Scalp ECG Fetal Heart Sound ?ACF ?Peak Detection?FHR Internal monitoring Intrauterine Probe 图 3.4.5 胎儿心电的监护 （a） 胸部导联心电（参考输入）；（b） 腹部导联的心电（原始输入）； (c) 自适应对消器的输出 ? 3) 在有他人讲话的背景中提取他人讲话 可将靠近所需要的讲话人的扩音器的输出送至原始输入端， 而将远离讲话人而靠近干扰源（其他人的讲话）的扩音器的输出送至参考输入端，则在抵消器的输出端可以提取消除了干扰的讲话。 4) 用来作为天线阵的自适应旁瓣抑制器 如图3.4.6所示，假定干扰的来向已定，首先将天线阵的主瓣对准拟接收的信号来源，此时，从点干扰源来的干扰通过天线方向图旁瓣与信号一起进入接收机，将它作为自适应抵消器的原始输入。然后，取天线阵中的若干根，组成的方向图正好在信号的入射方向有一个零响应谷点，这个组合只能得到干扰而得不到信号。于是，将后一组合的输出作为参考输入，根据抵消器的原理，就达到抑制干扰的目的。 图 3.4.6 自适应天线阵 图3.4.7 3元阵、 4元阵线天线阵方向图 图 3.2.11 LMS算法稳态误差 可以推出失调系数为 (3.2.77) 或者 M=μNPin (3.2.78) 式中，N是滤波器的阶数，Pin是输入信号功率。上式说明μ和输入功率加大都会增加失调系数。在保证收敛的情况下μ加大，会提高收敛速度，也说明为了减小失调系数，应该适当选择收敛速度，以保证收敛速度和失调系数都满足要求。 图3.2.12是一个LMS自适应滤波器的计算机仿真结果，阶数N=5，其输入是信号加白噪声，输入信号功率为1，中心频率是0.03fs(fs为采样频率)，噪声功率为0.5, 输入信号自相关函数的特征值为： 5.14、0.853、0.502、0.500、0.500，权系数初始值取0，μ=0.0065。图中画出了一条样本学习曲线和150条样本学习曲线的平均曲线。该