使用交叉熵最小化方法求解矩阵，求哪位lingo高手或matlab高手帮忙

svanscan

各位大神，我想要计算一个16*16的矩阵，本来想用lingo计算，在网上找了一个类似的程序但是运行出来，结过出错，希望给位走过路过的朋友，给予指点，用lingo 或者是matlab都可以，谢谢了，有些急。有关模型描述包括目标函数和约束条件都在附件里面。lingo的运算程序如下：
7 _4 E) C; a8 `4 `, U7 {( j$ n( M& ^model:!16家银行；sets：
1 L4 }( p1 p' ~3 P7 g. O' @A/1..16/:ai；
7 h/ k( H7 X  y0 g( V' f. M7 uB/1..16/:bj;
3 @# U6 Q5 X7 c- u- h$ ]! Ylinks(A,B):X,X0;
% y% }7 F/ T: C! I$ M!X1=X_*;
endsets
data:
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;
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;
enddata
0 b  L8 i0 ^4 J!目标函数：
5 t% [" Y' \/ S- C8 L% \[obj]min=@sum(links(i,j)/i#ne#j:X(i,j)*(@log(X(i,j))-@log(X0(i,j))));
aa=@sum(A(i):ai(i));
bb=@sum(b(j):bj(j));
!约束条件；
  K" ?* Q0 P9 E8 v2 n) q$ u+ y! Y@for(A(i):@sum(B(j):X(i,j))=X(i,i)+ai(i););
- `+ s3 x" y; ^  i6 o& D1 ]$ P@for(B(j):@sum(A(i):X(i,j))=X(j,j)+bj(j););
@for(links(i,j):X0(i,j)=@if(i#EQ#j,0,ai(i)*bj(j)));
end
! D! ]0 C5 {" C( g7 I1 F
, p( N. M3 c/ y8 s' H

4 }2 y- c" ]5 l# F

lxfdjt158

我现在也在做这个研究，能否指导一下
) F7 F/ {) p2 [  }& ~; q我还是没能用lingo

lxfdjt158

求指导。本人也想学习

lxfdjt158

希望给个联系方式，本人论文有涉及到，可是不懂，都快抓狂了。希望前辈指导一下，联系方式2385676071.电话18382425941希望尽快能联系到前辈
' s! [0 w& G4 H  u3 M

PL6160082

谢谢楼主  不明觉厉

dingxuemeiwang

看看。。。

hellowen

~\(≧▽≦)/~

yankeke

什么事 啊不懂。。。。。。。。。。。。。。。

svanscan

lingo运行程序中的约束条件应该为：
!约束条件；
  w6 q9 F: j& U2 {@for(A(i):@sum(B(j):X(i,j))=X(i,i)+ai(i););
* W$ V( {- j$ d0 L@for(B(j):@sum(A(i):X(i,j))=X(j,j)+bj(j););
+ k3 c% s" ~5 D9 B* N9 T/ Q@for(links(i,j):X0(i,j)=@if(i#EQ#j,0,ai(i)*bj(j)));
5 x% B4 C- t* C* V% B2 Xend
: |' J9 n) k3 U$ `
那个愤怒的表情是： @