hzcmaster 发表于 2015-5-23 13:05 
你没理解我的意思,你这里的h(a)记录的是每一次序列的最后一个值,这个值跟VAR没有任何关系,你应该写
v ...
所以我把同学你后半部分的var加进去之后,
clear
clc
mcmc=1000;
p=zeros(1,mcmc+1);
h=zeros(1,mcmc+1);
p(1)=6.1089;
alpha=0.05;
for a=1:mcmc
for k=1:mcmc
epsilon=normrnd(0,1);
p(k+1)=p(k)+p(k)*(1.94152E-05+0.000638662*epsilon);
end
p(1)=[];
lo=prctile(p,alpha);
up=prctile(p,alpha);
var=mean(p-quantile(p,alpha));//①
h(a)=var;
end
vars=sort(h);
varF=abs(vars(mcmc*0.05));//②
这样是不是有错啊?感觉1和2有冲突呢,也不知道我理解的对不对