《第9章维纳滤波ppt课件.ppt》由会员分享，可在线阅读，更多相关《第9章维纳滤波ppt课件.ppt（96页珍藏版）》请在三一办公上搜索。
1、2022/11/13,1,第九章维纳滤波(Wiener Filtering),2022/11/13,2,主要内容,9.1 概述9.2 波形线性均方估计的正交原理9.3 维纳-霍夫（Wiener-Horf)积分方程9.4 非因果的维纳滤波问题9.5 因果的维纳滤波器9.6 预测问题9.7 后验维纳滤波与互补维纳滤波9.8 维纳滤波器的应用,2022/11/13,3,在实际应用中，有用信号往往会受到一些外界干扰，我们实际观察到的是受到噪声干扰了的信号。如何最大限度地抑制噪声，并将有用信号分离出来，是信号处理中经常遇到的问题。,9.1 概述,2022/11/13,4,在传输或测量信号s(n)时，由于
2、信道噪声或者测量噪声w(n)，接收或测量到的数据x(n)将与s(n)不同。设噪声是加性的: 即:x(n)=s(n)+w(n),2022/11/13,5,如果s(n)和w(n)的频谱是分离的，那么设计一个具有恰当频率特性的线性滤波器即能有效地抑制噪声并提取信号，这就是前面经典数字信号处理理论中详细讨论过的数字滤波器的设计问题.但是如果s(n)和w(n)的频谱互相重叠,或者s(n)和w(n)是随机信号，它们的频谱根本就不存在，问题就要复杂得多，这就是本章要讨论的内容。,2022/11/13,6,随机性是生物医学信号的特点之一，在本章中主要讨论噪声中随机信号的线性估计问题。维纳滤波适用于平稳随机过程
3、。观察x(t)中既含有随机信号s(t)又含有噪声n(t)。,2022/11/13,7,经处理器处理后得一估计值 作为对所希望取得的信息d(t) 的估计值， d(t)可能是s(t)，也可能是预测值s(t+a)，导数ds(t)/dt等。估计的任务就是要求 与d(t) 的差值在一定判据意义下取得极小值。,2022/11/13,8,处理器,判据,极小,9-1估计原理方框图,2022/11/13,9,根据待估计量d(t)的形势，波形估计问题可分为三类：,滤波问题：由t0tf一段时期内的观察x(t), t0ttf，估计t=tf瞬间信号s(t)的值s(t)。即：d(t)=s(t)。预测问题：由t0tf一段时
4、期内的观察x(t), 估计ttf的某一时刻待估计信号的可能值。即：d(t)=s(t+a)，a0。平滑问题：由t0tf一段时期内的观察x(t), 估计t0ttf期内待估计信号的值.也叫内插问题。,2022/11/13,10,在本章应用中作以下假设：,x(t)中信号s(t)和噪声信号n(t)是加法地结合的，即：x(t)=s(t)+n(t)。处理器所采用的算法限于线性运算。即：如该处理器的冲击响应是h(t),则：所采用的判据是均方误差最小，即：,2022/11/13,11,9.2 波形线性均方估计的正交原理,设信号是随机过程s(t)，观察是t0tf期内测得的随机过程x(t), tt0，tf,为了简化
5、讨论，设s(t)和x(t)都是零均值的。采用最小均方误差作为估计判据。即：又限定估计 是由观察x(t)经线性滤波器h(t)得出的：,2022/11/13,12,最优线性均方估计的选取原则是使估计误差 与所有的观察值x(), t0，tf正交，也就是说，如果对每一个 t0，tf都有：,则均方误差最小，它等于：,2022/11/13,13,证明：,2022/11/13,14,例9.1 简单预测问题,设观察中没有噪声，即x(t)=s(t)，又待估计量是信号的预测值d(t)=s(t+a)，a0，设只用t时刻的观察值x(t)对d(t)作线性估计：按最小均方误差判据做估计，即求。,2022/11/13,15
6、,解：根据正交原理，估计误差：,应和观察值x(t)=s(t)正交，即：,2022/11/13,16,例 9.2,设观察中没有噪声，即x(t)=s(t)，又待估计量是信号的预测值d(t)=s(t+a)，a0，设估计算子采用：按最小均方误差判据做估计，即求估计系数a和b。,2022/11/13,17,解：此时正交原理表现为：,2022/11/13,18,由于Rss(t)是奇函数，所以Rss(0)=0把上式化简得到：,2022/11/13,19,代入均方误差式中，得到：,因为Rss(0)0,此式最后一项大于零，所以，它要比例9.1的LMS要小，主要是他利用了更多的测量信息，估计效果更好些。,2022
7、/11/13,20,9.3维纳-霍夫（Wiener-Horf)积分方程,观察x(t)由信号s(t)和噪声n(t)相加组成，观察时间t0tf，则： x(t)=s(t)+n(t), t0，tf 待估计过程是d(t)，x(t)经线性处理得到估计为：要求估计均方误差最小，试求h(t),2022/11/13,21,根据正交原理可知：,即t时刻的估计误差要和t0tf区间所有时刻的观察值x()正交，推得：,2022/11/13,22,这就是h(t)应满足的条件,称为维纳-霍夫积分方程，只要相关函数Rxd和Rxx已知，就可以由此解出h(t).而h(t)一经解出，就有：,问题是维纳-霍夫方程是一个积分方程，未必
8、能求出解析解答。,2022/11/13,23,9.4 非因果的维纳滤波问题,对于滤波问题，利用从t0=-一直到tf=t时刻为止的全部观察来估计t时刻的信号s(t)。此时有： t0=-,tf=t d(t)=s(t) 于是维纳-霍夫方程变为：,2022/11/13,24,做变量替换，t-=,t-=,得到：,或：此时：,2022/11/13,25,9.4.1 连续时间形式的解,如果不要求滤波器是因果的，可以把观察时间的上限tf扩展到-,也就是利用x(t)在全时间轴上的值来进行估计。此时维纳-霍夫方程为：,2022/11/13,26,把它做付氏变换有：,如果n(t)和s(t)统计独立，则有：,做反傅里
9、叶变换得到滤波器的冲击响应h(t),2022/11/13,27,例9.3 设信号的功率谱是：,噪声是白色的，其功率谱是常数而且噪声与信号统计独立，求维纳滤波器的频率特性和冲击响应。,2022/11/13,28,解：此时有,求付里叶反变换得到：它显然是非因果的，物理不可实现的,2022/11/13,29,在离散情况下，在不要求物理可实现的条件下。可以类似推出以下结论：维纳-霍夫原方程为：现在放宽为：,9.4.2 离散时间形式的解,2022/11/13,30,因此滤波器的频率特性是：,实际中一般采用Z变换的传递函数,将H(z)做反演Z变换得到冲击响应h(n),2022/11/13,31,可见H(e
10、j) 决定于信号与噪声的功率谱密度;当噪声为零时,即Snn(ej) =0; H(ej) =1，信号全部通过；当信号为零时,即Sss(ej) =0; H(ej) =0 ，噪声被全部抑制掉;因此维纳滤波器确有滤除噪声的能力。非因果维纳滤波器的幅频特性如下图所示。,2022/11/13,32,2022/11/13,33,0,非因果维纳滤波器的幅频特性,w,Sss(ej),Snn(ej),H(ej),1,2022/11/13,34,例9.4 设信号的自相关函数是：,噪声是白色的设计非因果的维纳滤波器,2022/11/13,35,传递函数：,解：此时有,2022/11/13,36,它是非因果的，而且是无
11、限长的，可以取短近似，如只取4项为：,2022/11/13,37,9.5因果的维纳滤波器,非因果维纳滤波器需要用全时间上的观察值来估计s(n),所以不能实时实现，即使采用把h(n)截短的近似估计，也必须延迟若干拍，待xn+k输入后（k是截短范围)才能做出本次估计。,维纳滤波器的时域解（Time domain solution of the Wiener filter）,2022/11/13,38,设计维纳滤波器的过程就是寻求在最小均方误差下滤波器的单位脉冲响应h(n)或传递函数H(Z)的表达式，其实质就是解维纳霍夫（WienerHorf）方程。我们从时域入手求最小均方误差下的h(n)。这里只讨
12、论因果可实现滤波器的设计，既：,2022/11/13,39,在此主要介绍两种方法：,FIR型：限制处理器的形式，只用最近的p+1个观察值x(n),x(n-1),x(n-p)来估计s(n),即：预白化处理，把观察序列值x(n)白化。变成白噪声w(n)，再对w(n)做可实现的最优滤波，如图：,X(n),d(n),预白化,H2,w(n)白色,2022/11/13,40,9.5.1 FIR型处理,由正交原理得：,2022/11/13,41,令:m=n-m,上式可改写为：,用矩阵表示： RxxH=G,2022/11/13,42,只要Rxx是非奇异的，就可以求得H:,H=Rxx-1GRxx是对称且Topl
13、itz型的这时的最小均方误差为：,2022/11/13,43,随机信号都可以看成是由一白色噪声W(n)激励一个物理可实现的系统或模型的响应，如图所示：,信号模型图,9.5.2 预白化处理,2022/11/13,44,预白化方法是基于如下事实：当x(n)是方差x2=1的白噪声时，有： Rxx(m-n)=1 当m=n =0 其它 所以上式直接就可以得出： h(m)=Rxs(m) m=0,离散维纳-霍夫方程为:,2022/11/13,45,对白色的x(n)而言，它的维纳解是Rxs(m)的 m0的部分也可以写成:h(m)=Rxs(m) u(m),它的付氏变换写作： Sxs(ej)+或用Z变换表示： H
14、(z)= Sxs(z)+符号 . +表示原函数m0部分对应的付氏变换和Z变换。,2022/11/13,46,X(n)S(n)+n(n),d(n)=s(n),预白化,H2(z),w(n)白色,H1(z),E.2,极小,优化,-,+,预白化算法处理框图,2022/11/13,47,复习一个Z变换的性质,如果h(n)=h(-n),则有:H(z)=H(z-1);那么，如果z1是()的极点，1/Z1一定是H(z)的极点;同样，零点也有这样的性质；还有如果h(n)是实函数，则(z)极点一定是共轭对称的,2022/11/13,48,1.预白化滤波器H1(z)的设计 对x(n)的可实现白化滤波器H1(z)可如
15、下求得:,2022/11/13,49,2.最优滤波器H2(z)的设计 因为W(n)是白色的，所以：,2022/11/13,50,最后得总滤波器:,2022/11/13,51,证明：,2022/11/13,52,例9.5 设x(n)=s(n)+n(n),s、n统计独立,且： 设计可实时实现的维纳滤波器。 解：,2022/11/13,53,因此得：,可见：,2022/11/13,54,括弧中的因子可按下式做部分分式分解：前一项对应与n0部分，后一项对应与n0 部分。分解得A=B=0.6。,2022/11/13,55,因此得:,2022/11/13,56,这种处理器有两种具体实现方法：,()把h(n
16、)截短成有限长，再与x(n) 做卷积：()按IIR方式递归实现：,2022/11/13,57,9.6 预测问题,预测的任务是根据观察x(n)对d(n)的未来值做预报。提法全然不变，只是待估过程d(n)的具体含义变了，因此估计结果也改变。在最常见的一步预测情况下： 因此：,2022/11/13,58,9.6.1一般解答,因为此时维纳-霍夫方程中的唯一变化是Rxd(m) 成为Rxs(m+1) ,所以一步预测时离散时间的维纳-霍夫方程是:,2022/11/13,59,(i)在非因果情况下，把上式放宽成： 从而得：,当s(n)和n(n)互相独立时，上式变成：,2022/11/13,60,（ii）在因果
17、情况下，解答是：例9.6 试为例9.5设计可实现的一步预测器。 因为：,2022/11/13,61,本例中：做部分分式分解，上式可分成:前一项对应于n0部分，故:,2022/11/13,62,因此:得:或用IIR方式实现：,2022/11/13,63,9.6.2 用有限项FIR滤波器实现,此时:由正交原理：得维纳-霍夫方程组:,2022/11/13,64,用矩阵形式表示:这里的自相关阵仍是对称且Toeplitz型的。,2022/11/13,65,例9.7 用p=2的FIR结构给例9.5设计维纳一 步预测器。,2022/11/13,66,故得:解得：故,2022/11/13,67,9.7 后验维
18、纳滤波与互补维纳滤波,9.7.1 后验维纳滤波维纳滤波是以信号和噪声的相关函数或功率谱已有先验知识为前提的。如果这些统计特性未知，就需先作出它们的估计，然后再据以设计维纳滤波器。但是从严格意义上说，这时所得结果并不是真正的维纳滤波器，所以称之为后验“维纳”滤波。,2022/11/13,68,从频域上应用后验维纳滤波的核心问题是由各次观察xi(n)中分解出信号的谱估计 和噪声的谱估计 。通常可采用下述步骤：1.先对各次观察求均值，设做N次观察： 式中，s(n)是确定性的诱发响应，ni(n)是第i次刺激后记录中的噪声。,2022/11/13,69,则平均诱发响应：然后求 的功率谱。如果s(n)和n
19、i(n) 统计独立，各次噪声也互相独立，则:,2022/11/13,70,2.再分别对每一次观察xi(n) 求功率谱:并求这些功率谱的平均值:,2022/11/13,71,3.联立解 1 ,2中最后两式，便可求得Sss(ej) 和Snn(ej) 的估计：,2022/11/13,72,据此，得后验维纳滤波器如下：（i）用于对单次观察进行滤波：（ii）用于对平均诱发响应x(n)进行滤波：,2022/11/13,73,许多研究者用这种滤波方法对各种平均诱发响应进行了滤波，但效果报道不一。有的效果较好，有的却不甚见效。其原因除了谱估计不是真实值，因此所得得H()只能是近似的估计外，还由于“维纳”滤波的
20、其它假设也未必能满足。其中： 过程不是平稳的;“信号和噪声是相加的”这一假设是一个有用模型，但刺激愈接近阈值不正确；信号与噪声未必统计独立。实际上刺激对作为噪声的自发活动往往也有一些作用。,2022/11/13,74,为了改进后验“维纳”滤波的效果。又做出了许多改进方案，介绍如下：1.交替集均法 此法除按前式求外,又按下式计算另一种平均值： 即：每当序号i 为偶数时，就将观察值取负号，通过这样的“相加”，S(n)将被平均掉，因此 的功率谱将只反映噪声：,2022/11/13,75,便可得维纳滤波器。 这个方法的优点式计算量大为下降：只需要求两次功率谱一次对 ，一次对 ，而采用前法却需要求N+1
21、次功率谱对每个xi(n) 求 ，还要对 求 。但理论分析可以证明，所得谱估计方差较大是这种方法的缺点。,2022/11/13,76,2.谱平滑：把前面得到的功率谱 和 加以平滑，然后再代维纳滤波公式，可以改善滤波效果。设施加在 和 上的谱窗口分别是1(ej)和2(ej)，则平滑后的谱分别是:,2022/11/13,77,然后把它们代入维纳滤波公式，得：窗口长度要比较数据长度短，其具体值要在方差和偏差之间取折中。时窗宽则谱窗窄，因此平滑作用小，偏差小，方差大。反之，时窗窄则谱窗宽，平滑作用显著，因此偏差大，方差小。,2022/11/13,78,9.7.2 互补维纳滤波,维纳滤波器得基本假设是信号
22、为随机的，但是实际工作中信号常有些确定性结构，并非纯粹随机，因此应用效果未必好。因此对它简单应用维纳滤波效果未必好，因为这样处理所得的 充其量也只是真实 S(t) 在最小均方误差意义下的逼近，不是真实S(t) 。对这类S(t) 是确定性信号的情况，采用互补维纳滤波可能是更合理的方案。,2022/11/13,79,以做两次观察为例。如果设计滤波器时，H1(z)和H2(z)是分别独立设计的,如图(a)然后再把处理结果相加，效果就未必好.,2022/11/13,80,如果设计时多引入一个限制条件，如： 效果就会好些，如图(b)。,1-H1(z),H1(z),s(n)+n2(n),s(n)+n1(n)
23、,+,+,x(n),x2(n),x1(n),2022/11/13,81,因为此时:,H1(z),s(n)+n2(n),s(n)+n1(n),+,+,x(n),+,n1(n)-n2(n),2022/11/13,82,可见处理结果中信号s(n)成分不变。H1(z)的任务是把n1-n2变成对n2的最优抵消。由于n2和n1-n2都是随机信号,所以应用维纳滤波的效果就比较好。,2022/11/13,83,9.8维纳滤波器的应用（Application of Wiener filter）要设计维纳滤波器必须知道观测信号和估计信号之间的相关函数，即先验知识。如果我们不知道它们之间的相关函数，就必须先对它们的
24、统计特性做估计，然后才能设计出维纳滤波器，这样设计出的滤波器被称为“后验维纳滤波器”。,2022/11/13,84,9.8.1.在生物医学信号处理中比较典型的应用就是关于诱发脑电信号的提取。大脑诱发电位（Evoked Potential，EP）指在外界刺激下，从头皮上记录到的特异电位，它反映了外周感觉神经、感觉通路及中枢神经系统中相关结构在特定刺激情况下的状态反应。在神经学研究以及临床诊断、手术监护中有重要意义。,2022/11/13,85,EP信号十分微弱，一般都淹没在自发脑电（EEG）之中，从EEG背景中提取诱发电位一直是个难题：EP的幅度比自发脑电低一个数量级，无法从一次观察中直接得到；
25、EP的频谱与自发脑电频谱完全重迭，使得频率滤波失效；在统计上EP是非平稳的、时变的脑诱发电位。通过多次刺激得到的脑电信号进行叠加来提取EP，这是现今最为广泛使用的EP提取方法。,2022/11/13,86,为了解决诱发电位提取问题，研究者利用维纳滤波来提高信噪比，先后有Walter、Doyle、Weerd等对维纳滤波方法进行了改进。在频域应用后验维纳滤波的核心就是由各次观察信号中分解出信号的谱估计和噪声的谱估计，通过设计出的滤波器来提高信噪比。,2022/11/13,87,9.8.2. 介绍时频平面的维纳滤波（timefrequency plane wiener filtering，简称TFP
26、W）在高分辨心电图（HRECG）中的应用。,2022/11/13,88,1.设一共有N次观测样本： 其中S(t)是周期确定的心电信号；Wi(t)是第i次记录时的噪声，包括肌电、测量仪器噪声等，假设每次记录的噪声之间互不相关；xi(t)是观测信号；信号和噪声相互独立。对每次观测用短时傅立叶变换求时频表示（TFR）：,2022/11/13,89,对N次观测的时频表示（TFR）求平均：样本平均为：样本平均的时频表示（TFR）为：,2022/11/13,90,可以得到一个基于样本平均的简单时-频平面后验维纳滤波器：,2022/11/13,91,.在时频域上对式进行修正，给出更实际的表示：,IF是干扰项；,2022/11/13,92,式中COV表示信号和噪声之间的互协方差，也就是考虑了信号和噪声并非相互独立；IF是干扰项；表示样本平均的噪声功率；表示样本噪声功率的平均。,2022/11/13,93,2022/11/13,94,TFPW的模拟实验结果原信号是线性调频信号，观测信号混有白噪声,2022/11/13,95,TFPW的模拟实验结果原信号是正弦波，观测信号混有白噪声,2022/11/13,96,TFPW滤波中由于有二次TFR中的相关噪声以及IF项，滤波器可能包含虚部，也就是包含信号的相位信息，直接在时频平面上考虑相位问题还需要进一步研究。,
