165

Analiza stabilności lokalnej i globalnej 165

if nargin<2

Eprim=1.2; %sem przejściowa GENERATORA if nargin<l

k=l; % wybór Dpu suwakiem

end

end

end

end

end

end

end

end

end

fd=fopen(¹stabout.m’,'wt'); % otwarcie pliku dla wyników sind=Pm*X/(Eprim*Us);

dO=asin(sind); % kat wirnika w stanie ustalonym Kl=Eprim*Us/X*cos(dO);

ws=100*pi; % prędkość kątową synchroniczna

wm=sqrt(KI*ws/Tm); % pulsacja własnych kołysań wirnika

fm=wm/2/pi; % częstotliwość własnych klysan wirnika

set(rys,'Enable',' on ' ) ;

if k==l

Dpusuw=get(vslider,'value');

set (vedit,'string',['Dpu=',num2str(Dpusuw)]);

Dpu=Dpusuw; end

if k= = 2 Dpu=Dpusuw;

psi=Dpu/2/wm/Tm; % współczynnik tłumienia kołysań

fprintf(fd, fprintf(fd, fprintf(fd, fprintf(fd, fprintf(fd, fprintf(fd, fprintf(fd, fprintf(fd,

\n Parametry wejściowe układu; GENERATOR - SYSTEM'); \n    Eprim    =    %8.4f    pu',Eprim);

\n    Us    =    %8.4f    pu',Us);

\n    X    =    %8.4f    pu',X);

\ n    Pm    =    % 8.4 f    pu',Pm);

\n    Tm    =    % 8.2 f    s ' ,Tm) ;

\n    Dpu    =    % 8.2 f    pu',Dpu) ;

. 2 f pu

dO =	%8	. 2 f	st^1	,ddOst)
KI =	%8	. 4 f	pu'	, KI) ;
l wm -	%8	. 2 f	rad/s',wm)
fm =	%8	. 2 f	Hz '	, fm) ;
psi =	%8	. 4 f	pu'	,psi);

\n Wyniki analizy kołysań wirnika GENEARATORA'); dd0st=2; dd=ddOst*pi/180; fprintf(fd,'\n Kat początkowy wirnika fprintf(fd,'\n Moc synchronizującą fprintf(fd,'\n Pulsacja kołysań wirnika wm fprintf(fd,'\n Czestotliowsc kołysań fprintf(fd,'\n Współczynnik tłumienia

%MACIERZ STANU

A=[0    1; -wm -2*psi*wm];

%Wartosci własne macierzy stanu A lambda=eig(A);

Rell=real(lambda(1)); Imll=imag(lambda(1)); Rel2=real(lambda(2)); Iml2=imag(lambda(2)); subplot (212) ;

plot(Rell,Imll,'bO',Re12,Iml2,'rO'); text(Rell+0.2,Imll+0.2,'lambda(1)' ) ; text (Rel2 + 0.2 , Im.12 + 0.2 , ' lambda (2 ) ' ) ;