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?
조회 수: 2 (최근 30일)
이전 댓글 표시
%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
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.
채택된 답변
Jan
2017년 5월 9일
편집: Jan
2017년 5월 9일
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
추가 답변 (3개)
Walter Roberson
2017년 5월 9일
Aeq = []; beq = [];
[x,fval,eflag] = intlinprog(-f, intcon, A, B, Aeq, beq, lb, ub)
댓글 수: 2
MUTHUSELVI G
2018년 7월 10일
Have you got the feasible solution..?If it is so please send the code that will be better for understanding.
참고 항목
카테고리
Help Center 및 File Exchange에서 GPU Computing에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!