comparing element of matrix

1 次查看（过去 30 天）

Mariam Gasra 2019-3-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/453336-comparing-element-of-matrix

评论： Jan 2019-5-3

clc;
unit=input('input the number of unit:');
D=input('input total load :'); dP=D;
Bdat1=input('transmission losses considered matrix:');
B=Bdat1;
dB=diag(B);
data1=input('input the data:n  a   b   c     min  max');
DA=array2table(data1,'V',{'Unit' 'a' 'b' 'c' 'Pl' 'Ph'});
x=max(DA.b);
n=input('insert number of iteration:')
for i=1:n
while abs(dP)>0.00001
    P=(x-DA.b)./(2*(DA.c+x*dB));
    P=min(P,DA.Ph); P=max(P,DA.Pl);
    dP=D+P'*B*P-sum(P);
    x=x+dP*2/(sum(1./DA.c));
end
end
P((P> DA.Pl) | (P < DA.Ph)) = 0;
C=DA.a+DA.b.*P+DA.c.*P.*P;
totalCost=sum(C); 
display(totalCost);
lamda=x; display(lamda); 
Loss=P'*B*P; display(Loss);
display(P);
display(C);

%% Economic load dispatch of Real power
% D  is total Load
% dP is change during iterations
% B  is loss coefficients matrix
% dB is loss vector
% DA is data of units a Table
% x  is lambda,
% P  is Real Power vector of all units
% C  is cost of production individual units
% totalCost for all units
% Loss in transmission
% Pl Ph are min and max Constraints