6_3经典功率谱估计.ppt

周期图法的改进 2. 平滑周期图 (Blackman-Tukey 法 ) B-T 法进行功率谱估计的主要步骤： (1) 利用观测数据估计自相关序列。 (2) 对自相关函数估计值加窗。 (3) 计算加窗后自相关函数的 DTFT 。 优点： P M ( ? ) 波动比 I N ( ? ) 小，可证是一致估计 缺点：降低了频率分辨率 经典功率谱估计 周期图法的改进 3. 平均周期图法 (Welch-Bartlett 法 ) 将随机序列 X [ k ] 的 N 个观测值分成 A 段 x [ k ] ? x [ iM ? k ]; i ? 0 , 1 , ? A ? 1 ; k ? 0 , 1 , ? M ? 1 i 第 i 段序列的周期图为 平均周期图 1 j ? 2 I ( ? ) ? X M ( e ) M i M 1 i P ( ? ) ? ? I M ( ? ) A i ? 0 A M 经典功率谱估计 A ? 1 周期图法的改进 3. 平均周期图法 (Welch-Bartlett 法 ) 平均周期图法估计质量 1 i var{ P ( ? )} ? var{ I M ( ? )} A A M A ?? , 方差为零，是一致估计 因为 所以 n ? j ? n bia [ I N ( ? )] ? ? ? R [ n ] e n ? ? ( N ? 1 ) N N ? 1 bia [ I ( ? )] ? bia [ I N ( ? )] 经典功率谱估计 i M 平均周期图方差减小的代价之一是偏差增大 周期图法的改进 4. 重叠平均周期图法 (Welch 法 ) 平均周期图法 优点 ：减小方差 缺点 ：增加估计的偏差， 降低了谱的分辨率 原因 ：分段即加窗，段越多，窗越短，主瓣宽度越大 解决方法 ：将各段数据有一定程度的重叠 i = 1 M 0 M i = 2 经典功率谱估计 i = 3 M i = A M N - 1 例：利用 Welch 法进行平稳高斯白噪声的谱估计 产生 30 组 512 点均值为零，方差为 1 的平稳高斯 白噪声，利用 Welch 法按照 50% 重叠分别将其分 成 A =3 ， 7 ， 15 ， 31 段，计算功率谱估计值，并 分析谱估计质量。 分析： (1) 对每组 512 点数据按各段数据重叠 50% 的方式 分成 3 段 256 点序列， 7 段 128 点序列， 15 段 64 点序列， 31 段 32 点序列。 (2) 求出每段数据的周期图 A ? 1 i I M [ m ] ? 2 1 i X M [ m ] M (3) 再取平均即得各组数据的功率谱估计 , 即 1 i P ( ? ) ? ? I M ( ? ) A i ? 0 A M 数字信号处理 ( Digital Signal Processing ) 信号与系统系列课程组 国家电工电子教学基地 离散随机序列的特征描述 平稳随机序列通过 LTI 系统 经典功率谱估计 现代功率谱估计 经典功率谱估计 谱估计的质量 相关法 ( 间接法 ) 周期图法 ( 直接法 ) 周期图法的改进 利用 MATLAB 实现功率谱估计 经典功率谱估计 谱估计的质量 通常以估计量的偏差、方差、均方差等来衡量 估计的质量 1 ．估计量的偏差 bia { ? ? } ? E { ? ? } ? ? 表示近似程度。 当 bia { ? ? } ? 0 时，称无偏估计； N 趋向无穷时 bia { ? ? } ? 0 ，称渐进无偏估计。 经典功率谱估计 谱估计的质量 2 ．估计量的方差 2 ? ? ? var{ ? } ? E { ( ? ? E { ? } ) } 越小，估计值离估计均值越小 ? } ? 0 , 若 bia { ? N ? ? ? } ? 0 , 则称 ? ? 为 ? 的一致估计 lim var{ ? 3 ．估计量的均方差 2 ? ? ? ? MSE { ? } ? E { ( ? ? ? ) } ? var{ ? } ? bia { ? } ? } ? 0 , 则称 ? ? 为 ? 的一致估计 若 MSE { ? 经典功率谱估计 相关法 ( 间接法 ) 进行功率谱估计 相关法的理论基础 自相关函数估计的计算