166

166 Rozdział 13

if Rell>0 | Rel2>0

text(0.2,0.2,'Utrata Stabilności lokalnej'); end

title('Wartosci własne macierzy stanu A. Wciśnij ENTER!'); xlabel('Real'); ylabel('Imag'); grid on; pause;

% Ustalenie czasu symulacji 10 = 0 ; tmax=10;

xO=[dO; 0]; utratastab=0;

[t,x]=ode4 5( 'stabdx' ,tO,tmax,x0);

fprintf(¹\n\n .trwa całkowanie met. ode45() ... CZEKAJ !\n\n’);

nt=length(t);

fprintf(fd,'\n t,s d,st w,rad/s');

for K=1:nt

fprintf(fd,'\n %8.4f %12.6f %12.6f %12.6f',t(k),x(k,1),x(k,2) );

end

subplot(212) ; plot( Ł , X ( : , 1 ) ) ;

grid on; xlabel(’czas t, s'), ylabel('kat wirnika, st' ) ; title('Kat wirnika w stopniach - zwarcie 3 fazowe'); if utratastab

text(tpzw+tzw,dO,'Utrata Stabilności - przerwano obliczenia’);

else

text(tpzw,0.5,'Zwarcie na generatorze');

end

grid on; end

fclose('all'); return function fig = stabsym

%interfejs graficzny do analizy małych kołysań global wm psi Dpusuw vslider vedit rys Dpu format compact

'Tag' , 'Fig2');

set(gcf, 'NumberTitle' , 'of f' , . . .

¹Name¹,'Wybór Dpu suwakiem');

hl = uicontrol('Parent',hO, ...

'Tag', 'Stat icText2') ;

% zmiana wartości tłumienia za pomocą suwaka

Dpusuw= 0;

vslider=uiccntrol('style','slider',...

'units','normalized',...

'positicn' , [ .05, .084, .035,0.75], . . .

'min' , - 0.5, 'max' ,5, . . .

' value' ,Dpusuw....