Thanks to visit codestin.com
Credit goes to www.scribd.com

Scribd
Upload
63 views27 pages

Assignment#08 Problem 3.10 From Allen J. Wood Book (2 Edition)

This document contains the results of a power system operation and control problem solved using MATLAB. It presents the parameters and equations for three generating units, runs an iterative optimization algorithm over 5 iterations to determine the optimal dispatch of the units to minimize total generation cost while meeting a load of 900 MW, and shows the outputs of each iteration. The optimal dispatch after 5 iterations was found to be 301.4411 MW for unit 1, 301.1275 MW for unit 2, and 297.4315 MW for unit 3.

Uploaded by

Asif Bhati
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
63 views27 pages

Assignment#08 Problem 3.10 From Allen J. Wood Book (2 Edition)

This document contains the results of a power system operation and control problem solved using MATLAB. It presents the parameters and equations for three generating units, runs an iterative optimization algorithm over 5 iterations to determine the optimal dispatch of the units to minimize total generation cost while meeting a load of 900 MW, and shows the outputs of each iteration. The optimal dispatch after 5 iterations was found to be 301.4411 MW for unit 1, 301.1275 MW for unit 2, and 297.4315 MW for unit 3.

Uploaded by

Asif Bhati
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 27

ASSIGNMENT#08

Problem 3.10 from Allen J. Wood book (2nd edition)


Subject:
Power System Operation & Control
Submitted To:
Dr. Shahzad
Submitted By:
M. Mohsin Raza (2017-EE-729)
M. Asif (2017-EE-715)
Abdul Satar (2017-EE-707)
Program:
B.Sc. Electrical Engineering

Electrical Engineering Department


Muhammad Nawaz Sharif
University of Engineering and Technology
Multan, Pakistan
Matlab Results
A1=561
B1=7.92
C1=0.001562
A2=310
B2=7.85
C2=0.00194
A3=78
B3=7.97
C3=0.00482
fprintf('P1=%d+%dP1+%dP1^2+%dP1^3',A1,B1,C1)
fprintf('P2=%d+%dP2+%dP2^2+%dP2^3',A2,B2,C3)
fprintf('P3=%d+%dP3+%dP3^2+%dP3^3',A3,B3,C3)
1=input('Enter of P1')
P2=input('Enter of P2')
P3=input('Enter of P3')
Pl=900
for a=1:5 % a is showing the no of itteration
C(1,a)=((2*C1*P1)+(B1))%cost
C(2,a)=((2*C2*P2)+(B2))%cost
C(3,a)=((2*C3*P3)+(B3))%cost
% For first Itteration
if(a==1)
lamda=(C(1,a)+C(2,a)+C(3,a))*1/3
P(4,1)=lamda
end
if(a>1)
X(1,a-1)=P1
X(2,a-1)=P2
X(3,a-1)=P3
X(4,a-1)=lamda
for i=1:3
L(i,a)=C(i,a)-lamda
end
P(1,a)=X(1,a-1)-L(1,a-1)
P(2,a)=X(2,a-1)-L(2,a-1)
P(3,a)=X(3,a-1)-L(3,a-1)
P(4,a)=Pl-(P1+P2+P3)
P1=P(1,a)
P2=P(2,a)
P3=P(3,a)
lamda=P(4,1)-P(4,a)
end
end

A1 =

561

B1 =

7.9200

C1 =

0.0016

A2 =

310

B2 =

7.8500

C2 =

0.0019

A3 =

78

B3 =

7.9700
C3 =

0.0048

P1=561+7.920000e+00P1+1.562000e-03P1^2+P2=310+7.850000e+00P2+4.820000e
-03P2^2+P3=78+7.970000e+00P3+4.820000e-03P3^2+Enter of P1300

P1 =

300

Enter of P2300

P2 =

300

Enter of P3300

P3 =

300

Pl =

900

C=

8.8572

C=

8.8572
9.0140

C=

8.8572
9.0140
10.8620

lamda =

9.5777

P=

0
0
0
9.5777

C=

8.8572 8.8572
9.0140 0
10.8620 0

C=

8.8572 8.8572
9.0140 9.0140
10.8620 0

C=

8.8572 8.8572
9.0140 9.0140
10.8620 10.8620

X=

300

X=
300
300

X=

300
300
300

X=

300.0000
300.0000
300.0000
9.5777

L=

0 -0.7205

L=

0 -0.7205
0 -0.5637

