From 7bc77cb1ed33745c720952c92b3b2747c5cbf2df Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Sat, 3 Feb 2018 11:01:52 +0530
Subject: Added new code

---
 3871/CH6/EX6.26/Ex6_26.sce | 35 +++++++++++++++++++++++++++++++++++
 1 file changed, 35 insertions(+)
 create mode 100644 3871/CH6/EX6.26/Ex6_26.sce

(limited to '3871/CH6/EX6.26/Ex6_26.sce')

diff --git a/3871/CH6/EX6.26/Ex6_26.sce b/3871/CH6/EX6.26/Ex6_26.sce
new file mode 100644
index 000000000..0d15d7ff9
--- /dev/null
+++ b/3871/CH6/EX6.26/Ex6_26.sce
@@ -0,0 +1,35 @@
+//===========================================================================
+//chapter 6 example 26
+clc;
+clear all;
+
+//variable declaration
+KT	= 199;		//turn ration
+Is	= 5;		//secondary current in A
+Im	= 7;		//magnetising component of exciting current in A
+Ie	= 4;		//cross-loss component of exciting current in A
+delta	=0;
+
+//calculations
+KN	= 1000/5             //nominal ratio
+alpha	=atan(Ie/Im)	//angle in ° 
+alpha1  = (alpha*180)/%pi     
+theta	= delta+(((acos(0.8))*180)/%pi)     //from figure taken the value of gamma
+z   = cos((theta*%pi)/180)
+z1  = sin(((theta)*%pi)/180)
+Kc	= KT+(((Ie*z)+(Im*z1))/Ie)  //actual current in A
+e	= ((Kn-Kc)/Kc)*100	//ratio error
+b	=(180/%pi)*(((Im*z)-(Ie*z1))/(KT*Is))
+theta1	= delta-(((acos(0.8))*180)/%pi)    //from figure taken the value of gamma
+z11   = cos((theta1*%pi)/180)
+z12  = sin(((theta1*%pi)/180))
+Kc1	= KT+(((Ie*z11)+(Im*z12))/Is)	//actual current in A
+e1	= ((Kn-Kc1)/Kc1)*100	//ratio error
+b1	=(180/%pi)*(((Im*z11)-(Ie*z12))/(KT*Is))
+
+
+//result
+mprintf("ratio error = %3.2f  percentage",e);
+mprintf("\nphase angle error = %3.1f °",b);
+mprintf("\nratio error = %3.2f percentage",e1);
+mprintf("\nphase angle error = %3.2f percentage°",b1);
-- 
cgit