From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Tue, 10 Oct 2017 12:27:19 +0530
Subject: initial commit / add all books

---
 147/CH3/EX3.13/Example3_13.sce | 17 +++++++++++++++++
 147/CH3/EX3.13/Result3_13.txt  |  4 ++++
 147/CH3/EX3.14/Example3_14.sce | 17 +++++++++++++++++
 147/CH3/EX3.14/Result3_14.txt  |  3 +++
 147/CH3/EX3.18/Example3_18.sce | 15 +++++++++++++++
 147/CH3/EX3.18/Result3_18.txt  |  4 ++++
 147/CH3/EX3.19/Example3_19.sce | 24 ++++++++++++++++++++++++
 147/CH3/EX3.19/Result3_19.txt  |  7 +++++++
 147/CH3/EX3.2/Example3_2.sce   | 12 ++++++++++++
 147/CH3/EX3.2/Result3_2.txt    |  3 +++
 147/CH3/EX3.3/Example3_3.sce   | 26 ++++++++++++++++++++++++++
 147/CH3/EX3.3/Result3_3.txt    |  1 +
 147/CH3/EX3.4/Example3_4.sce   | 12 ++++++++++++
 147/CH3/EX3.4/Result3_4.txt    |  1 +
 147/CH3/EX3.5/Example3_5.sce   | 22 ++++++++++++++++++++++
 147/CH3/EX3.5/Result3_5.txt    |  7 +++++++
 147/CH3/EX3.6/Example3_6.sce   | 14 ++++++++++++++
 147/CH3/EX3.6/Result3_6.txt    |  3 +++
 147/CH3/EX3.7/Example3_7.sce   | 41 +++++++++++++++++++++++++++++++++++++++++
 147/CH3/EX3.7/Result3_7.txt    |  7 +++++++
 147/CH3/EX3.8/Example3_8.sce   | 25 +++++++++++++++++++++++++
 147/CH3/EX3.8/Result3_8.txt    |  5 +++++
 147/CH3/EX3.9/Example3_9.sce   | 18 ++++++++++++++++++
 147/CH3/EX3.9/Result3_9.txt    |  1 +
 24 files changed, 289 insertions(+)
 create mode 100644 147/CH3/EX3.13/Example3_13.sce
 create mode 100644 147/CH3/EX3.13/Result3_13.txt
 create mode 100644 147/CH3/EX3.14/Example3_14.sce
 create mode 100644 147/CH3/EX3.14/Result3_14.txt
 create mode 100644 147/CH3/EX3.18/Example3_18.sce
 create mode 100644 147/CH3/EX3.18/Result3_18.txt
 create mode 100644 147/CH3/EX3.19/Example3_19.sce
 create mode 100644 147/CH3/EX3.19/Result3_19.txt
 create mode 100644 147/CH3/EX3.2/Example3_2.sce
 create mode 100644 147/CH3/EX3.2/Result3_2.txt
 create mode 100644 147/CH3/EX3.3/Example3_3.sce
 create mode 100644 147/CH3/EX3.3/Result3_3.txt
 create mode 100644 147/CH3/EX3.4/Example3_4.sce
 create mode 100644 147/CH3/EX3.4/Result3_4.txt
 create mode 100644 147/CH3/EX3.5/Example3_5.sce
 create mode 100644 147/CH3/EX3.5/Result3_5.txt
 create mode 100644 147/CH3/EX3.6/Example3_6.sce
 create mode 100644 147/CH3/EX3.6/Result3_6.txt
 create mode 100644 147/CH3/EX3.7/Example3_7.sce
 create mode 100644 147/CH3/EX3.7/Result3_7.txt
 create mode 100644 147/CH3/EX3.8/Example3_8.sce
 create mode 100644 147/CH3/EX3.8/Result3_8.txt
 create mode 100644 147/CH3/EX3.9/Example3_9.sce
 create mode 100644 147/CH3/EX3.9/Result3_9.txt

(limited to '147/CH3')

