clear all;
close all;

p=[2*exp(j*pi/2) 2*exp(-j*pi/2)];
n=[1 -1];
cp=poly(p);
cn=poly(n);
zplane(cn,cp);
title('rozlozenie polov prenosovej funkcie nestabilneho systemu');

figure;
freqz(cn,cp);
title('frekvencne charakteristiky nestabilneho systemu');
h=filter(cn,cp,[1 zeros(1,10)]);

figure;
bar(h);
title('impulzova charakteristika nestabilneho systemu');
ps=[(1/2)*exp(-j*pi/2) (1/2)*exp(j*pi/2)];
cps=4*poly(ps);

figure;
zplane(cn,cps);
title('rozlozenie polov prenosovej funkcie stabilneho systemu');

figure;
freqz(cn,cps);
title('frekvencne charakteristiky stabilneho systemu');
hs=filter(cn,cps,[1 zeros(1,100)]);

figure;
bar(hs);
title('impulzova charakteristika stabilneho systemu');