Markov 链法上证综合指数建模和仿真程序
%数据选择
A=[];
%数据前后做差
t=length(A);
u=2:1:t
TT(u-1,1)=A(u,1)-A(u-1,1);
m=length(TT);
%统计各状态出现频数
S0=0; S1=0; S2=0; S3=0;
for j=1:1:m
if TT(j)>=60
S0=S0+1;
elseif TT(j)<60 &&TT(j)>0
S1=S1+1;
elseif TT(j)>-60&& TT(j)<0
S2=S2+1;
elseif TT(j)<=-60
S3=S3+1;
end
End
%16种状态转移频数统计
p00=0;p01=0;p02=0;p03=0;
p10=0;p11=0;p12=0;p13=0;
p20=0;p21=0;p22=0;p23=0;
p30=0;p31=0;p32=0;p33=0;
%强涨至其他状态
for j=1:1:(m-1)
if TT(j)>=60&&TT(j+1)>=60
p00=p00+1;
elseif TT(j)>=60&&(TT(j+1)<60&&TT(j+1)>0)
p01=p01+1;
elseif TT(j)>=60&&(TT(j+1)<=0&&TT(j+1)>-60)
p02=p02+1;
elseif TT(j)>=60&&(TT(j+1)<=-60)
p03=p03+1;
%弱涨至其他状态
elseif (TT(j)<60&&TT(j)>=0)&&(TT(j+1)>=60)
p10=p10+1;
elseif (TT(j)<60&&TT(j)>=0)&&(TT(j+1)<60&&TT(j+1)>0)
p11=p11+1;
elseif (TT(j)<60&&TT(j)>=0)&&(TT(j+1)<0&&TT(j+1)>-60)
p12=p12+1;