Introduction to Feedback and Control Systems: Dynamic Systems Control
the typa guy to flex: "new plots"
Changed files
hw2/code_1-2.PNG
Binary files differ
hw2/code_1.PNG
Binary files differ
hw2/code_6.PNG
Binary files differ
hw2/plot_1-2.PNG
Binary files differ
hw2/plot_1.png
Binary files differ
hw2/plot_6-2.PNG
Binary files differ
hw2/plot_6-3.PNG
Binary files differ
hw2/plot_6.PNG
Binary files differ
hw2/prob_1.py
@@ -1,13 +1,31 @@
1
1
import numpy as np
2
2
import matplotlib.pyplot as plt
3
Added:
from scipy import signal
3
4
4
5
t = [0.001 * x for x in range(0, 141)]
5
6
y = [0.873
6
7
- 0.873 * np.exp(-46 * t) * np.cos(105.4 * t)
7
8
- 0.381 * np.exp(-46 * t) * np.sin(105.4 * t) for t in t]
8
9
10
Added:
# plt.plot(t, y)
11
Added:
# plt.xlabel('t (sec)')
12
Added:
# plt.ylabel('y(t)')
13
Added:
# plt.title('MAE 171A - Homework 2.1\n Output y(t)')
14
Added:
# plt.show()
15
Added:
16
Added:
# numH = []
17
Added:
# denH =
18
Added:
# sysH = signal.TransferFunction(numH, denH)
19
Added:
# impulse = signal.impulse(sysH)
20
Added:
21
Added:
system = ([11550], [1, 92, 13225, 0])
22
Added:
t, y = signal.impulse2(system)
9
23
plt.plot(t, y)
24
Added:
25
Added:
# y_s = [11550 / (t * (t**2 + 92 * t + 13225)) for t in t]
26
Added:
27
Added:
# plt.plot(impulse)
10
28
plt.xlabel('t (sec)')
11
Removed:
plt.ylabel('y(t)')
12
Removed:
plt.title('MAE 171A - Homework 2.1\n Output y(t)')
29
Added:
plt.ylabel('Amplitude')
30
Added:
plt.title('MAE 171A - Homework 2.1\n Impulse Response')
13
31
plt.show()
hw2/prob_5.py
@@ -0,0 +1,5 @@
1
Added:
import numpy as np
2
Added:
3
Added:
p = [1, 2, 3, 8, 8]
4
Added:
r = np.roots(p)
5
Added:
print(np.around(r, 2))
hw2/prob_6.m.txt
@@ -0,0 +1,4 @@
1
Added:
N=1;K=0;
2
Added:
D=[1 5 10 10 5 K];
3
Added:
sys=tf(N,D);
4
Added:
pzplot(sys)