echo on, clc, clear, close all
n=[2.5 0 1];
m=[3 4.2 2];
sys = tf(n,m)
pause

%transmitancja
disp('G(s)')
printsys(n,m), sys
pause

%funkcja tf2zip
[z,b,w] = tf2zp(n,m)
[nzp,mzp]=zp2tf(z,b,w)
printsys(nzp,mzp)
pause

%residue
[r,p,k]=residue(n,m)
[nres,mres]=residue(r,p,k)
printsys(nres,mres)
pause

%zestawienie
[z,b];
[r,p];
[w,k];
pause

%odp na step i impulse
t=[0:0.01:10];
yi = impulse(n,m,t);
ys = step(n,m,t);
ysres = step(nres,mres,t);
figure(1), plot(t,yi,t,ys), hold on,
grid, axis([0 10 -1.5 1.5])
legend('impulse()','step()')

echo off