diff options
Diffstat (limited to '1445/CH2/EX2.15/ch2_ex_15.sce')
-rw-r--r-- | 1445/CH2/EX2.15/ch2_ex_15.sce | 39 |
1 files changed, 39 insertions, 0 deletions
diff --git a/1445/CH2/EX2.15/ch2_ex_15.sce b/1445/CH2/EX2.15/ch2_ex_15.sce new file mode 100644 index 000000000..09182c5ed --- /dev/null +++ b/1445/CH2/EX2.15/ch2_ex_15.sce @@ -0,0 +1,39 @@ +//CHAPTER 2- STEADY-STATE ANALYSIS OF SINGLE-PHASE A.C. CIRCUIT +//Example 15 + +disp("CHAPTER 2"); +disp("EXAMPLE 15"); + +//VARIABLE INITIALIZATION +I=2; //in Amperes +angle_I=60; //in degrees +v1=200; //in Volts +f=50; //in Hertz + +//SOLUTION +z1=v1/I; +disp(sprintf("The impedance is %d Ω, %d degrees",z1,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; +[r,x1]=pol2rect(z1,angle_I); +disp(sprintf("The resistance is %d Ω",r)); +L=x1/(2*%pi*f); +disp(sprintf("The inductance is %f H",L)); + +v2=100; +f2=25; +x2=2*%pi*f2*L; +z2=sqrt((r^2)+(x2^2)); +angle=atan(x2/r); +I1=v2/z2; +p=v2*I1*cos(-angle); +disp(sprintf("The power consumed is %f W",p)); + +//Answer may be slightly different due to precision of floating point numbers + +//END + + |