求助，关于空间经济模型数值模拟问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

1. sigma=5;psi=0.94;Y0=0.7;T=1.2;miu1=0.083;miu2=0.088;miu3=0.094;alfa0=0.5;L=1;
   
2. alfa11=0.042;alfa12=0.062;alfa13=0.082;
   
3. alfa21=0.032;alfa22=0.042;alfa23=0.062;
   
4. alfa31=0.022;alfa32=0.032;alfa33=0.042;
   
5. 
   
6. lamdar1=linspace(0,1,100);
   
7. lamdar2=linspace(0,1,100);
   
8. lamdar3=linspace(0,1,100);
   
9. lamdar0=1-lamdar1-lamdar2-lamdar3;
   
10. lamdas1=linspace(0,1,100);
    
11. lamdas2=linspace(0,1,100);
    
12. lamdas3=linspace(0,1,100);
    
13. lamdas0=1-lamdas1-lamdas2-lamdas3;
    
14. lamdat1=linspace(0,1,100);
    
15. lamdat2=linspace(0,1,100);
    
16. lamdat3=linspace(0,1,100);
    
17. lamdat0=1-lamdat1-lamdat2-lamdat3;
    
18. 
    
19. len=length(lamdar1);
    
20. 
    
21. omegar1=1*ones(1,len);
    
22. omegar2=1*ones(1,len);
    
23. omegar3=1*ones(1,len);
    
24. omegas1=1*ones(1,len);
    
25. omegas2=1*ones(1,len);
    
26. omegas3=1*ones(1,len);
    
27. omegat1=1*ones(1,len);
    
28. omegat2=1*ones(1,len);
    
29. omegat3=1*ones(1,len);
    
30. 
    
31. omegar1old=0;
    
32. omegar2old=0;
    
33. omegar3old=0;
    
34. omegas1old=0;
    
35. omegas2old=0;
    
36. omegas3old=0;
    
37. omegat1old=0;
    
38. omegat2old=0;
    
39. omegat3old=0;
    
40. 
    
41. valr1=(omegar1-omegar1old)./omegar1>0.00001;
    
42. valr2=(omegar2-omegar2old)./omegar2>0.00001;
    
43. valr3=(omegar3-omegar3old)./omegar3>0.00001;
    
44. vals1=(omegas1-omegas1old)./omegas1>0.00001;
    
45. vals2=(omegas2-omegas2old)./omegas2>0.00001;
    
46. vals3=(omegas3-omegas3old)./omegas3>0.00001;
    
47. valt1=(omegat1-omegat1old)./omegat1>0.00001;
    
48. valt2=(omegat2-omegat2old)./omegat2>0.00001;
    
49. valt3=(omegat3-omegat3old)./omegat3>0.00001;
    
50. ii=1;
    
51. syms lnGr1 lnGr2 lnGr3
    
52. while(valr1(ii)&&valr2(ii)&&valr3(ii)&&...
    
53.         vals1(ii)&&vals2(ii)&&vals3(ii)&&...
    
54.         valt1(ii)&&valt2(ii)&&valt3(ii)||ii>length(valr1))
    
55. 
    
56. omegar1old=omegar1;
    
57. omegar2old=omegar2;
    
58. omegar3old=omegar3;
    
59. omegas1old=omegas1;
    
60. omegas2old=omegas2;
    
61. omegas3old=omegas3;
    
62. omegat1old=omegat1;
    
63. omegat2old=omegat2;
    
64. omegat3old=omegat3;
    
65. 
    
66. omegar0=-(1-lamdar1-lamdar2-lamdar3)^(psi-1);
    
67. omegas0=-(1-lamdas1-lamdas2-lamdas3)^(psi-1);
    
68. omegat0=-(1-lamdat1-lamdat2-lamdat3)^(psi-1);
    
69. 
    
70. Yr=(L.*omegar1.*lamdar1)+(L.*omegar2.*lamdar2)+(L.*omegar3.*lamdar3)+(L.*(lamdar0.^psi))./psi;
    
71. Ys=(L.*omegas1.*lamdas1)+(L.*omegas2.*lamdas2)+(L.*omegas3.*lamdas3)+(L.*(lamdas0.^psi))./psi;
    
72. Yt=(L.*omegat1.*lamdat1)+(L.*omegat2.*lamdat2)+(L.*omegat3.*lamdat3)+(L.*(lamdat0.^psi))./psi;
    
73. 
    
74. Er1=miu1.*(Yr-Y0)+(L.*lamdar1.*omegar1.*alfa11/alfa0)+(L.*lamdar2.*omegar2.*alfa12./alfa0)+(L.*lamdar3.*omegar3.*alfa13./alfa0);
    
