Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- from numpy import array, float64, complex64
- from numpy.fft import fft as fft_numpy, ifft as ifft_numpy
- from matplotlib import pylab as plt
- #Изначальный массив отсчетов
- arr_in = array([140, 32, 1, 50], dtype = complex64)
- #Прямое БПФ
- forward_fft_numpy_res = fft_numpy(arr_in)
- #Обратное преобразование Фурье
- backward_fft_numpy_res = ifft_numpy(forward_fft_numpy_res)
- #Вывод данных и графиков
- print (arr_in)
- print (forward_fft_numpy_res)
- print (backward_fft_numpy_res)
- plt.xlabel('j')
- plt.ylabel('u')
- plt.bar ([0,1,2,3], arr_in,width = 0.02)
- plt.show()
- plt.xlabel('k')
- plt.ylabel('U')
- plt.bar ([0,1,2,3], forward_fft_numpy_res,width = 0.02,color = 'red')
- plt.show()
- plt.xlabel('j')
- plt.ylabel('u back')
- plt.bar ([0,1,2,3], backward_fft_numpy_res,width = 0.02)
- plt.show()
Add Comment
Please, Sign In to add comment