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