Experiment 3 Economic Load Dispatch
Experiment 3 Economic Load Dispatch
Experiment 3 Economic Load Dispatch
Experiment No.-03
Aim: - To obtain economic load dispatch for a general power system using
the MATLAB program.
Since transmission losses are neglected, the total demand PD is equal to the sum of
all power generation. A cost function Ci is assumed to be known for each plant. The
problem is to find the real power generation for each plant such that the objective
function (i.e: total production cost) as defined by the equation
ng
Ct =∑ Ci
i=1
n
=∑ α i + β i Pi + γ i P2i
i=1 (4-1)
∑ Pi =PD
i =1 (4-2)
where Ct is the total production cost, Ci is the production cost of ith plant, Pi is the
A typical approach is to augment the constraints into objective function by using the
Lagrange multiplier
( )
ng
L=C t + λ P D−∑ P i
i=1 (4-3)
The minimum of L is found by equating the partials derivatives of the function with
respect to Pi and λ to zero.
∂L
=0
∂ Pi (4-4)
∂L
=0
∂λ (4-5)
Since
Ct =C1 +C 2 +. . .. .+C n
g
then
∂C t ∂ Ci
= =λ
∂ Pi ∂ P i
i=1,….,ng (4-6)
or
β i +2 γ i Pi =λ (4-7)
∑ Pi =PD
i =1 (4-8)
Equation (4-8) is precisely the power flow equality constraint. In summary, when
losses are neglected with no generator limits, for most economic operation, all plants
must operate at equal increment production cost, while satisfying the equality
constraint given by (4-8). In order to find the solution, (4-7) is solved for Pi as
λ−β i
Pi=
2γi i=1,….,ng (4.9)
The relations given by (4.9) are known as the coordination equations. They are
functions of λ. An analytical solution can be obtained for λ by substituting (4.9) in (4-
8), i.e
ng
λ−β i
∑ 2 γi
=PD
i=1 (4-10)
or
ng
βi
( PD + ∑ )
i=1 2γi
λ= ng
1
∑ 2γ
i=1 i (4-11)
The value of λ found from (4-11) is substituted in (4-7) to obtain the optimal
scheduling of generation. The solution for economic dispatch neglecting losses has
been found analytically.
f(λ) = PD (4-12)
Expanding the left-hand side of the above equation in Taylor’s series about an
operating point λ(k), and neglecting the higher order terms results in
(k) df ( λ ) ( k ) ( k )
f ( λ) +( ) Δλ =P D
dλ (4-13)
or
(k)
ΔP
Δλ( k )=
df ( λ ) ( k )
( )
dλ
ΔP( k )
¿ ∂ Pi
∑ ( ∂ λ )( k)
(4-14)
or
ΔP( k )
(k )
Δλ =
∑ 21γ
i (4-15)
and therefore,
where
ng
ΔP =P D−∑ P ( k )
(k)
i=1 i
(4-17)
PROCEDURE:-
Example:-
The Fuel cost function for three thermal power plants in Rs/hr is given by,
C1=500+5.3P1+0.004P2
C2=400+5.5P1+0.006P2
C3=200+5.8P1+0.009P2
Total power PD=800 neglecting line losses and no generator limit. Find optimal
dispatch and total cost in Rs/hr.
Manual Calculation:-
a) by analytical method using (4-11)
b) by graphical demonstration
Substituting for λ in the coordination equation given in (2-9), the optimal dispatch
is
8 .5−5. 3
P1= =400 MW
2∗0 . 004
8 .5−5 .5
P2 = =250 MW
2∗0 . 006
8 .5−5 .8
P3 = =150 MW
2∗0 . 009
b) From equation (4-6), the necessary conditions for optimal dispatch are
dC 1
=5 .3+ 0. 004 P1 =λ
dP1
dC 2
=5 . 5+0 . 012 P2 =λ
dP 2
dC 3
=5 . 3+0 . 018 P3 = λ
dP 3
subject to
P1+P2+P3=PD
To demonstrate the concept of equal increment cost for optimal dispatch, we can use
MATLAB plot command to plot the incremental cost of each plant on the same graph
as shown in Figure 4-2. To obtain a solution, various values of λ could be tried until
otherwise, if ∑ Pi >P D , we reduce λ. Therefore, the horizontal line shown in the graph
Solution:
% Iterative solution Using Newton method
alpha =[500; 400; 200];
beta = [5.3; 5.5; 5.8]; gama=[.004; .006; .009];
PD=800;
DelP = 10; % Error in DelP is set to a high value
lambda = input('Enter estimated value of Lambda = ');
fprintf('\n ')
disp([' Lambda P1 P2 P3 DP'...
' grad Delambda'])
iter = 0; % Iteration counter
while abs(DelP) >= 0.001 % Test for convergence
iter = iter + 1; % No. of iterations
P = (lambda - beta)./(2*gama);
DelP =PD - sum(P); % Residual
J = sum( ones(length(gama), 1)./(2*gama)); % Gradient sum
Delambda = DelP/J; % Change in variable
disp([lambda, P(1), P(2), P(3), DelP, J, Delambda])
lambda = lambda + Delambda; % Successive solution
end
totalcost = sum(alpha + beta.*P + gama.*P.^2)
Result:
Conclusion: