西电排队论大作业

一步转移概率矩阵的收敛特性。

陈灿枫 03124016

一步转移概率矩阵的特性应从以下两方面来分析：

第一：什么矩阵具有收敛特性即p^n=p^(n+1)。

第二：若一个转移矩阵（以下称一步转移概率矩阵为转移矩阵）有收敛性，那么其收敛的速度与什么有关呢？

对于一般的一步转移矩阵p 若有： an+1=anp=an

那么称该一步转移矩阵可收敛。anp=an

关于那些一步转移矩阵能够收敛我用matlab验证了几个比较具有代表性的矩阵：

1. 单位矩阵。

可以看到单位矩阵不具有收敛性。

2. 类单位矩阵。

类单位矩阵我们可以看到原本并非单位矩阵但是经过n次后也变为单位矩阵。由此可见此矩阵也不具有收敛特性。此类矩阵也易证明其不具有收敛性。

3. 一般一步转移概率矩阵（1）

我们可以看到经过18次后矩阵收敛到一个稳定的值。

4. 一般一步转移概率矩阵（2）

从这个矩阵我们可以看到该一步转移概率矩阵只经过了4次就趋于稳定收敛了。

有上述的四个例子我们能够总结：类单位矩阵单位矩阵是不具有收敛性的而一般的一步转移矩阵是有收敛性，而且收敛有快有慢。

那么是什么影响了一步转移概率矩阵的收敛的快慢呢？：

矩阵相乘行乘列的和列相同即行相加的和乘列行的和根据转移矩阵特性为1 所以也就收敛了。

若一开始的矩阵就是上面的转移矩阵那么他也就是收敛最快的因为他已经收敛了。我们再来对比（1）和（2）。不难发现矩阵（1）的列的差值比矩阵（2）的要大即矩阵（1）的方差要大的多。

那么我们就可以猜测是不是列的相似度越高其收敛的的速度也就越快呢。那么用什么指标去判断一个矩阵的列值得相似程度呢？

最先想到的就是矩阵的行列式的值，因为第一列为0的行列式值为0。不难看出矩阵收敛后的矩阵行列式值为0。

那么我们计算一下上述两个矩阵的行列式的值。

从上述的验证中可以看到矩阵1的行列式的绝对值为0.0255 而矩阵2的行列式绝对值为6*10-6远小于行列式1中的值而正好矩阵1的列值相似度要小于矩阵2。

上述只是总结性的验证，并没用理论的知识来证明该过程是否准确。那么行列式的值是否真的能刻画一步转移概率矩阵的收敛快慢呢？

我们先看类单位矩阵的行列式的值为1 而且不难证明所以得一步转移概率矩阵的行列式的值得绝对值都在[0,1]之间。假设一个n阶一步转移概率矩阵其行列式的表达式为：det(p)=a11*(-1)1+1det(c(11))+a12*

-1)1+2det(c(12))…a1n*(-1)1+ndet(c(1n))。

由上式可以看出若列值的差值越大那么行列式的值就取决于该列的值中的较大的值，，也会比较靠近0。

换个角度可以将单位矩阵看成1和很多无穷小ε组成。那么其决定因素就为1 那么其行列式的值就为1了。

所以我认为，利用一步转移概率矩阵的行列式的值来刻画矩阵的收敛快慢是可行的行列式的值越小其收敛的越快。

后记：到此也结束了由于这篇大作业总结是在较早时间完成的，：p=[0 1 0;0.

5 0 0.5; 0 1 0] 就是这个3阶的矩阵也是书上的一个例题的矩阵这个矩阵并不是上述我说的类单位矩阵或者是单位矩阵。而是一个一般的矩阵（就是有点对称）然而这个矩阵是没有办法收敛的其n次的值是在两个值之间循环跳动的。

我算了一下这个矩阵的det发现值为0 但是并没有上述验证中的列相同达到收敛的规律。但是其行列式的值也为0.之后我算了一下他的秩发现是2 也就是说秩的值小于阶的值而我之前举得例子中秩的值都是等于阶的值。

之后我又验证了一个矩阵p=[0.1 0.1 0.

1 0.7;0 0.2 0.

2 0.6;0 0 0.4 0.

6;0.1 0.1 0.

1 0.7] 。若不为满秩的话则可以算其各列的方差的平均值来进行比较单位矩阵的列平均方差为（n-1）/n 而其他的一步转移概率矩阵则介于0-(n-1)/n之间。