我要用MATLAB求解控制图中的链长,用到的是马尔可夫链。应用二维马尔科夫链[0,c]*[a,b]来构建模型,其中:0<=Z<=c;a<=T<=b.
把Z的受控区间划分为m1个状态,每个状态长度为f=2c/2m1-1,除去第一个状态区间为l/2;同样的,T的受控区间也划分为m2个状态,每个状态的长度为v=2(b-a)/2m2-1,除去第一个状态区间为v/2.
矩阵R中的每一个元素转移概率p(i,j)(k,l)可表示为:
p(i,j)(k,l)=pr{Z(t)处于状态k,T(t)处于状态l/ Z(t-1)处于状态i,T(t-1)处于状态j }
=pr{min[b2,max(a1,b1)]<X<max[b1,min(a2,b2)]}
k不等于0,l不等于0;
Pr{a1<X<max[a1,min(a2,b2)]}
k不等于0,l等于0;
Pr{b1<X<max[b1,min(a2,b2)]}
k等于0, l不等于0;
Pr{X<min(a2,b2)}
k等于0, l等于0;
其中,a1,a2,b1,b2是变量,涉及到状态空间长度等,公式比较麻烦,就不写了。
最后计算ARL=P(I-R)I
其中,P为全为1的列向量,I为单位阵。
请问有哪位高手能不能给指点如何用MATLAB编写这个程序,由于本人初步涉及,实在不知道如何下手,拜托了!!!
有没有相关方面的程序代码,我参考一下也行,拜托啦!!!!