现代数字信号处理(chap3确定性最小二乘)(修订版)ppt课件

1、主讲教师主讲教师: 何松华何松华 教授教授联络方式联络方式: 73132618 (0731) 82687718 73132618 现代数字信号处置/确定性最小二乘滤波器第三章第三章 确定性最小二乘滤波器确定性最小二乘滤波器本章的教学内容本章的教学内容正那么方程正那么方程滤波器的渐近性滤波器的渐近性最小二乘逆滤波器最小二乘逆滤波器白化滤波器白化滤波器前前 言言一、一、LP、HP、BP、BS 滤波器设计滤波器设计(本本科科)设计特点：在频域上给出容限图，期望能逼近理想设计特点：在频域上给出容限图，期望能逼近理想滤波器。滤波器设计与输入信号的特性关系不亲密。滤波器。滤波器设计与输入信号的特性关系不亲

2、密。以低通滤波器为例前前 言言二、最正确滤波器的概念二、最正确滤波器的概念期望输出期望输出实践输出实践输出 设计特点：在时域上，我们希望实践输出尽量逼近设计特点：在时域上，我们希望实践输出尽量逼近期望输出。与输入和输出信号的波形特性亲密相关。期望输出。与输入和输出信号的波形特性亲密相关。 输入和输出既可是随机信号，也可是确定信号。输入和输出既可是随机信号，也可是确定信号。 逼近准那么：最小均方误差、最小二乘逼近准那么：最小均方误差、最小二乘h(k)g(k)y(k)f(k)以信道平衡为例?前前 言言确定性最小二乘滤波器确定性最小二乘滤波器2( )min ( )( )( )h knf kg kh

3、k直接对样本数据g(k)进展处置,无须知道g(k)的统计特性,只对该样本现实最正确统计性最小二乘滤波器统计性最小二乘滤波器(最小均方误差最小均方误差)2( )min ( )( )( ) h kEf kg kh k需求知道g(k)的二阶统计特性,对具有该二阶统计特性的随机信号的一切样本现实,从平均意义上最正确第4章前前 言言滤波滤波预测预测平滑的概念平滑的概念期望输出期望输出实践输出实践输出预测滤波器(k)y(k)f(j+k)j0期望输出期望输出实践输出实践输出平滑滤波器(k)y(k)f(j+k)j0gg2( )min ()( )( )h knfjkg kh k( )h k( )h k第一节第一

