View raw

1 import numpy as np 2 import matplotlib.pyplot as plt 3 from scipy import signal 4 5 # numH = [] 6 # denH = 7 # sysH = signal.TransferFunction(numH, denH) 8 # impulse = signal.impulse(sysH) 9 10 11 num = 1 12 den = [1, 2, 2] 13 sys = signal.TransferFunction(num, den) 14 t = [0.0001 * t for t in range(60000)] 15 y = signal.step(sys, t) 16 17 18 plt.plot(t, y) 19 plt.xlabel('t (sec)') 20 plt.ylabel('y(t)') 21 plt.title('MAE 171A - Homework 3.3\n Step Response') 22 plt.show() 23