周期图法功率谱估计窗口效应.pptVIP

周期图法功率谱估计------窗口效应 专业： 学号： 姓名： 主要内容 概述 平均周期图法 窗函数法 总结 概述 功率谱估计就是用有限长的数据来估计信号的功率谱。 周期图法：利用FFT算法直接对采样数据序列进行傅立叶变换来估计功率谱的方法。 BT谱估计法：利用维纳相关法从采样数据序列的相关函数来得到功率谱估计的方法。 功率谱估计 BT谱估计法 周期图法 平均周期图法 平滑平均周期图法 二、周期图法 平均周期图法 巴特利特（ Bartlett ）提出了一种通过对若干独立估计值进行平均，来降低周期图估计量的方差的方法。 这种改进方法就是将数据分段，各分段数据之间没有重叠，分别求出每一段的功率谱，然后加以平均。 平滑平均周期图法 另一种方法是选择合适的窗函数 与周期图卷积而使周期图平滑，故称为窗函数法。 巴特利特平均周期图法的实现过程就是把数据分成互不重叠的L 段，每段有M个样本，即 其中 , 是长度为M的矩形窗。令 ,则第i段的周期图可有下面公式求得： 由上式可求得整个序列的谱估计： 巴特利特平均周期图法的仿真 在巴特利特平均周期图法中，对于固定长度的N，分段K数目愈大，方差愈小，但M愈小，会使偏倚增大，谱分辨率变得愈差。 因此在应用巴特利特平均周期图法时，需要在偏倚和估计方差之间进行权衡，根据实际要求选择适当的M和N。 三、窗函数法(平滑平均周期图法) 加窗处理的宗旨是减小频谱泄漏，但却降低了频谱分辨率。一般说来，减小泄漏和提高分辨率是互相矛盾的，但是如何使两者都取得最好的效果，就必须对窗函数的特点进行分析，然后根据实际要求来选择。其公式如下： 式中窗函数序列w(m)的长度为2M-1的有限宽窗序列，为了保证 是一个实、偶、非负函数， w(m)应是一个偶序列，并且满足条件 一种改进的窗处理是韦尔奇（Welch）法，因为窗在两边渐变为零，所以这种方法降低了由于重叠导致的段间统计相依的效应，而且，一个合适的非矩形窗可以减小“旁瓣效应”，即“谱泄漏”。 经常采用的窗函数有矩形窗、汉宁窗(Hanning)、海明窗(Hamming)、凯瑟窗(Kaiser)和切比雪夫窗(Chebwin)等，其中矩形窗、汉宁窗以及海明窗都是广义余弦窗的特殊形式，它们可以看作是频率为0、 ／(N一1)和 ／(N一1)的余弦曲线的线性组合，这里N为窗的长度，在MALAB中可以采用下面的命令来生成： 加窗后的窗口效应 矩形窗处理的谱估计的主瓣较窄，分辨率最好，同时其主瓣附近的衰减比其他两个窗函数小，旁瓣部分泄漏比较严重，所以其方差特性最差，噪声水平较高。 加海明窗后的谱估计效果 海明窗处理谱估计的旁瓣部分衰减较大，方差较小，噪声水平较低，性能良好，改善了由矩形窗处理的谱估计所产生的较大谱失真 加布莱克曼窗后的效果 布莱克曼窗处理谱估计的旁瓣部分衰减较大，方差较小，噪声水平较低，性能良好，改善了由矩形窗处理的谱估计所产生的较大谱失真 窗函数对估计结果的影响主要取决于其主瓣和旁瓣特性，所以窗函数评价指标有 。 窗函数评价指标 主瓣宽度 旁瓣大小 旁瓣衰减斜率 四、总结 例如测量震动物体的自震频率时，则可以选用主瓣宽度比较窄的Rectangular窗；当要求谱泄漏小时，则应选用旁瓣幅度较小的窗。 对于随时间衰减的信号如脉冲响应信号，可以采用指数窗来提高分析的信噪比。 对存在强干扰情况下的窄带信号，若干扰靠近信号，则可选用旁瓣幅度小的窗函数。 若离开通带较远，则可选用旁瓣衰减速度比较快的窗函数。 当要求分辨率高而不计及幅值的精度时，可选用主瓣窄而陡的窗。 谱分析时窗函数的选择 谢 谢 Bartlett 平均周期图的方法是将N 点的有限长序列x(n)分段求周期图再平均。平均周期图法使谱估计方差减少是用增加偏倚和降低谱分辨率为代价得来的。 xn=cos(2*pi*20*n)