dynare报错

var c_1  pi w_1 n_1 c_2 r lambda_2 w_2 n_2 Q h_1 h_2 b_2 lev_2 y k b_e q c_e lambda_e lev_e I x   e_a e_j e_2 e_e e_r;

varexo es_a es_j es_2 es_e es_r;

parameters eta beta_1 beta_2 beta_e m2_ss me_ss thet phi delt kappa alpha omega1 omega2 omega3 omega4 omega5 omega6 omega7
omega8 omega9 omega10 x_ss j omega11 omega12  rho_1 rho_pi rho_y rho_a rho_j rho_me rho_m2 rho_r sigma_a sigma_j sigma_me sigma_m2 sigma_r;

beta_1 = 0.09926;   //储蓄家庭主观贴现因子
beta_2 = 0.095;     //借贷家庭主观贴现因子
beta_e = 0.96;      //企业家主观贴现因子
eta = 0.5;          //劳动供给弹性
m2_ss = 0.7231;
me_ss = 0.8266;
thet = 0.75;
phi = 0.64;
delt = 0.025;
kappa = 10;
alpha = 0.33;
omega1 = 0.5192;    //储蓄家庭消费占比
omega2 = 0.3037;    //借贷家庭消费占比
omega3 = 0.1741;    //企业家消费占比
omega4 = 0.3022;    //投资占GDP比重
omega5 = 1-omega4;   //消费占GDP比重
omega6 = omega1*omega5; //储蓄家庭占GDP比重
omega7 = omega2*omega5; //借贷家庭占GDP比重
omega8 = omega3*omega5; //企业家消费占GDP比重
omega9 = omega4/delt; //资本存量占GDP比重
omega10 = 1-beta_2-(beta_1-beta_2)*m2_ss; //自定义参数
x_ss = (1/omega9)*(alpha*beta_e)*(1/(1-beta_e*(1-delt)-(1-delt)*(beta_1-beta_e)*me_ss));
j = (((((1-phi)*(1-alpha)/x_ss)/omega2)-1)*omega10)/(1-beta_1);
omega11 = j*omega7/omega10; //借贷家庭房地产占GDP比重
omega12 = j*omega6/(1-beta_1); //储蓄家庭房地产占GDP比重
rho_1 = 0.75;      //货币政策平滑系数
rho_pi = 1.6;      //利率对通胀的反应系数
rho_y = 0.63;      //利率对产出的反应系数


//parameters of shocks
rho_a = 0.7806;
rho_j = 0.8106;
rho_me = 0.8333;
rho_m2 = 0.7723;
rho_r = 0.0637;
sigma_a = 0.0066;
sigma_j = 0.0288;
sigma_me = 0.1729;
sigma_m2 = 0.175;
sigma_r = 0.0433;
model(linear);
//储蓄家庭
c_1 = r+c_1(+1)-pi(+1);
w_1 = c_1+(1+eta)*n_1;
Q-c_1 = (1-beta_1)*(e_j-h_1)+beta_1*(Q(+1)-c_1(+1));

//借贷家庭
-c_2 = (beta_2/beta_1)*(r-c_2(+1)-pi(+1))+(1-(beta_2/beta_1))*lambda_2;
w_2 = c_2+(1+eta)*n_2;
(1/omega10)*(Q-c_2)=e_j-h_2+(beta_2/omega10)*(Q(+1)-c_2(+1))+(((beta_1-beta_2)*m2_ss)/omega10)*(lambda_2+e_2+Q(+1)+pi(+1)-r);
b_2 = Q(+1)+e_2+h_2+pi(+1)-r;
omega7*c_2+(m2_ss*omega11)*(r(-1)+b_2(-1)-pi)+omega11*(h_2-h_2(-1))=(m2_ss*omega11*beta_1)*b_2+(((1-phi)*(1-alpha))/x_ss)*(w_2+n_2);
lev_2 = b_2-Q-h_2;

//企业家
y = e_a + alpha*k(-1)+(1-alpha)*phi*n_1+(1-alpha)*(1-phi)*n_2;
b_e = e_e+q(+1)+k+pi(+1)-r;
-c_e = (beta_e/beta_1)*(r-c_e(+1)-pi(+1))+(1-(beta_e/beta_1))*lambda_e;
w_1 = y - x - n_1;
w_2 = y - x - n_2;
((alpha*beta_e)/(omega9*x_ss))*(y-c_e-k(-1)-x)+(1-delt)*beta_e*(q-c_e)+(1-delt)*(beta_1-beta_e)*me_ss*(lambda_e(-1)+e_e(-1)+q+pi-r(-1));
omega8*c_e+((phi*(1-alpha))/x_ss)*(w_1+n_1)+((1-phi)*(1-alpha)/x_ss)*(w_2+n_2)+me_ss*(1-delt)*omega9*(r(-1)+b_e(-1)-pi)+delt*omega9*(q+k)=(1/x_ss)*(y-x)+me_ss*beta_1*(1-delt)*omega9*b_e;
lev_e = b_e - q - k;

//零售商
pi=beta_1*pi(+1)-(((1-beta_1*thet)*(1-thet))/thet)*x;

//资本品生产商
k=(1-delt)*k(-1)+delt*I;

1=(1+q)*(1-kappa*(I-I(-1)))+(beta_1+beta_1/c_1(+1)-beta_1/c_1)*(I(+1)-I);


// 市场出清
omega12*h_1+omega11*h_2=0;
y = omega6*c_1+omega7*c_2+omega8*c_e+omega4*I;


//货币政策
r = rho_1*r(-1)+(1-rho_1)*(rho_pi*pi+rho_y*y)+e_r;

//外生冲击
e_a = rho_a*e_a(-1)+es_a;
e_j = rho_j*e_j(-1)+es_j;
e_2 = rho_m2 *e_2(-1)+ es_2;
e_e = rho_me*e_e(-1)-es_e;
e_r = rho_r*e_r(-1)+es_r;

end;
steady;
check;

shocks;
var es_a = sigma_a^2;
var es_j = sigma_j^2;
var es_2 = sigma_m2^2;
var es_e = sigma_me^2;
var es_r = sigma_r^2;
end;

stoch_simul(irf = 40) y c_1 c_2 c_e Q q lev_2 lev_e I ;


出现错误：
错误使用 print_info (line 65)
The Jacobian contains NaNs

出错 check (line 76)
    print_info(info, 0, options);

出错 x2 (line 378)
oo_.dr.eigval = check(M_,options_,oo_);

出错 dynare (line 180)
evalin('base',fname) ;

>> 各位帮帮忙

已解决，对数线性化时，资本生产商错误