initial commit / add all books

24 files changed, 289 insertions, 0 deletions

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