Matlab中的FFT函数示例。

1、在命令行窗口，输入如下命令：load sunspot.datyear = sunspot(:,1);relNums = sunspot(:,2);plot(year,relNums)title('Sunspot Data')如图1所示。

3、以下是前50年的近况。在命令行窗口，输入如下命令：plot(year(1:50),relNums(1:50),'b.-');如图3所示。

5、Y的复震级平方称为功率，功率与频率的关系图称为“周期图”。在命令行窗口，输入如下命令：n = length(Y);power = abs(Y(1:floor(n/2))).^2;nyquist = 1/2;freq = (1:n/2)/(n/2)*nyquist;plot(freq,power)xlabel('cycles/year')title('Periodogram')如图5所示。

7、为了方便起见，我们绘制了功率与周期的关系图（其中period=1./freq）。正如预期的那样，有一个非常显著的周期，其长度约为11年。在命令行窗口，输入如下命令：period = 1./freq;plot(period,power);axis([0 40 0 2e+7]);ylabel('Power');xlabel('Period (Years/Cycle)');如图7所示。