034

34 Rozdział 3

d-j2)S ₂    0.5-j2)S

(1 -j 2 )S

(0.5-j2.5)S

Rys. 3.1. Liniowy obwód elektryczny z numerami węzłów i admitancjami gałęzi.

Węzeł o numerze 0 jest węzłem odniesienia

Do tworzenia macierzy admitancji węzłowych wykorzystać następujące instrukcje:

j=sqrt(-l); ygal=g+j*b; nbr=length(ygal); for k=l:nbr

i=wp(k); j=wk(k); y=ygal(k); if i*j > 0

Y(i, i)=Y(i,i)+y; Y(i,j)=-y;

Y(j, j)=Y(j,j)+y; Y(j,i)=-y;

end

if i==0 Y(j,j)=Y(j,j)+y; end if j==0 Y(i,i)=Y(i,i)+y; end

end

Utworzoną zespoloną macierz admitancji zapisać na dysku w pliku macY.mat.