Program Required Theory::: Matlab 7.0
Program Required Theory::: Matlab 7.0
Fundamentals of economic dispatch and problem solving using classical method with and
without line losses
Program Required:
MATLAB 7.0
Theory:
n a power s!stem" with negligible transmission loss and with # number of spinning thermal
generating units the total s!stem load $% at a particular interval can be met b! different sets of
generation schedules
$&'()* " $&+()* " ,,,,,,$&#(-* ./ ) 0 '"+",,..#1
2ut of these #1 set of generation schedules" the s!stem operator has to choose the set of
schedules"
3hich minimi4e the s!stem operating cost" which is essentiall! the sum of the production cost
of all the generating units. This economic dispatch problem is mathematicall! stated as an
optimi4ation problem
Given5 The number of available generating units #" their production cost functions" their
operating limits and the s!stem load $%"
To determine: The set of generation schedules"
$&i / i 0 '"+,,,# (1)
3hich minimi4e the total production cost"
#
Min / FT 0 Fi ($&i * (2)
i0'
and satisfies the power balance constraint
#
0 $&i 6$% 0 0 (3)
i0'
and the operating limits
$&i"min $&i $&i" "ma7 (4)
The units production cost function is usuall! appro7imated b! 8uadratic function
Fi ($&i* 0 ai $&+
i 9 bi $&i 9 ci / i 0 '"+",,.# (5)
where ai " bi and ci are constants
Necessary conditions for the existence of solution to ED problem
The :% problem given b! the e8uations ('* to (;*. B! omitting the ine8ualit! constraints
(;* tentativel!" the reduce :% problem ('*"(+* and (<* ma! be restated as an unconstrained
optimi4ation problem b! augmenting the ob=ective function ('* with the constraint multiplied
b! La&range multiplier" to obtained the La&range function" L as
# #
Min 5 L ($&' ,,..$&# " * 0 Fi($&i* > ? $&i 6 $%@ ()
The necessar! conditions for the e7istence of solution to (A* are given b!
L B $&i 0 0 0 dFi ($&i* B d$&i > / i 0 '" +",,..# (!)
#
L B 0 0 0 $&i 6 $% (")
i0'
The solution to :% problem can be obtained b! solving simultaneousl! the necessar!
conditions
(7* and (C* which state that the economic generation schedules not onl! satisf! the s!stem
power
balance e8uation (C* but also demand that the incremental cost rates of all the units be e8ual be
e8ual to which can be interpreted as Dincremental cost of received powerE.
3hen the ine8ualit! constraints(;* are included in the :% problem the necessar! condition (7*
gets
modified as
dFi ($&i* B d$&i 0 for $&i"min $&i $&i" "ma7
for $&i 0 $&i" "ma7
# for $&i 0 $&i" "mi FFFF(G*
:conomic 1chedule
$&i 0 ( >bi*B +ai / i0'"+,,,,,.# $1%)
ncremental fuel cost
# #
0 $% 9 ( biB+ai * B ('B+ai* (11)
PROCEDURE
'. :nter the command window of the MATLAB.
+. Hreate a new M 6 file b! selecting File > #ew 6 M 6 File
<. T!pe and save the program.
;. :7ecute the program b! either pressing Tools 6 Iun.
J. Kiew the results.
DATA-1
'. The fuel cost functions for three thermal plants in LBh are given b!
H' 0 J00 9 J.< $' 9 0.00; $'
+ / $' in M3
H+ 0 ;00 9 J.J $+ 9 0.00A $+
+ / $+ in M3
H< 0 +00 9J.C $< 9 0.00G $<
+ / $< in M3
The total load " $% is C00M3.#eglecting line losses and generator limits" find the optimal
dispatch and the total cost in LBh b! anal!tical method. Kerif! the result using MATLAB
program.
&'()'*+
alpha = [500; 400; 200];
beta = [5.3; 5.5; 5.8]; gamma = [0.004; 0.006; 0.009];
PD = 800;
DelP = 10;
lamda = input(!nte" e#timated $alue %& 'amda = (;
&p"int&( (
di#p(['amda P1 P2 P3 DP...
g"ad Delamda](
ite" = 0;
)hile ab#(DelP( *= 0.001
ite" = ite" + 1;
P = (lamda , beta(.-(2.gamma(;
DelP = PD , #um(P(;
/ = #um(%ne#(length(gamma(01(.-(2.gamma((;
Delamda = DelP-/;
di#p([lamda0P(1(0P(2(0P(3(0DelP0/0Delamda](
lamda = lamda + Delamda;
end
t%tal1%#t = #um(alpha + beta..P + gamma..P.22(
DATA-2
The fuel cost functions for three thermal plants in LBh are given b!
H' 0 J00 9 J.< $' 9 0.00; $'
+ / $' in M3
H+ 0 ;00 9 J.J $+ 9 0.00A $+
+ / $+ in M3
H< 0 +00 9J.C $< 9 0.00G $<
+ / $< in M3
The total load " $% is G7JM3.
Generation limit
+00 $' ;J0 M3
'J0 $+ <J0 M3
'00 $< ++J M3
Find the optimal dispatch and the total cost in LBh b! anal!tical method. Kerif! the result
using MATLAB program
PROGRA!
1%#t = [500 5.3 0.004
400 5.5 0.006
200 5.8 0.009];
m)limit# = [200 450
150 350
100 225];
Pdt = 935;
di#pat1h
gen1%#t
RE"U#T: