093

Współpraca z plikami zewnętrznymi 93

if j==0 Yzw(i,i)=Yzw(i,i)+yg; end

end

% macierz impedancji zwarciowych Zzw=inv(Yzw) ;

% Obliczenia zwarciowe z wyprowadzaniem wyników na ekran for i = l: n

k=wezly(i , 1) ; Un=wezly(i,2); fprintf('\n');

^ * \r"l****¹*^r*********'************ ' J .

fprintf(' \n    ZWARCIE 3-fazowe W WEZLE k = %d',k);

^    J .

p3 = sqrt(3.0) ;

fprintf('\n Napięcie znamionowe wezla Un = %8.2f kV',Un);

% Przeliczenie impedancji zastępczej % na napięcie znamionowe w p. zwarcia tn=wezly(i,3) ;

Zkk=Zzw ( i , i) /tn^/N2;

Rkk=real(Zkk); Xkk=imag(Zkk);

fprintf('\n Rezystancja własna wezla Rkk = %8.2f omy',Rkk); fprintf(¹\n Reaktancja własna wezla Xkk = %8.2f omy’,Xkk); %{zwarcie 3-fazowe}

c=l.1;

Ik=c*Un/p3/abs(Zkk);

fprintf('\n PRĄD ZWARCIOWY POCZĄTKOWY Ik = %8.2f kA',Ik'); Sk=p3 *Un*Ik;

fprintf('\n MOC ZWARCIOWA    Sk = %8.2f MVA',Sk');

% zapamiętanie wyników obliczeń w postaci wektorów

tnrw(i)=k; tUn(i)=Un;

tRkk(i)=Rkk; tXkk(i)=Xkk;

tlk(i)=lk; tSk(i)=Sk;

end

% zapisywanie wyników obliczeń w m-pliku o nadanej nazwie

zapwyn(tnrw,tUn,tRkk,tXkk,tik,tSk);

et=etime(clock,tO); % całkowity czas obliczeń fprintf(. . .

¹ \n= = = = = === = = = = = = = = = = = = = = = = = === = = = = = = = = = = = = = = = = = = = = = = = = = = = = ' ) ; fprintf(    '\n! Czas po zakończeniu analizy: %.2f sekund !', et);

fclose('a11 ¹ ) ;

fprintf('\n\n*** KONIEC OBLICZEŃ ZWARCIOWYCH ***\n\n'); return;

Wyniki obliczeń zostały zapisane w pliku zwarwyn.m o następującej postaci:

function [z3fwyniki]=zwarwyn

% Wyniki obliczeń zwarciowych w postaci funkcji z3 fwyniki=[

% nrw    Un kV    Rkk om    Xkk om    Ik kA    Sk MVA