通过Matlab 使用 FFT 分析周期性数据

傅里叶变换可以用来分析数据中的变化

前言：天文学家使用苏黎世太阳黑子相对数将几乎 300 年的太阳黑子的数量和大小制成表格。对大约 1700 至 2000 年间的苏黎世数绘图。

load sunspot.dat
year = sunspot(:,1);
relNums = sunspot(:,2);
plot(year,relNums)
xlabel('Year')
ylabel('Zurich Number')
title('Sunspot Data')

为了更详细地看太阳黑子活动的周期特性，将对前 50 年的数据绘图。

plot(year(1:50),relNums(1:50),'b.-');
xlabel('Year')
ylabel('Zurich Number')
title('Sunspot Data')

傅里叶变换是一种基础的信号处理工具，可确定数据中的频率分量。使用 fft 函数可获取苏黎世数据的傅里叶变换。删除存储数据总和的输出的第一个元素。绘制该输出的其余部分，其中包含复傅里叶系数关于实轴的镜像图像。

y = fft(relNums);
y(1) = [];
plot(y,'ro')
xlabel(&#