def invFourier(X,fs,N): # Approximate inverse Fourier transform on the interval (-T/2,T/2)
NFFT=len(X)
k=np.arange(NFFT)
x=fs*np.exp(1j*np.pi*(N/2.-k))*ifft(X*np.exp(-1j*np.pi*k*N/NFFT),NFFT)
x=x[:N]
return x
评论列表
文章目录