exa4_18.m

来自「新版Matlab 7辅助信号处理技术与应用」· M 代码 · 共 17 行

% exa4-18_ellipstop.m , for example 4-18
% using chebyshev II filter to design bandpass DF

clear all;
Wpl=650;
Wph=850;
Wp=[Wpl,Wph];
Rp=0.1;
Wsl=700;
Wsh=800;
Ws=[Wsl,Wsh];
Rs=50;
Fs=2000;
[N,Wn]=cheb2ord(Wp/(Fs/2),Ws/(Fs/2),Rp,Rs);
[b,a]=ellip(N,Rp,Rs,Wn,'stop');
freqz(b,a,512,1000)