维纳滤波器计算机实现.docxVIP

维纳滤波器的计算机实现 一、 实验目的 1、 利用计算机编程实现加性噪声信号的维纳滤波。 2、 将计算机模拟实验结果与理论分析结果相比较，分析影响维纳 滤波效果的各种因素，从而加深对维纳滤波的理解。 3、 利用维纳滤波一步纯预测方法实现对信号生成模型的参数估 计。 二、 实验原理 1、 维纳滤波器是一种从噪声中提取信号的最佳线性估计方法，假 定一个随机信号形式为： x(n)=s(n)+v(n), 其中 s(n)为有用信号， v(n)为噪声信号。而维纳滤波的作用就是让 x(n)通过一个系统 h(n)尽可能滤掉噪声，提取近似 s(n),h(n)的选择以最小均方误差 为准则。由维纳 -霍夫方程知，只要求出 φxx 及φxs 就可求出h（h=φ-1 xxφxs ）。但要求 h(n) 满足因果性要求，维纳 - 霍夫方 程便是一个难题， 这里利用最佳 FIR 维纳滤波方法求解 h(n) 的 近似，这也便于在计算机上实现，公式为： h =R-1 xx r xs 。实验中 s(n) 由 信 号 生 成 模 型 ： s(n)=as(n-1)+w(n) 确 定 ， 其 中 2 a=0.95,w(n) 是均值为 0，方差为 бw =1 的高斯白噪声， v(n) 为均值为 0，方差为 1 的高斯白噪声，且 s(n) 与 v(n) 不相关。实验中 s(n) 是已知的，但实际中如果 s(n) 已知，维纳滤波也就失去意义了，因此实验纯粹是为了理解维纳滤波原理而设计。 2、 维纳一步纯预测问题 S(n) 的生成模型：s(n)+a 1(n-1)+ ? +aps(n-p)=w(n), 已知 φxx（n）,利用 Yule-walker 方程即可得到信号生成模型参数 ai (i=1,2 ?p) 和 2 б w 。 三、 实验步骤及结果分析 1、 根据维纳滤波原理绘制程序流程图 开始 输入样本个数 L,FIR 滤波器阶数 N 产生 L个 v(n),w(n),s(n) 和 x(n) ， 利用 L 个 s(n)和 x(n) ，估计 RSS和 r xs N  检验产生序列 x(n 的自相 关和互相关函数 Y 在同一坐标内绘出 x(n) 自相关函数的 理论值和实际值 在同一坐标内绘出最后 100个 s(n) 和 调矩阵求逆子程序计算 , 将N个理想的 h(n) 和 估计的 h(n) 绘于同一坐标内 进行理想的维纳滤波得 L 个 SI (n),和最 后 100 个 s(n)绘制于同一坐标 对 x(n) 进行过滤得 L个 SR(n) ，和最后 100 个 s(n)和绘于同一坐标内 2、 根据流程图编写程序（见附录 2, 2 1）并分析运行结果： L 个 x(n),s(n), S I (n), 2 ,e SR(n), 统计 ex I ,e R L=5000,N=10 察并 、分析 果。 结束 1） 与 s(n)相比，信号 x(n)在维纳滤波前后效果比较： 最 后 100个 sn（黑色 ） 和由维纳 滤波 器 得到 的 sn （红色）的比较 最 后100个 s(n)（黑色）和 x(n)（ 红色 ） 2 w 5 1.5 4 1 3 0.5 2 0 1 -0.5 0 -1 -1 -1.5 -2 -2 -3 -2.5 10 20 30 40 50 60 70 80 90 100 -4 10 20 30 40 50 60 70 80 90 100 0 0 图1 图2 图 1 为维纳滤波后的 s(n)与最后 100 个 s(n)比较图 图 2 为未经维纳滤波的 x(n)与最后 100 个 s(n)比较图。分析：显然与 s(n)相比，x(n)在维纳滤波前与 s(n)相差很大，维纳滤波后较接近 s(n),可见滤波效果比较好。 2） 估计 (n)与理想 h(n)的比较： h(n) 估 计 值 （ 红 色 ） ， 与 真 值 （ 黑 色 ） 的 比 较 0.25 0.2 0.15 0.1 0.05 0 1 2 3 4 5 6 7 8 9 10 图 3 图 3 为估计 (n)与理想 h(n)的对比图。 分析：由图可见，二者近似程度除最后几个点外，其他近似度还是满高的，总体而言，近似效果不错。 3） 理想的维纳滤波与 FIR 维纳滤波效果对比： 最后100个sn（黑色）和由维纳滤波器得到的snw（红色）的比较 最后100个sn（黑色）和由 FIR滤波器得到的snw （红色）的比较 2 2 1.5 1.5 1 1 0.5 0.5 0 0 -0.5 -0.5 -1 -1 -1.5 -1.5 -2 -2 -2.5 10 20 30 40 50 60 70 80 90 100 -2.5 10 20 30 40 50 60 70 80 90 100 0 0 图4 图5 图 4 为理想维纳滤波效