L=

0 -0.7205
0 -0.5637
0 1.2843

P=

0 300.0000
0 0
0 0
9.5777 0
P=

0 300.0000
0 300.0000
0 0
9.5777 0

P=

0 300.0000
0 300.0000
0 300.0000
9.5777 0

P=

0 300.0000
0 300.0000
0 300.0000
9.5777 0

P1 =

300

P2 =

300

P3 =

300

lamda =

9.5777
C=

8.8572 8.8572 8.8572


9.0140 9.0140 0
10.8620 10.8620 0

C=

8.8572 8.8572 8.8572


9.0140 9.0140 9.0140
10.8620 10.8620 0

C=

8.8572 8.8572 8.8572


9.0140 9.0140 9.0140
10.8620 10.8620 10.8620

X=

300.0000 300.0000
300.0000 0
300.0000 0
9.5777 0

X=

300.0000 300.0000
300.0000 300.0000
300.0000 0
9.5777 0

X=

300.0000 300.0000
300.0000 300.0000
300.0000 300.0000
9.5777 0

X=

300.0000 300.0000
300.0000 300.0000
300.0000 300.0000
9.5777 9.5777

L=

0 -0.7205 -0.7205
0 -0.5637 0
0 1.2843 0

L=

0 -0.7205 -0.7205
0 -0.5637 -0.5637
0 1.2843 0

L=

0 -0.7205 -0.7205
0 -0.5637 -0.5637
0 1.2843 1.2843

P=

0 300.0000 300.7205
0 300.0000 0
0 300.0000 0
9.5777 0 0

P=

0 300.0000 300.7205
0 300.0000 300.5637
0 300.0000 0
9.5777 0 0

P=

0 300.0000 300.7205
0 300.0000 300.5637
0 300.0000 298.7157
9.5777 0 0

P=

0 300.0000 300.7205
0 300.0000 300.5637
0 300.0000 298.7157
9.5777 0 0

P1 =

300.7205

P2 =

300.5637

P3 =

298.7157

lamda =

9.5777

C=

8.8572 8.8572 8.8572 8.8595


9.0140 9.0140 9.0140 0
10.8620 10.8620 10.8620 0

C=

8.8572 8.8572 8.8572 8.8595


9.0140 9.0140 9.0140 9.0162
10.8620 10.8620 10.8620 0

C=

8.8572 8.8572 8.8572 8.8595


9.0140 9.0140 9.0140 9.0162
10.8620 10.8620 10.8620 10.8496

X=

300.0000 300.0000 300.7205


300.0000 300.0000 0
300.0000 300.0000 0
9.5777 9.5777 0

X=

300.0000 300.0000 300.7205


300.0000 300.0000 300.5637
300.0000 300.0000 0
9.5777 9.5777 0

X=

300.0000 300.0000 300.7205


300.0000 300.0000 300.5637
300.0000 300.0000 298.7157
9.5777 9.5777 0

X=

300.0000 300.0000 300.7205


300.0000 300.0000 300.5637
300.0000 300.0000 298.7157
9.5777 9.5777 9.5777

L=

0 -0.7205 -0.7205 -0.7183


0 -0.5637 -0.5637 0
0 1.2843 1.2843 0

L=

0 -0.7205 -0.7205 -0.7183


0 -0.5637 -0.5637 -0.5615
0 1.2843 1.2843 0

L=

0 -0.7205 -0.7205 -0.7183


0 -0.5637 -0.5637 -0.5615
0 1.2843 1.2843 1.2719

P=

0 300.0000 300.7205 301.4411


0 300.0000 300.5637 0
0 300.0000 298.7157 0
9.5777 0 0 0

P=

0 300.0000 300.7205 301.4411


0 300.0000 300.5637 301.1275
0 300.0000 298.7157 0
9.5777 0 0 0

P=
0 300.0000 300.7205 301.4411
0 300.0000 300.5637 301.1275
0 300.0000 298.7157 297.4315
9.5777 0 0 0

P=

0 300.0000 300.7205 301.4411


0 300.0000 300.5637 301.1275
0 300.0000 298.7157 297.4315
9.5777 0 0 0

P1 =

301.4411

P2 =

301.1275

P3 =

297.4315

lamda =

9.5777

C=

8.8572 8.8572 8.8572 8.8595 8.8617


9.0140 9.0140 9.0140 9.0162 0
10.8620 10.8620 10.8620 10.8496 0

