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计算结果
扫码加好友，拉您进群
全部版块
我的主页
收藏
