维纳滤波器.docxVIP

西安电子科技大学统计与自适应信号处理仿真 学院： 班级： 学号： 姓名： 2013年12月FIR维纳滤波器1维纳滤波原理概述维纳（Wiener）是用来解决从噪声中提取信号的一种过滤（或滤波）方法。这种线性滤波问题，可以看做是一种估计问题或一种线性估计问题。一个线性系统，如果它的单位样本响应为，当输入一个随机信号，且 （1）其中表示信号，)表示噪声，则输出为 （2）我们希望通过线性系统后得到的尽量接近于，因此称为的估计值，用表示，即 （3） 则维纳滤波器的输入—输出关系可用下面图1表示。图1 维纳滤波器的输入—输出关系实际上，式（2）所示的卷积形式可以理解为从当前和过去的观察值，，…，…来估计信号的当前值。因此，用进行过滤问题实际上是一种统计估计问题。一般地，从当前的和过去的观察值，，…估计当前的信号值成为过滤或滤波；从过去的观察值，估计当前的或者将来的信号值称为外推或预测；从过去的观察值，估计过去的信号值称为平滑或内插。因此维纳滤波器又常常被称为最佳线性过滤与预测或线性最优估计。这里所谓的最佳与最优是以最小均方误差为准则的。　如果我们分别以与表示信号的真实值与估计值，而用表示他们之间的误差，即 （4）显然可能是正值，也可能是负值，并且它是一个随机变量。因此，用它的均方误差来表达误差是合理的，所谓均方误差最小即它的平方的统计期望最小：采用最小均方误差准则作为最佳过滤准则的原因还在于它的理论分析比较简单，不要求对概率的描述。2维纳-霍夫方程的求解为了按（5）式所示的最小均方误差准则来确定维纳滤波器的冲激响应，令对的导数等于零，即可得 （6）式中，是与的互相关函数，是的自相关函数，分别定义为式（6）称为维纳滤波器的标准方程或维纳-霍夫（Wiener-Hopf）方程。如果已知和，那么解此方程即可求的维纳滤波器的冲激响应。式（6）所示标准方程右端的求和范围即的取值范围没有具体标明，实际上有三种情况：有限冲激响应（FIR）维纳滤波器，从到取得有限个整数值；非因果无限冲激响应（非因果IIR）维纳滤波器，从到取所有整数值；因果无限冲激响应（因果IIR）维纳滤波器，从到取正整数值。上述三种情况下标准方程的解法不同，本文将描述因果IIR维纳滤波器和FIR维纳滤波器的求解。3 FIR维纳滤波器的求解设滤波器冲激响应序列的长度为，冲激响应矢量为 （14）滤波器输入数据矢量为 （15）则滤波器的输出为 （16）这样，式（6）所示的维纳-霍夫方程可写成 或 （17）其中 （18）是与的互相关函数，它是一个维列矢量；是的自相关函数，是阶方阵 （19）利用求逆矩阵的方法直接求解式（10），得 （20）这里表示“最佳”，这就是FIR维纳滤波器的冲激响应。4 FIR维纳滤波器的matlab实现这是FIR维纳滤波器的原理框图：图2 FIR维纳滤波器的原理框图设计维纳滤波器的过程就是寻求在最小均方误差下滤波器的单位脉冲响应或传递函数的表达式，其实质就是解维纳－霍夫（Wiener－Hopf）方程。这是用matlab设计FIR维纳滤波器的流程图：图3 FIR维纳滤波器的流程图5 FIR维纳滤波器的matlab程序%************************************************************% wiener.m% FIR维纳滤波器的设计与仿真% ***********************************************************clear all;clc;%************************* 生成信号 *************************% 信号振幅A=1.0;% 初始的相位th0=0;N=input(请输入信号长度);vn=randn(1,N);%产生随机高斯白噪声%产生信号sn = zeros(1,N);k = 0:N-1;sn = A*sin(2*pi*k/N+th0)+2;xn=sn+vn; %信号中加入随机高斯白噪声%******************** 计算自相关函数Rxx **********************Rxx=xcorr(xn,xn); %得到混合信号的自相关函数M=input(输入滤波器阶数 );for i=1:M %得到混合信号的自相关矩阵 for j=1:M rxx(i,j)=Rxx(i-j+N); endend%******************** 计算互相关函数Rxy **********************Rxy=xcorr(xn,sn); %得到混合信号和原信号的互相关函数for i=1:M rxy(i)=Rxy(i+N-1);end