Power System Stability

and Analysis

Assignment: 1

Yashwant Singh
101804044
3EE3
 Aesmunn

Epman MV VP

I)
SoovSAv 9:6 AVA s . 302 2 I1.713

2i2eV -2 v 4 0.08 2

New ba
AgSov.s v S 4o -2pv lokVA 8. 138

).2l2o 7.2bv Jobva

dasd /30 V 6kVA -4084ab) 120v -bvz_ -92 4
1-u

in8 . v .2v
133 p4 + -6-655pT PT
A M - AMw-loae
H Joscin Canlo n

V Yoo V
 :
S AsemmLnk at
TEEE S bus data

shunt
0 02t 0 . S-
A -2
B 3 0-03 ti 2y 2S-3-2 0.o 23S
C -3 0-6t 0.r O077-SD) 2
D -y P7-S
E -S .o ti 0-12 S-7S
3- .o ti 0.03 lo- 3
003 j:24 .32S
 Page

GO wiHhou Shb Eleen

O Wihat shunb Elemanb

Yh

63513.1 -5t 1s -25+i3-1

-SHlS 1o.01131.4 o17 tj313 67+StiS
-1-2s H3-75 -o.9174322-1s137.$3 lot3o o
-1-67S lo+i3 1-16j28 =L23j8-23
21-S 12 i3-s 3.7 2
 ) ih Shun men
Shunt olOmenh a arlae aong bA_DDPago

Ybu 62Sj1372s - S I 25 3s
-5 +jls 31-27-0.) h3732 -6) tjs -
12s +3-7 0.o74 3-722 13-1S773747-btj3o
-167 +jS lo tj3o R-933- .3ti31
s+i1- -35t3-7 3.1s-
i wim Shun lenuno

p-) 3

Y
-

-2

13 m
9-3
Y

-1 3

(3,3) Yhu 3)
Yu (3,1)=y bu ()-.
Yh
y buo (33) =Ybu (33) tja

Vb (2) -Ybus (h2)th Y tn

hu (a,2)= Yhuu (3,2)J -
yb (33)-Ybus (3,3 Jn
 Page

V 6S-j1b:lo-St+io-2s 9-7s O
-S tlo 10-7-31.2-o.914913-614S -2St 1s
F-28 i3-25-0.4)ti:13 12-S7-j147lot3
-7+jS--+3 91639-49-I. s 3-2
2 i 37S 3-1s6 j3/
)Sigan hrana lor mauHon Mahad
ihau Shun
S
AA
C

Y S
1:23-j37s

2S-j7S
lo-j3
1-2-i3

AT YA
 Page

Aly 5-s s-3-73 O

-Sis 04073712 I-iS 73 o

-12+37s 3-182 O

-615S -o30 125-3-3

-4115 1253-7

bus

-SS 2s3:7S O
-StL lo.37-j312 09a H313167 7
-1.25 3-7 -OGs7F3.182 2-1S37.S3 i3o
fo j3 RI33--12j3-
-3-54 2413s 3.%-i1-2
w; n shun
3

A L

C0

-l

O.0 S5
SSo
4-
 Page

Y=
9 25-1375
O .9)-3.73)

1:-S
1-13
23-7s 0

9.02
. o SS
00SSo

O O

Vaun A A
 Page
T 5-Is 1-3ti3:7o
0.0)j3-14 -61-isSo-5
-25+13-75 -14|3.78 0

-.5+15 -2s375

.0SSo

j.y

C.2s13 735 -S EiS -1-2s4j3-13

-stjs lo-0713122 -.0 31%
3-2 0.967 313 12-597-j71-

-2s7-
tj
h9-A 386 2 3:2
3-7sJ).2)
 wh Mh Tnpaoho
bw B anel

A
-

0 183-
o- 0.907-3-7T22) 0
Q O-

lo3

10-
O 100Sj»
jo-as
0.0
 Yhia
