查看“︁轉移矩陣”︁的源代码

{{noteTA
|G1 = Math
|1 = zh-hans:概率; zh-hant:機率;
|2 = zh-hans:马尔可夫; zh-hant:馬可夫;
}}
{{distinguish|轉置矩陣}}
在[[数学]]中，'''随机矩阵'''（stochastic matrix）是用来描述一个[[马尔可夫链]]的转变的[[矩阵]]，亦称为'''概率矩阵'''（probability matrix）、'''转移矩阵'''（transition matrix）<ref>{{Cite book|url=https://link.springer.com/chapter/10.1007/0-387-21525-5_1|title=Applied Probability and Queues|date=2003|publisher=Springer, New York, NY|isbn=9780387002118|series=Stochastic Modelling and Applied Probability|pages=3–38|language=en|chapter=Markov Chains|doi=10.1007/0-387-21525-5_1|access-date=2018-03-03|archive-date=2019-06-29|archive-url=https://web.archive.org/web/20190629075122/https://link.springer.com/chapter/10.1007%2F0-387-21525-5_1}}</ref>、'''替代矩阵'''（substitution matrix）、'''马尔可夫矩阵'''（Markov matrix）或'''转移概率矩阵'''（transition probability matrix）。它的每一项都是一个表示[[概率]]的非负[[实数]]。它适用于[[概率论]]、[[统计学]]和[[线性代数]]，也在[[计算机科学]]和[[群体遗传学]]中使用。
有几种不同的定义和类型随机矩阵：

* '''右随机矩阵'''（right stochastic matrix）是实方阵，其中每一-{zh-cn:行;zh-hk:行;zh-mo:行;zh-sg:行;zh-tw:列}-求和为1。
* '''左随机矩阵'''（left stochastic matrix）是实方阵，其中每一-{zh-cn:列;zh-hk:列;zh-mo:列;zh-sg:列;zh-tw:行}-求和为1。
* '''双随机矩阵'''（doubly stochastic matrix）是非负实数方阵，每个行和列求和均为1。

同理，可以定义'''随机向量'''（也称为'''概率向量'''）为元素为非负实数且和为1的[[向量]]。因此，右随机矩阵的每一-{zh-cn:行;zh-hk:行;zh-mo:行;zh-sg:行;zh-tw:列}-（或左随机矩阵的每一-{zh-cn:列;zh-hk:列;zh-mo:列;zh-sg:列;zh-tw:行}-）都是一个随机向量。

在英语数学文献中的惯例是用概率的-{zh-cn:[[行向量]];zh-hk:[[行向量]];zh-mo:[[行向量]];zh-sg:[[行向量]];zh-tw:[[列向量]]}-和概率的右随机矩阵，而不用-{zh-cn:[[列向量]];zh-hk:[[列向量]];zh-mo:[[列向量]];zh-sg:[[列向量]];zh-tw:[[行向量]]}-和左随机矩阵，本文遵循此惯例。

==定义和性质==
随机矩阵描述了在一个[[有限集合|有限]][[概率空間|状态空间]] ''S'' 上的[[马尔可夫链]] <math>\boldsymbol{X}_{t}</math>。

如果在一个时间步长内从 <math>i</math> 到<math>j</math>移动的[[概率]]为 <math>\operatorname{Pr}(j|i)=P_{i,j}</math> ，随机矩阵 <math>P</math> 的第 <math>i</math> -{zh-cn:行;zh-hk:行;zh-mo:行;zh-sg:行;zh-tw:列}-，第 <math>j</math> -{zh-cn:列;zh-hk:列;zh-mo:列;zh-sg:列;zh-tw:行}-元素由 <math>P_{i,j}</math> 给出，例如，

:<math>P=\left[\begin{matrix}
P_{1,1}&P_{1,2}&\dots&P_{1,j}&\dots&P_{1,S}\\
P_{2,1}&P_{2,2}&\dots&P_{2,j}&\dots&P_{2,S}\\
\vdots&\vdots&\ddots&\vdots&\ddots&\vdots\\
P_{i,1}&P_{i,2}&\dots&P_{i,j}&\dots&P_{i,S}\\
\vdots&\vdots&\ddots&\vdots&\ddots&\vdots\\
P_{S,1}&P_{S,2}&\dots&P_{S,j}&\dots&P_{S,S}\\
\end{matrix}\right].</math>

由于从状态 <math>i</math> 到下一状态的概率总和必须是 1，这个矩阵是一个右随机矩阵，于是

:<math>\sum_j P_{i,j}=1 \,</math>

从 <math>i</math> 到 <math>j</math> 分两步转变的概率由然后由给定的 <math>P</math> 的平方矩阵的 <math>(i,j)</math> 号元素给出：

:<math>\left(P ^{2}\right)_{i,j} </math>

一般地，在由矩阵<math>P</math>给出的有限马尔可夫链上从任何状态转移到另一个状态的 ''k'' 步转移概率为 <math>P^k</math>。

初始分布为一个-{zh-cn:[[行向量]];zh-hk:[[行向量]];zh-mo:[[行向量]];zh-sg:[[行向量]];zh-tw:[[列向量]]}-。

平稳概率向量 <math>\boldsymbol{\pi}</math> 定义为不随转移矩阵的运用而变化的一个向量；也就是说，它定义为概率矩阵的左特征向量，其特征值为1：

:<math>\boldsymbol{\pi}P=\boldsymbol{\pi} </math>

{{Link-en|佩龙一弗罗宾尼斯定理|Perron–Frobenius theorem}}保证了每个随机矩阵都具有这样的向量，而特征值的最大绝对值始终为1。在一般情况下，可能有多个这样的向量。然而，对于具有严格正项的矩阵，该向量是唯一的，并可以观察到对任意 <math>i</math> 我们都有以下极限而求出，

:<math>\lim_{k\rightarrow\infty}\left(P^k \right)_{i,j}=\boldsymbol{\pi}_j </math> ， 

其中 <math>\boldsymbol{\pi}_{j}</math> 是-{zh-cn:行;zh-hk:行;zh-mo:行;zh-sg:行;zh-tw:列}-向量 <math>\boldsymbol{\pi}</math> 的第 <math>j</math> 个元素。在其他方面，这表示处在状态 <math>j</math> 下的长期概率与初始状态 <math>i</math> 是独立的。这两种计算得到相同的稳定向量是遍历定理的一种形式，在各种各样的[[耗散系統|耗散动力系统]]广泛成立：该系统随着时间演变到[[定态]]。

直观地看，随机矩阵表示一个马尔可夫链；对概率分布应用随机矩阵，就是将原始分布的概率质量进行重新分布，同时保持其总质量。如果反复应用此过程，分布就会收敛为马尔可夫链的平稳分布。

==應用==

轉移矩陣可用以表示[[機率]](或變化比率)，而矩陣相乘的結果可用以預測未來事件發生的[[機率]]。

==性質==

設<math>\mathbf{A}</math>、<math>\mathbf{B}</math>為二個''n''×''n''階轉移矩陣，則以下亦為轉移矩陣:

*<math>\mathbf{AB}</math>
*<math>\mathbf{A}^2</math>
*<math>\frac{1}{2}(\mathbf{A}+\mathbf{B})</math>

== 范例：猫和老鼠 ==

假设你有一个计时器和五个相邻的格子排成一行，零时刻有一只猫在第一个格子中，而一只老鼠在第五个格子中。在计时器增加的时候猫和老鼠都会随机跳到一个相邻的格子中。例如，如果猫在第二个格子，老鼠在第四个，''在计时器增加后，猫会出现在第一个格子'''且'''老鼠会出现在第五个格子''的概率为1/4。如果猫在第一个格子而老鼠在第五个，那么计时器增加后，猫会出现在第二个格子且老鼠会出现在第四个的概率为1。当它们处于同一个格子的时候，猫会吃掉老鼠，游戏结束。[[随机变量]] ''K'' 给出了老鼠仍留在游戏中的时间步长。

表示这个包含五种位置组合 (猫,鼠) 的状态的游戏的[[马尔可夫链]]为：

* 状态 1：(1,3)
* 状态 2：(1,5)
* 状态 3：(2,4)
* 状态 4：(3,5)
* 状态 5：游戏结束：(2,2), (3,3) & (4,4) 

我们使用一个随机矩阵来表示这个系统的[[马尔可夫链|转移概率]]（这个矩阵中的行和列用上面提到的可能状态来索引），

:<math> P =
\begin{bmatrix}
    0  & 0    & 1/2  &   0  & 1/2 \\
    0  & 0    & 1    &   0  & 0 \\
  1/4  & 1/4  & 0    & 1/4  & 1/4 \\
    0  & 0    & 1/2  &   0  & 1/2 \\
    0  & 0    & 0    &   0  & 1
\end{bmatrix} </math>

===长期平均===

无论初始状态是什么，猫最终都会抓到老鼠（概率为1），且极限为稳态 '''&pi;''' = (0,0,0,0,1)。要计算随机变量 Y 的长期平均或期望值。对每种状态 S<sub>j</sub> 和时间 t<sub>k</sub>，都有 Y<sub>j,k</sub>·P(S=S<sub>j</sub>,t=t<sub>k</sub>) 的贡献。生存与否可以视作一个二值变量，Y=1 代表生存状态而 Y=0 代表终止状态。Y=0 的状态不对长期平均有贡献。

===位相型表示===
[[Image:Mousesurvival.jpg|thumb|right|350px|老鼠的生存函数。老鼠至少在第一个时间步长存活。]]

由于状态 5 是一个吸收态，吸收对时间的分布为{{le|离散位相型分布|Discrete phase-type distribution}}。假设系统从状态 2 开始，表示为向量 <math>[0,1,0,0,0]</math>。老鼠死亡后的状态不会对生存平均产生影响，所以状态五可以忽略。初始状态和转移矩阵可以化简为，

:<math>\boldsymbol{\tau}=[0,1,0,0]</math>

以及，

:<math>T=\begin{bmatrix}
    0  & 0    & 1/2  &   0\\
    0  & 0    & 1    &   0\\
  1/4  & 1/4  &   0  & 1/4\\
    0  & 0    & 1/2  &   0
\end{bmatrix}</math>；而 <math>(I-T)^{-1}\boldsymbol{1}
=\begin{bmatrix}2.75 \\ 4.5 \\ 3.5 \\ 2.75\end{bmatrix}\,,</math>
其中 <math>I</math> 为[[单位矩阵]]，<math>\mathbf{1}</math> 表示全为1的-{zh-cn:列;zh-hk:列;zh-mo:列;zh-sg:列;zh-tw:列}-矩阵，进行状态的相加。

由于每个状态都占据一个时间步长，老鼠生存时间的期望就是在所有生存状态和时间步长中占据的概率之[[矩阵多项式|和]]，

:<math>E[K]=\boldsymbol{\tau}(I+T+T^2+\cdots)\boldsymbol{1}=\boldsymbol{\tau}(I-T)^{-1}\boldsymbol{1}=4.5 </math>

其高阶矩为

:<math>E[K(K-1)\dots(K-n+1)]=n!\boldsymbol{\tau}(I-{T})^{-n}{T}^{n-1}\mathbf{1}\ </math>

==参见==
* [[Muirhead's inequality]]
* {{Link-en|佩龙一弗罗宾尼斯定理|Perron–Frobenius theorem}}
* [[密度矩陣]]
* {{Link-en|双随机矩阵|Doubly stochastic matrix}}
* [[Discrete phase-type distribution]]
* {{Link-en|概率自动机|Probabilistic automaton}}
* [[Models of DNA evolution]]
* {{Link-en|马尔可夫核|Markov kernel}}，随机矩阵在连续状态空间的等价形式

==参考文献==
{{Reflist}}
*  G. Latouche, V. Ramaswami. ''Introduction to Matrix Analytic Methods in Stochastic Modeling'', 1st edition. Chapter 2: PH Distributions; ASA SIAM, 1999.

{{概率分布}}

[[Category:矩陣]]
[[Category:马尔可夫模型]]