Szabályozástechnika Matlab Gyakorlatok, Villamosmérnöki - Index of

Soros PID Kompenzáció Hetthéssy Jenő, Bars Ruth, Barta András, 2005
Írjunk Matlab programot (hozzunk létre egy m fájlt, ez egy szövegfájl m kiterjesztéssel). Nyissunk meg
egy új m-fajlt a Matlab file menűjéből. Írjuk a Matlab utasításokat az üres fájlba. Mentsük el a fájlt: Save
As, C:/Matlab/work/myfile.m
A program meghívásához egyszerűen írjuk be a program nevét kiterjesztés nélkül a Matlab parancs
ablakába.
» myfile
Az alábbi Matlab program szolgálja a szabályozó tervezést.
clear
s=tf('s');
P=1/((1+10*s)*(1+s)*(1+0.2*s
))
P=zpk(P)
Cp=1
Cpi=(1+10*s)/(10*s)
Cpd=(1+s)/(1+0.1*s)
Cpid=Cpi*Cpd
Cpid=zpk(Cpid)
[mag,phase,w]=bode(Cp*P);
kp=margin(mag,phase-60,w)
Cp=kp*Cp;
[mag,phase,w]=bode(Cpi*P);
kpi=margin(mag,phase-60,w);
Cpi=kpi*Cpi;
[mag,phase,w]=bode(Cpd*P);
kpd=margin(mag,phase-60,w);
Cpd=kpd*Cpd;
[mag,phase,w]=bode(Cpid*P);
kpid=margin(mag,phase-60,w);
Cpid=kpid*Cpid;
figure(1),step(Tp,'r',Tpi,'b'
,Tpd,'g',Tpid,'y')
figure(2),step(Up,'r',Upi,'b'
)
figure(3),step(Upd,'g',Upid,'
y')
t=0:0.05:10;
yp=step(Tp,t);
ypi=step(Tpi,t);
ypd=step(Tpd,t);
ypid=step(Tpid,t);
ysp=dcgain(Tp)
yspi=dcgain(Tpi)
yspd=dcgain(Tpd)
yspid=dcgain(Tpid)
ep=1-ysp
epi=1-yspi
epd=1-yspd
epid=1-yspid
ytp=(max(yp)-ysp)/ysp
ytpi=(max(ypi)-yspi)/yspi
ytpd=(max(ypd)-yspd)/yspd
ytpid=(max(ypid)-yspid)/yspid
Tp=feedback(Cp*P,1);
Tpi=feedback(Cpi*P,1);
Tpd=feedback(Cpd*P,1);
Tpid=feedback(Cpid*P,1);
Up=feedback(Cp,P);
Upi=feedback(Cpi,P);
( )
up=step(Up,t);
upi=step(Upi,t);
upd=step(Upd,t);
upid=step(Upid,t);
upim=max(upi)
updm=max(upd)
upidm=max(upid)
A szabályozási kör viselkedése vizsgálható Simulink blokk-diagram megépítésével és futtatásával az adott
szakasszal és a megtervezett szabályozókkal.
41