数学建模社区-数学中国

标题: 帮忙看看这段代码！老是报雅克比错误。 [打印本页]

作者: wahehehehehe 时间: 2010-11-20 08:31
标题: 帮忙看看这段代码！老是报雅克比错误。
model:
sets:
wh /1..256/:ver;
h  /1..16/;
circulation(h,wh):flow;
endsets

data:
ver=
  0  4  5  3  2  5  3  6  3  2  6  3  2  4  1  2
  3  0  3  5  2  2  4  5  2  5  2  3  4  6  4  0
  3  2  0  4  8  2  5  2  1  3  3  2  4  2  2  2
  0  3  5  0  5  2  6  3  4  4  1  2  5  4  2  4
  5  4  0  7  0  2  5  2  5  7  3  5  1  4  3  4
  5  3  7  4  2  0  3  2  2  4  3  4  2  4  2  1
  3  7  2  4  3  2  0  2  3  5  1  2  5  4  3  1
  5  2  2  3  3  2  1  0  3  5  3  5  1  3  6  2
  2  6  2  3  1  4  2  3  0  3  4  5  1  2  3  5
  7  5  2  2  3  1  2  0  5  0  3  6  1  5  4  6
  5  3  6  2  1  3  3  4  2  1  0  1  5  3  3  5
  2  2  4  1  4  4  3  3  4  7  6  0  2  7  3  4
  6  3  1  4  3  2  5  3  4  6  7  4  0  4  2  3
  2  2  5  5  5  3  1  2  1  2  5  5  4  0  3  1
  2  2  7  2  3  3  3  2  4  4  3  3  5  4  0  4
  8  0  4  2  4  2  3  7  6  3  3  2  3  2  2  0
;
enddata

[obj] max=@sum(circulation(i,j)|
(@mod(j-1,16) #ne# @floor((j-1)/16) #and# @mod(j-1,16) #ne# (i-1)):flow(i,j));

@for(circulation(i,j)|
(@mod(j-1,16) #eq# @floor((j-1)/16) #or# @mod(j-1,16) #eq# (i-1)):flow(i,j)=0);

@for(wh(j)| @mod(j-1,16) #ne# @floor((j-1)/16):
@sum(h(i)| @mod(j-1,16) #ne# (i-1) : flow(i,j)) <= ver(j));

@for(wh(j)| @mod(j-1,16) #ne# @floor((j-1)/16):
@sum(h(i)| @mod(j-1,16) #ne# (i-1) : flow(i,j)) =
   @sum(circulation(m,n)| (m-1) #eq# @mod(j-1,16) #and# @mod(n-1,16) #eq# @floor((j-1)/16):flow(m,n)) );

end

作者: wahehehehehe 时间: 2010-11-20 08:34
本帖最后由 wahehehehehe 于 2010-11-20 10:46 编辑

坛子里下载的11. 12版都不好用。
找了一个9.0版的不报雅克比错误。但是当问题规模增大，变量数达到24*(24-1)(24-1)=12696时，又开始报错：not enough virtual system memory.

郁闷！

作者: Greedy 时间: 2010-11-20 13:28
帮你顶一下~

作者: 柳亦非台 时间: 2011-4-27 23:55
唉，没人解释啊~

作者: shengshengchina 时间: 2011-9-21 08:49
你的软件问题，我的可以求解出：
  Global optimal solution found.
  Objective value:                            750.0000
  Infeasibilities:                            0.000000
  Total solver iterations:                         1108

                     Variable          Value       Reduced Cost
                     VER( 1)       0.000000          0.000000
                     VER( 2)       4.000000          0.000000
                     VER( 3)       5.000000          0.000000
                     VER( 4)       3.000000          0.000000
                     VER( 5)       2.000000          0.000000
                     VER( 6)       5.000000          0.000000
                     VER( 7)       3.000000          0.000000
                     VER( 8)       6.000000          0.000000
                     VER( 9)       3.000000          0.000000
                     VER( 10)       2.000000          0.000000

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)