数字信号处理实验报告一西交大殷.doc.docxVIP

数字信号处理Ⅱ实验报告

11896

1

1896

实验题目：维纳滤波器的计算机实现

姓名：

学号：

班级：

专业：

2

数字信号处理Ⅱ实验报告

一、实验目的

1.利用计算机编程实现加性噪声信号的维纳滤波。

2.将计算机模拟实验结果与理论分析结果相比较，分析影响维纳滤波效果的各种因素，从而加深对维纳滤波的理解。

3.利用维纳一步纯预测方法实现对信号生成模型的参数估计。

二、实验原理

维纳滤波是一种从噪声背景中提取信号的最佳线性滤波方法，假定一个随机信号x(n)具有以下形式：

x(n)=s(n)+v(n)1-1

其中，s(n)为有用信号，v(n)为噪声干扰，将其输入一个单位脉冲响应为h(n)的线性系统，其输出为

1-2

我们希望x(n)通过这个系统后得到的y(n)尽可能接近于s(n),因此，称y(n)为信号s(n)的估值。按照最小均方误差准则，b(n)应满足下面的正则方程：

1-3

这就是著名的维纳一霍夫方程，其中是φ(m)是x(n)的自相关函数，中(m)是x(n)和s(n)是的互相关函数。

在要求h(n)满足因果性的条件下，求解维纳-霍夫方程是一个典型的难题。虽然目前有几种求解h(n)的解析方法，但它们在计算机

数字信号处理Ⅱ实验报告

3

上实现起来非常困难。因此，本实验中，利用近似方法，即最佳

FIR维纳滤波方法，在计算机上实现随机信号的维纳滤波。设h(n)为一因果序列，其长度为N,则

同样利用最小均方误差准则，h(n)满足下面方程：

Rh=r

其中h=[h(O),h(1),..,h(N-D]

r=[?(0)…φ(N-1]

1-4

1-5

当R为满秩矩阵时，

1-6

由此可见，利用有限长的h(n)实现维纳滤波器，只要已知R和

rxs,就可以按上式解得满足因果性的h。只要N选择的足够大，它就可以很好地逼近理想无限长的维纳滤波器。这一点我们可以在下面实验中得到证实。

考虑维纳一步纯预测问题，假定s(n)的生成模型为

s(n)+a?s(n-1)+…+aps(n-p)=w(n)1-7

其中w(n)是均值为零，方差等于的高斯白噪声。在已知准确自

数字信号处理Ⅱ实验报告

4

相关函数φ,(n)的情况下，由下面Yule-Walker方程可以得到信号

生成模型参数a,(1,2,…,p)和9

1-8

其中R为(p+1)×(p+1)的自相关矩阵，A为(p+1)×1的系数列向量，定义为

其中e为(p+1)×1的单位列向量，除第一个元素等于1外，其余元素均为零，即

ε=[1,0,…,0]

三、实验内容和步骤：

3.1、实验流程图

5

数字信号处理Ⅱ实验报告

开始

输入样本个数L,FIR滤波器阶数N

产生L个v(n),w(n),s(n)和x(n),

利用L个s(n)和x(n),估计Rss和rxs

N

检验产生序列x(n的自相关和互相关函数

Y

在同一坐标内绘出x(n)自相关函数的理论值和实际值

在同一坐标内绘出最后100个s(n)和

调矩阵求逆子程序计算，将N个理想的h(n)和估计的h(n)绘于同一坐标内

进行理想的维纳滤波得L个S?(n),和最后100个s(n)绘制于同一坐标

对x(n)进行过滤得L个Sg(n),和最后100个s(n)和绘于同一坐标内

L个x(n),s(n),S?(n),Sg(n),统计

结束

数字信号处理Ⅱ实验报告

6

3.2、运行维纳滤波器程序(见附录程序1),N=10,L=5000,观察并记录实验结果：

1)与s(n)相比，信号x(n)在维纳滤波前后效果比较：

图1未经维纳滤波的x(n)与最后100个s(n)比较图

数字信号处理Ⅱ实验报告

7

图2维纳滤波后的s(n)与最后100个s(n)比较图

分析：显然与s(n)相比，x(n)在维纳滤波前与s(n)相差很大，维纳滤波后较接近s(n),可见滤波效果比较好。

2)估计h(n)与理想h(n)的比较：

图3为估计h(n)与理想h(n)的对比图

分析：由图可见，二者近似程度除最后几个点外，其他近似度还是满高的，总体而言，近似效果不错。

3)理想的维纳滤波与FIR维纳滤波效果对比：

数字信号处理Ⅱ实验报告

8

图4理想维纳滤波效果

图5FIR维纳滤波效果

分析：直接从图形观察，