clear all;
close all;
clc;

Rp=2; Rs=3; Wn=[500 700]; N=3;

[Ns,Ds]=cheby2(N, Rs, Wn,'s');
[Nsk,Dsk]=ellip(N, Rp, Rs, Wn,'s');

[T,w]=freqs(Ns,Ds);
[Tk,wk]=freqs(Nsk,Dsk);

m=20*log10(abs(T));
mk=20*log10(abs(Tk));
figure(1)

plot(w,m);
grid;
ylabel('AMPLITUDATA V dB');
xlabel('CHESTOTATA V rad/s');
title('NCHF PO CHERBISHEV');
figure(2)

plot(wk,mk);
grid;
ylabel('AMPLITUDATA V dB');
xlabel('CHESTOTATA V rad/s');
title('NCHF PO KAUER');

fi = angle(T);
fi1 = unwrap(angle(T));
fi2 = fi1*180/pi;

figure(3);

plot(w,fi2);
xlabel('CHESTOTA (rad/s)');
ylabel('FCHH(gradusi)');
title('FAZOVO-CHESTOTNA HARAKTERISTIKA CHERBISHEV');
grid;

figure(4);
plot(wk,fi2);
xlabel('CHESTOTA (rad/s)');
ylabel('FCHH(gradusi)');
title('FAZOVO-CHESTOTNA HARAKTERISTIKA KAUER');
grid;

figure(5);
pzmap(Ns,Ds);
title('PND CHERBISHEV');
grid;

figure(6);
pzmap(Nsk,Dsk);
title('PND CHERBISHEV');
grid;