数学建模社区-数学中国
标题:
使用交叉熵最小化方法求解矩阵，求哪位lingo高手或matlab高手帮忙
[打印本页]
作者:
svanscan
时间:
2013-11-8 16:54
标题:
使用交叉熵最小化方法求解矩阵，求哪位lingo高手或matlab高手帮忙
各位大神，
我想要计算一个
16*16
的矩阵，本来想用lingo计算，在网上找了一个类似的程序但是运行出来，结过出错，希望给位走过路过的朋友，给予指点，用lingo 或者是matlab都可以，谢谢了，有些急。有关模型描述包括目标函数和约束条件都在附件里面。lingo的运算程序如下：
- X3 v. p2 f. a* \
model:!16家银行；sets：
0 |5 e! J1 |! s- h& c6 d
A/1..16/:ai；
; ~" H- q* W; ]( @; A) I. u
B/1..16/:bj;
M- m# _( k/ X
links(A,B):X,X0;
& h* w. e7 \! I0 b% }0 H
!X1=X_*;
" J9 T {" i2 ?8 P# o2 _" J
endsets
+ O; V" U# d$ _; h- A# T
data:
' S. N, h6 w; D
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;
( r0 ?+ s) l, I) @
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;
) l, ], a0 z$ @ F( ~( s; T
enddata
3 E* W" n& x2 T
!目标函数：
# w, ~) V. c# L& y, S
[obj]min=@sum(links(i,j)/i#ne#j:X(i,j)*(@log(X(i,j))-@log(X0(i,j))));
1 E' ?1 r4 m% a( I+ @- `+ x
aa=@sum(A(i):ai(i));
7 o& y0 H6 a: q/ U0 A9 x8 x9 a0 G
bb=@sum(b(j):bj(j));
( F7 v$ z8 ?/ v( D
!约束条件；
+ } ]: \" B0 d( r& a% G
@for(A(i):@sum(B(j):X(i,j))=X(i,i)+ai(i););
: r: v2 V9 M$ f4 [' z
@for(B(j):@sum(A(i):X(i,j))=X(j,j)+bj(j););
& C- u8 x9 W5 N# U
@for(links(i,j):X0(i,j)=@if(i#EQ#j,0,ai(i)*bj(j)));
( ^. n4 `$ N6 k" ?7 A
end
9 E% N. I5 x d' x: J, E* V% P% u
; v6 l/ w2 U# y* v( I. Z
: M( B A3 g3 ?; W: ~1 H0 x
6 }3 g# F' y7 B& o
矩阵模型.doc
2013-11-8 16:50 上传
点击文件名下载附件
下载积分: 体力 -2 点
25 KB, 下载次数: 11, 下载积分: 体力 -2 点
所用的模型
作者:
svanscan
时间:
2013-11-8 16:58
lingo运行程序中的约束条件应该为：
# k: @8 V$ w: Y4 a% h
!约束条件；
' X) o, t5 G3 d( B6 x) Q0 h: F
@for(A(i):@sum(B(j):X(i,j))=X(i,i)+ai(i););
( V# n9 B9 r. p8 z, P0 A
@for(B(j):@sum(A(i):X(i,j))=X(j,j)+bj(j););
. Y% H# V' T9 }; E" r/ `& ~$ j
@for(links(i,j):X0(i,j)=@if(i#EQ#j,0,ai(i)*bj(j)));
0 U5 d$ V1 A5 W3 ~
end
; u4 ]! B9 \6 B) g5 Y4 Q/ L
! O% U& j3 T2 @, L. u
那个愤怒的表情是：
@
作者:
yankeke
时间:
2013-11-9 17:37
& w+ K D& [ Y, u5 d
什么事 啊不懂。。。。。。。。。。。。。。。
作者:
hellowen
时间:
2013-11-9 18:59
~\(≧▽≦)/~
作者:
dingxuemeiwang
时间:
2013-11-15 20:27
看看。。。
作者:
PL6160082
时间:
2013-12-9 15:48
谢谢楼主 不明觉厉
作者:
lxfdjt158
时间:
2015-1-25 15:24
希望给个联系方式，本人论文有涉及到，可是不懂，都快抓狂了。希望前辈指导一下，联系方式2385676071.电话18382425941希望尽快能联系到前辈
8 N( X! l' H& o# l$ w# v" k
作者:
lxfdjt158
时间:
2015-1-25 15:37
求指导。本人也想学习
. o6 ?- X7 M: }, Z/ l0 |6 m( w
作者:
lxfdjt158
时间:
2015-1-25 15:40
我现在也在做这个研究，能否指导一下
9 D8 q" t3 o& R4 N9 r
我还是没能用lingo
3 c; \( o1 j0 e+ I/ H1 h/ F
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