Ay-Ss 12-3- s
0.437-3-73 6IS_ S75

-Stj-S

j0.0SSo
aa
+1030
1S-1375
O 34j3.-7 jo
6s-137 -SJo
-o.91 i 73
-S+lo 1o-07-.32
o ) + 8-79 RIS7.47
j3o
-2Stj7S

- fo 3 >

191 23-9 - 3-7

- 1.2s 3-7S
 Page
K Alittanu hbiq

LO o3 jo.21

O2 %

.03 ojo2

0 2 . fa 2 jo 21
2(o) 0 24 0 )

Os2j3:6 -jo.1

o0011+j3-00
Zh
0 o.175 ho.%3
9 OD120 oy jo-13s
0-o7tjoos's 0.0392 2

Zh 0-0099 to-o 2
ti o-s3y 7SpoB
0-0 jo-o36 oYjo
oE.73
ol7s o-o
0O17sio-o53 o-0099 io-02 0.02S o.35
O 0 175 tja.S? s30.o75
t jo-

-.007 -j07 o.bo)

Zb Zu tio-
o.06-270-32-2
.3 B.1
a7-j021 a5377-ja
-a011-ja

017 0.9 o.O1nt.33) -9117 jo.o37 0-07 tjo-o

9.0 tjo.03 023) hj8:1 o a1i.99 9. ti
0-0127 t oo3 29jo.os
17t OOS)3 121i 3 .0S7jo

Zbua L o.613-ja.
co2 io2
.174jo.0s62 Doy3-0-314
-o2 tjo.072
O-0413-j./3Yf -9o4t3-jo.
 9

0169 4j0as 0.0123 E036 0-0174jp

-o123o3 o.p3 0 . 0 tjoo
O:O 3 i 34 9.0 o-0821 0 03 oo114 1SE0A
313 iS39 .013io13
Direct method

1.Without shunt

Code:

clc
clear
bus_system=5;
[line_data]=databus(bus_system);
nbr=size(line_data,1);
nbus=max(max(line_data(:,1)),max(line_data(:,2)));
z=line_data(:,3)+j*line_data(:,4);
y=1./z;
ysh=j*line_data(:,5);
ybus=zeros(nbus);
for k=1:nbr
m=line_data(k,1);
n=line_data(k,2);
ybus(m,n)=ybus(m,n)-y(k);
ybus(n,m)=ybus(m,n);
ybus(m,m)=ybus(m,m)+y(k)
ybus(n,n)=ybus(n,n)+y(k)
end
ybus
function [line_data]=databus(bus_system)
switch bus_system
case 5
line_data=[1 2 0.02 0.06 0;
1 3 0.08 0.24 0.025;
2 3 0.06 0.25 0.020;
2 4 0.06 0.18 0.020;
2 5 0.04 0.12 0.015;
3 4 0.01 0.03 0.010;
4 5 0.08 0.24 0.025]
end

Output:
2.With Shunt:

Code:

%direct method
clc
clear
bus_system=5;
[line_data]=databus(bus_system);
nbr=size(line_data,1);
nbus=max(max(line_data(:,1)),max(line_data(:,2)));
z=line_data(:,3)+j*line_data(:,4);
y=1./z;
ysh=j*line_data(:,5);
ybus=zeros(nbus);
for k=1:nbr
m=line_data(k,1);
n=line_data(k,2);
ybus(m,n)=ybus(m,n)-y(k);
ybus(n,m)=ybus(m,n);
ybus(m,m)=ybus(m,m)+y(k)+ysh(k);
ybus(n,n)=ybus(n,n)+y(k)+ysh(k);
end
ybus
function [line_data]=databus(bus_system)
switch bus_system
case 5
line_data=[1 2 0.02 0.06 0;
1 3 0.08 0.24 0.025;
2 3 0.06 0.25 0.020;
2 4 0.06 0.18 0.020;
2 5 0.04 0.12 0.015;
3 4 0.01 0.03 0.010;
4 5 0.08 0.24 0.025]
end

