[]信号处理实验.doc

生物医学信号处理实验报告 实验一 维纳滤波器的计算机实现 一、设计目的 1．利用计算机编程实现加性干扰信号的维纳滤波。 2．将计算机模拟实验结果与理论分析结果相比较，分析影响维纳滤波效果的各种因素，从而加深对维纳滤波的理解。 二、设计原理与方法 维纳滤波是一种从噪声背景中提取信号的最佳线性方法，假定一个随机信号x(n)具有以下形式： （1-1） 其中，s(n)为有用信号，v(n)为噪声干扰，将其输入一个单位脉冲响应为h(n)的线性系统，其输出为 （1-2） 我们希望x(n)通过这个系统后得到的y(n)尽可能接近于s(n)，因此，称y(n)为信号s(n)的估值。按照最小均方误差准则，h(n)应满足下面的正则方程： （1-3） 这就是著名的维纳－霍夫方程，其中 (1-4) 是x(n)与s(n)的互相关函数，定义为 (1-5) 这里，E[·]表示求数学期望，*表示取共轭。 在要求h(n)满足因果性的条件下，求解维纳-霍夫方程是一个典型的难题。虽然目前有几种求解h(n)的解析方法，但它们在计算机上实现起来非常困难。因此本实验中利用近似方法，即最佳FIR维纳滤波方法，在计算机上实现随机信号的维纳滤波。 设h(n)为一因果序列，其长度为N，则 (1-6) 同样利用最小均方误差准则，h(n)满足下面正则方程： (1-7) 其中 (1-8) (1-9) (1-10) 这里T表示转置运算。称为信号x(n)的N阶自相关矩阵，为x(n)与s(n)的互相关函数向量。当为满秩矩阵时，由公式(1-7)可得 (1-11) 由此可见，利用有限长的h(n)实现维纳滤波器，只要已知和，就可以按上式解得满足因果性的h。虽然它不同于真正的维纳滤波器，但是只要N选择的足够大，它就可以很好地逼近真正的维纳滤波器，这一点我们可以在下面实验中得到证实。 在本实验中，s(n)由下式来确定： (1-12) 其中, 是零均值方差为的均匀分布白噪声，v(n)是与s(n)互不相关的均匀分布白噪声，其均值为零，方差。 根据理论推导，此时维纳最佳滤波器为 (1-13) 单位脉冲响应为 (1-14) 由此可以实现对信号x(n)的最佳过滤，即 (1-15) 其中为s(n)的最佳估值。同时可以推出，经过理想维纳滤波后，均方误差应为 (1-16) 在实验中，我们利用下面公式来统计均方误差： (1-17) 其中L为维纳滤波数据长度。 实际中，一般很难确切地知道和，通常是利用有限个x(n)和s(n)的样本来估计它们 (1-18) (1-19) 三、设计内容与步骤 1． 仔细阅读维纳滤波原理，根据图1.1给出的框图编制维纳滤波程序。 维纳滤波程序使用实验指导书中的程序。 2． 运行维纳滤波程序，选择L=500,N