用马尔科夫链求解控制图的链长问题，
是求解EWMA控制图的平均链长，
其中转移概率P=Pr{1/r(LCL+(j一1)d一(1一r)(LCL+(i—o．5)d))≤
x≤{1/r(LCL+jd一(1一r)(LCL+(i一0·5)d))}
其中，x是服从标准正态分布。
具体代码如下：
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计算结果
想问一下，有关y=FEWMAARL(M0)是调用的什么函数啊？百度里面没有搜到，有没有专家给指点一下，谢谢啦！！！