diff options
Diffstat (limited to '3872/CH12/EX12.9')
| -rw-r--r-- | 3872/CH12/EX12.9/Ex12_9.JPG | bin | 0 -> 21453 bytes | |||
| -rw-r--r-- | 3872/CH12/EX12.9/Ex12_9.sce | 30 |
2 files changed, 30 insertions, 0 deletions
diff --git a/3872/CH12/EX12.9/Ex12_9.JPG b/3872/CH12/EX12.9/Ex12_9.JPG Binary files differnew file mode 100644 index 000000000..3a323446a --- /dev/null +++ b/3872/CH12/EX12.9/Ex12_9.JPG diff --git a/3872/CH12/EX12.9/Ex12_9.sce b/3872/CH12/EX12.9/Ex12_9.sce new file mode 100644 index 000000000..6cd8eba6e --- /dev/null +++ b/3872/CH12/EX12.9/Ex12_9.sce @@ -0,0 +1,30 @@ +//Book - Power System: Analysis & Design 5th Edition
+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye
+//Chapter - 12 ; Example 12.9
+//Scilab Version - 6.0.0 ; OS - Windows
+
+clc;
+clear;
+
+B11=1.5e-4; B12=2e-5; B22=3e-5; // Loss coefficients
+lamda=16; //Area incremental cost in $/MWhr
+
+e1=[20.8e-3 32e-5 6]; //Coefficients of incremental operating cost equation1
+e2=[32e-5 18.96e-3 8]; //Coefficients of incremental operating cost equation2
+
+P1=(e2(2)*e1(3)-e1(2)*e2(3))/(e2(2)*e1(1)-e1(2)*e2(1)); //Solution of P1 from incremental cost equations
+P2=(e2(1)*e1(3)-e1(1)*e2(3))/(e1(2)*e2(1)-e2(2)*e1(1)); //Solution of P2 from incremental cost equations
+
+Pl=B11*P1^2+B12*P1*P2+B22*P2^2; //Total losses
+
+Pt=P1+P2-Pl; //Total demand
+
+CT=10*P1+8e-3*P1^2+8*P2+9e-3*P2^2; // Cost equation taken from example 12.6
+
+printf('The output of each unit are given by P1=%d MW and P2=%d MW\n',P1,P2)
+printf('The total transmission loss is %.2f MW\n',Pl)
+printf('The total demand is %.2f MW\n',Pt)
+printf('The total operation cost is %.2f $/hr',CT)
+
+
+
|
