dynare 运行结果出问题了，求助

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

//封闭经济，无投资，无资本，无ZF支出，calvo定价，手工计算稳态，
var q Cj P i STj Wj MCh Yh yh ph PH U1 U2 BFkx;//q:债券价格，Cj:消费，P:价格指数，STj:劳动，Wj:工资，MCh：边际成本，Yh:总产出，yh:对商品的需求，ph:商品的价格，PH:本国商品价格指数（其实和P一回事，不过是为讨论开放经济做准备）,U1，U2:calvo定价中的辅助变量，BFkx:债券
varexo Mj;//货币供应量
parameters beta rou kai gamma niu you yeta n sigma theta alfa a lambdapai lambday kappa qbar Cjbar Mjbar Pbar ibar STjbar Wjbar MChbar Yhbar yhbar phbar PHbar U1bar U2bar BFkxbar;//很多参数在封闭经济下没有用，放在这里为开放经济讨论做准备
beta=0.99;
rou=3;//如果简化问题 rou=1
kai=8.2;
gamma=1;
niu=3;
you=0.1;
yeta=1;//封闭经济下
n=1;//封闭经济
sigma=10;
theta=1.9;
alfa=0.1;
a=0.3;
lambdapai=2;
lambday=2;
kappa=2;
qbar=beta;
yhbar=((sigma-1)/sigma)^(rou-niu);
Cjbar=yhbar;
Mjbar=3;
//以下手工计算的稳态
Wjbar=(Mjbar*((gamma*((sigma-1)/sigma)^(rou^2-rou*niu))/(1-q))^(-1/kai))*(sigma-1)/sigma;
phbar=(sigma/(sigma-1))*Wjbar;
PHbar=phbar;
Pbar=phbar;
ibar=1/qbar-1;
STjbar=yhbar;
MChbar=Wjbar;
Yhbar=yhbar;
U1bar=Wjbar*yhbar/(1-qbar*alfa);
U2bar=yhbar/(1-qbar*alfa);
BFkxbar=0;

model(linear);
//对数线性近似后的模型，14个内生变量，14个方程
q-rou*Cj=P-rou*Cj(+1)-P(+1); //本国欧拉函数
kai*Mj-kai*P=rou*Cj-i/(1+ibar);//货币需求
STj=Wj-P-rou*Cj;//劳动供给
yh=-sigma*(ph-PH)+Cj;//消费篮子分配于h商品产生的需求
Yh=yh; //封闭经济下本国出清
Yh=STj;//只有劳动这个生产要素
MCh=Wj;//只有劳动这个要素
//以下calvo定价
ph=U1-U2;
U1=(MChbar*yhbar*(yh+MCh)+alfa*qbar*U1bar*(U1(+1)+q))/(MChbar*yhbar+qbar*alfa*U1bar);
U2=(yhbar*yh+qbar*alfa*U2bar*(q+U2(+1)))/(yhbar+qbar*alfa*U2bar);
//以下价格指数
PH=((1-alfa)*phbar^(1-sigma)*ph+alfa*PHbar^(1-sigma)*PH(-1))/((1-alfa)*phbar^(1-sigma)+alfa*PHbar^(1-sigma));
P=PH;
//以下资源约束
phbar*yhbar*yh+alfa*phbar*yhbar*ph(-1)+(1-alfa)*yhbar*phbar*ph=(Pbar*Cjbar*(P+Cj)-BFkxbar*qbar*(BFkx+q)+BFkxbar*BFkx(-1))/(Pbar*Cjbar-BFkxbar*qbar+BFkxbar);

q=-ibar/(1+ibar)*i;//债券价格与名义利率关系
end;

initval;
q=qbar;
Mj=Mjbar;
Cj=Cjbar;
P=Pbar;
i=ibar;
STj=STjbar;
Wj=Wjbar;
MCh=MChbar;
Yh=Yhbar;
yh=yhbar;
ph=phbar;
PH=PHbar;
U1=U1bar;
U2=U2bar;
BFkx=BFkxbar;


