Fmincon problems with A and b matrices that differ depending on for loop

Question

Andrew Alkiviades 2011-6-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/8716-fmincon-problems-with-a-and-b-matrices-that-differ-depending-on-for-loop

Hi

Im doing a nonlinear optimization using fmincon that includes nonlinear equality constraints linear inequality constraints (These depend on values from two excel spreadsheets.

This is my code for the A and b matrices for j=1:3 for i=1:5

        L = xlsread('hourly electric demand.xlsx');
        Lth = xlsread('hourly thermal demand.xlsx');
        if xin(1)>= L(j,i) &&  xin(2)>=Lth(j,i);
             A(1,:)=[FIT RHI 0 0 -(FIT-Xe);-1 0 -1 0 0; 0 -1 0 -1 0];
             b = [max(Pe*L(j,i))+Pth*max(Lth(j,i))+ Operation; -L(j,i);...
                 -Lth(j,i)];
        elseif xin(1)<=L(j,i) &&  xin(2)<=Lth(j,i);
           A(2,:)=[FIT RHI -Pe -Pth 0;-1 0 -1 0 0; 0 -1 0 -1 0];
          b = [max(Pe*L(j,i)+Pth*Lth(j,i))+ Operation;-L(j,i); -Lth(j,i)];
        elseif xin(1)>=L(j,i) &&  xin(2)<=Lth(j,i);
            A(3,:)=[FIT RHI 0 -Pth -(FIT-Xe);-1 0 -1 0 0; 0 -1 0 -1 0];
            b = [max(Pe*L(j,i)+Pth*Lth(j,i))+ Operation; -L(j,i); -Lth(j,i)];
        elseif xin(1)<= L(j,i) &&  xin(2)>= Lth(j,i);
            A(4,:)=[FIT RHI -Pe -Pth 0;-1 0 -1 0 0; 0 -1 0 -1 0];
      b = [max(Pe*L(j,i)+Pth*Lth(j,i))+ Operation; -L(j,i); -Lth(j,i)];
        end
    end
end

and this is the code for the objective function

f = sum(((0.537*(L(j,i)-xin(1)))+(0.185*(Lth(j,i)-xin(2)))));

However I am getting the following error

??? Subscripted assignment dimension mismatch.

Error in ==> invoke at 42 A(1,:)=[FIT RHI 0 0 -(FIT-Xe);-1 0 -1 0 0; 0 -1 0 -1 0];

Do I need to make changes to the objective function since it only includes xin(1) and xin(2) whereas I am constraining to 5 difference xin?

thanks

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

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2011-6-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/8716-fmincon-problems-with-a-and-b-matrices-that-differ-depending-on-for-loop#answer_12016

[FIT RHI 0 0 -(FIT-Xe);-1 0 -1 0 0; 0 -1 0 -1 0] is clearly 2 dimensional, but you are trying to assign it to a vector A(1,:) . That isn't going to work.