马尔可夫链中转移矩阵计算的MATLAB程序

A=[124455433325322241];
B=unique（A）;
C=sort（B‘ascend‘）;
m=length（C）;
n=length（A）;
k=4;%转移步数
for i=0:1:k-1 
    
x（i+1）=A（n-i）;%取到A序列最后一个状态
end

for j=1:1:m
    L=find（A==C（j））;%找到A中独立状态e的位置
        a=0;
        b=0;
        c=0;
        d=0;
        e=0;
    for i=1:1:length（L）
        if L（i）+k>n 
         break
        elseif A（L（i）+k）==C（1）
            a=a+1;
        elseif A（L