Optimization giving unexpected results

2014 8 11
1 个回答
回答已采纳
更新时间：2014 8 11
4 次查看（30 天）

0 个投票

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 个评论
显示 1更早的评论 隐藏 1更早的评论

Gimpy
Gimpy 2014-8-11
maybe I would like to do it with simplex method(My benchmark in Excel) or wahatever method is the best... thank you!
Yu Jiang
Yu Jiang 2014-8-11
编辑：Yu Jiang 2014-8-11
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 2014-8-11
f = [1815;379]
A = [-1475 -11203]
b = [-700000]
ub= (ones(1,2)*80)'
lb = zeros(2,1)
x = linprog(f,A,b,[],[],lb,ub)

请先登录，再进行评论。

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

 采纳的回答

Yu Jiang
Yu Jiang 2014-8-11
编辑：Yu Jiang 2014-8-11

1 个投票

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 个评论
显示 4更早的评论 隐藏 4更早的评论

Gimpy
Gimpy 2014-8-11
tried that option, still not working
Yu Jiang
Yu Jiang 2014-8-11
What is the expected answer?
Yu Jiang
Yu Jiang 2014-8-11
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 2014-8-11
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

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Solver Outputs and Iterative Display 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by