我有一组数据,以3142x2双矩阵的形式给出,我想计算DFT(X[k] from k=0 to k=N-1
)。第1列表示时间,第2列表示太阳黑子数。
计算DFT后,我想绘制幅度函数abs(X[k])
。
以下是我尝试解决的问题,其中 ssd
是数据集:
y = fft(ssd(:,2)');
plot(ssd(:,1)',abs(y))
我很确定这不对,但我不明白错误。
如果文件中有任何用途,请输入带有数据的文件:https://dl.dropboxusercontent.com/u/15991564/sunspots.m