我按照GK2010的论文用dynare模拟结果，运行之后显示“The rank condition ISN'T verified!”以及“Blanchard Kahn conditions are not satisfied: indeterminacy due to rank failure”，请各位大佬指点！
代码如下：
var y k c z l uc r qi i phin ni nn n d qbar phii nu omegai omegan mui qn rkii rkin rkni rknn wuc wc wi womegai womegan wrkii wrkin wrkni wrknn a psi phic psig;
varexo ea epsi;
parameters alpha gamma beta pii pin chi epsilon theta xi rhoa rhopsi sigmaa sigmapsi delta sigma kappa C Y Z I Ni Nn K Qn RKin RKii N D Qbar MUi PHIi PHIn OMEGAi OMEGAn NU RKni RKnn etac etauc etai etaomegai etaomegan etarkii etarkin etarkni etarknn nug;
alpha = 0.33;
gamma = 0.5;
beta = 0.99;
pii = 0.25;
pin = 1-pii;
chi = 5.013;
epsilon = 0.1;
theta = 0.1288;
xi = 0.0019;
rhoa=0.95;
sigmaa=0.02;
rhopsi=0.66;
sigmapsi=-0.05;
delta = 0.025;
sigma = 0.972;
kappa = 1;
C = 3.2706;
Y = 4.1817;
Z = 0.0379;
I = 0.9110;
Ni = 2.2781;
Nn = 6.8323;
K = 36.4417;
Qn = 1.0094;
RKin = 1.0220;
RKii = 1.0129;
N = 9.1104;
D = 27.3313;
Qbar = 1.0070;
MUi = 0.0050;
PHIi = 4.2991;
PHIn = 3.9369;
OMEGAi = 0.5661;
OMEGAn = 0.5207;
NU = 0.5321;
RKni = 1.0034;
RKnn = 1.0125;
etac = 0;
etauc = 0;
etai = 0;
etaomegai = 0;
etaomegan = 0;
etarkii = 0;
etarkin = 0;
etarkni = 0;
etarknn = 0;
nug =0;

