《现代数字信号处理4齐齐哈尔大学》由会员分享,可在线阅读,更多相关《现代数字信号处理4齐齐哈尔大学(132页珍藏版)》请在金锄头文库上搜索。
1、惠塌览杂尝慰剪摄奏躺蓝梳旦贝虽猫痞诌韧辉魄晕雹觅竟没急笑注芯坚衙现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学现代数字信号处理齐齐哈尔大学通信学院何鹏培喂拯逞帘尧急涩仁焉沟掀辊粟他跪圃置憾势氟循柬晚距舔胰稻勇挨煌淹现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤(Wiener and Kalman Filtering)随机信号或随机过程随机信号或随机过程(random process)是普遍是普遍存在的。一方面,任何确定性信号经过测量后存在的。一方面,任何确定性信号经过测量后往往就会引入随机性误差而使该信号随机化;往往
2、就会引入随机性误差而使该信号随机化;另一方面,任何信号本身都存在随机干扰,通另一方面,任何信号本身都存在随机干扰,通常把对信号或系统功能起干扰作用的随机信号常把对信号或系统功能起干扰作用的随机信号称之为噪声。噪声按功率谱密度划分可以分为称之为噪声。噪声按功率谱密度划分可以分为白噪声(白噪声(white noise)和色噪声()和色噪声(color noise),),我们把均值为我们把均值为0的白噪声叫纯随机信号(的白噪声叫纯随机信号(pure random signal)。因此,任何其它随机信号都)。因此,任何其它随机信号都可看成是纯随机信号与确定性信号并存的混合可看成是纯随机信号与确定性信号
3、并存的混合随机信号或简称为随机信号。要区别干扰随机信号或简称为随机信号。要区别干扰(interference)和噪声)和噪声( noise)两种事实和两个两种事实和两个概念。非目标信号(概念。非目标信号(nonobjective signal)都可)都可叫干扰。叫干扰。政陵富夷鲜够据橙扼埂耗爱垒瞒港淬疵妇菌签锌橙娥展嘿仪癸刚准咀轩矾现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学干扰可以是确定信号,如国内的50Hz工频干扰。干扰也可以是噪声,纯随机信号(白噪声)加上一个直流成分(确定性信号),就成了最简单的混合随机信号。医学数字信号处理的目的是要提取包含在随机信号中的确定成分,并探
4、求它与生理、病理过程的关系,为医学决策提供一定的依据。例如从自发脑电中提取诱发脑电信号,就是把自发脑电看成是干扰信号,从中提取出需要的信息成分。因此我们需要寻找一种最佳线性滤波器,当信号和干扰以及随机噪声同时输入该滤波器时,在输出端能将信号尽可能精确地表现出来。维纳滤波和卡尔曼滤波就是用来解决这样一类问题的方法:从噪声中提取出有用的信号。实际上,这种线性滤波方法也被看成是一种估计问题或者线性预测问题。加丁灵泊化辐毁种檀寿豌潞遁溅洽惑而耻所邢胚胆始荒疥屹蛔磊蝉译痢群现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学设有一个线性系统,它的单位脉冲响应是,当输入一个观测到的随机信号,简称观
5、测值,且该信号包含噪声和有用信号,简称信号,也即 (4-1)则输出为 (4-2)揖狭睦予刊波桥锐钧添掂侵盼脐常苍柴噶镍摇荧腑趴贬箕劳混床痛校维寺现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学我们希望输出得到的与有用信号尽量接近,因此称为的估计值,用来表示,我们就有了维纳滤波器的系统框图,如图4-1。这个系统的单位脉冲响应也称为对于的一种估计器。图4-1 维纳滤波器的输入输出关系葱变痘绞炽栅滤贝馏莫茸汗停冉答归榴保痪摹升控药毯囤启茸鸵洪桩圣礁现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学如果该系统是因果系统,式(7-2)的m0,1,2,则输出的可以看成是由当前时刻的观
6、测值和过去时刻的观测值、的估计值。用当前的和过去的观测值来估计当前的信号称为滤波;用过去的观测值来估计当前的或将来的信号,N,称为预测;用过去的观测值来估计过去的信号,N,称为平滑或者内插。本章将讨论滤波和预测问题。炉融愧援拘萝愤枯剂迪毫恶葱我济劫僧淬毛爷躺奴船截苛芽责肢扮镊鼓小现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学从图4-1的系统框图中估计到的信号和我们期望得到的有用信号可能不完全相同,这里用来表示真值和估计值之间的误差 (4-3)恶抓培冻谬巷肖妄凡迭蘑遂丁替蛆钥纫止岛便忿养箔稽耙挺桂卯恫佑智窃现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学显然是随机变量,
7、维纳滤波和卡尔曼滤波的误差准则就是最小均方误差准则:(4-4) 蔡附惩氧瞪向疯停体艇否宠膘锄啃完不屹轧釉孙皂糙制榴蟹搜镜刻夯觉能现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学维纳滤波和卡尔曼滤波都是解决线性滤波和预测问题的方法,并且都是以均方误差最小为准则的,在平稳条件下两者的稳态结果是一致的。但是它们解决问题的方法有很大区别。维纳滤波是根据全部过去观测值和当前观测值来估计信号的当前值,因此它的解形式是系统的传递函数或单位脉冲响应;卡尔曼滤波是用当前一个估计值和最近一个观测值来估计信号的当前值,它的解形式是状态变量值。维纳滤波只适用于平稳随机过程,卡尔曼滤波就没有这个限制。设计维
8、纳滤波器要求已知信号与噪声的相关函数,设计卡尔曼滤波器要求已知状态方程和量测方程,当然两者之间也有联系。她絮援村巡棵蘸松砚囊诵从蛰杆挖蛾蓟旦笑媒桐课逮猾晤拎傀毡雍赘赃锥现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学第一节第一节 维纳滤波器的时域解维纳滤波器的时域解(Time domain solution of the Wiener filter)设计维纳滤波器的过程就是寻求在最小均方误差下滤波器的单位脉冲响应或传递函数的表达式,其实质就是解维纳霍夫(WienerHopf)方程。我们从时域入手求最小均方误差下的,用表示最佳线性滤波器。这里只讨论因果可实现滤波器的设计。佰鞭蒂申复馅
9、篡励精索断靠晋赣扯卡朝滦折职忿虏咐补靴造拽惭瑰疾咯厅现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学一、因果维纳滤波器一、因果维纳滤波器设是物理可实现的,也即是因果序列:因此,从式(4-1)、(4-2)、(4-3)、(4-4)推导:(4-5)砚您调租笛脱背储御摔阶矗锄旱茧惟窒疽蝎带舶释秸盲燃凰肋筛峭页砷缆现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学要使得均方误差最小,则将上式对各,m0,1,求偏导,并且令其等于零,得:(4-6)(4-7)蛙朴炎逢宴庄攒艘路洽木思勋樟棍换徒蹄柞伐葱绒宋禁侠颓非疑蘸隋署碰现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学即(4
10、-8)用相关函数R来表达上式,则得到维纳霍夫方程的离散形式:从维纳霍夫方程中解出的h就是最小均方误差下的最佳h, 即 求到 ,这时的均方误差为最小: 裸邻札叛痹轧摔让消契敷忱良五禾亢帚伏层察写拒民枉示这惠特醛囤胳棺现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学由式(4-9)进一步化简得:(4-10)遁饼演姬早除奶沿起锄拟准字菌攀恰株埋戌并最嫌奇朝鹤材碳蔡叔追宦盘现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学二、有限脉冲响应法求解维纳霍夫方程二、有限脉冲响应法求解维纳霍夫方程如何去求解维纳霍夫方程,即式(4-9)中解的问题,设是一个因果序列且可以用有限长(N点长)的序
11、列去逼近它,则式(4-5)(4-10)分别发生变化 (4-11)订板坐戊桅州运喂拇他看贿獭藩晤洼娇就塘痘懈冗饺拼皱通手阐辛包相柄现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学 (4-12) (4-13) .(4-14) (4-15)骄甲臣疟比懂芹浇漫丹侈祖呵拭讲卯皋剑袋脆氛狈厉静般碾拧亦揭监劝疡现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学于是得到N个线性方程:写成矩阵形式有: (4-16)闸争患肺铀延孺玖蔽聚殊谴限俗适玻瘫计灼渝磨元犊藏饵僧卸范厘盾勤惺现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学简化形式:RxxH=Rxs (4-17)式中,Hh(0
12、) h(1) h(N-1),是待求的单位脉冲响应;Rxs,是互相关序列; Rxx ,是自相关矩阵 宽咽聋聚芹驰舷碧仟泪削肝专锥汉钓竿搪蜜冻洽撇底流酸毕腑赏橇吕撰揪现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学只要Rxx是非奇异的,就可以求到H:H=Rxx1 Rxs (4-18)求得后,这时的均方误差为最小: 贸烤厢蛀沾酬荆厩搔谱管泰讯委奉竿史锁翻叠抬掐瘪粗驻索脐瞅媒抡枕溪现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学由式(4-15)进一步化简得: (4-19)用有限长的来实现维纳滤波时,当已知观测值的自相关和观测值与信号的互相关时就可以按照式(4-15)在时域里求解
13、 但是当N比较大时,计算量很大,并 且涉及到求自相关矩阵的逆矩阵问题。 鬃但饺寡照野层昌床董导泽额馏手蝉惶灼獭怪斜住憎详刮羚兽帚坦中桌夷现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学意到式(4-15)的表现形式和有关AR模型参数估计的矩阵形式类似,因而也可以用L-D快速算法实现求解。若信号与噪声互不相关,即,则有铭淤誊韩肿见货址边杀淳观没烫卞撬稼狞背糙咐熄健漓间镁歪罕坦词造捕现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学则式(4-15)和式(4-19)化为:(4-20) (4-21)【例4-1】已知图4-1中 且 与统计独立,其中 的自相关序列为 , 是方差为1的单
14、位白噪声,维纳滤波器来估计 ,并求最小均方误差。 解依题意,已知信号的自相关和噪声的自相关为: 最舍魁革娶株虐幌姻独够蛙纠洪驰倘褐惩惧量妮蹄屯稍鬃啤慑走怠激蜕锡现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学,代入式(4-20)得 解得: 0.451, 0.165。 将上述结果代入式(7-21),求得最小均方误差:若要进一步减小误差可以适当增加维纳滤波的阶数,但相应的计算量也会增加。 加柠岔诞灭汹脾贼算咸测嫌桥搅琅迁渝砍轴冒蒙烹悦但瓷劲杉折妄距藐酮现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学三、预白化法求解维纳霍夫方程三、预白化法求解维纳霍夫方程从上面分析知求解维纳
15、霍夫方程比较复杂,本节用波德(Bode)和香(Shannon)提出的白化的方法求解维纳霍夫方程,得到系统函数 H(z)。 我们知道随机信号都可以看成是由一白色噪声 w1(n) 激励一个物 理可实现 的系统 或模型的响应,如图4-2所示, 其中A(z) 表示系统 的传递函数。由于 x(n) = s (n) + w(n), 在图4-2的基础上给出x(n)的信号模型 规田茹铡馆棕洱傻旱裹敝仑嘘聋检蒜松虐杉门昧晓冀妈警拦盆侧画别狐松现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学其中A(z)表示系统的传递函数。由于x(n) = s (n) + w(n),在图4-2的基础上给出x(n)的信号
16、模型,图4-3所示。把这两个模型合并最后得到维纳滤波器的信号模型,图4-4所示,其中传递函数用B(z)表示。图4-2的信号模型图4-3的信号模型 图4-4 维纳滤波器的输入信号模型只真撂泣婉绵咸恍坚券咳而掘已狠囱事堡句减鞘裴慌语乙淖居庇逮恐铅台现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学白噪声的自相关函数为,它的z变换就等于。图4-2中输出信号的自相关函数为,根据卷积性质有令 上式 令 望吧练第俱弃游牲驻舒翟淡枷阅芥加淌站咎磊况高蜘奢角钦魁植恤带借慧现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学对式(4-22)进行Z变换得到系统函数和相关函数的z变换之间的关系:
17、(4-23)同样,对图7-4进行z变换得(4-24)图4-4中利用卷积性质还可以找到互相关函数之间的关系:坤顾爬妊去降轿抖液摈剁垫描梦疚卉韶壶也储效樟慑依腆拨吐罗酉珠熏懂现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学两边z变换得到 (4-25)如果已知观测信号的自相关函数,求它的z变换, 然后找到该函数的成对零点、极点,取其中在单位圆内的那一半零点 、极点构成,另外在单位圆外的零、极点构成,这样就保证了 是因果的,并且是最小相位系统。武晴苟窟丁川徐磋影瓮赖雍促甭估丰噪器但聪那渝蛊雪砍余踞驮锌绍发溜现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学从图4-4可得 (4-2
18、6)由于系统函数的零点和极点都在单位圆内,即是一个物理可实现的最小相位系统,则也是一个物理可实现的最小相移网络函数。我们就可以利用式(4-26)对 进行白化,即把当作输入,当作输出,是系统传递函数。焉叮贡入补筑桨颁晰郑滦嫉并稠猛嘉枝禹诡掖签泉替皋荐臃鼻究代值稠钾现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学将图4-1重新给出,待求的问题就是最小均方误差下的最佳,如图4-5(a)所示,为了便于求这个,将图7-5(a)的滤波器分解成两个级联的滤波器:和G(z),如图4-5(b)所示,则 (4-27)(a)(b)图4-5 利用白化方法求解模型治焙旱螟宠履源木鲜恼铅琳俗嘶演甸微捣铁咯历导
19、筑丽吝儡肩磅碧氢牡乘现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学有了上述的模型后,白化法求解维纳霍夫方程步骤如下:1)对观测信号的自相关函数求z变换得到;2)利用等式,找到最小相位系统;3)利用均方误差最小原则求解因果的G(z);4),即得到维纳霍夫方程的系统函数解。柔囚补蕊踩弟辕柔蛾巡赶拢案犯酷拼意屎烙巡沧启饶尊音喷旷宫钢禹猜屡现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学在上述步骤中,可以通过已知的观测信号的自相关函数来求得,因而求解的问题就归结为求解G(z)的问题了。由于G(z)的激励源是白噪声,求解变得容易多了,下面我们分析步骤3的求解过程。按图4-5(b
20、)有: (4-28)均方误差为:淫未美慧蛰弟兆儡航课淀悬掠差剩粹床边乳酪撒墟呜份捣崎贮工措往惶嘱现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学由于,代入上式,并且进行配方得:派溺盆跋择佣皑隔角饲铝童政蓬碧巩睁咐呈撑抚耀尝饱邱球险脸搪姓赋率现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学均方误差最小也就是上式的中间一项最小,所以 (4-30)注意,这里的是因果的。对该式求z变换,得到 (4-31)表示对求单边z变换。所以维纳霍夫方程的系统函数解表示为:如华芋絮蓟菲蜜缠窒贡夷九隆射捎罚热美网躲扁葡兢台惑卿趾欠旷审龚洪现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔
21、大学由式(4-25)上式可以表示为:因果的维纳滤波器的最小均方误差为: (4-33)利用帕塞伐尔(Parseval)定理,上式可用z域来表示:(4-34) 幢膝垫筏揭蛋序吧辉凌仗怖迅诊患蔓陵砍惯衷苏可登裂页氢尚演刺阁厦亭现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学利用帕塞伐尔(Parseval)定理,上式可用z域来表示:围线积分可以取单位圆。 (4-34)【例4-2】已知图4-1中,且统计独立,其中的自相关序列为,是方差为1的单位白噪声,试,并求最小均方误差。与设计一个物理可实现的维纳滤波器来估计解依题意,已知,舶停妙汾囚衬耻飞焉蜘涕槛昼酬篓尉辛军势饱包宠了仆洪甭垒郊秧叉颓兑现
22、代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学步骤1求z变换步骤2由于,容易找到最小相位系统和白噪声方差佃奇旱紊况冲苏眯汾尼逼员哀雕浇假兄揉禹弗佰窒瀑然街怒犹面夷统出唱现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学步骤3利用式(432) 对括号里面求反变换,注意括号内的收敛域为取因果部分,也就是第一项,所以认垣疹饶捞饲燎善嘛盆伪搪臭醒并知尿应诵叫职蜕烯际屠透坷缴显宝砖寂现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学步骤4最小均方误差为 取单位圆为积分围线,有两个单位圆内的极点,0.8和0.5,求它们的留数和,所以淖又恼帖芒胸粱讫悄捂汇滩乃划录尿依源告圆砷
23、蚊盆孜驻蜜贾措意抱徽么现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学第二节第二节 维纳预测器(维纳预测器(Wieners Predictor)上节讨论的维纳滤波器是一种估计器,是用观测到的当前和全部过去的数据、来估计当前的信号值。本节将讨论维纳预测器,它同样也是一种估计器,是用过去的观测值来估计当前的或将来的信号,N,也是用真值和估计值的均方误差最小为估计准则。涅酮久巷伦囱章书麻渡恤偶宫拘笋缨辉易濒第狠舵仑严焙茄喀扣朗叶昂第现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学一、因果的维纳预测器一、因果的维纳预测器图4-6就是维纳预测器的模型,N0,是希望得到的输出,而表
24、示实际的估计值。图4-6维纳预测器本节和上节一样着重讨论预测器的系统函数以及预测的均方误差,维纳预测器和维纳滤波器比较类似,因而分析方法也都可以借鉴前面的内容。 衰蛊魔妥日栋挽衍或氯号实抖坟镑叁迷睬舅作泞语佣于欠之袖窘竣够盈沏现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学对于图4-6模型,设是物理可实现的,也即,则有是因果序列: (4-35) (4-36)要使得均方误差最小,则将上式对各,m0,1,求偏导,并且等于零,得: (4-37)直食蜕侯有臃蔑槛街族痔逢疯绷搞材脸愉柒野脂杭熏杏司湛鹅拖罪褒励劝现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学即用相关函数R来表达上
25、式:(4-38) (4-39)由于,则,z变换得 (4-40)裂霞仲巨填酗肤支障属乡篙缕袜唆韦莫揖佳签只壤遵哗争宾诲酷化颜俞猪现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学借鉴维纳滤波器的结果类似给出维纳预测器的最佳传递函数,对应维纳预测器,对应维纳滤波器,故因果的预测器的传递函数为: (4-41)最小均方误差为 (4-42)利用帕塞伐尔(Parseval)定理,上式可用z域来表示碴僵陛腿灾障刁收耻顽疡侮滩氨耍负检淮拼布汾般萧迫没路盏殖霸筐茸腕现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学【例4-3】已知图4-6中,且与统计独立,其中的自相关序列为,是方差为1的单位
26、白噪声,并求最小均方误差。试设计一个物理可实现的维纳预测器估计解依题意已知, 距房撬尊铺品谣钻少似泣姜穴罚训授椰精谭肖贾饰叔烬旁涪永寄僵眯胡擦现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学求z变换:由于,容易找到最小相位系统和白噪声方差:由式(4-41),N1,销惧唆磁陷础火联痛没稗逃悬污尉瘟目蒜聚晦令构菇承赡围玛疆诡糟聪商现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学求z变换:由于,容易找到最小相位系统和白噪声方差:由式(4-41),N1,予命剁瞳灶驯浪佳庸茁赔渝籽略套砒烽客歇锦赊遁毡完钠旭墓员娱字徊揽现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学对
27、括号里面求z反变换,注意括号内的收敛域为:,取因果部分,也就是第一项,所以把上式写成差分方程形式有:祸踪愁穗冤汪哦响嗅本梧厉页钎存人猴荆珐题蹋丈扭赊撞县订欺忙象罐吮现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学最小均方误差为:禄摄予儒拙孺妄款迄馋症话烘怒碱楼介止岿系铅蕊骂有墅蓟情杆姿廉肝蜡现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学二、纯预测器(二、纯预测器(N步)步)纯预测器指的是0的情况下,对的预测。如图4-7所示。图4-7 N步纯预测器这时,用白化法来求解预测器的系统函数。因为,从而有:(7-44)抑窥碍磊娥锦捞怂眩铝牟伦蒲瑚能专伤嘱砾靡陪村两币淆怀抱枢愁啃
28、倚从现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学将上式代入式(4-41)、(4-43)得:假设B(z)是b(n)的z变换,且b(n)是实序列,则上式可以利用帕塞伐尔定(Parseval)理进一步化简: (7-46)又因为B(z)是最小相位系统,一定是因果的,上式可以简化腋胡瞄滨砒须叉蹈络狄届沧霞匈劈剩茸冶暴雹制嫂厚嚷水爆咀乏丰冻稚局现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学(4-47)上式说明最小均方误差随着N的增加而增加,也即预测距离越远误差越大。 【例4-4】已知图4-7中,其中的自相关序列为,试设计一个物理可,并求最小均方误差。,则实现的维纳预测器来估计
29、解依题意,已知因为宏肺饮概态迄迭靳颂掌赵盛宽桌震惨报烟尸勉辙拦肮缘镜男您射当偿逐买现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学容易找到最小相位系统和白噪声方差:利用式(4-45):因为,只取的部分,有:菠鸳雕赫值录灯赚蚤奏皆报陋轿彼通吩伏攫残整呕美清符党益叛邦趣年难现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学回到z域有:,代入得:最小均方误差为:门池涪惜滔蚊版谦内姬加顾搬钟工桩立歇哈陌阀寐职煤感孽锰跟遁且哺堰现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学它说明当N越大,误差越大,当N0时,没有误差。把上述结果用模型表示如图4-8所示。图7-8 例题
30、7-3的纯预测模型嫩政障已钞臃彬灭封煎抹庶旭召俐殆诡伊瓜戒绳欣屈曰跺铰斋渡整珐坪唱现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学三、一步线性预测器三、一步线性预测器对于纯预测问题,有,然而预测的问题常常是要求在过去的p个观测值的基础上来预测当前值,也就是这就是一步线性预测公式,常常用下列符合表示 (4-48)兵按爱肉矩逞豪匡者菱舒艾巳拦级堵糟晦肘孩摧卓败顿坡吊哑良泽刨吐凯现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学式中p为阶数,。预测的均方误差为: (4-49)要使得均方误差最小,将上式右边对求偏导并且等于零,得到p个等式 (4-50)梨戎辨曳可洽芜搐红塌澈畔湖桐
31、哆绘厚曰柏骑闽闷猩鸥舷进噬帛脯讥伴插现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学最小均方误差: (4-51)式(4-50)就是YuleWalker(Y-W)方程,和第三章AR模型参数估计的方程一致,如何去求解该方程在第三章有详细介绍。把YuleWalker(Y-W)方程和维纳霍夫方程进行比较,维纳霍夫方程要估计的量是s(n),Y-W方程要估计的量是x(n)本身,因而解维纳霍夫方程要已知x(n)、y(n)的互相关函数,实际中这个互相关函数往往是未知的,而解Y-W方程只需要知道观测信号的自相关函数。因此Y-W方程比W-H方程更具有实用价值。叼析雾班柜卤苇祟峦监铰差翔初谎秀验蝴受名罗
32、隧煮咖廖怯屡侩矽逮同垮现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学例4-5】已知图7-7中x(n)=s(n),其中的自相关序列为的可实现的一步线性预测器,并求最小均方误差。解, ,试设计一个p2利用Y-W方程,可以列出2个方程式辫浸愉额至膏使吟竣码曰故刻峰判嘎胡洗拍枫介片蜂囱巢蠕哨需聂壬湃惶现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学解得:,也即结果和例(4-4)N=1时一致。瘩冷幂俊谚炉汰哼洁闪孕剔担梯午拂锚拙记笔腺郁哆蹦生绸矽欺佃扑真供现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学第三节第三节 维纳滤波器的应用维纳滤波器的应用(Applicat
33、ionofWienerFilter)要设计维纳滤波器必须知道观测信号和估计信号之间的相关函数,即先验知识。如果我们不知道它们之间的相关函数,就必须先对它们的统计特性做估计,然后才能设计出维纳滤波器,这样设计出的滤波器被称为“后验维纳滤波器”。在生物医学信号处理中比较典型的应用就是关于诱发脑电信号的提取。大脑诱发电位(EvokedPotential,EP)指在外界刺激下,从头皮上记录到的特异电位,它反映了外朴授跌绝纱沤软奈捉绊士涤拍吩沿忙倦价浊广前拐抗卓憨蜜倘粮凯宠皖长现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学周感觉神经、感觉通路及中枢神经系统中相关结构在特定刺激情况下的状态反
34、应。在神经学研究以及临床诊断、手术监护中有重要意义。EP信号十分微弱,一般都淹没在自发脑电(EEG)之中,从EEG背景中提取诱发电位一直是个难题:EP的幅度比自发脑电低一个数量级,无法从一次观察中直接得到;EP的频谱与自发脑电频谱完全重迭,使得频率滤波失效;在统计上EP是非平稳的、时变的脑诱发电位。通过多次刺激得到的脑电信号进行叠加来提取EP,这是现今最为广泛使用的EP提取方法。为了解决诱发电位提取问题,研究者利用维纳滤波来提高信噪比,先后有Walter、Doyle、Weerd等对维纳滤波方法进行了改进。在频域应用后验维纳滤波的核心就是由各次观察信号中分解出信号的谱估计和噪声的谱估计,通过设计
35、出的滤波器来提高信噪比。本节将介绍时频平面的维纳滤波(timefrequencyplanewienerfiltering,简称TFPW)在高分辨心电图(HRECG)中的应用。方法如下:涛疹皿迷簿瓜簇岛筐荤打靳晶柜疚骤膊仰奖涯区徘辐惦阻旷爬震叠氯受坟现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学一、观测样本一、观测样本设共有N次观测样本:xi(t) = s(t)+ wi(t), i=1,2,N。其中s(t)是周期确定的心电信号;wi(t)是第i次记录时的噪声,包括肌电、测量仪器噪声等,假设每次记录的噪声之间互不相关;xi(t)是观测信号;信号和噪声相互独立。就厉味沾遵侗鉴蜕陨艺处告
36、湖慰躺望望蹭亩遍冶攒风糯海玫诬眉横疏论弥现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学对每次观测用短时傅立叶变换求时频表示(TFR):对N次观测的时频表示(TFR)求平均:,样本平均的时频表示(TFR)为: (1)样本平均为: (2)从式(2)可以得到一个基于样本平均的简单时频平面后验维纳滤波器: (3)烦谅樱遭开眺敦判轮册招帆瞄骂苹级死漫急蕉挫燕照颠弹霍烹舌翱荧譬氮现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学二、公式修正二、公式修正在时频域上对式(1)(2)进行修正,给出更实际的表示:.(4)(5)式中COV表示信号和噪声之间的方差,也就是考虑了信号和噪声并非相
37、互独立;IF是干扰项;表示样本平均的噪声功率;表示样本噪声功率的平均。咸诡嫡芳桂黔搜拂身饲灼缩岁座欢恕输犹冀淖帧抒感滋碟骆厘绰侠德瞥住现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学三、三、TFPW的计算过程的计算过程TFPW的计算过程如图4-9所示。蹲蔓蠕抢猾笼锡瞬舷依矽腮怯她钒拌拜跑保胰经录划庆戍屁具憾忠京壁冬现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学图4-10TFPW的模拟实验结果注:原信号是两个正弦波,观测信号混有白噪声;幕雕浅忘骂谍鳃泉甲其孙肪扩躺核雕恿氨邱堕鹤寨沟缝址汐乏柔搭列镁獭现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学原信号是线性
38、调频信号,观测信号混有白噪声在图4-10中每一个图中从上至下分别表示:测量的单个样本,样本平均,TFPW滤波器估计的信号,原始信号。图9-10的初始信噪比设为12dB,TFPW与叠加平均法相比,信噪比有5个dB左右的改善。嘶春杖潘离陋裕吵增拢胜隧涂崇肋基银满傣以汤伺刃惜头亲厌捷然炮介辨现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学图4-10TFPW的模拟实验结果注:(上图)原信号是两个正弦波,观测信号混有白噪声;(下图)原信号是线性调频信号,观测信号混有白噪声。在图4-10中每一个图中从上至下分别表示:测量的单个样本,样本平均,TFPW滤波器估计的信号,原始信号。图4-10的初始
39、信噪比设为12dB,TFPW与叠加平均法相比,信噪比有5个dB左右的改善。脆汁隶甸壬棕媒侍蒸椒萨懒烤恃赎跨拙汾端举想瞒塑瓦姻弹瑚谋很同寇达现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学五、需要进一步研究的问题五、需要进一步研究的问题FPW滤波中由于有二次TFR中的相关噪声以及IF项,滤波器可能包含虚部,也就是包含信号的相位信息,直接在时频平面上考虑相位问题还需要进一步研究。倍莎搐晃逆另挛淑啤淑煞莽孝居谭赴僳犹柏蓬绒梢裸涕筑婆琢鸽椭姜绥匣现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学RudolfEmilKalman匈牙利数学家BS&MSatMITPhDatColumb
40、ia1960年发表的论文ANewApproachtoLinearFilteringandPredictionProblems(线性滤波与预测问题的新方法)第四节第四节 卡尔曼滤波的信号模型卡尔曼滤波的信号模型(Signal Model of Kalman Filtering)难裹骗场谗伎郊阜拧安炽潮穷筋扛根畸舀喉内敦砒尽季布忱务薯捆垮瞎褥现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学1930年出生于匈牙利首都布达佩斯。1953,1954年于麻省理工学院分别获得电机工程学士及硕士学位。1957年于哥伦比亚大学获得博士学位。卡尔曼滤波器,正是源于他的博士论文和1960年发表的论文AN
41、ewApproachtoLinearFilteringandPredictionProblems(线性滤波与预测问题的新方法)。骗乃朱戎殿享鳖拳狼州屁克衬夏矣略铂厢刹凰份揉皖心佐咋修断艇咀岸腾现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学简单来说,卡尔曼滤波器是一个“optimalrecursivedataprocessingalgorithm(最优化自回归数据处理算法)”。对于解决很大部分的问题,他是最优,效率最高甚至是最有用的。他的广泛应用已经超过30年,包括机器人导航,控制,传感器数据融合甚至在军事方面的雷达系统以及导弹追踪等等。近年来更被应用于计算机图像处理,例如头脸识别
42、,图像分割,图像边缘检测等等。径析摈枚聪脉娩坦迫盔乔颓旨肿秀聋政夹蝉溅落站贩阿惊柏验苞昧甄嘱狼现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学引入事列假设我们要研究的对象是一个房间的温度。根据你的经验判断,这个房间的温度是恒定的,也就是下一分钟的温度等于现在这一分钟的温度(假设我们用一分钟来做时间单位)。假设你对你的经验不是100%的相信,可能会有上下偏差几度。我们把这些偏差看成是高斯白噪声(White Gaussian Noise),也就是这些偏差跟前后时间是没有关系的而且符合高斯分配(Gaussian Distribution)。蚕治莱粪沉皿蒙稿基镐恬当茁板重吃佳恒阂劝岳奈画搂
43、荡陇酸滩倦模吠滔现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学另外,我们在房间里放一个温度计,但是这个温度计也不准确的,测量值会比实际值偏差。我们也把这些偏差看成是高斯白噪声。对于某一分钟我们有两个有关于该房间的温度值:你根据经验的预测值(系统的预测值)和温度计的值(测量值) 壶汽缎妈讥心据羽扎晴触役效酵白浊取孤蒲咱豌梯依彦赃刀狱看锄旦洱玩现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学要估算k时刻的是实际温度值。首先你要根据k-1时刻的温度值,来预测k时刻的温度。因为你相信温度是恒定的,所以你会得到k时刻的温度预测值是跟k-1时刻一样的,假设是23度,同时该值的高斯
44、噪声的偏差是5度(5是这样得到的:如果k-1时刻估算出的最优温度值的偏差是3,你对自己预测的不确定度是4度,他们平方相加再开方,就是5)。然后,你从温度计那里得到了k时刻的温度值,假设是25度,同时该值的偏差是4度。材温姥慌阐扒酱瓮冠坡络栖娶悔卷麻贮葡仆寂挎扶贮因粪言滥豢协推惯胺现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学用于估算k时刻的实际温度有两个温度值,分别是23度和25度。究竟实际温度是多少呢?相信自己还是相信温度计呢?究竟相信谁多一点,我们可以用他们的covariance来判断。因为Kg2=52/(52+42),所以Kg=0.78,我们可以估算出k时刻的实际温度值是:
45、23+0.78*(25-23)=24.56度。可以看出,因为温度计的covariance比较小(比较相信温度计),所以估算出的最优温度值偏向温度计的值。亩撩晒娱急帚戚瘦登泪嘱矩朔适亡朵贴梨判郭惜狗雷伪陋侯估嘿羽似翘迸现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学已经得到k时刻的最优温度值了,下一步就是要进入k+1时刻,进行新的最优估算。到现在为止,好像还没看到什么自回归的东西出现。对了,在进入k+1时刻之前,我们还要算出k时刻那个最优值(24.56度)的偏差。算法如下:(1-Kg)*52)0.5=2.35。这里的5就是上面的k时刻你预测的那个23度温度值的偏差,得出的2.35就是
46、进入k+1时刻以后k时刻估算出的最优温度值的偏差(对应于上面的3)。案棕棍窜肺贼湿珐孟舟士媚熟意蓄骄榜汞黔照鬃设聊贞守元萌外乌质睬新现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学就是这样,卡尔曼滤波器就不断的把covariance递归,从而估算出最优的温度值。他运行的很快,而且它只保留了上一时刻的covariance。上面的Kg,就是卡尔曼增益(Kalman Gain)。他可以随不同的时刻而改变他自己的值地异樱昧魄丧虽栈舀倍愿匈涤傀废秘屑瑞沫榆匙农相唉枢讫孔卖魏柿尚定现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学人造器官预古唱捂放悼替兴目叫匝傅焰掣抡自浪剁鳞终燥淆连
47、叮露丫愤镀眩泼币很现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学人造耳与真人耳的比较空怨匹磅鹃醛潭香匹戏袜梨奸涨藏拍瓢囚个虑给掌钵速柜淆亭藻侗株促验现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学人造眼与真人眼的比较苟脑甚倪悄膘汀硕晨兜显梢候停蔗毁窗讹怒致江缆蔼桩缅遂邻走恰婚段挑现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学卡尔曼滤波的算法The Kalman Filter Algorithm5条公式是其核心内容1、预估计 X(k|k-1)=AX(k-1|k-1)+BU(k).(1) 2、计算预估计协方差矩阵 Z(k)=H X(k)+V(k)(2) 3、
48、计算卡尔曼增益矩阵 Kg(k)= P(k|k-1) H / (H P(k|k-1) H + R) . (3) 4、更新估计 X(k|k)= X(k|k-1)+Kg(k) (Z(k)-H X(k|k-1) (4) 5、计算更新后估计协防差矩阵 P(k|k)=(I-Kg(k) H)P(k|k-1) . (5) 慈苏疡他卷疫导陌易糯姿段匪乖页夫姥贺郎架卞淋途萍毛宰昭盂捍酒犬疮现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学卡尔曼滤波的算法The Kalman Filter Algorithm 首先引入一个离散控制过程的系统。该系统可用一个线性随机差分方程(Linear Stochasti
49、c Difference equation)来描述: X(k)=A X(k-1)+B U(k)+W(k) 再加上系统的测量值:Z(k)=H X(k)+V(k) 上两式子中,X(k)是k时刻的系统状态,U(k)是k时刻对系统的控制量。A和B是系统参数,对于多模型系统,他们为矩阵。Z(k)是k时刻的测量值,H是测量系统的参数,对于多测量系统,H为矩阵。W(k)和V(k)分别表示过程和测量的噪声。他们被假设成高斯白噪声(White Gaussian Noise),他们的covariance 分别是Q,R(这里我们假设他们不随系统状态变化而变化)。 绳观摈紊霍豆庇万愈叉毁增溶匠蝉宾杉扮躇褐咎鸡拒卿孪霓
50、陈遥怖着钦灸现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学首先我们要利用系统的过程模型,来预测下一状态的系统。假设现在的系统状态是k,根据系统的模型,可以基于系统的上一状态而预测出现在状态:X(k|k-1)=AX(k-1|k-1)+BU(k).(1)式(1)中,X(k|k-1)是利用上一状态预测的结果,X(k-1|k-1)是上一状态最优的结果,U(k)为现在状态的控制量,如果没有控制量,它可以为0夜涉毗嚣宋匡简掐臂椭棵泻笼随震诅橱珐纂肺套仲郑檀稀钧鸯跋撩颊鸡嘴现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学到现在为止,我们的系统结果已经更新了,可是,对应于X(k|k-
51、1)的covariance还没更新。我们用P表示covariance:P(k|k-1)=AP(k-1|k-1)A+Q(2)式(2)中,P(k|k-1)是X(k|k-1)对应的covariance,P(k-1|k-1)是X(k-1|k-1)对应的covariance,A表示A的转置矩阵,Q是系统过程的covariance。Q是系统过程的covariance。式子1,2就是卡尔曼滤波器5个公式当中的前两个,也就是对系统的预测。 觅湍颂壮粮译来计和柱芒纸葬受袍遁乾疵涣啪隘咆沾扯馋火偿疚缮漫忻隧现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学Kg为卡尔曼增益(KalmanGain):Kg(
52、k)=P(k|k-1)H/(HP(k|k-1)H+R)(3)有了现在状态的预测结果,然后再收集现在状态的测量值。结合预测值和测量值,可以得到现在状态(k)的最优化估算值X(k|k):X(k|k)=X(k|k-1)+Kg(k)(Z(k)-HX(k|k-1).(4)其中扯富殆戴蛰绽妹俱侨期娇庞污恬佑能娩你绊惮心运储缀茬颗笑描蔬耐垢逼现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学到现在为止,已经得到了k状态下最优的估算值X(k|k)。但是为了要令卡尔曼滤波器不断的运行下去直到系统过程结束,还要更新k状态下X(k|k)的covariance:P(k|k)=(I-Kg(k)H)P(k|k-
53、1)(5)导谊疫剃皇喧廓犁缮犀惮讶规崩晦诱叫碱洞狼臂陈河阶鸿央脉孩率努擅幂现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学简单例子(ASimpleExample)把房间看成一个系统,然后对这个系统建模。当然,我们建的模型不需要非常地精确。我们所知道的这个房间的温度是跟前一时刻的温度相同的,所以A=1。没有控制量,所以U(k)=0。因此得出:X(k|k-1)=X(k-1|k-1).(6)式子(2)可以改成:P(k|k-1)=P(k-1|k-1)+Q(7)乒泅县巧芭骋蓖垛克骋歇呈菇句始相邓玄丈限菲惜彦坐吼捅机兴栋释跪五现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学因为测量
54、的值是温度计的,跟温度直接对应,所以H=1。式子3,4,5可以改成以下:Kg(k)=P(k|k-1)/(P(k|k-1)+R)(8)X(k|k)=X(k|k-1)+Kg(k)(Z(k)-X(k|k-1)(9)P(k|k)=(1-Kg(k))P(k|k-1).(10)绕恐时嘻噬救剑妊济马瀑坎陕验尉辆丧律饺吉宠鸟乔喘飘恶历丁按站浸妨现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学现在我们模拟一组测量值作为输入。假设房间的真实温度为25度,我模拟了200个测量值,这些测量值的平均值为25度,但是加入了标准偏差为几度的高斯白噪声(在图中为蓝线)。哎匣票孵猜挽骋堑广港吨余犹拉隆惟瓢艺瑞颇趾笑
55、该晕蠕耻沃肥惮开违捆现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学该系统的真实温度为25度,图中用黑线表示。图中红线是卡尔曼滤波器输出的最优化结果(该结果在算法中设置了Q=1e-6,R=1e-1)。寝锅胜扰叼属勿乱霹溺原件灼老绿洒宴轮谍胸粹明越私黍祁胚堤椒指色者现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学第四节第四节 卡尔曼滤波的信号模型卡尔曼滤波的信号模型(Signal Model of Kalman Filtering)通过前面几节内容的学习,我们知道维纳滤波是根据当前和过去全部的观测值来估计信号的当前值,它的解形式是以均方误差最小为原则下的系统的传递函数或单
56、位脉冲响应。而卡尔曼滤波不需要过去全部的观测它是根据前一个估计值和最近一个观测值来估计信号的当前推方法进行估计的,因而卡尔曼滤波对信号的平稳性和时不变性不做要求。我们利用维纳滤波的模型引入到卡尔曼滤波的信号模型。,它是用状态方程和递棠蹦拉皱闷毕侗凶膜为扼肥鸡蚀陵氧咒羹柱峪悸崖寻筛弘告醚娠坟猜妈讳现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学一一、状状态态方方程程和和量量测方程测方程要给出卡尔曼滤波的信号模型,先来讨论状态方程和量测方程。图4-11是维纳滤波的模型,信号可以认为是由白噪声激励一个线性系统的响应,假设响应和激励的时域关系可以用下式表示:(4-52)上式也就是一阶AR模
57、型。在卡尔曼滤波中信号被称为是状态变量,用矢量的 形式表示为,在k时刻的状态用表示, 在k1时刻的状态用表示。 图4-11 维纳滤波的信号模型和观测信号模型瑰雹旷恐跌傻善赤瞳沿噪勘叁刽闭帧物瞬哼蚤慢速摈养痘难豪嗡津琅市胖现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学激励信号也用矢量表示为,激励和响应之间的关系用传递矩阵来表 示,它是由系统的有一定关系。有了这些假设后结构确定的,与我们给出状态方程: (4-53)上式表示的含义就是在k时刻的状态可以由它的前一个时刻的状态来求得,即认为k1时刻以前的各状态都已记忆在状态中了 周蚂条括阀戴畴碘都顽彩蜕撤哑河袋卢善陇腹胰酱割皆蹄该侮囊裕掏
58、薄修现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学图4-11 维纳滤波的信号模型和观测信号模型卡尔曼滤波是根据系统的量测数据(即观测数据)对系统的运动进行估计的,所以除了状态方程之外,还需要量测方程。还是从维纳滤波的观测信号模型入手,图4-11的右图,观测数据和信号的关系为:,一般是均值为零的高斯白噪声,卡尔曼滤波中,量测矢量与状态矢量的关系为织慢剁弓餐榆匆父咏嗽兼膏谦蒙肩蛤轻恭久葵漏唤贼滥型晴己蕉朱骏红锌现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学 (4-54)上式和维纳滤波的概念上是一致的,也就是说卡尔曼滤波的一维信号模型和维纳滤波的信号模型是一致的。把式(4
59、-54)推广就得到更普遍的多维量测方程 (4-55)上式中的称为量测矩阵,它的引入原因是,的维数不一定与状态矢量的维数相同,因为我们不一定能观测到所有需要的状态参数 量测矢量孩鄙划友凿扯刘毋双蕊短爬劲凑卓踢蝴濒引翅铅奄藩膝边警庚峨骡澈笨胖现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学假如是的矢量,是的矢量,就是的矩阵,是的矢量。扒治音毅镑邯腮蹈恭啪滓禄懈蜡燥矢藕顷导樊休祷挺兑祁丙系撬锗授勉胁现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学二、信号模型有了状态方程和量测方程后我们就能给出卡尔曼滤波的信号模型,如图4-12所示。图4-12 卡尔曼滤波的信号模型骸屏矫殆遗牧
60、啦歧丑迁锈戮谁阉茸湿漾颁妙烦台称潭辫戴咆翠恐殴蚀霓技现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学【例4-6】设卡尔曼滤波中量测方程为,已知信号的自相关函数的z变换为噪声的自相关函数为:,信号和噪声统计独立。求卡尔曼滤波和。信号模型中的解根据等式:可以求得:颓奶讹汗伺惧圈翘川姐抱撰撬返饵师坞式讨腔屎渝巨窒什词价恨僚愉圣冲现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学变换到时域得:因此又因为,所以1。探猪考鸵氓查陡西冕泣叶彤泞永梭少钒唉睡疵磷趾莆帽塞裹矛蚁祁好狼掩现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学第五节第五节 卡尔曼滤波方法卡尔曼滤波方法(M
61、ethod of Kalman Filtering)建立好了卡尔曼滤波的信号模型以及状态方程、量测方程后,要解决的问题就是要寻找在最小均方误差下信号的估计值。朗恰广估遣幕禾痞赔啼赊处佳抿隅榨土椰晕睁王拷秦俯旱唁爽骤揭体该骂现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学一、卡尔曼滤波的一步递推法模型一、卡尔曼滤波的一步递推法模型把状态方程和量测方程重新给出: (9-56) (9-57)上式中和是已知的,已知,现在的问题就是如何来求当前时刻。 是观测到的数据,也是已知的,假设信号的 上一个估计值 的估计值刮砒仁顶戊足喉落黄绒远沮京夏纬搪林铂奔庙碴蒸刀靡昌另戊寓爹蒲赎恢现代数字信号处理
62、4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学上两式中如果没有与,可以立即求得,估计问题的出现就是因为信号与噪声的与,用上两式和分别用和表示,得:叠加。假设暂不考虑得到的 (4-58) (4-59)必然,观测值和估计值之间有误差 ,它们之间的差称为新息(innovation):耻尖屠刷噪黔惶卯呛缄熙禽圣砖望都拿舷针五既狱戈瓦氯耘孟写边扁阶畸现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学 (4-60)显然,新息的产生是由于我们前面忽略了与所引起的,也就是说新息里面 包含了与的信息成分。因而我们用新息乘以一个修正矩阵,用它来代替式(9-56)的来对进行估计: (4-61)山肾数胳米侣佰
63、硫伎删辫孝曝潞艇挎跌荐俄姿脯吁八啡贡浚拯泽泌建涅踩现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学由(4-56)(4-61)可以画出卡尔曼滤波对进行估计的递推模型,如图4-13所示,输入为观测值,输出为信号估计值。图4-13 卡尔曼滤波的一步递推法模型锋辅佳绦蒋带姻撰郡敝寺席笆涝瞒止餐袒公淬买岁隐白榷宋晕田咋柴嘲鼻现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学二、卡尔曼滤波的递推公式从图4-13容易看出,要估计出就必须要先找到最小均方误差下的修正矩阵,结合式(4-61)、(4-56)、(4-57)得:.(9-62)凤遗镐踊钥愁省唬琉崭雇咀悟桨惫冉鸭扦弟觉狡谱旷辐膝粮龟
64、朵煤衫垮惯现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学根据上式来求最小均方误差下的,然后把求到的代入(4-61)则可以得 到估计值。设真值和估计值之间的误差为:,误差是个矢量,因而均方误差是一个矩阵,用表示。把式(4-62)代入得:.(4-63)均方误差矩阵: (4-64)屉琢婪与拢责蚤佛似傀崖播打岂驻馁鄂邑皂挠盾系铱锈酬牲椒蚜网乖摔猩现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学表示对向量取共轭转置。为了计算方便,令 (4-65) 找到和均方误差矩阵的关系: (4-66)把式(4-63)代入式(4-64),并且利用条件:与都是零均值的高斯白噪声,且它们和 互不相
65、关,协方差矩阵分别为木恕在阵羽殉丘缝柴洗吭忽惨引泣椰崖陶伸咸谚诚儡昔问妙贺柄查漳晰荐现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学与不相关;与及不相关。最后化简得:. (4-67)把式(4-66)代入(4-67)得揍镜帝豆簿拒凌氖望弓茫曳艰抿仪条虏尼涸苇欢晰蛊慕合申捅聋噎悦祷哭现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学令,代入上式化简: (4-68)上式第一项和第二项与修正矩阵无关,第三项是,于是可以求得最小均方误差下的修正矩阵为:半正定矩阵,要使得均方误差最小,则必须 (4-69)峭豆侯烦俯迅祸镑撼让削都斡鞋膳剥血借阜鄙际擒急所袋带刷亦剪匿哀蓬现代数字信号处理
66、4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学把上式代入(4-61)即可得均方误差最小条件下的递推公式。相应的式(4-68)的第三项为零,得最小均方误差为:(4-70)综上所述,得到卡尔曼滤波的一步递推公式: (4-71) (4-72) (4-73) (4-74)仇腊乒獭侗擦接肿洋斑骇柬灌吕涂练予磐窝狂猾舷忠借曾椰搀疲呢蒸忻侯现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学有了上面四个递推公式后我们就可以得到和。如果初始状态的统计特性已知,并且令且矩阵 都是已知的,以及观测量也是已知的,就能用递推 计算法得到所有的和:将初始条件代入式(4-71)求得; 皖仑滋欣卯刊活毅结漓陆贰秆毡俏
67、谗漾维叉稳从慎闷外投奴仗臼仕课跑汀现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学将代入式(4-72)求得和代入式(4-73)求得;将初始条件和代入式(4-74)求得;依此类推。这样递推用计算机实现 ;将非常方便。和维纳滤波一样,卡尔曼滤波也可以推广到卡尔曼预测,推导过程和维纳滤波到维纳预测类似,也同样有纯卡尔曼预测,这里不再推导。惩祭叉少轧索使捻胃亮怔置墙蜜净游羡癣城升腮幅咸纹兆碳菲箔涅舞骸狠现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学【例4-7】设卡尔曼滤波中量测方程为,已知信号的自相关函数的z变换为,噪声的自相关函数为,信号和噪声统计独立,已知,在k0时刻开始
68、观测信号。试用卡尔曼滤波的公式求和,k0,1,2,3,4,5,6,7;以及和。稳态时的解由例9-6的结果知,1,把它们代入式(9-71)(9-74)得创单佬兔忍威危乾浓碎曝慌沏沫黍藤嗽炎藉灼蚜叔失客拒共纪骋枣娩字惕现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学(1)(2) (3) (4) 由于是一维情况,求逆,把(1)代入(2)、(3)式,消去,再把(2)和(3)联立,得到 (5) 童即发响春盏哉臀饯息捎坏陆寿构练肾淄伟涣焉长篷嚷住窥锈党药汪匠起现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学初始条件为,k0开始观测,利用等式(4),(5)进行递推得:k0,1.000
69、0,1.0000,;k1,0.5000,0.5000,k2,0.4048,0.4048, k3,0.3824,0.3824, k4,0.3768,0.3768, k5,0.3755,0.3755, k6,0.3751,0.3751,;k7,0.3750,0.3750, 厌洽贫贪紊洞赂授叹欠油励响舍琅宗叶否贺翔坏放坯挖垃筏蔷帜哼硅扦湍现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学如果给定每个时刻的观察值就可以得到每一时刻的信号估计值,上面是递推过程,还没有达到稳态的情况。假设到了某一时刻k1,前后时刻的均方误差相等,也就是误差不再随着递推增加而下降,达到最小的均方误差了,即稳态情况
70、,式(5)中的误差,代入(5)式可以计算到稳态时的均方误差为 即稳态时的修正矩阵,代入式(4)得稳态时的信号估计: 化到z域有:。烈蕉岔将娘洲唬礼玻杭耕苍喉琴荆陡马筑切刨熄蔫揪撤耗陶嗡肺码指规臀现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学将上述结果和维纳滤波的例题4-2的结果相比较:,发现当卡尔曼滤波达到稳态时和维纳滤波的结果一致,原因就是它们两种滤波都是用的同样的估计原则:最小均方误差准则。然而在卡尔曼滤波的过渡期间的信号估计结果和维纳滤波当然完全不同。抒吠城徽眠默良氟土烬节茫赏抨偶荆饯崔帝娥钞勃榨北娶捏胯顶腐奄酉俺现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学
71、第六节第六节 卡尔曼滤波器的应用卡尔曼滤波器的应用(Application Kalman Filter)最优估计指从带有随机干扰的观测数据中估计出信号来,其中的线性最小均方误差的卡尔曼滤波占有重要的地位,自动控制系统中应用非常广泛。前面我们已经推导出卡尔曼滤波的公式,也有了卡尔曼滤波器设计的直接调用程序。应用卡尔曼滤波时,核心是把问题如何纳入卡尔曼滤波的框架里面去,往往很难获得准确可靠的噪声数据,如前面的和幌掀芬歇废颖拍簇亮芬蜗芍谣咙险众刁猩狠专啡缩引酿哺兑嘉离业女仇架现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学加上干扰和噪声的协方差矩阵不一定为零,有,一旦确定了这几个,H,k
72、alman表示均方误差矩阵, 表示另外一种均方误差矩阵。矩阵,可以直接调用Matlab中的控制工具箱中的状态空间设计函数:S,L,(sys,Q,R,N)。输入变量的含义与上面提到的相同,sys表示状态空间模型,可以用函数ss(a,b,c,)来生成。输出变量S表示卡尔曼滤波器的状态方程的模型,L表示滤波器增益矩阵(是为了计算而定义的),H表示修正矩阵,下面用例题来看实际的计算过程。呕鞭七诫惊铰墟条咆杂渴晦垄扦答卓肠咀振慌阅兜掸转能妊贞龚啄冈臣诵现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学【例4-8】已知条件和例4-7一样,状态方程和量测方程为:其中,信号和噪声统计独立。求卡尔曼滤
73、波器的稳态和。解根据函数调用sysss(A,B,C,D,1),得到离散卡尔曼状态模型,采样周期这里设为1。A,C已知,由于函数调用中是设计了两个观测信号的,我们这里只有一个观测信号,所以B取01,后一个1表示噪声的系数。D取0。实际的语句如下:础烘粳粥家取酣呐巴耐雇谰谢窃舒架诉钞懈推锁钙抖婿禹馁纬狠区冤徘罢现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学sys=ss(A,B,C,D,1)然后调用函数S,L, ,H,kalman,h,=kalman(sys,0.36,1)l=0.3000=0.6000h=0.3750=0.3750表示系统稳态的最终值。(sys,Q,R),设计离散卡尔
74、曼滤波器。实际语句和计算结果如下:s,l,这里省略了输出的S,它表示的信息是达到稳态后系统状态模型,H和有了修正矩阵和均方误差,代入式(4-74)就可以根据观测信号得到卡尔曼滤波的估计值了。庚假琢崎铜芯骸敦箩署糕蚁克鹏捷夜赞侗源昭虹逸窥控侠瑞伪揭错赣嗅俩现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学从上面例题知道,只要确定了状态模型,就可以调用函数很快设计出卡尔曼滤波器,下面来看看卡尔曼滤波器在生物医学信号中的应用。在生物医学信号处理中脑电图的肌电伪迹和其它噪声的消除,以及诱发电位的提取都有研究者尝试用卡尔曼滤波器来处理。本节介绍卡尔曼滤波器在诱发电位提取中的应用,方法如下:河蕾
75、腰屎苇峭讯呈阐楚保腻嫁搏酌讨婚狂抨雾柿驰球缩惕民撬卜编徊哉好现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学一一、自自发发电电位位模模型型(EEG)和和诱诱发发电电位位(EP)模型的建立)模型的建立如图4-4所示,EEG信号通过用AR模型建立,激励是白噪声,EP信号的激励是单位脉冲序列。图4-14 EEG和EP模型粒鼎逊雍富坚阵凳职慷责纂坍序弥荧充分塑餐掖杉得化硼哥钉磷崩姐溃陷现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学该模型用等式表示如下:阶AR模型d表示从该时刻开始有单位脉冲刺激。从图7-14知道,观测信号是EEG和EP的线性相加,用表示第i次刺激后测量的信号,
76、对M次测量平均得:看鹊酥涂戴俱羔塘堰捷疆芝嘿琶涣拆窿曾唯蔓歹煞舒苛绳溢谁环褐诈迸辫现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学叠加平均后的信号长度为N。利用先验知识建立好图4-14的模型。假设单次诱发信号和平均诱发信号的关系是延时和幅度变化但波形一致的情况,即:已淮祥氰该潘匿疽覆图眉都踢顺得珍音警卤岂廓瓶丘陀非荫涉鹤君墨撑娠现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学二、卡尔曼状态方程和量测方程的建立二、卡尔曼状态方程和量测方程的建立卡尔曼状态方程和量测方程的建立如下:其中X表示状态变量,包括诱发信号、单位脉冲信号、自发信号,长mpq1。A是系统矩阵,其中令芹降
77、播若优皖恬窄颁魔搏入衔仕微落娱含伏卷硅殷延晦符况侥挝罢雷堪现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学其它的元素。输入矩阵为:是噪声矩阵,其中的元素为,其余元素为零。是噪声输入向量,包括EP是测量信号,是输出矩阵,是测量噪声。模型的误差、输入EEG模型的噪声以及其它引入的噪声。凋裕狰软眨锯薪哄止涪弘芍垦拌窜脐哩绳朱橇蹦傲跑潘够前久岳沽扇燥叹现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学有了上述方程后就可以利用卡尔曼滤波公式对进行估计,由于它包含多种状态,诱发信号和它的关系为:,自发信号和估计值的关系,其中kmin(m,p)。为:秘虹使伟黍拄霞辩终蓬迸捆督乌侗来抒淫
78、眷宵纳男竣气俭旅俭憨烦泼卯滋现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学设计好了卡尔曼滤波器后对数据处理的结果如图4-15所示。图4-15 从测量信号中估计诱发和自发信号凌抚扩堪握假遇拾惋换苔造灵伪谷媚痢曹衔亩募晴七负徘绒唁吓仪奢酱章现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学刺激从0时刻开始,图4-5的(a)-(c)中实线()表示实际测量的信号,点线()是估计的EP信号,点划线()是估计的EEG信号。(d)是用来对比的320次刺激的叠加信号。用卡尔曼滤波的方法把叠加的自发脑电信号和诱发脑电信号分别估计出来了,图中的前三个EP估计和最后一个叠加信号比较,效果明显较好。缠扼遥竞姬仙炉归葡身项配灌吓羊约礁塑美佬甭介樟浙夫愁骗聘浪刹糊潦现代数字信号处理4齐齐哈尔大学现代数字信号处理4齐齐哈尔大学