QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 21888|回复: 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的运算程序如下:# u6 z( h( F% ]( W# c* ^6 p6 n1 b# r
    model:!16家银行;sets:
    ) P! H4 e$ K9 e  x! N6 H* RA/1..16/:ai;8 l' E& j3 X. i/ w
    B/1..16/:bj;' @8 p5 P/ b* P+ G
    links(A,B):X,X0;
    7 ~* m( b8 p" }9 d3 G2 V) D!X1=X_*;+ L" t2 L& I& d* m% X
    endsets
    9 D7 r) T. N/ ]6 f& W4 Sdata:
    1 z1 A/ q0 R9 Jai=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;
    ) P  Y/ m! Q3 zbj=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;9 K/ r( _* w4 m7 X, g! P. Q
    enddata% _" Y# @$ v7 u6 t2 X
    !目标函数:
      N/ o& C) L4 m0 Z+ l7 l[obj]min=@sum(links(i,j)/i#ne#j:X(i,j)*(@log(X(i,j))-@log(X0(i,j))));' v* j% @+ \( d- e
    aa=@sum(A(i):ai(i));
    1 p1 _& R, g5 u& h/ \bb=@sum(b(j):bj(j));
    - x8 e/ E2 `6 A1 Y: [6 _7 q/ h!约束条件;9 N% m! A( N( S
    @for(A(i):@sum(B(j):X(i,j))=X(i,i)+ai(i););- z0 b& g! m1 |. A; _! l8 \
    @for(B(j):@sum(A(i):X(i,j))=X(j,j)+bj(j););
    1 U: y& X2 a. z; ~! s  k! j@for(links(i,j):X0(i,j)=@if(i#EQ#j,0,ai(i)*bj(j)));0 w  E# N4 c) E3 ]4 J4 f
    end
    7 j& o3 b4 ^- D7 y9 n
    3 d3 [1 h9 ?1 |7 d- z* M7 ~! J1 u
    5 b4 w2 H$ H# F! B* w
    8 z8 ^/ J0 G0 `* ?

    矩阵模型.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运行程序中的约束条件应该为:
    1 p, ?3 t4 v7 A! j!约束条件;
    : F7 |/ F0 y3 b) F@for(A(i):@sum(B(j):X(i,j))=X(i,i)+ai(i););
    0 h+ B6 w) K/ x3 |@for(B(j):@sum(A(i):X(i,j))=X(j,j)+bj(j););
    & M4 h" O& _9 o6 x@for(links(i,j):X0(i,j)=@if(i#EQ#j,0,ai(i)*bj(j)));6 P3 Z( f7 u3 @( s, ?# Q5 c, q
    end
    4 M; @  _2 B/ ^, ^" Y. d
    ( b1 h( u- I, Y# Z2 f那个愤怒的表情是: @
    回复

    使用道具 举报

    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希望尽快能联系到前辈
      S& T* I% K. s& U+ x& f7 F
    回复

    使用道具 举报

    lxfdjt158        

    0

    主题

    2

    听众

    4

    积分

    升级  80%

    该用户从未签到

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

    使用道具 举报

    lxfdjt158        

    0

    主题

    2

    听众

    4

    积分

    升级  80%

    该用户从未签到

    自我介绍
    在大学参加过数学建模比赛,获得国家一等奖,省二等奖,研究生也参加了数学建模大赛,获得优秀参赛奖,现任学术部部长。
    我现在也在做这个研究,能否指导一下
    . g  S- U) d0 G+ t9 ^$ _  r我还是没能用lingo8 e% a4 I( l/ X# ~# l
    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-7-6 10:32 , Processed in 0.528285 second(s), 103 queries .

    回顶部