75. Er2=miu1.*(Yr-Y0)+(L.*lamdar1.*omegar1.*alfa21/alfa0)+(L.*lamdar2.*omegar2.*alfa22./alfa0)+(L.*lamdar3.*omegar3.*alfa23./alfa0);
    
76. Er3=miu1.*(Yr-Y0)+(L.*lamdar1.*omegar1.*alfa31/alfa0)+(L.*lamdar2.*omegar2.*alfa32./alfa0)+(L.*lamdar3.*omegar3.*alfa33./alfa0);
    
77. Es1=miu1.*(Yr-Y0)+(L.*lamdas1.*omegas1.*alfa11/alfa0)+(L.*lamdas2.*omegas2.*alfa12./alfa0)+(L.*lamdas3.*omegas3.*alfa13./alfa0);
    
78. Es2=miu1.*(Yr-Y0)+(L.*lamdas1.*omegas1.*alfa21/alfa0)+(L.*lamdas2.*omegas2.*alfa22./alfa0)+(L.*lamdas3.*omegas3.*alfa23./alfa0);
    
79. Es3=miu1.*(Yr-Y0)+(L.*lamdas1.*omegas1.*alfa31/alfa0)+(L.*lamdas2.*omegas2.*alfa32./alfa0)+(L.*lamdas3.*omegas3.*alfa33./alfa0);
    
80. Et1=miu1.*(Yr-Y0)+(L.*lamdat1.*omegat1.*alfa11/alfa0)+(L.*lamdat2.*omegat2.*alfa12./alfa0)+(L.*lamdat3.*omegat3.*alfa13./alfa0);
    
81. Et2=miu1.*(Yr-Y0)+(L.*lamdat1.*omegat1.*alfa21/alfa0)+(L.*lamdat2.*omegat2.*alfa22./alfa0)+(L.*lamdat3.*omegat3.*alfa23./alfa0);
    
82. Et3=miu1.*(Yr-Y0)+(L.*lamdat1.*omegat1.*alfa31/alfa0)+(L.*lamdat2.*omegat2.*alfa32./alfa0)+(L.*lamdat3.*omegat3.*alfa33./alfa0);
    
83. 
    
84. A1=L.*lamdar1.*(T.^(1-sigma)).*((omegar1).^(1-sigma.*alfa0))+L.*lamdas1.*(T.^(1-sigma)).*((omegas1).^(1-sigma.*alfa0))+(L.*lamdat1.*(T.^(1-sigma)).*((omegat1).^(1-sigma.*alfa0)));
    
85. A2=L.*lamdar2.*(T.^(1-sigma)).*((omegar2).^(1-sigma.*alfa0))+L.*lamdas2.*(T.^(1-sigma)).*((omegas2).^(1-sigma.*alfa0))+(L.*lamdat2.*(T.^(1-sigma)).*((omegat2).^(1-sigma.*alfa0)));
    
86. A3=L.*lamdar3.*(T.^(1-sigma)).*((omegar3).^(1-sigma.*alfa0))+L.*lamdas3.*(T.^(1-sigma)).*((omegas3).^(1-sigma.*alfa0))+(L.*lamdat3.*(T.^(1-sigma)).*((omegat3).^(1-sigma.*alfa0)));
    
87. lnA1=log(A1);lnA2=log(A2);lnA3=log(A3);
    
88. eqn1=((1-sigma)/sigma)*lnGr1-alfa11*lnGr1-alfa21*lnGr2-alfa31*lnGr3-lnA1==0;
    
89. eqn2=((1-sigma)/sigma)*lnGr2-alfa12*lnGr1-alfa22*lnGr2-alfa32*lnGr3-lnA2==0;
    
90. eqn3=((1-sigma)/sigma)*lnGr3-alfa13*lnGr1-alfa23*lnGr2-alfa33*lnGr3-lnA3==0;
    
91. [lnGr1,lnGr2,lnGr3]=solve([eqn1,eqn2,eqn3],[lnGr1, lnGr2, lnGr3]);
    
92. Gr1=exp(lnGr1);Gr2=exp(lnGr2);Gr3=exp(lnGr3);
    
93. Gs1=Gr1;Gt1=Gr1;Gs2=Gr2;Gt2=Gr2;Gs3=Gr3;Gt3=Gr3;
    
94. 
    
95. MGr1=(Gr1.^alfa11).*(Gr2.^alfa21).*(Gr3.^alfa31);
    
96. MGr2=(Gr1.^alfa12).*(Gr2.^alfa22).*(Gr3.^alfa32);
    
97. MGr3=(Gr1.^alfa13).*(Gr2.^alfa23).*(Gr3.^alfa33);
    
