diff options
Diffstat (limited to '3872/CH12/EX12.7/Ex12_7.sce')
| -rw-r--r-- | 3872/CH12/EX12.7/Ex12_7.sce | 46 |
1 files changed, 46 insertions, 0 deletions
diff --git a/3872/CH12/EX12.7/Ex12_7.sce b/3872/CH12/EX12.7/Ex12_7.sce new file mode 100644 index 000000000..6bf46aba1 --- /dev/null +++ b/3872/CH12/EX12.7/Ex12_7.sce @@ -0,0 +1,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)');
|