DSGE模型求解出现问题

本人在求解DSGE模型的时候遇到了困难，现已将模型设置好，但是Dynare代码出现错误如下，恳请前辈们慷慨相助。。。或者联系本人QQ350641328或791935570，定有重谢！！！本人的程序如下所示：
var a b y c  k kG I IG l lA lB m n i w wA wB r;
varexo epsm epsn epsi;
//-----------------------------------------------
// 2. Declaration of Parameters
//-----------------------------------------------
parameters
tauK
tauL
G
chi
beta
alphaA
betaA
alphaB
betaB
gamma
theta
delta
deltag
rhom
rhon
rhog
rhoi
sigmam
sigman
sigmai;
//-----------------------------------------------
// 3. Parameter Calibration
//-----------------------------------------------
tauK=0.002;
tauL=0.002;
G = 0.1;
beta = 0.9992;
alphaA = 0.21;
betaA = 0.79;
alphaB =0.31;
betaB =0.42;
gamma =0.27;
theta=0.5;
delta = 0.02;
deltag=0.05;
rhom = 0.9;
rhon =0.9;
rhog = 0.9;
rhoi=0.9;
chi = 1;
sigmam = 0.001;
sigman = 0.001;
sigmai=0.001;
// parameters used for initial value calculation
//-----------------------------------------------
// 4. The MODEL
//-----------------------------------------------
model;
(1/(c)) = beta*(1/(c(+1)))*(1+(1-tauK)*r-delta);
c = (1-tauL)*w*l +(1-tauK)*r*k-I;
chi*(c)/(1-l) = (1-tauL)*w;
a = exp(m)*(kG^alphaA)*(lA^betaA);
b = exp(n)*(k^alphaB)*(lB^betaB);
y = (theta*b^gamma+(1-theta)*a^gamma)^(1/gamma);
y=c+I+IG+G;
l =lA+lB;
I = k(+1)-(1-delta)*k;
IG = k(+1)-(1-deltag)*k;
wA*lA = betaA*a;
wB*lB = betaB*b;
w*l= wA*lA+wB*lB;
r*k= alphaB*b;
G + exp(i)*IG = w*tauL*l + r*tauK*k;
exp(m) = exp((m(-1))^rhom)*exp(epsm);
exp(i) = exp((i(-1))^rhoi)*exp(epsi);
exp(n) = exp((n(-1))^rhon)*exp(epsn);
end;
resid(1);
//-----------------------------------------------
// 5. Initial guesses for steady-state computation
//-----------------------------------------------
initval;
k=20;
kG=8;
I=0.3;
IG=0.1;
l=0.8;
lA=0.2;
lB=0.6;
m=0.001;
n=0.001;
i=0.001;
a=exp(m)*(kG^alphaA)*(lA^betaA);
b=exp(n)*(k^alphaB)*(lB^betaB);
y=(theta*b^gamma+(1-theta)*a^gamma)^(1/gamma);
c=(1-tauL)*(betaA*exp(m)*(kG^alphaA)*(lA^betaA)+betaB*exp(n)*(k^alphaB)*(lB^betaB))*(1-l)/l;
w=betaA*exp(m)*(kG^alphaA)*(lA^betaA)+betaB*exp(n)*(k^alphaB)*(lB^betaB);
wA=betaA*exp(m)*(kG^alphaA)*(lA^betaA)/lA;
wB=betaB*exp(n)*(k^alphaB)*(lB^betaB)/lB;
r=alphaB*exp(n)*(k^alphaB)*(lB^betaB)/k;
end;
//-----------------------------------------------
// 6. Specification of shocks
//-----------------------------------------------
shocks;
var epsm = sigmam^2;
var epsn = sigman^2;
var epsi = sigmai^2;
end;
steady;
solve_algo=1;

stoch_simul ( order=1);

本人的错误如下所示：??? Error using ==> print_info at 74
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

Error in ==> steady at 92
    print_info(info,options_.noprint, options_);

Error in ==> try3 at 277
steady;

Error in ==> dynare at 162
evalin('base',fname) ;
  

看报错应该是你的稳态值计算错误了

initval;是不是应该放到model;的前面啊。另外，代码中为什么不将大写字母换成小写呢

大力鼓励广大网友发送学术及技术问题来宏观经济学板块探讨。

richardgu26 发表于 2016-2-14 16:16
看报错应该是你的稳态值计算错误了

前辈您好，请问您怎么求稳态值呢？若您方便可否加QQ聊，我的是791935570

孤独的散步者翱 发表于 2016-2-14 16:20
initval;是不是应该放到model;的前面啊。另外，代码中为什么不将大写字母换成小写呢

为了方便识别，所以用大写了，感谢您回复，能否再费神诊断一下？不甚感激

791935570 发表于 2016-2-14 16:59
前辈您好，请问您怎么求稳态值呢？若您方便可否加QQ聊，我的是791935570

你把所有的一阶条件以及出清条件都放在一起，然后算他们的稳态值，应该会出现一组方程。其中，一些核心变量，就是你在编成一开始出现的那些，可以用结构参数来表示的，而剩下的一部分可以用其他核心变量和结构参数来表示。这样，你就可以求的初始值了。

还有，你的model里，有一些是exp(n)什么的。我不知道你的模型具体怎样，感觉应该是非线性化的。你再检查下你的模型推导吧。

richardgu26 发表于 2016-2-14 19:33
你把所有的一阶条件以及出清条件都放在一起，然后算他们的稳态值，应该会出现一组方程。其中，一些核心变 ...

前辈您好，请问这个核心变量怎样用参数表示，仅仅依靠原始方程无法表示呀，总会有变量