matlab - 如何将 fir1 阻带滤波器变成多阻带滤波器

Question

单个带阻滤波器可以通过fir1(N,[w1 w2],'stop')来制作，其中N是滤波器系数的个数，w1是开始停止频率，w2是结束停止频率，怎么转这变成一个多阻带滤波器？有没有直接的方法来做到这一点？

Answer 1

使用 conv 确定适当的系数：

b =  fir1(50,[0.2 0.3],'stop');
b2 =  fir1(50,[0.7 0.8],'stop');
freqz(conv(b,b2),1);

Answer 2

使用 'DC-0' 而不是 'stop' 并给它每个波段的开始和结束：

bd1start=1500/(fs/2);
bd1end=3500/(fs/2);

bd2start=6500/(fs/2);
bd2end=8500/(fs/2);
order=64;
bb=fir1( order, [bd1start bd1end bd2start bd2end ],'DC-1' );