filter1.m

来自「一些数字信号处理课程的小程序」· M 代码 · 共 4 行

[b,a]=ellip(4,0.1,40,[10,20]*2/fs);
[H,w]=freqz(b,a,512);
plot(w*fs/(2*pi),abs(H));
grid;