实验六功率谱估计.docVIP

实验六 功率谱估计 1、实验要求： （1）理解功率谱估计的基本概念； （2）掌握经典功率谱估计方法——直接法和间接法； （3）掌握改进的经典功率谱估计方法，例如Welch法。 2、实验原理： 功率谱估计就是基于有限的数据寻找信号、随机过程或系统的频率成分。它表示随机信号频率域的统计特性。随机信号是无始无终具有无限能量的，所以其傅立叶变换并不存在，因为它不满足绝对可积的条件。因此需要研究其在频率域上的功率分布情况，即功率谱密度或功率谱。 经典功率谱分析分为直接法和间接法。直接法，即周期图法，是由傅立叶变换得到的：将随机信号的N点样本值看作能量有限信号，取其傅立叶变换，得到；然后再取其幅值的平方，并除以N作为的真实功率谱的估计，即 间接法，又称为自相关法或BT法，是由随机信号N个观察值，估计出自相关函数，然后再求的傅立叶变换作为功率谱的估计： 直接法和间接法的方差性能很差，而且当数据长度太大时，谱曲线起伏加剧；若数据长度太小，则谱的分辨率又不好，所以需要改进。改进的直接谱估计方法由Bartlett法和Welch法。 Bartlett法将采样数据分成L段，每段的长度都是M，即N LM，对每段数据加矩形窗，再计算其各自的功率谱，把对应相加，再取平均，得到平均周期图。Welch法是对Bartlett法的改进：一，在对分段时，可允许每段数据有部分重叠；二，每段数据窗口可以不是矩形窗口，例如使用汉宁窗或哈明窗，记为。然后按Bartlett法求每一段的功率谱，记为 ，其中。平均后的功率谱为：（段数） 3.程序说明 （1）zhouqitu.m 最普通的周期图法 （2）bartlett.m 分段平均周期图法 Bartlett法 （3）welch.m 加窗平均周期图法 Welch 求信号的功率谱。其中， 0.05, 0.120,w n 为白噪声，信号长度N 1024。用Welch法求信号的功率谱密度估计。（注：本处采用2：1重叠分段） 分析：本处窗函数用的是hanning 窗对信号进行处理，这样可以减小“频谱泄漏”，同时增加频峰的宽度，从而提高了频谱分辨率。