diff options
| author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
| commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
| tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3872/CH12/EX12.7/Ex12_7.sce | |
| parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
| download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip | |
initial commit / add all books
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)');
|