Chapter+4中国科学院大学现代数字信号处理课程课件.pptxVIP

第四章 功 率 谱 估 计 4.1 引言 4.2 经典谱估计 4.3 现代谱估计中的参数建模 4.4 AR模型谱估计方法 4.5 最大熵谱估计方法4.6 特征分析法谱估计4.1 引 言功率谱定义功率谱估计的方法功率谱估计的应用1、功率谱的定义信号的功率谱和其自相关函数服从一对傅里叶变换关系2、功率谱估计的方法经典谱估计方法间接方法：BT法直接方法：周期图法现代谱估计方法参数法：ARMA模型法（AR模型、MA模型、ARMA模型）非参数法：Pisarenko谐波分解法、MUSIC方法3、功率谱估计的应用在信号处理的许多场所，要求预先知道信号的功率谱密度(或自相关函数)；常常利用功率谱估计来得到线性系统的参数估计；从宽带噪声中检测窄带信号。 4.2 经 典 谱 估 计BT法周期图法改进的周期图法4.2.1 BT法BT法是先估计自相关函数, 然后进行傅里叶变换得到功率谱。有偏自相关函数估计的误差相对较小，是一种渐近一致估计：4.2.2 周期图法 周期图法的定义如下: 1. 周期图与BT法的等价关系 令 m=k-n, 即k=m+n，则 利用有偏自相关函数的BT法和周期图法是等价的。 2. 周期图法谱估计质量分析 1) 周期图的偏移式中 上式在频域表示为： 式中 周期图的统计平均值等于它的真值卷积三角谱窗函数，因此周期图是有偏估计，但当N→∞时，wB(m)→1，三角谱窗函数趋近于δ函数，周期图的统计平均值趋于它的真值，因此周期图属于渐近无偏估计。 2) 周期图的方差 用这种方法估计的功率谱在σ2x附近起伏很大，故周期图是非一致估计。 白噪声的周期图 4.2.3 经典谱估计方法改进Bartlett平均周期图法窗口处理法平均周期图Welch法（修正的周期图求平均法）1. Bartlett平均周期图法 主要思想：对序列x(n)进行L次独立观测或将其分成L段，计算每组观测数据的周期图，再将L个周期图加和后求平均。 估计方法： 假设随机信号x(n)的观测数据区间为：0≤n≤N-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组的周期图用下式表示： 将得到的L个周期图进行平均，作为信号x(n)的功率谱估计， 公式如下： 估计效果分析： 平均周期图的估计方差是周期图的方差的1/L，L越大方差越小，功率谱越平滑；相应的，M越小，偏移越大，分辨率越低；估计的均方误差也减少； 以分辨率的降低换取了估计方差的减少，估计量的方差和分辨率是一对矛盾。平均周期图法 2、窗口处理法平均周期图主要思想：用一适当的功率谱窗函数W(ejω)与周期图进行卷积，来达到使周期图平滑的目的。 -(M-1)≤n≤M-1 估计方法：式中 那么 估计效果分析： 偏移分析： 又 可得 周期图的窗函数法仍然是有偏估计, 其偏移和wB(m)、w(m)两个窗函数有关。 如果w(m)窗的宽度比较窄，M比N小得多，这样|m|N，则wB(m)~1, 由于w(m) 比wB(m)窄， W(ejw) 的主瓣比WB(ejw)宽，故可以利用窗函数法进一步平滑周期图，减小估计方差；但相应的会增加偏移，降低频率分辨率。3. 修正的周期图求平均法（Welch法）主要思想：对Bartlett法进行修正，使之更适合FFT计算。选择适当的窗函数w(n)，并在周期图计算前直接加进去；在分段时，可使各段之间有重迭，这样将会使方差减小。 估计方法：首先把数据长度为N的信号x(n)分成L段，每一段数据长度为M，N=LM；然后把窗函数w(n)加到每一个数据段上， 求出每一段的周期图，形成修正的周期图；再对每一个修正的周期图进行平均。i=1, 2, 3, …，L 第i段的修正周期图为 式中 同样, 将每一段的修正的周期图之间近似看成互不相关, 最后功率谱估计为 估计效果分析：对上式求统计平均, 得到 式中 估计是渐近无偏的； 这种方法对窗函数没有限制，无论什么样的窗函数均可使谱估计非负； 分段时，为了减少因分段数增加给分辨率带来的影响，可使各段间有重叠，例如重叠50%。 4.3 现代谱估计中的参数建模存在问题：经典谱估计中，采用观测到的N个样本值估计功率谱，认为在此观察到的N个数据以外的x(n)=0，这与实际情况是不符合的，从而造成分辨率的降低。解决方法：根据已观察到的数据，选择一个正确的模型，认为x(n)是白噪声通过此模型产生的，那么就不必认为N个以外的数据全为零了。 方法步骤：（1）选择合适的信号模型；（2）根据x(n)有限的观测数据，或者它的有限个自相关函数估计值，估计模型参数;（3）基于此模型，计算输