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.21/Ex6_21.sce | 28 ++++++++++++++++++++++++++++
 1 file changed, 28 insertions(+)
 create mode 100644 3871/CH6/EX6.21/Ex6_21.sce

(limited to '3871/CH6/EX6.21/Ex6_21.sce')

diff --git a/3871/CH6/EX6.21/Ex6_21.sce b/3871/CH6/EX6.21/Ex6_21.sce
new file mode 100644
index 000000000..a2ab474e8
--- /dev/null
+++ b/3871/CH6/EX6.21/Ex6_21.sce
@@ -0,0 +1,28 @@
+//============================================================================
+//Chapter 6 Example 21
+
+
+clc;
+clear all;
+
+//variable declaration
+R       = 25;           //rate burden in VA
+Is      = 5;            //current in A
+r       = 6;            //Rs/Es ratio of resistance to reactance 
+IL      = 0.2;          //iron loss in W
+Im         = 1.5;       //magnetising compnent of current in A
+
+
+//calculations
+KT      = 100/5;            //turn ratio
+Es      = R/Is;             // Secondary rated voltage in V
+Zs      = Es/Is;            //total secondary impedance in Ω
+theta       = (atan(1/r))*180/%pi;        //angle in °
+Zs1          = (Zs*cos(theta*%pi/180))+(Zs*sin(theta*%pi/180))*%i;
+Ie          = KT*0.04;          //energy component of exciting current on primary side 
+r           = (((Im*sin(theta*%pi/180))+(Ie*cos(theta*%pi/180)))/((KT*Is)+(Ie*cos(theta*%pi/180))+(Im*sin(theta*%pi/180))))*100;      //percentage ratio error in %
+beta        = (180/%pi)*(((Im*cos(theta*%pi/180))-(Ie*sin(theta*%pi/180)))/(KT*Is));        //phase angle erro in  °
+
+//result
+mprintf("percentage ratio error  = -%3.1f percentage",r);
+mprintf("\nphase angle error = %3.4f °",beta);
-- 
cgit