现代信号处理研究生课程报告

1、华南师范大学现代信号处理课程设计课程名称：现代信号处理课程题目：wiener滤波器和kalman滤波器的原理分析及其matlab实现指导老师：李xx专业班级：2015级电路与系统姓名：xxxx学号：xxxxwiener滤波器和kalman滤波器的原理分析及matlab实现摘要：信号处理的实际问题，常常是要解决在噪声中提取信号的问题，因此，我们需要寻找一种所谓有最佳线性过滤特性的滤波器。这种滤波器当信号与噪声同时输入时，在输出端能将信号尽可能精确地重现出来，而噪声却受到最大抑制。Wiener滤波Kalman滤波就是用来解决这样一类从噪声中提取信号问题的一种过滤（或滤波）方法1。Wiener滤波与

2、Kalman滤波都是解决最佳线性过滤和预测问题，并且都是以均方误差最小为准则的。但与Wiener滤波器不同的是，Kalman滤波器是一种自适应滤波器，Kalman滤波器提供了推导称作递推最小二乘滤波器的一大类自适应滤波器的统一框架。关键词：Wiener滤波Kalman滤波均方误差最小自适应滤波器第一章绪论41.1 滤波器的发展历程41.2 现代信号处理的滤波器分类51.3 wiener和kalman滤波各自的运用领域61.3.1 wiener滤波的运用范围61.3.2 kalman滤波的运用范围6第二章wiener和kalman的各自的滤波原理72.1 wiener滤波器的原理分析72.2 2

3、维纳-霍夫方程92.3 kalman滤波的自适应原理分析112.4 wiener滤波和kalman滤波的区别与联系13第三章wiener和kalman滤波的matlab仿真实现143.1 FIR维纳滤波器的matlab实现143.2 kalman滤波器的matlab实现19第四章总结与展望23参考文献25第一章绪论1.1 滤波器的发展历程从滤波器的发展现状来看，滤波器从处理信号的类型可以分为模拟滤波器和数字滤波器，模拟滤波器可分为无源滤波器(Passivefilter)和有源滤波器(Activefilter),而数字滤波器已可用计算机软件实现，也可用大规模集成数字硬件实时实现。本文主要针对数字

4、滤波器。从形式上看，数字滤波有线性滤波和非线性滤波。线性滤波是指卷积滤波，又分为频域滤波和时域滤波，在实域中根据滤波方式又分为递归滤波和递归滤波。非线性滤波主要是指同态滤波，它是用取对数的方法将非线性问题线性化。近些年，线性滤波方法，如Wiener滤波、Kalman滤波和自适应滤波得到、2了广泛的研究和应用。同时一些非线性滤波方法，如小波滤波、同态滤波、中值滤波和形态滤波等都是现代信号处理的前沿课题，不但有重要的理论意义，而且有广阔的应用前景。Wiener滤波是最早提出的一种滤波方法，当信号混有白噪声时，可以在最小均方误差条件下得到信号的最佳估计。但是，由于求解Wiener-Hoff方程的复杂

5、性，使得Wiener滤波实际应用起来很困难，不过Wiener滤波在理论上的意义是非常重要的，利用Wiener滤波的纯一步预测，可以求解信号的模型参数，进而获得著名的Levinson算法。Kalman滤波是20世纪60年代初提出的一种滤波方法。与Wiener滤波相似，它同样可以在最小均方误差条件下给出信号的最佳估计。所不同的是，这种滤波技术在时域中采用递推方式进行，因此速度快，便于实时处理，从而得到了广泛的应用。Kalman滤波推广到二维，可以用于图象的去噪。当假设Wiener滤波器的单位脉冲响应为有限长时，可以采用自适应滤波的方法得到滤波器的最佳响应。由于它避开了求解Wiener-Hoff方程

