一、BT法(间接法)

BT法是先估计自相关函数,然后进行傅里叶变换得到功率谱。

有偏自相关函数估计的误差相对较小,是一种渐近一致估计(渐近无偏估计):

自相关图和偏自相关图怎么判断自回归性_数据

二、周期图法(直接法) 

整个计算过程中没有计算样本的自相关,拿到随机信号的样本后,直接对样本进行傅里叶变换,跳过了样本自相关计算的过程。

周期图法的定义如下:

自相关图和偏自相关图怎么判断自回归性_自相关图和偏自相关图怎么判断自回归性_02

1、周期图与BT法的等价关系

自相关图和偏自相关图怎么判断自回归性_java_03

 令m=k-n,即k=m+n,则

自相关图和偏自相关图怎么判断自回归性_java_04

 2、周期图法谱估计质量分析

周期图是有偏估计,当N→oo时,周期图的统计平均值趋于它的真值,因此周期图属于渐近无偏估计。

周期图估计的功率谱起伏很大,是一种很差的功率谱估计方法。

自相关图和偏自相关图怎么判断自回归性_java_05

 

自相关图和偏自相关图怎么判断自回归性_自相关图和偏自相关图怎么判断自回归性_06

三、经典谱估计方法改进 

存在问题:BT法和周期图法估计功率谱都不是一致估计,频率分辨率低。

解决方法:对周期图进行修正,使其满足一致估计条件。可以采用平滑处理的方法,使其方差减小。

Bartlett平均周期图法

窗口处理法平均周期图

Welch法(修正的周期图求平均法)

1、Bartlett平均周期图法(功率谱加三角窗)

主要思想:对序列x(n)进行L次独立观测或将其分成L段,计算每组观测数据的周期图,再将L个周期图加和后求平均。

自相关图和偏自相关图怎么判断自回归性_java_07

估计方法:假设随机信号x(n)的观测数据区间为:0≤n≤M-1,共进行了L次独立观测,得到L组记录数据,每一组记录数据用xi(n),i=1,2,3,…,L表示;

或对长为N的数据x(n)分成L段,每段有M个数据,N=LM,第i段数据表示为xi(n)=x(n+iM-M)。

第i组的周期图用下式表示: 

自相关图和偏自相关图怎么判断自回归性_数据_08

 将得到的L个周期图进行平均,作为信号x(n)的功率谱估计,公式如下:

自相关图和偏自相关图怎么判断自回归性_java_09

 

自相关图和偏自相关图怎么判断自回归性_开发语言_10

估计效果分析

偏移分析:

自相关图和偏自相关图怎么判断自回归性_数据_11

平均周期图仍然是有偏估计,偏移和每一段的数据个数M有关;

偏移的大小反映分辨率的高低。

方差分析:

自相关图和偏自相关图怎么判断自回归性_开发语言_12

平均周期图的估计方差是周期图的方差的1/L,L越大方差越小,功率谱越平滑;相应的,M越小,偏移越大,分辨率越低;

估计的均方误差也减少;

以分辨率的降低换取了估计方差的减少,估计量的方差和分辨率是一对矛盾。 

自相关图和偏自相关图怎么判断自回归性_开发语言_13

 

自相关图和偏自相关图怎么判断自回归性_开发语言_14

2、窗口处理法平均周期图

主要思想:用一适当的功率谱窗函数W(e^jw)与周期图进行卷积,来达到使周期图平滑的目的。

自相关图和偏自相关图怎么判断自回归性_java_15

自相关图和偏自相关图怎么判断自回归性_自相关图和偏自相关图怎么判断自回归性_16

估计效果分析: 

偏移分析:

自相关图和偏自相关图怎么判断自回归性_java_17

周期图的窗函数法仍然是有偏估计,其偏移和wB(m)、w(m)两个窗函数有关。 

 

自相关图和偏自相关图怎么判断自回归性_java_18

 

自相关图和偏自相关图怎么判断自回归性_java_19

 

自相关图和偏自相关图怎么判断自回归性_开发语言_20

 3、修正的周期图求平均法(Welch法:分段+窗函数)

主要思想:对Bartlett法进行修正,使之更适合FFT计算。

选择适当的窗函数w(n),并在周期图计算前直接加进去;

在分段时,可使各段之间有重迭,这样将会使方差减小。

估计方法:首先把数据长度为N的信号x(n)分成L段,每一段数据长度为M;

然后把窗函数w(n)加到每一个数据段上,求出每一段的周期图,形成修正的周期图;

再对每一个修正的周期图进行平均。

自相关图和偏自相关图怎么判断自回归性_java_21

 

自相关图和偏自相关图怎么判断自回归性_开发语言_22

 

自相关图和偏自相关图怎么判断自回归性_java_23

结论:传统的功率谱估计方法无论采取哪一种改进方法,总是以减少分辨率为代价,换取估计方差的减少,提高分辨率的问题无法根本解决。

自相关图和偏自相关图怎么判断自回归性_java_24

 

自相关图和偏自相关图怎么判断自回归性_开发语言_25

 

自相关图和偏自相关图怎么判断自回归性_开发语言_26

 

自相关图和偏自相关图怎么判断自回归性_自相关图和偏自相关图怎么判断自回归性_27