首先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滤波之后的图像