維纳滤波原理及其matlab实现(.docVIP

维纳滤波 滤波技术是信号分析、处理技术的重要分支，无论是信号的获取、传输，还是信号的处理和交换都离不开滤波技术，它对信号安全可靠和有效灵活地传递是至关重要的。信号分析检测与处理的一个十分重要的内容就是从噪声中提取信号，实现这种功能的有效手段之一是设计一种具有最佳线性过滤特性的滤波器，当伴有噪声的信号通过这种滤波器的时候，它可以将信号尽可能精确地重现或对信号做出尽可能精确的估计，而对所伴随噪声进行最大限度地抑制。维纳滤波器就是这种滤波器的典型代表之一。 1.维纳滤波概述 维纳（Wiener）是用来解决从噪声中提取信号的一种过滤（或滤波）方法。这种线性滤波问题，可以看做是一种估计问题或一种线性估计问题。 一个线性系统，如果它的单位样本响应为，当输入一个随机信号，且 （1） 其中表示信号，)表示噪声，则输出为 （2） 我们希望通过线性系统后得到的尽量接近于，因此称为的估计值，用表示，即 （3） 则维纳滤波器的输入—输出关系可用下面图1表示。 图1 实际上，式（2）所示的卷积形式可以理解为从当前和过去的观察值，，…，…来估计信号的当前值。因此，用进行过滤问题实际上是一种统计估计问题。 　　一般地，从当前的和过去的观察值，，…估计当前的信号值成为过滤或滤波；从过去的观察值，估计当前的或者将来的信号值称为外推或预测；从过去的观察值，估计过去的信号值称为平滑或内插。因此维纳滤波器又常常被称为最佳线性过滤与预测或线性最优估计。这里所谓的最佳与最优是以最小均方误差为准则的。 　　如果我们分别以与表示信号的真实值与估计值，而用表示他们之间的误差，即 （4） 显然可能是正值，也可能是负值，并且它是一个随机变量。因此，用它的均方误差来表达误差是合理的，所谓均方误差最小即它的平方的统计期望最小： （5） 采用最小均方误差准则作为最佳过滤准则的原因还在于它的理论分析比较简单，不要求对概率的描述。 2.维纳-霍夫方程的求解 　　为了按（5）式所示的最小均方误差准则来确定维纳滤波器的冲激响应，令对的导数等于零，即可得 （6） 式中，是与的互相关函数，是的自相关函数，分别定义为 式（6）称为维纳滤波器的标准方程或维纳-霍夫（Wiener-Hopf）方程。如果已知和，那么解此方程即可求的维纳滤波器的冲激响应。 　　式（6）所示标准方程右端的求和范围即的取值范围没有具体标明，实际上有三种情况： 有限冲激响应（FIR）维纳滤波器，从到取得有限个整数值； 非因果无限冲激响应（非因果IIR）维纳滤波器，从到取所有整数值； 因果无限冲激响应（因果IIR）维纳滤波器，从到取正整数值。 上述三种情况下标准方程的解法不同，本文只描述FIR维纳滤波器的求解。 　设滤波器冲激响应序列的长度为，冲激响应矢量为 (7) 滤波器输入数据矢量为 （8） 则滤波器的输出为 （9） 这样，式（6）所示的维纳-霍夫方程可写成 或 （10） 其中 （11） 是与的互相关函数，它是一个维列矢量；是的自相关函数，是阶方阵 （12） 利用求逆矩阵的方法直接求解式（10），得 （13） 这里表示“最佳”，这就是FIR维纳滤波器的冲激响应。