081

Rozwiązywanie zadań opisanych równaniami nieliniowymi 81

Układ równań napięciowo-mocowych można rozwiązać za pomocą fsolve. Podczas pisania programu wygodnie jest dokonać podstawień pozwalających wykorzystać standardowe sposoby wywołania funkcji fsolve.

W naszym przypadku będą to następujące podstawienia:

4D = e, 42) =/, p\ = G, pl = B, _P3 = P, pĄ = Q.

W celu rozwiązania zadania za pomocą funkcji fsolve utworzono trzy następujące funkcje:

•    rmf() - funkcja definiująca układ równań nieliniowych,

•    rmg() - funkcja obliczająca gradient równań węzłowych,

•    rmobl() - funkcja sterująca obliczeniami

function f = rmf (x, G, B, P, Q)

% równania wezlowe mocy

f(l)= G*x(l)~2 + G*x(2)~2 - G*x(1) - B*x(2) + P; f (2) =-B*x (1)^2 - B*x(2)^,'2 + B*x(l) - G*x(2) + Q ; return

function df = rmg (x, G, B, P, Q)

% gradient rownan węzłowych mocy odbioru % f(l)= G*x(l)'^s2 + G*x(2)^A2 - G*x (1) - B*x(2) + P;

% f(2)=-B*x(1)"2 - B*x(2)^2 + B*x(l) - G*x(2) + Q; macJac=[ 2*G*x(l)-G    2*G*x(2)-B

-2 *B*x(1)+B    -2 *B*x(2)-G] ;

df = macJac¹; %gradient równy transponowanej m. Jacobiego return

function [Q,U]    = rmobl (R, X, Podb, Qodb, Sb, Un)

% function [Q,U]    = rmroz(P,Q,Sb,Un)

do Qmax, kV

Mvar

% wyznacza charakterystykę Q-U mocy biernej % dostarczanej linia wysokiego napięcia % R, X - rezystancja i reaktancja linii w omach, % P - moc czynna odbioru, MW % Q - moc bierna odbioru, Mvar % Sb - moc bazowa, MVA % Un - napięcie znamionowe linii, kV % Q - wektor mocy biernej zmienianej od 0

U - wartość skuteczna napięcia odbioru,

if nargin if nargin if nargin if nargin if nargin if nargin < j=sqrt(-1) ; Zb=Un^2/Sb ;

Un=220;

Sb = 4 84; Podb=242; Qodb=121;

end

end

end

end

R= 10 ; X=2 0 ;

end

end

impedancja bazowa