k=0,110; %konstanta
n1=1.008*exp(-i*1.695); % nula
n2=1.267*exp(-i*2.884); % nula
p1=0.736*exp(-i*1.138); % pol
p2=0.263*exp(-i*-2.177); % 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);

subplot(327);
zplane((n1z+n2z),(p1z+p2z));