[ < ] [ > ]   [ << ] [] [ >> ]         [冒頭] [目次] [見出し] [ ? ]

72. simplex


[ < ] [ > ]   [ << ] [] [ >> ]         [冒頭] [目次] [見出し] [ ? ]

72.1 Introduction to simplex

simplexはシンプレックスアルゴリズムを使った線形最適化のパッケージです。

例:

 
(%i1) load("simplex")$
(%i2) minimize_lp(x+y, [3*x+2*y>2, x+4*y>3]);
                  9        7       1
(%o2)            [--, [y = --, x = -]]
                  10       10      5
·

Categories:  Numerical methods ·Optimization ·Share packages ·Package simplex


[ < ] [ > ]   [ << ] [] [ >> ]         [冒頭] [目次] [見出し] [ ? ]

72.2 Functions and Variables for simplex

オプション変数: epsilon_lp

デフォルト値: 10^-8

linear_programの数値計算で使われるイプシロン。

以下も参照してください: linear_program

·

Categories:  Package simplex

関数: linear_program (A, b, c)

linear_programはシンプレックスアルゴリズムの実装です。 linear_program(A, b, c)は、 A.x = bかつ x >= 0を満たすベクトルの中で c.xが可能な最小となるベクトル xを計算します。 引数 Aは行列で、引数 bcはリストです。

linear_programは、最小化ベクトル xと最小値 c.xを含むリストを返します。 もし問題が有界でないなら、 "Problem not bounded!"を返し、 もし問題が実現可能でないなら、 "Problem not feasible!"を返します。

この関数を使うためには、最初に load(simplex);simplexパッケージをロードしてください。

例:

 
(%i2) A: matrix([1,1,-1,0], [2,-3,0,-1], [4,-5,0,0])$
(%i3) b: [1,1,6]$
(%i4) c: [1,-2,0,0]$
(%i5) linear_program(A, b, c);
                   13     19        3
(%o5)            [[--, 4, --, 0], - -]
                   2      2         2

以下も参照してください: minimize_lp, scale_lp, epsilon_lp

関数: maximize_lp (obj, cond, [pos])

いくつかの線形制約 condに従う線形目標関数 objを最大化します。 引数と戻り値の詳細な記述に関しては、 minimize_lpを参照してください。

以下も参照してください: minimize_lp.

関数: minimize_lp (obj, cond, [pos])

いくつかの線形制約 condに従う線形目標関数 objを最小化します。 condは線形等式や不等式のリストです。 厳密な不等式では、 >>=に、 <<=に置き換えられます。 オプションの引数 posは正と仮定される決定変数のリストです。

もし最小が存在するなら、 minimize_lpは目標関数の最小値と最小が得られる決定変数値のリストを含むリストです。 もし問題が有界でないなら、 minimize_lpは "Problem not bounded!"を返し、 もし問題が実現可能でないなら、 "Ploblem not feasible!"を返します。

決定変数はデフォルトでは非負とは仮定されません。 もし決定変数すべてが正なら、 nonegative_lptrueに設定してください。 もし決定変数のいくつかだけが正なら、オプション引数 posの中でそれらをリストしてください。 (これは制約を足すより効率的だということに注意してください。)

minimize_lpは Maximaの linear_program関数で実装されたシンプレックスアルゴリズムを使います。

この関数を使うためには、最初に load(simplex);simplexパッケージをロードしてください。

例:

 
(%i1) minimize_lp(x+y, [3*x+y=0, x+2*y>2]);
                      4       6        2
(%o1)                [-, [y = -, x = - -]]
                      5       5        5
(%i2) minimize_lp(x+y, [3*x+y>0, x+2*y>2]), nonegative_lp=true;
(%o2)                [1, [y = 1, x = 0]]
(%i3) minimize_lp(x+y, [3*x+y=0, x+2*y>2]), nonegative_lp=true;
(%o3)                Problem not feasible!
(%i4) minimize_lp(x+y, [3*x+y>0]);
(%o4)                Problem not bounded!

いかも参照してください: maximize_lp, nonegative_lp, epsilon_lp

オプション変数: nonegative_lp

デフォルト値: false

もし nonegative_lpが trueなら、 minimize_lpmaximize_lpの決定変数すべては正と仮定されます。

以下も参照してください: minimize_lp

·

Categories:  Package simplex


[ << ] [ >> ]           [冒頭] [目次] [見出し] [ ? ]

この文書は市川 雄二によって2013年8月月24日texi2html 1.82を用いて生成されました。