fmincon should return a scalar value error

50 次查看（过去 30 天）

显示更早的评论

Hassan garni 2018-1-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/379051-fmincon-should-return-a-scalar-value-error

评论： Torsten 2022-8-30

采纳的回答： Nicolas Schmit

在 MATLAB Online 中打开

Hello,

I am trying to do an optimization for a large function using a "fmincon" function. I have got an error and couldn't find the problem. The main function has so many trigonometric functions (cos, sin. ...). The following is the error obtained :

Error using fmincon (line 609)
  Supplied objective function must return a scalar value.
    Error in Energy_test2 (line 74)
    x = fmincon(objective,x0,A,b,Aeq,beq,lb,ub,'flcon');

The following is the main script using "fmincon". I also have checked the matrices dimensions.

objective = @(x)( cos(x(1)).*sin(x(2))+ ....(cos(t).*cos((Ay)-(x(1)))).*...*Ar*e));
x0 = [2,42];
disp(['initial objective:' num2str(objective(x0))])
A = [];
b=[];
Aeq = [];
beq = [];
x0 = [1,4];
lb = [0,0]; 
ub = [90,90] ;
nonlinecon = @flcon;
x = fmincon(objective,x0,A,b,Aeq,beq,lb,ub,'flcon');

I appreciate if somebody help with this issue and how can be solved.

Thank you.

Hassan

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

AJEET Kumar 2022-8-30

Hello sir, Sir how to decide the initial condition in fmincon problem.because fir different initial condition we get different solution.so please ask me any formula for take initial condition.

Torsten 2022-8-30

Take the initial condition in the neighbourhood of the solution you want to get from fmincon :-)

请先登录，再进行评论。

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

采纳的回答

Nicolas Schmit 2018-1-26

3
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/379051-fmincon-should-return-a-scalar-value-error#answer_301819

在 MATLAB Online 中打开

fmincon() crashed because it expected the objective function to return a scalar value but got something that was not a scalar.

Test your objective function on the initial conditions and check that the output is a scalar.

objective(x0) % Scalar?

If you need further help, post the complete program.

13 个评论
显示 11更早的评论隐藏 11更早的评论

Sanyam Maheshwari 2020-7-16

在 MATLAB Online 中打开

I am also getting the same error. Please help.

Below is the Data File and Below that is the objective function file. Please help.

global SL V p r k m n C Z s 
SL = 0.75;
V = 94100;
p = [0.07,0.18,0.2,0.3];
r = [55,47,45,49];
k = [33,28,29,30];
m = 4;
n= 4;
C = [78,69,70,73;
    64,68,56,59;
    34,39,42,41;
    52,47,48,45];
Z =[250	250	250	250;
    320	320	320	320;
    440	440	440	440;
    350	350	350	350];
s = [110,95,99,100;
    110,95,99,100;
    110,95,99,100;
    110,95,99,100];
bguess = randi(4,4);
bguess = (bguess<2);
Q0 = [zeros(4,4) bguess];
% call the solver
xopt = fmincon(@EP, Q0,[],[],[],[],[],[]);

function G =  EP(Q,b)
global s r k C m n 
rng;
y = rand(4);
b = randi(4,4);
b = (b<2);
q = 0;
for i = 1:m
    for j = 1:n
        q = q + y(i,j)*Q(i,j);
    end
end
x = sym('x',[1,4]); %x as symbolic variable
f1 = 0;
for i = 1:m
    for j = i:n 
        f1  = f1 + s(i,j)*b(i,j)*x(j);
    end
end       
f2 = 0;
for i = 1:m
    for j = 1:n
        f2 = f2+C(i,j)*b(i,j)*y(i,j)*Q(i,j);
    end
end
f3 = 0;
        f3a = 0;
        for i = 1:m
            for j = 1:n
                f3a = f3a + y(i,j)*Q(i,j)-x(j);
            end
        end
for i  = 1:m
    for j = 1:n
        f3 = f3 + b(i,j)*r(j)*f3a;
    end
