数学建模社区-数学中国
标题:
使用交叉熵最小化方法求解矩阵，求哪位lingo高手或matlab高手帮忙
[打印本页]
作者:
svanscan
时间:
2013-11-8 16:54
标题:
使用交叉熵最小化方法求解矩阵，求哪位lingo高手或matlab高手帮忙
各位大神，
我想要计算一个
16*16
的矩阵，本来想用lingo计算，在网上找了一个类似的程序但是运行出来，结过出错，希望给位走过路过的朋友，给予指点，用lingo 或者是matlab都可以，谢谢了，有些急。有关模型描述包括目标函数和约束条件都在附件里面。lingo的运算程序如下：
) q) v0 n L8 [. b3 A3 t6 n
model:!16家银行；sets：
) t) O% n. \/ ]; q/ x
A/1..16/:ai；
8 z, V3 O) o2 [- W5 Y& T/ W
B/1..16/:bj;
- S" e8 I7 i/ V7 r1 P8 }4 k. u
links(A,B):X,X0;
5 i# I3 x1 b: K0 N& y* }7 P
!X1=X_*;
) M9 E+ H7 ]: L6 p
endsets
3 T9 p1 E9 d+ q' q
data:
! [, C0 m7 K, e$ p9 G
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;
0 Q' Y3 ]' V3 ?" L* V
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;
* o; b- o' F# e" H8 `% u
enddata
4 d4 R% l7 X+ h! p& e! N( w/ E5 ~
!目标函数：
0 n, r( u) ]' S. e7 e) Z7 ^6 M
[obj]min=@sum(links(i,j)/i#ne#j:X(i,j)*(@log(X(i,j))-@log(X0(i,j))));
q1 c ?$ F& q: R* L0 r' C0 y
aa=@sum(A(i):ai(i));
6 O: Q4 c2 f# }. ?; r
bb=@sum(b(j):bj(j));
9 H- [/ x6 n o" M, x; V: C: j
!约束条件；
9 U9 T0 K+ H' @8 o
@for(A(i):@sum(B(j):X(i,j))=X(i,i)+ai(i););
* ~$ m! L# W: f% b+ T! ?3 ]
@for(B(j):@sum(A(i):X(i,j))=X(j,j)+bj(j););
- e! o) {( F" P9 P
@for(links(i,j):X0(i,j)=@if(i#EQ#j,0,ai(i)*bj(j)));
0 i# Q! C( y4 ]% [
end
, E0 z& O* G E. Y5 K
. p* r7 s. i% s
1 H4 P% H0 T" Q5 ]7 b$ q
7 X: M1 U/ Z) c7 c( M7 e# q5 x
矩阵模型.doc
2013-11-8 16:50 上传
点击文件名下载附件
下载积分: 体力 -2 点
25 KB, 下载次数: 11, 下载积分: 体力 -2 点
所用的模型
作者:
svanscan
时间:
2013-11-8 16:58
lingo运行程序中的约束条件应该为：
' c7 K; [# T" v
!约束条件；
! e N' ^$ I/ f' J
@for(A(i):@sum(B(j):X(i,j))=X(i,i)+ai(i););
5 a- [- r' u P
@for(B(j):@sum(A(i):X(i,j))=X(j,j)+bj(j););
5 q2 ]$ D# W* E# C7 C
@for(links(i,j):X0(i,j)=@if(i#EQ#j,0,ai(i)*bj(j)));
' @! Y5 b* E+ i& b3 N
end
- ~7 L7 u0 d$ p) Z
, R3 g7 \# e2 p0 R) `: x3 r
那个愤怒的表情是：
@
作者:
yankeke
时间:
2013-11-9 17:37
" H) C' J0 R( n) H- M
什么事 啊不懂。。。。。。。。。。。。。。。
作者:
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希望尽快能联系到前辈
+ M2 k! ^5 d1 k5 J" f
作者:
lxfdjt158
时间:
2015-1-25 15:37
求指导。本人也想学习
( d) \! b" l- R# \& v2 q) p( F
作者:
lxfdjt158
时间:
2015-1-25 15:40
我现在也在做这个研究，能否指导一下
4 k. a! v* }8 \: ~" R
我还是没能用lingo
! e2 Z' |8 }0 ?$ Q5 u/ \: x- _. U
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
