Page 177 - 6734

y = A1*np.sin(omega1*x) + A2*np.sin(omega2*x) #
сигнал (лінійна комбінація синусоїдальних сигналів)
plt.plot(x,y,'k') # нарисувати сигнал
plt.xlabel('t');plt.ylabel('y');plt.grid();plt.show()
print "Рисунок - Сигнал"

yf = fft(y) # розрахувати спектр частот для сигналів
змінних у часі за допомогою дискретного перетворення
Фур'є. Декомпозиція сигналу на частоти і амплітуди.
xf = np.linspace(0.0, 1.0/(2.0*T), N/2) # частоти
yf_=2.0/N * np.abs(yf[0:N/2]) # амплітуди
plt.figure()
plt.plot(xf, yf_,'k') # нарисувати спектр
plt.xlabel('f');plt.ylabel('A');plt.grid();plt.show()
print "Рисунок - Амплітудно-частотна характеристика
сигналу"

Рисунок 38 – Сигнал

176