发帖
雷达卡
这是我的模型；
Y我是提前解出来的，所以我直接赋值。房地产总供给为1.
var c l d b c_1 l_1 d_1 b_1 y d_z b_z r w a b_z q c_z x pai h;
varexo etaa;
parameters beta_a beta_b beta_z varphi m_b m_z rhoa alpha omega1 omega2 omega3 omega4 omega5 D1 L1 C1 B1 D L C B Y DZ BZ CZ W sigmaa X Q R theta PAI v H rhor rhoy rhopai;
beta_a = 0.99;
beta_b = 0.98;
beta_z = 0.886;
varphi = 1.99;
m_b = 0.75;
m_z = 0.5;
v=0.7;
rhoa = 0.96;
alpha = 0.4;
theta = 0.75;
sigmaa = 0.01;
rhor = 0.8;
rhoy = 2.4;
rhopai = 0.9;
PAI=1;
R=1/beta_a;
X = 1.2;
omega1 = (1-beta_a);
omega2 = (1-(beta_a-beta_b)*m_b-beta_b);
omega3 = X/(1-alpha)*(1-(beta_a-beta_z)*m_z-beta_z);
omega4 = (1-v)*alpha/X/(1+(1-beta_a)*m_b/omega2);
omega5 = ((1-alpha)/X-(1-beta_a)*m_z/omega3);
H=1;
Y=1.35;
C1=omega4*Y;
CZ=omega5*Y;
C=Y-C1-CZ;
Q=(1/omega3+(1-omega4-omega5)/omega1+omega4/omega2)*Y;
DZ=Y/omega3/Q;
D1=C1/omega2/Q;
D=C/omega1/Q;
B1=m_b*beta_a*Q*D1;
BZ=m_z*beta_a*Q*DZ;
W=(v*alpha)^(v)*((1-v)*alpha)^(1-v)*(DZ/Y)^((1-alpha)/alpha)/X;
L=(W/C)^(1/varphi);
L1=(W/C1)^(1/varphi);
B=B1+BZ;
model;
(l)^(varphi)=w/c;
(1/d)-(q/c)+beta_a*q(+1)/c(+1)=0;
c(+1)/c=beta_a*r/pai(+1);
c_1+r(-1)*b_1(-1)/pai+q*(d_1-d_1(-1))=w*l_1+b_1;
b_1*r=m_b*q(+1)*d_1*pai(+1);
(l_1)^(varphi)=w/c_1;
(1/d_1)-(q/c_1)+(1/c_1/r)*m_b*q(+1)*pai(+1)-beta_b/c_1(+1)*m_b*q(+1)+beta_b*q(+1)/c_1(+1)=0;
y=a*(l)^(v*alpha)*(l_1)^((1-v)*alpha)*(d_z)^(0.6);
c_z+w*(l+l_1)+q*(d_z-d_z(-1))+r(-1)*b_z(-1)/pai=y/x+b_z;
b_z*r = m_z*q(+1)*d_z*pai(+1);
(((1-alpha)*y/d_z/x)-q)/c_z+(1/r/c_z-beta_z/c_z(+1)/pai(+1))*m_z*q(+1)*pai(+1)+beta_z*q(+1)/c_z(+1)=0;
v*alpha*y/x/l=w;
(1-v)*alpha*y/x/l_1=w;
pai=beta_a*pai(+1)-(1-theta)*(1-beta_a*theta)*x/theta;
r/R=(r(-1)/R)^(rhor)*((y/Y)^(rhoy)*(pai/PAI)^(rhopai))^(1-rhor);
y=c+c+1+c_z;
h=d+d_1+d_z;
b=b_1+b_z;
log(a)=rhoa*log(a(-1))-etaa;
end;
steady_state_model;
c=C;
d=D;
l=L;
b=B;
c_1=C1;
d_1=D1;
l_1=L1;
b_1=B1;
y=Y;
A=1;
r=R;
c_z=CZ;
d_z=DZ;
b_z=BZ;
q=Q;
x=X;
w=W;
pai=PAI;
h=H;
end;
shocks;
var etaa = sigmaa^2;
end;
stoch_simul(irf = 30, order = 1, periods = 0);
京公网安备 11010802022788号
