查看: 22026|回复: 8

[问题求助] 使用交叉熵最小化方法求解矩阵，求哪位lingo高手或matlab高手帮忙

字体大小: 正常放大

svanscan

6 主题	10 听众	42 积分

升级 38.95%

TA的每日心情

	难过 2014-3-27 10:56

签到天数: 8 天

[LV.3]偶尔看看II

自我介绍: master

电梯直达

1^#

发表于 2013-11-8 16:54 |只看该作者 |正序浏览

|招呼Ta 关注Ta

各位大神，我想要计算一个16*16的矩阵，本来想用lingo计算，在网上找了一个类似的程序但是运行出来，结过出错，希望给位走过路过的朋友，给予指点，用lingo 或者是matlab都可以，谢谢了，有些急。有关模型描述包括目标函数和约束条件都在附件里面。lingo的运算程序如下：
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

矩阵模型.doc

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

所用的模型

zan

matlab

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

相关帖子

使用道具举报

lxfdjt158

0 主题	2 听众	4 积分

升级 80%

该用户从未签到

自我介绍: 在大学参加过数学建模比赛，获得国家一等奖，省二等奖，研究生也参加了数学建模大赛，获得优秀参赛奖，现任学术部部长。

9^#

发表于 2015-1-25 15:40 |只看该作者 |招呼Ta 关注Ta

我现在也在做这个研究，能否指导一下
我还是没能用lingo

我也说一句

发表

使用道具举报

lxfdjt158

0 主题	2 听众	4 积分

升级 80%

该用户从未签到

自我介绍: 在大学参加过数学建模比赛，获得国家一等奖，省二等奖，研究生也参加了数学建模大赛，获得优秀参赛奖，现任学术部部长。

8^#

发表于 2015-1-25 15:37 |只看该作者 |招呼Ta 关注Ta

求指导。本人也想学习

我也说一句

发表

使用道具举报

lxfdjt158

0 主题	2 听众	4 积分

升级 80%

该用户从未签到

自我介绍: 在大学参加过数学建模比赛，获得国家一等奖，省二等奖，研究生也参加了数学建模大赛，获得优秀参赛奖，现任学术部部长。

7^#

发表于 2015-1-25 15:24 |只看该作者 |招呼Ta 关注Ta

希望给个联系方式，本人论文有涉及到，可是不懂，都快抓狂了。希望前辈指导一下，联系方式2385676071.电话18382425941希望尽快能联系到前辈

我也说一句

发表

使用道具举报

PL6160082

3 主题	7 听众	127 积分

升级 13.5%

TA的每日心情

	慵懒 2014-4-8 19:36

签到天数: 26 天

[LV.4]偶尔看看III

自我介绍: 大学二年级...

6^#

发表于 2013-12-9 15:48 |只看该作者 |招呼Ta 关注Ta

谢谢楼主不明觉厉

我也说一句

发表

使用道具举报

dingxuemeiwang

0 主题	4 听众	38 积分

升级 34.74%

TA的每日心情

	奋斗 2014-1-24 16:04

签到天数: 12 天

[LV.3]偶尔看看II

5^#

发表于 2013-11-15 20:27 |只看该作者 |招呼Ta 关注Ta

看看。。。

我也说一句

发表

使用道具举报

hellowen

6 主题	9 听众	69 积分

升级 67.37%

TA的每日心情

	开心 2014-12-26 15:12

签到天数: 23 天

[LV.4]偶尔看看III

自我介绍: 华工工科大学生

群组: 2013电工杯A题讨论群组

4^#

发表于 2013-11-9 18:59 |只看该作者 |招呼Ta 关注Ta

~\(≧▽≦)/~

我也说一句

发表

使用道具举报

yankeke

0 主题	1 听众	2 积分

升级 40%

该用户从未签到

自我介绍: 制作沙盘模型：15110202628

3^#

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

什么事啊不懂。。。。。。。。。。。。。。。

我也说一句

发表

使用道具举报

svanscan

6 主题	10 听众	42 积分

升级 38.95%

TA的每日心情

	难过 2014-3-27 10:56

签到天数: 8 天

[LV.3]偶尔看看II

自我介绍: master

2^#

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

lingo运行程序中的约束条件应该为：
!约束条件；
@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

那个愤怒的表情是： @

我也说一句

发表

使用道具举报

返回列表

帐号		密码		只需要一步，快速开始		注册地址	找回密码

[问题求助] 使用交叉熵最小化方法求解矩阵，求哪位lingo高手或matlab高手帮忙

相关帖子

浏览过的版块

QQ

电话咨询

关于我们| 联系我们| 诚征英才| 对外合作| 产品服务|