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

svanscan

各位大神，我想要计算一个16*16的矩阵，本来想用lingo计算，在网上找了一个类似的程序但是运行出来，结过出错，希望给位走过路过的朋友，给予指点，用lingo 或者是matlab都可以，谢谢了，有些急。有关模型描述包括目标函数和约束条件都在附件里面。lingo的运算程序如下：
) n# N. R2 S  `2 m3 Bmodel:!16家银行；sets：
) |  [. O9 Z$ S8 g/ g# bA/1..16/:ai；
  C9 i+ k9 B6 Q0 z, mB/1..16/:bj;
8 k2 X- B# q$ M/ Tlinks(A,B):X,X0;
0 C  `% |; j5 M6 R. n. |) r  N. c2 J!X1=X_*;
endsets
data:
: E: u/ |1 h( n$ m- U: [# Yai=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;
( i1 p9 E; S' {8 R& H4 g* jbj=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
!目标函数：
6 }+ `/ H8 o( W9 d6 K5 B[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));
" B. q, K. q! G! R) h/ g% ]bb=@sum(b(j):bj(j));
  o4 b# a- X; m4 j# s!约束条件；
" H  c! w6 f/ s1 o@for(A(i):@sum(B(j):X(i,j))=X(i,i)+ai(i););
+ m( T# v& a0 A: E3 j@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

7 w8 _* c+ S8 p
) m5 x$ p1 s, D- H; c8 @2 {

svanscan

lingo运行程序中的约束条件应该为：
!约束条件；
@for(A(i):@sum(B(j):X(i,j))=X(i,i)+ai(i););
2 b' o) {1 ~' L/ l@for(B(j):@sum(A(i):X(i,j))=X(j,j)+bj(j););
( H; S) ?$ ?. t2 U! w" E# \@for(links(i,j):X0(i,j)=@if(i#EQ#j,0,ai(i)*bj(j)));
+ t* X) h5 U" o) C, \end
9 i) v5 ]7 {: I' i% p2 W8 n
: D9 o; H0 x3 O# e: }* c+ I那个愤怒的表情是： @

yankeke

0 r4 N1 a4 B% n# _2 |. Q什么事 啊不懂。。。。。。。。。。。。。。。

hellowen

~\(≧▽≦)/~

dingxuemeiwang

看看。。。

PL6160082

谢谢楼主  不明觉厉

lxfdjt158

希望给个联系方式，本人论文有涉及到，可是不懂，都快抓狂了。希望前辈指导一下，联系方式2385676071.电话18382425941希望尽快能联系到前辈

lxfdjt158

求指导。本人也想学习
  Y, T0 m0 I( _& ~4 Z8 ~( g& \/ Q

lxfdjt158

我现在也在做这个研究，能否指导一下
% x1 J! x/ R. E/ [1 m. U' y' H& D% k我还是没能用lingo
, d# ^- Y5 W5 W* A* _