勘误和解释：《动态随机一般均衡(DSGE)模型：理论、方法和Dynare实践》2020.11第二次

ahnulxy 发表于 2020-1-7 02:54
《动态随机一般均衡(DSGE)模型：理论、方法和Dynare实践》一书自2018年12月出版以来，已经一年有余，收到 ...

tulipsliu 发表于 2020-1-17 21:21
已下载，会传播。李向阳教授的书籍我也买了的。

非常棒的书籍

感谢支持！

higgspaul1 发表于 2020-6-13 16:10
请问楼主找到Monetary policy, inflation and the Business Cycle的课后题答案了么  有偿求一个

这本书目前没标准答案，不过有人正在组织准备弄

维柠a 发表于 2020-7-5 11:04
李老师您好
我是一名硕士研究生，我购买了您的书籍《动态随机一般均衡模型》，在学习的过程中遇到了一些问 ...

论坛里就有！你搜一下。或扫描本书封底二维码也能看到

tulipsliu 发表于 2020-4-2 18:44
我们买的纸质版的书籍。这个是中文的书籍，暂时没有PDF 版本的吧。
150块钱，你想学买一本吧。

某多多不香吗，40多就可以了，盗版20多元。

飞马洪 发表于 2021-2-21 09:11
某多多不香吗，40多就可以了，盗版20多元。

李向阳老师的教材，呕心沥血；

是别人的科研成果；作者拥有微软软件工程师资格。虽然教材里，编程在我看来是入门级的。
我也会 Micro soft visual studio.

现在电脑里的 MSVC ,除了可以编译 C++ 语言，已经集成了  Intel 公司的  FORTRAN 语言到 Micro soft Visual studio 里；

作者的编程功底强大， C++ 语言的编程，在程序员里也是非常难的。我们知道现在的动态语言 Python 非常容易入手。

但是，专业的程序员，最难编程的是 C++ ，我已经在论坛发布一些 R 混编 C++ 的工具包。比如  rugarch  package  ,做 GARCH  模型的R 工具包，也混编 C ++ ,

国际上最前沿的  DSGE ，也有混编 C++ ，以及 FORTRAN 语言。

我已经在论坛分享了程序，  MATLAB 编译 FORTRAN 语言。通过MATLAB 的 mex 命令，编译了FORTRAN 程序，生成为 MATLAB 可以执行的 .mexw64 （我的电脑系统和 MATLAB 均为 64-bit）。

作者的教材非常好。不但帮助国内的学生从基础的 MATLAB 语法以及 Dynare 的应用入手，我可以用这个句话来说：作者用心良苦。

即使这样的基础教学，没有五年以上的教学带，是无法完成的。
我打开MATLAB ， 展示我的软件配置吧。

作者不但有可以从学生的基础入门编程（在我看来教材的编程非常入门，与作者的微软软件工程车资格不相符）。也有经济学国际上各种经济现象研究的深度理解。 是作者多年经济学研究呕心沥血之作。

这样的好教材， 150.00 元购买正版，一点也不亏。

飞马洪 发表于 2021-2-21 09:11
某多多不香吗，40多就可以了，盗版20多元。

这里，我先回复里，我的电脑里MATLAB 编译FORTRAN的一段代码。
下一帖我不回复你，你自己看我在这里的留言，一个 MATLAB 可以编译的 FORTRAN 语言代码， MEX 文件不论是 C++ 语言的，还是 FORTRAN 语言的，funcion 前面必须有关键字 ， mex
比如计算DSER的有一个叫  valueiterationfunction , 那么 FORTRAN 语言要被 MATLAB 能 mex 编译的话，必须有 mex 关键字开头。

mexValueIterationFunction.F90

我的编译结果如下：
>> mex policies_on_itc_mex.F90
Building with 'Intel Parallel Studio XE 2018 for Fortran with Microsoft Visual Studio 2017'.
MEX completed successfully.
>> dir

.                           myinterp4_mex.F90           policies_on_mex.F90         spcoefs_mex.F90            
..                          myinterp4_mex.mexw64        policies_on_mex.mexw64      spcoefs_mex.mexw64         
myinterp3_mex.F90           policies_on_itc_mex.F90     ppval_mex.F90               
myinterp3_mex.mexw64        policies_on_itc_mex.mexw64  ppval_mex.mexw64            

>>

一段 MATLAB 可以编译的高级程序语言代码；

