[P,F]=periodogram(x, window, nfft, Fs)

function [P,f]=fftpsd(x,Fs)
    m = length(x);          % 整个数据长度
    n = pow2(nextpow2(m));  % 扩展到合适长度，2的整数次幂
    y = fft(x,n);           
    f = (0:n-1)*(Fs/n);     %频率范围
    P = y.*conj(y)/n;   %功率谱密度

    f=f(1:n/2);
    P=P(1:n/2);
   
    %plot(f,10*log10(P)); %把数据转换为dB
    plot(f,P);
    grid;
    xlabel('Frequency (Hz)');
    ylabel('Power(dB)');
    title('{\bf Periodogram}');

ty = y(1:5:end); %将原来2000Hz的数据降采样到400Hz
ty = ty - mean(ty); %减去直流分量
[P,F]=periodogram(ty,[],[],400); %求功率谱
plot(F, P) %画图

696000楼: Originally posted by sudo at 2011-07-25 20:28:50
唉，可能解释起来一言难尽~建议楼主阅读一下信号处理的入门书嗯

首先，你师兄的图，应该是降采样过后的（或者是后来截取了一段，不过我倾向于认为是降采样），也就是说他用来做功率谱的数据采样率只有400Hz，你 ...

3楼: Originally posted by sudo at 2011-07-25 20:44:28
对了，还有一些想补充的~

首先就是直流分量的问题，如果不先对数据都减去一个均值，那么接近0Hz部分的功率会非常非常大~这是有害的，一般做数据处理的时候如果只关注交流量的话，就要先使得整个序列的均值为0，方 ...

3楼: Originally posted by sudo at 2011-07-25 20:44:28
对了，还有一些想补充的~

首先就是直流分量的问题，如果不先对数据都减去一个均值，那么接近0Hz部分的功率会非常非常大~这是有害的，一般做数据处理的时候如果只关注交流量的话，就要先使得整个序列的均值为0，方 ...

9楼: Originally posted by mallory19 at 2014-04-03 23:53:13
sudo，你好。
我在用Mathematica中的Periodogram和PeriodogramArray对同一个数据列表计算时，得到的结果却不同。我在网上查了一下，也有别人在问这个问题。但没有回答。
不知道你对Mathematica是否了解，是否可以 ...