blob: 6bf46aba17b3c242bd9a9c5458ebee8c5db082f3 (
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
42
43
44
45
46
|
//Book - Power System: Analysis & Design 5th Edition
//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye
//Chapter - 12 ; Example 12.7
//Scilab Version - 6.0.0 ; OS - Windows
clc;
clear;
C1=[8e-3 10 0] //Coefficients of cost equation for unit 1
C2=[9e-3 8 0] //Coefficients of cost equation for unit 2
dC1=[2*C1(1) C1(2)] //Coefficients of incremental cost equation for unit 1
dC2=[2*C2(1) C2(2)] //Coefficients of incremental cost equation for unit 2
P1lim=[100 600]; //Lower and upper generation limit for unit 1
P2lim=[400 1000]; //Lower and upper generation limit for unit 2
result=[];
for PT=[500 600 700 725 800 900 1000 1100 1200 1244 1300 1400 1500]
P1=(dC2(1)*PT+(dC2(2)-dC1(2)))/(dC2(1)+dC1(1));
P2=PT-P1;
dC1value=dC1(1)*P1+dC1(2);
dC2value=dC2(1)*P2+dC2(2);
if P1<P1lim(1) | P1>P1lim(2) //Checking for limits of P1
if P1<P1lim(1)
P1=P1lim(1)
else
P1=P1lim(2)
end
P2=PT-P1;
dC1value=dC2(1)*P2+dC2(2);
elseif P2<P2lim(1) | P2>P2lim(2) //Checking for limits of P2
if P2<P2lim(1)
P2=P2lim(1)
else
P2=P2lim(2)
end
P1=PT-P2;
dC1value=dC1(1)*P1+dC1(2);
end
CT=C1(1)*P1^2+C1(2)*P1+C1(3)+C2(1)*P2^2+C2(2)*P2+C2(3); //Total cost in $/hr
result=[result;PT P1 P2 dC1value CT]
end
disp(result,' PT(MW) P1(MW) P2(MW) dC/dP($/MWhr) CT($/hr)');
|