98. MGs1=(Gs1.^alfa11).*(Gs2.^alfa21).*(Gs3.^alfa31);
    
99. MGs2=(Gs1.^alfa12).*(Gs2.^alfa22).*(Gs3.^alfa32);
    
100. MGs3=(Gs1.^alfa13).*(Gs2.^alfa23).*(Gs3.^alfa33);
     
101. MGt1=(Gt1.^alfa11).*(Gt2.^alfa21).*(Gt3.^alfa31);
     
102. MGt2=(Gt1.^alfa12).*(Gt2.^alfa22).*(Gt3.^alfa32);
     
103. MGt3=(Gt1.^alfa13).*(Gt2.^alfa23).*(Gt3.^alfa33);
     
104. 
     
105. omegar1=((((alfa0.*((Er1.*(T./Gr1).^(1-sigma))+(Es1.*(T./Gs1).^(1-sigma))+(Et1.*(T./Gt1).^(1-sigma)))).^(1/sigma))./(MGr1))).^(1/alfa0);
     
106. omegar2=((((alfa0.*((Er2.*(T./Gr1).^(1-sigma))+(Es2.*(T./Gs1).^(1-sigma))+(Et2.*(T./Gt1).^(1-sigma)))).^(1/sigma))./(MGr2))).^(1/alfa0);
     
107. omegar3=((((alfa0.*((Er3.*(T./Gr1).^(1-sigma))+(Es3.*(T./Gs1).^(1-sigma))+(Et3.*(T./Gt1).^(1-sigma)))).^(1/sigma))./(MGr3))).^(1/alfa0);
     
108. omegas1=((((alfa0.*((Er1.*(T./Gr1).^(1-sigma))+(Es1.*(T./Gs1).^(1-sigma))+(Et1.*(T./Gt1).^(1-sigma)))).^(1/sigma))./(MGs1))).^(1/alfa0);
     
109. omegas2=((((alfa0.*((Er2.*(T./Gr1).^(1-sigma))+(Es2.*(T./Gs1).^(1-sigma))+(Et2.*(T./Gt1).^(1-sigma)))).^(1/sigma))./(MGs2))).^(1/alfa0);
     
110. omegas3=((((alfa0.*((Er3.*(T./Gr1).^(1-sigma))+(Es3.*(T./Gs1).^(1-sigma))+(Et3.*(T./Gt1).^(1-sigma)))).^(1/sigma))./(MGs3))).^(1/alfa0);
     
111. omegat1=((((alfa0.*((Er1.*(T./Gr1).^(1-sigma))+(Es1.*(T./Gs1).^(1-sigma))+(Et1.*(T./Gt1).^(1-sigma)))).^(1/sigma))./(MGt1))).^(1/alfa0);
     
112. omegat2=((((alfa0.*((Er2.*(T./Gr1).^(1-sigma))+(Es2.*(T./Gs1).^(1-sigma))+(Et2.*(T./Gt1).^(1-sigma)))).^(1/sigma))./(MGt2))).^(1/alfa0);
     
113. omegat3=((((alfa0.*((Er3.*(T./Gr1).^(1-sigma))+(Es3.*(T./Gs1).^(1-sigma))+(Et3.*(T./Gt1).^(1-sigma)))).^(1/sigma))./(MGt3))).^(1/alfa0);
     
114. 
     
115. valr1=(omegar1-omegar1old)./omegar1>0.00001;
     
116. valr2=(omegar2-omegar2old)./omegar2>0.00001;
     
117. valr3=(omegar3-omegar3old)./omegar3>0.00001;
     
118. vals1=(omegas1-omegas1old)./omegas1>0.00001;
     
119. vals2=(omegas2-omegas2old)./omegas2>0.00001;
     
120. vals3=(omegas3-omegas3old)./omegas3>0.00001;
     
121. valt1=(omegat1-omegat1old)./omegat1>0.00001;
     
122. valt2=(omegat2-omegat2old)./omegat2>0.00001;
     
123. valt3=(omegat3-omegat3old)./omegat3>0.00001;
     
124. ii=ii+1;
     
125. 
     
126. end
     
127. L=L+0.1;
     
128. plot(lamdar1,omegar1)

复制代码

我参照论坛里发的关于空间经济核心边缘模型的数值模拟代码进行修改，推广到三个地区、三种制造业的情况，修改出的程序如上述，目的是想绘制出参数L对lamda的影响，但是代码运行却没有结果，想知道问题出在哪里

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：经济模型 数值模拟 空间经济 length omega 空间经济 数值模拟 核心边缘模型 产业扩散模型 克鲁格曼