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);

m=20*log10(abs(T));

[Tk,wk]=freqs(Nsk,Dsk);

mk=20*log10(abs(Tk));

figure(1)

plot(w,m);
ylabel('AMPLITUDATA V dB');
xlabel('CHESTOTATA V rad/s');
title('NCHF PO CHERBISHEV(BLUE) I KAUER (RED)');

hold on;

plot(wk,mk,'r');
grid;


hold off;

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

fik = angle(Tk);
fi1k = unwrap(angle(Tk));
fi2k = fi1k*180/pi;

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


hold on

plot(wk,fi2k,'r');


hold off

figure(3)
subplot(1,2,1);
pzmap(Ns,Ds);
title('PND INV CHEBISHEV');

grid;
subplot(1,2,2);
pzmap(Nsk,Dsk);

title('PND KAUER');

grid;