diff options
Diffstat (limited to '2441/CH6/EX6.3/Ex6_3.sce')
-rwxr-xr-x | 2441/CH6/EX6.3/Ex6_3.sce | 73 |
1 files changed, 73 insertions, 0 deletions
diff --git a/2441/CH6/EX6.3/Ex6_3.sce b/2441/CH6/EX6.3/Ex6_3.sce new file mode 100755 index 000000000..ae154bdda --- /dev/null +++ b/2441/CH6/EX6.3/Ex6_3.sce @@ -0,0 +1,73 @@ +//exa 6.3
+clc;clear;close;
+format('v',7);
+l=350;//km(length of line)
+Z=180*expm(%i*75*%pi/180);//ohm/phase(Total)
+Y=1*10^-3*expm(%i*90*%pi/180);//Siemens/phase(Total)
+z=Z/l;//ohm/km
+y=Y/l;//Siemens/km
+re=l*sqrt(z*y);//
+Zc=sqrt(z/y);//ohm
+disp("Part(a) A,B,C,D parameters are : ");
+A=cosh(re);//unitless
+D=A;//unitless
+B=Zc*sinh(re);//ohm
+C=sinh(re)/Zc;//unitless
+A_mag=abs(A);//unitless
+A_angle=atand(imag(A)/real(A));//degree
+B_mag=abs(B);//ohm
+B_angle=atand(imag(B)/real(B));//degree
+C_mag=abs(C);//unitless
+C_angle=atand(imag(C)/real(C));//degree
+C_angle=C_angle+180;//degree(Converting -ve to +ve angle)
+D_mag=abs(D);//unitless
+D_angle=atand(imag(D)/real(D));//degree
+disp(A_mag,"Magnitude of A : ");
+format('v',5);
+disp(A_angle,"Angle of A(degree) : ");
+format('v',7);
+disp(B_mag,"Magnitude of B(ohm) : ");
+format('v',6);
+disp(B_angle,"Angle of B(degree) : ");
+format('v',8);
+disp(C_mag,"Magnitude of C : ");
+format('v',6);
+disp(C_angle,"Angle of C(degree) : ");
+format('v',7);
+disp(D_mag,"Magnitude of D : ");
+format('v',5);
+disp(D_angle,"Angle of D(degree) : ");
+//60% series compensation
+B=B-%i*60/100*abs(Z)*sind(atand(imag(Z),real(Z)));//ohm(considering series compensation=60%)
+//For Equivalent pi-circuit
+disp("Part(b) A,B,C,D parameters of compensated line are : ");
+Ydash=2/Zc*[(cosh(re)-1)/sinh(re)];//S
+A=1+B*Ydash/2;//unitless
+D=A;//unitless
+C=2*Ydash/2+B*(Ydash/2)^2;//unitless
+A_mag=abs(A);//unitless
+A_angle=atand(imag(A)/real(A));//degree
+B_mag=abs(B);//ohm
+B_angle=atand(imag(B)/real(B));//degree
+C_mag=abs(C);//unitless
+C_angle=atand(imag(C)/real(C));//degree
+C_angle=C_angle+180;//degree(Converting -ve to +ve angle)
+D_mag=abs(D);//unitless
+D_angle=atand(imag(D)/real(D));//degree
+format('v',4);
+disp(B_mag,"Magnitude of B(ohm) : ");
+format('v',6);
+disp(B_angle,"Angle of B(degree) : ");
+format('v',7);
+disp(A_mag,"Magnitude of A : ");
+format('v',5);
+disp(A_angle,"Angle of A(degree) : ");
+format('v',6);
+disp(C_mag,"Magnitude of C : ");
+format('v',5);
+disp(C_angle,"Angle of C(degree) : ");
+format('v',7);
+disp(D_mag,"Magnitude of D : ");
+format('v',5);
+disp(D_angle,"Angle of D(degree) : ");
+//Answer for some parts are not accurate in the textbook.
|