数字信号处理实验——维纳滤波器设计
最佳滤波器(维纳滤波,噪声抵消)
function [a err]=wienerhopf(x,d,p) %x 为带噪观测数据； %d 为期望信号； %p 为 FIR 维纳滤波器的阶数。 x=x(:); N=length(x); if p>=N error('Model order too large'); end X=convm(x,p+1); d(N+p)=0; a=X\d'; err=d*d'-d*X*a;
a. 编写设计 p 阶 FIR 维纳滤波器的 Matlab 程序，要求能返回滤波器系数和均方误差； b. 取 α = 0.8 ，分别计算 p = 1, 2, , 2 0 的均方误差，并比较这些误差的大小。得出什 么结论并解释原因。 c. 取 p = 10 ，画出均方误差在 α = 0.1, 0.2, , 0.9 时的值，解释图中的结果。
实验三
实验内容、步骤：
实验内容包括两个：
最佳滤波器
实验一、如果 d ( n) 是 1 阶的 AR 过程，其自相关序列为
rd (k ) = α ， 0 < α < 1 。观测
k
数据 x(n) 包含期望信号 d ( n) 和零均值的不相关白噪声，方差为
σ v2 ，
x = ( n) d ( n) + v ( n)
close all;clear;clc; %%%%%%%%%%%%% a %%%%%%%%%%%%% w=0.05*pi;alpha=0; N=500;n=0:N-1; phi=2*pi*(rand(1)-0.5); d=sin(n*w+phi); g=randn(size(n)); b=1;a=[1 -0.8]; v=filter(b,a,g)+alpha*d; x=d+g; subplot(4,1,1);plot(n,x,'b-'); hold on;plot(n,d,'r--'); %%%%%%%%%%%%% c %%%%%%%%%%%%% error=[]; for p=[2 4 6]
维纳滤波器设计
维纳滤波器设计(总10页) --本页仅作预览文档封面，使用时请删除本页--维纳滤波器设计随机信号处理学号：姓名：实验一 维纳滤波器设计1实验内容设计一个维纳滤波器：（1） 产生三组观测数据，首先根据()(1)()s n as n w n =-+产生信号()s n ，将其加噪，（信噪比分别为20,10,6dB dB dB ），得到观测数据123(),(),()x n x n x n 。
（2） 估计()i x n ，1,2,3i =的AR 模型参数。
假设信号长度为L ，AR 模型阶数为N ，分析实验结果，并讨论改变L ，N 对实验结果的影响。
2实验原理滤波目的是从被噪声污染的信号中分离出有用的信号来，最大限度地抑制噪声。
对信号进行滤波的实质就是对信号进行估计。
滤波问题就是设计一个线性滤波器，使得滤波器的输出信号()y n 是期望响应()s n 的一个估计值。
下图就是观测信号的组成和信号滤波的一般模型。
观测信号()()()x n s n v n =+ 信号滤波的一般模型维纳滤波解决从噪声中提取信号的滤波问题，并以估计的结果与真值之间的误差均方值最小作为最佳准则。
它根据()()(),1,,x n x n x n m --估计信号的当前值，它的解以系统的系统函数()H z 或单位脉冲()h n 形式给出，这种系统常称为最佳线性滤波器。
维纳滤波器设计的任务就是选择()h n ，使其输出信号()y n 与期望信号()d n 误差的均方值最小。
假设滤波系统()h n 是一个线性时不变系统，它的()h n 和输入信号都是复函数，设()()()h n a n jb n =+ 0,1,n=考虑系统的因果性，可得到滤波器的输出()()()()()0*m y n h n x n h m x n m +∞===-∑ 0,1,n=设期望信号()d n ，误差信号()e n 及其均方误差()2E e n ⎡⎤⎣⎦分别为()()()()()e n d n y n s n y n =-=-()()()()()()2220m E e n E d n y n E d n h m x n m ∞=⎡⎤⎡⎤⎡⎤=-=--⎢⎥⎣⎦⎣⎦⎢⎥⎣⎦∑要使均方误差为最小，需满足：()()20E e n h j ⎡⎤∂⎣⎦=∂ 整理得()()0E x n j e n *⎡⎤-=⎣⎦，等价于()()0E x n j e n *⎡⎤-=⎣⎦上式说明，均方误差达到最小值的充要条件使误差信号与任一进入估计的输入信号正交，这就是正交性原理。
维纳滤波器的设计及Matlab仿真实现
Wiener 滤波器的设计及Matlab 仿真实现1.实验原理在许多实际应用中，人们往往无法直接获得所需的有用信号，能够得到的是退化了或失真了的有用信号。
例如，在传输或测量信号s(n)时，由于存在信道噪声或测量噪声v(n)，接受或测量到的数据x(n)将与s(n)不同。
为了从x(n)中提取或恢复原始信号s(n)，需要设计一种滤波器，对x(n)进行滤波，使它的输出y(n)尽可能逼近s(n)，成为s(n)的最佳估计，即y(n) = )(ˆn s。
这种滤波器成为最优滤波器。
Wiener 滤波器是“理想”意义上的最优滤波器，有一个期望响应d(n)，滤波器系数的设计准则是使滤波器的输出y(n)(也常用)(ˆn d表示)是均方意义上对期望响应的最优线性估计。
Wiener 滤波器的目的是求最优滤波系数],,,,,,[,1,0,1, k o o o o w w w w w -=，从而使])(ˆ)([])([)(22n d n d E n e E n J -==最小。
通过正交性原理，导出)()(k r k i r w xd x i oi -=-∑∞-∞=， 2,1,0,1,-=k该式称为Wiener-Hopf 方程，解此方程，可得最优权系数},2,1,0,1,,{ -=i w oi 。
Wiener-Hopf 方程的矩阵形式为xd o x r w R =，解方程求得xd x o r R w 1-= 2.设计思路下面我们通过具体的例子来说明Wiener 滤波器的设计方法：考虑如下图所示的简单通信系统。
其中，产生信号S(n)所用的模型为)95.01/(1)(11-+=z z H ，激励信号为)3.0,0(~)(WGN n w 。
信号s(n)通过系统函数为)85.01/(1)(12--=z z H 的信道，并被加性噪声)1.0,0(~)(WGN n v 干扰，v(n)与w(n)不相关。
确定阶数M=2的最优FIR 滤波器，以从接收到的信号x(n) = z(n) + v(n)中尽可能恢复发送信号s(n)，并用MATLAB 进行仿真。
(完整word版)维纳滤波器设计(word文档良心出品)
1.设计要求Sequence s(n) of N=2000 points is generated by AR(1) model: s(n)=as(n-1)+w(n), in which a=0.8, w(n) is white noise sequence, the mean and variance of w(n) is 0w m =，20.36w σ=.The measurement model is x(n) =s(n) +v(n), in which white noise sequence v (n) andw (n) is not related, the mean and variance of v(n) is 0v m =,21mσ=. Requirements:(1)Design IIR causal Wiener filter , calculate the filtered sequence and mean square error;(2)Design FIR Wiener filter , calculate the filtered sequence and mean square error;(3)Display raw data , noise data and filtered data on the same graph , compare the mean square error between the two cases and draw a conclusion.2.设计原理2.1维纳滤波原理概述维纳（Wiener ）是用来解决从噪声中提取信号的一种过滤（或滤波）方法。
这种线性滤波问题，可以看做是一种估计问题或一种线性估计问题。
一个线性系统，如果它的单位样本响应为)(n h ，当输入一个随机信号)(n x ，且)()()(n v n s n x += （1） 其中)(n x 表示信号，)(n v )表示噪声，则输出)(n y 为∑-=mm n x m h n y )()()( （2）我们希望)(n x 通过线性系统)(n h 后得到的)(n y 尽量接近于)(n s ，因此称)(n y 为)(n s 的估计值，用^)(n s 表示，即^)()(n s n y = （3） 则维纳滤波器的输入—输出关系可用下面图1表示。
维纳滤波毕业设计
维纳滤波毕业设计维纳滤波毕业设计维纳滤波是一种常用于信号处理领域的滤波方法，其主要目的是通过对信号进行数学处理，去除其中的噪声成分，从而提取出有用的信息。
在我的毕业设计中，我选择了维纳滤波作为研究对象，旨在探究其在图像处理中的应用。
首先，我对维纳滤波的原理进行了深入的学习和理解。
维纳滤波是一种最小均方误差滤波器，其基本思想是通过最小化信号与噪声之间的均方误差，来实现对信号的优化处理。
在图像处理中，维纳滤波可以通过对图像进行频域变换，将信号和噪声分离，然后对信号进行加权平均，从而去除噪声的影响。
接着，我开始了实验部分的工作。
我选取了一些常见的图像，包括自然风景、人物肖像等，作为实验对象。
首先，我使用了一些常见的图像处理软件，如MATLAB、OpenCV等，对原始图像进行了预处理，包括去噪、平滑等操作，以便于后续的维纳滤波处理。
然后，我使用了维纳滤波算法对预处理后的图像进行了处理。
在实验过程中，我采用了不同的参数设置，如滤波器大小、信噪比等，以比较不同参数对滤波效果的影响。
通过对比实验结果，我发现在一定范围内，滤波器大小和信噪比对维纳滤波效果有一定的影响，但并不是绝对的，需要根据具体情况进行调整。
在实验过程中，我还发现维纳滤波在一些特定情况下可能会出现一些问题。
例如，在图像中存在边缘、纹理等细节信息时，维纳滤波可能会导致图像模糊，失去一些细节。
为了解决这个问题，我尝试了一些改进的方法，如结合边缘检测算法、纹理增强算法等，以提高滤波效果。
最后，我对实验结果进行了分析和总结。
通过对比不同图像的处理效果，我发现维纳滤波在去除噪声方面具有一定的优势，可以有效地提取出图像中的有用信息。
然而，在处理一些特殊情况下的图像时，维纳滤波的效果可能会受到一些限制，需要结合其他算法进行改进。
综上所述，我的毕业设计主题是维纳滤波的应用研究。
通过对维纳滤波原理的学习和实验的开展，我对维纳滤波的工作原理和应用场景有了更深入的了解。
维纳滤波复原实验报告
维纳滤波复原实验报告一、实验介绍维纳滤波是一种常用的图像复原技术，可以通过提供滤波器来降低图像的噪声和估计原始图像。
本次实验旨在通过使用维纳滤波器来复原被噪声污染的图像。
二、实验方法1. 实验准备首先需要准备一个带有噪声的图像作为输入图像，以及一个用作参考的干净图像。
通过加载图像，可以将两幅图像转换为灰度图像来简化处理。
2. 维纳滤波器的建立维纳滤波器可以通过以下公式来构建：H(u, v) = \frac{1}{H(u, v)} \cdot \frac{{ F(u, v) ^2}}{{ F(u, v) ^2 + S_n(u, v)}} 其中，H(u, v)是滤波器的频域函数，F(u, v)是输入图像的傅里叶变换，S_n(u, v)是噪声功率谱。
通过计算输入图像的傅里叶变换，以及噪声功率谱，可以根据上述公式来生成维纳滤波器。
3. 图像复原将输入图像通过傅里叶变换转换到频域，然后与维纳滤波器相乘，最后再进行傅里叶反变换，即可得到复原后的图像。
三、结果与讨论在实验中，我们使用了一幅被高斯噪声污染的图像作为输入图像，并使用了一个无噪声的参考图像。
通过对输入图像进行傅里叶变换，我们得到了输入图像的频域表示。
接着，根据输入图像和参考图像的功率谱，我们生成了对应的维纳滤波器。
最后，我们将输入图像通过傅里叶变换转换到频域，然后与维纳滤波器相乘，再进行傅里叶反变换，得到了复原后的图像。
实验结果显示，通过应用维纳滤波器，最终得到的复原图像与参考图像相比较为接近，且噪声得到了明显的减少。
这证明了维纳滤波的有效性和可行性。
然而，维纳滤波也存在一些限制。
由于维纳滤波是一种线性滤波方法，当输入图像中存在较大的模糊或失真时，滤波器可能无法恢复出清晰的图像。
此外，既有的维纳滤波器还无法处理复杂的噪声类型，如椒盐噪声或周期性噪声。
四、实验总结本次实验通过使用维纳滤波器来复原被噪声污染的图像，展示了维纳滤波的效果和限制。
维纳滤波是一种常用的图像复原技术，能够有效地降低图像噪声并估计原始图像。
中科院数字信号处理作业维纳滤波
中科院数字信号处理作业维纳滤波数字信号处理实验⼀⼀、实验⽬的1、掌握FIR 维纳滤波器的基本原理2、设计FIR 维纳滤波器，对加噪声的信号进⾏滤波处理；⼆、实验原理设期望信号为()S n ，噪声信号为()V n ，则观测信号为()()()X n S n V n =+，为了设计滤波器的最佳传递函数，要先求出期望信号与观测数据的互相关函数xs R 和观测数据的⾃相关函数xx R ，然后利⽤公式-1opt h xs xs R R =求出最佳传递函数，再对信号进⾏()X n 进⾏滤波处理。
三、实验内容及要求1、在直⾓坐标系产⽣⼀个围绕原点做半径为1的圆；2、产⽣均值为0，⽅差分别为0.05和0.06的⾼斯噪声500个；3、将两种杂声分别加在圆的x 、y ⽅向；4、设计FIR 滤波器，根据要求的误差范围不断改变滤波器的阶数，直⾄满⾜误差允许值；5、分别绘制出x ⽅向和y ⽅向的期望信号、噪声信号、观测信号、滤波后的信号、误差信号的曲线图；6、在同⼀幅图中绘制出期望信号、观测信号和滤波后点⽬标的运动轨迹四、实验结果及分析1、在滤波器阶数分别为100和400时，X 轴⽅向的期望信号、噪声信号、观测信号、滤波后信号和误差信号分别如图1和图2所⽰。
图1 X 轴⽅向各个信号滤波器阶数为100图2 X轴⽅向各个信号滤波器阶数为4002、在滤波器阶数分别为100和400时，Y轴防线的期望信号、噪声信号、观测信号、滤波后信号和误差信号如图2所⽰。
图3 Y轴⽅向各个信号滤波器阶数为100图4 Y轴⽅向各个信号滤波器阶数为4003、在滤波器阶数分别为100和400时，期望信号、观测信号和滤波信号在⼀张图中显⽰的结果如图5和图6所⽰，其中红⾊的线为期望信号、绿⾊的线为加噪后的信号、蓝⾊的线为滤波后的信号。
图5 滤波器阶数为100 图6 滤波器阶数为4004、滤波器的阶数与误差的关系表如表1所⽰表1 滤波器阶数与误差对照表滤波器阶数 50 100 200 300 400 x ⽅向误差 0.028512 0.014955 0.012172 0.008156 0.007988 y ⽅向误差 0.0291100.018449 0.005544 0.006908 0.006098 总误差 0.040747 0.023749 0.013375 0.010688 0.010050由实验结果可以看出，滤波器的阶数越⾼，滤波的效果越好。
(整理)维纳滤波器的设计
一、实验题目维纳滤波器的设计 (结维纳霍夫方程，最小的均方误差)维纳-霍夫方程矩阵形式已知期望信号与观测数据的互相关函数及观测数据的自相关函数二、实验要求：设计一维纳滤波器。
(1)产生三组观测数据，首先根据()(1)()s n as n w n =-+产生信号()s n ，将其加噪(信噪比分别为20,10,6dB dB dB )，得到观测数据123(),(),()x n x n x n 。
(2)估计()i x n ，1,2,3i =的AR 模型参数。
X （n ）=w(n)+…+bw(n-q)假设信号长度为L ，AR 模型阶数为N ，分析实验结果，并讨论改变L ，N 对实验结果的影响。
三、实验原理：维纳滤波器设计的任务就是选择()h n ，使其输出信号()y n 与期望信号()d n 误差的均方值最小，实质是解维纳－霍夫方程。
假设滤波系统是一个线性时不变系统，它的()h n 和输入信号都是复函数， 设 ()()()h n a n jb n =+ n=0,1,… 考虑系统的因果性，可得到滤波器的输出()()()()()0y n h n x n h m x n m m +∞=*=-∑= n=0,1,…设期望信号d (n )，误差信号()e n 及其均方误差()2E e n ⎡⎤⎢⎥⎣⎦分别为()()()()()e n d n y n s n y n =-=-()()()()()()2220E e n E d n yn E d n h m x n m m ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎣⎦+∞=-=--∑= 要使均方误差为最小，需满足：1xx xd -=h R R()20E e n h j⎡⎤⎢⎥⎣⎦∂=∂由上式可以推导得到()()0E x n j e n *⎡⎤⎢⎥⎣⎦-= 上式说明，均方误差达到最小值的充要条件是误差信号与任一进入估计的输入信号正交，这就是正交性原理。
将()()0E x n j e n *⎡⎤⎢⎥⎣⎦-=展开，得 ()()()0dx xx m r k h m r n m +∞*=-=-∑ k=0,1,…对两边取共轭，并利用相关函数的性质()()yx xy r k r k *-=，得()()()()()0xx xd xx m k m k m r k r h r h k +∞==-=*∑ k=0,1,…此式称为维纳-霍夫(Wiener-Hopf)方程。
维纳滤波器的设计
基于运动图像复原的维纳滤波器设计通信工程 09250104 张明哲指导教师陈海燕副教授摘要景物成像过程中可能会出现模糊、失真或混入噪声，最终导致图像质量下降。
这种质量的下降会造成图像中的目标很难识别或者图像中的特征无法提取，必须对其进行恢复,维纳滤波是一种常见的图像复原方法。
本设计主要对维纳滤波的基本原理进行研究，并结合MATLAB中的函数，设计相应的维纳滤波器，对运动模糊图像和它的加噪图像进行复原。
之后，对逆滤波和维纳滤波进行图像复原仿真实验，并对比它们的复原效果。
关键词：维纳滤波；图像恢复；退化模型AbstractImaging features may appear blurred, distorted or mixed with noise in the process of scene imaging. As a consequence, quality of images is lowered,which in digital images is likely to make it difficult to identify the target image or to extract the image features, images must be restored, then. Wiener filter is a common method for image restoration.This paper mainly introduces the basic principles of Wiener filtering, and function of MA TLAB are combined to design the corresponding Wiener filter to restore motion-blur images. And noise restoration are also taken into account in the design.Finally, the inverse filter and Wiener filter for image restoration simulation experiment be taken, and compared the differences between them.Key words: Wiener filter; image restoration; degraded image一，概述（一）图像复原的定义图像复原就是研究如何从所得的变质图像中复原出真实图像，或说是研究如何从获得的信息中反演出有关真实目标的信息，如果是因为在摄像时相机和被摄景物之间有相对运动而造成的图像模糊则称为运动模糊。
数字信号处理总结+维纳声音滤波器的设计
数字信号处理总结+维纳声音滤波器的设计第一篇：数字信号处理总结+维纳声音滤波器的设计现代数字信号处理大作业维纳声音滤波器的设计摘要在数字信号中往往存在很多扰动信号，即各种噪声。
在信噪比较低时，原始信号会变得难以分辨，因此需要对输入信号进行处理，以提取有用信号，即数字信号处理。
其中主要方法是数字滤波器的设计，数字滤波器主要分为两大类，有限脉冲滤波器（FIR）和无限脉冲响应滤波器（IIR）。
本文主要介绍有限长冲击响应（FIR）维纳滤波器的设计，采用MATLAB软件设计了一个FIR维纳滤波器，用以对有噪声干扰的鸟鸣声进行降噪处理。
最后用基于MATLAB函数设计的维纳滤波器对一段鸟鸣声进行滤波处理，通过滤波前后信号的波形图的对比，分析不同信噪比下的滤波效果。
关键词：维纳滤波器；降噪；信噪比现代数字信号处理大作业目录1课程总结..................................................................................................................1 2 绪论.. (1)2 2.1 数字滤波器简介............................................................................................12 2.1.1 数字滤波器概述.....................................................................................12 2.1.2 数字滤波器的优点.................................................................................12 2.2 设计内容........................................................................................................12 2.3 研究方法.. (12)3 数字滤波器...........................................................................................................14 3.1 数字滤波器原理.. (14)3.2 数字滤波器的分类........................................................................................14 3.3 数字滤波器的实现........................................................................................15 3.4 维纳滤波器....................................................................................................15 4 维纳滤波器设计 (16)4.1 维纳滤波器的设计原理................................................................................16 4.2 维纳滤波器的仿真结果................................................................................17 5总结 (2)1 参考文献 (22)现代数字信号处理大作业1课程总结信号是携带信息的工具，是信息的载体。
信号处理实验--FIR维纳滤波器
实验报告一、 实验题目假设一个点目标在x ，y 平面上绕单位圆做圆周运动，由于外界干扰，其运动轨迹发生了偏移。
其中，x 方向的干扰为均值为0，方差为0.05的高斯噪声；y 方向干扰为均值为0，方差为0.06的高斯噪声。
1. 产生满足要求的x 方向和y 方向随机噪声500个样本；2. 明确期望信号和观测信号；3. 试设计一FIR 维纳滤波器，确定最佳传递函数：1opt xx xs h R R -=，并用该滤波器处理观测信号，得到其最佳估计。
（注：自行设定误差判定阈值，根据阈值确定滤波器的阶数或传递函数的长度）。
4.要求3中，也可以选择Kalman 滤波器进行滤波处理，采用哪种滤波器可以自由选择。
5.分别绘制出x 方向和y 方向的期望信号、噪声信号、观测信号、滤波后信号、误差信号的曲线图； 6.在同一幅图中绘制出期望信号、观测信号和滤波后点目标的运动轨迹。
7. 实验报告要求：给出求解思路和结果分析，给出MATLAB 实现源程序和程序注解。
二、 解题思路目标信号在x 、y 方向均受到高斯噪声的干扰，本实验分别从两个方向分别进行滤波处理，最终得到整体滤波处理后的信号。
具体解题思路如下：8. 选用FIR 维纳滤波器进行滤波处理9. 产生x 、y 方向噪声信号xnoise 和 ynoise ，产生x 、y 方向期望信号xs 、ys ，产生x 、y 方向观测信号x 、y 。
10. 分别求出x 、y 方向上观测信号的自相关矩阵以及观测信号与期望信号的互相关矩阵，根据公式1opt xx xs h R R -=确定最佳传递函数。
11. 根据min 2)]([n E e =()opt T d h xd *2R -σ分别求出x 、y 方向上最小均方误差。
三、 MATLAB 实现源程序close all ;clear all ;%x 轴方向和y 轴方向期望信号x0=0;y0=0;r=1;theta=0:pi/249.5:2*pi;xs=x0+r*cos(theta);ys=y0+r*sin(theta);%产生x 轴方向和y 轴方向噪声n=500;%产生500个随机样本ux=0;vx=0.05;t=randn(1,n); %产生1*n个伪随机数t=t/std(t);t=t-mean(t);xnoise=ux+sqrt(vx)* t;uy=0;vy=0.06;g=randn(1,n);g=g/std(g);g=g-mean(g);ynoise=uy+sqrt(vy)* g;% x轴方向和y轴方向观测信号x=xs+xnoise;y=ys+ynoise;%x方向上观测信号的自相关矩阵rxx=xcorr(x);for i=1:100for j=1:100mrxx(i,j)=rxx(500-i+j);endend%x方向上观测信号与期望信号的互相关矩阵rxs=xcorr(x,xs);for i=1:100mrxs(i)=rxs(499+i);endhoptx=inv(mrxx)*mrxs';%求出x方向上最佳传递函数fx=conv(x,hoptx);%维纳滤波后x方向上的输出nx=sum(abs(xs).^2);eminx=nx-mrxs*hoptx;%求出x方向上的最小均方误差%y方向上观测信号的自相关矩阵ryy=xcorr(y);for i=1:100for j=1:100mryy(i,j)=ryy(500-i+j);endend%y方向上观测信号与期望信号的互相关矩阵rys=xcorr(y,ys);for i=1:100mrys(i)=rys(499+i);endhopty=inv(mryy)*mrys'; %求出y方向上最佳传递函数fy=conv(y,hopty); %维纳滤波后y方向上的输出ny=sum(abs(ys).^2);eminy=ny-mrys*hopty; %求出y方向上的最小均方误差plot(xs,ys,'-r',x0,y0,'.');hold on;plot(x,y,'b:');hold on;plot(fx,fy,'g-');legend('期望信号','圆周运动中心','观测信号','滤波后点目标的运动轨迹') axis square;figure;subplot(2,4,1)plot(xs);title('x方向期望信号');subplot(2,4,2)plot(xnoise);title('x方向噪声信号');subplot(2,4,3)plot(x);title('x方向观测信号');subplot(2,4,4)n=0:500;plot(n,eminx);title('x方向最小均方误差');subplot(2,4,5)plot(ys);title('y方向期望信号');subplot(2,4,6)plot(ynoise);title('y方向噪声信号');subplot(2,4,7)plot(y);title('y方向观测信号);subplot(2,4,8)plot(n,eminy);title('y方向最小均方误差');四、结果及其分析运行结果如下图：分析：1.运行结果中滤波后的到的信号与原始信号和噪声信号的对比可以看出，滤波后的结果与期望信号还是很接近的，整体上达到了最优滤波的效果。
改进的维纳滤波器的设计及应用
信息产业改进的维纳滤波器的设计及应用颜军田祎(陕西商洛学院计算机科学系，陕西商洛726000)1维纳滤波的基本原理设是某平稳随机过程的一个取样序列，该随机过程的自相关函数或功率谱是已知的或能够由估计得到。
在传输或测量时，由于存在信道噪声或测量噪声，使得接收或测量到的数据与不同，如果噪声是加性的，则。
为了能从中提取或恢复原始信号，需要设计一个滤波器对进行滤波，使滤波器的输出近可能逼近，成为的最佳估计，，这种滤波就是维纳滤波。
维纳滤波是一种最佳线性滤波，它根据混有噪声的观测数据，在最小均方误差准则下，得到对信号的最佳估计。
2维纳滤波器的设计2.1维纳滤波器的实现滤波器研究的一个基本问题就是：如何设计和建立最佳或最优的滤波器。
所谓最佳滤波器是指能够根据某一最佳准则进行设计的滤波器。
假定线性滤波器的输入为有用信号和噪声之和，两者均为广义平稳过程且已知它们的二阶统计特性，根据最小均方误差准则（滤波器的输出信号与期望信号之差的均方值最小），维纳求得了最佳线性滤波器的参数。
这种滤波器称为维纳滤波器。
设观测信号是由信号和噪声干扰组成的。
在图1中，代表误差信号，是理想信号，一般指期望输出信号。
当时，此时希望从中估计出。
由正交性原理，则：即：此方程称之为维纳霍甫夫方程。
对IIR 滤波器，由维纳霍甫夫方程有：。
对上式的两边进行Z 变换，有：2.2维纳滤波器的推广假设噪声为加性的，且与纯净信号不相关，则有，设为与的互功率谱密度，而对应为的功率谱密度，则：，得到非因果维纳滤波器的滤波因子为：在实际的操作过程中，还有经常用到如下的推广形式：3改进的维纳滤波器的算法实现假设噪声是加性的，有，这里是观测信号，是纯净信号，是噪声信号。
这样，由维纳滤波器的推广可得到：(1)是纯净信号的功率谱估计，是噪声的功率谱，用来调整滤波器的增益，是滤波器的功率，此处一般取为0.5。
这样可得到：(2)(3)是观测信号的傅立叶变换，是纯净信号的傅立叶变换的估计。
维纳滤波器实验昆明理工大学
昆明理工大学信息工程与自动化学院学生实验报告（2017 —2018 学年第二学期）课程名称：生物医学信号处理开课实验室：设备编号：实验日期：一、实验目的1、掌握维纳滤波器，并能用MATLAB实现维纳滤波器。
2、通过实验学习利用维纳滤波器从噪声中估计信号，以便在工程中加以应用。
3、能用MATLA的toolbar里SPT工具处理信号。
二、实验原理纳滤波是指从连续的（或离散的）输入信号中滤除噪声和干扰以提取有用信息的过程，而相应的装置称为滤波器。
维纳滤波器（Wiener filter）是一种以最小平方为最优准则的线性滤波器，在一定的约束条件下，其输出与一给定函数的差的平方达到最小，通过数学运算最终可变为一个托不利兹方程的求解问题。
维纳滤波器有被称为最小二乘法或最小平方法滤波器，目前是最基本的滤波方法之一。
从噪声中提取信号波形的各种估计方法中，维纳滤波是一种最基本的方法，适用于需要从最小噪声中分离出有用的信号是整个信号，而不是其中的几个参量。
二、实验原理1、设计维纳滤波器对噪声语言信号进行恢复。
噪声信号如下load chirp;s=y;M=length(s);N=200;n=randn(M,1)/10;Z=s+n;计算Wiener-Hof方程的最佳滤波系数，从滤波结果（滤波后波形和信噪比改善）分析滤波器的作用。
（1）源程序及其说明：（2）程序运行结果（3）其他滤波方法对噪声语音信号的处理使用SPTool工具，下图为chirp信号经pzip滤波处理的结果：下图为chirp的信号经FIRbp滤波器处理的结果。
下图为chirp信号经pzip滤波处理的结果：下图为chirp信号经维纳滤波处理结果：分析：比较维纳滤波器对于被污染的信号起到了更好的滤波作用，pzlp滤波器得到的的波形相比于firbp滤波器得到的的波形有一定的效果，但总体都没有达到我们要求的滤波效果。
2、GUI程序和图绘制输入信号s的波形function pushbutton1_Callback(hObject, eventdata, handles)S=evalin('base','S');axes(handles.axes1);plot(S);绘制白噪声n的波形function pushbutton2_Callback(hObject, eventdata, handles)n=evalin('base','n');axes(handles.axes2);plot(n);绘制观测信号Z的波形function pushbutton3_Callback(hObject, eventdata, handles)Z=evalin('base','Z');axes(handles.axes3);plot(Z);绘制输出信号S的波形function pushbutton4_Callback(hObject, eventdata, handles)s=evalin('base','s');axes(handles.axes4);plot(s);控制输入信号Z的声音function pushbutton5_Callback(hObject, eventdata, handles) Z=evalin('base','Z');sound(Z);控制输出信号S的声音function pushbutton7_Callback(hObject, eventdata, handles) S=evalin('base','S');sound(S);GUI运行图：。
维纳滤波器的设计与分析
（2）固 （3） （4） （5）
图2 维纳滤波前)(n x 和)(n s 图3 理想维纳滤波后)(n x 和)(n s
（6） （7）
f f 图6 4=N 估计∧
)(n h 与理想)(n h 图7 4=N FIR 维纳滤波后)(n s 与(x 图8 21=N 估计∧
)(n h 与理想)(n h 图9 21=N FIR 维纳滤波后)(n s 与x 固定11=N ，分别取5000010000、=L ，观察L 的大小对∧
)(n h 的估计和滤波效果的影响。
10000时，2395.0=f e ；当50000=L 时，2465.0=f e 。
从仿真结果可知，L 越大，
)(n 越接近，精度越高，滤波效果越好。
图10 10000=L 估计∧
)(n h 与理想(n h 图11 10000=L FIR 维纳滤波后)(n s 图12 50000=L 估计∧
)(n h 与理想)(n h 图13 50000=L FIR 维纳滤波后)(n s 与x 总结
从仿真结果可以看出，样本个数越大，参数精度越高。
影响维纳滤波效果的因
素包括样本个数L 、FIR 滤波阶数，且均成正比关系，L 越大或者FIR 滤波器的阶数越大则维纳滤波的效果越好。
维纳滤波器的优点是适应面较广，无论平稳随机过程是连续的还是离散的，是标量的还是向量的，都可应用。
对某些问题，还可求出滤波器传递函数的显式解，并进而采用由简单的物理元件组成的网络构成维纳滤波器。
维纳滤波器的缺点是，要求得到半无限时间区间内的全部观察数据的条件很难满足，同时它也不能用。
数字信号处理2自适应信号滤波
实验二自适应信号滤波实验报告班级：姓名：学号：实验二 自适应信号滤波一、实验目的1．利用自适应LMS 算法实现FIR 最佳维纳滤波器。
2．观察影响自适应 LMS 算法收敛性，收敛速度以及失调量的各种因素，领会自适应信号处理方法的优缺点。
3．通过实现AR 模型参数的自适应估计，了解自适应信号处理方法的应用。
二、实验原理如果信号是由有用信号和干扰信号组成，即n n n y x w =+（2－1）利用维纳滤波方法可以从信号中得到有用信号的最佳估计。
假如最佳维纳滤波器由一个FIR 滤波器所构成，则其最佳权系数向量h 可表示为1h R r -=（2－2） 其中0,1,....,Tm h h h h ⎡⎤=⎣⎦（2－3） TR E yy ⎡⎤=⎣⎦（2－4） [()]r E x n y =（2－5）[(),(1),...,()]T y y n y n y n M =--（2－6）但是实际中，一般很难知道准确的统计量R 和r ，因此，若设计一个维纳滤波器,事先要估计出R 和r 。
同时，当R 和r 改变时（如果信号或干扰是非平稳的），需要重新计算h ，这是非常不方便的。
虽然卡尔曼滤波方法无需事先知道R 和r ，但它必须知道系统的状态方程和噪声的统计特性，这在实际中也是很难办到的。
根据卡尔曼滤波的思想，Widrow 等提出了一种自适应最小均方误差算法(LMS)，这种算法不需要事先知道相关矩阵R 和r ，当得到一个观察值，滤波器自动“学习”所需要的相关函数，从而调整FIR 滤波器的权系数，并最终使之收敛于最佳值，即维纳解。
下面是自适应FIR 维纳滤波器的LMS 算法公式：^^()()()Mm m x n h n y n m ==-∑（2－7）^()()()e n x n x n =-（2－8）^^(1)()2()(),0,...,m m h n h n e n y n m m M μ+=+⋅-=（2－9）因此，给定初始值，每得到一个样本n y ，可以递归得到一组新的滤波器权系数，只要步长满足 max10μλ<<（2－10）其中为max λ矩阵R 的最大特征值，当n →∞时，()h n 收敛于维纳解。
数字信号处理II实验一 维纳滤波 基于MATLAB
数字信号处理II实验一实验一维纳滤波器的计算机实现一、实验目的1．利用计算机编程实现加性干扰信号的维纳滤波。
2．将计算机模拟实验结果与理论分析结果相比较，分析影响维纳滤波效果的各种因素，从而加深对维纳滤波的理解。
3．利用维纳预测方法实现对 AR模型的参数估计。
二、实验原理与方法维纳滤波是一种从噪声背景中提取信号的最佳线性滤波方法，假定一个随机信号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）这就是著名的维纳－霍夫方程，其中是 x(n)的自相关函数，定义为：（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。
对于 ARMA 模型或 MA 模型来说，虽然它不同于理想的维纳滤波器，但是只要N选择的足够大，它就可以很好地逼近理想的维纳滤波器，这一点我们可以在下面实验中得到证实。
维纳滤波器
随机信号分析实验报告维纳滤波器的设计与仿真①问题描述滤波技术是信号分析、处理技术的重要分支，无论是信号的获取、传输，还是信号的处理和交换都离不开滤波技术，它对信号安全可靠和有效灵活地传递是至关重要的。
信号分析检测与处理的一个十分重要的内容就是从噪声中提取信号，实现这种功能的有效手段之一是设计一种具有最佳线性过滤特性的滤波器，当伴有噪声的信号通过这种滤波器的时候，它可以将信号尽可能精确地重现或对信号做出尽可能精确的估计，而对所伴随噪声进行最大限度地抑制。
维纳滤波器就是这种滤波器的典型代表之一。
维纳滤波器具有以下一些特点：1.假设：信号以及附加噪声都是已知频谱特性或者自相关和互相关的随机过程2.性能标准：最小均方差3.能够用标量的方法找到最优滤波器4.维纳滤波器的设计目的是就是滤除按照统计方式干扰信号的噪声。
②模型建立假设维纳滤波器的输入信号是s(t)，叠加噪声n(t)。
输出信号x(t) 通过滤波器g(τ) 使用下面的卷积运算得到：x(t) = g(τ) * (s(t) + n(t))其中s(t) 是需要估计的原始信号n(t) 是噪声x(t) 是估计出的信号（我们希望它能等同于s(t)）g(τ) 是维纳滤波器误差是e(t) = s(t + d) − x(t) ，方差是e2(t) = s2(t + d) − 2s(t + d)x(t) + x2(t) 其中s(t + d) 是所期望的滤波器输出e(t) 是误差根据d 的不同，问题名称可以更换为：如果d > 0 那么问题是预测如果d = 0 那么问题是滤波如果d < 0 那么问题是平滑将x(t) 写成卷积积分：计算平方误差的均值，可得其中Rs 是s(t) 的自相关函数Rx 是x(t) 的自相关函数Rxs 是x(t) 和s(t) 的互相关函数如果信号s(t) 和噪声n(t) 是不相关的（例如，互相关是0）那么请注意这个的目的是求最优的g(t)，使得E(e2)最小。
维纳滤波器
维纳滤波实验——噪声中语音信号的恢复班级：学号：姓名：实验内容信号从发送者传送到接受者往往受到几种形式的变形而削弱，维纳滤波是一种从接收的原始信号中恢复信号的方法。
给定信号为chiep信号(鸟叫声)，数据文件为chirp.mat，用频域法实现维纳滤波。
步骤如下：1、产生信号s(n)和观测z(n),信号为chiep信号，观测为信号叠加上高斯白噪声。
2、估计信号s(n)和z(n)的功率谱，计算维纳滤波的传递函数。
3、计算输入z(t)和输出y(t)的频谱，并对输出信号求反变换得到时域的输出信号。
用到的函数：装入数据文件，Load傅里叶变换与反变换，fftn ifftn谱估计，periodogra一．实验结果用频域法在MATLAB仿真文件为：clear all;load('chirp','y','Fs');z=awgn(y,20); %加噪声subplot(1,3,1); %画语音信号y的时域图plot(y);title('语音信号y的时域图');subplot(1,3,2) %画加噪后的观测信号图plot(z);title('加噪后的观测信号');X=fftn(y);Px1=X.*conj(X)/50; %原信号频谱密度X=fftn(z);Px2=X.*conj(X)/50; %观测信号频谱密度H=Px1\Px2; %传递函数py=Px2.*H; %输出信号的频谱yn=ifftn(py); %输出信号subplot(1,3,3)plot(yn);title('输出信号');sound(y); %原信号sound(z); %加了噪声的测试信号 sound(yn) %输出信号得到的结果图：原语音信号的时域图语音加噪后的观测信号时域图输出信号的时域图二．实验结果分析及心得体会维纳（Wiener ）滤波是一种从噪声中提取引号波形最基本的方法，适用于需要从噪声中分离出的有用信号是整个信号（波形），而不只是它的几个参量。
维纳滤波器-维纳霍夫方程课程设计
电气工程学院信号处理课程设计报告设计题目：维纳滤波器专业：生物医学工程指导教师:叶立夏学生姓名: 叶立夏学号: 144起迄日期: 2013年12月20日—2014年1月15日如果有做课程设计的同学不懂的，可以联系我！Qq：2目录前言 (3)1 设计任务及指标 (3)1.1 课程设计的内容和要求 (4)1.2 对课程设计成果的要求 (4)2 设计思想 (4)2.1 概述 (4)2.2 主要仪器设备及耗材 (5)3 课程设计的组成部分 (5)3.1具体操作 (5)4 实验分析 (9)4.1 原始图像显示 (9)4.2 噪声的强度对维纳滤波器的影响 (13)4.3 阶数对滤波效果的影响 (17)4.4 数据长度对维纳滤波的影响 (21)5 设计总结 (26)5.1思考题 (26)5.2实验心得 (26)6 主要参考文献 (27)附录 (27)前言去除信号中的噪声影响是信号处理中的一个重要内容，而滤波则是实现这一功能的重要手段之一。
滤波器可以分为两类，及经典滤波器和现代滤波器经典滤波器是假定输入信号中嘚瑟有用成分和希望去除的成分各自占有不同的频带。
当输入信号通过一个滤波器是可将欲去除的成分有效的去除，如果信号和噪声的频谱相互重叠，那么经典滤波器将无能为力。
现代滤波器理论研究的主要内容是从含有噪声的数据记录中估计出信号的某些特征灬信号本身。
一旦信号被估计出，那么估计出的信号的信噪比将比原信号的高。
现代滤波器把信号和噪声都视为随机信号，利用它们的统计特征导出一套最佳的估值算法，然后用硬件或软件予以实现。
现代滤波器理论源于维纳在20世纪40年代及其以后的工作，因此维纳滤波器便是这一类滤波器的典型代表。
维纳滤波器，也是最小平方滤波器，其基本思路为：设计一个滤波器，使其与输入信号滤波后的输出与期望输出在最小平方意义下的最佳逼近。
寻求最小均方误差的实质就是解维纳-霍夫方程。
1.设计任务及指标1.1课程设计的内容和要求（包括原始数据、技术要求、工作要求等）：本设计的目的是产生用于信号滤波的维纳霍夫方程。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验一 维纳滤波1. 实验内容设计一个维纳滤波器：（1) 产生三组观测数据,首先根据()(1)()s n as n w n =-+产生信号()s n ，将其加噪，（信噪比分别为20,10,6dB dB dB ）,得到观测数据123(),(),()x n x n x n 。
（2） 估计()i x n ，1,2,3i =的AR 模型参数。
假设信号长度为Ｌ,AR 模型阶数为N,分析实验结果,并讨论改变L,N 对实验结果的影响。
2. 实验原理滤波目的是从被噪声污染的信号中分离出有用的信号来，最大限度地抑制噪声。
对信号进行滤波的实质就是对信号进行估计。
滤波问题就是设计一个线性滤波器,使得滤波器的输出信号()y n 是期望响应()s n 的一个估计值。
下图就是观测信号的组成和信号滤波的一般模型。
观测信号()()()x n s n v n =+ 信号滤波的一般模型 维纳滤波解决从噪声中提取信号的滤波问题,并以估计的结果与真值之间的误差均方值最小作为最佳准则。
它根据()()(),1,,x n x n x n m --估计信号的当前值,它的解以系统的系统函数()H z 或单位脉冲()h n 形式给出，这种系统常称为最佳线性滤波器。
维纳滤波器设计的任务就是选择()h n ,使其输出信号()y n 与期望信号()d n 误差的均方值最小。
假设滤波系统()h n 是一个线性时不变系统，它的()h n 和输入信号都是复函数，设()()()h n a n jb n =+ 0,1,n=考虑系统的因果性,可得到滤波器的输出()()()()()0*m y n h n x n h m x n m +∞===-∑ 0,1,n=设期望信号()d n ，误差信号()e n 及其均方误差()2E e n ⎡⎤⎣⎦分别为()()()()()e n d n y n s n y n =-=-()()()()()()2220m E e n E d n y n E d n h m x n m ∞=⎡⎤⎡⎤⎡⎤=-=--⎢⎥⎣⎦⎣⎦⎢⎥⎣⎦∑ 要使均方误差为最小,需满足:()()20E e n h j ⎡⎤∂⎣⎦=∂ 整理得()()0E x n j e n *⎡⎤-=⎣⎦，等价于()()0E x n j e n *⎡⎤-=⎣⎦上式说明,均方误差达到最小值的充要条件使误差信号与任一进入估计的输入信号正交，这就是正交性原理。
将()()0E x n j e n *⎡⎤-=⎣⎦展开，得()()()()00m E x n k d n h m x m +∞***=⎡⎤⎛⎫--=⎢⎥ ⎪⎝⎭⎣⎦∑整理得 ()()()0dx xx m r k h m r m k +∞*=-=-∑ 0,1,2,k= 等价于()()()()()0dx xx xx m r k h m r k m h k r k +∞==-=*∑ 0,1,2,k=此式称为维纳-霍夫（W ｉen ｅr －Ｈolf)方程。
解此方程可得到最优权系数012,,,h h h ，此式是Wiener 滤波器的一般方程。
定义121M M h h h h ⨯⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦()()()1011xd xd xd xd M r r R r M ⨯⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥-⎢⎥⎣⎦()()()()()()()()()011102120xx xx xx xxxx xx xx xx xx xx M Mr r r M r r r M R r M r M r ***⨯⎡⎤-⎢⎥-⎢⎥=⎢⎥⎢⎥--⎢⎥⎣⎦则维纳-霍夫方程可写成矩阵形式xd xx R R h =求逆,得1xx xd h R R -=此式表明，已知期望信号与观测数据的互相关函数及观测数据的自相关函数时,可以通过矩阵求逆运算，得到维纳滤波器的最佳解。
3. 实验结果及分析（1)当L=2００,N=6信噪比为20ｄB 的滤波效果信噪比为１０dＢ的滤波效果信噪比为6dB的滤波效果（2)当L＝2０0，N=60信噪比为20ｄB的滤波效果信噪比为1０dＢ的滤波效果信噪比为6dB的滤波效果（３)当L=６00，Ｎ＝6信噪比为20dB的滤波效果信噪比为10dＢ的滤波效果信噪比为６dB的滤波效果实验分析:别取信号长度为200、600，滤波器长度为6、６0,加噪信噪比为20dB、10dB、6dB，组合进行实验。
每组实验得到的最小均方误差统计如下表。
由此表可以看出，信号长度越长，最小均方误差（绝对值）越大，精度越差；在信噪比较大（误差影响较小）的滤波过程中,滤波器长度约长,最小均方误差（绝对值)越小,精度越好。
表1 最小均方误差统计表加噪信号信噪比L=2０0N＝6L=200Ｎ＝6０L=600N＝61()x n20dB －178５．３971-１574.１1０7-6８8８.41032()x n10dB -845.５714 -１727.5260-1021９．33973()x n６ｄB -1213.636７-2１3５.459６－77４3．9３5８对于相同信号和滤波器(这里取L=２00,N=6),信噪比越大，最小距离误差约小;而当信噪比较小时,信号与噪声值接近,导致滤波效果受到影响,最小距离误差变大。
4. 源代码clｅaｒ；clc；％初始化变量L＝200；%信号长度N=6;%滤波器的阶次a=０.96;ｗn=randｎ（L,1）;%wn为用于生成信号的噪声信号,随机生成一个Ｌ*1矩阵,生成矩阵的元素值在％区间（０.０，１．0）之间sn=ｚｅroｓ(L,１）;%sｎ为信号，生成一个Ｌ＊1的零矩阵ｈn=zeros（N,1)；%hｎ为系统单位脉冲响应生成一个N*１的零矩阵rxx=zｅros(N，１）;%rxx为自相关函数,生成一个N＊1的零矩阵rｘd=zeroｓ(Ｎ，1);％ｒｘd为互相关函数,生成一个N＊1的零矩阵ｙn＝zeros(L,1）;%ｙn为输出信号,生成一个L*１的零矩阵xt=ｚerｏs（L＋N,1）；%生成一个（L+Ｎ)＊１的零矩阵gn=zｅｒos（L,1);%gn为yn与sn最小距离误差信号，生成一个L*1的零矩阵％根据给定公式ｓ（ｎ）=aｓ(n-１）+w(n）,生成理想信号for i=2:Ｌsn(ｉ,１)=a＊ｓn(i-1,1)+ｗn（i,1);ｅｎdsn(1，1)=ｗn(1,１);suｂpｌot(2,2，1)；pｌot（ｓn,'r＇）,ａxｉs([0 2０0 -10 10])，ｘlabeｌ('时间'),ｙlaｂe ｌ('幅度'),ｔｉｔle（＇sn')%生成期望信号方差ｃdcd=(ｖaｒ（wn))/(１－a＾2）；%对信号加噪ｘ1=aｗgn(sｎ,20);x2=ａｗgｎ(sn，１0)；x３=awgn(sn,6）；subploｔ(２,2，2)plot（ｘ3,＇g')，axiｓ(［0 ２00 -1０ 10]）,xｌａbeｌ('时间＇）,ｙｌab ｅl(＇幅度')，titｌe('x3')；%生成输入信号与理想信号的互相关函数,此处x1为输入信号，ｓn为期望信号for i=1:N，for ｍ=i:1:L,rxd(i，1)＝rxd(i,1)+x３（m,１）＊ｓn(m-ｉ＋1，1);enｄｅnd％生成输入信号的自相关函数foｒ i=1：N,for ｍ=i：１:L,rxx(i,１)=rxx(i，1)+x３（m，1)*x3(m-ｉ+1，1);ｅｎdend%将自相关函数生成托普勒斯矩阵rxx1=ｔoepｌｉｔｚ(rxx）;%生成逆矩阵iｒｘx=inｖ(rｘx1）;%生成滤波器系数ｈ(n)hn=irxx*rxd;ｆｏr i=1:Lxt(ｉ＋N,１)=ｘ3（i,１);eｎｄ%实现滤波for i＝１：L,for m=１:Ｎ，ｙn(ｉ,1）=ｙn(ｉ,1)＋ｘt(i+N+１-m,1)*ｈn(m,1);eｎｄend％计算最小均方误差信号eｎeｎ=０;eｎ=cｄ-(rxｄ＇)*hn;%生成最小距离误差信号gnｇn＝yｎ-sｎ；%画出滤波后的信号时域图sｕｂpｌｏt(2,2,3);plot（yn),aｘiｓ（[0 2００-1０ 1０]),xｌabｅｌ('时间＇),ｙｌabｅl（'幅度＇),titｌe('yn＇);%画出理想信号与输出信号对比图ｓｕbplot(2,2，4）;ｐlot(sｎ,'r')，aｘis([0 200 -１０ 10]),ｘlaｂｅl('时间'）,ｙｌabel （'幅度＇),title('ｓn与ｙn对比')；hold ｏn;pｌot(yｎ，＇b')，axiｓ(［0 2０0 -１0 10］)；hｏld oｆｆ;%画出最小距离误差信号图ｆigure；plｏt（gｎ）,ａxiｓ([０２0０－2 2]),ｘlaｂｅl('时间'),yｌabel('幅度'),titlｅ(＇gn'）;。