求助

kenstorm

MODEL:
Titel Location problem;
sets:
      demand/1..6/:a,b,d;
      supply/1..2/:x,y,e;
      Link(demand,supply):c;
endsets
data:
!Locations for the demand(需求点的位置);
a=1.25,8.75,0.5,5.75,3,7.25;
b=1.25,0.75,4.75,5,6.5,7.75;
!quantites of the demand and supply(供需量);
d=3,5,4,7,6,11;e=20,20;
enddata
init:
!initial Locations for the supply(初始点);
x,y=5 1 2 7;
endinit
!Objective function(目标);
[OBJ]min=@sum(Link(i,j):c(i,j)*((x(j)-a(i))^2+(y(j)-b(i))^2)^(1/2));
!demand constraints(需求约束);
@for(demand(i):[DEMAND_CON]@sum(supply(j):c(i,j))=d(i););
!supply constraints(供应约束);
@for(supply(i):[SUPPLY_CON]@sum(demand(j):c(i,j))<=e(i););
@for(supplyfree(X);@free(Y););
END





model:
sets:
    demand/1..6/:a,b,d;
    supply/1..2/:x,y,e;
    link(demand,supply):c;
endsets
data:
a=1.25,8.75,0.5,5.75,3,7.25;
b=1.25,0.75,4.75,5,6.5,7.75;
d=3,5,4,7,6,11;e=20,20;
enddata
init:
x,y=5,1,2,7;
endinit
[OBJ]min=@sum(link(i,j):c(i,j)*((x(j)-a(i))^2+(y(j)-b(i))^2)^(1/2));
@for(demand(i):[demand_con]@sum(supply(j):c(i,j))=d(i););
@for(supply(i):[supply_con]@sum(supply(j):c(i,j))<=e(i););
@for(supplyfree(x);@free(y););
end
以上两段代码，前者无法运算，但是我看不出有什么错误的地方，请高手指点，还有，后面一个如果加上Titel Location problem;也不能求解，请问加标题正确的代码是怎样的，谢谢啦
不明白为什么要强迫我选好书共享，我是求助帖啊，汗。。。。。。。

hubulwm

第一段错误之处：
1.英语单词拼错：Titel Location problem;!语法错误，Titel应改为Title;
第二段你加上了错误的语句当然不能运行.加标题正确的代码是
MODEL:

TITLE 标题;

2.@for(supply(i):[SUPPLY_CON]@sum(demand(j):c(i,j))<=e(i););!集合属性错误：应该为c(j,i);
3.@for(supplyfree(X);@free(Y););!语法错误：supplyfree(X)应改为supply:@free(X);
完整的程序如下：
MODEL:
Title Location problem;
sets:
      demand/1..6/:a,b,d;
      supply/1..2/:x,y,e;
      Link(demand,supply):c;
endsets
data:
!Locations for the demand(需求点的位置);
a=1.25,8.75,0.5,5.75,3,7.25;
b=1.25,0.75,4.75,5,6.5,7.75;
!quantites of the demand and supply(供需量);
d=3,5,4,7,6,11;
e=20,20;
enddata
init:
!initial Locations for the supply(初始点);
x,y=5 1 2 7;
endinit
!Objective function(目标);
[OBJ]min=@sum(Link(i,j):c(i,j)*((x(j)-a(i))^2+(y(j)-b(i))^2)^(1/2));
!demand constraints(需求约束);
@for(demand(i):[DEMAND_CON]@sum(supply(j):c(i,j))=d(i););
!supply constraints(供应约束);
@for(supply(i):[SUPPLY_CON]@sum(demand(j):c(j,i))<=e(i););
@for(supply:@free(X);@free(Y););
END

运行结果：
  Local optimal solution found.
  Objective value:                              85.26604
  Total solver iterations:                            18


  Model Title: Location problem

                       Variable           Value        Reduced Cost
                          A( 1)        1.250000            0.000000
                          A( 2)        8.750000            0.000000
                          A( 3)       0.5000000            0.000000
                          A( 4)        5.750000            0.000000
                          A( 5)        3.000000            0.000000
                          A( 6)        7.250000            0.000000
                          B( 1)        1.250000            0.000000
                          B( 2)       0.7500000            0.000000
                          B( 3)        4.750000            0.000000
                          B( 4)        5.000000            0.000000
                          B( 5)        6.500000            0.000000
                          B( 6)        7.750000            0.000000
                          D( 1)        3.000000            0.000000
                          D( 2)        5.000000            0.000000
                          D( 3)        4.000000            0.000000
                          D( 4)        7.000000            0.000000
                          D( 5)        6.000000            0.000000
                          D( 6)        11.00000            0.000000
                          X( 1)        3.254883            0.000000
                          X( 2)        7.250000           0.8084079E-07
                          Y( 1)        5.652332            0.000000
                          Y( 2)        7.750000           0.2675276E-06
                          E( 1)        20.00000            0.000000
                          E( 2)        20.00000            0.000000
                       C( 1, 1)        3.000000            0.000000
                       C( 1, 2)        0.000000            4.008540
                       C( 2, 1)        0.000000           0.2051358
                       C( 2, 2)        5.000000            0.000000
                       C( 3, 1)        4.000000            0.000000
                       C( 3, 2)        0.000000            4.487750
                       C( 4, 1)        7.000000            0.000000
                       C( 4, 2)        0.000000           0.5535090
                       C( 5, 1)        6.000000            0.000000
                       C( 5, 2)        0.000000            3.544853
                       C( 6, 1)        0.000000            4.512336
                       C( 6, 2)        11.00000            0.000000

                            Row    Slack or Surplus      Dual Price
                            OBJ        85.26604           -1.000000
                 DEMAND_CON( 1)        0.000000           -4.837363
                 DEMAND_CON( 2)        0.000000           -7.158911
                 DEMAND_CON( 3)        0.000000           -2.898893
                 DEMAND_CON( 4)        0.000000           -2.578982
                 DEMAND_CON( 5)        0.000000          -0.8851584
                 DEMAND_CON( 6)        0.000000            0.000000
                 SUPPLY_CON( 1)        0.000000            0.000000
                 SUPPLY_CON( 2)        4.000000            0.000000


若有其他问题发到
hubulwm@qq.com

treefly001

支持一下。。。。。