summaryrefslogtreecommitdiff
path: root/1445/CH2/EX2.17
diff options
context:
space:
mode:
Diffstat (limited to '1445/CH2/EX2.17')
-rw-r--r--1445/CH2/EX2.17/ch2_ex_17.sce51
1 files changed, 51 insertions, 0 deletions
diff --git a/1445/CH2/EX2.17/ch2_ex_17.sce b/1445/CH2/EX2.17/ch2_ex_17.sce
new file mode 100644
index 000000000..ba82080a2
--- /dev/null
+++ b/1445/CH2/EX2.17/ch2_ex_17.sce
@@ -0,0 +1,51 @@
+//CHAPTER 2- STEADY-STATE ANALYSIS OF SINGLE-PHASE A.C. CIRCUIT
+//Example 17
+
+disp("CHAPTER 2");
+disp("EXAMPLE 17");
+
+//VARIABLE INITIALIZATION
+e=141.4; //in Volts
+E=141.4/sqrt(2); //in Volts
+angle_E=0; //in degrees
+//i(t)=(14.14<0)+(7.07<120)
+i1=14.14; //in Amperes
+angle_i1=0; //in degrees
+i2=7.07; //in Amperes
+angle_i2=120; //in degrees
+
+//SOLUTION
+//function to convert from polar form to rectangular form
+function [x,y]=pol2rect(mag,angle);
+x=mag*cos(angle*(%pi/180)); //to convert the angle from degrees to radians
+y=mag*sin(angle*(%pi/180));
+endfunction;
+[i1_x,i1_y]=pol2rect(i1,angle_i1);
+[i2_x,i2_y]=pol2rect(i2,angle_i2);
+i=(i1_x+i2_x)+(%i*(i1_y+i2_y));
+//function to convert from rectangular form to polar form
+function [mag,angle]=rect2pol(x,y);
+mag=sqrt((x^2)+(y^2));
+angle=atan(y/x)*(180/%pi); //to convert the angle from radians to degrees
+endfunction;
+[I,angle_I]=rect2pol((i1_x+i2_x),(i1_y+i2_y));
+I=I/sqrt(2);
+
+//solution (i)
+z=E/I;
+angle_z=angle_E-angle_I;
+[r,xc]=pol2rect(z,angle_z);
+f=50;
+c=1/(2*%pi*f*(-xc));
+disp(sprintf("(i) The value of resistance is %f Ω",r));
+disp(sprintf(" The value of capacitance is %f μF",c*10^6));
+
+//solution (ii)
+pf=cos(angle_z*(%pi/180));
+disp(sprintf("(ii) The power factor is %f ",pf));
+
+//solution (iii)
+p=E*I*pf;
+disp(sprintf("(iii) The power absorbed by the source is %f W",p));
+
+//END