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

svanscan

各位大神，我想要计算一个16*16的矩阵，本来想用lingo计算，在网上找了一个类似的程序但是运行出来，结过出错，希望给位走过路过的朋友，给予指点，用lingo 或者是matlab都可以，谢谢了，有些急。有关模型描述包括目标函数和约束条件都在附件里面。lingo的运算程序如下：
model:!16家银行；sets：
A/1..16/:ai；
4 Q6 U$ Z$ ^' D( \. N+ O# ]1 t: aB/1..16/:bj;
links(A,B):X,X0;
! M3 j- r, ], t% V$ }0 q!X1=X_*;
endsets
% p. }' c6 w2 D1 m/ Pdata:
. u) E" ~4 V( a* O6 F+ P  L+ A" gai=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;
9 v5 V& x+ }( pbj=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
2 d: y8 Z: X/ }1 A( Y" K0 s& X!目标函数：
' N# Y/ `* D3 D! f. \+ P6 G# X[obj]min=@sum(links(i,j)/i#ne#j:X(i,j)*(@log(X(i,j))-@log(X0(i,j))));
' _* {0 l9 C5 m; S" R( Laa=@sum(A(i):ai(i));
bb=@sum(b(j):bj(j));
!约束条件；
@for(A(i):@sum(B(j):X(i,j))=X(i,i)+ai(i););
. F- s' n; ~4 e- ?8 z3 X@for(B(j):@sum(A(i):X(i,j))=X(j,j)+bj(j););
- U3 H) `9 V, f+ K% ?' s@for(links(i,j):X0(i,j)=@if(i#EQ#j,0,ai(i)*bj(j)));
, h# |  U# m! ?" N( }; ~end

8 @/ B! S+ s" c& Z. o, ^+ |% V9 I$ B

svanscan

lingo运行程序中的约束条件应该为：
4 e! q& O# q6 E9 E/ ^3 o!约束条件；
@for(A(i):@sum(B(j):X(i,j))=X(i,i)+ai(i););
1 b# Y  A: C' p0 b@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
9 h- z2 U3 M. @1 e! H. `
& b' n  d, I% a* b# F  j那个愤怒的表情是： @

yankeke

$ k5 ]9 u8 p! y4 M什么事 啊不懂。。。。。。。。。。。。。。。

hellowen

~\(≧▽≦)/~

dingxuemeiwang

看看。。。

PL6160082

谢谢楼主  不明觉厉

lxfdjt158

希望给个联系方式，本人论文有涉及到，可是不懂，都快抓狂了。希望前辈指导一下，联系方式2385676071.电话18382425941希望尽快能联系到前辈
8 e8 [1 N: }% \7 d4 G8 o

lxfdjt158

求指导。本人也想学习

lxfdjt158

我现在也在做这个研究，能否指导一下
我还是没能用lingo
8 w# v* S! @7 W5 m" L: C! R