Introduction to Feedback and Control Systems: Dynamic Systems Control
homework submission
Changed files
hw1/code_2.PNG
Binary files differ
hw1/code_5.PNG
Binary files differ
hw1/plot_2-a.png
Binary files differ
hw1/plot_5.png
Binary files differ
hw1/prob_2.py
@@ -0,0 +1,30 @@
1
Added:
import numpy as np
2
Added:
import matplotlib.pyplot as plt
3
Added:
from math import pi
4
Added:
5
Added:
N = 300
6
Added:
fr = np.logspace(-3, 3, N)
7
Added:
w = 2 * pi * fr
8
Added:
s = w * 1j
9
Added:
10
Added:
R1 = 1e6
11
Added:
R2 = 1e4
12
Added:
C = 1.6e-6 / 2 / pi
13
Added:
14
Added:
num = R1 * R2 * C * s + R2
15
Added:
den = R1 * R2 * C * s + R1 + R2
16
Added:
H = num / den
17
Added:
18
Added:
plt.figure(1)
19
Added:
plt.subplot(211)
20
Added:
plt.loglog(fr, abs(H))
21
Added:
plt.xlabel('Frequency (Hz)')
22
Added:
plt.ylabel('Magnitude (V/V)')
23
Added:
plt.title('R1 = 1MΩ, R2 = 1kΩ,C = 1.6μF, Rf = 100kΩ')
24
Added:
25
Added:
plt.subplot(212)
26
Added:
plt.semilogx(fr, np.angle(H) * 180 / pi)
27
Added:
plt.xlabel('Frequency (Hz)')
28
Added:
plt.ylabel('Phase (deg)')
29
Added:
30
Added:
plt.show()
hw1/prob_5.py
@@ -0,0 +1,11 @@
1
Added:
import numpy as np
2
Added:
import matplotlib.pyplot as plt
3
Added:
4
Added:
t = [0.01 * x for x in range(0, 1001)]
5
Added:
y = [-(t + 1) * np.exp(-t) for t in t]
6
Added:
7
Added:
plt.plot(t, y)
8
Added:
plt.xlabel('t (sec)')
9
Added:
plt.ylabel('y(t)')
10
Added:
plt.title('MAE 171A - Homework 1.5\n Step Response')
11
Added:
plt.show()