请问我的模型出现了Blanchard Kahn conditions are not satisfied: indeterminacy due to rank failure，用判错命令告诉我是无明显错误，所以我判定是出现了[size=15.008px]timing error，但我自己目前看不出具体错误出现在哪，求助论[size=15.008px]坛里的各位！
model部分如下：
model;
%(1)
lab_e^phi=w_e/(c_e+pi*c_ge);

%(2)
j/h_e=1/(c_e+pi*c_ge)*((1+tau_p)*theta_e*p_h+(1-theta_e)*R_he)-beta_e*((1-theta_e)*R_he(+1)+theta_e*p_h(+1))/(c_e(+1)+pi*c_ge(+1));

%(3)
(c_e(+1)+pi*c_ge(+1))/(c_e+pi*c_ge)=beta_e*(r_ce(+1)+1-delt);

%(4)
(c_e(+1)+pi*c_ge(+1))/(c_e+pi*c_ge)=beta_e*(r_h(+1)+1-delt);

%(5)
lab_l^phi=w_l/(c_l+pi*c_gl);

%(6)
j/h_l=1/(c_l+pi*c_gl)*((1+tau_p)*theta_l*p_h+(1-theta_l)*R_hl)-beta_l*((1-theta_l)*R_hl(+1)+theta_l*p_h(+1))/(c_l(+1)+pi*c_gl(+1));

%(7)
(c_l(+1)+pi*c_gl(+1))/(c_l+pi*c_gl)=beta_l*(r_cl(+1)+1-delt);

%(8)
y_c=a*k_ce^(alph*alph_c)*k_cl^(alph_c*(1-alph))*lab_ce^((1-alph_c)*alph)*lab_cl^((1-alph_c)*(1-alph));

%(9)
r_ce=alph_c*alph*y_c/k_ce;

%(10)
r_cl=alph_c*(1-alph)*y_c/k_cl;

%(11)
w_e=(1-alph_c)*alph*y_c/lab_ce;

%(12)
w_l=(1-alph_c)*(1-alph)*y_c/lab_cl;

%(13)
y_h=L*k_h^alph_h*lab_he^((1-alph_h)*alph)*lab_hl^((1-alph_h)*(1-alph));

%(14)
p_h=xi/(xi-1)*mc;

%(15)
r_h=mc*alph_h*y_h/k_h;

%(16)
lab_hl/lab_cl=gamma_s;

%(17)
c_e/c_l= epsilo;

%(18)
d_h=p_h*y_h-mc*y_h;

%(19)
y_c=c_e+c_l+c_ge+c_gl+omega2*y_c;

%(20)
lab_e=lab_ce+lab_he;

%(21)
lab_l=lab_cl+lab_hl;

%(22)
R_he=gamm*w_e;

%(23)
R_hl=gamm*w_l;

%(24)
k_ce=(1-delt)*k_ce(-1)+dinv_ce;

%(25)
k_cl=(1-delt)*k_cl(-1)+dinv_cl;

%(26)
k_h=(1-delt)*k_h(-1)+dinv_he;

%(27)
c_ge+c_gl=omeg*(y_c+y_h);

%(28)
b=v*(y_c+y_h);

%(29)
log(j)=rho_j*log(j(-1))+ej;

%(30)
log(a)=rho_a*log(a(-1))+ea;

%(31)
log(L)=rho_L*log(L(-1))+eL;

%(32)
log(v)=rho_v*log(v(-1))+ev;

%(33)
y_c/y_h=mu;

%(34)
lab_he/lab_ce=gamma_s;

%(35)
lab_e/lab_l=eta;