Как добавить минимальное временное ограничение для домашней нагрузки в смешанное целочисленное линейное программирование в Matlab? - PullRequest
Новая
       13

Как добавить минимальное временное ограничение для домашней нагрузки в смешанное целочисленное линейное программирование в Matlab?

0 голосов
/ 18 мая 2018

Я написал код для работы приборов в течение минимального времени работы в проблеме MILP в MATLAB.

Код, который я написал, основан на строках «оптимальной работы тепловых установок»."учебное видео по проблеме MILP. Ниже приводится ссылка на него -

https://in.mathworks.com/videos/mixed-integer-linear-programming-in-matlab-91541.html

Ниже приведен мой код - 

powerprob=optimproblem;
E=[3 2 3 5 2 3 ]; % Cost of electricity
P=[0.5 1 0.5]; % power rating of appliances
D=[2 2 3]; % Total duration of operation
% For making the appliances operate within a fixed start and finish time

V=zeros(3,3);
TstartM=[2 3 1]; % Start time of each appliances
TfinishM=[5 5 4 ]; % Finish time
for jj=1:3 
  for kk=1:3 % total appliances
      DurationM=[1 1 1];
    if kk<TstartM(jj) 
       V(kk,jj)=0;
    elseif kk>TfinishM(jj)
       V(kk,jj)=0;
    elseif kk>=TstartM(jj) && kk<=TfinishM(jj)
       V(kk,jj)=DurationM(jj);
    end
  end
end
V
VV=cat(1,V,V) % Start and finish time window for appliances
Pgrid=optimvar('Pgrid',6,1,'LowerBound',0,'UpperBound',7); % Power Imported from the grid
A=optimvar('A',6,3,'Type','integer','LowerBound',0,'UpperBound',1); % ON/Off schedule of the appliances

HLD=((A.*VV)*P'); % Hourly Load

gridcost=sum((Pgrid).*E');

powerprob.Objective=gridcost;
% min up time for the appliances
 powerprob.Constraints.D1=optimconstr(6,3);

for jj=1:3
   for kk=1:6 % based on possible startups; no penalty at end for running over

       MinUpTimeM=[2 1 1];

     if kk>6-MinUpTimeM(jj)
       sumidx=kk:6;
     elseif kk<=6-MinUpTimeM(jj)
       sumidx=kk:kk+MinUpTimeM(jj)-1;
     end
powerprob.Constraints.D1(kk,jj)=...
A(kk,jj) - sum(A(sumidx,jj)/length(sumidx))<=0;

showconstr(powerprob.Constraints.D1) ;
   end
end
powerprob.Constraints.C1=sum(A.*VV)==D; % For limiting the appliances operation for a particular number of hours

showconstr(powerprob.Constraints.C1(:,1:3))

powerprob.Constraints.C3=Pgrid==HLD;

% Options for the optimization algorithm, here we set the max time it can run for
options=optimoptions('intlinprog','Maxtime',10);

% Call the optimization solver to find the best solution

[sol,TotalCost,exitflag,output]=solve(powerprob,options);

sol
struct2table(sol)
A=sol.A

Хотя код генерирует решение,Но решение неверно, поскольку оно нарушает ограничение продолжительности работы приборов.Для первого столбца переменной «A» в решении устройство работает в течение 3 часов, но его продолжительность составляет только 2 часа.

Любая помощь будет принята с благодарностью.Спасибо

...