Output:
3.With Mutual Impedance

Code:
clc

clear

[line_data]=databus();
nbr=size(line_data,1);
nbus=max(max(line_data(:,1)),max(line_data(:,2)));
z=line_data(:,3)+j*line_data(:,4);
y=1./z;
ysh=j*line_data(:,5);
ybus=zeros(nbus);
l=0;
for k=1:nbr
if line_data(k,6)~=0
l=k;
end
m=line_data(k,1);
n=line_data(k,2);
ybus(m,n)=ybus(m,n)-y(k);
ybus(n,m)=ybus(m,n);
ybus(m,m)=ybus(m,m)+y(k)+ysh(k);
ybus(n,n)=ybus(n,n)+y(k)+ysh(k);
end
if l~=0
l=l-1;
p=line_data(l,1);
q=line_data(l,2);
r=line_data(l+1,1);
s=line_data(l+1,2);
ym1=-j*1/line_data(l,6);
ym2=-j*1/line_data(l+1,6);
h=[p q r s];
adsh=[1 -1 0 0;
0 0 1 -1;];
yb=adsh'*[y(l) ym1;
ym2 y(l+1);]*adsh;
for x=1:2
for y=3:4
ybus(h(x),h(y))=ybus(h(x),h(y))+yb(x,y);
ybus(h(y),h(x))=ybus(h(y),h(x))+yb(y,x) ;
end
end
end
ybus
function [line_data]=databus()
line_data=[1 2 0.02 0.06 0.025 0;
1 3 0.08 0.24 0.025 0.2;
2 3 0.06 0.25 0.020 0.2;
2 4 0.06 0.18 0.020 0;
2 5 0.04 0.12 0.015 0;
 3 4 0.01 0.03 0.010 0;
4 5 0.08 0.24 0.025 0;];

end
Output:

Singular Transformation Method

1.Without Shunt

Code:

%incidence matrix method

clc
clear
bus_system=5;
[line_data]=databus(bus_system);
nbr=size(line_data,1);
nbus=max(max(line_data(:,1)),max(line_data(:,2)));
z=line_data(:,3)+j*line_data(:,4);
y=1./z;
disp(y)
ysh=j*line_data(:,5);
ybus=zeros(nbus);
a=zeros(nbr,nbus);
ybus=zeros(nbus,nbus);
ypr=zeros(nbr,nbr);
zpr=zeros(nbr)
%ybus without shunt elements
for k=1:nbr
m=line_data(k,1);
n=line_data(k,2);
a(k,m)=1;
a(k,n)=-1;
ypr(k,k)=y(k);
zpr(k,k)=z(k);
end
ybus=a'*ypr*a
function [line_data]=databus(bus_system)
switch bus_system
case 5
line_data=[1 2 0.02 0.06 0;
1 3 0.08 0.24 0.025;
2 3 0.06 0.25 0.020;
2 4 0.06 0.18 0.020;
 2 5 0.04 0.12 0.015;
3 4 0.01 0.03 0.010;
4 5 0.08 0.24 0.025]

end
end
Output:

2.With Shunt

Code:

%incidence matrix method

clc
clear
bus_system=5;
[line_data]=databus(bus_system);
nbr=size(line_data,1);
nbus=max(max(line_data(:,1)),max(line_data(:,2)));
z=line_data(:,3)+j*line_data(:,4);
y=1./z;
yshc=zeros(nbus,nbus);
ysh=j*line_data(:,5);
for i=1:nbus
for z=1:nbr
if line_data(z,1)==i||line_data(z,2)==i
yshc(i,i)=yshc(i,i)+ysh(z);
end
end
end
ybus=zeros(nbus);
a=zeros(nbr+nbus,nbus);
ybus=zeros(nbus,nbus);
ypr=zeros(nbr+nbus,nbr+nbus);
ym=0;
k=0;
for k=1:nbr+nbus
if k<=nbr
m=line_data(k,1);
n=line_data(k,2);
a(k,m)=1;
a(k,n)=-1;
ypr(k,k)=y(k);
else
a(k,k-nbr)=1;
ypr(k,k)=yshc(k-nbr,k-nbr);
end
end
ybus=a'*ypr*a