文件的后缀名是 .F90 格式。文件名 ： policies_on_itc_mex.F90

1. !
   
2. !   myinterp3_mex.F90
   
3. !   
   
4. !
   
5. !   Created by Volker Tjaden on 23.08.2013.
   
6. !   Copyright 2013 __MyCompanyName__. All rights reserved.
   
7. !
   
8. #include "fintrf.h"
   
9. 
   
10. !-----------------------------------------------------------------------
    
11. ! Gateway routine
    
12. !-----------------------------------------------------------------------
    
13. 
    
14. subroutine mexFunction(nlhs,plhs,nrhs,prhs)
    
15.             implicit none
    
16.         ! Arguments in mexFunction
    
17.             mwPointer  :: plhs(*), prhs(*)
    
18.             integer    :: nlhs, nrhs
    
19. 
    
20.             ! Function declarations
    
21.             mwPointer   ::  mxGetPr, mxCreateNumericMatrix
    
22.             integer(4)  ::  mxClassIDFromClassName
    
23.             integer(4)  ::  mxGetString
    
24.             mwPointer   ::  mxGetN , mxGetM
    
25. 
    
26.             ! Pointers to input/output mxArrays:
    
27.             ! Input
    
28.             mwPointer  :: prices,EV,EV_down,EV_up,kgrid,q,rho,w,gam,del,b,beta,xibar
    
29.                 mwPointer  :: tcred,itc
    
30.             ! Output
    
31.             mwPointer  :: kopt,xihat,up,down,uncont,kopt_ind
    
32.        
    
33.             ! Array size information
    
34.             mwSize             :: nks, neps, mwOne
    
35.             integer(4)   :: classid
    
36.             integer(4)   :: complexflag
    
37. 
    
38. !-----------------------------------------------------------------------
    
39.         ! Copy pointers
    
40.         prices        = mxGetPr(prhs(1))
    
41.         EV        = mxGetPr(prhs(2))
    
42.         EV_down = mxGetPr(prhs(3))
    
43.         EV_up   = mxGetPr(prhs(4))
    
44.         kgrid        = mxGetPr(prhs(5))
    
45.         q        = mxGetPr(prhs(6))
    
46.         rho        = mxGetPr(prhs(7))
    
47.         w        = mxGetPr(prhs(8))
    
48.         gam        = mxGetPr(prhs(9))
    
49.         del        = mxGetPr(prhs(10))
    
50.         b        = mxGetPr(prhs(11))
    
51.         beta        = mxGetPr(prhs(12))
    
52.         xibar   = mxGetPr(prhs(13))
    
53.     tcred   = mxGetPr(prhs(14))
    
54.         itc  = mxGetPr(prhs(15))
    
55. 
    
56.         ! Retrieve size information
    
57.         nks        = mxGetM(prhs(2))
    
58.         neps    = mxGetN(prhs(2))       
    
59. 
    
60.         ! Create return arguments
    
61.         complexflag = 0  
    
62.         mwOne=1       
    
63. 
    
64.         classid = mxClassIDFromClassName('double')      
    
65.             plhs(1) = mxCreateNumericMatrix(mwOne,neps,classid, complexflag)
    
66.             kopt = mxGetPr(plhs(1))
    
67.         plhs(2) = mxCreateNumericMatrix(nks,neps,classid, complexflag)       
    
68.         xihat = mxGetPr(plhs(2))
    
69. 
    
70.         classid = mxClassIDFromClassName('logical')  
    
71.         plhs(3) = mxCreateNumericMatrix(nks,neps,classid, complexflag)
    
72.             up = mxGetPr(plhs(3))
    
73.         plhs(4) = mxCreateNumericMatrix(nks,neps,classid, complexflag)
    
74.             down = mxGetPr(plhs(4))
    
75.         plhs(5) = mxCreateNumericMatrix(nks,neps,classid, complexflag)
    
76.             uncont = mxGetPr(plhs(5))
    
77.         classid = mxClassIDFromClassName('double')  
    
78.           plhs(6) = mxCreateNumericMatrix(mwOne,neps,classid, complexflag)
    
79.             kopt_ind = mxGetPr(plhs(6))
    
80.        
    
81.         ! call subroutine
    
82.         call policies_on(%val(kopt),%val(xihat),%val(up),%val(down),&
    
83.                 %val(uncont),%val(kopt_ind),%val(EV),%val(EV_down),&
    
84.                 %val(EV_up),%val(kgrid),%val(q),%val(rho),%val(w),&
    
85.                 %val(gam),%val(del),%val(b),%val(beta),%val(xibar),&
    
86.                 %val(tcred),%val(itc),nks,neps)
    
