Scribd
Upload
19 views

WITHLOSSS1

The document describes solving an economic load dispatch problem with losses and no generator limits using MATLAB. The problem involves minimizing total generation cost to meet a load of 640.82 MW by optimally scheduling generation across two power plants, taking into account each plant's cost function and transmission losses.

Uploaded by

Kamal Chapagain
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
19 views

WITHLOSSS1

The document describes solving an economic load dispatch problem with losses and no generator limits using MATLAB. The problem involves minimizing total generation cost to meet a load of 640.82 MW by optimally scheduling generation across two power plants, taking into account each plant's cost function and transmission losses.

Uploaded by

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

EXAMPLE: The fuel cost functions in Rs.

/hr for TWO


thermal plants are given by
C1=420+9.2 P +0.004 P
1
2
1

C2 =350+8.5 P +0.0029 P
2
2
2

Where P1, P2, are in MW. Determine the optimal


scheduling of generation if the total load is
640.82MW. Estimate value of LAMBDA =12
Rs./MWh. The transmission power loss is given by
the expression
PL=0.0346 P +0.00643 P
2
1( PU )
2
2 (PU )

Solution:
%MATLAB PROGRAM FOR ECONOMIC LOAD DISPATCH WITH LOSSES AND NO GENERATOR LIMITS
clc;
clear;
% uno d b a
costdata=[1 420 9.2 0.004;
2 350 8.5 0.0029];

Ng=length(costdata(:,1));
for i=1:Ng
uno(i)=costdata(i,1);
d(i)=costdata(i,2);
b(i)=costdata(i,3);
a(i)=costdata(i,4);
end
lambda=12;
pd=640.82;
delp=0.1;
dellambda=0;
lossdata=[0.0346 0.00643];
totgencost=0;
for i=1:Ng
B(i)=lossdata(1,i);
end
iter=0;
disp('iter lambda pg1 pg2 ploss');
while(abs(delp)>=0.001)
iter=iter+1;
lambda=lambda+dellambda;
pl=0;
sum=0;
delpla=0;
for i=1:Ng
den=2*(a(i)+lambda*B(i)*0.01);
p(i)=(lambda-b(i))/den;
pl=pl+(B(i)*0.01*p(i)*p(i));
sum=sum+p(i);
end
delp=pd+pl-sum;
for i=1:Ng
den=2*(a(i)+lambda*B(i)*0.01)^2;
delpla=delpla+(a(i)+B(i)*0.01*b(i))/den;
end
dellambda=delp/delpla;
[iter;lambda;p(1);p(2);pl]'
end
for i=1:Ng
den=1-(B(i)*p(i)*2*0.01);
l(i)=1/den;
end
totgencost=0;
for i=1:Ng
totgencost=totgencost+(d(i)+b(i)*p(i)+a(i)*p(i)*p(i));
ifc(i)=2*a(i)*p(i)+b(i);
end
disp('FINAL OUTPUT OF MATLAB PROGRAM dispatch3.m');
lambda
disp('GENERATING UNIT OPTIMAL GENRATION(MW)');
[uno; p]'
disp('INCREMENTAL FUEL COST(Rs./MWhr.) AND PENALTY FACTOR ARE');
disp('UNIT NO. IFC L');
[uno; ifc; l]'
disp('CHECK LAMBDA=IFC*L');
disp('UNIT NO. LAMBDA');
[uno; ifc.*l]'
disp('TOTAL GENERATION COST(Rs./hr.)');
totgencost

You might also like