6、，为某些问题的解决带来了极大的方便阔。小波滤波就是利用信号和噪声的目的。同态滤波主要用于解决信号和噪声之间不是相加而是相乘关系时滤波问题。另外，当信号和噪声之间为卷积关系的时候，在一定条件下可以利用同态滤波把信号有效地分离开来，由同态滤波理论引申出的复时谱也成为现代信号处理中极为重要的概念.Wiener滤波、Kalman滤波等自适应滤波都是线性滤波，线性滤波的最大缺点就是在消除噪声的同时，会造成信号边缘的模糊。中值滤波是20世纪70年代提出的一种非线性滤波方法，它可以在最小绝对误差条件下，给出信号的最佳估计。这种滤波方法的优点，就是能够保持信号的边缘不模糊。另外它对脉冲噪声也有良好的消除作用。

7、形态滤波是建立在集合运算上的一种非线性滤波方法，它除了用于滤除信号中的噪声外，还在图象分析中发挥了重要的作用。1.2 现代信号处理的滤波器分类数字滤波分空间域和频率域的方法。空间域的滤波处理，是根据平滑窗口内的统计值或自适应参数进行处理，很难达到在消除相干斑噪声的同时又能很好地保留边缘和纹理细节的理想状态。一般只能在相干斑噪声消除和细节信息保留两个方面进行折衷，综合这两个方面的较好效果。频率域的傅立叶变换能够进行高频或低频的带通滤波，但不能区分噪声和信息相近的频率。基于小波分析的方法由于具有多分辨率和时频联合分析的特征，使得频率域的去噪有了更好的途径。空间域的几种著名滤波器可分为以下两类：传统

8、方法、局域统计自适应滤波方法。均值滤波器和中值滤波器属于经典传统滤波器范畴。传统方法在对SA就像进行滤波时，对噪声和边缘信息是不加区分的。为了解决传统方法存在的问题，人们提出了各种形式的自适应滤波器，自适应滤波器一般通过局域统计参数的调节，对噪声进行较强的平滑，而对边缘则尽量予以保留3o基于频率域的滤波方法有Fourier变换滤波方法以及近年兴起的小波变换滤波的方法，传统的建立在傅里叶变换基础上的频率域滤波方法在提高信噪比和提高空间分辨率两项指标上存在矛盾。低通滤波能较好地平滑抑制噪声，但同时也模糊了图像的边缘。高通滤波可以使边缘更加陡峭，但背景噪声同时也被加强。此外相干平均也是滤除噪声常用的

9、手段，但需时间较长，不能作动态提取，而且当各次纪录中的信号没有对齐时处理结果也会产生低通模糊。与之相比，基于小波变换的多分辨率滤波技术有明显优点。小波分析最大的特点在于具有极敏感的变焦特征，在不同的分辨率下，反映出不同的图像结构特征，使其在处理突变信息方面具有特殊的能力，利于噪声的滤除和边缘的保留。1.3 wiener和kalman滤波各自的运用领域Weiner滤波是一种线性滤波的方法，而kalman是对Weiner滤波理论的发展，即与weiner滤波器不同的是，kalman滤波器是一种自适应滤波器，kalman滤波器提供了推到称作递推最小二乘滤波器的一大类自适应滤波器的统一框架。1.3.1

10、wiener滤波的运用范围维纳滤波通常用于深层地震勘探数据处理，特别在于深层石油天然气勘探开发中有诸多优点，然而当我们将维纳滤波的思路与原理应用于浅层地震勘探的数据处理中时，发现维纳滤波仍然能够较好地滤除噪音信号，提高信号的信噪比。微地震勘探技术中，信噪比是衡量地震资料好坏的一个重要指标，信噪比越高，则地震资料质量越好，处理结果就越可信。所以，信噪比的估值无论是对处理资料还是对地质解释都有一定的参考价值。微地震勘探中，噪声是不可避免的，提高信噪比是地震资料数据处理中一项最基本的任务。目前数字滤波技术的应用是提高信噪比最常用的方法。维纳滤波是数字信号处理中滤波技术研究的一个主要内容。作为最佳滤波

11、器，其最优化的准则是使均方误差最小。1.3.2 kalman滤波的运用范围简单来说，卡尔曼滤波器是一个最优化自回归数据处理。对于解决很大部分的问题，他是最优，效率最高甚至是最有用的。他的广泛应用已经超过30年,包括机器人导航，控制，传感器数据融合甚至在军事方面的雷达系统以及导弹追踪等等。近来更被应用于计算机图像处理，例如头脸识别，图像分割，图像边缘检测等等。卡尔曼滤波是根据前一个估计值和最近一个观察数据来估计当前值，用状态方程和递推方法估计，解是以估计值形式给出。卡尔曼滤波不要求信号和噪声都是平稳过程的假设条件。对于每个时刻的系统扰动和观测误差(即噪声)，只要对它们的统计性质作某些适当的假定，

12、通过对含有噪声的观测信号进行处理，就能在平均的意义上，求得误差为最小的真实信号的估计值。因此，自从卡尔曼滤波理论问世以来，在通信系统、电力系统、航空航天、环境污染控制、工业控制、雷达信号处理等许多部门都得到了应用，取得了许多成功应用的成果。例如在图像处理方面，应用卡尔曼滤波对由于某些噪声影响而造成模糊的图像进行复原。在对噪声作了某些统计性质的假定后，就可以用卡尔曼的算法以递推的方式从模糊图像中得到均方差最小的真实图像，使模糊的图像得到复原。第二章wiener和kalman的各自的滤波原理2.1wiener滤波器的原理分析维纳(Wiener)是用来解决从噪声中提取信号的一种过滤(或滤波)方法。这

13、种线性滤波问题，可以看做是一种估计问题或一种线性估计问题。一个线性系统，如果它的单位样本响应为h(n),当输入一个随机信号x(n),且x(n)=s(n)+v(n)(1)其中x(n)表示信号，v(n)表示噪声，则输出y(n)为(2)我们希望x(n)通过线Tt系统h(n)后得到的y(n)尽量接近于s(n),因此称y(n)为(3)s(n)的估计值，用s(n)表示，即y(n)=s(n)则维纳滤波器的输入一输出关系可用下面图1表示。=昌(冗)+乂月)*细)图1实际上，式(2)所示的卷积形式可以理解为从当前和过去的观察值x(n),x(n-1),x(n-2),x(n-m),来估计信号的当前值s(n)。因此，

14、用h(n)进行过滤问题实际上是一种统计估计问题。一般地，从当前的和过去的观察值x(n),x(n-1),x(n-2),估计当前的信号值y(n)=s(n)成为过滤或滤波；从过去的观察值，估计当前的或者将来的信号值y(n)=s(n+N)(N20)称为外推或预测；从过去的观察值，估计过去的信号值y(n)=s(n-N)(N1)称为平滑或内插。因此维纳滤波器又常常被称为最佳线性过滤与预测或线性最优估计。这里所谓的最佳与最优是以最小均方误差为准则的。e(n)=s(n)-s(n)(4)显然e(n)可能是正值，也可能是负值，并且它是一个随机变量。因此，用它的均方误差来表达误差是合理的，所谓均方误差最小即它的平方

15、的统计期望最小：(5)(n)=Ee2(n)=min采用最小均方误差准则作为最佳过滤准则的原因还在于它的理论分析比较简单,不要求对概率的描述，维纳滤波基本原理框图如下所示：轼入2.2维纳-霍夫方程设计维纳滤波器的过程就是寻求在最小均方误差下滤波器的单位脉冲响应h(n)或传递函数H(z)的表达式，其实质就是解维纳-霍夫(Wiener-Hopf)方程。我们从时域入手求最小均方误差下的h(n),用hopt(n)表示最佳线性滤波器。这里只讨论因果可实现滤波器的设计。因果的维纳滤波器，设h(n)是物理可实现的，也即是因果序列：h(n)=0,当n0因此，从式上式中可推导：-bey(n)=?(n)=、h(m)