4、节 正那么方程正那么方程一、最正确滤波的引出一、最正确滤波的引出y = f = g* h)()()(zHzGzF)()()(zGzFzH假设输入信号假设输入信号g g为单位脉冲信号为单位脉冲信号u(k)u(k)经过线性系统经过线性系统G(z)G(z)产生产生1( )( )( )F zh kZG z问题的处理似乎很容易,真有这么简单吗?第一节第一节 正那么方程正那么方程在因果、稳定的要求下，在因果、稳定的要求下，H(z)H(z)无法物理实现的缘由：无法物理实现的缘由： ()y = g* hff最优准那么：最优准那么：LMSLMS准那么、最小准那么、最小二乘二乘 的零点不一定在单位圆内的零点不一定

5、在单位圆内, , 假设假设 不是最小相不是最小相位的，那么位的，那么 就不是稳定的；就不是稳定的； 对应的能够是一个因果稳定的对应的能够是一个因果稳定的IIRIIR滤滤 波波器，而所设计的器，而所设计的 要求是一个有限阶的要求是一个有限阶的FIRFIR滤滤 波器。波器。 假 设假 设 不 是 因 果 的不 是 因 果 的 ( ( 平 滑 滤 波平 滑 滤 波 ) ) ， 那 么， 那 么 不是因果的不是因果的)(zG1/( )G z)(/ )(zGzF)(zH( )F z)(/ )(zGzF物理可实现的处理方案：设计因果、稳定或物理可实现的处理方案：设计因果、稳定或FIRFIR的的)(zH)(

6、zG最正确线性滤波第一节第一节 正那么方程正那么方程二、正那么方程二、正那么方程22( )lf k f( )0,0 ()h kkor kne = f - y = f - g* h误差序列误差序列222( )( )( )( )kVf kg kh khf - yf - gh21( )()()()2jjjVF eG eH edh22( )lg k g因果输入序列因果输入序列期望输出期望输出n阶阶FIR滤波器滤波器(不限定不限定)表示平方可积(和)序列帕斯瓦尔定理误差总能量第一节第一节 正那么方程正那么方程最小二乘准那么最小二乘准那么 (LS : Least-Square) (LS : Least-S

7、quare) min( )Vhh220nkkmV hf ky kf kh m g km( )( )( )( )() ()20( )2( )( ) ()nkmkfkh mf k g km00( ) ( )() ()nnlmkh m h lg km g kl要求h(n)因果、FIR积分的乘积变成二维积分，采用两个积分变量第一节第一节 正那么方程正那么方程22( ) |kfkf与与h h无关的常数项无关的常数项 ： 期望输出信号的能量期望输出信号的能量 线性项的系数线性项的系数 ：输入与输出的相互关函数：输入与输出的相互关函数 0( )( ) ()( ) ()( ) ()kkkq mf k g km

8、g k f kmg k f km1()()2jjjmF eG eed0nm输入序列是因果的卷积定义？( )()f mgm() (g k mgm k 第一节第一节 正那么方程正那么方程二次项的系数二次项的系数 输入信号自相关输入信号自相关 ( , )() ()kr m lg km g kl00( ) ()( ) ()kkg k g kmlg k g klm 2()1()2jj m lG eed()()r mlr lm积分变量置换kmk积分变量置换klk 0nm0nl 参见前面求和表达式( ) ()g mgml第一节第一节 正那么方程正那么方程(0)(1)( )Thhh nh(0)(1)( )Tq

9、qq nq(0)(1)( )(1)(1)( )(1)(0)rrr nrRrr nrrToeplitz矩阵矩阵 (n+1)(n+1)维的输入自相关矩阵，第i行第j列元素值为r(|i-j|)n+1维滤波矢量令n+1维输入输出相互关矢量0( )( ) ()kq mg k f km0( )( ) (|)kr mg k g kmn阶滤波器第一节第一节 正那么方程正那么方程2TTV(h)q h+ h Rh20( )2( )( ) ()nkmkfkh mf k g kmV(h)00( ) ( )() ()nnlmkh m h lg km g kl( )q m()r ml上式用矢量、矩阵方式表示为上式用矢量、

10、矩阵方式表示为练习:反过来验证第一节第一节 正那么方程正那么方程2 TTV(h)q h+ h Rhmin( )0Vhh标量V对矢量h求偏导零矢量2 0 q+ 2Rh当R为对称矩阵时1 = hR q下面引见另外一种相对复杂的推导方法(展开法)第一节第一节 正那么方程正那么方程( )00,( )Vlnh lh0( )( ) ()()00nkmf kh m g km g klln 0(|) ( )( )0nmr lm h mq lln 估计误差序列与估计误差序列与输入序列输入序列(l=0)及及其平移序列正交其平移序列正交2 ( ) ()00ke k g klln 220nkkme kf kh m g

11、 km( )( )() ()V(h)参见r()、q()的定义及性质第一节第一节 正那么方程正那么方程(0)(1)( )(0)(0)(1)(1)( )(1)(0)( )( )rrr nhqrrr nrrh nq nRhq1Rhq0( )( ) ()kr mg k g km0( )( ) ()kq mg k f km最小二乘滤波器的正那么方程第一节第一节 正那么方程正那么方程三、误差分析三、误差分析R hq( )2TTVRhq hhhmin( )2TTVhq hh qTTRq hhh200( ) ( )() ()( )nnlmkkh m h lg km g klyk222( )( )( )0kkk

12、e kfkyk 即：期望输出能量一定大于最小二乘滤波器实践输即：期望输出能量一定大于最小二乘滤波器实践输出能量。类似于勾股定理。出能量。类似于勾股定理。误差能量等于期望输出信号能量-实践输出信号能量练习:根据定义以及可分别的二维求和性质第一节第一节 正那么方程正那么方程相对误差能量：相对误差能量：2( )( )( )kV hV hfk2minmin22( )( )11( )( )TkkkykV hq hfkfk 10minmin0期望输出期望输出 f f 与实践输出与实践输出 y y 完全一致完全一致1min期望输出期望输出 f f与实践输出与实践输出 y y完全不一致完全不一致输入与期望输出

13、完全不相关。输入与期望输出完全不相关。 第一节第一节 正那么方程正那么方程000hqy0g1gfyfeff输入序列长度为2情况下的举例，大于2时，引入超平面子空间概念期望输出与输入不相关期望输出与输入完全相关为输入的线性组合，必在输入信号构成的子空间内yf 窄带干扰信号y(k)有长的相关长度，自相关在范围 上有较大的值；干扰数目、频率等参数未知。例：窄带干扰例：窄带干扰(NBI)消除消除( )( )( )( )x ks ky kv k( )0sr llD0,lMMDs(k)、y(k)和v(k)相互之间不相关；热噪声v(k)是白色的；有用信号s(k)是宽带的，因此有短的相关长度，即 第一节第一节

14、 正那么方程正那么方程超宽带雷达(探地、救灾)信号处置,如何根据观测信号x(k)得到回波信号s(k)的最正确估计,或如何消除干扰y(k)例如各种窄带的通讯信号由于由于y(k)和和x(k)是相关的，可以经过是相关的，可以经过x(k)用最正确线用最正确线性估计方法得到一个对性估计方法得到一个对NBI的估计：的估计： *yx h( )( )y ky k假设假设 那那么么( )( )( )( )x ky ks kv k第一节第一节 正那么方程正那么方程线性预测器可以利用最小二乘滤波器来实现。滤波线性预测器可以利用最小二乘滤波器来实现。滤波器期望输出为器期望输出为y(k)，输入为，输入为x(k) Rhq

15、问题：只需观测数据问题：只需观测数据 x(k) ，y(k)不可得。不可得。序列的卷积如何计算相互关?可计算自相关?思绪：最正确滤波器的设计需求的是期望输出与输入思绪：最正确滤波器的设计需求的是期望输出与输入 之间的相互关函数之间的相互关函数q(m)|m=0,1,n,而不是而不是 期望输出期望输出f(k);进一步的问题是进一步的问题是: 依然无法计算依然无法计算 y(k)与与x(k)的相互关的相互关思索：思索：x(k)由由s(k)、y(k)、v(k)三部分组成，三部分组成， s(k)的相的相 关长度小于关长度小于D， y(k)的相关长度大于的相关长度大于D， v(k) 相关长度为零，那么相关长度

16、为零，那么x(k)中的中的s(k) 与与x(k-D)中的中的 s(k-D)、y(k-D)、v(k-D)三部分都不相关三部分都不相关;同理同理, x(k)中的中的v(k) 也与也与x(k-D)不相关不相关,那么那么 x(k)与与x(k-D) 的自相关就是期望输出的自相关就是期望输出y(k)与输入与输入 x(k-D)的互的互 相关相关第一节第一节 正那么方程正那么方程思绪：以思绪：以x(k-D)作为输入，以作为输入，以y(k)作为期望输出；那么作为期望输出；那么y(k)与与x(k-D) 的相互关函数是可计算的的相互关函数是可计算的转化为转化为x(k)的自相的自相关关第一节第一节 正那么方程正那么方

17、程1100( )( ) ()( ) ()kkq mx k y kmx k x km1( )()()()()x kx kDs kDy kDv kD( )y nZ-D前向线性预测+-x(k)e(k)性质:1100( ) ()( ) ()()()kkx k x kmx k s kmy kmv km证:第一节第一节 正那么方程正那么方程11000( )( ) ()( ) () ( )() ()( ) ()kkk Dkq mx k x kmx k y kmq mx kD x kmx k x kDmm0时，时， 与与 互不相关，那么互不相关，那么 与与()s kD()s km()s km 互不相关，同理，

18、互不相关，同理， 与与 互不相关，那么互不相关，那么1( )x k()v km1( )x ky(k)未知情况下依然可以经过x(k)及其延时序列x1(k)估计y(k)与输入序列的相互关函数11100( )( )() () () ( ) ()( )xkk Dxkrmx k x kmx kD x kDmx k x kmr m第一节第一节 正那么方程正那么方程10( )( ) () (0,1,., )N D mkq mx k x kDmmn设设x(k)的实践长度范围为的实践长度范围为x(k)|k=0,1,N-1，滤波器，滤波器的长度范围为的长度范围为h(k)|k=0,1,n (nN-D)；计算；计算1

19、0( )( ) () (0,1,., )N mkr mx k x kmmn构造如下矢量与矩阵(0)(1)( )Thhh nh(0)(1)( )Tqqq nq第一节第一节 正那么方程正那么方程那么有：(0)(1)( )(1)(1)( )(1)(0)rrr nrRrr nrr1R qhmax(0,)( )( ) () (0,1,.,1)kmk ny kx m h kmkN( )( ) (0,1,.,1)x ky kkN去除窄带干扰后的信号为第一节第一节 正那么方程正那么方程(2)/(2)/1( ) (08)frtnsttnstg ttnsee 2.3,0.4,2,2rfstns tnsfGHz(

20、)( /)(1)/ (0,1,2,.,15)sss kg kfg kfk141( )sin()iiiiy kAk5 (96)( )( )96,97,.,111( )( )( )96,111s ky kv kkx ky kv kk( ) (0, 0.1)v kWGN举例(0,1,2,.,199)k 采样间隔0.5ns第一节第一节 正那么方程正那么方程第一节第一节 正那么方程正那么方程14个干扰频率课程上机实验4：确定性最小二乘滤波器的实现上述举例中，其他参数不变，干扰信号参数设置为 A1=A2=A14=2; 1=2= 14=0 1=0.06,2=0.10,3=0.18,4=0.21 5=0.30

21、,6=0.48,7=0.52,8=0.57 9=0.61,10=0.64,11=0.67,12=0.70 13=0.78,14=0.94采用Matlab言语编程(1) 参照实验3的方法产生均值为0，根方差为0.1的正态白噪声第一节第一节 正那么方程正那么方程 ( )|1,2,.,200v kk 数据点数从实验3的100000变成200(2) 按设定的参数产生干扰信号数据第一节第一节 正那么方程正那么方程141( )sin(1) (1,2,.,200)iiiiy kAkk针对Matlab的数组下标只能从1开场的调整(3) k0=96，产生信号数据2.30.5(1) 2/0.42.30.5(1)

22、2/21( ) (1,2,.,15,16)kkg kkee( )( )(1) (2,3,.,15,16)s kg kg kk(1)(1)sg第一节第一节 正那么方程正那么方程5 (96)97,98,.,112( )01,2,.,96;113,114,.,200s kks kk画出 以及 的波形图(4) 画出含有噪声以及干扰的合成信号x(k)的波形图( )( )( )( ) (1,2,.,200)x ks ky kv kk( )s k( )( )s kv k(5) 给定参数D=16，n=12,N=200；求FIR线性最小二乘滤 波器 h(k)|k=1,2,n+1数组下标从1开场,对应实践的h(0

23、)到h(n)第一节第一节 正那么方程正那么方程11( )( ) (1) (1,2,.,1)N D mkq mx k x kDmmn11( )( ) (1) (1,2,.,1)N mkr mx k x kmmn(1)(2)(1)Tqqq nq(1)(2)(1)(2)(2)(1)(2)(1)rrr nrRrr nrr( , )(| 1)R i jr ij第一节第一节 正那么方程正那么方程(6) 利用线性最小二乘滤波器对观测数据x(k)进展滤波处置，得到窄带干扰信号序列的估计值序列(7) 利用滤波器的输出数据序列对观测信号中的干扰信号进展对消处置；画出对消后信号 的波形图；察看其与 的类似性及差别性

24、max(1,)( )( ) (1) (1,2,.,200)kmk ny kx m h kmk( )( )( ) (1,2,.,200)s kx ky kk( )s k( )( )s kv k1(1)(2)(1)Thhh nR qh利用Matlab的矩阵运算函数第二节第二节 最小二乘滤波器的渐近性最小二乘滤波器的渐近性22( )lf k f( )0,0h kk22( )lg k g输入序列，因果，不一输入序列，因果，不一定最小相位定最小相位期望输出序列，不一定因果期望输出序列，不一定因果因果因果IIR滤波器滤波器(n)下面讨论最小二乘滤波器的阶n足够大或趋于的情况220nkkme kf kh m

25、 g km( )( )() ()V(h)问题：n时，V(h)可以下降到多少？能否可下降到0？期望输出实践输出y(k)数据平滑运用第二节第二节 最小二乘滤波器的渐近性最小二乘滤波器的渐近性一、最优因果一、最优因果IIR LS 滤波器的求解问题滤波器的求解问题1min0( )( )()G zGzzz01z min( )Gz最小相位最小相位假设输入信号g(k)是因果、稳定的但非最小相位的，G(z)存在一个单位圆外的零点1/z0110min0101( )( )()1z zG zGz zzz z作等值变换矩阵求解?不行Rhq问题转化: 序列所对应的全通函数序列的概念第二节第二节 滤波器的渐近性滤波器的渐

26、近性10min0( )( ) (1)G zGzz z最小相位最小相位1010()( )1zzD zz z全通函数全通函数1*10min00*10( )(1)( )( )1zzGzz zG z D zz z练习3.1：令z=ej,z0=|z0|ej利用复数运算性质证明|()| 1jD e0( )( )( )G zG z D z对于G(z)有多个单位圆外零点的情况，采用同样方法可以得到最小相位全通函数稳定因果第二节第二节 滤波器的渐近性滤波器的渐近性1110()( )1IiiizzD zz z全通函数全通函数|()| 1jD e|0,1,.,1iziIG(z)的一切单位圆外的零点的一切单位圆外的零

27、点1( )( )d kZTD z设设0, ( )0kd k全通函数对应的全通函数对应的序列为因果的序列为因果的输入序列输入序列 g(k) 知的情况下知的情况下, d(k) 也是知的也是知的2201( )|()|12jkdkD ed序列总能量序列总能量分子分母的阶一样 第二节第二节 滤波器的渐近性滤波器的渐近性000( )( )1( )( )( )( )( )( )F zF zH zHzG zG z D zG z000001( )( )( ) ( )( )( )( ) ( )( ) * * G z H zG z D zHzD z HzF zG zg hdhf01( )G z0( )Hz( )H

28、zgfd0h0( )( )( )F zHzD z现实情况下，要求现实情况下，要求理想情况下相当于以d(k)为输入,f(k)为期望输出( )G z( )D z( )F z第二节第二节 滤波器的渐近性滤波器的渐近性22( )()minVfV00hfg hdhhdR 0dhq0(|)( ) (|)dkrmld k d kml2|1()()2jj m lD eedml0dhq00( )( ) ()() ( )0kk mh md k f kmd km f kmRd、h0、qd为无穷阶的即即RI有有的解h0满足根据功率谱定义以及全通函数的性质即即 Rd的第m行第l列元素值第二节第二节 滤波器的渐近性滤波器

29、的渐近性22min0( )( )( )TTkmVfkqmd0dddhqhqq20( )() ( )() ( )kmk ml mfkd km f kd lm f l 2000( )( ) ( )() ()kklmfkf k f ld km d lm二、误差分析二、误差分析根据第一节(三)2000( )( ) ()( ) ()kmklfkf k d kmf l d lm 定义i1G zz 试设计一个n阶的FIR滤波器，满足 ( )1/( )H zG z解：解： ) 1()(kkkg2( )( ) ()1010kr mg k g kmmmothers (0)10( )00gmq mm实数实数( )(

30、 )( )g kh kk第三节第三节 最小二乘逆滤波器最小二乘逆滤波器22(0)110(1)0( )001hhh n 2(1) (0)(1)1hh2(1)(1) ( )0h nh n对于此特殊情况,不一定要经过矩阵求逆解方程组第0个方程2(1)(1) ( )(1)0h kh kh k第k个方程第n个方程第三节第三节 最小二乘逆滤波器最小二乘逆滤波器2(1)(1) ( )(1)0h kh kh k12( )kkh kcc通解通解0kn12(1)0zz 差分方程有两个根1z21/z根据实际，差分方程的解通式思索到( )0 0 or h kkkn一切方程满足第三节第三节 最小二乘逆滤波器最小二乘逆滤

31、波器解出系数解出系数c1、c2 1(1)3(1)0( )0nknknnknh kothers min( )1(0) (0)1(0)Vghh h1(1)313(1)3(1)1nnnnnnnn 练习3.2：将上式代入方程组中的第1,n个方程求得：第三节第三节 最小二乘逆滤波器最小二乘逆滤波器渐近性：渐近性：123lim ( )0nkknnh kk 21121111lim( )1( )nH zzzG z11211111 11lim( ) ( )( )1nzzH z G zD zzz min21lim( )1 lim (0)1nnVh h全通函数第三节第三节 最小二乘逆滤波器最小二乘逆滤波器min,

32、lim( )1nV h min1, lim( )0nV h11011( )()( )( )| 1zG zzD z G zz11( )( )(1)kkd ku ku k 1(0)d22211101nkV hdkd minlim ( )( )( )也可利用渐近性进展误差能量求解：也可利用渐近性进展误差能量求解：根据Z变换性质例例221( )4G zz解：解： (1) L = 12第三节第三节 最小二乘逆滤波器最小二乘逆滤波器试设计一个10阶的FIR滤波器，使得 ( )LF zz(1) L = 12； (2) L = 0 009( )110mq mm1.06250( )0.2520mr mmothe

33、rwise 1( )( )(2)4g kkk( )(12)f kk练习3.3练习3.3续g(k)非最小相位序列10阶滤波器,只须计算0至10的q(m),虽然q(12)=1/4第三节第三节 最小二乘逆滤波器最小二乘逆滤波器4322 9.155 10 ,0, 3.891 10 ,0, 1.562 10 ,0, 6.25 10 ,0,-0.25,0,-1.0h ygh= -2.28910-4,0,-5.72210-5,0,-1.43110-5,0, -3.57610-6,0,-8.94110-7,0,-2.23610-7,0,18( )5.588 10V h(2) L = 01/40( )01,2,

34、10kq kk练习3.3续解正那么方程组练习3.3续：计算两个有限长度序列的卷积k=12时 最大误差很小( )y k( )(12)y kk( )( )f kk第三节第三节 最小二乘逆滤波器最小二乘逆滤波器ygh=1/16,0,-0.234,0,-0.0586,0,-0.0146,0,-0.00366, 0,-0.00092,0,-0.00023 ( )0.9373V h h =1/4,0,6.2510-2,0,-1.56210-2,0, 3.90410-3, 0,9.72410-4,0,2.28810-4y(k)与(k)之间的误差很大两种情况的比较1( )( )(2)4g kkk输入期望输出(

35、 )(12)f kk群时延为正，用因果滤波器实现时，误差较小练习3.3续第三节第三节 最小二乘逆滤波器最小二乘逆滤波器例例3 输入输入:14 . 011)(zzG期望输出期望输出:14 . 01)(zzF求四阶求四阶( ) ( )0,04optHzh kkor k解：解：-2( )( )1-0.16( ) ( )optF zH zzHzG zmin( )0Vh1( )( )(2)4g kkk输入期望输出( )( )f kk群时延为负，且输入序列为非最小相位的，用因果滤波器呼应在输入之后无法实现相位的补偿，误差较大显然，输入为最小相位的，输出为因果的；n足够大时，存在准确的解第四节第四节 确定性

36、白化滤波器确定性白化滤波器生成模型生成模型白化模型白化模型一、白化滤波器引出一、白化滤波器引出假设对假设对H(z)=1/G(z)无任何限制，平稳随机信号经过无任何限制，平稳随机信号经过该系统后一定能得到完全白化；该系统后一定能得到完全白化；假设假设H(z)为因果、稳定为因果、稳定IIR系统，那么系统，那么G(z)必需为必需为最小相位才干得到完全白化；最小相位才干得到完全白化；假设要求假设要求H(z)为为n阶阶FIR系统，此时系统，此时 H(z)是是G(z)的的n阶最小二乘逆滤波器阶最小二乘逆滤波器, 不一定能完全白化。不一定能完全白化。G(z)(n)nx不思索物理可实现性1/G(z)nx(n)

37、二、第一种解法二、第一种解法 逆滤波器逆滤波器2( ) |minVhg h(0)(1)( )(0)(0)(1)(1)0(1)( )(1)(0)( )0rrr nhgrhrr nrrh n1( )( )H zG z第四节第四节 确定性白化滤波器确定性白化滤波器min( )1(0) (0)Vgh h正那么方程正那么方程三、第二种解法三、第二种解法白化滤波器白化滤波器 与逆滤波器本质上一致与逆滤波器本质上一致第四节第四节 确定性白化滤波器确定性白化滤波器用归一化的n阶因果滤波器a (k)a(0)=1(k=0,1,n)对信号进展滤波处置，期望输出为(k) (未知,白化,但1)2min220222022

38、( )|( )*| (0) (0)( ) 2(0) (0) (0)( ) 2(0) |*|kkVkgagy kagyy kggaaa逆滤波器的h(0)不一定为12|minag那么( )minVa2( ) |minVaag(0)1a220( )*( ) ()nkmVa m g km aag( )0 1,2,( )Vlna la00( ) () ()( )() ()0nnkmmka m g km g kla mg km g kl 第四节第四节 确定性白化滤波器确定性白化滤波器(|)r ml正二次型,具有最小值点0( ) (|)0 1,2,(0)1nma m r mllna第四节第四节 确定性白化滤波器确定性白化滤波器(0)(1)(1)(1)(1)(1)(2)(2)(1)(1)(1)(0)( )( )rrr narrarrr nrra nr n1( ) (|)(|)( ) 1,2,nma m r mlrlr lln 根据上述方程组可以求得a(1),a(2),a(n)第四节第四节 确定性白化滤波器确定性白化滤波器00000( )( )() ()( )(