你好，我也遇到了这种情况，我的内生变量是31个，但是方程是32个，要怎么处理才好呢


// 内生变量 31
var c k d i ni nj wi wj rd rk ki kj gs yi a kg gi ri li yj rj lj f l t e g y n lambda mu;
// 外生变量
varexo eps_a eps_ig eps_sg eps_f eps_e ;
// 参数
parameters beta thet1 thet2 alpha1 alpha2 delta_k delta_g tau s rho_a rho_ig rho_sg rho_f rho_e
A LI_L LJ_L RK RD MU RI RJ F YI_Y YJ_Y GS_G GI_G C_Y LI_YI LJ_YJ LI_LJ YI_YJ YJ_YI KI_KJ KI_K KJ_K KI_YI KJ_YJ WINI_YI WJNJ_YJ K_Y I_Y G_Y T_Y LI_Y E_Y NI_NJ NI_N NJ_N;
beta=0.99;thet1=0.3;thet2=0.2;alpha1=0.45;alpha2=0.2;delta_k=0.04;delta_g=0.04;tau=0.09;s=0.15;
rho_a=0.75;rho_ig=0.35;rho_sg=0.35;rho_f=0.9;rho_e=0.9;
A=1;
LI_L=0.09;
LJ_L=1-LI_L;
RK=1/beta+delta_k-1;
RD=1/beta;
MU=(s-RD)/(s-1);
RI=MU/1.5;
RJ=MU;
F=0.5*RI;
GS_G=0.2;
GI_G=0.8;
C_Y=0.384;

LI_YI=(1-tau)/RI;
LJ_YJ=(1-tau)/RJ;
LI_LJ=LI_L/(1-LI_L);
YI_YJ=LI_LJ*(RI/RJ);
YJ_YI=1/YI_YJ;
YI_Y=YI_YJ/(YI_YJ+1);
YJ_Y=1-YI_Y;

KI_KJ=YI_YJ*(RJ/RI);
KI_K=KI_KJ/(KI_KJ+1);
KJ_K=1-KI_K;
KI_YI=(1-tau)*alpha1/(RI*(RK+delta_k));
KJ_YJ=(1-tau)*alpha1/(RJ*(RK+delta_k));

WINI_YI=(1-tau)*(1-alpha)/RI;
WJNJ_YJ=(1-tau)*(1-alpha)/RJ;
K_Y=KJ_YJ*((1+KI_KJ)/(1+YI_YJ));
I_Y=delta_k*K_Y;
G_Y=1-C_Y-I_Y;
T_Y=tau;
LI_Y=LI_YI/(1+YJ_YI);
E_Y=G_Y+F*LI_Y-T_Y;

NI_NJ=((YI_YJ)*(KI_KJ^alpha1))^(1/(1-alpha1));
NI_N=NI_NJ/(NI_NJ+1);
NJ_N=1-NI_N;

// 27+5 Equations of the model
predetermined_variables k ;
model(linear);
//家庭
k(+1)=delta_k*i+(1-delta_k)*k;
lambda+c=0;
wi=c+thet1*ni;
wj=c+thet2*nj;
c(+1)=c+beta*(1/beta-1+delta_k)*rk(+1);
c(+1)=c+rd;
//企业
kg(+1)=delta_g*gi+(1-delta_g)*kg;
yi=a+alpha1*ki+alpha2*kg+(1-alpha1)*ni;
LI_YI*li=WINI_YI*(wi+ni)+RK*KI_YI*(rk+ki)+delta_k*KI_YI*ki;
yi=ki+ni+wi;
yi+ki+ri+rk*RK/(RK+delta_k);
yj=a+alpha1*kj+alpha2*kg+(1-alpha1)*nj;
LJ_YJ*lj=WJNJ_YJ*(wj+nj)+RK*KJ_YJ*(rk+kj)+delta_k*KJ_YJ*kj;
yj=kj+nj+wj;
yj+kj+rj+rk*RK/(RK+delta_k);
y=YI_Y*yi+YJ_Y*yj;
l=LI_L*li+LJ_L*lj;
//银行
l=d;
MU*mu=RI*ri+F*f;
mu=rj;
mu=rd;
//政府
T_Y*t+E_Y*e=G_Y*g+F*LI_Y*(f+li);
t=y;
g=GS_G*gs+GI_G*gi;
//市场出清
y=C_Y*c+I_Y*i+G_Y*g;
k=KI_K*ki+KJ_K*kj;
n=NI_N*ni+NJ_N*nj;
//冲击
a=rho_a*a(-1)+eps_a;
gi=rho_ig*gi(-1)+eps_ig;
gs=rho_sg*gs(-1)+eps_sg;
f=rho_f*f(-1)+eps_f;
e=rho_e*e(-1)+eps_e;

end;