QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 21495|回复: 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的运算程序如下:
    . A( S- G; |& ]. S4 Kmodel:!16家银行;sets:
    1 A% z* n  w6 i9 v) @  NA/1..16/:ai;
    6 N) w: t! I, v- Z2 k- G' PB/1..16/:bj;3 J9 D" D% K  @: A, ]
    links(A,B):X,X0;
    0 W! Q" D+ r: Q8 k! Y!X1=X_*;
    ) h5 L9 {& U: x/ ^: J# Zendsets4 M2 ~* v! E4 c- k/ D3 I
    data:
    ' j5 x8 `+ r& uai=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;
    ! _$ `) [# s$ }+ cbj=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;' A3 j# A/ d' q1 t$ C
    enddata% O5 L0 |) F- `+ q) N& z/ e
    !目标函数:# c* t: S/ U2 {7 W; @& J2 x; R8 A/ C
    [obj]min=@sum(links(i,j)/i#ne#j:X(i,j)*(@log(X(i,j))-@log(X0(i,j))));
    ' d! ]' z6 G5 B( {7 _aa=@sum(A(i):ai(i));
    % F# h4 c- O. u: ?bb=@sum(b(j):bj(j));
    2 x3 ]6 b' t3 d2 l9 e+ k* S!约束条件;
    # Q, Y  P( U2 X1 U: s1 ]@for(A(i):@sum(B(j):X(i,j))=X(i,i)+ai(i););
    , \# {+ m3 S  P& e9 Q@for(B(j):@sum(A(i):X(i,j))=X(j,j)+bj(j););
    1 f) I9 V% p* i/ ~: F! G1 u@for(links(i,j):X0(i,j)=@if(i#EQ#j,0,ai(i)*bj(j)));  j* U5 [7 x0 ^& \2 @# ~/ j
    end
    4 q( h9 u6 _! F3 y7 {9 w  I3 W( e; C
    + s6 T. o7 |; D9 e% l  J3 u$ b& @, {, w, U. q" ^0 j

    ; _; ]5 }3 m' }5 N$ @; }$ h

    矩阵模型.doc

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

    所用的模型

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

    6

    主题

    10

    听众

    42

    积分

    升级  38.95%

  • TA的每日心情
    难过
    2014-3-27 10:56
  • 签到天数: 8 天

    [LV.3]偶尔看看II

    自我介绍
    master
    lingo运行程序中的约束条件应该为:% p- y+ e. t2 A
    !约束条件;* C* p: @# P5 ~" L  x8 t% ]
    @for(A(i):@sum(B(j):X(i,j))=X(i,i)+ai(i););
    2 p- w' J' i8 X$ X6 `4 P@for(B(j):@sum(A(i):X(i,j))=X(j,j)+bj(j););: I5 a! O8 P9 v  t
    @for(links(i,j):X0(i,j)=@if(i#EQ#j,0,ai(i)*bj(j)));
    9 @$ ~! M8 x6 z+ U1 @end, `5 b0 ~2 q4 v% M  q$ x
    ) T5 j7 D% h( f; E
    那个愤怒的表情是: @
    回复

    使用道具 举报

    yankeke        

    0

    主题

    1

    听众

    2

    积分

    升级  40%

    该用户从未签到

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

    使用道具 举报

    hellowen        

    6

    主题

    9

    听众

    69

    积分

    升级  67.37%

  • TA的每日心情
    开心
    2014-12-26 15:12
  • 签到天数: 23 天

    [LV.4]偶尔看看III

    自我介绍
    华工工科大学生

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

    回复

    使用道具 举报

    0

    主题

    4

    听众

    38

    积分

    升级  34.74%

  • TA的每日心情
    奋斗
    2014-1-24 16:04
  • 签到天数: 12 天

    [LV.3]偶尔看看II

    回复

    使用道具 举报

    PL6160082        

    3

    主题

    7

    听众

    127

    积分

    升级  13.5%

  • TA的每日心情
    慵懒
    2014-4-8 19:36
  • 签到天数: 26 天

    [LV.4]偶尔看看III

    自我介绍
    大学二年级...
    回复

    使用道具 举报

    lxfdjt158        

    0

    主题

    2

    听众

    4

    积分

    升级  80%

    该用户从未签到

    自我介绍
    在大学参加过数学建模比赛,获得国家一等奖,省二等奖,研究生也参加了数学建模大赛,获得优秀参赛奖,现任学术部部长。
    希望给个联系方式,本人论文有涉及到,可是不懂,都快抓狂了。希望前辈指导一下,联系方式2385676071.电话18382425941希望尽快能联系到前辈# Z; V4 F+ E; u
    回复

    使用道具 举报

    lxfdjt158        

    0

    主题

    2

    听众

    4

    积分

    升级  80%

    该用户从未签到

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

    使用道具 举报

    lxfdjt158        

    0

    主题

    2

    听众

    4

    积分

    升级  80%

    该用户从未签到

    自我介绍
    在大学参加过数学建模比赛,获得国家一等奖,省二等奖,研究生也参加了数学建模大赛,获得优秀参赛奖,现任学术部部长。
    我现在也在做这个研究,能否指导一下
    7 Y( o- G, f$ I我还是没能用lingo7 d, a% {* d9 d' {8 g5 r
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-4-21 00:42 , Processed in 0.776472 second(s), 104 queries .

    回顶部