Optimization giving unexpected results

1 visualización (últimos 30 días)
Gimpy
Gimpy el 11 de Ag. de 2014
Comentada: Gimpy el 11 de Ag. de 2014
Hi, I'm trying to optimize a function which I know the results but matlab is giving me weird results. Here's what I'm trying to do:
max: f(x)= 1815x1+379x2
subject to: -1475x1-112013x2 >= -700000
x1,,x2<=80
x1,x2>=0
How would you do it?
  3 comentarios
Yu Jiang
Yu Jiang el 11 de Ag. de 2014
Editada: Yu Jiang el 11 de Ag. de 2014
Can you share your code? Which built-in function are you using?
Bear in mind that both functions fminsearch and fmincon search for the minimum not the maximum.
-Yu
Gimpy
Gimpy el 11 de Ag. de 2014
f = [1815;379]
A = [-1475 -11203]
b = [-700000]
ub= (ones(1,2)*80)'
lb = zeros(2,1)
x = linprog(f,A,b,[],[],lb,ub)

Iniciar sesión para comentar.

Respuesta aceptada

Yu Jiang
Yu Jiang el 11 de Ag. de 2014
Editada: Yu Jiang el 11 de Ag. de 2014
Hi Gimpy
Since linprog finds the minimum for the function f(x) but the problem you are trying to solve is a maximization problem, I suggest you change f = [1815;379] to f = -[1815;379]. That is to say, you solve the following problem instead
min: f(x)= - 1815x1 -379x2
subject to: -1475x1-112013x2 >= -700000
x1,,x2<=80
x1,x2>=0
There will only be a sign difference on the objective function. The optimal solutions x will be the same.
-Yu
  6 comentarios
Yu Jiang
Yu Jiang el 11 de Ag. de 2014
Hi Gimpy
The answer seems not correct, since -1475*40 -11203*80 = -955240 <=700000, meaning the first inequality constraint is not satisfied.
Also, linprog treats inequality constraints as "<=" by default, not ">=". So your A and b matrices should be A = [1475 11203] and b = 700000
-Yu
Gimpy
Gimpy el 11 de Ag. de 2014
Yes after a double check I realized it's my model that was mispecification not the matlab function that I don't understand. My understanding was correct. thank you

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Get Started with Optimization Toolbox en Help Center y File Exchange.

Etiquetas

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by