Economic Load Dispatch In Power System With Considering Losses Using Direct Method

ECONOMIC LOAD DISPATCH IN POWER SYSTEM WITH CONSIDERING LOSSES USING DIRECT METHOD

PROGRAM:

clc;

clear all;

warning off;

a=[.658;.586;.789];

b=[5.3;5.5;5.8];

c=[500;400;200];

pd=800;

delp=10;

noofunits=length(a)

lambda=input('enter estimated value of lam=');

fprintf('\n')

disp(['lambda p1 p2 p3 delta data lambda'])

iter=0;

while abs(delp)>=0.001

iter=iter+1;

    p=(lambda-b)./(2*a);

delp=pd-sum(p);

    j=sum(ones(length(a),1)./(2*a));

delambda=delp/j;

disp([lambda,p(1),p(2),p(3),delp,delambda])

lambda=lambda+delambda;

end

lambda

p

totalcost=sum(c+b.*p+a.*p.^2)