一、结构陷波滤波器
1.代码
kesi1_f=0.0002;%开口大小
kesi=40;%深度
num1=2*pi*400;%陷波频率
den1=num1;
G_filter_y1=tf([1/num1/num1 2*kesi1_f/num1 1],[1/den1/den1 2*kesi*kesi1_f/den1 1]);
bode(G_filter_y1)
G_filter_y1z=c2d(G_filter_y1,0.0002,'prewarp',num1)
Epf11=G_filter_y1z.num{1}
Epf12=G_filter_y1z.den{1}
2. 输出结果
3. 滤波器bode图
二、陷波器
1.代码
w1=2,wj=2;%陷波器中心频率w1=wj
ksai1=100;
ksai2=500;%ksai1/ksai2越小,深度越大;ksai1越大谷底越平坦
Gf=tf([1/w1^2 2*ksai1/w1 1],[1/wj^2 2*ksai2/wj 1]);
bode(Gf)
2. 输出
3. bode图
三、切比雪夫二型滤波器(cheby2)
1. 代码
[bb3 aa3]=cheby2(4,50,0.04);%阶数、深度、起始频率=(采样频率/2)*X
GG=tf(bb3,aa3,0.0002)
bode(GG)
EG11=GG.num{1}
EG12=GG.den{1}
2. 输出
3. bode图
四、filterdesigner设计滤波器
参考:https://blog.csdn.net/weixin_43778388/article/details/123952909