2011年3月27日 星期日

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

作者: chungyuandye (養花種魚數月亮賞星星)
標題: Re: [其他] mathematica求整數解
時間: Sun Mar 27 19:29:38 2011

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

5.2的話應該只是把RegionPlot改成InequalityPlot,
不過在6.0之前,這些事額外的套件,必須帶入,所以要先導入InequalityGraphics
我手上沒有5.2,可以的話跟我講一下∼科科∼感恩喔!

<<Graphics`InequalityGraphics`
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[InequalityPlot[cons,{x,-1,3},{y,-,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://cydye1069.blogspot.com
http://tinyurl.com/25nedr2

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.173.128.199

沒有留言:

張貼留言