„SzabtechLabZH” változatai közötti eltérés

s = zpk('s); <br />

P = 2 / ( (1+s)*(1+4*s)*(1+8*s)  ); <br />

Ts = 0.5; ''%ez a mintavételezési idő, adott a feladat szövegében'' <br />

sysd = c2d(P, Ts, 'zoh'); ''%ez alakítja át diszkrétre. A 'zoh' adja meg, hogy zérusrendű a tartószerv.'' <br />

s=zpk('s');<br />

z=zpk('z',0.5);<br />

P=2/((1+s)*(1+4*s)*(1+8*s)) ''% átviteli függvény'' <br />

ts=0.5; ''% mintaveteli ido'' <br />

Pz=c2d(P,ts); ''% diszkrét átviteli függvény'' <br />

Pz ''% zérus-pólus alak'' <br />

kc=1; ''%egység körerősítés először'' <br />

[zd,pd,kd] = zpkdata(Pz,'v'); ''%diszkrét folyamat zérusai pólusai'' <br />

Cz=kc*(((z-pd(1))*(z-pd(2)))/((z-1)*z)); ''%diszkrét szab. átv. fv.'' <br />

Lz = Cz*Pz<br />

Lz=minreal(Lz,0.001);<br />

w=logspace(-1,1,500); ''%saját logaritmikus skála beállítása'' <br />

[mag,phase]=bode(Lz,w); ''%diszkrétből számolt bode diagram'' <br />

[kc,pm,wg,wc]=margin(mag,phase-60,w); ''%60 fokos fázistöbblet'' <br />

Cz=Cz*kc;<br />

Lz=Cz*Pz; ''%a nyitott kör diszkrét átv. fv.'' <br />

Lz=minreal(Lz,0.001);<br />

Tz= Lz/(1+Lz); ''%diszkrét zárt rendszer átv. fv.'' <br />

Ts=d2c(Tz,'tustin') ''%folytonos zárt rendszer átv. fv.'' <br />

figure(2),step(Ts),grid ''%zart rendszer kimeneti jele'' <br />

y=step(Ts);<br />

ymax = max(y) ''%a kimenet maximális értéke'' <br />

ys=dcgain(Ts) ''%állandósult érték'' <br />

yt=(ymax-ys)/ys*100 ''%túllövés értéke százalékban'' <br />

es=1-ys ''%az előírt 1-tól való eltérés'' <br />

Uz=Cz/(1+Lz); ''%zárt rendszer beavatkozó jele'' <br />

Uz=minreal(Uz,0.001); ''%egyszerűsítések elvégzése'' <br />

figure(3),step(Uz),grid ''%beavatkozó jel kirajzolása'' <br />

u=step(Uz);<br />

umax = max(u) ''%maximális érték kiírása'' <br />