diff --git a/147/CH3/EX3.13/Example3_13.sce b/147/CH3/EX3.13/Example3_13.sce
new file mode 100644
index 000000000..a43ff963a
--- /dev/null
+++ b/147/CH3/EX3.13/Example3_13.sce
@@ -0,0 +1,17 @@
+close();
+clear;
+clc;
+//resistance 'R', inductance 'L', angular frequency 'w', phase difference 'theta'
+R = 20; //ohm
+L = 15*10^(-3); //mH
+w = 1000; //rad/s
+theta = 45; //degree
+//net reactance of parallel LC combination = 20 ohm or -20 ohm
+C1 = poly(0,'C1');
+pc1 = 1/(%i*w*L) - 1/(R*%i) + %i*w*C1;
+C1 = roots(pc1);
+C2 = poly(0,'C2');
+pc2 = 1/(%i*w*L) + 1/(R*%i) + %i*w*C2;
+C2 = roots(pc2);
+mprintf("C1 = %0.2f uF\nC2 = %0.1f uF\n\n",C1*10^6, C2*10^6);
+mprintf("For smaller (larger) capacitance, I lags (leads) V by 45 degree"); 
\ No newline at end of file
diff --git a/147/CH3/EX3.13/Result3_13.txt b/147/CH3/EX3.13/Result3_13.txt
new file mode 100644
index 000000000..201c6fd76
--- /dev/null
+++ b/147/CH3/EX3.13/Result3_13.txt
@@ -0,0 +1,4 @@
+C1 = 16.67 uF
+C2 = 116.7 uF
+
+For smaller (larger) capacitance, I lags (leads) V by 45 degree
\ No newline at end of file
diff --git a/147/CH3/EX3.14/Example3_14.sce b/147/CH3/EX3.14/Example3_14.sce
new file mode 100644
index 000000000..58268bee5
--- /dev/null
+++ b/147/CH3/EX3.14/Example3_14.sce
@@ -0,0 +1,17 @@
+//Inductance L, Resistance of coil Rl, Voltage V, Frequency f
+close();
+clear;
+clc;
+L = 0.046;//H
+Rl = 10;//ohm
+V = 100;//V
+f = 60;//Hz
+omega = 2*%pi*f;
+Zl = Rl + %i*omega*L;
+Il = V/Zl;
+Ilr = polar(Il);
+Ilarg = atan(imag(Il),real(Il))*180/%pi;
+Power_factor = cos(Ilarg*%pi/180);
+//For unity power factor imaginary part of admittance must be zero hence Capacitance across the coil
+C = -imag(1/Zl)/omega;
+mprintf('Current drawn I = %0.1f arg(%0.0f degree)\nPower factor = %0.1f lagging \nCapacitance that must be connected across coil in order to make the power factor unity C = %0.0f micro F',Ilr,Ilarg,Power_factor,C*10^6);
\ No newline at end of file
diff --git a/147/CH3/EX3.14/Result3_14.txt b/147/CH3/EX3.14/Result3_14.txt
new file mode 100644
index 000000000..141dd6483
--- /dev/null
+++ b/147/CH3/EX3.14/Result3_14.txt
@@ -0,0 +1,3 @@
+Current drawn I = 5.0 arg(-60 degree)
+Power factor = 0.5 lagging 
+Capacitance that must be connected across coil in order to make the power factor unity C = 115 micro F
\ No newline at end of file
diff --git a/147/CH3/EX3.18/Example3_18.sce b/147/CH3/EX3.18/Example3_18.sce
new file mode 100644
index 000000000..ae7d91ad8
--- /dev/null
+++ b/147/CH3/EX3.18/Example3_18.sce
@@ -0,0 +1,15 @@
+//Phase voltage Vp, Impedance Z
+close();
+clear;
+clc;
+Vp = 120;//V
+Z = 36 + %i*48;//ohm
+//Line Voltage 'Vl'
+Vl = 3^(1/2)*Vp;
+Ip = Vp/polar(Z);
+Il = Ip;
+R = real(Z);
+//Power factor 'Pf'
+Pf = R/polar(Z);
+P = 3^(1/2)*Vl*Il*Pf;
+mprintf('Line voltage Vl = %0.1f V\nLine current Il = %0.0f A\nPower factor = %0.1f lagging \nPower P = %0.0f W',Vl,Il,Pf,P);
\ No newline at end of file
diff --git a/147/CH3/EX3.18/Result3_18.txt b/147/CH3/EX3.18/Result3_18.txt
new file mode 100644
index 000000000..125d432ba
--- /dev/null
+++ b/147/CH3/EX3.18/Result3_18.txt
@@ -0,0 +1,4 @@
+Line voltage Vl = 207.8 V
+Line current Il = 2 A
+Power factor = 0.6 lagging 
+Power P = 432 W
\ No newline at end of file
diff --git a/147/CH3/EX3.19/Example3_19.sce b/147/CH3/EX3.19/Example3_19.sce
new file mode 100644
index 000000000..e3c756f19
--- /dev/null
+++ b/147/CH3/EX3.19/Example3_19.sce
@@ -0,0 +1,24 @@
+close();
+clear;
+clc;
+//line voltage 'Vl', resistance 'R', reactance 'X'
+Vl = 207.8; 
+R = 36;
+X = 48;
+//inpedance
+Z = sqrt(R^2 + X^2);
+//(a)phase current 'Ip'
+Ip = Vl/Z; //A
+mprintf("Phase current, Ip = %0.2f A\n\n",Ip);
+
+//(b)line current 'Il'
+Il = sqrt(3)*Ip; //A
+mprintf("Line current, Il = %d A\n\n",round(Il));
+
+//(c)power factor 'pf'
+pf = R/Z;
+mprintf("Power factor, pf = %0.1f lagging\n\n",pf);
+
+//(d)total power 'P'
+P = sqrt(3)*Vl*Il*pf; //W
+mprintf("Total power, P = %d W",round(P));
diff --git a/147/CH3/EX3.19/Result3_19.txt b/147/CH3/EX3.19/Result3_19.txt
new file mode 100644
index 000000000..4980692d7
--- /dev/null
+++ b/147/CH3/EX3.19/Result3_19.txt
@@ -0,0 +1,7 @@
+Phase current, Ip = 3.46 A
+
+Line current, Il = 6 A
+
+Power factor, pf = 0.6 lagging
+
+Total power, P = 1295 W
\ No newline at end of file
diff --git a/147/CH3/EX3.2/Example3_2.sce b/147/CH3/EX3.2/Example3_2.sce
new file mode 100644
index 000000000..6bc54c7d6
--- /dev/null
+++ b/147/CH3/EX3.2/Example3_2.sce
@@ -0,0 +1,12 @@
+//V = Vmax*sin(omega*t) i = imax*sin(omega*t-phi)
+close();
+clear;
+clc;
+Vmax = 155.6;
+omega = 377;//Deg
+imax = 7.07;
+phi = 36.87;
+f = omega/(2*%pi);
+T = 1/f;
+phase_angle = phi*%pi/180;
+mprintf('frequency f = %0.0f Hz\nT = %0.4f s \nphase angle betwwen v and i = %0.2f rad',f,T,phase_angle);
\ No newline at end of file
diff --git a/147/CH3/EX3.2/Result3_2.txt b/147/CH3/EX3.2/Result3_2.txt
new file mode 100644
index 000000000..418829921
--- /dev/null
+++ b/147/CH3/EX3.2/Result3_2.txt
@@ -0,0 +1,3 @@
+frequency f = 60 Hz
+T = 0.0167 s 
+phase angle betwwen v and i = 0.64 rad
\ No newline at end of file
diff --git a/147/CH3/EX3.3/Example3_3.sce b/147/CH3/EX3.3/Example3_3.sce
new file mode 100644
index 000000000..c81b4d6c2
--- /dev/null
+++ b/147/CH3/EX3.3/Example3_3.sce
@@ -0,0 +1,26 @@
+close();
+clear;
+clc;
+//resistance 'R', rate of electrical energy 'rate', time of operation 't', time period 'T'
+T = 0.01; //s
+R = 20; //ohm
+rate = 0.06; //$ per kWh
+time = 24; //hours
+//from graph, voltage as function of time 'v'
+function V = v(t)
+  V = 10^4 * t; //V
+  endfunction
+//current as a function of time 'i'
+function I = i(t)
+  I = v(t)/R; //A
+endfunction
+//power as afunction of time 'p'
+function power = p(t)
+  power = v(t)*i(t); //W
+  endfunction
+//power 'P'
+P = 1/T * intg(0,T,p)/1000; //kW
+energy = P*time; //kWh
+//cost of electrical energy 'cost'
+cost = rate*energy; //$
+mprintf("Electrical energy costs $ %0.2f for %d hours",cost,time);
diff --git a/147/CH3/EX3.3/Result3_3.txt b/147/CH3/EX3.3/Result3_3.txt
new file mode 100644
index 000000000..f702a98f5
--- /dev/null
+++ b/147/CH3/EX3.3/Result3_3.txt
@@ -0,0 +1 @@
+Electrical energy costs $ 0.24 for 24 hours
\ No newline at end of file
diff --git a/147/CH3/EX3.4/Example3_4.sce b/147/CH3/EX3.4/Example3_4.sce
new file mode 100644
index 000000000..6acae06aa
--- /dev/null
+++ b/147/CH3/EX3.4/Example3_4.sce
@@ -0,0 +1,12 @@
+close();
+clear;
+clc;
+//time period 'T'
+T = 0.01; //s
+//voltage as a function of time 'v' 
+function volt = v(t)
+  volt = (10^4*t)^2;
+endfunction
+//RMS value of voltage 'V'
+V = sqrt(1/T * intg(0,T,v)); //V
+mprintf("RMS value of voltage, V = %f V",V);
\ No newline at end of file
diff --git a/147/CH3/EX3.4/Result3_4.txt b/147/CH3/EX3.4/Result3_4.txt
new file mode 100644
index 000000000..c018c147d
--- /dev/null
+++ b/147/CH3/EX3.4/Result3_4.txt
@@ -0,0 +1 @@
+RMS value of voltage, V = 57.735027 V
\ No newline at end of file
diff --git a/147/CH3/EX3.5/Example3_5.sce b/147/CH3/EX3.5/Example3_5.sce
new file mode 100644
index 000000000..021d0f13e
--- /dev/null
+++ b/147/CH3/EX3.5/Example3_5.sce
@@ -0,0 +1,22 @@
+close();
+clear;
+clc;
+Vm = 200;
+w = 377;
+Im = 8; //A
+theta = 30; //degree
+//(a) power factor 'pf'
+pf = cos(theta*(%pi)/180);
+mprintf("Power factor, pf = %0.3f\n\n",pf);
+
+//(b) true power 'S'
+S = Vm*Im/2 * pf; //W
+mprintf("True Power, S = %0.1f W\n\n",S);
+
+//(c)Apparant power 'Q'
+Q = Vm*Im/2; //VA
+mprintf("Apparant power, Q = %d VA\n\n",round(Q));
+
+//(d)Reactive Power 'P'
+P = Vm*Im*sqrt(1-pf^2)/2; //VAr
+mprintf("Reactive power, P = %d VAr",round(P));
\ No newline at end of file
diff --git a/147/CH3/EX3.5/Result3_5.txt b/147/CH3/EX3.5/Result3_5.txt
new file mode 100644
index 000000000..32d0f1fff
--- /dev/null
+++ b/147/CH3/EX3.5/Result3_5.txt
@@ -0,0 +1,7 @@
+Power factor, pf = 0.866
+
+True Power, S = 692.8 W
+
+Apparant power, Q = 800 VA
+
+Reactive power, P = 400 VAr
\ No newline at end of file
diff --git a/147/CH3/EX3.6/Example3_6.sce b/147/CH3/EX3.6/Example3_6.sce
new file mode 100644
index 000000000..686ecae74
--- /dev/null
+++ b/147/CH3/EX3.6/Example3_6.sce
@@ -0,0 +1,14 @@
+//Resistance of coil R, Current I, Voltage V, Frequency f
+close();
+clear;
+clc;
+R = 10;//ohm
+I = 5;//A
+V = 100;//V
+f = 60;//Hz
+Z = V/I;
+omega = 2*%pi*f;
+L = (Z^2-R^2)^(1/2)/omega*1000;
+Power_factor = R/Z;
+Q = V*I*sin(acos(R/Z));
+mprintf('Inductance of coil L = %0.2f mH\nPower Factor = %0.1f \nReactive Power = %0.0f var',L,Power_factor,Q);
\ No newline at end of file
diff --git a/147/CH3/EX3.6/Result3_6.txt b/147/CH3/EX3.6/Result3_6.txt
new file mode 100644
index 000000000..592db7a6c
--- /dev/null
+++ b/147/CH3/EX3.6/Result3_6.txt
@@ -0,0 +1,3 @@
+Inductance of coil L = 45.94 mH
+Power Factor = 0.5 
+Reactive Power = 433 var
\ No newline at end of file
diff --git a/147/CH3/EX3.7/Example3_7.sce b/147/CH3/EX3.7/Example3_7.sce
new file mode 100644
index 000000000..55860a4ec
--- /dev/null
+++ b/147/CH3/EX3.7/Example3_7.sce
@@ -0,0 +1,41 @@
+close();
+clear;
+clc;
+//voltage source 'V', frequency of source 'f', resistance 'R', inductance 'L', capacitance 'C'
+V = 100;
+V_arg = 0;
+f = 79.6; //Hz
+R = 100; //ohm
+L = 1; //H
+C = 5*10^(-6); //F
+
+//(a) 
+//angular frequency 'w'
+w = 2*%pi*f; //rad/s
+//inductive reactance 'Xl'
+Xl = w*L; //ohm
+Xl_arg = %pi/2; //rad
+//capacitive reactance 'Xc'
+Xc = 1/(w*C);
+Xc_arg = -%pi/2; //rad
+//impedance 'Z'
+Z = sqrt(R^2 + (Xl-Xc)^2);
+Z_arg = atan((Xl-Xc)/R); //rad
+//input current magnitude 'I' and argument 'I_arg'
+I = V/Z;
+I_arg = V_arg-Z_arg; //rad
+mprintf("input current, I = %0.3f arg(%d degree )\n\n",I,I_arg*180/%pi);
+
+//(b)
+//voltage across R 'Vr'
+Vr = R*I;
+Vr_arg = I_arg;
+//volatge across L 'Vl'
+Vl = Xl*I;
+Vl_arg = Xl_arg+I_arg;
+//voltage across C 'Vc' 
+Vc = Xc*I;
+Vc_arg = Xc_arg+I_arg;
+mprintf("Voltage across resistance, Vr = %0.1f arg(%d degree )\n\n",Vr, round(Vr_arg*180/%pi));
+mprintf("Voltage across inductor, Vl = %0.1f arg(%d degree )\n\n",Vl, round(Vl_arg*180/%pi));
+mprintf("Voltage across capacitor, Vc = %0.1f arg(%d degree )",Vc, round(Vc_arg*180/%pi));
\ No newline at end of file
diff --git a/147/CH3/EX3.7/Result3_7.txt b/147/CH3/EX3.7/Result3_7.txt
new file mode 100644
index 000000000..47106bdae
--- /dev/null
+++ b/147/CH3/EX3.7/Result3_7.txt
@@ -0,0 +1,7 @@
+input current, I = 0.706 arg(-45 degree )
+
+Voltage across resistance, Vr = 70.6 arg(-45 degree )
+
+Voltage across inductor, Vl = 353.2 arg(45 degree )
+
+Voltage across capacitor, Vc = 282.4 arg(-135 degree )
\ No newline at end of file
diff --git a/147/CH3/EX3.8/Example3_8.sce b/147/CH3/EX3.8/Example3_8.sce
new file mode 100644
index 000000000..a9cdcab70
--- /dev/null
+++ b/147/CH3/EX3.8/Example3_8.sce
@@ -0,0 +1,25 @@
+//Voltage V, Resistances R1 and R2, Capacitive reactance Xc
+close();
+clear;
+clc;
+R1 = 10;//ohm
+R2 = 20;
+Xc = -11.55*%i;
+V = 173.2;//V
+//By Nodal analysis
+V1 = V/(R1*((1/R1)+(1/R2)+(1/Xc)));
+V1r = polar(V1);
+V1arg = atan(imag(V1),real(V1))*180/%pi;
+V2 = V-V1;
+V2r = polar(V2);
+V2arg = atan(imag(V2),real(V2))*180/%pi;
+I = V2/R1;
+Ir = polar(I);
+Iarg =atan(imag(I),real(I))*180/%pi;
+I1 = V1/R2;
+I1r = polar(I1);
+I1arg =atan(imag(I1),real(I1))*180/%pi;
+I2 = V1/Xc;
+I2r = polar(I2);
+I2arg =atan(imag(I2),real(I2))*180/%pi;
+mprintf('Voltage across R2 and Xc = %0.0f arg(%0.0f degree) V \nVoltage across R1 = %0.0f arg(%0.0f degree) V \nCurrent through R1, I = %0.0f arg(%0.0f degree) A \nCurrent through R2, I1 = %0.0f arg(%0.0f degree) A \nCurrent through capacitor, I2 = %0.2f arg(%0.0f degree) A',V1r,V1arg,V2r,V2arg,Ir,Iarg,I1r,I1arg,I2r,I2arg);
\ No newline at end of file
diff --git a/147/CH3/EX3.8/Result3_8.txt b/147/CH3/EX3.8/Result3_8.txt
new file mode 100644
index 000000000..706e90368
--- /dev/null
+++ b/147/CH3/EX3.8/Result3_8.txt
@@ -0,0 +1,5 @@
+Voltage across R2 and Xc = 100 arg(-30 degree) V 
+Voltage across R1 = 100 arg(30 degree) V 
+Current through R1, I = 10 arg(30 degree) A 
+Current through R2, I1 = 5 arg(-30 degree) A 
+Current through capacitor, I2 = 8.66 arg(60 degree) A
\ No newline at end of file
diff --git a/147/CH3/EX3.9/Example3_9.sce b/147/CH3/EX3.9/Example3_9.sce
new file mode 100644
index 000000000..9851e8674
--- /dev/null
+++ b/147/CH3/EX3.9/Example3_9.sce
@@ -0,0 +1,18 @@
+close();
+clear;
+clc;
+//at node I, 
+R1 = 50; //ohm
+R2 = 50; //ohm
+X = 25*(%i); //ohm
+V2 = 100; //V
+V2_arg = 0;
+V2 = V2*cos(V2_arg) + (%i)*sin(V2_arg);
+I1 = 2.83;
+I1_arg = 45*%pi/180;
+I1 = I1*cos(I1_arg) + (%i)*I1*sin(I1_arg);
+V1 = poly(0,'V1');
+pv = I1 - V1/R1 - V1/X - (V1-V2)/R2;
+V1 = roots(pv);
+I = (V2-V1)/R2;
+disp(I,"Current from the voltage source, I = ");
\ No newline at end of file
diff --git a/147/CH3/EX3.9/Result3_9.txt b/147/CH3/EX3.9/Result3_9.txt
new file mode 100644
index 000000000..8541f18e0
--- /dev/null
+++ b/147/CH3/EX3.9/Result3_9.txt
@@ -0,0 +1 @@
+Current from the voltage source, I = 1.5 - 1.5005561i  
-- 
cgit