12 / 2 页下一页

查看: 17907|回复: 12

求解答，我用lingo计算最小交叉熵的解，老是提示错误，怎么解决？

[复制链接]

字体大小: 正常放大

svanscan

6 主题	10 听众	42 积分

升级 38.95%

TA的每日心情

	难过 2014-3-27 10:56

签到天数: 8 天

[LV.3]偶尔看看II

自我介绍: master

电梯直达

1^#

发表于 2013-11-6 17:12 |只看该作者 |正序浏览

|招呼Ta 关注Ta

本帖最后由 svanscan 于 2013-11-6 17:12 编辑

model:!16家银行；
sets：
A/1..16/:ai；
B/1..16/:bj;
links(A,B):X,X0;
!X1=X_*;
endsets
data:
ai=0.049105298,0.19706518,0.079929381,0.040404067,0.00307861,0.029535259,0.030970463,0.081245849,
0.053001907,0.122173474,0.000229636,0.067807575,0.359324615,0.051351782,0.004097295,0.008038271;
bj=0.06029401,0.017726814,0.034382538,0.041643595,0.006514545,0.075761924,0.011835925,0.017374484,
0.11483671,0.17727039,0.029125245,0.127436517,0.242528843,0.025613802,0.002559618,0.015095039;
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

file:///C:\Users\lenovo\AppData\Roaming\Tencent\Users\605193618\QQ\WinTemp\RichOle\U)R3W8`~KYA$65%A~QPA6G7.jpg[img=0,1]file:///C:\Users\lenovo\AppData\Roaming\Tencent\Users\605193618\QQ\WinTemp\RichOle\L0OIYB[O6@$12DM)ZFW}1MR.jpg[/img][img=0,1]file:///C:\Users\lenovo\AppData\Roaming\Tencent\Users\605193618\QQ\WinTemp\RichOle\L0OIYB[O6@$12DM)ZFW}1MR.jpg[/img]file:///C:\Users\lenovo\AppData\Roaming\Tencent\Users\605193618\QQ\WinTemp\RichOle\U)R3W8`~KYA$65%A~QPA6G7.jpg
提示错误如下
error code 1017
error text
unexpected jacobian overflow. Contact Lingo Systems for patch availability(1).

zan

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

	13^# 无效楼层，该帖已经被删除

	12^# 无效楼层，该帖已经被删除

svanscan

6 主题	10 听众	42 积分

升级 38.95%

TA的每日心情

	难过 2014-3-27 10:56

签到天数: 8 天

[LV.3]偶尔看看II

自我介绍: master

11^#

发表于 2013-11-8 17:01 |只看该作者 |招呼Ta 关注Ta

wujianjack2 发表于 2013-11-8 12:52
匆匆写了几句，仍无解，可能有哪个地方写错了，你自己看看吧。
MODEL:
SETS:

朋友，附件里是我计算的模型，有些急，还请你费心看一下，指点迷津，到底是哪点出来问题

矩阵模型.doc

25 KB, 下载次数: 7, 下载积分: 体力 -2 点

我要计算的模型

我也说一句

发表

使用道具举报

svanscan

6 主题	10 听众	42 积分

升级 38.95%

TA的每日心情

	难过 2014-3-27 10:56

签到天数: 8 天

[LV.3]偶尔看看II

自我介绍: master

10^#

发表于 2013-11-8 13:12 |只看该作者 |招呼Ta 关注Ta

wujianjack2 发表于 2013-11-8 12:52
匆匆写了几句，仍无解，可能有哪个地方写错了，你自己看看吧。
MODEL:
SETS:

万分感谢呐！！！

我也说一句

发表

使用道具举报

wujianjack2

21 主题	97 听众	3110 积分

TA的每日心情

	奋斗 2014-3-2 00:26

签到天数: 243 天

[LV.8]以坛为家I

9^#

发表于 2013-11-8 12:52 |只看该作者 |招呼Ta 关注Ta

匆匆写了几句，仍无解，可能有哪个地方写错了，你自己看看吧。
MODEL:
SETS:
ASET:AI; !SIXTEEN BANKS;
BSET:BJ; !SIXTEEN BANKS;
LINKS(ASET,BSET):XZERO,XSTAR,RATIO_LOG_STAR_ZERO;
ENDSETS
DATA:
AI=0.049105298,0.19706518,0.079929381,0.040404067,0.00307861,0.029535259,0.030970463,0.081245849,
0.053001907,0.122173474,0.000229636,0.067807575,0.359324615,0.051351782,0.004097295,0.008038271;
BJ=0.06029401,0.017726814,0.034382538,0.041643595,0.006514545,0.075761924,0.011835925,0.017374484,
0.11483671,0.17727039,0.029125245,0.127436517,0.242528843,0.025613802,0.002559618,0.015095039;
ENDDATA
[OBJ] MIN=@SUM(LINKS(I,J):XSTAR(I,J)*RATIO_LOG_STAR_ZERO(I,J));
@FOR(LINKS(I,J):RATIO_LOG_STAR_ZERO(I,J)=@IF(XSTAR(I,J) #NE# 0 #AND# XZERO(I,J) #NE# 0,@LOG(XSTAR(I,J))-@LOG(XZERO(I,J)),0));
@FOR(BSET(J):@SUM(ASET(I):XSTAR(I,J))=BJ(J));
@FOR(ASET(I):@SUM(BSET(J):XSTAR(I,J))=AI(I));
@FOR(LINKS(I,J)|I #NE# J:XZERO(I,J)=AI(I)*BJ(J));
@FOR(LINKS(I,J)|I #EQ# J:XZERO(I,J)=0);
@SUM(LINKS(I,J):XZERO(I,J)+XSTAR(I,J))=1;
END