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

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

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

%dunkcja 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

%pd[pwoedzo ma step i pulse
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),
grid, axis([0 10 -1.5 1.5])
legend('impulse()','step()')
echo off

/////////////////////////////////////////////////

echo on
clc
clear
close all

%--------------------------------------
%dane
n1 = [4 5 6]; m1 = [1 5 4]
printsys(n1,m1)
n2=1;m2=[1 2.5];
printsys(n2,m2);
pause

%--------------------------------------
%połączenie szeregowe 
[nsz,msz]= series9n1,m1,n2,m2);
disp('gszer(s)')
printsys(nsz,msz)
pause 

%--------------------------------------
%połaczenie rownoległe
[nr,mr]= parallel(n1,m1,n2,m2);
disp('Grow(s)')
printsys(nr,mr)
pause

%--------------------------------------
%uklad z ujemnym sprzezeniem zwrotnym

[nusz,musz]=feedback(n1,m1,n2,m2);
disp('Gusz(s)')
printsys(nusz,musz)
pause

%--------------------------------------
%residua
[rsz,bsz,wsz] = residue(nsz,msz);
[rr,br,wr] = residue(nr,mr);
[rusz,busz,wusz]= residue(nusz,muszz);
pause

%--------------------------------------
%tf2zp
[zsz,psz,ksz] = tf2zp(nsz,msz);
[zr,pr,kr] = tf2zp(nr,mr);
[zusz,pusz,kusz]= tf2zp(nusz,muszz);
pause

tf2zp
%odpowiedz na step
t=[0:.01:6];
ysz=step(nsz,msz,t);
yr=step(nr,mr,t);
yusz=step(nusz,musz,t);
figure(1)
plot(t,ysz,t,yr,t,yusz) ,grid, axis ([ 0 6 0 2.5])
legend('ysz','yr','yusz')