082

82 Rozdział 7

Z=(R+j*X)/Zb; Y=1/Z;

G=real(Y); B=imag(Y); % parametry linii w jedn. względnych

P=Podb/Sb; Q=Qodb/Sb; % przeliczenie mocy na jednostki względne

xO=[l 0]; x=x0; % punkt startowy obliczeń

options=foptions; % opcje standardowe

% Wyznaczenie napięcia odbioru

x=fsolve('rmf',x0,options,'rmg',G,B,P,Q);

u=sqrt(x(1).~2 + x(2)~2); % wyznaczone napięcie odbioru

Urm = u*Un; % napięcie odbioru w kV

Qrm=Q*Sb; % moc bierna odbioru w Mvar

disp('Napięcie odbioru w kV wynosi:');

disp(Urm);

% Wyznaczanie charakterystyki napięciowej Q-U i=0; Q=0;

U=[]; u0=l; U(1)=u0; u=l;

Qodb=[]; Qodb(1)=Q;

Ulawina=0.65*u0; % przybliżony punkt lawiny napięcia while u > Ulawina % obliczenia az do lawiny napięcia i=i+l; % kolejny krok charakterystyki Q-U Q=Q+0.01; % przyrost mocy biernej x=fsolve('rmf¹,x0,options,'rmg',G,B,P,Q);

u=sqrt (x (1) . ~2 + x(2)^/'2); % napięcie Qodb(i)=Q*Sb; % moc bierna w Mvar Uodb(i)=u*Un; % napięcie w kV Urmi(i)=Urm; % napięcie odbioru

end

nroz=i; Qlawina=Q*Sb; Ulawina=Ulawina*Un; disp(' Liczba wyznaczonych punktów charakterystyki napięciowej ' ) ; disp(nroz);

plot(Qodb,Uodb,'b-',Qodb,Urmi,'r--');

title('Charakterystyka Q-U mocy biernej dostarczanej linia 220 kV' ) ;

xlabel('Q, Mvar'); ylabel('U, kV' ) ; legend('Krzywa Q-U','Napięcie odbioru'); hołd on;

plot(Qrm,Urm,'O');

text(Qrm,Urm+2,'Punkt pracy układu przesyłowego'); plot(Qlawina,Ulawina,'0');

text(Qlawina-150,Ulawina+1,'Punkt lawiny napięcia');

grid on;

disp('KONIEC');

return

Przebieg obliczeń numerycznych po wywołaniu funkcji rmobl jest następujący:

Napięcie odbioru w kV wynosi:

194.2949