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

实验一 维纳滤波器的计算机实现 一、实验目的 1 ．利用计算机编程实现加性噪声信号的维纳滤波。 2 ．将计算机模拟实验结果与理论分析结果相比较，分析影响维纳滤波效果的各种因素，从 而加深对维纳滤波的理解。 3 ．利用维纳一步纯预测方法实现对信号生成模型的参数估计。 二 实验原理与方法 维纳滤波是一种从噪声背景中提取信号的最佳线性滤波方法，假定一个随机信号 x(n) 具有以下形式： x (n) s(n) +v(n) (1-1) 其中，s(n)为有用信号，v(n)为噪声干扰，将其输入一个单位脉冲响应为 h(n) 的线性系统， 其输出为 ∞ y (n) ∑h(m)x (n −m) (1-2) m −∞ 我们希望 x(n)通过这个系统后得到的 y(n)尽可能接近于 s(n)，因此，称 y(n)为信号 s(n) 的估值。按照最小均方误差准则，h(n)应满足下面的正则方程： ∞ φ ( k ) ∑h ( m )φ ( k − m ) (1-3) xs xx m −∞ 这就是著名的维纳－霍夫方程，其中φ (m) 是 x(n)的自相关函数，定义为 xx φ (m) E[x(n)x(n +m)] (1-4) xx φ (m) 是 x(n)与 s(n)的互相关函数，定义为 xs φ (m) E[x(n)s(n +m)] (1-5) xs 这里，E[·]表示求数学期望。 在要求 h(n)满足因果性的条件下，求解维纳-霍夫方程是一个典型的难题。虽然目前 有几种求解 h(n)的解析方法，但它们在计算机上实现起来非常困难。因此，本实验中， 利用近似方法，即最佳 FIR 维纳滤波方法，在计算机上实现随机信号的维纳滤波。 设h(n)为一因果序列，其长度为 N，则 N −1 y (n) ∑h(m)x (n −m) (1-6) m 0 同样利用最小均方误差准则，h(n)满足下面方程： R h r (1-7) xx xs 其中 h [h(0), h(1),..., h(N −1)]T (1-8) φ φ φ N