Лаба 4 / main

import matplotlib.pyplot as plt
import numpy as np
import scipy as sp

freq = np.random.uniform(low=10, high=20, size=(1, ))
fs = 150

ampl = np.random.uniform(low=0.5, high=1, size=(1, ))

addition = np.random.uniform(low=-0.001,high=0.001,size=(1000,))

sig = ampl*np.sin(2 * np.pi * np.arange(1000) * freq / fs) + addition


#plt.plot(sig)
#plt.show()


fft_sig = sp.fftpack.fft(sig)

# plt.plot(abs(fft_sig))

# plt.show()

# plt.yscale('log')

# plt.plot(abs(fft_sig))

# plt.show()


#plt.plot(sig_window)
#plt.show()

hamming_window = np.hamming(1000)
sig_window = fft_sig*hamming_window
plt.yscale('log')
plt.plot(abs(sig_window))
plt.show()