175

Analiza stabilności lokalnej i globalnej 175

Analiza stabilności lokalnej i globalnej 175

sEqb	= l/Tdob*(Eqp -	(Eqb - Id*(xdp	- xdb))	)
sEdb	= l/Tqob*(Edp -	(Edb + Iq*(xqp	- xqb))	)
sEqp	= 1/Tdop*(Ef	(Eqp - Id*(xd	- xdp))	)
sEdp	= l/Tqop*(0	(Edp + Iq*(xq	- xqp))	)
1 równania ARN i wzbudnicy
Uf =	x3+x4;

% ograniczenie napięcia wzbudnicy Uf if Uf < VRmin Uf = VRmin; end if Uf > VRmax Uf = VRmax; end % ARN

SXl = 1 / TU 1 * ( KU 1 * uh - xl);

SX2 = 1/TU2*(KU2*xl - x2); sx3 = 1/TU3*(KU3*x2 - x3);

SX4 = 1/TU4*(KU4*KU2/TU2*xl - KU4/TU2*x2 - x4);

% równanie różniczkowe sem wzbudzenia Ef SdEf = 1/TE*(KE*Uf - dEf);

%wektor pochodnych zmiennych stanu

sx= [sdr;sdw;sEqb;sEdb;sEqp;sEdp;sdEf;sxl;sx2;sx3;sx4;] tis=t; % wyprowadzanie kontrolnych wyników na ekran eps=ie-4 ;

2)<eps)|(abs(tis- 3)<eps) 5)<eps)](abs(tis- 6)<eps) 8)<eps)j(abs(tis- 9)<eps)

CZEKAJ ... t = %9.4f s'

tis) ;

warunek = tis==0 ... |(abs(tis- 1)<eps)|(abs(tis-|(abs(tis 4)<eps)|(abs(tis-|(abs(tis- 7)<eps)|(abs(tis-if warunek

fprintf('\n .trwa całkowanie fprintf ( . . .

' \ndr=%5.2 frad.

dw=%7.4 frad/s,Eqb=%7.4 f,Edb=%7.4f,Eqp=%7.4f,Edp=%7. 4f ‘.... dr, dw,Eqb,Edb,Eqp,Edp); fprintf(...

' \nlg=%7.4f,Ug=%7.4f,uh=%9.6f,Ef=%7.4f,Pm=%7.4f,Pg=%7.4f,Qg=%7.4f•...

, Ig,Ug,uh,Ef,PmO, P, Q) ;

end

return

Funkcja czytająca dane, sterującą obliczeniami i wyprowadzaniem wyników ma następującą postać:

function gs

global ws xq xd xqb xdb xqp xdp Tqob Tdob Tqop Tdop Tm

global KU1 TUI KU2 TU2 KU3 TU3 KU4 TU4 KE TE

global Uref Efmin Efmax VRmin VRmax

global R X Us PmO EfO

global tzwpocz tzwkonc

fd = fopen(’gsout.m','wt');

% czytanie danych generatora asynchronicznego wdold = cd ;

[fname, ścieżka]~u iget f i 1e('*.*', . ..

'WYBIERZ, m plik z danymi GEN+ARN gsdat.m'); eval( [’cd(’ ’ ' ,sc i ezka, ' ' ' ) ’ ] ) ; pli kdat = strtok(fname [Us,R , X , P, Q,gen,arn]= feval (pli kdat);

% GEN