在C++中使用MATLAB的FFT

首先FFT后的数据是频域方面的数据，所以X轴一定是固定的（或者手动变动）不会是一直变动的，那么就需要和Hz对应上。
所以这里说一下，关于得出X轴数据的方法。

绘制出随频率变化的振幅

fftdata = fft(data);
N=length(fftdata);
n=0:N-1;
f=n*fs/N;
FFTy = abs(fftdata(1:N));
FFTx = f(1:N); // 1:N这个和上面的一定要对应上负责会有问题
figure;
plot(FFTx ,FFTy );

绘制Nyquist频率之前随频率变化的振幅

fftdata = fft(data);
N=length(fftdata);
n=0:N-1;
f=n*fs/N;
FFTy = abs(fftdata(1:N/2));
FFTx = f(1:N/2); // 1:N这个和上面的一定要对应上负责会有问题
figure;
plot(FFTx ,FFTy );

两个差别不是很大只要显示部分中没有重叠区域就可以，两种计算出来的结果，分别是

我这些图像都是1-40Hz滤波之后的图像