function [line_data]=databus(bus_system)
switch bus_system
case 5
line_data=[1 2 0.02 0.06 0 0;
1 3 0.08 0.24 0.025 0.2;
2 3 0.06 0.25 0.020 0.2;
2 4 0.06 0.18 0.020 0;
2 5 0.04 0.12 0.015 0;
3 4 0.01 0.03 0.010 0;
4 5 0.08 0.24 0.025 0];

end
end

Output:

3.With Mutual Impedance

Code:

%incidence matrix method

clc
clear
bus_system=5;
[line_data]=databus(bus_system);
nbr=size(line_data,1);
nbus=max(max(line_data(:,1)),max(line_data(:,2)));
z=line_data(:,3)+j*line_data(:,4);
y=1./z;
yshc=zeros(nbus,nbus);
ysh=j*line_data(:,5);
for i=1:nbus
for z=1:nbr
if line_data(z,1)==i||line_data(z,2)==i
yshc(i,i)=yshc(i,i)+ysh(z);
end
end
end
ybus=zeros(nbus);
a=zeros(nbr+nbus,nbus);
ybus=zeros(nbus,nbus);
ypr=zeros(nbr+nbus,nbr+nbus);
ym=0;
k=0;
for k=1:nbr+nbus
if k<=nbr
m=line_data(k,1);
n=line_data(k,2);
a(k,m)=1;
a(k,n)=-1;
ypr(k,k)=y(k);

if line_data(k,6)~=0
flag=k-1;
end
else
a(k,k-nbr)=1;
ypr(k,k)=yshc(k-nbr,k-nbr);
end
end
ym=1/j*inv(line_data(flag,6));
ypr(flag,flag+1)=ypr(flag,flag+1)+ym;
ypr(flag+1,flag)=ypr(flag+1,flag)+ym;
line_data()
ybus=a'*ypr*a

function [line_data]=databus(bus_system)
switch bus_system
case 5
line_data=[1 2 0.02 0.06 0 0;
1 3 0.08 0.24 0.025 0.2;
2 3 0.06 0.25 0.020 0.2;
2 4 0.06 0.18 0.020 0;
2 5 0.04 0.12 0.015 0;
3 4 0.01 0.03 0.010 0;
4 5 0.08 0.24 0.025 0];

end
end
Output:
Zbus Formulation

Code:

clc
clear
[line_data]=databus()
nbr=size(line_data,1);
zbus=[0];
currentbusno=0;
for i=1:7
from=line_data(i,2);
to=line_data(i,3);
value=line_data(i,4)+j*line_data(i,5);
dim=size(zbus,1);
newbus=max(from,to);
ref=min(from,to);
if newbus>currentbusno && ref==1
zbus=[zbus zeros(dim,1);
zeros(1,dim) value;]
currentbusno=newbus;
continue
end
if newbus>currentbusno && ref~=1
zbus=[zbus zbus(:,ref);
zbus(ref,:) zbus(ref,ref)+value;]
currentbusno=newbus;
continue
end
if newbus<=currentbusno && ref==1
x=zbus(:,newbus)
y=zbus(newbus,:)
z=(zbus(newbus,newbus)+value)
z=1/z
zbus=zbus-x*z*y
continue
end
if newbus<=currentbusno && ref~=1
x=zbus(:,ref)-zbus(:,newbus)
y=zbus(ref,:)-zbus(newbus,:)
z=(zbus(ref,ref)+zbus(newbus,newbus)-(2*zbus(ref,newbus))+value);
z=1/z
zbus=zbus-x*z*y
continue
end
end
Output: