From feaa9144eb28dd263ef79577ab369174b8c6c83b Mon Sep 17 00:00:00 2001 From: Marius Peter Date: Thu, 11 Jul 2019 17:14:32 -0700 Subject: the typa guy to flex: "new plots" --- hw2/prob_1.py | 22 ++++++++++++++++++++-- 1 file changed, 20 insertions(+), 2 deletions(-) (limited to 'hw2/prob_1.py') 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() -- cgit v1.2.3