Vaja št.: 1
Datum: 12.11.2001
Besedilo naloge:ČASOVNI ODZIV ELEKTRIČNIH VEZIJ
Za električni vezji na sliki a in b izračunajte časovne odzive za podana območja spreminjanja upornosti R oz. R1 ter ugotovite pri kateri vrednosti R oz R1 ( spodnji ali zgornji) se vezje približuje meji stabilnosti. Primerjajte časovne odzive.
PROGRAM
% Štumberger Andrej 12.11.2001; Regulacijska tehnika_Rač_Vaja
% izris lege korenov
% R=4,R=10, R=14, R=20, R=40
L=10e-3 % induktivnost [H]
C=100e-6 % kapacitivnost [F]
st=1 %definicija stevca
im=[L*C R*C 1] %definicija imenovalca
printsys(st,im,'s') %izpiše sistem
damp(im) %izračuna dušenje
figure(1),hold on %nariše graf in zadrži vrednosti
axis([-1100,-100,-1000,1000]) %določimo osi
pzmap(st,im) %nariše pole
%R nastavljamo v delovnem oknu
A=100*exp((-2e-1*pi)/(sqrt(1-(2e-1))))
A =49.5355
% Štumberger Andrej 12.11.2001; Regulacijska tehnika_Rač_Vaja
% R=4,R=10, R=14, R=20, R=40
L=10e-3 % induktivnost [H]
C=100e-6 % kapacitivnost [F]
st=1 %definicija stevca
im=[L*C R*C 1] %definicija imenovalca
printsys(st,im,'s') %izpiše sistem
damp(im) %izračuna dušenje
figure(1),hold on %nariše graf in zadrži vrednosti
axis([-1100,-100,-1000,1000]) %določimo osi
pzmap(st,im) %nariše pole
%R nastavljamo v delovnem oknu
figure(2),hold on
step(st,im)
% Regulacijska tehnika_Rač_Vaja
% R=4,R=10, R=14, R=20, R=40
% izris lege korenov vezja b
R2=100e3 %definicija R2
%definicija prenosne funkcije operacijskega ojačevalnika
cl1=[1e-1 1] % 10e-1*s+1
cl2=[1e-5 1] % 10e-5*s+1
cl3=[1e-6 1] % 10e-6*s+1
clp=conv(cl1,cl2) %cl1*cl2
imA=conv(cl3,clp) %cl1*cl2*cl3
stA=5e5 %5 *e5
%definicija prenosne funkcije v povratni vezavi
imR=1
stR=R1/(R1+R2)
%izračun in izpis prenosne funkcije vezja b iz blokovne sheme
[stz,imz]=feedback(stA,imA,stR,imR)>
printsys(stz,imz,'s')
%izris polov prenosne funkcije
figure(3)
axis([-1.1e6,-1e4,-0.3e6,0.3e6])
pzmap(stz,imz)
hold on
% Regulacijska tehnika_Rač_Vaja
% R =100, R=1k, R=2,5k, R=11k, R=20k
% izris lege korenov vezja b
R2=100e3 %definicija R2
% definicija prenosne funkcije operacijskega ojačevalnika
cl1=[1e-1 1] % 10e-1*s+1
cl2=[1e-5 1] % 10e-5*s+1
cl3=[1e-6 1] % 10e-6*s+1
clp=conv(cl1,cl2) % cl1*cl2
imA=conv(cl3,clp) % cl1*cl2*cl3
stA=5e5 % 5*e5
% definicija prenosne funkcije v povratni vezavi
imR=1
stR=R1/(R1+R2)
% izračun in izpis prenosne funkcije vezja b iz blokovne sheme
[stz,imz]=feedback(stA,imA,stR,imR)
printsys(stz,imz,'s')
% izris polov prenosne funkcije
figure(3)
axis([-1.1e6,-1e4,-0.3e6,0.3e6])
pzmap(stz,imz)
hold on
t=0:1e-6:0.7e-3
figure(4)
hold on
step(stz,imz,t)
% Regulacijska tehnika_Rač_Vaja
% R =100, R=1k, R=2,5k, R=11k, R=20k
R2=100e3 % definicija R2
% Definicija prenosne funkcije operacijskega ojačevalnika
cl1=[1e-1 1] % 10e-1*s+1
cl2=[1e-5 1] % 10e-5*s+1
cl3=[1e-6 1] % 10e-6*s+1
clp=conv(cl1,cl2) %cl1*cl2
imA=conv(cl3,clp) %cl1*cl2*cl3
stA=5e5 %5*e5
% Definicija prenosne funkcije v povratni vezavi
imR=1
stR=R1/(R1+R2)
% Izračun in izpis prenosne funkcije vezja b iz blokovne sheme
[stz,imz]=feedback(stA,imA,stR,imR)
printsys(stz,imz,'s')
% Izris polov prenosne funkcije
t=[0:0.1:4]
u=sin(8*t) % sinw*t in w=8
figure(1)
lsim(stz,imz,u,t)
grid
hold on
figure(2)
bode(stz,imz)
hold on
Pri R1=1000 Ω
| Avtor: Štumberger Andrej |
|
|
