数学建模社区-数学中国

标题: 没有最优解，是不是程序错了， [打印本页]

作者: svanscan 时间: 2014-3-27 10:58
标题: 没有最优解，是不是程序错了，
model:!16家银行;
sets:
A/1..16/:ai;
B/1..16/:bj;
links(A,B):X,X0;
!X1=X_*;
endsets
data:
ai=0.005858522,0.005087844,0.062226485,0.027219144,0.081445428,0.014868129,0.003513967,0.026944443,
0.030721218,0.03729378,0.044584185,0.102144315,0.0431425,0.05324312,0.399638873,0.062142449;
bj=0.019169918,0.004728038,0.08866658,0.024510633,0.033035937,0.02180133,0.010255439,0.06842088,
0.019047496,0.06073349,0.068611704,0.281815485,0.043940778,0.055856946,0.176789794,0.022615551;
enddata
!目标函数;
[obj]min=@sum(links(i,j)|i#ne#j:X(i,j)*(@log(X(i,j))-@log(X0(i,j))));
aa=@sum(A(i):ai(i));
bb=@sum(b(j):bj(j));
!约束条件;
@for(A(i):@sum(B(j):X(i,j))=X(i,i)+ai(i););
@for(B(j):@sum(A(i):X(i,j))=X(j,j)+bj(j););
@for(links(i,j):X0(i,j)=@if(i#EQ#j,0,ai(i)*bj(j)));
end

不知道是哪里出错了，请走过路过的高手给予指点

作者: madio 时间: 2014-3-30 11:29
我这里lingo11开了全局求解器算了一下，没有可行解，可能是你的约束条件有矛盾的地方。

作者: wujianjack2 时间: 2014-3-30 12:55
如果可以的话请附上模型的公式图片，方便对照。如站长所言，你的程序我用LINGO 14 Extended，使用默认设置也报告无解，poorly scaled.

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