Scribd
Upload
27 views

Compact Short

The document contains code for calculating the initial state of a power system using a Newton-Raphson load flow method. It defines line impedances and admittances, builds an admittance matrix, specifies an initial solution, calculates mismatches, computes corrections using the Jacobian, and iterates until mismatches are small. The final section calculates power flows between buses.

Uploaded by

Jser Misael Csar
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
27 views

Compact Short

The document contains code for calculating the initial state of a power system using a Newton-Raphson load flow method. It defines line impedances and admittances, builds an admittance matrix, specifies an initial solution, calculates mismatches, computes corrections using the Jacobian, and iterates until mismatches are small. The final section calculates power flows between buses.

Uploaded by

Jser Misael Csar
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 4

clc,clear,format compact ,format short

Sb=400;
Vb=100;
zb=(Vb^2)/Sb
disp('impedancias pu de las lineas')
z12=(0.0066+0.0321*i)
z13=(0.0052+0.0257*i)
z15=(0.0035+0.0171*i)
z24=(0.0061+0.0300*i)
z34=(0.0087+0.0428*i)
z45=(0.0035+0.0171*i)

disp('admitancias pu de las lineas')


Y12=1/z12
Y13=1/z13;
Y15=1/z15;
Y24=1/z24;
Y34=1/z34;
Y45=1/z45;
Y12SH=i*0.4007
Y13SH=i*0.3278
Y15SH=i*0.2185
Y24SH=i*0.3824
Y34SH=i*0.5463
Y45SH=i*0.2185

disp('Ybarra cartesiana')
ybarra=[Y12+Y13+Y15+Y12SH+Y13SH+Y15SH -Y12 -Y13 0 -Y15;-Y12
Y12+Y24+Y12SH+Y24SH 0 -Y24 0;-Y13 0 Y13+Y34+Y13SH+Y34SH -Y34 0;0 -Y24 -
Y34 Y24+Y34+Y45+Y24SH+Y34SH+Y45SH -Y45;-Y15 0 0 -Y45 Y15+Y45+Y15SH+Y45SH]
Y=ybarra;
disp('Ybarra polar')
YP=[abs(Y(1,1)) angle(Y(1,1))*180/pi abs(Y(1,2)) angle(Y(1,2))*180/pi
abs(Y(1,3)) angle(Y(1,3))*180/pi abs(Y(1,4))]
disp('SOLUSION INICIAL')
v1=1.03;d1=0;
v4=1.02;d4i=0;
v2i=1;d2i=0;
v3i=1;d3i=0;
v5i=1;d5i=0;
Sb=100;
Vb=400;
Sesp2=-2.4-1.96*j
P2in=real(Sesp2)
Q2in=imag(Sesp2)
Sesp3=-1.9-1.25*j
P3in=real(Sesp3)
Q3in=imag(Sesp3)
Sesp4=2.83
Sesp5=1.24-0.72*j
P5in=real(Sesp5)
Q5in=imag(Sesp5)
iteracion=input('ingresa tu iteracion : ')
for n=1:iteracion

syms v2 d2 v3 d3 d4 v5 d5
p2=v2*v1*abs(ybarra(2,1))*cos(angle(ybarra(2,1))-
d2+d1)+v2*v2*abs(ybarra(2,2))*cos(angle(ybarra(2,2)))+v2*v4*abs(ybarra(2,
4))*cos(angle(ybarra(2,4))-d2+d4);
p3=v3*v1*abs(ybarra(3,1))*cos(angle(ybarra(3,1))-
d3+d1)+v3*v4*abs(ybarra(3,4))*cos(angle(ybarra(3,4))-
d3+d4)+v3*v3*abs(ybarra(3,3))*cos(angle(ybarra(3,3)));
p4=v4*v3*abs(ybarra(4,3))*cos(angle(ybarra(4,3))-
d4+d3)+v4*v5*abs(ybarra(4,5))*cos(angle(ybarra(4,5))-
d4+d5)+v4*v4*abs(ybarra(4,4))*cos(angle(ybarra(4,4)))+v4*v2*abs(ybarra(4,
2))*cos(angle(ybarra(4,2))-d4+d2);
p5=v5*v1*abs(ybarra(5,1))*cos(angle(ybarra(5,1))-
d5+d1)+v5*v4*abs(ybarra(5,4))*cos(angle(ybarra(5,4))-
d5+d4)+v5*v5*abs(ybarra(5,5))*cos(angle(ybarra(5,5)));

