Osnove elektrotehnike II, udžbenik - Ljubomir Malešević

OSNOVE ELEKTROTEHNIKE II 123
Fazorski dijagram u uvjetima paralelne rezonancije prikazan je na Slici 6.17.
Slika 6.17 – Fazorski dijagrami – paralelna rezonancija
Rezonancijske krivulje paralelnoga titrajnog kruga prikazane su na sljedećim grafovima. Simulacija je
izvedena pomoću programskog paketa MATLAB. Kôd programa je:
% FREKVENCIJSKE KRIVULJE PARALELNOGA REZONANCIJSKOG KRUGA
format short
% izlazni format brojeva
R=1000;
% djelatni otpor
L=2.5*10^-3;
% induktivnost svitka
C=0.1*10^-6;
% kapacitivnost kondenzatora
U=100;
% narinuti napon
f=0:100:2*10^4;
% frekvencijski opseg 0–10 KHz
G=1/R;
% konduktancija
BL=j*1./(2*pi*f*L);
% induktivna susceptancija
BC=j*2*pi*f*C;
% kapacitivna susceptancija
Y=G-BL+BC;
% admitancija paralelnoga R-L-C kruga
Z=1./Y;
% impedancija paralelnoga R-L-C kruga
I=U.*Y;
% ukupna struja
IL=BL*U;
% struja kroz induktivnu susceptanciju
IC=BC*U;
% struja kroz kapacitivnu susceptanciju
IG=G*U;
% struja kroz konduktanciju
IX=I-IG;
% struja kroz ukupnu susceptanciju
Faza_rad=atan((abs(BL)-abs(BC))./G); % fazni kut u radijanima
Faza_stup=Faza_rad*180/pi; % fazni kut u stupnjevima
% Ovisnost reaktancija o frekvenciji
figure (1)
Nula=0*f;
% os apscisa
G=plot(f,Nula,'b',f,-abs(BC),f,abs(BL),f,abs(BL)-abs(BC)); % grafovi BL, BC, B
set(G,'linewidth',2);
% debljina krivulje
axis([0,(2*10^4),-20*10^-3,20*10^-3]);
% koordinatne osi
title('OVISNOST SUSCEPTANCIJA O FREKVENCIJI') % naslov
xlabel('Frekvencija f[Hz]')
% oznaka X-osi