求救，本人在编写DSGE模型的dynare文件，编完之后报错，具体文件内容如下：
var m pi Y M C I Rd Rs Rk Q MC W L K P picb Kcb phicb Ncb pisb Ksb phisb Nsb A S;




varexo ea em;




parameters gamma eta chi beta alpha delta tau thetaL thetaH sigma kappa rhom phipi phiy xip pB pG;




gamma=0.1986;
eta=3.4;
chi=0.25;
beta=0.99;
alpha=0.33;
delta=0.025;
tau=0.15;
thetaL=0.66;
thetaH=1.026;
sigma=0.95;
kappa=0.65;
rhom=0.7938;
phipi=3.65;
phiy=7.5861;
xip=0.75;
pB=0.66;
pG=1;




model;


(M/P)^(-1/gamma)+1/(C*Rd(+1))-1/C=0;
Rd=Rs;
W/(C*P)=chi*L^(eta);
beta*((C*P)/(C(+1)*P(+1)))*Rd(+1)=1;
pi=log(P/P(+1));
Y=A*(K^alpha)*(L^(1-alpha));
log(A)=log(A(-1))+ea;
K(-1)/L=(alpha/(1-alpha))*W/(Q*P);
L=(1-alpha)*(1/alpha)*K*((W/P)^(-1));
Rk(+1)*Q=P(+1)*alpha*(Y(+1)/K(+1))+(1-delta)*Q(+1);
P*(1-alpha)*(Y/L)=W;
picb=(Rk(+1)*((pG-pB)*thetaH+(pB-pG)*thetaL))/tau;
Q*Kcb=phicb*Ncb;
phicb=1/(1-(thetaL*Rk(+1))/(Rd(+1)));
Ncb=sigma*(Kcb/K(-1)-thetaL)*Rk*Q(-1)*Kcb(-1)+(1-sigma)*Ncb(-1);
Q*Ksb=Nsb+S;
phisb=Rs(+1)/(Rs(+1)-(Ksb(+1)/K*Rk(+1))+tau*(pisb^2));
Nsb=sigma*(Ksb/K(-1)*(Q+(1-delta)*Q)*Ksb(-1)-Rs*S(-1))+(1-sigma)*Nsb(-1);
m=log(M/M(-1));
m=rhom*m(-1)-(1-rhom)*(phipi*pi+phiy*Y)+em;
Y=C+I*(1+(I/I(-1)-1)*ln(I/I(-1)));
Q=1+(I/I(-1)-1)*ln(I/I(-1))+(ln(I/I(-1))+(I/I(-1)-1)*I/I(-1))*(I/I(-1)-(1/Rd(+1))*((I/I(-1))^2));
MC=(Q^alpha)*((W/P)^(1-alpha))*(alpha^(-alpha))*((1-alpha)^(alpha-1));
pi=beta/(1+beta)*pi(+1)+1/(1+beta)*pi(-1)-log(MC)*(1-xip)*(1-beta*xip)/((1+beta)*xip);
pisb=((Rs(+1)-tau*kappa-Rk(+1)*(pB*thetaH+(1-pB)*thetaL))+((Rs(+1)-tau*kappa-Rk(+1)*(pB*thetaH+(1-pB)*thetaL))^2+kappa*(Rk(+1)*(thetaH*(pG-pB)+thetaL*(pB-pG)))^2)^(1/2))/(Rk(+1)*((pG-pB)*thetaH+(pB-pG)*thetaL));
end;






initval;


m=1;
pi=0;
Y=0.7955;
M=2.1979;
C=0.5273;
I=0.2682;
Rd=1.0101;
Rs=1.0101;
Rk=1.0161;
Q=1;
MC=1;
W=0.388;
L=1.3737;
K=6.3872;
P=1;
picb=0.843;
Kcb=1.7267;
phicb=2.9755;
Ncb=0.5803;
pisb=0.7792;
Ksb=4.6605;
phisb=11.6672;
Nsb=0.3995;
A=1;
S=4.261;
end;






steady;


check;




shocks;


var ea;


stderr 0.048;


end;






shocks;


var em;


stderr 0.07;


end;






stoch_simul(order=1, periods=0,irf=20, conditional_variance_decomposition =[1:20]);




跑完之后显示结果：

错误使用 print_info (line 83)
Impossible to find the steady state. Either the model doesn't have a steady state, there are an infinity of steady states, or the guess
values are too far from the solution


出错 steady (line 104)
    print_info(info,options_.noprint, options_);


出错 lucky (line 302)
steady;


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


但是稳态求解好像是没有什么问题，求问大家这种情况怎么解决？