用马尔科夫链求解控制图的链长问题,
是求解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)是调用的什么函数啊?百度里面没有搜到,有没有专家给指点一下,谢谢啦!!!


雷达卡


京公网安备 11010802022788号







