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?

 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

Answer 1

[Alan Weiss]

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

[Scott Sanders]

I made the adjustments you suggested however it is still having trouble running the code. an error message is coming up with "Assignment has more non-singleton rhs dimensions than non-singleton subscripts"

[Walter Roberson]

Instead of x = zeros(10,10) make it zeros(10,10,3) and instead of assigning to x(X, Y) assign to x(X, Y, :)

[Scott Sanders]

That has almost completely worked. It is now giving me 100 values of each of the variables being optimized but is still overwriting the actually objective function value.

[Walter Roberson]

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]

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

Thank you all for your help!

Answer 2

[Walter Roberson]

Instead of assigning to x assign to x(X, Y)

[Scott Sanders]

which x in that code? I have tried this and couldn't get it to work.