功率谱估计方法介绍【转】(2009-01-20 11:12:25)标签:谱估计 功率谱 周期图 matlab 傅里叶 文化 分类:学习
功率谱估计在现代信号处理中是一个很重要的课题,涉及的问题很多。在这里,结合matlab,我做一个粗略介绍。功率谱估计可以分为经典谱估计方法与现代谱估计方法。经典谱估计中最简单的就是周期图法,又分为直接法与间接法。直接法先取N点数据的傅里叶变换(即频谱),然后取频谱与其共轭的乘积,就得到功率谱的估计;间接法先计算N点样本数据的自相关函数,然后取自相关函数的傅里叶变换,即得到功率谱的估计.都可以编程实现,很简单。在matlab中,周期图法可以用函数periodogram实现。
但是周期图法估计出的功率谱不够精细,分辨率比较低。因此需要对周期图法进行修正,可以将信号序列x(n)分为n个不相重叠的小段,分别用周期图法进行谱估计,然后将这n段数据估计的结果的平均值作为整段数据功率谱估计的结果。还可以将信号序列x(n)重叠分段,分别计算功率谱,再计算平均值作为整段数据的功率谱估计。
这2种称为分段平均周期图法,一般后者比前者效果好。加窗平均周期图法是对分段平均周期图法的改进,即在数据分段后,对每段数据加一个非矩形窗进行预处理,然后在按分段平均周期图法估计功率谱。相对于分段平均周期图法,加窗平均周期图法可以减小频率泄漏,增加频峰的宽度。welch法就是利用改进的平均周期图法估计估计随机信号的功率谱,它采用信号分段重叠,加窗,FFT等技术来计算功率谱。与周期图法比较,welch法可以改善估计谱曲线的光滑性,大大提高谱估计的分辨率。matlab中,welch法用函数psd实现。调用格式如下:
[Pxx,F] = PSD(X,NFFT,Fs,WINDOW,NOVERLAP)
X:输入样本数据
NFFT:FFT点数
Fs:采样率
WINDOW:窗类型
NOVERLAP,重叠长度
现代谱估计主要针对经典谱估计分辨率低和方差性不好提出的,可以极大的提高估计的分辨率和平滑性。可以分为参数模型谱估计和非参数模型谱估计。参数模型谱估计有AR模型,MA模型,ARMA模型等;非参数模型谱估计有最小方差法和MUSIC法等。由于涉及的问题太多,这里不再详述,可以参考有关资料。
matlab中,现代谱估计的很多方法都可以实现。music方法用pmusic命令实现;pburg函数利用burg法实现功率谱估计;pyulear函数利用yule-walker算法实现功率谱估计等等。
另外,sptool工具箱也具有功率谱估计的功能。窗口化的操作界面很方便,而且有多种方法可以选择。
分享到:
相关推荐
关于功率谱的详细说明,各种功率谱估计方法的描述
参数功率谱估计方法参数功率谱估计方法参数功率谱估计方法
经典功率谱估计与现代功率谱估计的对比,有程序,并加详细的分析,还有语谱图的分析
经典功率谱估计 BT法(间接法)
Welch法的功率谱估计实现,工具matlab。
经典功率谱估计 Bartlett法 自己编写的程序
利用小波变换的原理进行功率谱估计,包括经典功率谱估计和现代功率谱估计的主要方法MATLAB仿真,小波变换的分解重构和小波包变换的内容。
MUSIC算法估计功率谱
Welch法的功率谱估计实现,工具matlab。
信号为两个正弦信号加高斯白噪声,各正弦信号的信噪比均为10dB,长度为N,信号频率分别为f1,f2;初始相位均为0,取f1/fs=0.2,f2/fs取不同数值:0.3、0.25,,fs...利用经典功率谱估计和现代功率估计方法进行功率谱估计
功率谱估计及其MATLAB仿真 功率谱估计及其MATLAB仿真
其中包括经典功率谱估计和现代功率谱估计方法 自相关法,周期图法,batlet法,Welch法,Yule-walk法和Burg法 并附有matlab程序
介绍了现代功率谱估计中AR模型参数的几种典型求解算法。并借助matlab平台对各种算法的功率谱进行了仿真
用ls算法实现功率谱估计,现代信号处理。可直接运行。
本文主要介绍对比了两种经典谱估计方法:周期图法与自相关法和三种参数谱估计法:直接求解yule-walker方程,Levinson-Durbin快速递推法与Borg算法的原理、编程过程、优劣分析,并比较了三种参数谱方法的估计性能,还...
经典功率谱估计 Welch法 自己编写的程序,可以直接使用
经典功率谱估计,周期图法功率谱估计(直接法),没用MATLAB自带的函数,自己编写的函数,运行结果和MATLAB自带函数一样
BT法和周期图做功率谱估计 算法流程参考 何子述的现在数字信号处理及应用
capon用于功率谱估计的程序。现代谱估计的方法之一,程序可以运行,很好用。
从介绍功率谱的估计原理入手分析了经典谱估计和现代谱估计两类估计方法的原理、各自特点及在Matlab 中的实现 方法。