initial commit / add all books

24 files changed, 472 insertions, 0 deletions

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)");