张欣《可计算一般均衡模型的基本原理与编程》的课后习题答案

display Q.l parameter use2(i) /sec1 0 sec2 200 sec3 0/; display use2 equation commodityequi2(i); commodityequi2(i).. sum(j,a(i,j)*Q(j))+use2(i)=e=Q(i); model simiomodel /commodityequi2/ solve simiomodel using mcp; display 'Q.1,为求解的结果 ', Q.l
1
第四章 SAM 表
支出 1 活动 2 商品 市场销售 产出 中间 投入 增值 居民要素 收入 企业要素 收入 生产 税增 值税 销售税， 关税，出 口税 要素税， 政府要素 收入 直接 税，收 入税 居民储 蓄 进口 对国外要 素的支付 企业直接 税，企业向 政府缴纳盈 余 企业储蓄 企业向国外 支付盈余 政府储 蓄 政府对 国外的 居民之 间转移 支付 企业对居民 转移支付 政府对 居民转 移支付 政府对 企业的 转移支 付 交易成本 3 要素 4 居民 居民自 产自销 市场销 售的私 人消费 5 企业 6 政府 7 储蓄 —投 资账 户 8 国外 汇总 总产 出 政府消 费 投资 出口 国外要 素收入 国外对 居民转 移支付 国外对 企业的 转移支 付 国外对 政府的 转移支 付 国外净 储蓄 总需 求 要素 收入 居民 总收 入 企业 总收 入 政府 总收 入 总储 蓄 外汇 支出
张欣《可计算一般均衡模型的基本原理与编程》 的课后习题答案
第 2 章 投入产出表和投入产出模型
1、将下面投入产出表填满： 农业 制造业 中间投入合 计 370 460 最终使用 380 630 总产出 750 1090
160 210 农业 140 320 制造业 300 530 中间投入合 计 450 560 增加值 750 1090 总投入合计 计算过程： 第一行，160+210=370，750-370=380。 第二行，1090-630=460,460-140=320。 第三行，160+140=300,210+320=530 。 第五行，总投入=总产出 第四行，750-300=450,1090-530=560,450+560=1010 2.a 行模型为： 0.818 0.098 0.084 880 480 880 0.159 0.209 0.159 1530 990 1530 0.091 0.098 0.234 1070 590 1070
1 2 3 4
活动 商品 要素 居民
收 入
5
企业
பைடு நூலகம்
6
政府 储蓄 —投 资账 户 国外
7
8
汇总
总投 入
总供应
要素支出
居民支 出
企业支出
支付 政府支 出
总投 资
外汇收 入
注：表中标注黄色的为习题上要求归类的变量
第五章 SAM 表平衡
1.(a)最小二乘法平衡的结果是： 470 160 150 90 140 320 170 910 80 150 250 610 Q 320 360 430 170 560 150 1110 880 GAMS 程序为： $title===least square balance method===== set ac /sec1,sec2,sec3,lab,cap,hh,total/; set i(ac) /sec1,sec2,sec3,lab,cap,hh/; alias(ac,acp); alias(i,j); table sam(*, *) sec1 sec2 sec3 lab cap hh total sec1 160 150 90 540 940 sec2 140 320 170 910 1540 sec3 80 150 250 610 1000 lab 320 360 400 1080 cap 170 550 150 870 hh 1050 860 1910 total 870 1530 1060 1050 860 2060 ; parameters Q0(i,j) initial value; Q0(i,j)=sam(i,j); variables Q(i,j) z;
parameter pincrease(ii) increase in price in various sectors ; pincrease(ii)=p.l(ii)-1 display p.l,pincrease; parameter p20(s) p2fx(ss) a2(s,ss) a22(ss) b2(g,ss) ; p20(s)=1; p2fx(ss)=1.15; a2(s,ss)=i0(s,ss)/i0("total",ss); a22(ss)=i0("sec1",ss)/i0("total",ss); b2(g,ss)=i0(g,ss)/i0("total",ss); display a2,b2,p2fx; variable p2(s); equation sim2priceequ(ss); sim2priceequ(ss).. sum(s,a2(s,ss)*p2(s))+a22(ss)*p2fx(ss)+sum(g,b2(g,ss)*w(g))=e=p2(ss); p2.l(s)=p20(s); model sim2pricemodel /sim2priceequ/; solve sim2pricemodel using mcp ; parameter pincrease2(ss) increase in price in various sectors ;
w0(g) use(i) lab(ii) cap(ii) q0(i) ; p0(i)=1; w0(g)=1; a(i,ii)=(p0(i)*i0(i,ii))/(p0(ii)*i0("total",ii)); b(g,ii)=(w0(g)*i0(g,ii))/(p0(ii)*i0("total",ii)); use(i)=i0(i,"finaluse")/p0(i); lab(ii)=i0("labor",ii)/w0("labor"); cap(ii)=i0("capital",ii)/w0("capital"); q0(i)=i0("total",i)/p0(i); variable p(i) w(g) q(i) ; display a,b,p0,w0,use,lab,cap,q0; equation priceequ(ii); priceequ(ii).. sum(i,a(i,ii)*p(i))+sum(g,b(g,ii)*w(g))=e=p(ii); p.l(i)=p0(i); w.fx(g)=w0(g); model iopricemodel /all/; solve iopricemodel using mcp ; display p.l,w.l; parameter wl1 labor price wk1 capital cost ; wl1=1; wk1=1.1; display wl1,wk1;
positive variables Q(i,j); equations sumsquare balance ; sumsquare.. z=e=sum((i,j)$sam(i,j),(Q(i,j)-sam(i,j))**2); balance(i).. sum(j,Q(i,j))=e=sum(j,Q(j,i)); Q.l(i,j)=Q0(i,j); model sambal /all/; solve sambal using nlp minimizing z; display Q.l; (b)假设居民消费不可靠，得到的结果为： 540 160 80 90 140 320 170 910 80 150 250 610 Q 320 360 430 170 630 150 1110 950 在程序中加入限制条件，只需改 equation 步： equations sumsquare balance item ; sumsquare.. z=e=sum((i,j)$sam(i,j),(Q(i,j)-sam(i,j))**2); balance(i).. sum(j,Q(i,j))=e=sum(j,Q(j,i)); item(i).. Q(i,'hh')=e=Q0(i,'hh'); Q.l(i,j)=Q0(i,j); (d)RAS 法平衡的结果是： sec1 sec1 sec2 sec3 lab cap hh sec2 sec3 99.163 188.524 490.299 990.397 231.617 2000.000 900.000 lab cap hh 616.125 1045.024 1238.851 total 920.000 1540.000 2000.000 2000.000 900.000 2900.000 tartot 920.000 1540.000 2000.000 2000.000 900.000 2900.000
equation sim1priceequ(ii); sim1priceequ(ii).. sum(i,a(i,ii)*p(i))+b("capital",ii)*wk1+b("labor",ii)*wl1=e=p(ii); p.l(i)=p0(i); model simpricemodel /sim1priceequ/; solve simpricemodel using mcp ;
pincrease2(ss)=p2.l(ss)-1 display p2.l,pincrease2; *end 1.投入产出模型为：
0.159 0.091 0.205 *1.1 0.040 1 0.182 1 ΔP (I A' ) A x wx 0.098 1 0.209 0.098 0.366 *1.1 0.055 0.195 1 0.234 0.084 0.140 *1.1 0.034 2.价格模型为 P A' P A t wt A k wk 计算结果显示：制造业价格增加 2.1%，服务业增加 2.1%。
1010
Q1 39.062 b.计算结果： Q2 286.618 Q3 41.301 程序
$title input-output model set ac /sec1,sec2,sec3,labor,finaluse,fotal/; set i(ac) /sec1,sec2,sec3/; alias(i,j); table IO(*,*) sec1 sec2 sec3 finaluse sec1 160 150 100 480 sec2 140 320 530 900 sec3 80 150 250 590 labor 320 350 410 capital 180 560 150 total 880 1530 1070 ; parameter int(i,j) use(i) lab(i) cap(i) a(i,j) Q0(i) ; int(i,j)=IO(i,j); use(i)=IO(i,"finaluse"); lab(i)=IO("labor",i); cap(i)=IO("capital",i); Q0(i)=IO("total",i); a(i,j)=int(i,j)/Q0(j); display int,use,lab,Q0,a; variable Q(i); equation commodityequi(i); commodityequi(i).. sum(j,a(i,j)*Q(j))+use(i)=e=Q(i); Q.l(i)=Q0(i); model IOmodel /all/; solve IOmodel using mcp ; Total 880 1530 1070
第三章 投入产出中的价格关系
程序： set ac /sec1,sec2,sec3,labor,capital,finaluse,total/; set i(ac) /sec1,sec2,sec3/; set g(ac) /labor,capital/; set j(ac) /sec1,sec2,sec3,finaluse,total/; set s(ac) /sec2,sec3/; alias(i,ii); alias(s,ss); table i0(*,*) sec1 sec2 sec3 finaluse sec1 160 150 90 480 sec2 140 320 170 900 sec3 80 150 250 590 labor 320 350 410 capital 180 560 150 total 880 1530 1070 ; parameter a(i,ii) b(g,ii) *int(i,ii) p0(i) Total 880 1530 1070