详细请查阅《Why use fftshift(fft(fftshift(x))) instead of fft(x) in Matlab?》
使用fftshift(fft(fftshift(x)))后的效果:
1.不改变频谱的幅度和相位
2.使得频谱的范围为-Fs/2到Fs/2,中心频率为0
准确的频率是:
int n=N%2==0?N/2:((N-1)/2);
最左边频率 -Fs/N*n
最右边频率 Fs/N*(N-n-1)
其中N是系列长度,Fs是采样频率
3.频谱的实部全为正数,震荡性消失
使用注意:
为使得输出结果与理论分析相同,需要乘以系数
FFT: fftshift(fft(fftshift(x))) * dt
IFFT: fftshift(ifft(fftshift(X))) * fs
其次:
MATLAB中实现的IFFT已经除以序列长度L,FFTW中没有除以序列长度,需要自己补上