Can I get the code for Economic Load Dispatch using Gradient Descent (matlab)?

11 次查看（过去 30 天）

Suzzane M 2020-6-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/553564-can-i-get-the-code-for-economic-load-dispatch-using-gradient-descent-matlab

编辑： mahesh kumar 2020-6-29

I'm not getting how to proceed after this!!

clc
a=[500;400;200];
b=[5.31;5.5;5.8];
g=[0.004;0.005;0.009];
lambda=input('enter vvalue of lambda');
iter=0;
Pd=800;
n=length(a);
delP=5;tcost=0;
tolerance=0.0001;
while abs(delP)>tolerance
    iter=iter+1;
    x=0; 
    for i=1:n
        p(i)=(lambda-b(i))/(2*g(i));
        x=x+(1/(2*g(i)));
        cost=(a(i)+(b(i)*p(i))+(g(i)*p(i)*p(i)));
  

11 个评论
显示 9更早的评论隐藏 9更早的评论

mahesh kumar 2020-6-27

编辑：mahesh kumar 2020-6-27

在 MATLAB Online 中打开

try this::

B  = [.00003 0 0;0 .00009 0 ; 0 0 .00012]; dB=diag(B); load=850;
x=max(b); dP=1;i=0;
z=randperm(10000,length(b))'; P=load*z/sum(z);
while abs(dP)>0.000001
    i=i+1; disp(i);
    dP=load+P'*B*P-sum(P); % loss=P'*B*P;
    x=x+dP*2/(sum(1./a));    
    P=(x-b-2*(B*P-dB.*P))./(2*a+2*x*dB);
%     plot(x,P,'x'); hold on; pause(1);
%     C=c+b.*P+a.*P.*P;
%     plot(x,sum(C)/100,'o');
end