From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 2441/CH5/EX5.10/Ex5_10.sce | 30 ++++++++++++++++++++++++++++++
 1 file changed, 30 insertions(+)
 create mode 100755 2441/CH5/EX5.10/Ex5_10.sce

(limited to '2441/CH5/EX5.10')

diff --git a/2441/CH5/EX5.10/Ex5_10.sce b/2441/CH5/EX5.10/Ex5_10.sce
new file mode 100755
index 000000000..0248833b8
--- /dev/null
+++ b/2441/CH5/EX5.10/Ex5_10.sce
@@ -0,0 +1,30 @@
+//Example 5.10
+clc;clear;close;
+format('v',7);
+s=poly(0,'s');//for transfer function
+Tg=0.2;//sec///time constant of governing system
+Tt=2;//sec///time constant of turbine
+Gr=1/(1+Tg*s);//Transfer function of governer
+Gt=1/(1+Tt*s);//Transfer function of turbine
+C=1500;//MW
+f=50;//Hz
+R=4;//%////Speed regulation constant
+H=5;//sec////Inertia constant
+delPL=1;//%////change in load
+delf=1;//%////change in frequency
+disp("Part(a)");
+R=R/100*f;//z/p.u.MW
+D=delPL/delf*C/f;//MW/Hz
+D=D/C;//p.u.MW/Hz
+Kp=1/D;//Hz/p.u.MW
+Tp=2*H/f/D;//sec
+Gp=Kp/(1+Tp*s);//Transfer function of power system
+delFs=-Gp/(1+Gr*Gt*Gp/R);
+disp(delFs,"delFs = M/s*");
+disp("Part(b)");
+delf0_by_M=-Kp/(1+Kp/R);//Hz
+delf0=delf/100*f;//Hz
+M=delf0/delf0_by_M;//p.u.MW
+M=M*C;//MW
+disp(M,"Largest step change(MW)");
+//Transfer functions multiplication Gr*Gt*Gp is calculated & it is not possible to show together without calculated as in the book.
-- 
cgit