请问书中第5.6节直接交叉熵法，用张老师提供的代码运行和书中贴的表不一样

如果用GAMS Studio编译文件，打开Edit---Encoding---Reload with----system。

没有完整的GAMS代码，很难查出原因。如果是教科书里面的代码，请告诉我是哪一个？

请将NLP的求解器设定为PATHNLP，或在Solve之前插入代码，"Option nlp=pathnlp;"。

谢谢老师

请问老师，社会核算矩阵里各部门的进口关税是如何计算出的？统计年鉴只有总的进口关税。

应该是资本总量QKS增加10%吧

是不是参数校准出错了？参数校准后能复制出原来平衡的模型吗？

大家好，我也正在学习张欣老师的这本CGE模型，网络上的资源很好，我希望能够建一个CGE学习的交流群，如果大家想加入的话，可以加我qq：1078130537

大家好，我现在也在学习张欣老师的这本CGE模型的书，做练习的时候遇到问题，网上资源也比较少，希望能建立一个学习交流群（扣扣：一零七八一三零五三七 ）

这个应该是根据某些已有数据，自己按比例拆分吧

各位同学好，我也在学习张欣老师的CGE模型这本书，希望能够加您一起组建一个CGE的交流群

guan_xinjia@163.com

1078130537@qq.com

请问你是Mac版吗？我也出现了乱码，请问您最后如何解决的呢？？

老师好，同学们好！我用直接交叉熵第14章例题14-7-2的Sam表进行调平，却出现的如下的错误

**** Exec Error at line 83: log: FUNC SINGULAR: x = 0


代码是：

*编制GAMS程序语言，用直接交叉熵平衡SAM表。
*第十四章例题7.2数据
set ac   /sec1,sec2soe,sec2pri,sec2fdi,sec3,com1,com2,com3,lab,cap,hhr,hhu,ent,
        gov,vatl,vatk,bustax,tariff,invsav,row,total/;
set a(ac)  /sec1,sec2soe,sec2pri,sec2fdi,sec3/;
set c(ac)  /com1,com2,com3/;
set f(ac)  /lab,cap/;
set hoh(ac)  /hhr,hhu/;
set vat(ac) /vatl,vatk/;
*set acnt(ac)  除了总数之外的所有账户;
set i(ac) /sec1,sec2soe,sec2pri,sec2fdi,sec3,com1,com2,com3,lab,cap,hhr,hhu,ent
        ,gov,vatl,vatk,bustax,tariff,invsav,row/;

*acnt(ac)=YES; acnt('total')=NO;
*i(ac)=YES; i('total')=NO;

alias (ac,acp),(a,ap),(c,cp),(f,fp),(hoh,hp),(vat,vatp);
alias (i,j);

table sam(ac,acp)
$INCLUDE C:\Users\phili\Desktop\GAMs\必读文件夹\张欣 第二版所有CGE模型的GAMS程序 解压密码 gzcbs2017\chap14QXSAM.inc

*display sam

*这里SAM表是从外面用txt文件引进


parameters
    Q0(i,j)    initial value SAM表各个初始流量
    H0         sum of all transaction flows （初始流量总数）;

*Assignment for parameters
    Q0(i,j)=sam(i,j);
    H0=sum((i,j),sam(i,j));

display H0,sam;

Variables
    Q(i,j)       要调整的SAM表中的各个数值
    H            调整SAM表的总值
    Hratio       调整和原始两个总数的比例
    z            目标函数的数值即预期熵值;

*nonneg           每个变量必须是非负数
Positive variable Q(i,j);

equations
    totalsum         被调整的总数
    directentropy    目标函数 预期交叉熵
    balance          各个账户的平衡限制条件
    Hratiodef        Hratio 的定义和范围;

    totalsum..       H =e= sum((i,j),Q(i,j));
    Hratiodef..      Hratio =e= H/H0;
    directentropy..  z =e= sum((i,j)$sam(i,j),(1/H)*Q(i,j)*log(Q(i,j)/sam(i,j)))
                       -log(Hratio);
*    directentropy..  z =e= sum((i,j),(1/H)*Q(i,j)*log(Q(i,j)/sam(i,j)))-log(Hratio);
    balance(i)..     sum(j$sam(i,j),Q(i,j)) =e= sum(j,Q(j,i));

*对变量初始值赋值。
*这里对Hratio特别限制了范围。因为不限制的话，目标函数最小化时H会趋向无穷大
*得不到真实结果
    Q.l(i,j)=Q0(i,j);
    H.l=H0;
    Hratio.lo=0.5;
    Hratio.up=2;

model sambal  /all/;
option nlp=pathnlp;
solve sambal using nlp minimizing z;

display Q.l,H.l,Hratio.l;
*end 结束

感谢！

老师好，大家好！我用直接交叉熵的方法对14章的例14-7-2的Sam表调平，出现的下面的错误
**** Exec Error at line 82: log: FUNC SINGULAR: x = 0

谢谢了！

我用的代码是：
*第十四章例题7.2数据
set ac   /sec1,sec2soe,sec2pri,sec2fdi,sec3,com1,com2,com3,lab,cap,hhr,hhu,ent,
        gov,vatl,vatk,bustax,tariff,invsav,row,total/;
set a(ac)  /sec1,sec2soe,sec2pri,sec2fdi,sec3/;
set c(ac)  /com1,com2,com3/;
set f(ac)  /lab,cap/;
set hoh(ac)  /hhr,hhu/;
set vat(ac) /vatl,vatk/;
*set acnt(ac)  除了总数之外的所有账户;
set i(ac) /sec1,sec2soe,sec2pri,sec2fdi,sec3,com1,com2,com3,lab,cap,hhr,hhu,ent
        ,gov,vatl,vatk,bustax,tariff,invsav,row/;

