小白发代码，求怎么加入利率规则

var a y c i g n k w rk tau_k tau_w tau_c b r tw tk tc pi mc i_g k_g;
varexo epsilon_a epsilon_i_g epsilon_g epsilon_tau_w epsilon_tau_k epsilon_tau_c epsilon_r epsilon_tw epsilon_tk epsilon_tc;
parameters A Y C I G N K W RK TAU_K TAU_W TAU_C R B PI I_G K_G;
parameters alpha beta theta phi_n delta_k C_div_Y delta_kg eta;
PI       =1;
A        =1;
N        =1/3;
alpha    =0.54;
beta     =0.9804;
theta    =0.9;//可调整
TAU_K    =0.3;
TAU_W    =0.1;
TAU_C    =0.15;
delta_k  =0.1;
delta_kg =0.05;
eta      =0.06;
C_div_Y  =0.417;
I_G_div_Y=0.128;  
phi_n    =(1-N)/N*(1-alpha)/0.417;
R        =1/beta;
RK       =(R-1+delta_k)/(1-TAU_K);        
//Y        =((alpha*(1-TAU_K)/(R-1+delta_k))^alpha*N^(1-alpha))^(1/(1-alpha));//用N和参数解出Y
Y        =((alpha*(1-TAU_K)/(R-1+delta_k))^alpha*N^(1-alpha)*(0.128/delta_kg)^eta)^(1/(1-alpha-eta));//用N和参数解出Y
W        =(1-alpha)*Y/N;//根据不同的N的取值，可以得到不同的W值
C        =0.417*Y;
I_G      =0.128*Y;
I        =delta_k*alpha*(1-TAU_K)/(R-1+delta_k)*Y;
G        =Y-C-I-I_G;
B        =0.129*Y;
K        =alpha*(1-TAU_K)/(R-1+delta_k)*Y;
K_G      =0.128/delta_kg*Y;
parameters rho_a rho_g rho_i_g phi_yig rho_tau_w rho_tau_c rho_tau_k rho_r phi_yg phi_wb rho_tw phi_kb rho_tk phi_cb rho_tc;
rho_a    =0.9;
rho_g    =0.7;
rho_i_g  =0.7;
rho_tau_w=0.5;
rho_tau_c=0.5;
rho_tau_k=0.5;
rho_r    =0.5;
phi_yg   =0.8;//可调
phi_yig  =0.8;//可调
phi_wb   =0.5;
rho_tw   =0.5;
phi_kb   =0.5;
rho_tk   =0.5;
phi_cb   =0.5;
rho_tc   =0.5;
model;  //模型
//一、居民部门
//1.居民跨期消费欧拉方程
//c(+1)-c+TAU_C/(1+TAU_C)*(tau_c(+1)-tau_c)=beta*(1-TAU_K)*RK*(-TAU_K/(1-TAU_K)*tau_k(+1)+rk(+1));   
c(+1)-c+TAU_C/(1+TAU_C)*(tau_c(+1)-tau_c)=r-pi(+1);
//2.居民用于投资的回报等于购买债券的回报，这个回报用于下一期消费的效用的折现值等于现期消费的效用
(1-TAU_K)*RK*(-TAU_K/(1-TAU_K)*tau_k(+1)+rk(+1))=R/PI*(r-pi(+1));//这里的rk是实际变量，已经除以P了
//3.居民选择最优劳动时，消费、劳动与工资的关系
c+TAU_C/(1+TAU_C)*tau_c-N/(1-N)*n=-TAU_W/(1-TAU_W)*tau_w+w;//w也是除以P以后的实际变量了
//二、厂商部门
//4.技术冲击
a=rho_a*a(-1)+epsilon_a;
//5.生产函数
y =a+alpha*k(-1)+(1-alpha)*n;
//6.资本形成方程
k=(1-delta_k)*k(-1)+delta_k*i;
//7.公共资本形成方程
k_g=(1-delta_kg)*k_g(-1)+delta_kg*i_g;
//8.实际资本价格
rk=a+(alpha-1)*k(-1)+(1-alpha)*n;
//9.实际劳动价格
w=a+alpha*k(-1)-alpha*n;
//10.实际边际成本
mc=-a+(1-alpha)*w+alpha*rk-eta*k_g;
//11.新凯恩斯菲利普斯曲线NKPC
pi=beta*pi(+1)+(1-beta*theta)*(1-theta)/theta*mc;
//三、ZF部门
//12.ZF预算约束
G*g+I_G*i_g+B*R*(b(-1)+r(-1))=TAU_C*C*(tau_c+c)+TAU_K*RK*K*(tau_k+rk+k(-1))+TAU_W*W*N*(tau_w+w+n)+B*b;
//13.ZF消费
g=rho_g*g(-1)+(1-rho_g)*phi_yg*y(-1)+epsilon_g;
//14.ZF投资
i_g=rho_i_g*i_g(-1)+(1-rho_i_g)*phi_yig*y(-1)+epsilon_i_g;
//15.16劳动税税率规则
tau_w=phi_wb*b(-1)+tw;
tw=rho_tw*tw(-1)-epsilon_tw;
//17.18资本税税率规则
tau_k=phi_kb*b(-1)+tk;
tk=rho_tk*tk(-1)-epsilon_tk;
//19.20消费税税率规则
tau_c=phi_cb*b(-1)+tc;
tc=rho_tc*tc(-1)-epsilon_tc;
//四、市场出清条件
//21.产品市场出清
Y*y=C*c+I*i+G*g+I_G*i_g;
//22.利率规则
//r=rho_r*r(-1)+(1-rho_r)*y+x;
//x=rho_x*x(-1)+epsilon_x;
end;
initval;
a      =0;
tau_c  =0;
tau_k  =0;
tau_w  =0;
c      =0;
n      =0;
w      =0;
rk     =0;
y      =0;
k      =0;
i      =0;
g      =0;
i_g    =0;
b      =0;
r      =0;
tw     =0;
tk     =0;
tc     =0;
end;
steady;
//resid(1);
check;
shocks;
//var epsilon_a;
//stderr 0.01;
var epsilon_g;
stderr 0.01;
//var epsilon_i_g;
//stderr 0.01;
//var epsilon_tw;
//stderr 0.01;
//var epsilon_tk;
//stderr 0.01;
//var epsilon_tc;
//stderr 0.01;
end;
stoch_simul(irf=50);
最后的利率规则是屏蔽的，加进去方程数就超过变量数了，请问大家，我原来的方程组里是不是有多余的，ZF预算约束需要进入方程组吗？

