%% charakterystyki czestotliwosciowe
figure(5)
nyquist(obiekt)
ylim([-0.9 0.9])
xlim([-0.9 0.9])
title 'Obiekt'
figure(6)
nyquist(modelkupf)
ylim([-0.9 0.9])
xlim([-0.9 0.9])
title 'Model'
figure(7)
bode(obiekt)
hold on
bode(modelkupf)
legend('Obiekt', 'Model')
figure(8)
margin(obiekt)
figure(9)
margin(modelkupf)
[ZapasModuluModelu, ZapasFazyModelu] =  margin(modelkupf)
[ZapasModuluObiektu, ZapasFazyObiektu] = margin(obiekt)
PasmoObiektu=bandwidth(obiekt)
PasmoModelu=bandwidth(modelkupf)
%% bieguny i zera
figure(10)
pzplot(obiekt)
title 'Obiekt'
xlim([-40 5])
figure(11)
pzplot(modelkupf)
title 'Model'