2011年3月27日 星期日

[自動轉寄] Re: [其他] mathematica求整數解

作者: chungyuandye (養花種魚數月亮賞星星) 看板: Math
標題: Re: [其他] mathematica求整數解
時間: Sun Mar 27 17:50:57 2011

※ 引述《iamcrazyforu (愛宅)》之銘言:
: 用mathematica解線性規畫
: 是否能夠限制解必須為整數呢?
: 另外如果要輸出圖的話
: 有沒有辦法把最佳整數解標示在圖上呢?

方法很多,LinearProgramming, Minimize, FindMinimum, NMinimize
都可以用來求解線性規劃問題,當然包含整數規劃。

min x+2y
st 3x+4y>=5
x>=0,y>=0

obj=x+2y;
cons=And[3x+4y>=5,x>=0,y>=0];
ansR=Minimize[{obj,cons},{x,y}]
ansI=Minimize[{obj,cons,Element[{x,y},Integers]},{x,y}]
Show[RegionPlot[cons,{x,0,3},{y,0,3}],
ContourPlot[obj==ansR[[1]],{x,-1,3},{y,-1,3},ContourStyle->{Red,Thickness[0.01]
ContourPlot[obj==ansI[[1]],{x,-1,3},{y,-1,3},ContourStyle->{Blue,Thickness[0.01
Epilog->{Green,PointSize[0.05],Point[{x,y}/.ansR[[2]]],Point[{x,y}/.ansI[[2]]]}
PlotRange->{{-1,3},{-1,3}}
]

--
我打研究室走過 那獨坐電腦前的容顏如苦瓜的糾結
靈感不來 長壽的煙霧不散
研究室如小小的寂寞的城 恰如商管的電梯向晚

http://chungyuandye.twbbs.org

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.173.240.190
※ 編輯: chungyuandye 來自: 218.173.240.190 (03/27 17:55)

沒有留言:

張貼留言