我要用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编写这个程序，由于本人初步涉及，实在不知道如何下手，拜托了！！！
有没有相关方面的程序代码，我参考一下也行，拜托啦！！！！