-208-
}|{},,,|{
1111
ijPiiijP
nmnnnnmn
===
+−−+
L (2)
则称
},2,1,{ L=n
n
为一个马尔可夫链(简称马氏链),(2)式称为马氏性。
事实上,可以证明若等式(2)对 于
1
m
成立,则它对于任意的正整数
m
也成立。
因此,只要当
1=m 时(2)式成立,就可以称随机序列 },2,1,{ L
n
n
具有马氏性,
即
},2,1,{ L=n
n
是一个马尔可夫链。
定义 3 设
},2,1,{ L=n
n
是一个马氏链。如果等式(2)右边的条件概率与 n 无
关,即
)(}|{ mpijP
ijnmn
=
+
(3)
则称
},2,1,{ L=n
n
为时齐的马氏链。称
)(mp
ij
为系统由状态 i 经过 m 个时间间隔(或
m 步)转移到状态
的转移概率。(3)称为时齐性。它的含义是:系统由状态 i 到状态
的转移概率只依赖于时间间隔的长短,与起始的时刻无关。本章介绍的马氏链假定都
是时齐的,因此省略“时齐”二字。
2.2 转移概率矩阵及柯尔莫哥洛夫定理
对于一个马尔可夫链
},2,1,{ L=n
n
,称以 m 步转移概率
)(mp
ij
为元素的矩阵
))(()( mpmP
ij
= 为马尔可夫链的 m 步转移矩阵。当 1
m 时,记 PP
)1( 称为马尔可
夫链的一步转移矩阵,或简称转移矩阵。它们具有下列三个基本性质:
(i)对一切
i ∈, , 1)(0
≤ mp
ij
;
(ii)对一切
i
,
∑
∈
=
Ej
ij
mp 1)( ;
(iii)对一切
i ∈, ,
⎩
⎨
⎧
≠
=
==
时当
时当
, ji
ji
p
ijij
0
,1
)0(
δ
。
当实际问题可以用马尔可夫链来描述时,首先要确定它的状态空间及参数集合,然
后确定它的一步转移概率。关于这一概率的确定,可以由问题的内在规律得到,也可以
由过去经验给出,还可以根据观测数据来估计。
例 4 某计算机机房的一台计算机经常出故障,研究者每隔 15 分钟观察一次计算
机的运行状态,收集了 24 小时的数据(共作 97 次观察)。用 1 表示正常状态,用 0 表
示不正常状态,所得的数据序列如下:
1110010011111110011110111111001111111110001101101
111011011010111101110111101111110011011111100111
解 设
)97,,1( L
nX
n
为第 n 个时段的计算机状态,可以认为它是一个时齐马氏
链,状态空间
}1,0{=E ,编写如下 Matlab 程序:
a1='1110010011111110011110111111001111111110001101101';
a2='111011011010111101110111101111110011011111100111';
a=[a1 a2];
f00=length(findstr('00',a))
f01=length(findstr('01',a))
f10=length(findstr('10',a))
f11=length(findstr('11',a))
或者把上述数据序列保存到纯文本文件data1.txt中,存放在Matlab下的work
子目录中,编写程序如下:
clc,clear