matlab dynare问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

###Ramsey Model,BK's Method
bk< - function(mA,mB,mD,nk);
{
n=nrow(mA);
B=solve(mA)%*%mB;## B:B~
D=solve(mA,mD);##D:D~
B11=B[1:nk,1:nk];

Lamda=diag(eigen(B)Svalues[n:1]);
Lamda1=Lamda[1:nk,1:nk];
Lamda2=Lamda[(nk+1):n,(nk+1):n];

P=eigen(B)$vectors[ ,n:1];
Pi= solve(P);

J=Pi%*%B;
M=Pi%*%D;
P11=Pi[1:nk,1:nk];
P22=Pi[(nk+1):n,(nk+1):n];
P12=Pi[1:nk,(nk+1):n]
P21=Pi[(nk+1):n,1:nk]


F=matrix(NA,n,nk);
G=matrix(NA,n,ncol(D));

F[(nk+1):n,]=-solve(P22)%*%P21;
G[(nk+1):n,]=-solve(P22)%*%solve(Lamda2)%*%M[(nk+1):n,];

F[1:nk,]=B11-B12%*%solve(P22)%*%P21;
G[1:nk,]=D[1:nk,]-B12%*%solve(P22)%*%solve(lamda2)%*%M[(nk+1):n,];
out<-NULL;
out % F=F;
out % G=G;
return(out);
}
alpha=0.7;beta=0.95;sigma=2;
Kbar=(beta*(1-alpha))^(1/alpha);##Kbar=0.1664
Ybar=(Kbar)^(1-alpha)-Kbar;##Cbar=0.4175

mA=matrix(c(Kbar/Ybar,alpha,0,sigma),2,2);
mB=matrix(c(1-alpha,0,-Cbar/Ybar,sigma)2,2);
nk=1;
bk(mA,mB,mD,nk=1);


请问为什么会出现
??? Error using ==> dynare at 114
DYNARE: preprocessing failed
以后在编写程序时如何避免这种问题的出现？？？
谢谢！！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：MATLAB dynare matla atlab ARE matlab

人工置顶，求初学求大神

lyngy1314520 发表于 2013-7-17 19:40
人工置顶，求初学求大神

这是Matlab程序，不在dynare环境下运行。