bandpass.asv

来自「音频数据采集与IIR处理例程」· ASV 代码 · 共 12 行

close all;
Fp1=150;Fp2;3400;Fs1=50;Fs2=4000;
Wp1=Fp1/24000;Wp2=Fp1/24000;
Ws1=Fs1/24000;Ws2=Fs2/24000;
Wp=[Wp1 Wp2];
Ws=[Ws1 Ws2];
Rs=1;Rs=40;
[N,Wn]=ellipord(Wp,Ws,Rp,Rs);
[B,A]=ellip(N,1,40,Wn);
[Z,P,K]=tf2zp(B,A);
coeff=zp2sos(Z,P,K);