1 files changed, 38 insertions, 0 deletions
diff --git a/147/CH13/EX13.8/Example13_8.sce b/147/CH13/EX13.8/Example13_8.sce
new file mode 100644
index 000000000..571864d8c
--- /dev/null
+++ b/147/CH13/EX13.8/Example13_8.sce
@@ -0,0 +1,38 @@
+//Rated Power Pr, Resistance of primary winding R1,Resistance of secondary winding R2
+//Primary leakage reactance X1, Secondary leakage reactance X2
+//Magnetizing reactance Xm, Resistance accounting for core loss Rc
+//Secondary terminal voltage V2
+close();
+clear;
+clc;
+Po = 150;//kVA
+R1 = 0.2;//ohm
+R2 = 2/1000;
+X1 = 0.45;//ohm
+X2 = 4.5/1000;
+Rc = 10000;
+Xm = 1550;
+V2 = 240;
+a = 10;
+Pf =  0.8;
+theta2 = -acos(Pf);
+I2r = Po*1000/V2;
+I2 = I2r*cos(theta2)+I2r*sin(theta2)*%i;
+E1 = a*V2+ I2/a*(R1+%i*X1);
+E1r = polar(E1);
+E1_arg = atan(imag(E1),real(E1))*180/%pi;
+Im = E1/(%i*Xm);
+Ic = E1/Rc;
+Icr = polar(Ic);
+Io = Ic+Im;
+I1 = Io+I2/a;
+I1r = polar(I1);
+I1_arg = atan(imag(I1),real(I1))*180/%pi;
+V1 = E1 + I1*(R1+%i*X1);
+V1r = polar(V1);
+V1_arg = atan(imag(V1),real(V1))*180/%pi;
+per_regulation = (V1r - a*V2)/(a*V2)*100;
+output = Po*1000*Pf;
+losses = I1r^2*R1+Icr^2*Rc+I2r^2*R2;
+n = output/(output+losses)*100;
+mprintf('Voltage regulation = %0.1f %%\nEfficiency = %0.1f %%',per_regulation,n);