filter1.m

一些数字信号处理课程的小程序

[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;