View raw

1 import numpy as np 2 import matplotlib.pyplot as plt 3 from scipy import signal 4 5 t = [0.001 * x for x in range(0, 141)] 6 y = [0.873 7 - 0.873 * np.exp(-46 * t) * np.cos(105.4 * t) 8 - 0.381 * np.exp(-46 * t) * np.sin(105.4 * t) for t in t] 9 10 # plt.plot(t, y) 11 # plt.xlabel('t (sec)') 12 # plt.ylabel('y(t)') 13 # plt.title('MAE 171A - Homework 2.1\n Output y(t)') 14 # plt.show() 15 16 # numH = [] 17 # denH = 18 # sysH = signal.TransferFunction(numH, denH) 19 # impulse = signal.impulse(sysH) 20 21 system = ([11550], [1, 92, 13225, 0]) 22 t, y = signal.impulse2(system) 23 plt.plot(t, y) 24 25 # y_s = [11550 / (t * (t**2 + 92 * t + 13225)) for t in t] 26 27 # plt.plot(impulse) 28 plt.xlabel('t (sec)') 29 plt.ylabel('Amplitude') 30 plt.title('MAE 171A - Homework 2.1\n Impulse Response') 31 plt.show() 32