matlab 频谱不对称_Matlab中fft与fwelch有什么区别？如何用fft求功率谱？

讲这个话题，就要先搞清楚频谱、功率谱的概念，可参考我的另一篇文章

信号的频谱 频谱密度 功率谱密度 能量谱密度

做信号处理的朋友应该都会fft比较熟悉，就是求傅里叶变换。我在这里也不再去讲这个函数了，但需要注意的一点：实信号的频谱关于0频对称，是偶函数，如果st = cos(2pif0*t)+1;  t的长度为4000，那么0频的位置在第一个点，做fftshift后，0频的位置在低2001个点的位置，fft后的信号关于第2001个点对称，而不是4000个点左右对称。

pwelch是用来求功率谱的，采用Welch平均周期法对信号进行谱估计，它通过分段选取数据进行加窗求功率，再进行平均，pwelch函数的使用方式为：

pxx = pwelch(x,window,noverlap,nfft)
[pxx,f] = pwelch(x,window,noverlap,f,fs)

其中，

• X表示输入序列；

• window：当window是一个数值时，表示窗函数长度，即分段长度L，默认的窗函数为hamming窗；当window是一个序列时，表示窗函数序列；

• NFFT表示FFT的点数