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
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
|
# %%
import matplotlib.pyplot as plt
import scipy.signal
import numpy as np
import matplotlib
from matplotlib.patches import Rectangle
def ellip_filter(N):
order = N
passband_ripple_db = 3
stopband_attenuation_db = 20
omega_c = 1000
a, b = scipy.signal.ellip(
order,
passband_ripple_db,
stopband_attenuation_db,
omega_c,
btype='low',
analog=True,
output='ba',
fs=None
)
w, mag_db, phase = scipy.signal.bode((a, b), w=np.linspace(0*omega_c,2*omega_c, 4000))
mag = 10**(mag_db/20)
passband_ripple = 10**(-passband_ripple_db/20)
epsilon2 = (1/passband_ripple)**2 - 1
FN2 = ((1/mag**2) - 1)
return w/omega_c, FN2 / epsilon2
plt.figure(figsize=(4,2.5))
for N in [5]:
w, FN2 = ellip_filter(N)
plt.semilogy(w, FN2, label=f"$N={N}$")
plt.gca().add_patch(Rectangle(
(0, 0),
1, 1,
fc ='green',
alpha=0.2,
lw = 10,
))
plt.gca().add_patch(Rectangle(
(1, 1),
0.01, 1e2-1,
fc ='green',
alpha=0.2,
lw = 10,
))
plt.gca().add_patch(Rectangle(
(1.01, 100),
1, 1e6,
fc ='green',
alpha=0.2,
lw = 10,
))
plt.xlim([0,2])
plt.ylim([1e-4,1e6])
plt.grid()
plt.xlabel("$w$")
plt.ylabel("$F^2_N(w)$")
plt.legend()
plt.savefig("F_N_elliptic.pdf")
plt.show()
|