plt.subplot(2, 1, 1)
plt.plot(np.arange(0, sampling_period, interval), ori_func, 'black')
plt.xlabel...interval)
nfft = abs(ft[range(int(n / 2))] / n )
plt.plot(frequency, nfft, 'red')
plt.xlabel...fft_shift = np.fft.fftshift(fft)
plt.imshow(np.log(np.abs(fft_shift)))
# matplotlib.image.AxesImage...((2 * m - 1, int(n / 2)))
fft_shift = np.concatenate((b, fft_shift, b), axis = 0)
fft_shift = np.concatenate...# (633L, 1321L)
ifft = np.real(ifft)
plt.imshow(ifft, cmap = plt.get_cmap('gray'))
# matplotlib.image.AxesImage