blob: 20585c887e8d9c98a8768118c8e9ef684e1dc6c7 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
|
import numpy as np
import matplotlib.pyplot as plt
from math import pi
N = 300
fr = np.logspace(-3, 3, N)
w = 2 * pi * fr
s = w * 1j
R1 = 1e6
R2 = 1e4
C = 1.6e-6 / 2 / pi
num = R1 * R2 * C * s + R2
den = R1 * R2 * C * s + R1 + R2
H = num / den
plt.figure(1)
plt.subplot(211)
plt.loglog(fr, abs(H))
plt.xlabel('Frequency (Hz)')
plt.ylabel('Magnitude (V/V)')
plt.title('R1 = 1MΩ, R2 = 1kΩ,C = 1.6μF, Rf = 100kΩ')
plt.subplot(212)
plt.semilogx(fr, np.angle(H) * 180 / pi)
plt.xlabel('Frequency (Hz)')
plt.ylabel('Phase (deg)')
plt.show()
|
