在做毕业论文的过程中,用dynare运行代码,遇到这个问题该怎么解决。我做的应该算是RBC模型~
Residuals of the static equations:
Equation number 1 : NaN : M
Equation number 2 : NaN : K
Equation number 3 : NaN : 3
Equation number 4 : NaN : 4
Equation number 5 : NaN : 5
Equation number 6 : -0.00045517 : rb
Equation number 7 : NaN : 7
Equation number 8 : NaN : Y
Equation number 9 : NaN : 9
Equation number 10 : NaN : W
Equation number 11 : NaN : R
Equation number 12 : NaN : D
Equation number 13 : NaN : B
Equation number 14 : NaN : 14
错误使用 print_info (line 32)
The steady state has NaNs or Inf.
出错 steady (line 102)
print_info(info,options_.noprint, options_);
出错 RBC.driver (line 289)
steady;
以下是源代码
var M DM C I Y A K L W R P D B h;
varexo e_a e_h;
parameters alpha beta delta rhoa rhoh sigma eta theta
phi phid m rb rd;
alpha=0.45;
beta=0.9926;
delta=0.025;
rhoa=0.9;
rhoh=0.9;
sigma=2;
eta=0.1;
theta=0.75;
rs=1.008375;
rb=1.007;
rd=1.01225;
phi=1.5;
phid=0.0037;
m=1;
model;
//家庭
M=DM+(m*eta^2)/2;
//capital accumulation;
K=(1-delta)*K(-1)+I;
C^(sigma)*L^(sigma)=W/P;
(C(+1)/C)^sigma=beta*(1-delta+R(+1)/P(+1));
(2*M-m*eta^(2))=C^(-sigma)/(2*P);
rb=((P(+1)/P(+2))*(C(+1)/C)^(sigma))/beta;
log(h)=rhoh*log(h(-1))+e_h;
//厂商
Y=A*K^(alpha)*L^(1-alpha);
//technology shock
log(A)=rhoa*log(A(-1))+e_a;
//3 labor demand劳动需求方程
W=(1-alpha)*A*K(-1)^(alpha)*L^(-alpha);
//4 capital accumulation
R=alpha*A*K(-1)^(alpha-1)*L^(1-alpha);
//商业银行
D=W*L+R*K;
B=D;
//市场出清
Y = C+I+phid*D;
end;
initval;
A=1;
L=1/3;
C=0.76;
K=9;
e_h=0;
e_a=0;
h=0.2;
end;
shocks;
var e_h=.01^2;
end;
steady;
check;
stoch_simul(order=1);
出错 dynare (line 278)
evalin('base',[fname '.driver']);