diff options
Diffstat (limited to '3802/CH12/EX12.9/Ex12_9.sce')
-rw-r--r-- | 3802/CH12/EX12.9/Ex12_9.sce | 72 |
1 files changed, 72 insertions, 0 deletions
diff --git a/3802/CH12/EX12.9/Ex12_9.sce b/3802/CH12/EX12.9/Ex12_9.sce new file mode 100644 index 000000000..e8f84ce8a --- /dev/null +++ b/3802/CH12/EX12.9/Ex12_9.sce @@ -0,0 +1,72 @@ +//Book Name:Fundamentals of Electrical Engineering
+//Author:Rajendra Prasad
+//Publisher: PHI Learning Private Limited
+//Edition:Third ,2014
+
+//Ex12_9.sce
+
+clc;
+clear;
+z=10;
+ang1=0;
+ang2=37;
+ang3=-53;
+Zr=complex(z*cosd(ang1),z*sind(ang1));
+Zy=complex(z*cosd(ang2),z*sind(ang2));
+Zb=complex(z*cosd(ang3),z*sind(ang3));
+
+printf("\n (a)For phase sequence RYB:\n")
+V=220;
+theta1=0;
+theta2=-120;
+theta3=120;
+Vr=complex(V*cosd(theta1),V*sind(theta1));
+Vy=complex(V*cosd(theta2),V*sind(theta2));
+Vb=complex(V*cosd(theta3),V*sind(theta3));
+
+Ir=Vr/Zr;
+Iy=Vy/Zy;
+Ib=Vb/Zb;
+In=Ir+Iy+Ib;
+In_mag=sqrt(real(In)^2+imag(In)^2);
+In_angle=atand(imag(In)/real(In));
+printf("\n The current through the neutral wire,\n -In=\tMagnitude\tAngle(deg) \n\t %2.2f \t %2.2f \n",In_mag,In_angle)
+
+Ir_mag=sqrt(real(Ir)^2+imag(Ir)^2);
+Iy_mag=sqrt(real(Iy)^2+imag(Iy)^2);
+Ib_mag=sqrt(real(Ib)^2+imag(Ib)^2);
+Pr=(Ir_mag)^2*real(Zr);
+Py=(Iy_mag)^2*real(Zy);
+Pb=(Ib_mag)^2*real(Zb);
+printf("\n Power taken by each load:\n\t Pr=%d W \n\t Py=%4.1f W \n\t Pb=%4.1f W \n", Pr, Py, Pb)
+
+
+
+printf("\n\n (b)For phase sequence RBY:\n")
+V=220;
+theta1=0;
+theta2=120;
+theta3=-120;
+Vr=complex(V*cosd(theta1),V*sind(theta1));
+Vy=complex(V*cosd(theta2),V*sind(theta2));
+Vb=complex(V*cosd(theta3),V*sind(theta3));
+
+Ir=Vr/Zr;
+Iy=Vy/Zy;
+Ib=Vb/Zb;
+In=Ir+Iy+Ib;
+In_mag=sqrt(real(In)^2+imag(In)^2);
+In_angle=atand(imag(In)/real(In));
+printf("\n The current through the neutral wire,\n In=\tMagnitude\tAngle(deg) \n\t %2.2f \t %2.2f \n",In_mag,In_angle)
+
+Ir_mag=sqrt(real(Ir)^2+imag(Ir)^2);
+Iy_mag=sqrt(real(Iy)^2+imag(Iy)^2);
+Ib_mag=sqrt(real(Ib)^2+imag(Ib)^2);
+Pr=(Ir_mag)^2*real(Zr);
+Py=(Iy_mag)^2*real(Zy);
+Pb=(Ib_mag)^2*real(Zb);
+printf("\n Power taken by each load:\n\t Pr=%d W \n\t Py=%4.1f W \n\t Pb=%4.1f W \n", Pr, Py, Pb)
+
+
+
+
|