马尔可夫链的重要性质以及两个例题如下:
再看另一个例题:
注意:
- 例5中有几个地方需要注意:
(1)为什么 P 11 = P 22 = p q + ( 1 − p ) ( 1 − q ) ,而 P 33 = p q + ( 1 − p ) P_{11} = P_{22} = pq + (1-p)(1-q),而P_{33} = pq + (1-p) P11=P22=pq+(1−p)(1−q),而P33=pq+(1−p)?这里要看到上述题目中的限制条件,“一个需要服务的顾客到达系统时发现系统内有3个顾客则该顾客即离去”。因此,当顾客数是3时,不会发生q或者1-q。
(2)为什么 P i j = 0 ( ∣ i − j ∣ ≥ 2 ) P_{ij} = 0(|i - j| \ge 2) Pij=0(∣i−j∣≥2)?请看下面第3点的论述。
(3)为什么 Δ t → 0 \Delta t \to 0 Δt→0的真是含义是什么?当 Δ t → 0 \Delta t \to 0 Δt→0时,在此时间间隔内,事件发生两次的概率极限为0。 - 思考一下,例6中哪里体现了马尔可夫过程的齐次性呢?所有的条件转移概率都用 P i j ( 1 ) P_{ij}(1) Pij(1)来代替,既体现了马尔可夫链的齐次性原理。
- 第一张图中的(1.7)和 (1.7)'阐述了一个重要的结论,“马氏链在任一时刻的一维分布由初始分布p(0)和n步转移概率矩阵所确定”,该结论结合了线性代数中的矩阵乘法和矩阵正交化、矩阵分解等知识,将马尔可夫过程中的概率问题转换为线性代数和高等数学问题来求解。
- 对第3点的进一步分析。结合如下两个矩阵A与B,考虑矩阵乘法的含义。
A B = [ x 11 x 12 x 13 ] [ y 11 y 12 y 13 y 21 y 22 y 23 y 31 y 32 y 33 ] = [ x 11 y 11 + x 12 y 21 + x 13 y 31 x 11 y 12 + x 12 y 22 + x 13 y 32 x 11 y 13 + x 12 y 23 + x 13 y 33 ] = C AB = \begin {bmatrix} x_{11}&x_{12} &x_{13} \end{bmatrix} \begin {bmatrix} y_{11}&y_{12}&y_{13}\\y_{21}&y_{22}&y_{23} \\y_{31}&y_{32}&y_{33}\end{bmatrix} = \\ \\ \begin {bmatrix} x_{11}y_{11} + x_{12}y_{21} + x_{13}y_{31} &x_{11}y_{12}+x_{12}y_{22} + x_{13}y_{32} &x_{11}y_{13}+x_{12}y_{23} + x_{13}y_{33} \end{bmatrix} = C AB=[x11x12x13] y11y21y31y12y22y32y13y23y33 =[x11y11+x12y21+x13y31x11y12+x12y22+x13y32x11y13+x12y23+x13y33]=C
首先,
A
代表初始的行矩阵,
B
矩阵中的向量元素是这样排列组合的:每一个元素的行列坐标,
表示从行序号转移到列序号的转移概率,代表的是变量在随机过程所有的两个元素的之间的转移概率,
x
和
y
的矩阵乘法运算,可以算出从初始位置到不同终止位置的转移概率
:
\color{red}首先,A代表初始的行矩阵,B矩阵中的向量元素是这样排列组合的:每一个元素的行列坐标,\\表示从行序号转移到列序号的转移概率,代表的是变量在随机过程所有的两个元素的之间的转移概率,\\x和y的矩阵乘法运算,可以算出从初始位置到不同终止位置的转移概率:
首先,A代表初始的行矩阵,B矩阵中的向量元素是这样排列组合的:每一个元素的行列坐标,表示从行序号转移到列序号的转移概率,代表的是变量在随机过程所有的两个元素的之间的转移概率,x和y的矩阵乘法运算,可以算出从初始位置到不同终止位置的转移概率:
[
A
]
[
B
]
=
[
x
1
x
2
x
3
]
[
y
11
y
12
y
13
y
21
y
22
y
23
y
31
y
32
y
33
]
=
[
x
1
y
11
+
x
2
y
21
+
x
3
y
31
x
1
y
12
+
x
2
y
22
+
x
3
y
32
x
1
y
13
+
x
2
y
23
+
x
3
y
33
]
\color{red} [A][B] = \begin {bmatrix} x_{1}&x_{2} & x_{3} \end{bmatrix} \begin {bmatrix} y_{11}&y_{12}&y_{13}\\y_{21}&y_{22}&y_{23} \\y_{31}&y_{32}&y_{33}\end{bmatrix} = \\ \begin {bmatrix} x_{1}y_{11} + x_{2}y_{21} + x_{3}y_{31} &x_{1}y_{12}+x_{2}y_{22} + x_{3}y_{32} &x_{1}y_{13}+x_{2}y_{23} + x_{3}y_{33} \end{bmatrix}
[A][B]=[x1x2x3]
y11y21y31y12y22y32y13y23y33
=[x1y11+x2y21+x3y31x1y12+x2y22+x3y32x1y13+x2y23+x3y33]
结合上图,即可明白矩阵乘法和马尔可夫转移概率矩阵的关系。也就是说,马尔可夫初始向量乘以马尔可夫一步转移矩阵后,所得到的新矩阵元素,正好是马尔可夫一步转移概率。因此,乘以两次马尔可夫一步转移矩阵后,正好是马尔可夫2步转移概率矩阵。
从C矩阵的元素下标可以看出,x与y中元素的关系是,x元素的列坐标等于y的行坐标。
P
{
X
n
=
a
j
}
=
∑
i
=
1
+
∞
P
{
X
n
=
a
j
,
X
0
=
a
i
}
=
∑
i
=
1
+
∞
P
{
X
n
=
a
j
∣
X
0
=
a
i
}
P
{
X
0
=
a
i
}
=
∑
i
=
1
+
∞
P
i
(
0
)
P
i
j
(
n
)
P \{X_n = a_j\} =\sum_{i =1} ^{+\infty} P \{X_n = a_j,X_0 = a_i\} = \\ \sum_{i =1} ^{+\infty} P \{X_n = a_j|X_0 = a_i\} P\{ X_0 = a_i\} = \sum_{i =1} ^{+\infty} P_{i}(0) P_{ij}(n)
P{Xn=aj}=i=1∑+∞P{Xn=aj,X0=ai}=i=1∑+∞P{Xn=aj∣X0=ai}P{X0=ai}=i=1∑+∞Pi(0)Pij(n)
结合此公式可以看出,因为
P
i
P_i
Pi是行向量,故
P
i
P_i
Pi中省略了行号而只有列号,用初始的
P
i
P_i
Pi向量乘以马尔可夫转移矩阵后,所得到的行向量,其元素的列坐标即为马尔可夫转移概率的目标下标。
5. 考虑马尔可夫转移概率的步长。
A
B
B
=
[
x
11
x
12
x
13
]
[
y
11
y
12
y
13
y
21
y
22
y
23
y
31
y
32
y
33
]
[
y
11
y
12
y
13
y
21
y
22
y
23
y
31
y
32
y
33
]
=
[
x
11
y
11
+
x
12
y
21
+
x
13
y
31
x
11
y
12
+
x
12
y
22
+
x
13
y
32
x
11
y
13
+
x
12
y
23
+
x
13
y
33
]
[
y
11
y
12
y
13
y
21
y
22
y
23
y
31
y
32
y
33
]
=
[
x
11
y
11
y
11
+
x
12
y
21
y
11
+
x
13
y
31
y
11
+
x
11
y
12
y
21
+
x
12
y
22
y
21
+
x
13
y
32
y
21
+
x
11
y
13
y
31
+
x
12
y
23
y
31
+
x
13
y
33
y
31
x
11
y
11
y
12
+
x
12
y
21
y
12
+
x
13
y
31
y
12
+
x
11
y
12
y
22
+
x
12
y
22
y
22
+
x
13
y
32
y
22
+
x
11
y
13
y
32
+
x
12
y
23
y
32
+
x
13
y
33
y
32
x
11
y
11
y
13
+
x
12
y
21
y
13
+
x
13
y
31
y
13
+
x
11
y
12
y
23
+
x
12
y
22
y
23
+
x
13
y
32
y
23
+
x
11
y
13
y
33
+
x
12
y
23
y
33
+
x
13
y
33
y
33
]
T
=
D
ABB = \begin {bmatrix} x_{11}&x_{12} &x_{13} \end{bmatrix} \begin {bmatrix} y_{11}&y_{12}&y_{13}\\y_{21}&y_{22}&y_{23} \\y_{31}&y_{32}&y_{33}\end{bmatrix} \begin {bmatrix} y_{11}&y_{12}&y_{13}\\y_{21}&y_{22}&y_{23} \\y_{31}&y_{32}&y_{33}\end{bmatrix} = \\ \begin {bmatrix} x_{11}y_{11} + x_{12}y_{21} + x_{13}y_{31} &x_{11}y_{12}+x_{12}y_{22} + x_{13}y_{32} &x_{11}y_{13}+x_{12}y_{23} + x_{13}y_{33} \end{bmatrix} \begin {bmatrix} y_{11}&y_{12}&y_{13}\\y_{21}&y_{22}&y_{23} \\y_{31}&y_{32}&y_{33}\end{bmatrix} = \\ \begin {bmatrix} x_{11}y_{11}y_{11} + x_{12}y_{21}y_{11} + x_{13}y_{31}y_{11} + x_{11}y_{12}y_{21} +x_{12}y_{22}y_{21} + x_{13}y_{32}y_{21} + x_{11}y_{13}y_{31} +x_{12}y_{23}y_{31} + x_{13}y_{33} y_{31} \\ x_{11}y_{11}y_{12} + x_{12}y_{21}y_{12} + x_{13}y_{31}y_{12} + x_{11}y_{12}y_{22} +x_{12}y_{22}y_{22} + x_{13}y_{32}y_{22} + x_{11}y_{13}y_{32} +x_{12}y_{23}y_{32} + x_{13}y_{33} y_{32}\\ x_{11}y_{11}y_{13} + x_{12}y_{21}y_{13} + x_{13}y_{31}y_{13} + x_{11}y_{12}y_{23} +x_{12}y_{22}y_{23} + x_{13}y_{32}y_{23} + x_{11}y_{13}y_{33} +x_{12}y_{23}y_{33} + x_{13}y_{33} y_{33} \end{bmatrix}^T = D
ABB=[x11x12x13]
y11y21y31y12y22y32y13y23y33
y11y21y31y12y22y32y13y23y33
=[x11y11+x12y21+x13y31x11y12+x12y22+x13y32x11y13+x12y23+x13y33]
y11y21y31y12y22y32y13y23y33
=
x11y11y11+x12y21y11+x13y31y11+x11y12y21+x12y22y21+x13y32y21+x11y13y31+x12y23y31+x13y33y31x11y11y12+x12y21y12+x13y31y12+x11y12y22+x12y22y22+x13y32y22+x11y13y32+x12y23y32+x13y33y32x11y11y13+x12y21y13+x13y31y13+x11y12y23+x12y22y23+x13y32y23+x11y13y33+x12y23y33+x13y33y33
T=D
从上式可以看出,马尔可夫转移概率其步长,即为矩阵乘法的次数。