用Matlab画经典RBC模型的相位图

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

% This Matlab m file actually plot the phase diagram for classical RBC model.
% This code is written by Xiangyang Li, Ph.D in SUFE, School of Finance
% xli21@nd.edu
clear;
%set in quarterly frequency
alpha=.36;
beta =.99;
sigma =2;
delta=.025;
%steady states and the golden capital stock.
ks = (alpha/(1/beta - 1 + delta))^(1/(1-alpha));
kgold = (alpha/delta)^(1/(1-alpha));
cs = ks^alpha - delta*ks;
%draw the isocline
kv=0:.1:200;
cv=0:.1:6;
% the capital isocline
plot(kv,kv.^alpha-delta*kv,'-r','LineWidth',2);
hold on;
% the consumption isocline
plot(ks*ones(size(cv)),cv,'-b','LineWidth',2);
plot(kgold*ones(size(cv)),cv,'-g');
%1. suppose we start from the point c_t = 2 and k_t =20, to see what happend
T = 52;
C=zeros(1,T);
K=zeros(1,T);
C(1)=2;
K(1) = 20;
for ii = 1:T
    %resources constraint
    K(ii+1) = K(ii)^alpha + (1 -delta)*K(ii)  - C(ii);
    %FOC of the problem
    C(ii+1) =( C(ii)^(-sigma)/(alpha*K(ii+1)^(alpha -1) + (1-delta) ))^(-1/sigma);
end
plot(K,C,'-k');
%2. suppose we start from the point c_t = 1 and k_t =20, to see what happend
T = 120;
C=zeros(1,T);
K=zeros(1,T);
C(1)=1;
K(1) = 20;
for ii = 1:T
    K(ii+1) = K(ii)^alpha + (1 -delta)*K(ii)  - C(ii);
    C(ii+1) =( C(ii)^(-sigma)/(alpha*K(ii+1)^(alpha -1) + (1-delta) ))^(-1/sigma);
end
plot(K,C,'-m');
%3. suppose we start from the point c_t = 4.5 and k_t =70, to see what happend
T = 70;
C=zeros(1,T);
K=zeros(1,T);
C(1)=4.5;
K(1) = 120;
for ii = 1:T
    K(ii+1) = K(ii)^alpha + (1 -delta)*K(ii)  - C(ii);
    C(ii+1) =( C(ii)^(-sigma)/(alpha*K(ii+1)^(alpha -1) + (1-delta) ))^(-1/sigma);
end
plot(K,C,'-c');
%4. suppose we start from the point c_t = 3.5 and k_t =90, to see what happend
T = 370;
C=zeros(1,T);
K=zeros(1,T);
C(1)=3.1;
K(1) = 80;
for ii = 1:T
    K(ii+1) = K(ii)^alpha + (1 -delta)*K(ii)  - C(ii);
    C(ii+1) =( C(ii)^(-sigma)/(alpha*K(ii+1)^(alpha -1) + (1-delta) ))^(-1/sigma);
end
plot(K,C,'-r');

%5. plot the policy function, policy function coefficient is obtained
%by solving the RBC model.
pol=0.462886778502290;
K=10:0.1:110;
C=pol*K*cs/ks + cs*(1-pol);
plot(K,C,'-k','LineWidth',2);
hold off;
legend('captial isocline', 'consumption isocline', ...
    'capital gold','explosive path 1', 'explosive path 2',...
'explosive path 3', 'explosive path 4','Saddle Path');

扫码加我 拉你入群

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

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

分享2 收藏17 回帖

关键词：MATLAB matla atlab Mat Atl 经典 模型 RBC

thank you!

非常感谢分享

谢谢楼主，太精彩了

真是很厉害哦！