-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathBW_IIR_filter.m
64 lines (60 loc) · 964 Bytes
/
BW_IIR_filter.m
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
%IIR BUTTERWORTH FILTER
syms s
kp=-2;
ks=-10;
wp=20;
ws=30;
j=sqrt(-1);
N=log((10^(-kp/10)-1)/(10^(-ks/10)-1))/(2*log(wp/ws))
N=round(N)+1
wc=wp/(10^(-kp/10)-1)^(1/(2*N))
H1=1/(s+1);
H2=1/(s^2+1.414*s+1);
H3=1/((s^2+s+1)*(s+1));
H4=1/((s^2+0.765368*s+1)*(s^2+1.84776*s+1));
H5=1/((s+1)*(s^2+0.6180*s+1)*(s^2+1.6180*s+1));
if N==1
H(s)=H1;
elseif N==2
H(s)=H2;
elseif N==3
H(s)=H3;
elseif N==4
H(s)=H4;
elseif N==5
H(s)=H5;
end
H(s)=H(s/wc)
H(s)=H(j*wp)
n=abs(H(s))
m=20*log(H(s))
%IIR BUTTERWORTH FILTER
syms s
kp=-2;
ks=-10;
wp=20;
ws=30;
j=sqrt(-1);
N=log((10^(-kp/10)-1)/(10^(-ks/10)-1))/(2*log(wp/ws))
N=round(N)+1
wc=wp/(10^(-kp/10)-1)^(1/(2*N))
H1=1/(s+1);
H2=1/(s^2+1.414*s+1);
H3=1/((s^2+s+1)*(s+1));
H4=1/((s^2+0.765368*s+1)*(s^2+1.84776*s+1));
H5=1/((s+1)*(s^2+0.6180*s+1)*(s^2+1.6180*s+1));
if N==1
H(s)=H1;
elseif N==2
H(s)=H2;
elseif N==3
H(s)=H3;
elseif N==4
H(s)=H4;
elseif N==5
H(s)=H5;
end
H(s)=H(s/wc)
H(s)=H(j*wp)
n=abs(H(s))
m=20*log(H(s))