EXP 8: ECONOMIC DISPATCH IN POWER
SYSTEMS
AIM:
To develop a MATLAB program to solve the economic dispatch problem for
a given system.
APPARATUS:
[Link] Software Version
1. MATLAB R2022a
ALGORITHM:
LAMBDA-ITERATION METHOD:
• Set value of λ initially.
• Solve the below equation to obtain power generated by the units, Pi:
𝐝𝐅𝐢 (𝐏𝐢 )
=𝛌
𝐏𝐢
Where 𝐹𝑖 (𝑃𝑖 ) is the fuel cost function for units 𝑖 = 1 𝑡𝑜 𝑁; this is the function
that must be minimized.
• If PLoad is the load demand that must be satisfied, the difference between
this demand the and sum of the generated powers is found.
𝐍
∈ = 𝐏𝐥𝐨𝐚𝐝 − ∑ 𝐏𝐢
𝐢=𝟏
• If ∈ > 0, then load is greater than the generation. Therefore, increase the λ
value and repeat steps 2-3.
• If ∈ < 0, then generation is greater than load. Therefore, decrease the λ
value and repeat steps 2-3.
• Stop the iteration when ∈ falls below a particular threshold.
• This process is illustrated below as a flow chart.
MATLAB CODE:
%% Economic load dispatch of Real power
% NOT considering transmission losses
clear all;
clc;
% n a b*cost c*cost min max
d = [1 225*0.8 8.4*0.8 0.0025*0.8 45 350;
2 729*1.02 6.3*1.02 0.0081*1.02 45 350;
3 400*0.9 7.5*0.9 0.0025*0.9 47.5 450];
D = 450;
b = d(:,3);
c = d(:,4);
Pmin = d(:,5);
Pmax = d(:,6);
dP = D;
x = 6; % assume lambda
while abs(dP)>0.00001
P = (x-b)./(2*c);
P = min(P,Pmax);
P = max(P,Pmin);
dP = D - sum(P);
x = x + dP*0.0025;
end
fprintf("Final value of lambda is: %2.3f R/MWh", x)
Power = table(d(:,1),P,'V',{'Unit' 'Power'})
OUTPUT:
Final value of lambda is: 7.544 R/MWh
Power =
3×2 table
Unit Power
____ ______
1 205.95
2 67.647
3 176.4
RESULT:
Thus, the economic dispatch problem was solved using MATLAB and the
results were verified.