《可计算一般均衡模型的基本原理与编程》的CGE模型讨论

张欣genechang 发表于 2011-9-10 23:38
有同学照抄书中的GAMS程序仍然不能运行，很久也没有找到其中的差错。这里有书中第13章189-201页中的GAMS程序 ...

谢谢，有没有第十四章的程序，我照抄都不运行，郁闷啊

请教：张老师、各位同学，我的CGE模型校调正确运行，但是模拟部分不执行，出现以下错误：rPOWER:FUNC DOMAIN:X**Y,X=0,Y<0.是什么原因，我查询SAM表，没有X=0的情况，是什么原因造成的

我心飞翔2012 发表于 2012-5-19 11:04
谢谢，有没有第十四章的程序，我照抄都不运行，郁闷啊

你要用英文的txt文字编辑软件打印。注意字母l 和 o 与数字 1 和 0 的区别。然后，要把文章中的SAM表和文件路径修改成符合你电脑中的实际情况。
我来试试看能不能上载这个文件。

这里是第14章p.217-p.228的程序。运行前，仍然要
1、另行制作文章中的SAM表，
2、将程序中文件"$INCLUDE..."的路径修改成符合你电脑中SAM表存储的实际情况。

我心飞翔2012 发表于 2012-5-26 14:43
请教：张老师、各位同学，我的CGE模型校调正确运行，但是模拟部分不执行，出现以下错误：rPOWER:FUNC DOMAI ...

它认为你一个幂函数的底数成零，又是幂小于零，这是GAMS优化程序不允许的。
是不是SAM表的空格被程序认作零？

请教：我利用张老师教材上的代码，编制产业结构的程序，能运行一部分，但是嵌套函数出现如下错误：**** ERRORS/WARNINGS IN EQUATION QAfn(RACT2)
     24 error(s): vcPower: FUNC DOMAIN: x**c, x < 0 (RETURNED        0)

**** ERRORS/WARNINGS IN EQUATION QAfn(RACT3)
     12 error(s): vcPower: FUNC DOMAIN: x**c, x < 0 (RETURNED        0)

**** ERRORS/WARNINGS IN EQUATION QAfn(RACT4)
     4 error(s): vcPower: FUNC DOMAIN: x**c, x < 0 (RETURNED        0)

**** ERRORS/WARNINGS IN EQUATION QAfn(RACT5)
     36 error(s): vcPower: FUNC DOMAIN: x**c, x < 0 (RETURNED        0)

**** ERRORS/WARNINGS IN EQUATION QAfn(RACT6)
     44 error(s): vcPower: FUNC DOMAIN: x**c, x < 0 (RETURNED        0)
等等，我检查SAM表和数据，没有发现X<0的情况，不知道是什么原因造成的？

我心飞翔2012 发表于 2012-6-13 15:31
请教：我利用张老师教材上的代码，编制产业结构的程序，能运行一部分，但是嵌套函数出现如下错误：**** ERR ...

可能程序在计算或者优化运行过程中，X到了负值。可以试试在程序中加个限制条件指令
Positive variable X(..)
注意先把这个X找出来，确定是哪一个变量出了问题，然后正确写出该变量字母。
其它同学对此有没有idea和办法?

张欣genechang 发表于 2012-6-14 12:46
可能程序在计算或者优化运行过程中，X到了负值。可以试试在程序中加个限制条件指令
Positive variable X ...

谢谢张老师，我把程序分成多个模块进行，把方程调平衡，再核查吧

张老师，我将程序各个模块进行检验，发现方程都是平衡的，但是一给政策冲击，就出现三组方程不平衡，他们是：QAFOC(a)..
PVA(a)/((1+ltbus(a)+lgfsa(a)+cgfsa(a)+ctbus(a))*PINTA(a))=e=(deltaAa(a)/(1-deltaAa(a)))*(QINTA(a)/QVA(a))**(1-rhoAa(a));

PAeq(a)..
PA(a)*QA(a)=e=(1+ltbus(a)+lgfsa(a)+cgfsa(a)+ctbus(a))*PINTA(a)*QINTA(a)+PVA(a)*QVA(a);
YGLeq..
YGL=e=sum(h,tihl(h)*YH(h))+tilEnt*YENT+csc*CSCS0+sum(a,lgfsa(a)*PINTA(a)*QINTA(a))+
sum(a,ltbus(a)*PINTA(a)*QINTA(a))+sum(a,ltval(a)*WL*QLD(a))+sum(a,ltvak(a)*WK*QKD(a));
现在很困惑，不知道是什么原因？

正在学习中，觉得挺复杂的。希望大家啊能分享一些有用的资料。