C=

8.8572 8.8572 8.8572 8.8595 8.8617


9.0140 9.0140 9.0140 9.0162 9.0184
10.8620 10.8620 10.8620 10.8496 0

C=

8.8572 8.8572 8.8572 8.8595 8.8617


9.0140 9.0140 9.0140 9.0162 9.0184
10.8620 10.8620 10.8620 10.8496 10.8372

X=

300.0000 300.0000 300.7205 301.4411


300.0000 300.0000 300.5637 0
300.0000 300.0000 298.7157 0
9.5777 9.5777 9.5777 0

X=

300.0000 300.0000 300.7205 301.4411


300.0000 300.0000 300.5637 301.1275
300.0000 300.0000 298.7157 0
9.5777 9.5777 9.5777 0

X=

300.0000 300.0000 300.7205 301.4411


300.0000 300.0000 300.5637 301.1275
300.0000 300.0000 298.7157 297.4315
9.5777 9.5777 9.5777 0

X=

300.0000 300.0000 300.7205 301.4411


300.0000 300.0000 300.5637 301.1275
300.0000 300.0000 298.7157 297.4315
9.5777 9.5777 9.5777 9.5777

L=
0 -0.7205 -0.7205 -0.7183 -0.7160
0 -0.5637 -0.5637 -0.5615 0
0 1.2843 1.2843 1.2719 0

L=

0 -0.7205 -0.7205 -0.7183 -0.7160


0 -0.5637 -0.5637 -0.5615 -0.5594
0 1.2843 1.2843 1.2719 0

L=

0 -0.7205 -0.7205 -0.7183 -0.7160


0 -0.5637 -0.5637 -0.5615 -0.5594
0 1.2843 1.2843 1.2719 1.2595

P=

0 300.0000 300.7205 301.4411 302.1593


0 300.0000 300.5637 301.1275 0
0 300.0000 298.7157 297.4315 0
9.5777 0 0 0 0

P=

0 300.0000 300.7205 301.4411 302.1593


0 300.0000 300.5637 301.1275 301.6890
0 300.0000 298.7157 297.4315 0
9.5777 0 0 0 0

P=

0 300.0000 300.7205 301.4411 302.1593


0 300.0000 300.5637 301.1275 301.6890
0 300.0000 298.7157 297.4315 296.1596
9.5777 0 0 0 0
P=

0 300.0000 300.7205 301.4411 302.1593


0 300.0000 300.5637 301.1275 301.6890
0 300.0000 298.7157 297.4315 296.1596
9.5777 0 0 0 0

P1 =

302.1593

P2 =

301.6890

P3 =

296.1596

lamda =

9.5777
Matlab Results
alpha=[561;310;78]
beta=[792/100;785/100;797/100]
gamma=[0.001562;0.00194;0.00482]
Lambda=0;
Pload=900;
Pinitial(1)=300;
Pinitial(2)=300;
Pinitial(3)=300; Psum=sum(Pinitial(1)+Pinitial(2)+Pinitial(3));
%......Formation of Hessian Matrix......
for i=1:3
L(i)=beta(i)+2*gamma(i)*Pinitial(i);
end
dellambda=[L(1)-Lambda;L(2)-Lambda;L(3)-Lambda;Pload-Psum]
G=2*gamma;
H= [G(1) 0 0 -1;0 G(2) 0 -1;0 0 G(3) -1;-1 -1 -1 0]
M= inv(H)
delpPower=-[M*dellambda]
Lambdaoptimal=[delpPower(4)]
Poptimal=[Pinitial(1);Pinitial(2);Pinitial(3);Lambda]+delpPower

alpha =

561
310
78

beta =

7.9200
7.8500
7.9700

gamma =

0.0016
0.0019
0.0048

dellambda =
8.8572
9.0140
10.8620
0

H=

0.0031 0 0 -1.0000
0 0.0039 0 -1.0000
0 0 0.0096 -1.0000
-1.0000 -1.0000 -1.0000 0

M=

169.7646 -121.0452 -48.7194 -0.4697


-121.0452 160.2719 -39.2267 -0.3781
-48.7194 -39.2267 87.9461 -0.1522
-0.4697 -0.3781 -0.1522 -0.0015

delpPower =

116.6526
53.5110
-170.1636
9.2216

Lambdaoptimal =

9.2216

Poptimal =

416.6526
353.5110
129.8364
9.2216

You might also like