Rigorous Method
R=input('Resistance :');
XL=input('Inductive Reactance :');
XC=input('Capacitive Reactance :');
G=input('Conductance :');
length=input('Length of Transmission Line :');
f=input('Frequency :');
z= (R+j*XL)*length;
y=(G+j*XC)*length;
gm=sqrt(z*y);
zc=sqrt(y/z);
A=cosh(gm);
B=zc*sinh(gm);
C=(1/zc)*sinh(gm);
D=A;
TM=[A B;C D]; Z=B;
Y=(2/zc)*tanh(gm/2);
VRL=input('ENTER RECEIVEING END VOLTAGE :');
VRP=VRL*10^(3)/(sqrt(3));
PR = input('ENTER RECEIVING END LOAD IN MW :');
%%PR = PR*10^(6);
Pf=input('ENTER THE RECEIVING END LOAD POWER FACTOR :');
IR=PR*10^(6)/(sqrt(3)*VRL*10^(3)*Pf);
VS=VRP*A+IR*B;
IS1=VRP*B+IR*A;
%%h=acos(Pf);
%%SR=PR/Pf;
%%SR=SR*(cos(h)+j*sin(h));
%%QR=imag(SR);
%%SM=TM*[VRP;IR];
%%VS=SM(1,1);
%%IS=SM(2,1);
Pfs=cos(angle(VS)-angle(IS1));
%%VSA=angle(VS)*(180/pi);
%%ISA=angle(IS)*(180/pi);
VS=sqrt(3)*VS/1000;
IS=abs(IS1)/1000;
SS=3*VS*IS;
VREG=(((VS*sqrt(3)/A)-VRL)/VRL)*100;
PS=real(SS);
QS=imag(SS);
eff=(PR/PS)*100;
PL=PS-PR;
QL=QS-QR;
z
y
zc
Z
Y
TM
fprintf('receiving END current %g at %g degrees \n',IR);
fprintf('SENDING END VOLTAGE %g at %g degrees \n',VS,VSA);
fprintf('SENDING END CURRENT %g at %g degrees \n',IS1,ISA);
fprintf('SENDING END POWER FACTOR %g\n',Pfs);
fprintf('SENDING END REAL POWER %g\n',PS);
fprintf('SENDING END REACTIVE POWER %g\n',QS);
fprintf('PERCENTAGE VOLTAGE REGULATION %g\n',VREG);
fprintf('REAL POWER LOSS %g\n',PL);
fprintf('REACTIVE POWER LOSS %g\n',QL);
fprintf('EFFICIENCY %G', eff);
ABCD constant T method
R=input('Resistance :');
XL=input('Inductive Reactance :');
XC=input('Capacitive Reactance :');
G=input('Conductance :');
length=input('Length of Transmission Line :');
f=input('Frequency :');
z= (R+j*XL)*length;
y=(G+j*XC)*length;
gm=sqrt(z*y)/2;
zc=sqrt(y/z);
A=1+gm;
B=z*(1+(gm/2));
C=y;
D=A;
f=input('Frequency :');
VRL=input('ENTER RECEIVEING END VOLTAGE :');
VRP=VRL*10^(3)/(sqrt(3));
PR = input('ENTER RECEIVING END LOAD IN MW :');
%%PR = PR*10^(6);
Pf=input('ENTER THE RECEIVING END LOAD POWER FACTOR :');
IR=PR*10^(6)/(sqrt(3)*VRL*10^(3)*Pf);
VS=VRP*A+IR*B;
IS1=VRP*B+IR*A;
VS=sqrt(3)*abs(VS)/1000;
IS=abs(IS)*1000;
VREG=((VS/(abs(TM(1,1)))-VRL)/VRL)*100;
PS=real(SS);
QS=imag(SS);
eff=PR/PS*100;
PL=PS-PR;
QL=QS-QR;
A
B
C
D
fprintf('receiving END current %g at %g degrees \n',IR);
fprintf('SENDING END VOLTAGE %g at %g degrees \n',VS,VSA);
fprintf('SENDING END CURRENT %g at %g degrees \n',IS,ISA);
fprintf('SENDING END POWER FACTOR %g\n',Pfs);
fprintf('SENDING END REAL POWER %g\n',PS);
fprintf('SENDING END REACTIVE POWER %g\n',QS);
fprintf('PERCENTAGE VOLTAGE REGULATION %g\n',VREG);
fprintf('REAL POWER LOSS %g\n',PL);
fprintf('REACTIVE POWER LOSS %g\n',QL);
fprintf('EFFICIENCY %G', eff);