var y m  pi r u a e z ;
varexo eps_r eps_a eps_z eps_e;
parameters gamma  r_ss omega1 omega2 pi_ss psi rho_r rho_y rho_pi rho_u rho_a rho_e rho_z ;
  gamma= 0.16;
  omega1=0.6813;
  omega2=-0.3187;
  pi_ss=1.0122;
  r_ss=1.0544;
  psi=0.0034;

  rho_r=0.9144;
  rho_y=0.4312;
  rho_pi=0.3286;
  rho_u=0.0633;

  rho_a=0.24;
  rho_e=1.32;
  rho_z=0.209;

model(linear);
m=y-gamma/(r_ss-1)*r+e;//eq.7
y=y(+1)-omega1*(r-pi(+1))+omega2*(m(+1)-m)+1/(gamma-1)*omega2*(e(+1)-e)+omega1*(a(+1)-a); //eq.8
pi=pi_ss/r_ss*pi(+1)+psi*(y/omega1+omega2*m/omega1+omega2*e/omega1/(gamma-1)-a-z); //eq.9
u+m(-1)=m+pi;//eq.10
r=rho_r*r(-1)+(1-rho_r)*(rho_y*y(-1)+rho_pi*pi(-1)+rho_u*u(-1))+eps_r;//eq.11
a=rho_a*a(-1)+eps_a; //eq.12
e=rho_e*e(-1)+eps_e;//eq.13
z=rho_z*z(-1)+eps_z;//eq.14
end;

steady;
check;
resid;


shocks;
var eps_a;
stderr 0.1485;
var eps_e;
stderr 0.0634;
var eps_z;
stderr 0.4986;
var eps_r;
stderr 0.0068;
end;




stoch_simul(periods=1200,order=1,irf=60) y pi ;



就是想做一个脉冲响应分析，上述是代码，模型已经对数线性化，不知道哪里出了错，显示BK条件不满足，求助