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.1/Ex5_1.sce   | 15 +++++++++++++++
 2441/CH5/EX5.10/Ex5_10.sce | 30 ++++++++++++++++++++++++++++++
 2441/CH5/EX5.11/Ex5_11.sce | 22 ++++++++++++++++++++++
 2441/CH5/EX5.12/Ex5_12.sce | 34 ++++++++++++++++++++++++++++++++++
 2441/CH5/EX5.13/Ex5_13.sce | 21 +++++++++++++++++++++
 2441/CH5/EX5.14/Ex5_14.sce | 18 ++++++++++++++++++
 2441/CH5/EX5.15/Ex5_15.sce | 14 ++++++++++++++
 2441/CH5/EX5.16/Ex5_16.sce | 15 +++++++++++++++
 2441/CH5/EX5.17/Ex5_17.sce | 11 +++++++++++
 2441/CH5/EX5.18/Ex5_18.sce | 18 ++++++++++++++++++
 2441/CH5/EX5.19/Ex5_19.sce | 16 ++++++++++++++++
 2441/CH5/EX5.2/Ex5_2.sce   | 15 +++++++++++++++
 2441/CH5/EX5.20/Ex5_20.sce | 21 +++++++++++++++++++++
 2441/CH5/EX5.21/Ex5_21.sce | 15 +++++++++++++++
 2441/CH5/EX5.22/Ex5_22.sce | 20 ++++++++++++++++++++
 2441/CH5/EX5.23/Ex5_23.sce | 13 +++++++++++++
 2441/CH5/EX5.24/Ex5_24.sce | 28 ++++++++++++++++++++++++++++
 2441/CH5/EX5.3/Ex5_3.sce   | 15 +++++++++++++++
 2441/CH5/EX5.4/Ex5_4.sce   | 23 +++++++++++++++++++++++
 2441/CH5/EX5.5/Ex5_5.sce   | 18 ++++++++++++++++++
 2441/CH5/EX5.6/Ex5_6.sce   | 24 ++++++++++++++++++++++++
 2441/CH5/EX5.7/Ex5_7.sce   | 15 +++++++++++++++
 2441/CH5/EX5.8/Ex5_8.sce   | 24 ++++++++++++++++++++++++
 2441/CH5/EX5.9/Ex5_9.sce   | 27 +++++++++++++++++++++++++++
 24 files changed, 472 insertions(+)
 create mode 100755 2441/CH5/EX5.1/Ex5_1.sce
 create mode 100755 2441/CH5/EX5.10/Ex5_10.sce
 create mode 100755 2441/CH5/EX5.11/Ex5_11.sce
 create mode 100755 2441/CH5/EX5.12/Ex5_12.sce
 create mode 100755 2441/CH5/EX5.13/Ex5_13.sce
 create mode 100755 2441/CH5/EX5.14/Ex5_14.sce
 create mode 100755 2441/CH5/EX5.15/Ex5_15.sce
 create mode 100755 2441/CH5/EX5.16/Ex5_16.sce
 create mode 100755 2441/CH5/EX5.17/Ex5_17.sce
 create mode 100755 2441/CH5/EX5.18/Ex5_18.sce
 create mode 100755 2441/CH5/EX5.19/Ex5_19.sce
 create mode 100755 2441/CH5/EX5.2/Ex5_2.sce
 create mode 100755 2441/CH5/EX5.20/Ex5_20.sce
 create mode 100755 2441/CH5/EX5.21/Ex5_21.sce
 create mode 100755 2441/CH5/EX5.22/Ex5_22.sce
 create mode 100755 2441/CH5/EX5.23/Ex5_23.sce
 create mode 100755 2441/CH5/EX5.24/Ex5_24.sce
 create mode 100755 2441/CH5/EX5.3/Ex5_3.sce
 create mode 100755 2441/CH5/EX5.4/Ex5_4.sce
 create mode 100755 2441/CH5/EX5.5/Ex5_5.sce
 create mode 100755 2441/CH5/EX5.6/Ex5_6.sce
 create mode 100755 2441/CH5/EX5.7/Ex5_7.sce
 create mode 100755 2441/CH5/EX5.8/Ex5_8.sce
 create mode 100755 2441/CH5/EX5.9/Ex5_9.sce

