1. c) Computation of Parameters and Modelling of Transmission Lines in MATLAB | PSS LAB | EEE

Modelling and performance analysis of transmission line ABCD parameters in MATLAB PROBLEM 1(c): A 3 phase overhead line 200km long, R=0.16Ω and conductor diameter of 2cm with spacing 4m,5m and 6m transposed. Find A,B,C,D constants, sending end voltage and current, power factor, power when line is developing full load of 50MW at 132KV,0.8 p.f  lagging, transmission efficiency, receiving end voltage and regulation. Theoretical Calculation:                                                                                PROGRAM 1(c): clear; clear all; l=200000; PR=50000000; Pfr=0.8; Vrl=132000; R=32; D=0.02; r=D/2; Dab=4; Dbc=5; Dca=6; DM=(Dab*Dbc*Dca)^(1/3); GMR=r*exp(-1/4); L=(2*10^(-7)*(log(DM)-log(GMR))); c=(2*pi*8.854*10^(-12))/(log(DM)-log(r)); X=2*pi*50*L*l; Z=R+1j*X; Y=1j*(2*pi*50*c*l); Z1=abs(Z); Z2=angle(Z); A=1+((Y*Z)/2); D=1+((Y*Z)/2); B=Z; C=Y*(1+(Y*Z)/4); Vr=Vrl/(3^(1/2)); Ir1=(PR/(3*Vr*Pfr)); Ir2=sin(acos(Pfr)); Ir=Ir1*(Pfr-1j*Ir2); Vs=(A*Vr)+(B*Ir); Vs1=abs(Vs); Vs2=angle(Vs); Is=(C*Vr)+(D*Ir); Is1=abs(Is); Is2=angle(Is); e=Vs2-Is2; T=cos(e); Ps=3*Vs*Is*T; Ps1=abs(Ps); U=(PR/Ps1)*100; A1=abs(A); A2=angle(A); reg=(((Vs1/A1)-Vr)/Vr)*100; Ind = ['Value of Inductance               :' num2str(L)]; disp(Ind); Cap = ['Value of Capacitance              :' num2str(c)]; disp(Cap); z = ['Value of Impedance                :' num2str(Z)]; disp(z); y  = ['Value of Admittance               :' num2str(Y)]; disp(y); a = ['Value of Constant A               : ' num2str(A)]; disp(a); b = ['Value of Constant B               : ' num2str(B)]; disp(b); c1 = ['Value of Constant C               :' num2str(C)]; disp(c1); d = ['Value of Constant D               :' num2str(D)]; disp(d); vr = ['Value of Receving end Voltage Vr  :' num2str(Vr)]; disp(vr); ir = ['Value of Receving end Current Ir  :' num2str(Ir)]; disp(ir); vs = ['Value of Sending end Voltage Vr   :' num2str(Vs)]; disp(vs); is = ['Value of Sending end Current Is   :' num2str(Is)]; disp(is); u = ['Value of Eficency                 :' num2str(U)]; disp(u); Re = ['Value of %Regulation:             :' num2str(reg)]; disp(Re); t = ['Value of Sending end P.f          :' num2str(T)]; disp(t); ps1 = ['Value of Sending end Power Ps     :' num2str(Ps1)]; disp(ps1);                                         MATLAB Output:    YouTube Video of My Channel on the Problem you can use this also ( If you have any problem in doing this you can comment below  , I will help as soon as possible) 💗  Thanks For Visiting  💗