1 files changed, 45 insertions, 0 deletions
diff --git a/2273/CH5/EX5.10/ex5_10.sce b/2273/CH5/EX5.10/ex5_10.sce
new file mode 100755
index 000000000..77f788bdc
--- /dev/null
+++ b/2273/CH5/EX5.10/ex5_10.sce
@@ -0,0 +1,45 @@
+//Find regulation and charging current using nominal T method
+clear;
+clc;
+//soltion
+//FUNCTIONS
+function [z]=rxr(A,B)//Function for the multiplication of rectangular
+    z(1)=A(1)*B(1)
+    z(2)=A(2)+B(2)
+    endfunction
+
+function [a]=r2p(z)//Function for rectangular to polar
+    a=z(1)*complex(cosd(z(2)),sind(z(2)))
+endfunction
+
+//given
+P=50*10^6;//M watts
+Vrl=132*1000;//V
+pf=0.8//power factir
+Vr=[floor(Vrl/sqrt(3)) 0];
+Ir=[floor(P/(sqrt(3)*Vrl*pf)) -acosd(pf)];
+A=[0.95 1.4];
+B=[96 78];
+C=[0.0015 90];
+D=A;
+Z1=rxr(A,Vr);
+Z2=rxr(B,Ir);
+AV=r2p(Z1);
+BI=r2p(Z2);
+Vs=AV+BI;
+theta1=atand((imag(Vs)/real(Vs)));
+printf("Sending end voltage= %.0f∠%.2f° Volts\n",abs(Vs),theta1);
+Y1=rxr(C,Vr);
+Y2=rxr(D,Ir);
+CV=r2p(Y1);
+DI=r2p(Y2);
+Is=CV+DI;
+Ira=r2p(Ir);
+Ic=Is-Ira;
+theta2=atand(imag(Ic)/real(Ic));
+Ic_=sqrt(round(imag(Ic))^2+round(real(Ic))^2);
+printf("Charging current= %.1f∠%f° A\n",Ic_,theta2);
+Pr=((abs(Vs)/A(1)-Vr)*100)/Vr;
+printf("Percentage regulation= %.0f percent\n",Pr);
+//1. The Magnitude of Sending end voltage is 94066, it is due to rounding some of the values
+//2. The angle in the book is 93.1° in charging current