k=1.55; %konstanta
n1=1*exp(-i*pi*1/15); % nula
n2=1*exp(-i*pi*1/25); % nula
p1=0.9*exp(-i*pi*1/8); % pol
p2=0.71*exp(-i*pi*1/99); % pol
n1z=conj(n1);
n2z=conj(n2);
p1z=conj(p1);
p2z=conj(p2);

subplot (321);
zplane(n1,p1);
subplot(322);
zplane(n2,p2);
subplot(323);
zplane(n1z,p1z);
subplot(324);
zplane(n2z,p2z);


nuly=[n1 n2 n1z n2z];
pol =[p1 p2 p1z p2z];
nuly_polynom=real(poly(nuly));
pol_polynom=real(poly(pol));

ft=abs(fft(nuly_polynom,1000))./(k*abs(fft(pol_polynom,1000)));
subplot(325)
plot([0:500]*pi/500,ft(1:501)); grid on% vykreslenie magnitudy do pi!

n = 15;
[h,t] = impz(nuly_polynom,pol_polynom,n);
subplot(326);
stem(t,h);
