一维纳滤波器的计算机实现

《一维纳滤波器的计算机实现》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、实验一维纳滤波器的计算机实现专业：信息与通信工程班级：电信2119学号：3112313051姓名：胡道满指导老师：殷勤业2012年11月27日实验一维纳滤波器的计算机实现一、实验目的1．利用计算机编程实现加性噪声信号的维纳滤波。2．将计算机模拟实验结果与理论分析结果相比较，分析影响维纳滤波效果的各种因素，从而加深对维纳滤波的理解。3．利用维纳一步纯预测方法实现对信号生成模型的参数估计。二、实验原理1．维纳滤波维纳滤波是一种从噪声背景中提取信号的最佳线性滤波方法，假定随机信号表示为x(n)=s(n)+v(n)（1-1）其中s(n)是有用信号

2、，v(n)是加入的噪声，将其输入一个单位脉冲响应为h(n)的线性系统，其输出为（1-2）我们希望x(n)通过这个系统后得到的y(n)尽可能接近于s(n)，称y(n)为信号s(n)的估值。按照最小均方误差准则，h(n)应满足下面的正则方程：（1-3）这就是维纳－霍夫方程其中是x(n)的自相关函数，是x(n)与y(n)的互相关函数。在实验中，利用近似方法，即最佳FIR维纳滤波方法，在计算机上实现随机信号的维纳滤波。设h(n)为一因果序列，其长度为N，则（1-4）利用最小均方误差准则，h(n)满足下面方程：(1-5)其中(1-6)（1-7）（1

3、-8）称为信号x(n)的N阶自相关矩阵，为x(n)与s(n)的互相关函数向量。当为满秩矩阵时，由公式(1-5)可得（1-9）设信号的生成模型为（1-10）其中,w(n)是零均值方差为的高斯白噪声，v(n)是与s(n)互不相关的高斯白噪声，其均值为零，方差为。根据理论推导，此时维纳最佳滤波器为(1-11)其可以实现对信号x(n)的最佳过滤，即（1-12）经过理想维纳滤波后，均方误差为（1-13）在试验中，我们利用下面公式来统计均方误差：（1-14）其中L为维纳滤波数据长度。在实际中，我们用有限个x(n)和s(n)的样本来估计x(n)的自相关

4、函数和s(n)与x(n)的互相关函数：（1-15）（1-16）我们可以采用下式检验实际中某次产生序列的自相关特性与理论值的近似程度：（1-17）（1-18）该式表示了自相关函数的理论值与某次实现的实际值的相对平方误差。实验中为了得到与自相关特性理论值相符的观测序列，往往需要多次产生序列，直到两者的相对平方误差ρ足够小。本实验中，我们取K=50，并认为且的序列才是满足要求的。2．维纳一步纯预测假定s(n)的生成模型为（1-19）在已知准确自相关函数的情况下，由下面Yule-Walker方程可以得到信号生成模型参数和（1-20）其中为的系数列

5、向量，定义为（1-21）为的单位列向量，除第一个元素等于1外，其余元素均为零，即（1-22）三、实验步骤维纳滤波结果及分析1、L=5000，N=10a)最后100个x(n)和s(n)见图1图1最后100个sn（红）和xn（蓝）由图1可以看出，滤波后s(n)基本保留了滤波前x(n)的所有信息，但是幅度明显变小，并且有些点处的方向变为反方向，这是因为噪声干扰被滤除，随机信号能量变小的原因。b)估计出的和理想的比较图2估计FIR滤波器系数（红）和理想滤波器系数（蓝）图2是估计FIR滤波器系数和理想滤波器系数的比较，由图看出，估计的很接近理想的。

6、c)理想维纳滤波和FIR维纳滤波比较图3理想滤波最后100个sn（红）和（蓝）图4估计FIR维纳滤波最后100个s(n)（红）和（蓝）图3为理想维纳滤波的最后100个s(n)和，图4为估计FIR维纳滤波的最后100个s(n)和，可以看出，二者滤波后都能很好的保留输入信号的信息，并且都有很好的滤波效果，计算得出理想维纳滤波的均方误差为0.2482，估计FIR维纳滤波器的均方误差为0.2493，可见，理想滤波的误差更小。所以，理想维纳滤波的效果要优于FIR维纳滤波。d)若去掉判断自相关和互相关特性的步骤,理想滤波器滤波效果可能不如估计FIR滤

7、波器,原因在于产生的随机序列的自相关和互相关特性可能偏离理想滤波器,而于实际的估计FIR滤波器接近,这样滤波后的效果肯定是估计FIR滤波器的好.2、N的大小对的估计和滤波效果的影响a)L=5000,N=3图5N=3估计FIR滤波器系数（红）和理想滤波器系数（蓝）理想维纳滤波均方误差EI2=0.3422，估计FIR维纳滤波均方误差为ER2=0.2988b)L=5000,N=20图6N=20估计FIR滤波器系数（红）和理想滤波器系数（蓝）理想维纳滤波均方误差EI2=0.2283，估计FIR维纳滤波均方误差为ER2=0.2293图5为N=3时估

8、计的和理想的比较，其理想维纳滤波均方误差EI2=0.3422，估计FIR维纳滤波均方误差为ER2=0.2988；图6为N=20时估计的和理想的比较，其理想维纳滤波均方误差EI2=0.2283，