end;

resid(1);
steady;
solve algo = 1;
check;
shocks;
var Mj = 2;//什么意思
end;
stoch_simul(periods=2100);



运行的结果：

Residuals of the static equations:
Equation number 1 : 0.99Equation number 2 : -2.99Equation number 3 : 4.3Equation number 4 : 0Equation number 5 : 0Equation number 6 : 0Equation number 7 : 0Equation number 8 : 1.1132Equation number 9 : -0.73171Equation number 10 : 0.00099Equation number 11 : 0Equation number 12 : 0Equation number 13 : 8Equation number 14 : 0.9901

Warning: Matrix is singular to working precision.> In evaluate_steady_state at 85  In steady_ at 54  In steady at 81  In question at 294  In dynare at 180



Residuals of the static equations:
Equation number 1 : NaNEquation number 2 : NaNEquation number 3 : NaNEquation number 4 : NaNEquation number 5 : NaNEquation number 6 : NaNEquation number 7 : NaNEquation number 8 : NaNEquation number 9 : NaNEquation number 10 : NaNEquation number 11 : NaNEquation number 12 : NaNEquation number 13 : NaNEquation number 14 : NaN

??? Error using ==> print_infoThe steady state contains NaN or Inf
Error in ==> steady at 92    print_info(info,options_.noprint, options_);
Error in ==> question at 294steady;
Error in ==> dynare at 180evalin('base',fname) ;

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏3 回帖

关键词：dynare ARE parameters Equations Residuals sigma kappa 成本 手工

出现奇异矩阵了，在网上搜也没有搜到什么结果，是我对理论的理解有问题，模型建错了，还是什么问题，请求帮助，比较急，谢谢！

我发的贴有问题吗，似乎大家都看不到？

Equation number 1 : NaNEquation number 2 : NaNEquation number 3 : NaNEquation number 4 : NaNEquation number 5 : NaNEquation number 6 : NaNEquation number 7 : NaNEquation number 8 : NaNEquation number 9 : NaNEquation number 10 : NaNEquation number 11 : NaNEquation number 12 : NaNEquation number 13 : NaNEquation number 14 : NaN

说明你的参数设定可能有问题。

请问您的问题解决了吗？我也碰到了同样的问题，能告诉我该怎么解决吗？谢谢

榄外人L 发表于 2014-9-26 21:44
请问您的问题解决了吗？我也碰到了同样的问题，能告诉我该怎么解决吗？谢谢

说明你的稳态求解有问题。请将稳态带入方程确保每一个方程恒等于0

ahnulxy 发表于 2014-9-26 22:50
说明你的稳态求解有问题。请将稳态带入方程确保每一个方程恒等于0

对数线性化之后，哪些稳态不是都可以设为零的么？还需要再去重新求稳态吗？那我要是求出稳态之后该怎么用呢？我是菜鸟，正在试着学习，请您多多指导

Dynare forum,去看看吧，不过是英文的。这种情况下，很大可能是计算错误（稳态值），还有可能是timing写错了），如果不是计算错误什么的，那就是整个模型的结构错了。

榄外人L 发表于 2014-9-27 15:34
对数线性化之后，哪些稳态不是都可以设为零的么？还需要再去重新求稳态吗？那我要是求出稳态之后该怎么用 ...

您需要做的是先求出稳态，然而直接书写FOC条件，不需要线性化，最后在初始赋值模块里，将稳态值赋予给变量即可。自己线性化，容易出错，而且工作量极大，没有必要，让dynare自己帮你完成这个工作吧

榄外人L 发表于 2014-9-27 15:34
对数线性化之后，哪些稳态不是都可以设为零的么？还需要再去重新求稳态吗？那我要是求出稳态之后该怎么用 ...

在最后的stoch_simul 命令中 写order=1 就是线性化了