diff options
Diffstat (limited to '2441/CH2/EX2.24/Ex2_24.sce')
| -rwxr-xr-x | 2441/CH2/EX2.24/Ex2_24.sce | 24 |
1 files changed, 24 insertions, 0 deletions
diff --git a/2441/CH2/EX2.24/Ex2_24.sce b/2441/CH2/EX2.24/Ex2_24.sce new file mode 100755 index 000000000..7312c3b8c --- /dev/null +++ b/2441/CH2/EX2.24/Ex2_24.sce @@ -0,0 +1,24 @@ +//exa 2.24
+clc;clear;close;
+format('v',7);
+//dC1/dP1=0.02*P1+16;//
+//dC2/dP2=0.04*P2+20;//
+PL=10;//MW
+P1=100;//MW
+lambda=25;//Rs. per MWh
+B11=PL/P1^2;B22=0;B12=0;//Coefficient Loss
+//dF2/dP2*L2=lambda
+P2=poly(0,'P2');
+L2=1;//penalty factor
+eqn=(0.04*P2+20)*L2-lambda;
+P2=roots(eqn);//MW
+//PL=B11*P1^2
+P1=poly(0,'P1');
+dPLbydP1=2*B11*P1;
+L1=1/(1-dPLbydP1);//penalty factor
+eqn=(0.02*P1+16)-lambda/L1
+P1=roots(numer(eqn));//MW
+disp(P2,P1,"Generation P1 & P2 in MW are ")
+PL=B11*P1^2;//MW
+LD=P1-PL+P2;//MW
+disp(LD,"Load Demand(MW) :");
|