《《维纳滤波的使用》PPT课件.ppt》由会员分享，可在线阅读，更多相关《《维纳滤波的使用》PPT课件.ppt（79页珍藏版）》请在三一办公上搜索。
1、第五章维纳滤波（Wiener Filtering）,王 静,生物医学工程系,本章内容：,5.1 维纳滤波器的时域解5.2 维纳预测器5.3 维纳滤波器的应用,引言,（1）维纳,维纳于1894年生于美国密苏里州哥伦比亚市的一个犹太人的家庭中。他18岁就获得哈佛大学数学和哲学两个博士学位，是信息论的前驱和控制论的奠基人。,维纳：,开创维纳信息论 他从带直流电流或者至少可看作直流电流的电路出发来研究信息论，将统计方法引入通讯工程，奠定了信息论的理论基础。创立控制论 1947年10月，维纳写出了划时代的著作控制论，1948年出版后，立即风行世界。它揭示了机器中的通信和控制机能与人的神经、感觉机能的共同
2、规律；为现代科学技术研究提供了崭新的科学方法。,引言,在前面章节中已经介绍了噪声中确定性信号的检测和估计；随机性是生物医学信号的特点之一，因此，讨论噪声中随机信号的估计具有现实意义。维纳滤波器正是解决该问题。维纳滤波器的限制：要求被估计的随机信号是平稳的。,（2）维纳滤波技术：从噪声中提取有 用的平稳随机信号。,基础知识点回顾,1.卷积运算2.相关运算3.信号与系统,设有一个线性系统，它的单位脉冲响应是h(n)，当输入一个观测到的随机信号x(n)，简称观测值，且该信号包含噪声w(n)和有用信号s(n)，也即 x(n)=s(n)+w(n)则输出y(n)为 y(n)=x(n)*h(n)=,我们希望
3、输出得到的y(n)与有用信号s(n)尽量接近，因此称y(n)为s(n)的估计值，用来表示y(n)，我们就有了维纳滤波器的系统框图这个系统的单位脉冲响应也称为对于s(n)的一种估计器。,维纳滤波技术可应用于以下三个方面：,滤波：用当前的和过去的观测值来估计当前的信 号，称为滤波；预测：用过去的观测值来估计当前的或将来的信号，称为预测；平滑或内插：用过去的观测值来估计过去的信号，称为平滑或者内插。,系统框图中估计到的信号和我们期望得到的有用信号s(n)不可能完全相同，这里用e(n)来表示真值和估计值之间的误差,显然e(n)是随机变量，维纳滤波和卡尔曼滤波的误差准则就是最小均方误差准则,.维纳滤波器
4、的时域解（Time domain solution of the Wiener filter）,设计维纳滤波器的过程就是寻求在最小均方误差下滤波器的单位脉冲响应 h(n)或传递函数H(z)的表达式，其实质就是解维纳霍夫（WienerHopf）方程。,5.1.1 因果的维纳滤波器，即维纳-霍夫方程,5.1.2 有限脉冲响应法求解维纳-霍夫方程,5.1.3 预白化法求解维纳-霍夫方程,因果的维纳滤波器,设h(n)是物理可实现的，也即是因果序列：h(n)=0,当n0,要使得均方误差最小，则将上式对各h(m),m0，1，求偏导，并且等于零，得：(5-7)即(5-8)(5-9),从维纳霍夫方程中解出的h
5、就是最小均方误差下的最佳h，hopt(n)。求到hopt(n)，这时的均方误差为最小：,有限脉冲响应法求解维纳霍夫方程,设h(n)是一个因果序列且可以用有限长（N点长）的序列去逼进它，则式(5-5)(5-10)分别发生变化：(5-11)(5-12)(5-13),(5-14)(5-15)于是得到N个线性方程：,写成矩阵形式有：(5-16)简化形式：RxxH=Rxs(5-17)式中，Hh(0)h(1)h(N-1)是待求的单位脉冲响应:,RxxH=Rxs 只要Rxx是非奇异的，就可以求到H：H=Rxx1Rxs(5-18)求得hopt(n)后，这时的均方误差为最小：,若信号与噪声互不相关，即，Rsw(
6、m)=Rws(m)=0 则有Rxs(m)=Ex(n)s(n+m)=E(s(n)+w(n)s(n+m)=Es(n)s(n+m)+w(n)s(n+m)=Rss(m)Rxx(m)=E(s(n)+w(n)(s(n+m)+w(n+m)=Rss(m)+Rww(m),则式（515）和式（519）化为：(5-20)(5-21)【例5-1】如图，x(n)=s(n)+w(n),且s(n)与w(n)统计独立，其中s(n)的自相关序列,w(n)是方差为1的单位白噪声，试设计一个N2的维纳滤波器来估计s(n)，并求最小均方误差。,解：依题意，已知信号的自相关和噪声的自相关为：代入式（520）得:解得：h(0)0.451
7、，h(1)0.165。将上述结果代入式（521），求得最小均方误差：,5.1.1 因果的维纳滤波器，即维纳-霍夫方程,5.1.2 有限脉冲响应法求解维纳-霍夫方程,5.1.3 预白化法求解维纳-霍夫方程,预白化法求解维纳霍夫方程,随机信号都可以看成是由一白色噪声w1(n)激励一个物理可实现的系统或模型的响应，如图5.2所示 图5.2 s(n)的信号模型,由于x(n)=s(n)+w(n)，在图5.2的基础上给出x(n)的信号模型，图5.3所示。把这两个模型合并最后得到维纳滤波器的信号模型，图5.4所示，其中传递函数用B（z）表示。图5.3 x(n)的信号模型,维纳滤波器的输入信号模型,白噪声的自
8、相关函数为它的z变换就等于。图5.2中输出信号的自相关函数为，根据卷积性质有 令l=r-k，上式,令 代入上式得(5-22)对式（522）进行Z变换得到系统函数和相关函数的z变换之间的关系：(5-23)同样，对图5.4进行z变换得(5-24),图5.4中利用卷积性质还可以找到互相关函数之间的关系：两边z变换得到(5-25),如果已知观测信号的自相关函数，求它的z变换，然后找到该函数的成对零点、极点，取其中在单位圆内的那一半零点、极点构成B(z),另外在单位圆外的零、极点构成B(z-1)，这样就保证了B(z)是因果的，并且是最小相位系统。从图5.4可得(5-26),由于系统函数B(z)的零点和极
9、点都在单位圆内，即是一个物理可实现的最小相位系统，则1/B(z)也是一个物理可实现的最小相移网络函数。我们就可以利用式（525）对x(n)进行白化，即把x(n)当作输入，w1(n)当作输出，1/B(z)系统传递函数。,H(z)=G(z)/B(z),白化法求解维纳霍夫方程步骤如下：）对观测信号x(n)的自相关函数Rxx(m)求z变换得到 Rxx(z)）利用等式找到最小相位系统B(z)）利用均方误差最小原则求解因果的G（z）H(z)=G(z)/B(z),即得到维纳霍夫方程的系统函数解,步骤3的求解过程 按图5.5（b）有(5-28)均方误差为,由于代入上式，并且进行配方得(5-29)均方误差最小也
10、就是上式的中间一项最小，所以(5-30),注意，这里的g(m)是因果的。对该式求z变换,得到(5-31)表示 对求单边z变换。所以维纳霍夫方程的系统函数解表示为,由式(5-32)因果的维纳滤波器的最小均方误差为：(5-33),利用帕塞伐尔定理，上式可用z域来表示(5-34)例52。见书,5.2.1 因果的维纳预测器5.2.2 纯预测器（N步）5.2.3 一步线性预测器,5.2 维纳预测器,因果的维纳预测器,图5.6就是维纳预测器的模型，N0，yd(n)是希望得到的输出，而y(n)表示实际的估计值。图5.6维纳预测器,设h(n)是物理可实现的，也即是因果序列：h(n)=0,当n0，则有(5-35
11、)(5-36)要使得均方误差最小，则将上式对各h(m),m0，1，求偏导，并且等于零，得：,(5-37)即用相关函数R来表达上式：(5-39)由yd(n)=s(n+N)，则,z变换得(5-40)因果的预测器的传递函数为：(5-41)最小均方误差为(5-42),利用帕塞伐尔定理，上式可用z域来表示(5-43),例5-3 已知图5.6中，x(n)=s(n)+w(n),且s(n)与w(n)统计独立，其中s(n)的自相关序列为，w(n)是方差为1的单位白噪声，试设计一个物理可实现的维纳预测器估s(n+1)，并求最小均方误差。解：依题意已知,求z变换 由于，容易找到最小相位系统和白噪声方差,由式（541
12、），N1，对括号里面求z反变换，注意括号内的收敛域为,取因果部分，也就是第一项，所以 把上式写成差分方程形式有：,最小均方误差为：,纯预测器（N步）,纯预测器指的是w(n)0的情况下，对s(n+N)的预测。如图5.7所示。图5.7 N步纯预测器 这时，用白化法来求解预测器的系统函数。,因为，从而有(5-44)将上式代入式（541）、（543）得(5-45),假设B（z）是b（n）的z变换，且b（n）是实序列，则上式可以利用帕塞伐尔定理进一步化简：上式说明最小均方误差随着N的增加而增加，也即预测距离越远误差越大。,例5-4 已知图5.7中，x(n)=s(n)其中s(n)的自相关序列为，试设计一个
13、物理可实现的维纳预测器来估计，并求最小均方误差。解：依题意，已知，则,因为 容易找到最小相位系统和白噪声方差利用式（545）,因为 只取的部分，有代入 得:最小均方误差为：,它说明当N越大，误差越大，当N0时，没有误差 图5.8 例题5-3的纯预测模型,一步线性预测器,对于纯预测问题，有 然而预测的问题常常是要求在过去的p个观测值的基础上来预测当前值，也就是 这就是一步线性预测公式,一步线性预测公式，常常用下列符合表示(5-48)式中p为阶数，。预测的均方误差为(5-49),要使得均方误差最小，将上式右边对 求偏导并且等于零，得到p个等式：(5-50)最小均方误差：(5-51)式（550）就是
14、YuleWalker（Y-W）方程.,例5-5 已知图5.7中，x(n)=s(n),其中s(n)的自相关序列为，试设计一个p2的可实现的一步线性预测器，并求最小均方误差。,解：利用Y-W方程可以列出2个方程式,解得：结果和例（5-4）N=1时一致,维纳滤波器的应用（Application of Wiener filter）,要设计维纳滤波器必须知道观测信号和估计信号之间的相关函数，即先验知识。如果我们不知道它们之间的相关函数，就必须先对它们的统计特性做估计，然后才能设计出维纳滤波器，这样设计出的滤波器被称为“后验维纳滤波器”。,在生物医学信号处理中比较典型的应用就是关于诱发脑电信号的提取。大脑
15、诱发电位（Evoked Potential，EP）指在外界刺激下，从头皮上记录到的特异电位，它反映了外周感觉神经、感觉通路及中枢神经系统中相关结构在特定刺激情况下的状态反应。在神经学研究以及临床诊断、手术监护中有重要意义。EP信号十分微弱，一般都淹没在自发脑电（EEG）之中，从EEG背景中提取诱发电位一直是个难题：EP的幅度比自发脑电低一个数量级，无法从一次观察中直接得到；EP的频谱与自发脑电频谱完全重迭，使得频率滤波失效；在统计上EP是非平稳的、时变的脑诱发电位。通过多次刺激得到的脑电信号进行叠加来提取EP，这是现今最为广泛使用的EP提取方法,为了解决诱发电位提取问题，研究者利用维纳滤波来提
16、高信噪比，先后有Walter、Doyle、Weerd等对维纳滤波方法进行了改进，如时变维纳滤波。在频域应用后验维纳滤波的核心就是由各次观察信号中分解出EP信号的谱估计和噪声的谱估计，通过设计出的滤波器来提高信噪比,从而来提取出诱发电位。下面用实际例子说明如何用维纳滤波方法进行诱发电位的提取。（见参考文献）,下面将介绍时频平面的维纳滤波（timefrequency plane wiener filtering，简称TFPW）在高分辨心电图（HRECG）中的应用。,设一共有N次观测样本：其中s(t)是周期确定的心电信号；是第i次记录时的噪声，包括肌电、测量仪器噪声等，假设每次记录的噪声之间互不相关
17、；是观测信号；信号和噪声相互独立。,对每次观测用短时傅立叶变换求时频表示（TFR）：对N次观测的时频表示（TFR）求平均：样本平均为：,样本平均的时频表示（TFR）为：从上式可以得到一个基于样本平均的简单时频平面后验维纳滤波器：,.在时频域上对式（1）（2）进行修正，给出更实际的表示：式中COV表示信号和噪声之间的方差，也就是考虑了信号和噪声并非相互独立；IF是干扰项；表示样本平均的噪声功率；表示样本噪声功率的平均。,.TFPW的计算过程如图5.9所示。,.TFPW的模拟实验结果如图5.13所示。,图5.13（上图）原信号是两个正弦波，观测信号混有白噪声（下图）原信号是线性调频信号，观测信号混有白噪声.TFPW滤波中由于有二次TFR中的相关噪声以及IF项，滤波器可能包含虚部，也就是包含信号的相位信息，直接在时频平面上考虑相位问题还需要进一步研究。,thank you for your attention!,
