交叉熵平衡SAM表中标红的那句怎么不对了。。大家能不能帮我看看。

0关注
3粉丝

已卖：39份资源

初中生

66%

还不是VIP/贵宾

-

0%

威望: 0 级
论坛币: 66506 个
通用积分: 0.1200
学术水平: 2 点
热心指数: 2 点
信用等级: 2 点
经验: 1295 点
帖子: 8
精华: 0
在线时间: 21 小时
注册时间: 2014-12-7
最后登录: 2016-3-5

楼主

TeamoHe 发表于 2014-12-8 16:40:24 |AI写论文

是否 +2 论坛币

k人参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

*交叉熵平衡SAM表
set ac /sec1,sec2,sec3,sec4,total/;
set i(ac) /sec1*sec4/;
alias(i,j);
table sam(*,*)
      sec1 sec2 sec3    sec4       total
sec1    52    45             150          247
sec2    95    48             90          233
sec3    120    89                            209
sec4                   192                      192
total    267 182    192       240;
parameters
q0(i,j) SAM表各个初始流量,
h0  原始sum of all transaction flows;
q0(i,j)=sam(i,j);
h0=sum((i,j),q0(i,j));
display h0,sam;
variables
q(i,j) 要调整的SAM表中的各个数值，
h    调整SAM表的总值，
hratio 调整和原始两个总数的比例，
z 目标函数的数值，即预期熵值;
positive variable q(i,j);
equations
totalsum 被调整的总数,
hratiodef  hratio的定义和范围,
directentropy 目标函数预期交叉熵，
balance 各个账户平衡的限制条件;

totalsum..h=e=sum((i,j),q(i,j));
hratiodef..hratio=e=h/h0;
directentropy..z=e=sum((i,j),(1/h)*q(i,j)*log(q(i,j)/sam(i,j))-log(hratio));
balance(i)..sum(j,q(i,j))=e=sum(j,q(j,i));
q.l(i,j)=q0(i,j);
h.l=h0;
hratio.lo=0.5;
hratio.up=2;
model sambal /all/;
solve sambal useing nlp minimizing z;
display q.l,h.l,hratio.l;

图片里的****一直提示错误~~~