MATLAB中MVDR谱估计式,基于Welch算法的功率谱估计的实现

第 25 卷 第 3 期 2007 年 5 月 北京工商大学学报(自然科学版) Journal of Beijing Technology and Business University(Natural Science Edition) Vol25 No13 May 2007 　　文章编号 :167121513(2007) 030058202 基于 Welch 算法的功率谱估计的实现 罗 　敏1 , 　刘 　嵩2 (1. 湖北民族学院 理学院 , 湖北 恩施 445000 ; 2. 湖北民族学院 信息工程学院 , 湖北 恩施 445000) 摘 　要 : 在经典谱估计领域 ,由 Welch 提出的修正周期图法由于其简单 ,易于理解 ,便于计算等优点 ,获得了有效的应用 ,在多数情况下 ,其频率分辨率、估计方差能满足实际应用. 运用 MATLAB 软件进行 Welch 算法仿真 , 并重点分析了窗函数对谱估计的影响. 关键词 : 功率谱估计 ; 周期图 ; 窗函数 ; 分辨率 中图分类号: TN91117 　　　　　文献标识码 : A 收稿日期: 2007 01 03 基金项目: 国家民委自然科学基金项目(MZZ04004) 作者简介: 罗 　敏(1978 —) ,女 ,湖北随州人 ,硕士研究生 ,研究方向为模式识别与 MATLAB 软件应用研究. 　　功率谱估计应用范围很广 ,日益受到各学科和应用领域的极大重视. 它是用有限长的数据来估计信号的功率谱 ,是数字信号处理的重要研究内容之一[1]. 而 MATLAB 是一个强大的数值计算软件 ,程序设计自由度大 ,程序的可移植性好并具有强大的图形处理功能 ,因此可以利用 MATLAB 来实现功率谱估计. 介绍了用 MATLAB 实现基于 Welch 算法的功率谱估计. 1 　功率谱估计的 Welch 算法 经典谱估计方法在工程中都是以离散傅立叶变换为基础的 ,它隐含着对无限长数据加窗处理 ,所以经典谱估计有着分辩率不高、能量泄漏的固有缺点[2]. 为了克服这些缺点 ,人们做过长期的努力 ,提出了各种算法. 目前经典谱估计算法有周期图法、 Barlett 算法、Welch 算法、Nattall 算法等 ,Welch 算法是由 Welch 提出的对周期图的修正算法 ,是经典谱估计中获得有效应用的一种算法[3]. Welch 算法谱估计采取数据分段加窗处理再求平均的办法 ,先分别求出每段的谱估计 ,然后进行总平均[4]. 根据概率统计理论 ,若将原长度为 N 的数据分成 L 段 ,每段长度取 M = N/ L ,且各段数据互为独立 ,则估计的方差将只有原来不分段的 1/ L ,达到一致估计的目的. 其算法如式(1) . GL M (ω) = 1 L ∑ L l =1 GM , l (ω) . (1) 但若 L 增加 , M 减小 ,则分辨率下降. 相反 ,若L 减小 , M 增加 ,虽然估计的偏差减小 ,但估计的方差增大. 所以在实际中必须兼顾分辨率与方差的要求 ,适当选取 L 和 M 的值. 在分段时为了减小因分段数增加给分辨率带来的影响 ,采取各段数据有一定重叠的方法[5]. 2 　用 MATLAB实现 Welch 算法 211 　软件实现 Welch 法谱估计流程如图 1. 图 1 　Welch 法谱估计流程 在 MATLAB 中 ,psd 函数 与 pwelch 函数都可以实现 Welch 法的功率谱估计 ,其方法是一样的 ,只是部分参数设置有所不同. 利用 MATLAB 产生一随机信号 ,然后用