For each fmincon it overwrites the previous values, meaning at the end there is only 1 set of values when there should be 100. The code is running 100 times and displays all sets of values but not storing all sets. How do I resolve this?

Question

Scott Sanders 2018-4-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/392100-for-each-fmincon-it-overwrites-the-previous-values-meaning-at-the-end-there-is-only-1-set-of-values

评论： Scott Sanders 2018-4-3

 for X=1:1:10
    for Y=1:1:10
    ObjFcn = @myObjective;
    x0 = [10 0.001 7];
    LB = [0 0 0];
    UB = [50 0.5 7.5];
    ConsFcn = @(x)myConstraints(x,X,Y);
    [x] = fmincon(ObjFcn,x0,[],[],[],[],LB,UB,ConsFcn);
    disp(x)
    disp(myObjective(x))
    end
 end

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

请先登录，再进行评论。

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

Answer 1

Alan Weiss 2018-4-2

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/392100-for-each-fmincon-it-overwrites-the-previous-values-meaning-at-the-end-there-is-only-1-set-of-values#answer_313097

在 MATLAB Online 中打开

Your code is inefficient in that you define many things inside the loop that should be assigned outside the loop. Using Walter's suggestion, try this:

ObjFcn = @myObjective;
x0 = [10 0.001 7];
LB = [0 0 0];
UB = [50 0.5 7.5];
x = zeros(10,10); % initialization
for X=1:10
  for Y=1:10
    ConsFcn = @(x)myConstraints(x,X,Y);
    x(X,Y) = fmincon(ObjFcn,x0,[],[],[],[],LB,UB,ConsFcn);
    disp(x(X,Y))
    disp(myObjective(x(X,Y)))
  end
end

5 个评论
显示 3更早的评论隐藏 3更早的评论

Walter Roberson 2018-4-2

在 MATLAB Online 中打开

ObjFcn = @myObjective;
x0 = [10 0.001 7];
LB = [0 0 0];
UB = [50 0.5 7.5];
x = zeros(10,10,3 ); % initialization
fvals = inf(10,10); % initialization
for X=1:10
  for Y=1:10
    ConsFcn = @(x)myConstraints(x,X,Y);
    [x(X,Y,:), fvals(X,Y)] = fmincon(ObjFcn,x0,[],[],[],[],LB,UB,ConsFcn);
    disp(x(X,Y))
    disp(myObjective(x(X,Y)))
  end
end

Scott Sanders 2018-4-3

在 MATLAB Online 中打开

That seems to have done the trick. The final code that works is:

 ObjFcn = @myObjective;
 x0 = [10 0.001 7];
 LB = [0 0 0];
 UB = [100 0.3 7.5];
 x=zeros(10,10,3);
 fvals=inf(10,10);
 for X=1:10
     for Y=1:10
     ConsFcn = @(x)myConstraints(x,X,Y);
     [x(X,Y,:),fvals(X,Y)] = fmincon(ObjFcn,x0,[],[],[],[],LB,UB,ConsFcn);
     disp(x(X,Y,:))
     disp(myObjective(x(X,Y,:)))
     end
 end