From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 1076/CH3/EX3.18/3_18.sce | 50 ++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 50 insertions(+)
 create mode 100755 1076/CH3/EX3.18/3_18.sce

(limited to '1076/CH3/EX3.18/3_18.sce')

diff --git a/1076/CH3/EX3.18/3_18.sce b/1076/CH3/EX3.18/3_18.sce
new file mode 100755
index 000000000..9d3a2973c
--- /dev/null
+++ b/1076/CH3/EX3.18/3_18.sce
@@ -0,0 +1,50 @@
+clear;
+clc;
+
+A1=complex(.98 * cosd(2), .98* sind(2));
+B1=complex(28 * cosd(69), 28* sind(69));
+C1=complex(.0002 * cosd(80), .0002* sind(80));
+D1=A1;
+
+A2=complex(.95 * cosd(3), .95* sind(3));
+B2=complex(40 * cosd(85), 40* sind(85));
+C2=complex(.0004 * cosd(90), .0004* sind(90));
+D2=A2;
+
+//(a)
+
+A= (A1* A2) + (B1* C2);
+B= (A1* B2) + (B1* D2);
+C= (C1* A2) + (D1* C2);
+D= (C1* B2) + (D1* D2);
+
+
+mprintf("(a)");
+mprintf("\nA= %s", string(round(abs(A)*1000)/1000) +'/_'+ string(round(atand(imag(A)/real(A))*10)/10) )
+mprintf("\nB= %s", string(round(abs(B)*100)/100) +'/_'+ string(round(atand(imag(B)/real(B))*100)/100) )
+mprintf("\nC= %s *1e-4", string(round(abs(C)*100000)/10) +'/_'+ string(round(atand(imag(C)/real(C))*10)/10 ) )
+mprintf("\nD= %s", string(round(abs(D)*1000)/1000) +'/_'+ string(round(atand(imag(D)/real(D))*10)/10) )
+
+
+//(b)
+
+Vr=110e3/sqrt(3);
+pf=.95;
+Irm=200
+pfa=-1* acos(pf);
+Ir=complex(Irm *pf, Irm * sin(pfa));
+
+Vs=(A*Vr)+(B*Ir);
+V=abs(Vs)
+vs=V*sqrt(3)*1e-3;
+mprintf("\n\n\n(b)\nSending End Voltage= %.2f kV", vs)
+
+Is=(C*Vr)+(D*Ir);
+I=abs(Is)
+mprintf("\nSending End Current= %.1f A", I)
+
+phi1=atan(imag(Vs)/real(Vs))*(180/%pi);
+phi2=atan(imag(Is)/real(Is))*(180/%pi);
+phi=phi1-phi2;
+pfs=cosd(phi);
+mprintf("\nSending End pf= %.2f ",pfs)
-- 
cgit