*acnt(ac)=YES; acnt('total')=NO;
*i(ac)=YES; i('total')=NO;

alias (ac,acp),(a,ap),(c,cp),(f,fp),(hoh,hp),(vat,vatp);
alias (i,j);

table sam(ac,acp)
$INCLUDE C:\Users\phili\Desktop\GAMs\必读文件夹\张欣 第二版所有CGE模型的GAMS程序 解压密码 gzcbs2017\chap14QXSAM.inc
*display sa
*这里SAM表是从外面用txt文件引进
parameters
    Q0(i,j)    initial value SAM表各个初始流量
    H0         sum of all transaction flows （初始流量总数）;
*Assignment for parameters
    Q0(i,j)=sam(i,j);
    H0=sum((i,j),sam(i,j));
display H0,sam;

Variables
    Q(i,j)       要调整的SAM表中的各个数值
    H            调整SAM表的总值
    Hratio       调整和原始两个总数的比例
    z            目标函数的数值即预期熵值;

*nonneg           每个变量必须是非负数
Positive variable Q(i,j);

equations
    totalsum         被调整的总数
    directentropy    目标函数 预期交叉熵
    balance          各个账户的平衡限制条件
    Hratiodef        Hratio 的定义和范围;

    totalsum..       H =e= sum((i,j),Q(i,j));
    Hratiodef..      Hratio =e= H/H0;
    directentropy..  z =e= sum((i,j)$sam(i,j),(1/H)*Q(i,j)*log(Q(i,j)/sam(i,j)))
                       -log(Hratio);
*    directentropy..  z =e= sum((i,j),(1/H)*Q(i,j)*log(Q(i,j)/sam(i,j)))-log(Hratio);
    balance(i)..     sum(j$sam(i,j),Q(i,j)) =e= sum(j,Q(j,i));

*对变量初始值赋值。
*这里对Hratio特别限制了范围。因为不限制的话，目标函数最小化时H会趋向无穷大
*得不到真实结果
    Q.l(i,j)=Q0(i,j);
    H.l=H0;
    Hratio.lo=0.5;
    Hratio.up=2;

model sambal  /all/;
option nlp=pathnlp;
solve sambal using nlp minimizing z;

display Q.l,H.l,Hratio.l;
*end 结束

老师好，大家好！我用直接交叉熵法对书中例题14-7-2中的Sam表进行调平，却出现了如下的错误
**** Exec Error at line 82: log: FUNC SINGULAR: x = 0
请问这个问题怎么解决呢？
我用的代码在附件中。

学习中，谢谢楼主

各位同学好，我想请教一个问题，如果我想要把劳动力部门拆分成各个独立的劳动子部门，且要素市场独立，将QLD（a）进一步用ces函数来刻画，取要素为各劳动力子市场的投入值，除了生产函数/价格比值等式/价格关系等式以外，再加上劳动力子市场供给外生以外，还需要添加其他的约束条件吗？

各位老师同学好，大家有没有遇到出现optimal ，瓦尔拉斯为0，一致性检验通过，但是齐次性检验不通过的情况？
可以加个qq吗？
QQ：3501180048

您好，目前 更新的网址也打不开，您能再更新一下吗，谢谢您

您好，目前的网址也打不开，您能再更新一下？谢谢您

好贴就点赞，一起拿积分

   1848896750

yesyes

请问张欣老师书中13.4.1的程序，在加入劳动禀赋增加后的影响中，发现walras不充分为0，这是为什么？

《可计算一般均衡模型的基本原理与编程》（第二版）是张欣教授针对可计算一般均衡(CGE)模型的学习和应用提供的一部重要参考书。此书的更新版本不仅修正了第一版中的错误，还通过加强变量符号和定义的统一性来帮助读者更系统地理解和学习CGE模型的基本原理。

书中重新组织并优化了CGE模型程序、演示与讲解的内容，使之更加连贯且详细，方便初学者逐步掌握。特别是在宏闭合(macro closure)方面，作者先从一个通用的程序入手，然后通过不同的限制条件来展示如何构成不同形式的宏观闭合，以加深读者对这一概念的理解。

新版还新增了关于CGE模型设计与修改、可靠性稳健性敏感性评估等内容，并提供多个模型和编程实例。这些实例包括最简单的CGE模型示例以及在多部门多群体情况下植入LES函数(线性支出系统)的模型等，有助于加深理解并实操演练。

CGE模型广泛应用于财政政策分析、国际贸易研究以及气候变化经济学等领域。例如，今年诺贝尔经济学奖得主之一威廉·诺德豪斯教授开发的DICE（Dynamic Integrated Climate-Economy）模型就是基于动态CGE原理构建而成，用于评估和预测气候变化对经济的影响及其应对策略。

学习和应用CGE建模确实需要掌握相当数量的专业知识，并且涉及到GAMS等编程语言的学习与使用。对于初学者而言，可能会感到有一定难度，但通过系统性地阅读书籍、参与讨论以及实践编写代码等方式，可以逐步克服这些困难并最终掌握这一重要模型的运用技巧。

为了更好地讨论CGE建模与编程过程中的具体问题，建议建立一个学习小组或在线论坛，让大家能够分享经验、交流心得并解决遇到的实际难题。同时，利用书中提供的丰富实例和程序进行实践操作也是提高技能的有效途径之一。

此文本由CAIE学术大模型生成，添加下方二维码，优先体验功能试用

我也发现结果不一样，你解决了嘛