ttcs_vong.m

solve power flows by gauss seidel method ans newton Rapson

function ttcs_vong(R1,R2,R3,X1,X2,X3,B1,B2,B3,P1,P2,Q1,Q2,U)
Z=R1+R2+R3+j*(X1+X2+X3)
Y=1/Z
G=real(Z)/((real(Z)^2)+(imag(Z)^2))
B=imag(Z)/((real(Z)^2)+(imag(Z)^2))
R1B=R2+R3
X1B=X2+X3
R2B=R3
X2B=X3
S_1phay=P1+j*Q1-j*(B1*10^(-3)/2)*U^2-j*(B2*10^(-3)/2)*U^2
S_2phay=P2+j*Q2-j*(B2*10^(-3)/2)*U^2-j*(B3*10^(-3)/2)*U^2
M=(real(S_1phay)*R1B+imag(S_1phay)*X1B)+(real(S_2phay)*R2B+imag(S_2phay)*X2B)
N=(real(S_1phay)*X1B-imag(S_1phay)*R1B)+(real(S_2phay)*X2B-imag(S_2phay)*R2B)
PA1=G*M+B*N
QA1=-G*N+B*M
R1A=R1
X1A=X1
R2A=R1+R2
X2A=X1+X2
M_phay=(real(S_1phay)*R1A+imag(S_1phay)*X1A)+(real(S_2phay)*R2A+imag(S_2phay)*X2A)
N_phay=(real(S_1phay)*X1A-imag(S_1phay)*R1A)+(real(S_2phay)*X2A-imag(S_2phay)*R2A)
PB2=G*M_phay+B*N_phay
QB2=-G*N_phay+B*M_phay
SA1=PA1+j*QA1
SB2=PB2+j*QB2