QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 21496|回复: 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的运算程序如下:
    # a5 O; T7 W/ Y% Z( m4 p5 P) emodel:!16家银行;sets:
    9 t# F9 s! B. }$ M/ }A/1..16/:ai;
    ) Q+ O3 i/ M% D+ |+ C6 ]$ j, ]( ~6 gB/1..16/:bj;8 I4 ]/ N) k4 B* t3 Z0 x
    links(A,B):X,X0;- w+ P) l# s4 d5 K1 i
    !X1=X_*;6 l( ~0 x: o% }3 o6 W2 r8 h% w# y
    endsets
    : a) @& ~6 c$ |- H2 @data:( o. C0 s' W" f' ]9 H3 e- [
    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;
      d9 d1 b( i. c; N% Nbj=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;
    : X" G/ a# Q/ z! ?" j& l. denddata1 n4 v6 Y# m" R- Y) d9 k/ j- H
    !目标函数:
    3 L5 q* G$ H" o6 P[obj]min=@sum(links(i,j)/i#ne#j:X(i,j)*(@log(X(i,j))-@log(X0(i,j))));
    ! A9 q1 ^$ N+ m3 [1 l- L0 D  baa=@sum(A(i):ai(i));
    % ^! `" b, g4 L  @# Zbb=@sum(b(j):bj(j));
    # t/ u& I7 e, X5 O!约束条件;1 G4 C" c! r- o8 l$ u
    @for(A(i):@sum(B(j):X(i,j))=X(i,i)+ai(i););
    & [3 }. h( X! v2 [% w@for(B(j):@sum(A(i):X(i,j))=X(j,j)+bj(j););
    0 s2 Q/ E  M/ x5 @' v6 G: V@for(links(i,j):X0(i,j)=@if(i#EQ#j,0,ai(i)*bj(j)));9 \1 ^3 B( R$ Q- o' E6 ]
    end, l( r! R# k$ t& o/ h- g

      o" g+ y) C% |! N% T* U8 b6 |7 e# l. e

    & y! e+ M/ V% W0 O' }

    矩阵模型.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运行程序中的约束条件应该为:- h* m/ Q3 _8 Z& Q  H
    !约束条件;4 h4 _: U/ C: e  e0 I
    @for(A(i):@sum(B(j):X(i,j))=X(i,i)+ai(i););
    * `% D5 ]5 y5 {@for(B(j):@sum(A(i):X(i,j))=X(j,j)+bj(j););. Q$ \9 w5 z1 ]$ m, n
    @for(links(i,j):X0(i,j)=@if(i#EQ#j,0,ai(i)*bj(j)));
    ' E8 `7 `2 H# P* ^* _end9 C/ T' F" r, y) l( d0 E; S
    8 _- S  {4 @: k) k* J/ t+ U: y
    那个愤怒的表情是: @
    回复

    使用道具 举报

    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希望尽快能联系到前辈
    + B8 F; T% S* ?( S9 K, b
    回复

    使用道具 举报

    lxfdjt158        

    0

    主题

    2

    听众

    4

    积分

    升级  80%

    该用户从未签到

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

    使用道具 举报

    lxfdjt158        

    0

    主题

    2

    听众

    4

    积分

    升级  80%

    该用户从未签到

    自我介绍
    在大学参加过数学建模比赛,获得国家一等奖,省二等奖,研究生也参加了数学建模大赛,获得优秀参赛奖,现任学术部部长。
    我现在也在做这个研究,能否指导一下
    ! B9 @0 b' I+ P7 K( l( e$ v- G我还是没能用lingo) X. t1 T, V& |* p4 h
    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-4-21 04:48 , Processed in 0.707590 second(s), 104 queries .

    回顶部