- 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')