ttcs_tia.m

solve power flows by gauss seidel method ans newton Rapson

function ttcs_tia(P1,Q1,R,X,B,U)
S1=P1+j*Q1
S2=S1-j*(B*10^(-3)/2)*(U^2)
Delta_U=(real(S2)*R+imag(S2)*X)/U
Phan_tram_U=Delta_U*100/U
Delta_P=(real(S2)^2+imag(S2)^2)*R/U^2
Delta_Q=(real(S2)^2+imag(S2)^2)*X/U^2
S3=S2+Delta_P+j*Delta_Q
S4=S3-j*(B*10^(-3)/2)*(U^2)