20_经典谱估计与参数模型法

《20_经典谱估计与参数模型法》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第四章功率谱估计§4.1引言研究信号与系统有两类基本方法:时域法和频域法.由于随机信号不存在傅氏变换,因此在频域主要研究它的功率谱.功率谱的两种定义:（1）定义之一——维纳-辛钦定理根据维纳-辛钦定理,信号的功率谱和它的自相关函数构成一对傅氏变换:正变换(4.1.1)逆变换(4.1.2)该式建立了平稳随机序列的自相关函数与功率谱密度之间的关系.分别由Wiener(1930)和Khinchin(1934)建立.（2）定义之二式(4.1.1)中的通常根据自相关函数定义计算,设为复信号,则(4.1.3)对于平稳随机信号

2、,由于服从各态历经定理,因而集合平均可以用时间平均代替.因此上式可表示为(4.1.4)将上式代入式(4.1.1),可得的另一定义:(4.1.5)通常可能获得的观测数据为有限个,因此只能计算自相关函数和功率谱的估值,这就是功率谱估计问题.功率谱估计的两种方法:功率谱的估计方法有多种,一般可分为两大类:经典谱估计和现代谱估计.§4.2经典谱估计经典谱估计建立在传统的傅里叶变换基础之上.通常有两种方法:BT法和周期图法.4.2.1BT法1958年由布莱克曼(R.Blackman)和杜基(J.Tukey)提出.BT法先是

3、根据有限个观测数据估计自相关函数,然后再利用傅氏变换得到功率谱估值.通常采用有偏自相关函数估计(方差较小),公式为,(4.2.1)称为取样自相关函数,这种估计是渐近一致估计.对上式进行傅氏变换,得BT法的功率谱估计值为,(4.2.2)为了进一步减小估计方差,还可对上式进行加窗处理.可以证明,按上式估计是非一致估计,因而不是好的估计.取二人姓氏的第1个字母,故称BT法.又称自相关法.该式在无偏估计的基础上通过三角加窗得到的.4.2.2周期图法1.周期图定义周期图法是通过对随机数据直接进行傅里叶变换的一种功率谱估计方

4、法.1965年Cooley和Tukey提出FFT以后,使DFT的运算量显著降低,从而为周期图法的应用创造了条件.根据定义,的自相关函数可表示为与的卷积,即据此,式(4.2.1)可写成,上式的傅氏变换用符号表示,则(4.2.3)其中,是的点DFT:(4.2.4)由于是周期函数,所以用式(4.2.3)估计,称为“周期图法”,可用FFT实现.实际上,根据功率谱的第二种定义式:若忽略求统计平均的运算,并设观测数据为:,所得结果与式(4.2.3)相同.2.周期图与BT法的等价性BT估计算法可归结为:由随机序列,首先求得有偏

5、自相关函数估计,然后求的傅氏变换:周期图法可归结为:由随机序列,首先求的傅氏变换再求周期图即取多个样本傅氏变换幅值模的平方,再求平均值作为功率谱估计.利用以上关系式可以证明,周期图就是BT法中有偏自相关估值所得的功率谱估计.●证明如下:由周期图表达式令,即,则由于令m=k-n:(1)k=0m=-nn=0N-1m=-(N-1)0(2)k=N-1m=N-1-nn=0N-1m=0(N-1)综合(1),(2)可得m的数值范围是:-(N-1)≤m≤N-1有偏估计所以周期图法不需要估计自相关函数,可以利用FF

6、T进行计算,因此比BT法简单,在出现FFT以后,应用较为广泛.即同期图与采用有偏自相关估计的BT法是一致的.4.2.3经典谱估计的问题与改进1.主要问题经典谱估计的基础是傅氏变换.因傅里叶变换域为无穷大,而观测数据是有限的,对观测不到的数据实际上都被强制地当作0处理,这相当于无限长样本用矩形窗加以截断.由于窗口以外的数据仍有很强的相关性,因此,用有限长样本序列估计出来的功率谱必然存在很大的偏差.有限长序列的傅氏变换,等于无限长序列与矩形窗函数各自傅氏变换的卷积.矩形窗函数的频谱是sinc函数,卷积结果带来的影响是

7、:（1）sinc函数主瓣的影响——分辨率下降sinc函数的主瓣将引起信号功率谱向邻近频域扩展(即“频谱泄漏”),使谱的主瓣展宽,谱的分辨率下降.（2）sinc函数旁瓣的影响——谱间干扰谱间干扰的后果是,强信号的功率谱旁瓣会影响弱信号功率谱的检测，造成谱的模糊甚至畸变,严重时可能掩盖弱信号,或者把旁瓣误认为是信号,以至造成假信号.由于在采用自相关函数的有偏估计时,周期图法即是BT法,因此,当序列为有限长时,BT法与周期图法及其改进方法,都具有相同的问题.2.改进方法（1）平均周期图法对随机信号进行组独立观测,一次观

8、测得到个记录数据.共有个观测数据.其中第组的周期图为(4.2.5)的功率谱是个独立周期图的平均,即(4.2.6)改进效果:的估计方差是的.（2）窗函数法选择适当的窗函数对周期图进行平滑.设点序列周期图为,窗口的谱函数为,加窗后的功率谱是与的卷积,即(4.2.7)L个周期图之间近似互不相关.其中具有低通特性.上式等效于对周期图进行频域修正,即让周期图通过一个线性非频变系统,