Juni 2021

Elektrotehnički
fakultet
Univerzitet Sarajevo

SREES
Laboratorijska vježba 3

Student Asja Avdić, 17458

1.1
clc
clearvars
sb = 10 e6 ;
v 1 b l l = 110 e3 ;
v 2 b l l = 20 e3 ;
v 3 b l l = 0 . 4 e3 ;
bVals1 = g e t B a s e V a l u e s ( sb , v 1 b l l )

bVals1 = s t r u c t with f i e l d s :
v b l n : 6 . 3 5 0 9 e+04
ib : 157.4592
zb : 4 0 3 . 3 3 3 3

bVals2 = g e t B a s e V a l u e s ( sb , v 2 b l l )
bVals2 = s t r u c t with f i e l d s :
v b l n : 1 . 1 5 4 7 e+04
ib : 866.0254
zb : 1 3 . 3 3 3 3

bVals3 = g e t B a s e V a l u e s ( sb , v 3 b l l )
bVals3 = s t r u c t with f i e l d s :
vb ln : 230.9401
i b : 4 . 3 3 0 1 e+04
zb : 0 . 0 0 5 3

1.2
Model generator u p.u.
e g l l = 112 e3 ;
f i = 0;
r g = 0 . 1 e −3;
xg = 0 . 5 e −3;
[ Yg , Ig , Eg ] = g e n e r a t o r M o d e l ( e g l l , f i , rg , xg , bVals1
)

1
Yg = 3 3 complex
105
1.5513 − 7.7564 i 0.0000 + 0.0000 i 0.0000 + 0.0000 i
I g = 3 1 complex
105
1.5795 − 7.8974 i
−7.6291 + 2 . 5 8 0 8 i 6 . 0 4 9 6 + 5 . 3 1 6 6 i
Eg = 3 1 complex 1 . 0 1 8 2 + 0 . 0 0 0 0 i −0.5091 − 0 . 8 8 1 8 i
−0.5091 + 0 . 8 8 1 8 i

