FMCW_phase_error_correction.py 文件源码

def Fourier(x,fs,NFFT):     # Approximate Fourier transform on the interval (-fs/2,fs/2)
    k=np.arange(NFFT)
    N=len(x)
    T=N/fs
    f=(-.5+k/float(NFFT))*fs
    n=np.arange(N)
    return f, T/N*np.exp(1j*pi*f*T)*fft(x*(-1.)**n,NFFT)