16、x(n-m)m=0Ee2(n)LE|rs(n)-fh(m)x(n-m)i|IVm=0/要使得均方误差最小，则将上式对各h(m),m=0,1,，求偏导，并且等于零，得0=0j=0,1，即*、n)-hopt(m)x(n-m)x(n-j)m卫ES(n)x(n一j)=h0Pt(m)Ek(nm)x(nj)j0。m卫用相关函数R来表达上式，则得到维纳-霍夫方程的离散形式:-boRx八hopt(m)Rxx(j-m)j_0mz0由上式进一步化简得：Ee2(n)Lin=Rxx(0)-xhopt(m)Rxx(m)m-0为了按(5)式所示的最小均方误差准则来确定维纳滤波器的冲激响应h(n),令n)对h(j)的导数等

17、于零，即可得Rxs(m)=Lh(i)Rxx(m-i),vm(6)式中，Rxs(m)是s(n)与x(n)的互相关函数，Rxx(m)是x(n)的自相关函数，分别定义为Rxs=Ex(n)s(nm)Rxx=Ex(n)x(nm)式(6)称为维纳滤波器的标准方程或维纳-霍夫(Wiener-Hopf)方程。如果已知Rxs(m)和Rxx(m),那么解此方程即可求的维纳滤波器的冲激响应。式(6)所示标准方程右端的求和范围即i的取值范围没有具体标明，实际上有三种情况：(1)有限冲激响应(FIR)维纳滤波器，i从0到N-1取得有限个整数值；(2)非因果无限冲激响应(非因果IIR)维纳滤波器，i从-8到+如取所有整数

18、值；(3)因果无限冲激响应(因果IIR)维纳滤波器，i从0到+8取正整数值。上述三种情况下标准方程的解法不同，本文只描述FIR维纳滤波器的求解。设滤波器冲激响应序列的长度为N,冲激响应矢量为h=h(0)h.h(N-1)T(7)滤波器输入数据矢量为：x(n)=x(n)x(n-1).x(n一N+1)T(8)则滤波器的输出为：y(n)=s(n)=xT(n)h=hTx(n)(9)这样，式(6)所示的维纳-霍夫方程可写成PT=卜*或P=Rh(10)其中P=Ex(n)s(n)(11)是s(n)与x(n)的互相关函数，它是一个N维列矢量；R是x(n)的自相关函数,是N阶方阵R=Ex(n)xT(n)(12)利

19、用求逆矩阵的方法直接求解式(10),得1hpt=RP(13)这里opt表示“最佳”，这就是FIR维纳滤波器的冲激响应。2.2 kalman滤波的自适应原理分析卡尔曼滤波的含义是现时刻的最佳估计为在前一时刻的最佳估计的基础上根据现时刻的观测值作线性修正。卡尔曼滤波在数学上是一种线性最小方差统计估算方法，它是通过处理一系列带有误差的实际测量数据而得到物理参数的最佳估算。其实质要解决的问题是要寻找在最小均方误差下Xk的估计值Xknn15:_二!|三二匚-二103：040C600600EJmar：Ffer-2imulztisnfileEd*VinffirtwrtTgl，DrH叩VWndowHelp士i

20、d-,ljrn-第四章总结与展望信号处理的实际问题，常常是要解决在噪声中提取信号的问题，因此，我们需要寻找一种所谓有最佳线性过滤特性的滤波器。这种滤波器当信号与噪声同时输入时，在输出端能将信号尽可能精确地重现出来，而噪声却受到最大抑制。维纳(Wiener)滤波与卡尔曼(Kalman)滤波就是用来解决这样一类从噪声中提取信号问题的一种过滤(或滤波)方法。实际上这种线性滤波问题，可以看成是一种估计问题或一种线性估计问题。维纳过滤与卡尔曼过滤都是解决最佳线性过滤和预测问题，并且都是以均方误差最小为准则的。因此在平稳条件下，它们所得到的稳态结果是一致的。然而，它们解决的方法有很大区别。维纳过滤是根据全部过去的和当前的观察数据式门一1)：犬5-2).来估计信号的当前值，它的解是以均方误差最小条件下所得到的系统的传递函数H(z)或单位样本响应h(n)的形式给出的,因此更常称这种系统为最佳线性过滤器或滤波器。而卡尔曼过滤是用前一个估计值和最近一个观察数据来估计