1

Mathprog メソッドを使用して問題を解決しようとしていますが、R を使用して txt ファイルを読み取ると問題が発生します。

どこに問題があるかはわかっていると思いますが、それを解決する方法がわかりません。

これがコード全体ですが、 param a および param tc 3 変数で使用すると、問題は最終的に発生します。

set I;
/* plants*/

set J;
/* customers*/

set T;
/* types*/

set S;
/* modes*/

param pc{i in I, s in S};
/* variable production cost */

param tc{i in I, j in J, t in T};
/* transportation cost */

param a{i in I, t in T, s in S};
/* quantity of product produced */

param b{j in J, t in T};
/* transportation cost */

var z{i in I, s in S} >= 0;
/* intensity of the plant */

var u{i in I, j in J, t in T} >= 0;
/* amount of product shipped from plant to costumer*/

minimize cost: sum{i in I, s in S} pc[i,s] * z[i,s] + sum{t in T, j in J, i in I} tc[i,j,t] * u[i,j,t];

/* production transportation problem */

s.t. first{i in I, j in J, t in T}: sum{i in I} u[i,j,t] = b[j,t];
/* First Constraint */

s.t. second{i in I, j in J, t in T, s in S}: sum{j in J} u[i,j,t] <= sum{s in S} a[i,t,s] * z[i,s];
/* Second Constraint */

s.t. third{i in I, s in S}: sum{s in S} z[i,s] <= 1;
/* Third Constraint */

s.t. fourth{i in I, s in S}: z[i,s] >= 0;
/* Fourth Constraint */

s.t. fifth{i in I, j in J, t in T}: u[i,j,t] >= 0;
/* Fifth Constraint */

data;

set I := P1 P2;

set J := C1 C2 C3 C4 C5;

set T := T1 T2 T3 T4;

set S := S1;

param pc : S1 :=
P1     5000
P2     7000;

param tc : P1T1 P1T2 P1T3 P1T4 P2T1 P2T2 P2T3 P2T4:= 
C1    10 15 10 20 10 10 12 12
C2    20 10 15 10 12 12 10 15
C3    20 15 10 10 15 10 12 18
C4    10 10 10 10 10 10 10 10
C5    5  10 15 20 20 23 21 15;

param a : S1P1 S1P2 := 
T1    50 100
T2    150 50
T3    200 80
T4    150 200;

param b : T1 T2 T3 T4 := 
C1    80 100 80 120
C2    50 100 105 100
C3    80 105 100 100
C4    80 50 15 100
C5    50  100 150 200;

end;

これらの関数と一緒に使用して、R で読み取ります。

model.lp <- Rglpk_read_file ("mathprog_file.txt", type = "MathProg", verbose = F)
4

1 に答える 1

0

おそらく、Rglpk は正常に動作します。私が見ているのは、間違った数学プログラムの構文です。

制約では、同じインデックスを 2 回宣言しています。

s.t. first{i in I, j in J, t in T}: sum{i in I} u[i,j,t] = b[j,t];

i左の制約ステートメントで既に定義されていますが、sum ステートメントでも再度使用しiています。合計のように別の表記法でこれを修正して{ii in I}から、適切なiまたはiiを参照しながら使用できますu。

三次元パラメータの定義にも誤りがあります。ここで例を見てください: 例のデータブロックの説明とその方法の説明。

于 2016-12-06T07:12:01.553 に答える