the typa guy to flex: "new plots"

1 files changed, 20 insertions, 2 deletions

diff --git a/hw2/prob_1.py b/hw2/prob_1.py
index 615d855..e26a5fc 100644
--- a/hw2/prob_1.py
+++ b/hw2/prob_1.py
@@ -1,13 +1,31 @@
 import numpy as np
 import matplotlib.pyplot as plt
+from scipy import signal
 
 t = [0.001 * x for x in range(0, 141)]
 y = [0.873
      - 0.873 * np.exp(-46 * t) * np.cos(105.4 * t)
      - 0.381 * np.exp(-46 * t) * np.sin(105.4 * t) for t in t]
 
+# plt.plot(t, y)
+# plt.xlabel('t (sec)')
+# plt.ylabel('y(t)')
+# plt.title('MAE 171A - Homework 2.1\n Output y(t)')
+# plt.show()
+
+# numH = []
+# denH =
+# sysH = signal.TransferFunction(numH, denH)
+# impulse = signal.impulse(sysH)
+
+system = ([11550], [1, 92, 13225, 0])
+t, y = signal.impulse2(system)
 plt.plot(t, y)
+
+# y_s = [11550 / (t * (t**2 + 92 * t + 13225)) for t in t]
+
+# plt.plot(impulse)
 plt.xlabel('t (sec)')
-plt.ylabel('y(t)')
-plt.title('MAE 171A - Homework 2.1\n Output y(t)')
+plt.ylabel('Amplitude')
+plt.title('MAE 171A - Homework 2.1\n Impulse Response')
 plt.show()