1.3
Model transformatora 1-2 u p.u.
Sn = 10 e6 ;
Un1 = 110 e3 ; % Line v o l t a g e Un2 = 20 e3 ; % Line v o l t a g e
% Zero s e q u e n c e
uk 0 = 1 3 . 2 ; Pcu 0 = 2 . 7 5 e3 ; I 0 0 = 0 . 0 8 3 ; P f e 0 =
1 . 8 3 3 e3 ;
% P o s i t i v e sequence
uk 1 = 1 2 ; Pcu 1 = 3 . 3 e3 ; I 0 1 = 0 . 1 ; P f e 1 = 2 . 6 4 e3 ;
c = 5;
zg1 = 1 0 0 ;
zg2 = 0 ;
n = Un1/Un2 ; c o n n e c t i o n = ’Yd ’ ;
[ z s 0 , y s h 0 ] = xfmrSequenceParams ( Un1 , Sn , uk 0 ,
Pcu 0 , I 0 0 , Pfe 0 , bVals1 . zb ) ; [ z s 1 , y s h 1 ] =
xfmrSequenceParams ( Un1 , Sn , uk 1 , Pcu 1 , I 0 1 ,
Pfe 1 , bVals1 . zb ) ;
[ Y12 0 , Y12 1 , Y12 2 ] = xfmrModelSymm ( z s 0 , ysh 0 , z s 1
, ysh 1 , zg1 / bVals1 . zb , zg2 / bVals1 . z [ Y12 ] =
xfmrModel ( Y12 0 , Y12 1 , Y12 2 )
Y12 = 6 6 complex
0.354080744661644 − 2.037531452102654 i
0.346397855772755 + 0.740396575963023 i
0.346397855772755 + 0.740396575963023 i
0.354080744661644 − 2.037531452102654 i

2
 0.346397855772755 + 0.740396575963023 i
0.346397855772755 + 0.740396575963023 i
0.004410314556310 − 1.603744683554981 i
−0.004410314556309 + 1 . 6 0 3 7 4 4 6 8 3 5 5 4 9 8 0 i
0.000000000000000 + 0.000000000000000 i
0.004410314556309 − 1.603744683554981 i
−0.004410314556310 + 1 . 6 0 3 7 4 4 6 8 3 5 5 4 9 8 1 i
−0.000000000000000 + 0 . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 i

1.4
Model transformatora 4-5 u p.u.
Sn = 10 e6 ;
Un1 = 110 e3 ; % Line v o l t a g e Un2 = 20 e3 ; % Line v o l t a g e
% Zero s e q u e n c e
uk 0 = 1 3 . 2 ; Pcu 0 = 2 . 7 5 e3 ; I 0 0 = 0 . 0 8 3 ; P f e 0 =
1 . 8 3 3 e3 ;
% P o s i t i v e sequence
uk 1 = 1 2 ; Pcu 1 = 3 . 3 e3 ; I 0 1 = 0 . 1 ; P f e 1 = 2 . 6 4 e3 ;
c = 5;
zg1 = 1 0 0 ;
zg2 = 0 ;
n = Un1/Un2 ; c o n n e c t i o n = ’Yd ’ ;
[ z s 0 , y s h 0 ] = xfmrSequenceParams ( Un1 , Sn , uk 0 ,
Pcu 0 , I 0 0 , Pfe 0 , bVals1 . zb ) ; [ z s 1 , y s h 1 ] =
xfmrSequenceParams ( Un1 , Sn , uk 1 , Pcu 1 , I 0 1 ,
Pfe 1 , bVals1 . zb ) ;
[ Y12 0 , Y12 1 , Y12 2 ] = xfmrModelSymm ( z s 0 , ysh 0 , z s 1
, ysh 1 , zg1 / bVals1 . zb , zg2 / bVals1 . z [ Y12 ] =
xfmrModel ( Y12 0 , Y12 1 , Y12 2 )
Y12 = 6 6 complex
0.354080744661644 − 2.037531452102654 i
0.346397855772755 + 0.740396575963023 i
0.346397855772755 + 0.740396575963023 i
0.354080744661644 − 2.037531452102654 i
0.346397855772755 + 0.740396575963023 i

3
 0.346397855772755 + 0.740396575963023 i
0.004410314556310 − 1.603744683554981 i
−0.004410314556309 + 1 . 6 0 3 7 4 4 6 8 3 5 5 4 9 8 0 i
0.000000000000000 + 0.000000000000000 i
0.004410314556309 − 1.603744683554981 i
−0.004410314556310 + 1 . 6 0 3 7 4 4 6 8 3 5 5 4 9 8 1 i
−0.000000000000000 + 0 . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 i

1.5
Model linija
r 0 = 0 . 1 8 9 5 ; x0 = 0 . 3 6 8 3 ; c0 = 1 . 5 e −9; r 1 = 0 . 1 2 1 8 ; x1
= 0 . 1 8 9 5 ; c1 = 2 . 3 e −9; l 2 3 = 6 . 5 ; l 2 4 = 1 2 ; l 3 4 = 9 ;
[ Y23 0 , Y23 1 , Y23 2 ] = lineModelSymm ( r0 , x0 , c0 , r1 ,
x1 , c1 , l 2 3 , bVals2 ) ; [ Y23 ] = l i n e M o d e l P h a s e ( r0 , x0 ,
c0 , r1 , x1 , c1 , l 2 3 , bVals2 )
Y23 = 6 6 complex
4.037625791956102 − 6.574656515053729 i
−0.885888669056174 + 1 . 0 8 5 4 5 9 8 7 7 8 5 9 5 7 4 i
−0.885888669056175 + 1 . 0 8 5 4 5 9 8 7 7 8 5 9 5 7 4 i
4.037625791956102 − 6.574656515053729 i
−0.885888669056174 + 1 . 0 8 5 4 5 9 8 7 7 8 5 9 5 7 3 i
−0.885888669056175 + 1 . 0 8 5 4 5 9 8 7 7 8 5 9 5 7 3 i
−4.037625791956102 + 6 . 5 7 4 6 8 4 1 9 5 9 7 5 6 6 5 i
0.885888669056174 − 1.085463508144418 i
0.885888669056175 − 1.085463508144418 i
−4.037625791956102 + 6 . 5 7 4 6 8 4 1 9 5 9 7 5 6 6 6 i
0.885888669056174 − 1.085463508144417 i
0.885888669056175 − 1.085463508144417 i

[ Y24 0 , Y24 1 , Y24 2 ] = lineModelSymm ( r0 , x0 , c0 , r1 ,

x1 , c1 , l 2 4 , bVals2 ) ; [ Y24 ] = l i n e M o d e l P h a s e ( r0 , x0 ,
c0 , r1 , x1 , c1 , l 2 4 , bVals2 )
Y24 = 6 6 complex
2.187047303976222 − 3.561236169579655 i
−0.479856362405427 + 0 . 5 8 7 9 5 2 6 9 8 1 8 0 5 6 6 i

4
−0.479856362405428 + 0 . 5 8 7 9 5 2 6 9 8 1 8 0 5 6 6 i
2.187047303976222 − 3.561236169579655 i
−0.479856362405428 + 0 . 5 8 7 9 5 2 6 9 8 1 8 0 5 6 6 i
−0.479856362405428 + 0 . 5 8 7 9 5 2 6 9 8 1 8 0 5 6 5 i
−2.187047303976222 + 3 . 5 6 1 2 8 7 2 7 2 8 2 0 1 5 2 i
0.479856362405427 − 0.587959400244893 i
0.479856362405428 − 0.587959400244893 i
−2.187047303976222 + 3 . 5 6 1 2 8 7 2 7 2 8 2 0 1 5 3 i
0.479856362405428 − 0.587959400244893 i
0.479856362405428 − 0.587959400244893 i

[ Y34 0 , Y34 1 , Y34 2 ] = lineModelSymm ( r0 , x0 , c0 , r1 ,

x1 , c1 , l 3 4 , bVals2 ) ; [ Y34 ] = l i n e M o d e l P h a s e ( r0 , x0 ,
c0 , r1 , x1 , c1 , l 3 4 , bVals2 )
Y34 = 6 6 complex
2.916063071968295 − 4.748344702996495 i
−0.639808483207236 + 0 . 7 8 3 9 4 0 8 4 0 4 4 4 9 4 5 i
−0.639808483207237 + 0 . 7 8 3 9 4 0 8 4 0 4 4 4 9 4 5 i
2.916063071968295 − 4.748344702996496 i
−0.639808483207236 + 0 . 7 8 3 9 4 0 8 4 0 4 4 4 9 4 4 i
−0.639808483207237 + 0 . 7 8 3 9 4 0 8 4 0 4 4 4 9 4 4 i
−2.916063071968295 + 4 . 7 4 8 3 8 3 0 3 0 4 2 6 8 6 9 i
0.639808483207236 − 0.783945866993191 i
0.639808483207237 − 0.783945866993191 i
−2.916063071968295 + 4 . 7 4 8 3 8 3 0 3 0 4 2 6 8 7 0 i
0.639808483207237 − 0.783945866993190 i
0.639808483207237 − 0.783945866993190 i

1.6
Model potrosaca
Ip3 = 1/ bVals2 . i b ∗ [
1 5 3 . 7 5 4 1 7 ∗ exp(−1 i ∗ 1 7 4 . 0 3 ∗ p i /180) ; 1 6 1 . 9 8 9 7 9 ∗ exp
( 1 i ∗ 6 9 . 0 9 ∗ p i /180) ; 1 6 5 . 6 2 3 8 9 ∗ exp(−1 i ∗ 5 4 . 8 3 ∗
p i /180) ] ;
Ip5 = 1/ bVals3 . i b ∗ [

5
1 0 0 3 . 8 6 7 2 0 1 ∗ exp(−1 i ∗ 2 3 . 0 6 ∗ p i /180) ; 8 9 2 . 4 8 0 1 6 ∗
exp(−1 i ∗ 1 4 3 . 9 5 ∗ p i /180) ; 1 0 3 8 . 1 4 4 3 4 ∗ exp ( 1 i ∗
8 4 . 3 8 ∗ p i /180) ] ;

1.7
Model kondenzatorske baterije
c = 0 . 2 e −6;
conn = ’Y ’ ;
Y40 = shuntModel ( c , conn , bVals2 )
Y40 = 3 3 complex 10−3
0.000000000000000 + 0.837758040957278 i
0.000000000000000 + 0.000000000000000 i
0.000000000000000 + 0.000000000000000 i
0.000000000000000 + 0.837758040957278 i
0.000000000000000 + 0.000000000000000 i
0.000000000000000 + 0.000000000000000 i

1.8
Rjesenje u faznim koordinatama
Zm = z e r o s ( 3 ) ; Zv = z e r o s ( 3 , 1 ) ;
Y = zeros (15) ;
Y = s e t S u b m a t r i x (Y, 1 , 1 , Yg + getSubmatrix ( Y12 , 1 , 1 ) )
;
Y = s e t S u b m a t r i x (Y, 1 , 2 , getSubmatrix ( Y12 , 1 , 2 ) ) ;
Y = s e t S u b m a t r i x (Y, 2 , 1 , getSubmatrix ( Y12 , 2 , 1 ) ) ;
Y = s e t S u b m a t r i x (Y, 2 , 2 , getSubmatrix ( Y12 , 2 , 2 ) +
getSubmatrix ( Y23 , 1 , 1 ) + getSubmatrix ( Y24 Y =
s e t S u b m a t r i x (Y, 2 , 3 , getSubmatrix ( Y23 , 1 , 2 ) ) ;
Y = s e t S u b m a t r i x (Y, 2 , 4 , getSubmatrix ( Y24 , 1 , 2 ) ) ;
Y = s e t S u b m a t r i x (Y, 3 , 2 , getSubmatrix ( Y23 , 2 , 1 ) ) ;
Y = s e t S u b m a t r i x (Y, 3 , 3 , getSubmatrix ( Y23 , 2 , 2 ) +
getSubmatrix ( Y34 , 1 , 1 ) ) ;

6
Y = s e t S u b m a t r i x (Y, 3 , 4 , getSubmatrix ( Y34 , 1 , 2 ) ) ;
Y = s e t S u b m a t r i x (Y, 4 , 2 , getSubmatrix ( Y24 , 2 , 1 ) ) ;
Y = s e t S u b m a t r i x (Y, 4 , 3 , getSubmatrix ( Y34 , 2 , 1 ) ) ;
Y = s e t S u b m a t r i x (Y, 4 , 4 , Y40 + getSubmatrix ( Y24 , 2 , 2 )
+ getSubmatrix ( Y34 , 2 , 2 ) + getSubmatr Y =
s e t S u b m a t r i x (Y, 4 , 5 , getSubmatrix ( Y45 , 1 , 2 ) ) ;
Y = s e t S u b m a t r i x (Y, 5 , 4 , getSubmatrix ( Y45 , Y =
s e t S u b m a t r i x (Y, 5 , 5 , getSubmatrix ( Y45 ,
I = [ I g ; Zv ; −Ip3 ; Zv ; −Ip5 ] ;
V = Y\ I ;
V1 = V( 1 : 3 ) ; V2 = V( 4 : 6 ) ; V3 = V( 7 : 9 ) ; V4 = V( 1 0 : 1 2 ) ;
V5 = V( 1 3 : 1 5 ) ;
p o l a r P r i n t (V1)
ans = 3 1 c e l l a r r a y { ’ 1 . 0 1 8 1 8 2 < −0.000012 ’ } { ’
1 . 0 1 8 1 8 2 < −120.000012 ’ } { ’ 1 . 0 1 8 1 8 2 < 1 1 9 . 9 9 9 9 8 7 ’ }
p o l a r P r i n t (V2)
ans = 3 1 c e l l a r r a y { ’ 0 . 9 2 5 9 1 5 < −159.407030 ’ } { ’
1 . 1 0 6 2 7 9 < 8 5 . 1 9 0 7 1 1 ’ } { ’ 0 . 9 4 9 6 1 0 < −27.647865 ’ }
p o l a r P r i n t (V3)
ans = 3 1 c e l l a r r a y { ’ 0 . 9 1 4 2 7 6 < −160.084166 ’ } { ’
1 . 0 9 3 6 3 8 < 8 4 . 6 1 5 0 7 2 ’ } { ’ 0 . 9 3 4 8 5 6 < −28.169555 ’ }
p o l a r P r i n t (V4)

ans = 3 1 c e l l a r r a y { ’ 0 . 9 1 7 8 7 5 < −159.873928 ’ } { ’

1 . 0 9 7 4 4 0 < 8 4 . 8 0 0 2 9 8 ’ } { ’ 0 . 9 3 9 5 0 7 < −27.995419 ’ }
p o l a r P r i n t (V5)
ans = 3 1 c e l l a r r a y { ’ 0 . 9 5 5 8 6 7 < −9.020881 ’ } { ’
0 . 9 6 2 6 0 9 < −128.426289 ’ } { ’ 0 . 9 3 8 1 7 0 < 1 1 0 . 9 4 0 4 0 7 ’ }

Ipom = Y12 ∗ [ V1 ; V2 ] ; I 1 2 = Ipom ( 1 : 3 ) ;

I 2 1 = Ipom ( 4 : 6 ) ;
Ipom = Y23 ∗ [ V2 ; V3 ] ; I 2 3 = Ipom ( 1 : 3 ) ;
I 3 2 = Ipom ( 4 : 6 ) ;
Ipom = Y24 ∗ [ V2 ; V4 ] ; I 2 4 = Ipom ( 1 : 3 ) ;
I 4 2 = Ipom ( 4 : 6 ) ;
Ipom = Y34 ∗ [ V3 ; V4 ] ; I 3 4 = Ipom ( 1 : 3 ) ;
I 4 3 = Ipom ( 4 : 6 ) ;
polarPrint ( I12 )

7
ans = 3 1 c e l l a r r a y { ’ 0 . 2 0 4 9 5 2 < −25.688449 ’ } { ’
0 . 2 0 2 2 4 0 < −141.804973 ’ } { ’ 0 . 2 1 5 4 4 0 < 9 6 . 8 6 5 1 8 5 ’ }
polarPrint ( I21 )
ans = 3 1 c e l l a r r a y
{ ’ 0 . 1 9 9 2 9 5 < 6 . 0 9 1 8 5 4 ’ } { ’ 0 . 2 1 0 0 6 8 < −111.377381 ’ } { ’
0.212660 < 124.875488 ’ }
polarPrint ( I23 )
ans = 3 1 c e l l a r r a y { ’ 0 . 1 4 5 0 6 5 < −173.950988 ’ } { ’
0 . 1 5 2 8 7 6 < 6 8 . 8 1 6 1 4 1 ’ } { ’ 0 . 1 5 5 3 8 0 < −55.000552 ’ }
polarPrint ( I32 )
ans = 3 1 c e l l a r r a y { ’ 0 . 1 4 5 0 8 1 < 6 . 0 2 6 6 6 1 ’ } { ’
0 . 1 5 2 8 9 5 < −111.207742 ’ } { ’ 0 . 1 5 5 4 0 5 < 1 2 4 . 9 7 9 4 4 8 ’ }
polarPrint ( I24 )
ans = 3 1 c e l l a r r a y { ’ 0 . 0 5 4 2 3 0 < −173.793544 ’ } { ’
0 . 0 5 7 1 9 5 < 6 8 . 1 0 5 3 4 6 ’ } { ’ 0 . 0 5 7 2 8 2 < −55.460762 ’ }
polarPrint ( I42 )
ans = 3 1 c e l l a r r a y { ’ 0 . 0 5 4 2 6 1 < 6 . 0 9 5 8 7 2 ’ } { ’
0 . 0 5 7 2 3 1 < −112.012161 ’ } { ’ 0 . 0 5 7 3 2 9 < 1 2 4 . 4 3 9 3 5 7 ’ }
polarPrint ( I34 )
ans = 3 1 c e l l a r r a y { ’ 0 . 0 3 2 4 5 9 < 5 . 7 1 6 7 4 6 ’ } { ’
0 . 0 3 4 1 6 6 < −109.577481 ’ } { ’ 0 . 0 3 5 8 4 6 < 1 2 5 . 9 9 6 1 3 7 ’ }
polarPrint ( I43 )
ans = 3 1 c e l l a r r a y { ’ 0 . 0 3 2 4 3 6 < −174.145515 ’ } { ’
0 . 0 3 4 1 4 3 < 7 0 . 5 7 1 2 7 4 ’ } { ’ 0 . 0 3 5 8 1 3 < −53.883198 ’ }

1.9
Rjesenje u simetricnim komponentama
Egsim = phase to symm ( Eg , ’V ’ ) ; I g s i m = phase to symm (
Ig , ’ I ’ ) ; Ygsim = phase to symm (Yg , ’Y ’ ) ; Ip3sim =
phase to symm ( Ip3 , ’ I ’ ) ; Ip5sim = phase to symm ( Ip5 ,
’ I ’ ) ; Y40sim = phase to symm ( Y40 , ’Y ’ ) ;
Y1 = [ Ygsim ( 2 , 2 )+Y12 1 ( 1 , 1 ) Y12 1 ( 1 , 2 ) 0 0 0 ;
Y12 1 ( 2 , 1 ) Y12 1 ( 2 , 2 )+Y24 1 ( 1 , 1 )+Y23 1 ( 1 , 1 ) Y23 1
( 1 , 2 ) Y24 1 ( 1 , 2 ) 0 ; 0 Y23 1 ( 2 , 1 ) Y23 1 ( 2 , 2 )+Y34 1
( 1 , 1 ) Y34 1 ( 1 , 2 ) 0 ;

8
0 Y24 1 ( 2 , 1 ) Y34 1 ( 2 , 1 ) Y40sim ( 2 , 2 )+Y24 1 ( 2 , 2 )+
Y34 1 ( 2 , 2 )+Y45 1 ( 1 , 1 ) Y45 1 ( 1 , 2 ) ; 0 0 0 Y45 1 ( 2 , 1 )
Y45 1 ( 2 , 2 ) ] ;
Y2 = [ Ygsim ( 3 , 3 )+Y12 2 ( 1 , 1 ) Y12 2 ( 1 , 2 ) 0 0 0 ;
Y12 2 ( 2 , 1 ) Y12 2 ( 2 , 2 )+Y24 2 ( 1 , 1 )+Y23 2 ( 1 , 1 ) Y23 2
( 1 , 2 ) Y24 2 ( 1 , 2 ) 0 ;
0 Y23 2 ( 2 , 1 ) Y23 2 ( 2 , 2 )+Y34 2 ( 1 , 1 ) Y34 2 ( 1 , 2 ) 0 ;
0 Y24 2 ( 2 , 1 ) Y34 2 ( 2 , 1 ) Y40sim ( 3 , 3 )+Y24 2 ( 2 , 2 )+
Y34 2 ( 2 , 2 )+Y45 2 ( 1 , 1 ) Y45 2 ( 1 , 2 ) ; 0 0 0 Y45 2 ( 2 , 1 )
Y45 2 ( 2 , 2 ) ] ;
I 1 = [ I g s i m ( 2 ) ; 0 ; −Ip3sim ( 2 ) ; 0 ; −Ip5sim ( 2 ) ] ; I 2 = [
I g s i m ( 3 ) ; 0 ; −Ip3sim ( 3 ) ; 0 ; −Ip5sim ( 3 ) ] ;
%p o z i t i v n a s e k v e n c a
V1 = Y1\ I 1 ;
Ipom = Y12 1 ∗ [ V1 ( 1 ) ; V1 ( 2 ) ] ; I 1 2 1 = Ipom ( 1 ) ;
I 2 1 1 = Ipom ( 2 ) ;
Ipom = Y23 1 ∗ [ V1 ( 2 ) ; V1 ( 3 ) ] ; I 2 3 1 = Ipom ( 1 ) ;
I 3 2 1 = Ipom ( 2 ) ;
Ipom = Y24 1 ∗ [ V1 ( 2 ) ; V1 ( 4 ) ] ; I 2 4 1 = Ipom ( 1 ) ;
I 4 2 1 = Ipom ( 2 ) ;
Ipom = Y34 1 ∗ [ V1 ( 3 ) ; V1 ( 4 ) ] ; I 3 4 1 = Ipom ( 1 ) ;
I 4 3 1 = Ipom ( 2 ) ;
Ipom = Y45 1 ∗ [ V1 ( 4 ) ; V1 ( 5 ) ] ; I 4 5 1 = Ipom ( 1 ) ;
I 5 4 1 = Ipom ( 2 ) ;
I 4 0 1 = Y40sim ( 2 ) ∗V1 ( 4 ) ;
%n e g a t i v n a s e k v e n c a
V2 = Y2\ I 2 ;
Ipom = Y12 2 ∗ [ V2 ( 1 ) ; V2 ( 2 ) ] ; I 1 2 2 = Ipom ( 1 ) ;
I 2 1 2 = Ipom ( 2 ) ;
Ipom = Y23 2 ∗ [ V2 ( 2 ) ; V2 ( 3 ) ] ; I 2 3 2 = Ipom ( 1 ) ;
I 3 2 2 = Ipom ( 2 ) ;
Ipom = Y24 2 ∗ [ V2 ( 2 ) ; V2 ( 4 ) ] ; I 2 4 2 = Ipom ( 1 ) ;
I 4 2 2 = Ipom ( 2 ) ;
Ipom = Y34 2 ∗ [ V2 ( 3 ) ; V2 ( 4 ) ] ; I 3 4 2 = Ipom ( 1 ) ;
I 4 3 2 = Ipom ( 2 ) ;
Ipom = Y45 2 ∗ [ V2 ( 4 ) ; V2 ( 5 ) ] ; I 4 5 2 = Ipom ( 1 ) ;
I 5 4 2 = Ipom ( 2 ) ;
I 4 0 2 = Y40sim ( 3 ) ∗V2 ( 4 ) ;

9
%n u l t a s e k v e n c a : z a h t i j e v a r a z b i j a n j e g r a f a na povezane
komponente
V1 0 = I g s i m ( 1 ) / ( Ygsim ( 1 ) + Y12 0 ( 1 , 1 ) ) ; V5 0 =
Ip5sim ( 1 ) / Y45 0 ( 2 , 2 ) ;
Y234 0 = [ Y23 0 ( 1 , 1 )+Y24 0 ( 1 , 1 ) Y23 0 ( 1 , 2 ) Y24 0 ( 1 ,
2 ) ; Y23 0 ( 2 , 1 ) Y23 0 ( 2 , 2 )+Y34 0 ( 1 , 1 ) Y34 0 ( 1 , 2 ) ;
Y24 0 ( 2 , 1 ) Y34 0 ( 2 , 1 ) Y40sim ( 1 )+Y24 0 ( 2 , 2 )+Y34 0 ( 2 ,
2) ] ;
I 2 3 4 0 = [ 0 ; −Ip3sim ( 1 ) ; 0 ] ; V234 0 = Y234 0 \ I 2 3 4 0 ;
Ipom = Y12 0 ∗ [ V1 0 ; V234 0 ( 1 ) ] ; I 1 2 0 = Ipom ( 1 ) ;
I 2 1 0 = Ipom ( 2 ) ;
Ipom = Y23 0 ∗ [ V234 0 ( 1 ) ; V234 0 ( 2 ) ] ; I 2 3 0 = Ipom ( 1 ) ;
I 3 2 0 = Ipom ( 2 ) ;
Ipom = Y24 0 ∗ [ V234 0 ( 1 ) ; V234 0 ( 3 ) ] ; I 2 4 0 = Ipom ( 1 ) ;
I 4 2 0 = Ipom ( 2 ) ;
Ipom = Y34 0 ∗ [ V234 0 ( 2 ) ; V234 0 ( 3 ) ] ; I 3 4 0 = Ipom ( 1 ) ;
I 4 3 0 = Ipom ( 2 ) ;
Ipom = Y45 0 ∗ [ V234 0 ( 3 ) ; V5 0 ] ; I 4 5 0 = Ipom ( 1 ) ;
I 5 4 0 = Ipom ( 2 ) ;

V1sim = [ V1 0 ; V1 ( 1 ) ; V2 ( 1 ) ] ; p o l a r P r i n t ( V1sim )
ans = 3 1 c e l l a r r a y { ’ 0 . 0 0 0 0 0 0 < 1 2 3 . 6 9 0 0 0 1 ’ } { ’
1 . 0 1 8 1 8 2 < −0.000012 ’ } { ’ 0 . 0 0 0 0 0 0 < 1 2 6 . 0 5 4 8 6 7 ’ }

V2sim = [ V234 0 ( 1 ) ; V1 ( 2 ) ; V2 ( 2 ) ] ; p o l a r P r i n t ( V2sim )

ans = 3 1 c e l l a r r a y { ’ 0 . 1 1 4 2 3 7 < 7 8 . 6 9 5 6 3 3 ’ } { ’
0 . 9 9 0 5 8 2 < −153.957799 ’ } { ’ 0 . 0 0 2 9 2 2 < −72.792802 ’ }
V3sim = [ V234 0 ( 2 ) ; V1 ( 3 ) ; V2 ( 3 ) ] ; p o l a r P r i n t ( V3sim )
ans = 3 1 c e l l a r r a y { ’ 0 . 1 1 4 2 3 0 < 7 8 . 6 9 7 4 3 4 ’ } { ’
0 . 9 7 7 5 3 4 < −154.554259 ’ } { ’ 0 . 0 0 3 5 2 4 < −78.193365 ’ }
V4sim = [ V234 0 ( 3 ) ; V1 ( 4 ) ; V2 ( 4 ) ] ; p o l a r P r i n t ( V4sim )
ans = 3 1 c e l l a r r a y { ’ 0 . 1 1 4 2 4 8 < 7 8 . 6 9 2 8 7 7 ’ } { ’
0 . 9 8 1 5 5 6 < −154.362548 ’ } { ’ 0 . 0 0 3 2 3 9 < −77.640569 ’ }
V5sim = [ V5 0 ; V1 ( 5 ) ; V2 ( 5 ) ] ; p o l a r P r i n t ( V5sim )
ans = 3 1 c e l l a r r a y { ’ 0 . 0 0 9 8 9 8 < 1 1 0 . 6 6 0 4 1 9 ’ } { ’
0 . 9 5 2 2 0 3 < −8.833235 ’ } { ’ 0 . 0 0 5 6 1 8 < 9 3 . 6 4 9 4 0 9 ’ }
I12sim = [ I12 0 ; I12 1 ; I12 2 ] ; pol arP rin t ( I12sim )
ans = 3 1 c e l l a r r a y { ’ 0 . 0 0 0 0 0 0 < 9 5 . 6 8 5 5 5 5 ’ } { ’

10
 0 . 2 0 7 4 6 3 < −23.543364 ’ } { ’ 0 . 0 0 8 1 1 8 < −132.635200 ’ }
I21sim = [ I21 0 ; I21 1 ; I21 2 ] ; pol arP rin t ( I21sim )
ans = 3 1 c e l l a r r a y { ’ 0 . 0 0 0 0 0 0 < 0 . 0 0 0 0 0 0 ’ } { ’
0 . 2 0 7 2 6 3 < 6 . 5 3 0 6 0 4 ’ } { ’ 0 . 0 0 8 1 1 8 < −162.634301 ’ }
I23sim = [ I23 0 ; I23 1 ; I23 2 ] ; pol arP rin t ( I23sim )
ans = 3 1 c e l l a r r a y { ’ 0 . 0 0 0 0 3 7 < −11.306917 ’ } { ’
0 . 1 5 1 0 4 6 < −173.377662 ’ } { ’ 0 . 0 0 6 1 3 1 < 2 0 . 4 2 5 4 2 0 ’ }
I32sim = [ I32 0 ; I32 1 ; I32 2 ] ; pol arP rin t ( I32sim )
ans = 3 1 c e l l a r r a y { ’ 0 . 0 0 0 0 4 1 < 1 6 8 . 6 9 3 4 7 4 ’ } { ’
0 . 1 5 1 0 6 6 < 6 . 6 0 0 2 5 6 ’ } { ’ 0 . 0 0 6 1 3 1 < −159.574377 ’ }
I24sim = [ I24 0 ; I24 1 ; I24 2 ] ; pol arP rin t ( I24sim )
ans = 3 1 c e l l a r r a y { ’ 0 . 0 0 0 0 3 7 < 1 6 8 . 6 9 3 0 8 3 ’ } { ’
0 . 0 5 6 2 1 7 < −173.715872 ’ } { ’ 0 . 0 0 2 0 2 3 < 8 . 0 5 4 2 1 2 ’ }
I42sim = [ I42 0 ; I42 1 ; I42 2 ] ; pol arP rin t ( I42sim )
ans = 3 1 c e l l a r r a y { ’ 0 . 0 0 0 0 2 8 < −11.307278 ’ } { ’
0 . 0 5 6 2 5 5 < 6 . 1 7 4 7 1 8 ’ } { ’ 0 . 0 0 2 0 2 2 < −171.946947 ’ }
I34sim = [ I34 0 ; I34 1 ; I34 2 ] ; pol arP rin t ( I34sim )
ans = 3 1 c e l l a r r a y { ’ 0 . 0 0 0 0 7 4 < 1 6 8 . 6 9 3 1 3 3 ’ } { ’
0 . 0 3 4 1 3 3 < 7 . 3 8 3 1 3 5 ’ } { ’ 0 . 0 0 1 8 8 5 < −141.721249 ’ }
I43sim = [ I43 0 ; I43 1 ; I43 2 ] ; pol arP rin t ( I43sim )
ans = 3 1 c e l l a r r a y { ’ 0 . 0 0 0 0 6 8 < −11.307059 ’ } { ’
0 . 0 3 4 1 0 7 < −172.481290 ’ } { ’ 0 . 0 0 1 8 8 5 < 3 8 . 2 7 4 8 1 4 ’ }
I45sim = [ I45 0 ; I45 1 ; I45 2 ] ; pol arP rin t ( I45sim )
ans = 3 1 c e l l a r r a y { ’ 0 . 0 0 0 0 0 0 < 0 . 0 0 0 0 0 0 ’ } { ’
0 . 0 2 2 4 8 7 < −177.842290 ’ } { ’ 0 . 0 0 1 0 2 6 < −59.561640 ’ }
I54sim = [ I54 0 ; I54 1 ; I54 2 ] ; pol arP rin t ( I54sim )
ans = 3 1 c e l l a r r a y { ’ 0 . 0 0 2 4 9 9 < 2 0 . 6 9 2 1 0 6 ’ } { ’
0 . 0 2 2 4 7 3 < 1 5 2 . 2 3 1 9 4 5 ’ } { ’ 0 . 0 0 1 0 2 6 < −29.555555 ’ }

11