line.py 文件源码

def swap_Nq(fft_y, fu, fft_x, N):
    f = fu[:, 0].copy()
    fft_x[0] = f[0].real
    fft_x[1:N//2] = 0.5*(f[1:N//2] + np.conj(f[:N//2:-1]))
    fft_x[N//2] = f[N//2].real
    fu[:N//2+1, 0] = fft_x[:N//2+1]
    fu[N//2+1:, 0] = np.conj(fft_x[(N//2-1):0:-1])

    fft_y[0] = f[0].imag
    fft_y[1:N//2] = -0.5*1j*(f[1:N//2] - np.conj(f[:N//2:-1]))
    fft_y[N//2] = f[N//2].imag

    fft_y[N//2+1:] = np.conj(fft_y[(N//2-1):0:-1])
    return fft_y