解决赠送50个论坛币！

稳态含有复数或者无穷值，mod文件发上来看看

var
jch
MH
acloeb
acloee
cb
ce
ch
d
he
hh
loe
nh
q
re
rh
wh
y  ;

varexo jh,mh;

parameters BETAH BETAE BETAB TAUH NU RHO_J MN RHO_M GAMMAE FILOEE FILOEB LOE_SS;
BETAH   = 0.9925;   
BETAE   = 0.94;   
BETAB   = 0.945;   
TAUH    = 2;
NU      = 0.05;
RHO_J   = 0.1;
MN      = 1;
RHO_M   = 0.9;
GAMMAE  = 0.9;
FILOEE  = 0.25;
FILOEB  = 0.25;
LOE_SS  = 0.483;

phi = 0.1;

//Equations

model;

//% Households
ch+d+q*(hh-hh(-1))=rh(-1)*d(-1)+wh*nh;

1/ch=BETAH*rh*(1/(ch(+1)));

wh/ch=(TAUH/(1-nh));

q/ch=jch/hh+BETAH*q(+1)/(ch(+1));



//% Banks
cb+rh(-1)*d(-1)+loe=d+re*loe(-1);

d=GAMMAE*(loe);

(1-GAMMAE+acloeb)*(cb(+1))=BETAB*(re(+1)-GAMMAE*rh)*(cb);



//% Entrepreneurs
ce+q*(he-he(-1))+re*(loe(-1))+wh*nh=y;

y=((he(-1))^NU)*(nh^(1-NU));

loe=MH*(q(+1)/(re(+1)))*he-MN*wh*nh;

(q-(1-acloee)*((MH*q(+1))/(re(+1))))*(1/ce)=
BETAE*(q(+1)*(1-MH)+NU*(y(+1)/he))*(1/(ce(+1)));


(1-NU)*y=wh*nh*(1+MN*(1-acloee-BETAE*re(+1)*ce/ce(+1)));

he+hh=1;

jch=0.075+RHO_J*jch(-1)+jh;
MH=0.9+RHO_M*MH(-1)+mh;

acloeb = FILOEB*((loe)-(loe(-1)))/LOE_SS ;
acloee = FILOEE*((loe)-(loe(-1)))/LOE_SS;
end;

initval;
ch      = 0.5;
d      = 0.5;
wh      = 1;
nh      = 0.3436;
rh     = 1.0075;
re     = 1.0104;
he     = 0.2639;
hh     = 0.8;
q      = 3;
loe      = 0.5;
y      = 0.5;
ce     = 0.01;
cb     = 0.01;
jch    = 0;
MH     = 0;
acloee=0;
acloeb=0;
end;

steady;

check;

shocks;
var jh; stderr 0.009;
var mh; stderr 0.009;
var jh, mh = phi*0.009*0.009;
end;

stoch_simul;