求功率谱周期

flyapple1992

还是matlab强大

Obtain the periodogram of the Wolf (relative sunspot) number data sampled yearly between 1700 and 1987.
Load the relative sunspot number data. Obtain the periodogram using the default rectangular window and number of DFT points (512 in this example).
The sampling rate for these data is 1 sample/year.
Plot the periodogram.
load sunspot.dat
relNums=sunspot(:,2);
[pxx,f] = periodogram(relNums,[],[],1);
plot(f,10*log10(pxx))
xlabel('Cycles/Year'); ylabel('dB');
title('Periodogram of Relative Sunspot Number Data');
You see in the preceding figure that there is a peak in the
periodogram at approximately 0.1 cycles/year, which indicates a period
of approximately 10 years.

好方便！

大力水手

有图吗？？

kingjames

没有具体实际例子  看起来这个matlab还是很简单的  mark  感谢楼主的分享  学习一番 和我自己的貌似差异还是很大的