q2=-(v2*v1*abs(ybarra(2,1))*sin(angle(ybarra(2,1))-
d2+d1)+v2*v2*abs(ybarra(2,2))*sin(angle(ybarra(2,2)))+v2*v4*abs(ybarra(2,
4))*sin(angle(ybarra(2,4))-d2+d4));
q3=-(v3*v1*abs(ybarra(3,1))*sin(angle(ybarra(3,1))-
d3+d1)+v3*v4*abs(ybarra(3,4))*sin(angle(ybarra(3,4))-
d3+d4)+v3*v3*abs(ybarra(3,3))*sin(angle(ybarra(3,3))));
q5=-(v5*v1*abs(ybarra(5,1))*sin(angle(ybarra(5,1))-
d5+d1)+v5*v4*abs(ybarra(5,4))*sin(angle(ybarra(5,4))-
d5+d4)+v5*v5*abs(ybarra(5,5))*sin(angle(ybarra(5,5))));

F=[p2 p3 p4 p5 q2 q3 q5];incog=[d2 d3 d4 d5 v2 v3 v5];


j=jacobian(F,incog);
d2=d2i;
d3=d3i;
d4=d4i;
d5=d5i;
v2=v2i;
v3=v3i;
v5=v5i;
p2e=eval(p2);
p3e=eval(p3);
p4e=eval(p4);
p5e=eval(p5);
q2e=eval(q2);
q3e=eval(q3);
q5e=eval(q5);
DeVe=[d2;d3;d4;d5;v2;v3;v5];

dp2=P2in-p2e;
dp3=P3in-p3e;
dp4=Sesp4-p4e;
dp5=P5in-p5e;
dq2=Q2in-q2e;
dq3=Q3in-q3e;
dq5=Q5in-q5e;
if abs(dp2)<=0.00001 & abs(dp3)<=0.00001 & abs(dp4)<=0.00001 &
abs(dp5)<=0.00001 & abs(dq2)<=0.00001 & abs(dq3)<=0.00001 &...
abs(dq5)<=0.00001;
break
end
res=[dp2;dp3;dp4;dp5;dq2;dq3;dq5];
je=eval(j);
D=je\res
D1V1=DeVe+D;
d2i=D1V1(1);
d3i=D1V1(2);
d4i=D1V1(3);
d5i=D1V1(4);
v2i=D1V1(5);
v3i=D1V1(6);
v5i=D1V1(7);
salida(n,:)=[n v2 v3 v5 d2*180/pi d3*180/pi d4*180/pi d5*180/pi];
end
disp(' iter v2 v3 v5 d2 d3 d4
d5')
salida
V1=v1
V2=v2*(cos(d2)+i*sin(d2))
V3=v3*(cos(d3)+i*sin(d3))
V4=v4*(cos(d4)+i*sin(d4))
V5=v5*(cos(d5)+i*sin(d5))
y1i=i*(0.4007+0.3278+0.2185);
y2i=i*(0.4007+0.3824);
y3i=i*(0.3278+0.5463);
y4i=i*(0.3824+0.5463+0.2185);
y5i=i*(0.2185+0.2185);
i12=-(V1-V2)*ybarra(1,2)+y1i*V1;
i21=-(V2-V1)*ybarra(1,2)+y2i*V2;
s12=V1*conj(i12)
s21=V2*conj(i21)
SL12=s12+s21
i13=-(V1-V3)*ybarra(1,3)+y1i*V1;
i31=-(V3-V1)*ybarra(1,3)+y3i*V3;
s13=V1*conj(i13)
s31=V3*conj(i31)
SL13=s13+s31
i15=-(V1-V5)*ybarra(1,5)+y1i*V1;
i51=-(V5-V1)*ybarra(1,5)+y5i*V5;
s15=V1*conj(i15)
s51=V5*conj(i51)
SL15=s15+s51
i34=-(V3-V4)*ybarra(4,3)+y3i*V3;
i43=-(V4-V3)*ybarra(4,3)+y4i*V4;
s34=V3*conj(i34)
s43=V4*conj(i43)
SL34=s34+s43
i24=-(V2-V4)*ybarra(2,4)+y2i*V2;
i42=-(V4-V2)*ybarra(4,2)+y4i*V4;
s24=V2*conj(i24)
s42=V4*conj(i42)
SL24=s24+s42

You might also like