急求用MATLAB编程模拟马尔可夫链得到ARL值

本人写毕业论文用ARL特性分析小波动控制图的灵敏度，现在卡在MATLAB编程得到ARL值，请高人指点并传授经验，用其他模拟方法也行，急用，非常感谢！

function y=FEWMAARL(M0)
L=MO(1)；％控制限参数
r=MO(2)；％平滑系数
k=M0(3)；％间隔划分数，k为偶数
UCL=MO(4)；％上控制限
LCL=MO(5)；％下控制限
d=(UCL—LCL)／k；
m1=0：
％计算转移概率矩阵
for i=-k/2+1：1： k/2
k1=1：
ml=ml+1：
for j=-m+1:1:m
temp=((j*d)一(1一r)*(i一0．5)*d)／(r)；
T1=normcdf(temp，0，1)；
temp=((j一1)*d一(1一r)*(i--0．5)*d)／(r)；
T2=normcdf(temp，0，1)；
R(ml，k1)=T1一T2：
kl=kl+1：
end；
end； ·
PM=0*ones(1，2*m)；
PM(k／2)=1：
SS=PM*inv(eye(2*m)一R)*ones(2*m，1)；
y=SS；％返回ARL计算结果

楼主，请问一楼给的ARL代码是正确的吗？我也是这个的急求者。谢谢了！

围观中~