summaryrefslogtreecommitdiff
path: root/1445/CH2/EX2.47/Ex2_47.sce
diff options
context:
space:
mode:
Diffstat (limited to '1445/CH2/EX2.47/Ex2_47.sce')
-rw-r--r--1445/CH2/EX2.47/Ex2_47.sce48
1 files changed, 48 insertions, 0 deletions
diff --git a/1445/CH2/EX2.47/Ex2_47.sce b/1445/CH2/EX2.47/Ex2_47.sce
new file mode 100644
index 000000000..1e5f6185b
--- /dev/null
+++ b/1445/CH2/EX2.47/Ex2_47.sce
@@ -0,0 +1,48 @@
+//CHAPTER 2- STEADY-STATE ANALYSIS OF SINGLE-PHASE A.C. CIRCUIT
+//Example 47
+
+disp("CHAPTER 2");
+disp("EXAMPLE 47");
+
+//given
+//current in the circuit is 5-j.10 A
+//VARIABLE INITIALIZATION
+r=10; //in Ohms
+xl=8.66; //in Ohms
+I=5-(%i*10); //in Amperes
+
+//SOLUTION
+z=r+(%i*(xl));
+//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;
+[z,angle_z]=rect2pol(real(z),imag(z));
+[I,angle_I]=rect2pol(real(I),imag(I));
+
+//solution(i)
+v=I*z;
+angle_v=angle_I+angle_z;
+disp(sprintf("(i) The applied voltage is %.1f V, %.2f degrees",v,angle_v));
+
+//solution (ii)
+phi=angle_I-angle_v;
+pf=cos(phi*(%pi/180));
+disp(sprintf("(ii) The power factor is %.3f (lagging)",pf));
+
+//solution(iii)
+//S=phasor voltageXconjugate of phasor current
+//Converting v and I from polar to rectangular form
+s=v*I;
+angle_s=angle_v-angle_I;
+//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;
+[p,q]=pol2rect(s,angle_s);
+disp(sprintf("(iii) The active power is %.2f W",p));
+disp(sprintf(" The reactive power is %.2f VAR",q));
+
+//END