blob: ec18124240a03eb35db34f210b2f425178f89133 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
|
//Electric Power Generation, Transmission and Distribution by S.N.Singh
//Publisher:PHI Learning Private Limited
//Year: 2012 ; Edition - 2
//Example 4.6
//Scilab Version : 6.0.0 ; OS : Windows
clc;
clear;
k=0.6;
t=1.503032/0.6; //Time in hours
Df=3; //Density factor
P=30- 8*sin(k*t)+0.325*t; //Load variation at a power supply station in MW
i=1;
n=1;
while (t(i)<24)
t(i+1)=(2*n*%pi-1.503032)/0.6;
t(i+2)=(2*n*%pi+1.503032)/0.6;
if (t(i+1)<24)&(t(i+2)<24) then
i=i+2;
else
t(i+1)=25;
i=i+1;
end
n=n+1;
end
P=30- 8*sin(k*t)+0.325*t;
Max_demand=max(P); //Maximum demand on the system in MW
Avg_load=(1/24)*(30*24+(8/0.6)*(cosd(0.6*24)-cosd(0.6*0))+0.325*24^(2)/2); //Applying integration for power equation
Lf=Avg_load/Max_demand; //Load factor of the system
Total_load=Max_demand*Df; //Total installed load of the system in MW
printf("\nMaximum demand on the system is %.3f MW",Max_demand);
printf("\nLoad factor of the system %.3f",Lf);
printf("\nTotal installed load is %.3f MW",Total_load);
//different t values cannot be obtained for section d so cannot be solved in scilab
//Variation present in result due to wrong calculation of Avg load value in the book
|