model(linear);
y = a + alpha*k +(1-alpha)*l;  // prod f
k(+1) - psi(+1) = (1-delta)*k+delta*i;  // k accum
y -C/Y*c -I/Y*i=0;  // mkt clearing
z -a - (1-alpha)*l+(1-alpha)*k=0;  // k gross return
l - 1/(epsilon + alpha)*a - alpha/(epsilon + alpha)*k - 1/(epsilon + alpha)*uc=0;  // labor
(1-beta*gamma)*(1-gamma)*uc + (1+beta*gamma^2)*c - (beta*gamma)*wc = gamma*c(-1);  // mg U of C
r - uc + wuc = 0;  // risk-free rate
qi -(1-beta)*kappa*i - beta*kappa*wi = -kappa*i(-1); // qi
qi + delta/(delta+pii*(1-delta))*i + pii*(1-delta)/(delta+pii*(1-delta))*k - phic- ni=0;  // investment demand
qn + k -phin - nn =0; // leverage of n-banks
ni - (sigma+xi)*pii*K*Z/Ni*z - (sigma+xi)*pii*K*(1-delta)/Ni*qi -(sigma+xi)*pii*K*RKii/Ni*k(-1) = -sigma*pii*D/(beta*Ni)*r - sigma*pii*D/(beta*Ni)*d(-1); // net worth of i-banks
nn -(sigma+xi)*pin*K*Z/Nn*z -(sigma+xi)*pin*K*(1-delta)*Qn/Nn*qn -(sigma+xi)*pin*K*RKin/Nn*k(-1) = -sigma*pin*D/(beta*Nn)*r - sigma*pin*D/(beta*Nn)*d(-1);  // net worth of n-banks
n -Ni/N*ni-Nn/N*nn = 0;  // Aggregate net worth
D*d -Ni*PHIi*phic - (PHIi-1)*Ni*ni - PHIn*Nn*phin - (PHIn-1)*Nn*nn = 0;   // debt
qbar -pii/Qbar*qi - pin*Qn/Qbar*qn = 0;  // qbar
phii - nu - MUi/(theta-MUi)*mui = 0;  // phii (leverage of i-banks)
(pii*OMEGAi+pin*OMEGAn)*nu - pii*OMEGAi*womegai - pin*OMEGAn*womegan = 0;  // nu
omegai -(sigma*theta*PHIi)/(1-sigma+sigma*theta*PHIi)*phii = 0;   // omega_i
omegan -sigma*NU/(1-sigma+sigma*NU)*nu = 0;  // omega_n
MUi*mui -pii*OMEGAi*beta*RKii*wrkii-pin*OMEGAn*beta*RKin +beta*(pii*OMEGAi*RKii+pin*OMEGAn*RKin)*r -pii*OMEGAi*(beta*RKii-1)*womegai -pin*OMEGAn*(beta*RKin-1)*womegan =0;  // mui
pii*OMEGAi*beta*RKni*wrkni +pin*OMEGAn*beta*RKnn*wrknn -(pii*OMEGAi*beta*RKni+pin*OMEGAn*beta*RKnn)*r + pii*OMEGAi*(beta*RKni-1)*womegai +pin*OMEGAn*(beta*RKnn-1)*womegan = 0;  // mun=0
rkii(+1) - psi(+1) -Z/(Z+1-delta)*z(+1) - (1-delta)/(Z+1-delta)*qi(+1) = -qi;   // rkii
rkin(+1) - psi(+1) -Z/(Z+(1-delta)*Qn)*z(+1) -(1-delta)*Qn/(Z+(1-delta)*Qn)*qn(+1) = -qi;  // rkin
rkni(+1) - psi(+1) -Z/(Z+1-delta)*z(+1) -(1-delta)/(Z+1-delta)*qi(+1) = -qn;  // rkni
rknn(+1) - psi(+1) -Z/(Z+(1-delta)*Qn)*z(+1) -(1-delta)/(Z+(1-delta)*Qn)*qn(+1) = -qn;  // rknn
c = wc(-1) + etac;  // wc
uc = wuc(-1) +etauc;  // wuc
i = wi(-1) +etai;  // wi
omegai = womegai(-1) +etaomegai;  // womegai
omegan = womegan(-1) +etaomegan;  // womegan
rkii(+1) = wrkii +etarkii(+1);  // wrkii
rkin(+1) = wrkin +etarkin(+1);  // wrkin
rkni(+1) = wrkni +etarkni(+1);  // wrkni
rknn(+1) = wrknn +etarknn(+1);  // wrknn
a = rhoa*a(-1) +sigmaa*ea;  //a
psi = rhopsi*psi(-1)+sigmapsi*epsi;  //psi
phic - phii -psig = 0;  // effective leverage of investing banks
psig -nug*pii*RKii*wrkii -nug*pin*RKin*wrkin + nug/beta*r = 0;  // policy response to spread
end;
initval;
y = 0;
k = 0;
c = 0;
z = 0;
l = 0;
uc = 0;
r = 0;
qi = 0;
i = 0;
phin = 0;
ni = 0;
nn = 0;
n = 0;
d = 0;
qbar = 0;
phii = 0;
nu = 0;
omegai = 0;
omegan = 0;
mui = 0;
qn = 0;
rkii = 0;
rkin = 0;
rkni = 0;
rknn = 0;
wuc = 0;
wc = 0;
wi = 0;
womegai = 0;
womegan = 0;
wrkii = 0;
wrkin = 0;
wrkni = 0;
wrknn = 0;
a = 0;
psi = 0;
phic = 0;
psig = 0;
end;
steady;
check;
shocks;
var ea = sigmaa^2;
var epsi = sigmapsi^2;
end;

stoch_simul;