功率谱估计方法的比较

《功率谱估计方法的比较》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、.-功率谱估计方法的比拟摘要：本文归纳了信号处理中关键的一种分析方法,即谱估计方法。概述了频谱估计中的周期图法、修正的协方差法和伯格递推法的原理，并且对此三种方法通过仿真做出了比照。关键词：功率谱估计；AR模型；参数引言：谱估计是指用已观测到的一定数量的样本数据估计一个平稳随机信号的谱。由于谱中包含了信号的很多频率信息，所以分析谱、对谱进展估计是信号处理的重要容。谱估计技术开展渊源很长，它的应用领域十分广泛，普及雷达、声纳、通信、地质勘探、天文、生物医学工程等众多领域，其容、方法都在不断更新，是一个具有强大生

2、命力的研究领域。谱估计的理论和方法是伴随着随机信号统计量及其谱的开展而开展起来的，最早的谱估计方法是建立在基于二阶统计量,即自相关函数的功率谱估计的方法上。功率谱估计的方法经历了经典谱估计法和现代谱估计法两个研究历程，在过去及现在相当长一段时间里，功率谱估计一直占据着谱估计理论里的核心位置。经典谱估计也成为线性谱估计，包括BT法、周期图法。现代谱估计法也称为非线性普估计，包括自相关法、修正的协方差法、伯格〔Burg〕递推法、特征分解法等等。原理：经典谱估计方法计算简单，其主要特点是谱估计与任何模型参数无关，是

3、一类非参数化的方法。它的主要问题是：由于假定信号的自相关函数在数据的观测区间以外等于零，因此估计出来的功率谱很难与信号的真实功率谱相匹配。在一般情况下，经典法的渐进性能无法给出实际功率谱的一个满意的近似，因而是一种低分辨率的谱估计方法。现代谱估计方法使用参数化的模型，他们统称为参数化功率谱估计，由于这类方法能够给出比经典法高得多的频率分辨率，故又称为高分辨率方法。下面分别介绍周期图法、修正的协方差法和伯格递推法。修正的协方差法和伯格递推法采用的模型均为AR模型。..word.zl-.-〔1〕周期图法周期图法是

4、先估计自相关函数,然后进展傅里叶变换得到功率谱。假设随机信号x(n)只观测到一段样本数据，n=0,1,2,…,N-1。根据这一段样本数据估计自相关函数，如公式〔1〕〔1〕对〔1〕式进展傅里叶变换得到〔2〕式。〔3〕〔2〕如果忽略上式中求统计平均的运算，假设观测数据为:x(n)0≤n≤N-1，便得到周期图法的定义式〔3〕:〔3〕〔2〕修正的协方差法修正协方差法使用前向和后向预测误差平均值最小的方法，估计AR模型的参数，从而得到信号的功率谱。信号的前向和后向预测分别公式〔4〕，〔5〕：〔4〕..word.zl-.

5、-〔5〕式中apk是AR模型的参数。前向和后向预测误差功率ρpe、ρpb分别用〔6〕，〔7〕式表示〔6〕〔7〕〔8〕最小预测误差平均功率是模型输入白噪声的方差，即ρp=σ2w，前、后向预测误差平均功率为式〔8〕为了使预测误差平均功率最小，求ρp对apk(k=1,2,3,…,p)的微分,或者用复梯度法求，得到式〔9〕〔9〕..word.zl-.-化简并写成矩阵形式为式〔10〕：〔10〕协方差函数白噪声的方差估计值为式〔11〕：〔11〕..word.zl-.-观测数据x(n)(n=0,1,2,…,N-1)，利用上

6、面公式可以求出模型的参数：{api(i=1,2,3,…,p);w}。式中的协方差函数cxx(j,k)，有两个变量，因此也适合于非平稳随机信号。〔3〕伯格递推法设信号x(n)观测数据区间为：0≤n≤N-1，前向、后向预测误差功率分别用ρp,e和ρp,b表示，预测误差平均功率用ρp表示，公式分别为〔12〕，〔13〕，〔14〕〔12〕〔17〕〔13〕〔14〕前向、后向预测误差递推公式如式〔15〕：〔15〕将式〔15〕带入〔12〕，〔13〕，〔14〕公式中，得到式〔16〕..word.zl-.-〔16〕求预测误差平均

7、功率ρp最小时的反射系数kp，令式〔17〕为零〔17〕得式〔18〕就是利用伯格递推法求第p个反射系数的公式〔18〕实验：对高斯噪声中的信号x(t)=cos(2π*60*t)+cos(2π*63*t)进展功率谱估计，采样频率为fs=1000Hz,分别用周期图法、改良的协方差法和伯格递推法估计序列的功率谱，且AR模型的阶数取30与50两种情况，分别对应图1，图2，图3，图4，图5。实验结果：..word.zl-.-图1周期图法估计功率谱图230介改良的协方差法估计功率谱..word.zl-.-图350介改良的协方

8、差法估计功率谱图430介伯格递推法法估计功率谱..word.zl-.-图550介伯格递推法法估计功率谱实验结论：周期图法功率谱估计的分辨率低于修正的协方差法和伯格递推法，修正的协方差法估计功率谱分辨率与伯格递推法估计功率谱分辨率相当，但是修正的协方差法需要先由信号的观测数据估计自相关函数，这是它的缺点，而伯格递推法那么由信号的观测数据直接计算AR模型参数。而且增加AR模型的阶数可以提高分辨率。完毕语