end
f4 = 0;
for i = 1:m
    for j = 1:n
        f4 = f4+ s(j)*b(i,j)*y(i,j)*Q(i,j);
    end
end
f5 = 0;
        f5a = 0;
        for i = 1:m
            for j = 1:n
                f5a = f5a + x(j) -y(i,j)*Q(i,j);
            end
        end
    
for i  = 1:m
    for j = 1:n
        f5 = f5 + b(i,j)*k(j)*f5a;
    end
end
F1 = int(((f1 - f2 + f3).*normpdf(x, 400, 100)),0,q);
F2 = int(((f4 - f2 - f5).*normpdf(x, 400, 100)),q,Inf);
G = -(F1 +F2);
end

Walter Roberson 2020-7-16

在 MATLAB Online 中打开

x = sym('x',[1,4]); %x as symbolic variable

x is length 4

F1 = int(((f1 - f2 + f3).*normpdf(x, 400, 100)),0,q);

x is length 4, so normpdf(x, 400, 100) is a vector of length 4. You integrate an expression involving the vector of length 4, from 0 to q, where q is 0, and integral of anything from 0 to 0 is 0, so you get back a vector of 4 zeros.

F2 = int(((f4 - f2 - f5).*normpdf(x, 400, 100)),q,Inf);

Likewise, the normpdf is a vector of length 4. You integrate an expression involving the vector of length 4, from q to inf, where q is 0. You do not specify a variable of integration, so it does symvar() and picks the first one, which will be x1. So you integrate the expression of length 4 with respect to x1. As the expression involves x1, x2, x3, x4, the result is going to be of length 4 involving x2, x3, x4 (the x1 vanishes because it is a definite integration with respect to x1.) Because the result involves unresolved variables x2, x3, x4, there is no way to convert the result into a numeric scalar such as needed by fmincon -- for example -sum(F1+F2) will not work because although you would have reduced the length 4 to a scalar, it would still involve the unresolved variables x2, x3, x4.

Now, I could suggest that perhaps you wanted something like

syms X
F2 = int(((f4 - f2 - f5).*normpdf(X, 400, 100)),X,q,Inf);

which would give you a scalar result. But the f4-f2-f5 expression involves the unresolved variables x1, x2, x3, x4, so again the result would not be something you could convert to numeric.

Note: if you do manage to get down to a scalar without any unresolved variables, then beware that fmincon will not accept symbolic numbers: you must convert the symbolic expression into double()

Sanyam Maheshwari 2020-7-17

Hi Walter

I just posted in order to get my answer as fast as possible.

As suggested by you I tried using syms x.

That again resulted in different errors.

If possible if you could suggest me to use any other solver. I am new to MATLAB So I don't have much experience.

Here x represents demand of 4 different products.

I am solving an optimization problem in which is somewhere linked with newspaper vendor problem of operations management.

Also can you please explain how I have to use double() with an example.

It will be really helpfull

Walter Roberson 2020-7-17

在 MATLAB Online 中打开

Here x represents demand of 4 different products.

The demands of 4 different products can be optimized in more than 16 different ways. If we call them P1, P2, P3, P4, then one way would be to find the mimimum of P1, ignoring P2, P3, P4. Another would be to find the mnimum of P2, ignoring P1, P3, P4. Another would be to find the product of P2 and P3 ignoring P1 and P4. Other way would be to find the sum of P1+P4, ignoring P2 and P3. And so on. You might even ask about min() of the P1, P2, P3, P4 values --- but if you use min() then the function is not differentiable and fmincon() cannot be used.

Each of the products is an "objective", and fmincon() can only optimize one objective. When you have four objectives, you need to either use a different routine such as gamultiobj() or else you need to find a function that takes the four individual objectives and forms a unified score.

For example, suppose you could lower P1, P2, P3 by 0.02 each, but only at the cost of increasing P4 by 0.01: should you do that? Suppose the mininum for P1 is at its lower bound, but that forces P3 to be relatively large: should you do that?

请先登录，再进行评论。