Wp = 0.2;
Ws = 0.4;
Rp = 0.5;
Rs = 60;

[n,Wn] = ellipord(Wp,Ws,Rp,Rs)
[b,a] = ellip(n,Rp,Rs,Wp)
figure;freqz(b,a);
figure;zplane(b,a);


% [n,Wn] = buttord(Wp,Ws,Rp,Rs)
% [b,a] = butter(n,Wn)
% figure;zplane(b,a);
% 
% [n,Wn] = cheb1ord(Wp,Ws,Rp,Rs)
% [b,a] = cheby1(n,Rp,Wp)
% figure;zplane(b,a);
% 
% [n,Wn] = cheb2ord(Wp,Ws,Rp,Rs)
% [b,a] = cheby2(n,Rp,Wp)
% figure;zplane(b,a);



Wp = 0.8;
Ws = 0.6;
Rp = 0.5;
Rs = 60;

[n,Wn] = ellipord(Wp,Ws,Rp,Rs)
[b,a] = ellip(n,Rp,Rs,Wp, 'high')
figure;freqz(b,a);
figure;zplane(b,a);


Wp = [0.2,0.4];
Ws = [0.1,0.5];
Rp = 0.5;
Rs = 60;

[n,Wn] = ellipord(Wp,Ws,Rp,Rs)
[b,a] = ellip(n,Rp,Rs,Wp)
figure;freqz(b,a);
figure;zplane(b,a);


Wp = [0.4,0.6];
Ws = [0.3,0.7];
Rp = 0.5;
Rs = 60;

[n,Wn] = ellipord(Wp,Ws,Rp,Rs)
[b,a] = ellip(n,Rp,Rs,Wp)
figure;freqz(b,a);
figure;zplane(b,a);


Ws = [0.4,0.6];
Wp = [0.3,0.7];
Rp = 0.5;
Rs = 60;

[n,Wn] = ellipord(Wp,Ws,Rp,Rs,'s')
[b,a] = ellip(n,Rp,Rs,Wp,'stop')
figure;freqz(b,a);
figure;zplane(b,a);