下面是我的dynare代码  最近刚学习  不太知道哪里错误
var c l d b c_1 l_1 d_1 b_1 y d_z b_z r w a b_z q c_z x pai;
varexo etaa;
parameters beta_a beta_b beta_z varphi m_b m_z rhoa alpha omega1 omega2 omega3 omega4 D1 L1 C1 B1 D L C B Y DZ BZ CZ W sigmaa X Q R theta PAI;
beta_a = 0.979;
beta_b = 0.97;
beta_z = 0.98;
varphi = 0.01;
m_b =  0.75;
m_z = 0.5;
rhoa = 0.96;
alpha = 0.4;
theta = 0.75;
sigmaa = 0.01;
PAI=1;
R=1/beta_a;
X = 1.2;
omega1 = (1-beta_a);
omega2 = (1-(beta_a-beta_b)*m_b-beta_b);
omega3 = X/(1-alpha)*(1-(beta_a-beta_z)*m_z-beta_z);
omega5 = (omega2+(1-beta_a)*m_b)*X/(1-0.7)/alpha;
omega6 = omega2/omega5;
omega7 = ((1-alpha)/X)-(1-beta_a)*m_z/omega3;
omega4 = (0.7*alpha/X)^(0.7)*((0.3)*alpha/X)^(0.3)*(1/omega3)^((1-alpha)/alpha);
Y=1;
C1 = omega6*Y;
CZ = omega7*Y;
C = (1-omega6-omega7)*Y;
Q = (1/omega3+omega6/omega2+(1-omega6-omega7)/omega1)*Y;
D1 = omega6/omega2/Q;
DZ = 1/omega3/Q;
D = (1-omega6-omega7)/omega1/Q;
B1 = (m_b*beta_a*omega6/omega2)*Y;
BZ = (m_z*beta_a/omega3)*Y;
B = B1+BZ;
W = omega4*(Q)^(-(1/alpha));
L = (W/C)^(1/varphi);
L1 = (W/C1)^(1/varphi);
model;
(l)^(varphi)=w/c;
(1/d)-(q/c)+beta_a*q(+1)/c(+1)=0;
c(+1)/c=beta_a*r/pai(+1);
c_1+r(-1)*b_1(-1)/pai+q*(d_1-d_1(-1))=w*l_1+b_1;
b_1*r=m_b*q(+1)*d_1*pai(+1);
(l_1)^(varphi)=w/c_1;
(1/d_1)-(q/c_1)+(1/c_1/r)*m_b*q(+1)*pai(+1)-beta_b/c_1(+1)*m_b*q(+1)+beta_b*q(+1)/c_1(+1)=0;
y=a*(l)^(0.28)*(l_1)^(0.12)*(d_z)^(0.6);
c_z+w*(l+l_1)+q*(d_z-d_z(-1))+r(-1)*b_z(-1)/pai=y/x+b_z;
b_z*r = m_z*q(+1)*d_z*pai(+1);
(((1-alpha)*y/d_z/x)-q)/c_z+(1/r/c_z-beta_z/c_z(+1)/pai(+1))*m_z*q(+1)*pai(+1)+beta_z*q(+1)/c_z(+1)=0;
0.7*alpha*y/x/l=w;
(1-0.7)*alpha*y/x/l_1=w;
pai=beta_a*pai(+1)-(1-theta)*(1-beta_a*theta)*x/theta;
y=c+c+1+c_z;
1=d+d_1+d_z;
b=b_1+b_z;
log(a)=rhoa*log(a(-1))-etaa;
end;
steady_state_model;
c=C;
d=D;
l=L;
b=B;
c_1=C1;
d_1=D1;
l_1=L1;
b_1=B1;
y=Y;
A=1;
r=R;
c_z=CZ;
d_z=DZ;
b_z=BZ;
q=Q;
x=X;
w=W;
pai=PAI;
end;
shocks;
var etaa = sigmaa^2;
end;
stoch_simul(irf = 30, order = 1, periods = 0);