var y pi m r i oh c q of h w rk k u x a z p gama_w gama_p v;
varexo e_z e_x e_a e_p e_gama_w e_gama_p e_v;
parameters sigma_m sigma_c kesi beta sigma fy xi_w sigma_L alpha_1 alpha_2
alpha_3 delta xi_p omiga_r omiga_pi omiga_y omiga_u ep_c ep_i ep_o ro_a ro_x
ro_z ro_p ro_gama_p ro_gama_w ro_v sigma_a sigma_x sigma_gama_w sigma_p sigma_z
sigma_gama_p sigma_v ;

sigma_m=1.4;
sigma_c=2;
kesi=0.98;
beta=0.992;
sigma=0.03;
fy=4;
xi_w=0.5;
sigma_L=1;
alpha_1=0.6;
alpha_2=0.21;
alpha_3=0.19;
delta=0.5;
xi_p=0.5;
omiga_u=1;
omiga_pi=2;
omiga_y=2;
omiga_u=1;
ep_c=0.3649;
ep_i=0.4729;
ep_m=0.1622;
ro_a=0.5;
ro_x=0.5;
ro_p=0.5;
ro_gama_w=0.5;
ro_gama_p=0.5;
sigma_a=0.1;
sigma_x=0.1;
sigma_gama_w=0.1;
sigma_p=0.1;
sigma_z=0.1;
sigma_gama_p=0.1;
sigma_y=0.1;

model(linear);
sigma_m*m=(sigma_c-sigma_c*kesi)*oh+(sigma_c*kesi-kesi+1)*c-beta/(1-beta)*r+beta/(1-beta)*pi(+1);
p=(1-kesi)*c-kesi*oh(-1);
r=pi(+1)+a-a(+1)+(sigma_c*kesi-sigma_c)*(oh-oh(-1))+(kesi-sigma_c*kesi-1)*(c-c(-1));
q=beta*(1-sigma)*q(+1)-r+pi(+1)+(1-beta*(1-sigma))*rk(+1);
i=beta/(1+beta)*i(+1)+1/(1+beta)*i+fy/(1+beta)*q+1/(1+beta)*(beta*x(+1)-x);
k=(1-sigma)*k(-1)+sigma*i;
w=k(-1)+rk-h;
rk=of+p-k(-1);
w=1/(1+beta)*(pi(-1)+w(-1))-pi+beta/(1+beta)*(pi(+1)+w(+1))+(1-beta*xi_w)*(1-xi_w)*(-w+sigma_L*h+(sigma_c-sigma_c*kesi)*oh+(sigma_c*kesi-kesi+1)*c+gama_w);
y=z+alpha_1*h+alpha_2*k(-1)+alpha_3*of;
pi=(1-delta)*pi(-1)+beta*(pi(+1)-(1-delta)*pi)+(1-beta*xi_p)*(1-xi_p)/xi_p*(-z+alpha_1*w+alpha_2*rk+alpha_3*p+gama_p);
u=m-m(-1)+pi;
r=omiga_r*r(-1)+(1-omiga_r)*(omiga_pi*pi+omiga_y*y+omiga_u*u)+v;
y=ep_c*c+ep_i*i+ep_o*(oh+of+p);
a=sigma_a*e_a+ro_a*a(-1);
x=sigma_x*e_x+ro_x*x(-1);
z=sigma_z*e_z+ro_z*z(-1);
p=sigma_p*e_p+ro_p*p(-1);
gama_p=sigma_gama_p*e_gama_p+ro_gama_p*gama_p(-1);
gama_w=sigma_gama_w*e_gama_w+ro_gama_w*gama_w(-1);
v=sigma_v*e_v+ro_v*v(-1);
end;

initval;
y=0.0;
pi=0.0;
m=0.0;
r=0.0;
i=0.0;
oh=0.0;
c=0.0;
q=0.0;
of=0.0;
h=0.0;
w=0.0;
rk=0.0;
k=0.0;
u=0.0;
x=0.0;
a=0.0;
z=0.0;
p=0.0;
v=0.0;
gama_w=0.0;
gama_p=0.0;
end;

steady;
check;

shocks;
var e_p;
stderr 1;
var e_z;
stderr 1;
var e_a;
stderr 1;
var e_x;
stderr 1;
var e_gama_p;
stderr 1;
Var e_gama_w;
stderr 1;
var e_v;
stderr 1;
end;

estimated_params;
sigma_c,normal_pdf,2,0.25;
kesi,normal_pdf,0.98,0.01;
fy,normal_pdf,4,1.5;
xi_p,beta_pdf,0.5,0.1;
xi_w,beta_pdf,0.5,0.1;
delta,beta_pdf,0.5,0.015;
omiga_pi,normal_pdf,-0.5,0.1;
omiga_y,normal_pdf,-0.5,0.1;
omiga_u,normal_pdf,-0.5,0.1;
ro_a,beta_pdf,0.5,0.2;
ro_x,beta_pdf,0.5,0.2;
ro_z,beta_pdf,0.5,0.2;
ro_p,beta_pdf,0.5,0.2;
ro_gama_p,beta_pdf,0.5,0.2;
ro_gama_w,beta_pdf,0.5,0.2;
ro_v,beta_pdf,0.5,0.2;
sigma_a,inv_gamma_pdf,0.1,2;
sigma_x,inv_gamma_pdf,0.1,2;
sigma_z,inv_gamma_pdf,0.1,2;
sigma_p,inv_gamma_pdf,0.1,2;
sigma_v,inv_gamma_pdf,0.1,2;
sigma_gama_p,inv_gamma_pdf,0.1,2;
sigma_gama_w,inv_gamma_pdf,0.1,2;
end;
varobs y pi m r i oh c q of h w ;
estimation(datafile=datass,mh_replic=10000);
stoch_simul(periods=2100);

2楼附有代码还有mod文件，跪求各位大神指教啊

找到高人了么，我的mod也跑不出来

你们dsge都用什么软件跑的啊

啊坟贴勿回