87. 
    
88.         return
    
89. end subroutine
    
90. 
    
91. subroutine policies_on(kopt,xihat,up,down,uncont,kopt_ind,&
    
92.         EV,EV_down,EV_up,kgrid,q,rho,w,&
    
93.         gam,del,b,beta,xibar,tcred,tstate,nks,neps)
    
94.         implicit none
    
95.         real(8), intent(out)        :: kopt(neps)
    
96.         real(8), intent(out)        :: xihat(nks,neps)
    
97.         !logical*1, intent(out)        :: up(nks,neps)
    
98.         !logical*1, intent(out)        :: uncont(nks,neps)
    
99.         !logical*1, intent(out)        :: down(nks,neps)
    
100.         integer*1, intent(out)        :: up(nks,neps)
     
101.         integer*1, intent(out)        :: uncont(nks,neps)
     
102.         integer*1, intent(out)        :: down(nks,neps)
     
103.         real(8), intent(out)        :: kopt_ind(neps)
     
104.         real(8), intent(in)        :: EV(nks,neps)
     
105.         real(8), intent(in)        :: EV_down(nks,neps)
     
106.         real(8), intent(in)        :: EV_up(nks,neps)
     
107.         real(8), intent(in)        :: kgrid(nks)
     
108.         real(8), intent(in)        :: q,rho,w,gam,del,b,beta,xibar
     
109.         real(8), intent(in) :: tcred,tstate
     
110.         mwSize, intent(in)        :: nks,neps
     
111.         ! local variables       
     
112.         integer(4)                :: i1,ind,pos,i2
     
113.         real(8)                        :: E,EC,vbest,vcand
     
114.        
     
115.         ! Loop over idiosyncratic productivity states
     
116.         ind=1       
     
117.         do i1=1,neps
     
118.                 ! 1. Find capital choice
     
119.                 pos=0
     
120.                 vbest=-huge(vbest)
     
121.                 if (i1>1) then
     
122.                         ind=int(kopt_ind(i1-1),4)
     
123.                 end if
     
124.                 do while (ind<=nks)
     
125.                         vcand=-gam*kgrid(ind)*q*rho*(1.d0-tstate*tcred)&
     
126.                         +beta*EV(ind,i1)
     
127.                         if (vcand>vbest) then
     
128.                                 vbest=vcand
     
129.                                 pos=ind
     
130.                         else
     
131.                                 exit
     
132.                         end if
     
133.                         ind=ind+1
     
134.                 end do
     
135.                 E=vbest
     
136.                 kopt_ind(i1)=real(pos,8)
     
137.                 kopt(i1)=kgrid(pos)
     
138. 
     
139.                 ! 2. Check whether optimal choice is within constrained
     
140.                 !    set and if not, in what direction constrain
     
141.                 do i2=1,nks
     
142.                         ! (EC-repmat(E,mpar.nksim,1))./(rho*w);
     
143.                         if (((1.d0-del+b)/gam)*kgrid(i2)<kopt(i1)) then
     
144.                                 up(i2,i1)=1 !.true.
     
145.                                 EC=-(1.d0-del+b)*kgrid(i2)*q*rho*(1.d0-tstate*tcred)+&
     
146.                                         beta*EV_up(i2,i1)
     
147.                                 xihat(i2,i1)=-(EC-E)/(rho*w)
     
148.                                 xihat(i2,i1)=min(xibar,max(0.d0,xihat(i2,i1)))
     
149.                         elseif(((1.d0-del-b)/gam)*kgrid(i2)>kopt(i1)) then
     
150.                                 down(i2,i1)=1 !.true.
     
151.                                 EC=-(1.d0-del-b)*kgrid(i2)*q*rho*(1.d0-tstate*tcred)+&
     
152.                                         beta*EV_down(i2,i1)
     
153.                                 xihat(i2,i1)=-(EC-E)/(rho*w)
     
154.                                 xihat(i2,i1)=min(xibar,max(0.d0,xihat(i2,i1)))
     
155.                         else
     
156.                                 uncont(i2,i1)=1!.true.
     
157.                                 xihat(i2,i1)=0.d0
     
158.                         end if
     
159.                 end do
     
160.         end do                       
     
161. 
     
162. end subroutine

复制代码

支持