为什么使用 fftshift(fft(fftshift(x))) 而不是 fft(x)

详细请查阅《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中没有除以序列长度，需要自己补上