1、1第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波第1页，共137页。2第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波 信道信道s(n)x(n)()ns(n)：原始输入（发射）信号，随机平稳：原始输入（发射）信号，随机平稳信道噪声（测量噪声）信道噪声（测量噪声）x(n)：接收（测量）信号，：接收（测量）信号，随机平稳随机平稳问题提出：问题提出：第2页，共137页。3第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波准则：最大后验准则。均方准则，最大似然准则，准则：最大后验准则。均方准则，最大似然准则，滤波器滤波器 h(n)()()y ns nx(n)s(n)FIR,IIR逼近
2、（准则）逼近（准则）已知：已知：s(n),(n)的统计特性，要求：的统计特性，要求：设计线性移不变滤波器设计线性移不变滤波器 h(n),从从x(n)中恢复中恢复s(n)线性均方准则（最小二乘滤波）线性均方准则（最小二乘滤波）第3页，共137页。4第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波()()()()T11nu nu nu nM轾=-+臌uL()nu()u n2.2.横向滤波器结构横向滤波器结构第4页，共137页。5第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波iw*wT011Mwww-=wL()d n()d n()e n()()()e nd nd n=-()()()()
3、1HT0Miid nw u ninn-*=-=w uuw第5页，共137页。6第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波假设 由信号 与噪声 组成 如果 ，上图的系统称为滤波（filtering）；如果 ，上图的系统称为预测（prediction）；如果 ，上图的系统称为平滑 （smoothing）。()()()unsnvn=+()u n()s n()v n()()dns n=()()00,0d ns nnn=+()()00,0d ns nnn=+第6页，共137页。7第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波4.2 4.2 维纳滤波器的离散形式维纳滤波器的离散形式时域
4、解时域解 22()()()()()()()()()()()(),()()()ix ns nv ny ns nh i x nix nh ne ns ns nE e nE s ns n误差，平稳随机）均方误差：(n)最小均方误差准则：最小均方误差准则：2op寻求h(n)使(n)=Ee(n)取最小值（加性干扰）（加性干扰）d(n)=s(n)第7页，共137页。8第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波为此令为此令2()()2()()0()()()()()()()(),()()()()0,(inE e nEe ne nh jh jh je ns nh j x nix njjh jh jE
5、 e n x njj 正交标量方程组）一、维纳一、维纳霍夫方程霍夫方程 正交性原理：正交性原理：使代价函数最小化的充要条件是使代价函数最小化的充要条件是 n 时刻的时刻的最优估计误差正交于最优估计误差正交于n 时刻滤波器的每个输入值，或者时刻滤波器的每个输入值，或者说正交于说正交于n时刻的输入信号空间。时刻的输入信号空间。第8页，共137页。9第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波推论：推论：n时刻的最优估计误差正交于时刻的最优估计误差正交于n时刻滤时刻滤波器的最优输出值波器的最优输出值0)()(00nenyE000)()(kkknxhny一、维纳一、维纳霍夫方程霍夫方程第9页
6、，共137页。10第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波由正交方程可得：()()()()()()()()()()()0iiE e n x njEs nh i x nix njE s n x njE x ni x nj h i()()()()()iE s n x njh i E x ni x nj即：一、维纳一、维纳霍夫方程霍夫方程第10页，共137页。11第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波定义定义()()()(),()()()()()sxxxRmE s n x nms n x nRmE x n x nmx n的互相关的自相关可得可得(),sxxxiRmh i
7、Rmim（）（）维纳维纳霍夫（霍夫（Wiener-Holf)方程或标准方程方程或标准方程求和范围（求和范围（i）随滤波器的不同取不同区间）随滤波器的不同取不同区间:0 1:0 FIRMIIR（非因果）；（因果）一、维纳一、维纳霍夫方程霍夫方程第11页，共137页。12第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波FIR 维纳滤波器维纳滤波器令令(0)(1)(1)()()(1)(1)Hhhh NX nx nx nx nN10()()()()()()()()MkTy ns nx nh nh k x nkH X nXn H()()()()()()()()()Te hs ns ns ny ns
8、 nXn Hs nH X n对对FIR结构，假设其长度为结构，假设其长度为N,期望信号为期望信号为s(n)第12页，共137页。13第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波22()()()2 ()()()()ThE e nE snE s n Xn HH E X n Xn H22()()()()()()2()()()()Te ns nXn Hs nH X nsns n Xn HH X n Xn H第13页，共137页。14第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波()()()(1)()(1)(1)()(1)(1)(1)(1)(1)()(1)(1)(1)(1)ExnxnE
9、xnxnExnxn NExnxnExnxnExnxn NExn NxnExn NxnExn Nxn N ToeplitzToeplitz矩阵矩阵,N,NXN N对称半正定对称半正定()()RE X n Xn(0)(1)(1)(1)(0)(2)(1)(2)(0)xxxxxxxxxxxxxxxxxxrrrNrrrNrNrNr令令第14页，共137页。15第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波令令22()sE sn()()(0)()()(1)(1)()(1)(1)()sxsxsxPE X n s nrE x n s nrE x ns nrNE x nNs n2()2(shP HH R
10、H 二次型问题）代入可得：代入可得：第15页，共137页。16第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波则则222d P Hd H RHPRHHdHdH （）（）令令0H(0)(1)(1)(0)(0)(1)(0)(2)(1)(1)(1)(2)(0)(1)(1)xxxxxxsxxxxxxxsxxxxxxxsxrrr Nrhrrr Nrhr Nr Nrr Nh NRHP（维纳-霍夫方程）1OPHR P其解为：，且第16页，共137页。17第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波22122()222,()OPOPOPOPOPOPOPOPOPOPOPTTsTTsTTsTTTs
11、HP HH RHP HH R R PP HH PH PP HH P可以看出，均方误差与滤波器的单位脉冲响应是一个二次函数可以看出，均方误差与滤波器的单位脉冲响应是一个二次函数关系。由于单位脉冲响应关系。由于单位脉冲响应h(n)为为M维向量，因此均方误差维向量，因此均方误差是一个超椭圆抛物形曲面，该曲面有极小点存在。当滤是一个超椭圆抛物形曲面，该曲面有极小点存在。当滤波器工作于最佳状态时，均方误差取得最小值。波器工作于最佳状态时，均方误差取得最小值。第17页，共137页。18第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波 上式表明已知期望信号与观测数据的互相关函数及观上式表明已知期望信号
12、与观测数据的互相关函数及观测数据的自相关函数时，可以通过矩阵求逆运算，得到测数据的自相关函数时，可以通过矩阵求逆运算，得到维纳滤波器的最佳解。同时可以看到，直接从时域求解维纳滤波器的最佳解。同时可以看到，直接从时域求解因果的维纳滤波器，当选择的滤波器的长度因果的维纳滤波器，当选择的滤波器的长度M较大时，较大时，计算工作量很大，并且需要计算计算工作量很大，并且需要计算Rxx 的逆矩阵，从而要的逆矩阵，从而要求的存贮量也很大。此外，在具体实现时，滤波器的长求的存贮量也很大。此外，在具体实现时，滤波器的长度是由实验来确定的，如果想通过增加长度提高逼近的度是由实验来确定的，如果想通过增加长度提高逼近的
13、精度，就需要在新精度，就需要在新N基础上重新进行计算。因此，从时域基础上重新进行计算。因此，从时域求解维纳滤波器，并不是一个有效的方法。求解维纳滤波器，并不是一个有效的方法。12()OPOPOPTsHRPHH P即即第18页，共137页。19第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波图图 期望信号、期望信号、估计值与误差信号的几何关系估计值与误差信号的几何关系()()()ooe nd nd n=-由正交方程可知：误差与输入信号矢量正交，可推得其与估计值也正交，用下图表示。几何解释：几何解释：第19页，共137页。20第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波图表明在滤波
14、器处于最佳工作状态时，估计值加上估计偏差等于期望信号，即 d nynenoptopt()()()注意我们所研究的是随机信号，图中各矢量的几何表示应理解为相应量的统计平均或者是数学期望。再从能量的角度来看，假定输入信号和期望信号都是零均值,应用正交性原理，则,因此在滤波器处于最佳状态时，估计值的能量总是小于等于期望信号的能量。dyE eopt222opt|第20页，共137页。21第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波非因果性考虑非因果性考虑 可以证明：非因果可以证明：非因果WienerWiener滤波器的性能滤波器的性能(误差方差性误差方差性能能)要优于因果要优于因果Wiene
15、rWiener滤波器滤波器(参见郑南宁编参见郑南宁编数字信号数字信号处理处理)。所以，在实际。所以，在实际FIRFIR滤波器中，常用时延方法用滤波器中，常用时延方法用可实現的因果系统逼近非因果系统。可实現的因果系统逼近非因果系统。MMkkknxwny)()(0令令Mki第21页，共137页。22第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波MiiMiiiMnxwMinxwny200200)()()(MiiinxwMny200)()(n 因果系统因果系统n时刻的输出可以逼近非因果系统时刻的输出可以逼近非因果系统(n-M)的的输出。输出。第22页，共137页。23第四章第四章 维纳滤波和卡
16、尔曼滤维纳滤波和卡尔曼滤波波2),(8()(),()x nv nv nns(n)=sin(周期为 的确定信号）4nsin(（为方差的白噪声）4例例1 1：设计设计N=4N=4的的FIRFIR最佳滤波器最佳滤波器已知：已知：第23页，共137页。24第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波解：解：2222111022 22 2(0)(1)(2)(3)1110(1)(0)(1)(2)22 22 2(2)(1)(0)(1)11102(3)(2)(1)(0)2 22 2111022 22 2xxxxxxxxxxxxxxxxrrrrrrrrrrrrrrrr()RE XnT（n)X82222
17、181121(0)()sin()8821221(1)()(1sin()sin(1)08882 2xnxnrE xnnE vrE x n x nnn注：，第24页，共137页。25第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波12 2PE X121（n)s(n)=，2 20222min222min2min2min11()22(1)2(1)0,0,(1126114TXopE e nH P确定），(0分贝信噪比），同样：同样：121112 2opHRP1212 20第25页，共137页。26第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波 例例 设设y(n)=x(n)+v2(n)，v2
18、(n)是一白噪声，方差是一白噪声，方差22=0.1。期望信号。期望信号x1(n)的信号模型如图的信号模型如图(a)所示，其中白所示，其中白噪声噪声v1(n)的方差的方差21=0.27，且，且b0=0.8458。x(n)的信号模型的信号模型如图（如图（b）所示，）所示，b1=0.9458。假定。假定v1(n)与与v2(n)、x1(n)与与y(n)不相关，并都是实信号。设计一个维纳滤波器，不相关，并都是实信号。设计一个维纳滤波器，得到该信号的最佳估计，要求滤波器是一长度为得到该信号的最佳估计，要求滤波器是一长度为 2 的的FIR滤波器。滤波器。第26页，共137页。27第四章第四章 维纳滤波和卡尔
19、曼滤维纳滤波和卡尔曼滤波波图图 输入信号与观测数据的模型输入信号与观测数据的模型 解解 这个问题属于直接应用维纳这个问题属于直接应用维纳-霍夫方程的典型问题，其霍夫方程的典型问题，其关键在于求出观测信号的自相关函数和观测信号与期望信号关键在于求出观测信号的自相关函数和观测信号与期望信号的互相关函数。的互相关函数。图图 维纳滤波器的框图维纳滤波器的框图 第27页，共137页。28第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波根据题意，画出这个维纳滤波器的框图，如图所示。根据题意，画出这个维纳滤波器的框图，如图所示。用用H1(z)和和H2(z)分别表示分别表示x1(n)和和x(n)的信号模
20、型，那么滤波器的输的信号模型，那么滤波器的输入信号入信号x(n)可以看作是可以看作是v1(n)通过通过H1(z)和和H(z)级联后的输出，级联后的输出，H1(z)和和H(z)级联后的等效系统用级联后的等效系统用H(z)表示，输出信号表示，输出信号y(n)就就等于等于x(n)和和v2(n)之和。因此求出输出信号的自相关函数矩阵之和。因此求出输出信号的自相关函数矩阵Ryy和输出信号与期望信号的互相关矩阵和输出信号与期望信号的互相关矩阵Ryd是解决问题的关是解决问题的关键。相关函数矩阵由相关函数值组成，已知键。相关函数矩阵由相关函数值组成，已知x(n)与与v2(n)不相关，不相关，那么那么 2 2(
21、)()()yyxxv vrmrmrm第28页，共137页。29第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波(1)求出期望信号的方差。根据图求出期望信号的方差。根据图(a)，期望信号的时间序列，期望信号的时间序列模型所对应的差分方程为模型所对应的差分方程为 x1(n)=v1(n)-b0 x1(n-1)这里，这里，b0=0.8458,由于由于x1(n)的均值为零，其方差与自相关函数的均值为零，其方差与自相关函数在零点的值相等。在零点的值相等。11122122210 111122210(0)()()2()(1)(1)xxxRE x nE v nb v n x nb x nb 1222122
22、00.270.948611(0.8458)dxb 第29页，共137页。30第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波(2)计算输入信号和输出信号的自相关函数矩阵。根据自计算输入信号和输出信号的自相关函数矩阵。根据自相关函数、功率谱密度和时间序列信号模型的等价关系，已知相关函数、功率谱密度和时间序列信号模型的等价关系，已知时间序列信号模型，就可以求出自相关函数。这里，信号的模时间序列信号模型，就可以求出自相关函数。这里，信号的模型型H(z)可以通过计算得到。可以通过计算得到。12111()()()(10.8458)(10.9458)H zH z H zzz 这是一个二阶系统，所对应
23、的差分方程为这是一个二阶系统，所对应的差分方程为 x(n)+a1x(n-1)+a2x(n-2)=v1(n)式中，式中，a1=-0.1，a2=-0.8。由于。由于v1(n)、v2(n)的均值为零，因此，的均值为零，因此，x(n)的均值为的均值为0。将方程两边同乘以。将方程两边同乘以x*(n-m)，并取数学期望，得，并取数学期望，得 rxx(m)+a1rxx(m-1)+a2rxx(m-2)=0 m0(1)rxx(0)+a1rxx(1)+a2rxx(2)=21 m=0(2)第30页，共137页。31第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波对方程（对方程（1）取）取m=1,2，得到，得到
24、 rxx(1)+a1rxx(0)+a2rxx(1)=0(3)rxx(2)+a1rxx(1)+a2rxx(0)=0(4)方程方程(2)、(3)、(4)联立求解，得联立求解，得 222122222211211 0.80.27(0)11(1)1 0.8(1 0.8)(0.1)0.1(1)0.511 0.8xxxxxaraaaara 至此，至此，输入信号的自相关矩阵输入信号的自相关矩阵Rxx可以写出：可以写出：(0)(1)1 0.5(1)(0)0.5 1xxxxxxxxxxrrRrr 第31页，共137页。32第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波v2(n)是一个零均值的白噪声，它的自
25、相关函数矩阵是一个零均值的白噪声，它的自相关函数矩阵呈对角形，呈对角形，且且 ，2 222(0)v vr 2 22 22 22 22 2(0)(1)0.1 0(1)(0)00.1v vv vv vv vv vrrRrr 因此，输出信号的自相关因此，输出信号的自相关Ryy为为 2 22 22 22 2(0)(1)(1)(0)(1)1.1 0.5(1)(0)(1)(0)(0)0.5 1.1xxv vxxv vyyyyyyyyyyxxv vxxv vrrrrrrRrrrrrr 第32页，共137页。33第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波(3)计算输出信号与期望信号的互相关函数矩阵
26、。计算输出信号与期望信号的互相关函数矩阵。由由于两个信号都是实信号，故于两个信号都是实信号，故 ryd(m)=Ey(n)d(n-m)=Ey(n)x1(n-m)=E(x(n)+v2(n)x1(n-m)=Ex(n)x1(n-m)m=0,1 根据系统根据系统H2(z)的输入与输出的关系，有的输入与输出的关系，有 x1(n)-b1x(n-1)=x(n)推出推出 x1(n)=x(n)+b1x(n-1)这样这样 ryd(m)=Ex(n)x1(n-m)=Ex(n)(x(n-m)+b1x(n-1-m)=rxx(m)+b1rxx(m-1)第33页，共137页。34第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼
27、滤波波将将m=0,m=1代入上式代入上式,得得 ryd(0)=rxx(0)+b1rxx(-1)=1-0.94580.5=0.5272ryd(1)=rxx(1)+b1rxx(0)=0.5-0.94581=-0.4458 因此，输出信号与期望信号的互相关因此，输出信号与期望信号的互相关Ryd为为(0)0.5272(1)0.4458ydydydrRr 求出输出信号自相关的逆矩阵求出输出信号自相关的逆矩阵,并乘以并乘以Ryd,就可以得到就可以得到维纳滤波器的最佳解维纳滤波器的最佳解Wopt:1122(0)(1)(0)(1)1(0)(1)(1)(0)(1)(0)1.14560.52080.5208 1.
28、1456yyrrrrRrrrrrr 第34页，共137页。35第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波11.14560.5208 0.52720.83600.4580.78530.54281.1456optyyydWR R 可以计算出该维纳滤波达到最佳状态最小值可以计算出该维纳滤波达到最佳状态最小值E|e(n)|2min：22*Tminopt|()|()0.83600.94860.52720.44580.75830.1579dydE e nRW 第35页，共137页。36第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波2.3 离散维纳滤波器的离散维纳滤波器的 z 域解域解标
29、准方程：标准方程：双边双边Z Z变换变换一、非因果一、非因果IIR 维纳滤波器维纳滤波器第36页，共137页。37第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波假设信号和噪声不相关，即假设信号和噪声不相关，即rsv(m)=0，则，则 Sxs(z)=Sss(z)Sxx(z)=Sss(z)+Svv(z)则有则有 opt()()()()()()xsssxxssvvSzSzHzSzSzSz 当噪声为当噪声为0时，信号全部通过；当信号为时，信号全部通过；当信号为0时，时，噪声全部噪声全部被抑制掉，因此维纳滤波确有滤除噪声的能力。把信号的被抑制掉，因此维纳滤波确有滤除噪声的能力。把信号的频谱用频谱
30、用Pss(ej)表示，噪声的频谱用表示，噪声的频谱用Pvv(ej)表示，那么非表示，那么非因果的维纳滤波器的传输函数因果的维纳滤波器的传输函数Hopt(ej)的幅频特性的幅频特性如图所如图所示示。第37页，共137页。38第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波jopt1(e)10H Pss(ej)0,Pvv(ej)=0 Pss(ej)0,Pvv(ej)0 Pss(ej)=0,Pvv(ej)0 图图 非因果维纳滤波器的传输函数的幅频特性非因果维纳滤波器的传输函数的幅频特性 第38页，共137页。39第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波 然而实际的系统都是因果的。
31、对于一个因然而实际的系统都是因果的。对于一个因果系统，不能直接转入频域求解的原因是由于果系统，不能直接转入频域求解的原因是由于输入信号与期望信号的互相关序列是一个因果输入信号与期望信号的互相关序列是一个因果序列，如果能够把因果维纳滤波器的求解问题序列，如果能够把因果维纳滤波器的求解问题转化为非因果问题，求解方法将大大简化。那转化为非因果问题，求解方法将大大简化。那么怎样把一个因果序列转化为一个非因果序列么怎样把一个因果序列转化为一个非因果序列呢？呢？因果情况处理思路：因果情况处理思路：第39页，共137页。40第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波标准方程：标准方程：v(n)x
32、(n)S(n)信道信道s(n)h(n)s(n)IIR一般情况下直接求解比较困难。但如果滤波器输入一般情况下直接求解比较困难。但如果滤波器输入是白噪声，则维纳是白噪声，则维纳-霍夫方程容易求解；而任何平稳霍夫方程容易求解；而任何平稳随机信号可变换为等效的白噪声过程，故借助谱分解随机信号可变换为等效的白噪声过程，故借助谱分解定理可找到一种简单解决方法。定理可找到一种简单解决方法。二、因果二、因果IIR IIR 维纳滤波器维纳滤波器第40页，共137页。41第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波 回顾前面讲到的时间序列信号模型，假设回顾前面讲到的时间序列信号模型，假设x(n)的信的信
33、号模型号模型B(z)已知（如图已知（如图(a)所示），求出信号模型的逆系所示），求出信号模型的逆系统统B-1(z)，并将并将x(n)作为输入，那么逆系统作为输入，那么逆系统B-1(z)的输出的输出(n)为白噪声，白化滤波器（如图为白噪声，白化滤波器（如图(b)所示）。所示）。图图 x(n)的时间序列信号模型及其白化滤波器的时间序列信号模型及其白化滤波器(a)(b)第41页，共137页。42第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波（1 1）若滤波器的输入是白噪声时）若滤波器的输入是白噪声时20()():x nw n2()()()xxwwwrmrmm 2200()()()()()()
34、,0sxwwwwiiR mhi Rm ihim ih mm 21()(),0swwh nRnn对应传递函数：对应传递函数：()()swswSzSz表示取的因果部分（极点在单位园内）21()()swwH zSz第42页，共137页。43第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波（2 2）若滤波器的输入是平稳随机信号时）若滤波器的输入是平稳随机信号时则则x(n)x(n)可看作是由白噪声可看作是由白噪声w(n)w(n)激励一个线性移不变系统激励一个线性移不变系统的输出的输出)zB zzN（D（1()B z()Gz白化白化信号模型信号模型滤波器滤波器w(n)w(n)y(n)y(n)w(n)
35、w(n)x(n)x(n)b(n)b(n)20,s(n)s(n)1()()()CHzGzB zB(z)B(z)为有理分式，为有理分式，N(z),D(z)N(z),D(z)为最小为最小相位多项式相位多项式问题转化为求问题转化为求 的问题的问题(),()B zGz第43页，共137页。44第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波111()()()()()()()()()()()()()()()()swswsxsxRns nwnSzS z W zRns nxns nwnbnSzS z W zB z 1()()()sxswSzSzB z即由（由（1 1）221()11()()()sxsww
36、wSzGzSzB z21()1()()()sxCwSzHzB zB z11()()2ncCh nHz zdzj 其中其中211(),()()()()()wxxsxsxB zSzSzSzB zB z可由的谱分解得到可由进行因果、逆因果分解得到第44页，共137页。45第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波 具体思路如图所示。用白噪声作为待求的维纳滤波器具体思路如图所示。用白噪声作为待求的维纳滤波器的输入，设定的输入，设定1/B(z)为信号为信号x(n)的白化滤波器的传输函数，的白化滤波器的传输函数，那么维纳滤波器的传输函数那么维纳滤波器的传输函数G(z)的关系为的关系为()()(
37、)G zH zB z 因此，维纳滤波器的传输函数因此，维纳滤波器的传输函数H(z)的求解转化为的求解转化为G(z)的求解。的求解。图图 维纳滤波解题思路维纳滤波解题思路 2、第45页，共137页。46第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波假设待求维纳滤波器的单位脉冲响应为假设待求维纳滤波器的单位脉冲响应为g(n)，期望信号，期望信号d(n)=s(n)，系统的输出信号，系统的输出信号y(n)=s(n)，g(n)是是G(z)的逆的逆Z变换，变换，则输出信号则输出信号可表示为可表示为()()()()()()ky ns nng ng knk 第46页，共137页。47第四章第四章 维纳
38、滤波和卡尔曼滤维纳滤波和卡尔曼滤波波222*22*ss|()|()()()|()|()()()()()()()()()()(0)|()|()()(kkkkkskE e nE s ng knkE s nEg k g rnknrEg knk s ng knk s nrg kg k rkg k *22ss2)()()|(0)()skksskkrkrkrrg k 第47页，共137页。48第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波 可以看出，均方误差的第一项和第三项都是非负数可以看出，均方误差的第一项和第三项都是非负数,要使均方误差为最小，当且仅当要使均方误差为最小，当且仅当 ()()0s
39、rkg k -k 因此因此g(n)的最佳值为的最佳值为 opt2()()srkgk -k 对上式两边同时做对上式两边同时做Z变换，得到变换，得到 opt2()()sSzGz 第48页，共137页。49第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波这样，非因果维纳滤波器的最佳解为这样，非因果维纳滤波器的最佳解为 optopt2()()1()()()sGzSzHzB zB z 因为因为s(n)=s(n)*(n)，且，且x(n)=(n)*b(n)，根据相关卷，根据相关卷积定理，得到积定理，得到 rxs(m)=rs(m)*b(-m)对上式两边做对上式两边做Z变换，得到变换，得到 Sxs(z)=
40、Ss(z)B(z-1)因此因此 1()()()xssSzSzB z 第49页，共137页。50第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波根据根据x(n)的信号模型，得到非因果的维纳滤波器的复的信号模型，得到非因果的维纳滤波器的复频域最佳解的一般表达式频域最佳解的一般表达式 opt21()()11()()()()xsxsxxSzSzHzB z B zSz 假定信号与噪声不相关，即当假定信号与噪声不相关，即当Es(n)v(n)=0时，有时，有 rxs(m)=E(s(n)+v(n)*s(n+m)=rss(m)rxx(m)=E(s(n)+v(n)*(s(n+m)+v(n+m)=rss(m)
41、+rvv(m)对上边两式做对上边两式做Z变换，得到变换，得到 Sxs(z)=Sss(z)Sxx(z)=Sss(z)+Svv(z)第50页，共137页。51第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波1()()()sssSzSzB z 信号和噪声不相关时，非因果维纳滤波器的复频域信号和噪声不相关时，非因果维纳滤波器的复频域最佳解和频率响应分别为最佳解和频率响应分别为 xsssxxssvvSzSzHzSzSzSzopt()()()()()()jjoptjj(e)()(e)(e)(e)()()ssssssvvssvvSPHSSPP 第51页，共137页。52第四章第四章 维纳滤波和卡尔曼滤
42、维纳滤波和卡尔曼滤波波滤波器的最小均方误差滤波器的最小均方误差E|e(n)|2min的计算，的计算，22min2|()|()|(0)ssskrkE e nr 根据围线积分法求逆根据围线积分法求逆Z变换的公式变换的公式，rss(m)用下式表示用下式表示:11()()d2jmssssCrmSz zz 得出得出 11(0)()d2jssssCrSz zz 第52页，共137页。53第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波由复卷积定理由复卷积定理*11d()()()2jCnzx n y nX z Yzz 取取y(n)=x(n)，有有 211d|()|()2jCnzx nX z X zz
43、因此因此 211d|()|()()2jsssCnzrkSz Szz 21min211d|()|()()()2jsssCzE e nSzSz Szz 第53页，共137页。54第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波12min211()()11d|()|()2j()()1d()()()2jsssssCssoptssCSz SzzE e nSzB zB zzzSzHz Szz 因为实信号的自相关函数是偶函数，即因为实信号的自相关函数是偶函数，即rss(m)=rss(-m)，因此因此 Sss(z)=Sss(z-1)假定信号与噪声不相关，假定信号与噪声不相关，Es(n)v(n)=0，则则
44、 21min()1d|()|()()2j()()()1d2j()ssssssCxxssvvCxxSzzE e nSzSzSzzSz SzzSzz 第54页，共137页。55第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波 若维纳滤波器是一个因果滤波器，若维纳滤波器是一个因果滤波器，要求要求 g(n)=0 n0 则滤波器的输出信号则滤波器的输出信号 0()()()()()()ky ns nng ng knk 估计误差的均方值估计误差的均方值 E|e(n)|2=E|s(n)-y(n)|2 得到得到 222200()1|()|(0)()|()|sssskkrkE e nrg krk 因果维纳滤
45、波器的求解因果维纳滤波器的求解第55页，共137页。56第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波要使均方误差取得最小值，要使均方误差取得最小值，当且仅当当且仅当 2opt2()0()00()()ssrnngnnrnu n 令令 0optopt2()()()()1()ZT()()nnsssnnsSzrn u n zrn zGzgnSz 第56页，共137页。57第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波opt21()1()()xsSzGzB z 所以因果维纳滤波器的复频域最佳解为所以因果维纳滤波器的复频域最佳解为 optopt21()()11()()()()xsGzSz
46、HzB zB zB z 第57页，共137页。58第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波维纳滤波的最小均方误差为维纳滤波的最小均方误差为 22min20*212121|()|()|(0)1(0)()()()11d(0)()()2j()()11d()2j()()1()2jsssksssskssssCxsxsssCssCrkE e nrrrk u k rkzrSzSzzSzSzzSzB zB zzSz 1optd()()xszHz Szz 第58页，共137页。59第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波 比较可以看出因果维纳滤波器的最小均方误差与非因果维比较可以看出
47、因果维纳滤波器的最小均方误差与非因果维纳滤波器的最小均方误差的形式相同，但公式中的纳滤波器的最小均方误差的形式相同，但公式中的Hopt(z)的的表达式不同。表达式不同。前面已经导出，前面已经导出，对于非因果情况，对于非因果情况，22min2|()|()|(0)ssskrkE e nr 对于因果情况，对于因果情况，22min20|()|()|(0)ssskrkE e nr 比较两式，它们的第二项求和域不同，因为因果情况比较两式，它们的第二项求和域不同，因为因果情况下下,k=0+，因此可以说明非因果情况的因此可以说明非因果情况的E|e(n)|2min一定一定小于等于因果情况小于等于因果情况E|e(
48、n)|2min。在具体计算时。在具体计算时,可以选择可以选择单位圆作为积分曲线，单位圆作为积分曲线，应用留数定理，应用留数定理，计算积分函数在单计算积分函数在单位圆内的极点的留数来得到。位圆内的极点的留数来得到。第59页，共137页。60第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波因果维纳滤波器设计的一般步骤：因果维纳滤波器设计的一般步骤：(1)根据观测信号根据观测信号x(n)的功率谱求出它所对应的信号模型的功率谱求出它所对应的信号模型的传输函数，即采用谱分解的方法得到的传输函数，即采用谱分解的方法得到B(z)。具体方法为。具体方法为Sxx(z)=2B(z)B(z-1)，把单位圆内的
49、零极点分配给，把单位圆内的零极点分配给B(z)，单位，单位圆外的零极点分配给圆外的零极点分配给B(z-1)，系数分配给，系数分配给2。(2)求求 的的Z反变换，取其因果部分再做反变换，取其因果部分再做Z变换，即变换，即舍掉单位圆外的极点，得舍掉单位圆外的极点，得 (3)计算计算Hopt(z)，E|e(n)|2min。xsSzB z1()()xsSzB z1()()第60页，共137页。61第四章第四章 维纳滤波和卡尔曼滤维纳滤波和卡尔曼滤波波例例 已知已知 10.36()(10.8)(10.8)ssSzzz 信号和噪声不相关，即信号和噪声不相关，即rsv(m)=0，噪声，噪声v(n)是零均值、
50、单是零均值、单位功率的白噪声（位功率的白噪声（2v=1，mv=0)，求，求Hopt(z)和和E|e(n)|2min。解解 根据白噪声的特点得出根据白噪声的特点得出Svv(z)=1，由噪声和信号不，由噪声和信号不相关，相关，得到得到 rxx(m)=rss(m)+rvv(m)。对上式两边做对上式两边做Z变换，并代入已知条件，对变换，并代入已知条件，对x(n)进进行功率谱分解行功率谱分解:11211()()()0.361(10.8)(10.8)1.6(10.5)(10.5)()()(10.8)(10.8)xxssvvSzSzSzzzzzB z B zzz 第61页，共137页。62第四章第四章 维纳
