I am writing a code on mimimising hourly load using mixed integer linear programming. Error returned is that number of cols in Aeq must equal number of elements in f, but actually they are both equal to 168. Kindly explain why this error?

Question

nishtha mehta 2017-5-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/339220-i-am-writing-a-code-on-mimimising-hourly-load-using-mixed-integer-linear-programming-error-returned

回答： gurijala sai harsha 2021-12-14

%total 168 decision variables
%%5*24 as power shiftable variables for power distribution over time
%%%remaining 2*24 for time shiftable appliances to determine the time shift
%%%power pattern will remain same for them
x1=zeros(1,24); x1(1,19)=1;x1(1,20)=1;
x2=zeros(1,24); x2(1,3:5)=1; x2(1,21:22)=1;
x3=ones(1,24);
x4=ones(1,24);
x5=zeros(1,24); x5(1,1:8)=1; x5(1,20:24)=1;
p6=[1 0.5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]';
P6(:,1)= p6(:); 
for k=1:23
P6(:,k+1)=circshift(p6,k);
end 
s6=sum(P6);
p7=[0.8 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]';
   P7(:,1)=p7(:);
   for k=1:23 
P7(:,k+1)=circshift(p7,k);
   end 
s7=sum(P7);
f=[x1 x2 x3 x4 x5 s6 s7];
%bounds 
lb1=ones(1,24); ub1=Inf(1,24);
lb2=ones(1,24); ub2=Inf(1,24);
lb3(:,1:24)=0.12; ub3=Inf(1,24);
lb4=zeros(1,24); ub4(:,1:24)=1.5;
lb5(:,1:24)=0.1; ub5(:,1:24)=3;
lbs6=zeros(1,24); ubs6=ones(1,24);
lbs7=zeros(1,24); ubs7=ones(1,24);
lbi=[lb1 lb2 lb3 lb4 lb5 lbs6 lbs7]; ubi=[ub1 ub2 ub3 ub4 ub5 ubs6 ubs7];
lb=lbi'; ub=ubi';
%constraints
 Ai=zeros(7,168); Ai(1,1:24)=1; Ai(2,25:48)=1; Ai(3,49:72)=1; Ai(4,73:96)=1; Ai(5,97:120)=1; Ai(6,121:144)=1; Ai(7,145:168)=1;
   A=Ai; 
B=[1;3;2.88;3;5;1;1]; 
intcon=(121:168); 
size(A)
   size(f)
[x,fval,eflag] = intlinprog(-f,intcon,A,B,lb,ub)

6 个评论
显示 4更早的评论隐藏 4更早的评论

Walter Roberson 2017-5-9

Okay.. but your problem was that you passed lower bound and upper bound in the position reserved for equality constraints.

nishtha mehta 2017-5-9

Okay Sir, understood

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Jan 2017-5-9

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/339220-i-am-writing-a-code-on-mimimising-hourly-load-using-mixed-integer-linear-programming-error-returned#answer_266302

编辑：Jan 2017-5-9

在 MATLAB Online 中打开

Although you assume that the sizes of the arrays are matching, they do not match. You can use the debugger to find out wh. Type this in the command winodw:

dbstop if error

Now run your code again until it stops in the failing line. Now inspect the values of the used variables to find out, which array has the unexpected size. In addition you can set a breakpoint in the code and step through the program line by line to observe the reasons of the unexpected size.

Please format your code and post a copy of the complete error message. Currently we have to guess, where the error occurres, and this is very hard when the code is not readable.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Walter Roberson 2017-5-9

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/339220-i-am-writing-a-code-on-mimimising-hourly-load-using-mixed-integer-linear-programming-error-returned#answer_266313

在 MATLAB Online 中打开

Aeq = []; beq = [];
[x,fval,eflag] = intlinprog(-f, intcon, A, B, Aeq, beq, lb, ub)

2 个评论
显示无隐藏无

nishtha mehta 2017-5-9

Thank you Sir, it did work. Though it is returning no feasible point that satisfies the constraints and eflag = -2 , but the error that I was encountering previously has disappeared. I will try to reform my code again so as to get a feasible point.

MUTHUSELVI G 2018-7-10