计算功率谱密度的两种方法,直接计算和periodogram法求解

点击联系发帖人 时间：2017-10-05 12:07

直接法又称周期图法它是把随機序列x(n)的N个观测数据视为一能量有限的序列，直接计算x(n)的离散傅立叶变换得X(k)，然后再取其幅值的平方并除以N，作为序列x(n)真实功率谱的估计

间接法先由序列x(n)估计出自相关函数R(n)，然后对R(n)进行傅立叶变换便得到x(n)的功率谱估计。

 

 3.改进的直接法：
对于直接法的功率谱估计当數据长度N太大时，谱曲线起伏加剧若N太小，谱的分辨率又不好因此需要改进。

 

 Bartlett平均周期图的方法是将N点的有限长序列x(n)分段求周期图再岼均
Matlab代码示例：

 

 Welch法对Bartlett法进行了两方面的修正，一是选择适当的窗函数w(n)并再周期图计算前直接加进去，加窗的优点是无论什么样的窗函數均可使谱估计非负二是在分段时，可使各段之间有重叠这样会使方差减小。
Matlab代码示例：

}

f=x; %在信号中加入白噪声 y=fft(f,1024); %对原始信号進行离散傅里叶变换参加DFT采样点的个数为1024 %用周期图法估计功率谱密度 ff=1)/1024; %计算变换后不同点对应的频率值 title('功率谱密度（周期图法）');

求FFT变换时，纵坐标和横坐标的问题

}