加入利率规则，删去ZF预算约束后，BK条件不满足了，有6个大于1的特征值，却只有五个跳跃变量。真的不想延期毕业啊，求指导拍砖。

likejia0015 发表于 2016-9-10 12:32
加入利率规则，删去ZF预算约束后，BK条件不满足了，有6个大于1的特征值，却只有五个跳跃变量。真的不想延期 ...

EIGENVALUES:
         Modulus             Real        Imaginary

      4.522e-016      -4.522e-016                0
      5.718e-016       5.718e-016                0
             0.5              0.5                0
             0.5              0.5                0
             0.5              0.5                0
             0.5              0.5                0
             0.7              0.7                0
             0.7              0.7                0
          0.8023           0.8023                0
             0.9              0.9                0
            0.95             0.95                0
            1.02             1.02                0
           2.239            2.085           0.8155
           2.239            2.085          -0.8155
      4.593e+016       4.593e+016                0
      1.734e+034      -1.734e+034                0
             Inf              Inf                0


There are 6 eigenvalue(s) larger than 1 in modulus
for 5 forward-looking variable(s)

The rank condition ISN'T verified!

我把五个反馈规则的系数都设置为0，报错QZ分解失败。现在我也不知道是模型的方程组写的有问题还是参数赋值有问题了。有没有大神帮帮忙？先行谢过了