匿名用户不能发表回复！|
每天回帖即可获得10分可用分！小技巧：
你还可以输入10000个字符
(Ctrl+Enter)
请遵守CSDN，不得违反国家法律法规。
转载文章请注明出自“CSDN（www.csdn.net）”。如是商业用途请联系原作者。苹果/安卓/wp
积分 57, 距离下一级还需 28 积分
道具: 彩虹炫, 涂鸦板, 雷达卡, 热点灯, 金钱卡, 显身卡下一级可获得
权限: 自定义头衔
购买后可立即获得
权限: 隐身
道具: 金钱卡, 彩虹炫, 雷达卡, 热点灯, 涂鸦板
开心签到天数: 12 天连续签到: 1 天[LV.3]偶尔看看II
最近在看RBC之ABC，尝试将第8章CIA模型编程并用DYNARE模拟，程序运行报错，
There are 6 eigenvalue(s) larger than 1 in modulus
for 5 forward-looking variable(s)
The rank condition ISN'T verified!
错误使用 print_info (line 42)
Blanchard Kahn conditions are not satisfied: no stable equilibrium
出错 stoch_simul (line 98)
& & print_info(info, options_.noprint, options_);
出错 CIA (line 165)
info = stoch_simul(var_list_);
出错 dynare (line 180)
evalin('base',fname) ;
求助各位，不知如何修改，以下是code
var Y,C,K,H,W,R,M,p,g,
varexo e,u;
parameters&&beta, delta, B, theta, gamma,
beta = 0.99;
delta =0.025;
B=-2.5805;
theta =0.36;
gamma =0.95;
alpha=0.95;
1/beta=W/W(-1)*((1-delta)+R(+1));
B/(W*p)=-beta*(1/(p(+1)*C(+1)*g(+1)));
p*C=(M(-1)+g-1)/g;
K(+1)+M/p=(1-delta)*K+W*H+R*K;
W=(1-theta)*rho*(K/H)^
R=theta*rho*(K/H)^(theta-1);
Y=rho*(K^theta)*(H^(1-theta));
Y=C+K(+1)-(1-delta)*K;
log(rho)=gamma*log(rho(-1))+e;
log(g)=alpha*log(g(-1))+u;
Y = 1.2231;
C = 0.9095;
K = 12.544;
H = 0.3302;
W = (1-alpha)*Y/H;
R = alpha*Y/K;
p = 1.0995;
stderr 0.01;
stderr 0.01;
支持楼主：、
购买后，论坛将把您花费的资金全部奖励给楼主，以表示您对TA发好贴的支持
载入中......
总评分:&论坛币 + 2&
我也遇到相同的问题，楼主怎么解决的？
这种很难区分！有可能是参数取值导致的，也很有可能是模型输入错误导致的bk不满足。
个人方法是，先检查模型是否有问题，再调整参数，建议用其它文献常用得校准值，如果还是特征根问题那很有可能是模型问题。
无限扩大经管职场人脉圈！每天抽选10位免费名额，现在就扫& 论坛VIP& 贵宾会员& 可免费加入
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
如有投资本站或合作意向，请联系（010-）；
邮箱：service@pinggu.org
投诉或不良信息处理：（010-）
京ICP证090565号
论坛法律顾问：王进律师苹果/安卓/wp
积分 8, 距离下一级还需 2 积分
道具: 彩虹炫, 涂鸦板, 雷达卡, 热点灯下一级可获得
道具: 金钱卡
购买后可立即获得
权限: 隐身
道具: 金钱卡, 彩虹炫, 雷达卡, 热点灯, 涂鸦板
//封闭经济，无投资，无资本，无ZF支出，calvo定价，手工计算稳态，
var q Cj P i STj Wj MCh Yh yh ph PH U1 U2 BF//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 BF//很多参数在封闭经济下没有用，放在这里为开放经济讨论做准备
beta=0.99;
rou=3;//如果简化问题 rou=1
yeta=1;//封闭经济下
n=1;//封闭经济
theta=1.9;
lambdapai=2;
lambday=2;
yhbar=((sigma-1)/sigma)^(rou-niu);
//以下手工计算的稳态
Wjbar=(Mjbar*((gamma*((sigma-1)/sigma)^(rou^2-rou*niu))/(1-q))^(-1/kai))*(sigma-1)/
phbar=(sigma/(sigma-1))*W
ibar=1/qbar-1;
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=STj;//只有劳动这个生产要素
MCh=Wj;//只有劳动这个要素
//以下calvo定价
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));
//以下资源约束
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;//债券价格与名义利率关系
solve algo = 1;
var Mj = 2;//什么意思
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 294
Error in ==& dynare at 180evalin('base',fname) ;
支持楼主：、
购买后，论坛将奖励 10 元论坛资金给楼主，以表示您对TA发好贴的支持
载入中......
21:50:10 上传
出现奇异矩阵了，在网上搜也没有搜到什么结果，是我对理论的理解有问题，模型建错了，还是什么问题，请求帮助，比较急，谢谢！
我发的贴有问题吗，似乎大家都看不到？
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 发表于
请问您的问题解决了吗？我也碰到了同样的问题，能告诉我该怎么解决吗？谢谢说明你的稳态求解有问题。请将稳态带入方程确保每一个方程恒等于0
ahnulxy 发表于
说明你的稳态求解有问题。请将稳态带入方程确保每一个方程恒等于0对数线性化之后，哪些稳态不是都可以设为零的么？还需要再去重新求稳态吗？那我要是求出稳态之后该怎么用呢？我是菜鸟，正在试着学习，请您多多指导
Dynare forum,去看看吧，不过是英文的。这种情况下，很大可能是计算错误（稳态值），还有可能是timing写错了），如果不是计算错误什么的，那就是整个模型的结构错了。
榄外人L 发表于
对数线性化之后，哪些稳态不是都可以设为零的么？还需要再去重新求稳态吗？那我要是求出稳态之后该怎么用 ...您需要做的是先求出稳态，然而直接书写FOC条件，不需要线性化，最后在初始赋值模块里，将稳态值赋予给变量即可。自己线性化，容易出错，而且工作量极大，没有必要，让dynare自己帮你完成这个工作吧
榄外人L 发表于
对数线性化之后，哪些稳态不是都可以设为零的么？还需要再去重新求稳态吗？那我要是求出稳态之后该怎么用 ...在最后的stoch_simul 命令中 写order=1 就是线性化了
无限扩大经管职场人脉圈！每天抽选10位免费名额，现在就扫& 论坛VIP& 贵宾会员& 可免费加入
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
如有投资本站或合作意向，请联系（010-）；
邮箱：service@pinggu.org
投诉或不良信息处理：（010-）
京ICP证090565号
论坛法律顾问：王进律师