clc

//Variables de entrada//
open=input('presiona enter')

S=100                             //Potencia aparente del sistema//
V=110                             //Voltaje de la línea//
f=50                              //Frecuencia del sistema//
l=400                             //Largo de la línea//
cir=2                             //Número de circuitos//
seccion=765                       //Sección del conductor en mm2//
n=4                               //Número de subconductores//
G=0                               //Conductancia de la línea//
d=0.45                            //Distancia entre conductores//
D12=4.70744                        //Distancia entre conductores//
D23=4.70744                       //Distancia entre conductores//
D13=9.2                       //Distancia entre conductores//
RL=0.010172                        //Resistencia de la línea por subconductor//
fp=0.95
         //Corriente en pu del consumo//
r=(sqrt(seccion/%pi))*(0.001)     //Radio del conductor//
W=(2*%pi*f)                       //Omega//
E0=8.8542*10^(-12)                //Epsilon en vacío//

//Condición dependiente del n //
if n==1 then reql=(%e^-0.25)*r
   elseif n==2 then reql=sqrt((%e^-0.25)*r*d)
   elseif n==3 then reql=((%e^-0.25)*r*(d^2))^(1/3)
   elseif n==4 then reql=1.09*((((%e)^-0.25)*r*(d^3))^(1/4))
   else reql=r
end

if n==1 then reqc=r
elseif n==2 then reqc=sqrt(r*d)
elseif n==3 then reqc= (r*(d^2))^(1/3)
elseif n==4 then reqc=1.09*((r*(d^3))^(1/4))
else reqc=r
end

Deq=(D12*D13*D23)^(1/3)            //Distancia media entre conductores//
Zb=(V*V)/S                         //Cálculo de la impedancia base de la línea//

//Ecuaciones de inductancia//

L=(2*10^(-7))*(log(Deq/reql))*1000 //Inductancia en H/KM//
XL=W*L                             //Reactancia en ohm/KM//
aproxZ=(complex(RL,XL))*l                   //Impedancia de la línea//

//Ecuaciones de conductancia//

c=((2*%pi*E0)/(log(Deq/reqc)))*1000//Capacitancia en F/KM//
aproxY=W*c*l                              //Admitancia en S/KM//
aproxYc=complex(G,aproxY)                    //Conductancia de la línea//

//Ecuaciones para la matriz//

A2=((aproxZ*aproxYc)/2)+1
B2=(aproxZ/Zb)/2
C2=aproxYc*(1+((aproxZ*aproxYc)/4))*Zb*2
D2=A2
aproxT=[A2,B2;C2,D2]
aproxF=1/(sqrt(real(A2)^2+(imag(A2)^2)))    //Efecto ferranti en la línea//
fprad=acos(fp)

//Cálculo de los parámetros del generador//

I2=complex(1*cos(fprad),1*sin(-fprad))
V2=1
V1=(A2*V2)+(B2*I2)
I1=(C2*V2)+(D2*I2)
S=V1*I1

//Cálculo de la potencia del generador//

phi1=(atan((imag(V1))/(real(V1))))  //Ángulo phi del voltaje 1//
phi2=(atan((imag(I1))/(real(I1))))  //Ángulo phi de la corriente 1//

phi12=phi1-phi2                     //Diferencia de ángulos entre V1 e I1//
fp1=cos(phi12)                      //Factor de potencia del generador//

I1r=real(I1)                        //Valor real de la corriente 1//
I1i=imag(I1)                        //Valor imaginario de la corriente 1//

P1=real(V1*complex(I1r,-I1i))

//Cálculos de la potencia de salida//

phi3=atan(imag(V2)/real(V2))       //Ángulo phi del voltaje 2//
phi4=atan(imag(I2)/real(I2))       //Ángulo phi de la corriente 2//
phi34=phi3-phi4                    //Diferencia de ángulos entre V2 e I2//
fp2=cos(phi34)                     //Factor de potencia del consumo//
I2r=real(I2)                       //Valor real de la corriente 2//
I2i=imag(I2)                       //Valor imaginario de la corriente //
P2=real(V2*complex(I2r,-I2i))      //Potencia real del consumo//

aproxRend=P2/P1                         //Rendimiento del sistema//

aproxReg=(((sqrt(imag(V1)^2+real(V1)^2))-(sqrt(imag(V2)^2+real(V2)^2)))/(sqrt((imag(V1)^2)+(real(V1))^2)))*100

disp(aproxRend,'El Rendimiento del Sistema es:')
disp(aproxReg,'La Regulación es:')
disp(aproxF,'El efecto Ferranti es:')
disp(aproxT,'La matriz de la línea es:'