自适应滤波器
自适应滤波器是能够根据输入信号自动调整性能进行数字信号处理的数字滤波器。作为对比,非自适应滤波器有静态的滤波器系数,这些静态系数一起组成传递函数。
对于一些应用来说,由于事先并不知道所需要进行操作的参数,例如一些噪声信号的特性,所以要求使用自适应的系数进行处理。在这种情况下,通常使用自适应滤波器,自适应滤波器使用反馈来调整滤波器系数以及频率响应。
总的来说,自适应的过程涉及到将代价函数用于确定如何更改滤波器系数从而减小下一次迭代过程成本的算法。价值函数是滤波器最佳性能的判断准则,比如减小输入信号中的噪声成分的能力。
随着数字信号处理器性能的增强,自适应滤波器的应用越来越常见,时至今日它们已经广泛地用于手机以及其它通信设备、数码录像机和数码照相机以及医疗监测设备中。
结构框图
编辑下面图示的框图是最小均方滤波器(LMS)和递归最小平方(RLS)这些特殊自适应滤波器实现的基础。框图的理论基础是可变滤波器能够得到所要信号的估计。
在开始讨论结构框图之前,我们做以下假设:
- 输入信号是所要信号 和干扰噪声 之和
- 可变滤波器有有限脉冲响应结构,这样结构的脉冲响应等于滤波器系数。 阶滤波器的系数定义为
- .
- 误差信号或者叫作代价函数,是所要信号与估计信号之差
可变滤波器通过将输入信号与脉冲响应作卷积估计所要信号,用向量表示为
其中
是输入信号向量。另外,可变滤波器每次都会马上改变滤波器系数
其中 是滤波器系数的校正因子。自适应算法根据输入信号与误差信号生成这个校正因子,LMS 和 RLS 是两种不同的系数更新算法。
例子
编辑假设医院正在监测一个患者的心脏跳动,即心电图,这个信号受到 50 Hz (许多国家供电所用频率)噪声的干扰
剔除这个噪声的方法之一就是使用 50Hz 的点阻滤波器对信号进行滤波。但是,由于医院的电力供应会有少许波动,所以我们假设真正的电力供应可能会在 47Hz 到 53Hz 之间波动。为了剔除 47 到 53Hz 之间的频率的静态滤波器将会大幅度地降低心电图的质量,这是因为在这个阻带之内很有可能就有心脏跳动的频率分量。
为了避免这种可能的信息丢失,可以使用自适应滤波器。自适应滤波器将患者的信号与电力供应信号直接作为输入信号,动态地跟踪噪声波动的频率。这样的自适应滤波器通常阻带宽度更小,这就意味着这种情况下用于医疗诊断的输出信号就更加准确。
自适应滤波器的应用
编辑伪代码
编辑以下伪代码利用二阶自适应横向滤波器对一个加上噪声的正弦讯号进行滤波,可以用 Python 或 Matlab 实现。
t ← (0 TO 1 BY 0.0001) # 生成sin訊號 s ← sin(t) size_s = LEN(s) # 生成雜訊 n ← randn(size_s) # 將雜訊疊加至原sin訊號上 x ← s + n w ← [0, 0.5] u ← 0.00026 for i (1 TO 9999): y(i + 1) ← n(i: i+1) * w' e(i + 1) ← x(i + 1) - y(i + 1) w ← w + 2 * u * e(i + 1) * n(i: i + 1) ENDFOR plt.subplot(3, 1, 1) plt.plot(t, x) plt.title("sin + noise") plt.subplot(3, 1, 1) plt.plot(t, x) plt.title("sin") plt.subplot(3, 1, 1) plt.plot(t, x) plt.title("Filtered result") plt.show()
滤波器实现
编辑参考文献
编辑- Monson H. Hayes Statistical Digital Signal Processing and Modeling, Wiley, 1996, ISBN 0-471-59431-8
- Simon Haykin Adaptive Filter Theory, Prentice Hall, 2002, ISBN 0-13-048434-2