(limited to '2441/CH5')

diff --git a/2441/CH5/EX5.1/Ex5_1.sce b/2441/CH5/EX5.1/Ex5_1.sce
new file mode 100755
index 000000000..d876f6ae0
--- /dev/null
+++ b/2441/CH5/EX5.1/Ex5_1.sce
@@ -0,0 +1,15 @@
+//Example 5.1
+clc;clear;close;
+P=100;//MW
+drop=4;//%(No load to full load drop)
+f=50;//Hz
+disp("Part(i)");
+p=1;//MW(For calculating per unit MW)
+R=(drop/100)*f/p;//Hz/p.u.MW
+disp(R,"Speed regulation in Hz/p.u.MW");
+R=(drop/100)*f/P;//Hz/MW
+disp(R,"Speed regulation in Hz/MW");
+disp("Part(ii)");
+del_f=-0.1;//Hz(Frequency drop)
+delP=-1/R*del_f;//MW(Change in power output)
+disp(delP,"Change in power output(MW)");
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.
diff --git a/2441/CH5/EX5.11/Ex5_11.sce b/2441/CH5/EX5.11/Ex5_11.sce
new file mode 100755
index 000000000..3cbf38e28
--- /dev/null
+++ b/2441/CH5/EX5.11/Ex5_11.sce
@@ -0,0 +1,22 @@
+//Example 5.11
+clc;clear;close;
+format('v',8);
+GA=5000;//MW
+GB=10000;//MW
+R=2;//Hz/p.u.MW////Speed regulation constant
+D=0.01;//p.u.MW/Hz
+Ls=100;//MW//Load increase
+RA=R*GB/GA;//Hz/p.u.MW
+DA=D*GA/GB;//p.u.MW/Hz
+RB=R;//Hz/p.u.MW
+DB=D;//p.u.MW/Hz
+Beta_A=DA+1/RA;//p.u.MW/Hz
+Beta_B=DB+1/RB;//p.u.MW/Hz
+MA=0;//Load increase
+MB=Ls/GB;//p.u.MW
+delf0=-MB/(Beta_A+Beta_B);//Hz
+disp(delf0,"Static frequency drop(Hz)");
+format('v',6);
+delPAB=Beta_A*MB/(Beta_A+Beta_B);//p.u.MW
+delPAB=delPAB*GB;//MW
+disp(delPAB,"Change in tie line power(MW)");
diff --git a/2441/CH5/EX5.12/Ex5_12.sce b/2441/CH5/EX5.12/Ex5_12.sce
new file mode 100755
index 000000000..61f7c848f
--- /dev/null
+++ b/2441/CH5/EX5.12/Ex5_12.sce
@@ -0,0 +1,34 @@
+//Example 5.12
+clc;clear;close;
+format('v',8);
+GA=500;//MW
+GB=2000;//MW
+RA=2.5;//Hz/p.u.MW////Speed regulation constant
+RB=2;//Hz/p.u.MW////Speed regulation constant
+Ls=20;//MW//Load increase
+f=50;//Hz
+delL=1;//%////change in load
+delf=1;//%////change in frequency
+DA=delL/delf*GA/f;//MW/Hz
+DA=DA/GB;//p.u.MW/Hz
+DB=delL/delf*GB/f;//MW/Hz
+DB=DB/GB;//p.u.MW/Hz
+RA=RA*GB/GA;//Hz/p.u.MW
+Beta_A=DA+1/RA;//p.u.MW/Hz
+Beta_B=DB+1/RB;//p.u.MW/Hz
+disp("Part(a)");
+MA=Ls/GB;//unitless
+MB=0;//unitless
+delf0=-MA/(Beta_A+Beta_B);//Hz
+disp(delf0,"Change in frequency(Hz)");
+delPAB=-Beta_B*MA/(Beta_B+Beta_A);//p.u.MW
+delPAB=delPAB*GB;//MW
+disp(delPAB,"Change in tie line power(MW)");
+disp("Part(b)");
+MB=Ls/GB;//unitless
+MA=0;//unitless
+delf0=-MB/(Beta_A+Beta_B);//Hz
+disp(delf0,"Change in frequency(Hz)");
+delPAB=Beta_A*MB/(Beta_B+Beta_A);//p.u.MW
+delPAB=delPAB*GB;//MW
+disp(delPAB,"Change in tie line power(MW)");
diff --git a/2441/CH5/EX5.13/Ex5_13.sce b/2441/CH5/EX5.13/Ex5_13.sce
new file mode 100755
index 000000000..e15752a36
--- /dev/null
+++ b/2441/CH5/EX5.13/Ex5_13.sce
@@ -0,0 +1,21 @@
+//Example 5.13
+clc;clear;close;
+format('v',5);
+G=4000;//MW
+R=2;//Hz/p.u.MW////Speed regulation constant
+H=5;//sec
+C=600;//MW//Capacity
+theta=40;//degree///Power angle
+f=50;//Hz
+disp("Part(a)");
+T=C/G*cosd(theta);//sec
+omega0=sqrt([2*%pi*f*T/H-(f/4/R/H)^2]);//radian/sec
+disp(omega0,"Frequency of oscillation(radian/sec)");
+disp("Part(b)");
+delLB=100;//MW//change in load in area B
+delPAB=delLB/2;//MW//because Beta_A=Beta_B
+disp(delPAB,"Change in tie line power(MW)");
+disp("Part(c)");
+format('v',6);
+omega0=sqrt([2*%pi*f*T/H]);//radian/sec
+disp(omega0,"Frequency of oscillation(radian/sec)");
diff --git a/2441/CH5/EX5.14/Ex5_14.sce b/2441/CH5/EX5.14/Ex5_14.sce
new file mode 100755
index 000000000..ea9c0f992
--- /dev/null
+++ b/2441/CH5/EX5.14/Ex5_14.sce
@@ -0,0 +1,18 @@
+//Example 5.14
+clc;clear;close;
+format('v',6);
+C1=300;//MW
+C2=400;//MW
+G1=4;//%//droop characteristics of governer
+G2=5;//%//droop characteristics of governer
+L=600;//MW
+f=50;//Hz
+//Load on first generator =L1
+//Load on second generator =L-L1
+//f-G1*f/100*(L1/C1)=f-G2*f/100*(L2/C2)
+L1=G2*L/C2/(G1/C1+G2/C2);//MW
+L2=L-L1;//MW
+disp(L1,"Load on first generator(MW)");
+disp(L2,"Load on second generator(MW)");
+fLoad=f*(1-L1/C1*G1/100);//Hz
+disp(fLoad,"Frequency at load(Hz)");
diff --git a/2441/CH5/EX5.15/Ex5_15.sce b/2441/CH5/EX5.15/Ex5_15.sce
new file mode 100755
index 000000000..70dcdecff
--- /dev/null
+++ b/2441/CH5/EX5.15/Ex5_15.sce
@@ -0,0 +1,14 @@
+//Example 5.15 
+clc;clear;close;
+format('v',6);
+G=100;//MVA
+f=50;//Hz
+delL=50;//MW
+Tc=0.4;//sec
+H=5;///kWs/kVA
+KE=G*1000*H;//kWs
+delKE=delL*1000*Tc;////kWs///due to decrease in load
+fnew=sqrt((KE+delKE)/KE) *f;//Hz
+fdev=(fnew-f)/f*100;//%
+disp(fnew,"New frequency(Hz)");
+disp(fdev,"Frequency deviation(%)");
diff --git a/2441/CH5/EX5.16/Ex5_16.sce b/2441/CH5/EX5.16/Ex5_16.sce
new file mode 100755
index 000000000..d43b1991d
--- /dev/null
+++ b/2441/CH5/EX5.16/Ex5_16.sce
@@ -0,0 +1,15 @@
+//Example 5.16
+clc;clear;close;
+format('v',7);
+G=100;//MVA
+f=50;//Hz
+delL=60;//MW
+Tc=0.35;//sec
+H=5;///kWs/kVA
+KE=G*1000*H;//kWs
+delKE=(G-delL)*1000*Tc;////kWs///due to decrease in load
+fnew=sqrt((KE+delKE)/KE) *f;//Hz
+fdev=(fnew-f)/f*100;//%
+disp(fnew,"New frequency(Hz)");
+format('v',6);
+disp(fdev,"Frequency deviation(%)");
diff --git a/2441/CH5/EX5.17/Ex5_17.sce b/2441/CH5/EX5.17/Ex5_17.sce
new file mode 100755
index 000000000..32be3be11
--- /dev/null
+++ b/2441/CH5/EX5.17/Ex5_17.sce
@@ -0,0 +1,11 @@
+//Example 5.17
+clc;clear;close;
+format('v',6);
+KE=1500;//MJ
+Pin=5;//MW
+f=50;//Hz
+t=1;//sec
+delKE=Pin*t;////MJ///due to power inputs
+fnew=sqrt((KE+delKE)/KE) *f;//Hz
+delf=fnew-f;///Hz/second
+disp(delf,"Frequency increase rate(Hz/sec)");
diff --git a/2441/CH5/EX5.18/Ex5_18.sce b/2441/CH5/EX5.18/Ex5_18.sce
new file mode 100755
index 000000000..ae6e885b3
--- /dev/null
+++ b/2441/CH5/EX5.18/Ex5_18.sce
@@ -0,0 +1,18 @@
+//Example 5.18
+clc;clear;close;
+format('v',6);
+C=2000;//MW///Capacity
+L=1000;//MW//Load
+H=5;//kWs/KVA
+R=2.4;//Hz/puMW//Regulation
+f=50;//Hz
+delL=1;//%////change in load
+delf=1;//%////change in frequency
+D=delL/delf*L/f;//MW/Hz
+D=D/C;//p.u.MW/Hz
+Kp=1/D;//Hz/p.u.MW
+Tp=2*H/f/D;//sec
+disp("Primary ALFC loop parameters are : ");
+disp(D,"D(p.u.MW/Hz)");
+disp(Kp,"Kp(Hz/p.u.MW)");
+disp(Tp,"Tp(sec)");
diff --git a/2441/CH5/EX5.19/Ex5_19.sce b/2441/CH5/EX5.19/Ex5_19.sce
new file mode 100755
index 000000000..019ba7a74
--- /dev/null
+++ b/2441/CH5/EX5.19/Ex5_19.sce
@@ -0,0 +1,16 @@
+//Example 5.19
+clc;clear;close;
+format('v',6);
+Tp=10;//sec
+Tg=0;//sec
+Tt=0;//sec
+Kp=100;//Hz/p.u.MW
+R=3;///Hz/CuMW
+delPD=0.1;//p.u.
+Ki=0.1;//constant
+f=50;//Hz
+s=poly(0,'s');
+delFs=-Kp/Tp*[delPD/(s^2+s*{(1+Kp/R)/Tp})+Ki*Kp/Tp];
+n=1;//cycle
+time_error=n/f;//sec
+disp(time_error,"Total time error(sec)");
diff --git a/2441/CH5/EX5.2/Ex5_2.sce b/2441/CH5/EX5.2/Ex5_2.sce
new file mode 100755
index 000000000..6dd02222b
--- /dev/null
+++ b/2441/CH5/EX5.2/Ex5_2.sce
@@ -0,0 +1,15 @@
+//Example 5.2
+clc;clear;close;
+format('v',6);
+P=100;//MVA
+f=50;//Hz
+H=5;//kW-sec/kVA(Constant)
+delP=50;//MW(Increased Load)
+td=0.5;//s(Time delay)
+P=P/1000;//kVA
+KE=P*H;//kW-sec
+delP=delP/1000;//kW(Increased Load)
+KE_loss=delP*td;//kW-s
+f_new=sqrt((KE-KE_loss)/KE)*f;//Hz
+f_dev=(f-f_new)/f*100;//%(Frequency deviation)
+disp(f_dev,"Frequency deviation(%)");
diff --git a/2441/CH5/EX5.20/Ex5_20.sce b/2441/CH5/EX5.20/Ex5_20.sce
new file mode 100755
index 000000000..bf7d7895d
--- /dev/null
+++ b/2441/CH5/EX5.20/Ex5_20.sce
@@ -0,0 +1,21 @@
+//Example 5.20
+clc;clear;close;
+format('v',6);
+L=14;//MW//Total Load
+C1=15;//MW
+R1=3;//%//speed regulation
+C2=4;//MW
+R2=4;//%//speed regulation
+LB=4;//MW//Load on bus bar
+LA=10;//MW///Load on bus bar
+f=50;//Hz
+//Load on station A= L1 MW
+//Load on station B= L-L1 MW
+//f-C1*f/100*(L1/C1)=f-C2*f/100*(L2/C2)
+L1=R2*L/C2/(R1/C1+R2/C2);//MW
+L2=L-L1;//MW
+disp(L1,"Load generation at station A(MW)");
+disp(L2,"Load generation at station B(MW)");
+Pt=L1-LA;//MW//Power transmitted A to B
+f_oper=f-R1/100/C1*(L1)*f;//Hz
+disp(f_oper,"Operating Frequency(Hz)");
diff --git a/2441/CH5/EX5.21/Ex5_21.sce b/2441/CH5/EX5.21/Ex5_21.sce
new file mode 100755
index 000000000..91d72357b
--- /dev/null
+++ b/2441/CH5/EX5.21/Ex5_21.sce
@@ -0,0 +1,15 @@
+//Example 5.21
+clc;clear;close;
+format('v',6);
+C1=300;//MW
+C2=400;//MW
+G1=4;//%//droop characteristics of governer
+G2=6;//%//droop characteristics of governer
+L=400;//MW
+f=50;//Hz
+L1=C1*L/(C1+C2);//MW//Load on 300 MW generator
+L2=L*C2/(C1+C2);//MW//Load on 400 MW generator
+f01=f*(C1)/(C1-G1/100*L1);//Hz///No load frequency
+disp(f01,"No load frequency of 300 MW generator(Hz)");
+f02=f*(C2)/(C2-G2/100*L2);//Hz///No load frequency
+disp(f02,"No load frequency of 400 MW generator(Hz)");
diff --git a/2441/CH5/EX5.22/Ex5_22.sce b/2441/CH5/EX5.22/Ex5_22.sce
new file mode 100755
index 000000000..f5e545c14
--- /dev/null
+++ b/2441/CH5/EX5.22/Ex5_22.sce
@@ -0,0 +1,20 @@
+//Example 5.22
+clc;clear;close;
+format('v',6);
+C1=200;//MW
+C2=100;//MW
+R1=1.5;//%//speed regulation
+R2=3;//%//speed regulation
+L=100;//MW///Load on each bus
+f=50;//Hz
+RA=R1/100*f/C1;//Hz/MW
+RB=R2/100*f/C2;//Hz/MW
+//Let PA= generation at plant A
+//PB=2*L-PA will be generation at plant B
+//RA*PA=RB*PB
+PA=RB*2*L/(RA+RB);//MW
+PB=2*L-PA;//MW
+disp(PA,"Load generation at plant A(MW)");
+disp(PB,"Load generation at plant B(MW)");
+Pt=PA-L;//MW///Power transfer
+disp(Pt,"Power transfer from A to B(MW)");
diff --git a/2441/CH5/EX5.23/Ex5_23.sce b/2441/CH5/EX5.23/Ex5_23.sce
new file mode 100755
index 000000000..067acca6d
--- /dev/null
+++ b/2441/CH5/EX5.23/Ex5_23.sce
@@ -0,0 +1,13 @@
+//Example 5.23
+clc;clear;close;
+format('v',7);
+Z=1.5+%i*2.5;//ohm
+V=11;//kV
+P=20;//MW
+pf=0.8;//power factor
+theta=acosd(pf);
+I=P*1000/sqrt(3)/V/pf;//
+I=I*expm(%i*-theta*%pi/180);//A
+Vdrop=I*Z;//V
+Vboost=Vdrop;//V
+disp(Vboost,"Voltage boost needed at station A(V)");
diff --git a/2441/CH5/EX5.24/Ex5_24.sce b/2441/CH5/EX5.24/Ex5_24.sce
new file mode 100755
index 000000000..ac0085ca4
--- /dev/null
+++ b/2441/CH5/EX5.24/Ex5_24.sce
@@ -0,0 +1,28 @@
+//Example 5.24
+clc;clear;close;
+format('v',6);
+Z=3+%i*9;//%///impedence
+Z=Z/100;//p.u.///Impedence
+I=1;//p.u.
+IZ=Z;//p.u.
+disp("Part(a)");
+//2*I^2-2*cos(del)=[abs(IZ)]^2
+cos_del=acosd((2*I^2-[abs(IZ)]^2)/2);//degree
+disp(cos_del,"Phase angle between two station(degree)");
+angle_abc=87.277;///degree
+theta=180-angle_abc-atand(imag(IZ)/real(IZ));//degree
+Preal=I^2*cosd(theta);//p.u.
+disp(Preal,"Real power transfer(p.u.)");
+Preactive=I^2*sind(theta);//p.u.
+disp(Preactive,"Reactive power transfer(p.u.)");
+disp("Part(b)");
+//1.05^2+1^2-2*1.05*cos(del)=[abs(IZ)]^2
+cos_del=acosd((1.05^2+1^2-[abs(IZ)]^2)/2/1.05);//degree
+disp(cos_del,"Phase angle between two station(degree)");
+angle_dbc=60.53;///degree
+theta=atand(imag(IZ)/real(IZ))-angle_dbc//degree
+Preal=I^2*cosd(theta);//p.u.
+disp(Preal,"Real power transfer(p.u.)");
+Preactive=I^2*sind(theta);//p.u.
+disp(Preactive,"Reactive power transfer(p.u.)");
+//Answer in the textbook is not accurate.
diff --git a/2441/CH5/EX5.3/Ex5_3.sce b/2441/CH5/EX5.3/Ex5_3.sce
new file mode 100755
index 000000000..a8f494730
--- /dev/null
+++ b/2441/CH5/EX5.3/Ex5_3.sce
@@ -0,0 +1,15 @@
+//Example 5.3
+clc;clear;close;
+format('v',7);
+P1=500;//MW
+P2=200;//MW
+f=50;//Hz
+delP=140;//MW(System load increase)
+f_new=49.5;//Hz(Frequency after drop)
+delP1=delP*P1/(P1+P2);//MW
+delP2=delP*P2/(P1+P2);//MW
+f_dev=f_new-f;//Hz
+//For delPdash=0, R1 &R2 can be calculated as :
+R1=-1/delP1*f_dev;//Hz/MW
+R2=-1/delP2*f_dev;//Hz/MW
+disp(R2,R1,"Value of R for unit 1 & 2(Hz/MW)");
diff --git a/2441/CH5/EX5.4/Ex5_4.sce b/2441/CH5/EX5.4/Ex5_4.sce
new file mode 100755
index 000000000..bb242095e
--- /dev/null
+++ b/2441/CH5/EX5.4/Ex5_4.sce
@@ -0,0 +1,23 @@
+//Example 5.4
+clc;clear;close;
+format('v',8);
+f=50;//Hz
+R=2;//Hz/pu MW
+Pr=10000;//MW(Rated Capacity)
+P=Pr/2;//MW(Operating Power)
+delP=2;//%(Load Increase)
+del_f=f*1/100;//Hz(1% change in frequency)
+del_PL=P*1/100;//MW(1% change in load)
+//Rate of change of load with frequency :
+D=del_PL/del_f;//MW/Hz
+D=D/Pr;//p.u. MW/Hz
+//Frequency response characteristic : 
+Beta=D+1/R;//p.u. MW/Hz
+M=delP/100*P;//MW
+M=M/Pr;//p.u. MW
+del_fo=-M/Beta;//Hz
+disp(del_fo,"Static frequency drop(Hz)")
+R=%inf;
+Beta=D+1/R;//p.u. MW/Hz
+del_fo=-M/Beta;//Hz
+disp(del_fo,"If speed governer loop is open, frequency drop(Hz)")
diff --git a/2441/CH5/EX5.5/Ex5_5.sce b/2441/CH5/EX5.5/Ex5_5.sce
new file mode 100755
index 000000000..ec2d952b6
--- /dev/null
+++ b/2441/CH5/EX5.5/Ex5_5.sce
@@ -0,0 +1,18 @@
+//Example 5.5
+clc;clear;close;
+format('v',7);
+C=10000;//MW(Control area capacity)
+P=5000;//MW
+H=5;//s
+R=3;//Hz/pu MW
+f=50;//Hz
+del_f=f*1/100;//Hz
+del_PL=P*1/100;//MW
+D=del_PL/del_f;//MW/Hz
+D=D/C;//p.u. MW/Hz
+//Primary ALFC loop parameters :
+Kp=1/D;//Hz/p.u. MW
+Tp=2*H/f/D;//s
+disp("Primary ALFC loop parameters :")
+disp(Kp,"Kp(Hz/p.u. MW)");
+disp(Tp,"Tp(seconds)");
diff --git a/2441/CH5/EX5.6/Ex5_6.sce b/2441/CH5/EX5.6/Ex5_6.sce
new file mode 100755
index 000000000..1f0349eac
--- /dev/null
+++ b/2441/CH5/EX5.6/Ex5_6.sce
@@ -0,0 +1,24 @@
+//Example 5.6
+clc;clear;close;
+format('v',6);
+f=50;//Hz
+R=2;//Hz/pu MW
+Pr=10000;//MW(Rated Capacity)
+P=Pr/2;//MW(Operating Power)
+delP=2;//%(Load Increase)
+del_f=f*1/100;//Hz(1% change in frequency)
+del_PL=P*1/100;//MW(1% change in load)
+//Rate of change of load with frequency :
+D=del_PL/del_f;//MW/Hz
+D=D/Pr;//p.u. MW/Hz
+//Frequency response characteristic : 
+Beta=D+1/R;//p.u. MW/Hz
+M=delP/100*P;//MW
+M=M/Pr;//p.u. MW
+del_fo=-M/Beta;//Hz
+disp("Frequency drop contribution to increase in load(MW) : ");
+delP_fo=-del_fo*(D*Pr);//MW
+disp(delP_fo);
+disp("Increase in generation to meet the increase load(MW) ");
+delP_gen=-del_fo/R*Pr;//MW
+disp(delP_gen);
diff --git a/2441/CH5/EX5.7/Ex5_7.sce b/2441/CH5/EX5.7/Ex5_7.sce
new file mode 100755
index 000000000..3bb22d797
--- /dev/null
+++ b/2441/CH5/EX5.7/Ex5_7.sce
@@ -0,0 +1,15 @@
+//Example 5.7
+clc;clear;close;
+format('v',5);
+G=100;//MVA
+f=50;//Hz
+n=3000;//rpm
+L=25;//MW//Load
+td=0.5;//sec
+H=4.5;//MW-sec/MVA
+//Calculation
+KE=H*G;//MW-sec////at no load
+KE_Loss=L*td;//MW-sec///due to increase in load
+f_new=sqrt((KE-KE_Loss)/KE)*f;//Hz
+delF=(f-f_new)/f*100;//%////frequency deviation
+disp(delF,"Frequency deviation(%)");
diff --git a/2441/CH5/EX5.8/Ex5_8.sce b/2441/CH5/EX5.8/Ex5_8.sce
new file mode 100755
index 000000000..9194eff06
--- /dev/null
+++ b/2441/CH5/EX5.8/Ex5_8.sce
@@ -0,0 +1,24 @@
+//Example 5.8
+clc;clear;close;
+format('v',6);
+C=4000;//MW
+f=50;//Hz
+L=2500;//MW//Load
+R=2;//Hz/p.u.MW////Speed regulation constant
+H=5;//sec////Inertia constant
+delPL=2;//%////change in load
+delf=1;//%////change in frequency
+disp("Part(a)");
+D=delPL/delf*L/f;//MW/Hz
+D=D/C;//p.u.MW/Hz
+Beta=D+1/R;//p.u.MW/Hz
+delf0=-0.2;//Hz
+M=-(delf0)*Beta;//p.u.MW
+M=M*C;//MW
+disp(M,"Largest change in step load(MW)");
+disp("Part(b)");
+Kp=1/D;//Hz/p.u.MW
+Tp=2*H/f/D;//sec
+Tdash=(R+Kp)/R/Tp;//sec
+disp(Tdash,"(R+Kp)/(R*Tp) in seconds = ");
+printf('Change in frequency as a funtion of time, \ndelf(t) = -0.2*(1-epsilon^(-%f*t))',Tdash);
diff --git a/2441/CH5/EX5.9/Ex5_9.sce b/2441/CH5/EX5.9/Ex5_9.sce
new file mode 100755
index 000000000..ef755fe69
--- /dev/null
+++ b/2441/CH5/EX5.9/Ex5_9.sce
@@ -0,0 +1,27 @@
+//Example 5.9
+clc;clear;close;
+format('v',7);
+C=4000;//MW
+f=50;//Hz
+L=C;//MW//Load
+R=2.5;//%////Speed regulation constant
+H=5;//sec////Inertia constant
+delPL=1;//%////change in load
+delf=1;//%////change in frequency
+disp("Part(a)");
+Ls=80;//MW;//increase in step to load
+R=R/100*f;//z/p.u.MW
+D=delPL/delf*L/f;//MW/Hz
+D=D/C;//p.u.MW/Hz
+M=Ls/L;//unitless//for given step load
+Kp=1/D;//Hz/p.u.MW
+Tp=2*H/f/D;//sec
+Tdash1=(R+Kp)/R/Tp;//sec
+disp(Tdash1,"(R+Kp)/(R*Tp) in seconds = ");
+Tdash2=(R*Kp*M)/(R+Kp);//sec
+disp(Tdash2,"(R*Kp*M)/(R+Kp) in seconds = ");
+delf0=-Tdash2;//Hz////Static frequency error
+disp(delf0,"Static frequency error(Hz)");
+disp("Part(b)");
+Ki=(1+Kp/R)^2/4/Tp/Kp;//p.u.MW/Hz
+disp(Ki,"Critical value of Ki(p.u.MW/Hz)");
-- 
cgit