当前位置:文档之家› 第1章离散时间的马尔可夫链

第1章离散时间的马尔可夫链

第1章离散时间的马尔可夫链
第1章离散时间的马尔可夫链

第1章 离散时间的马尔可夫链

§1 随机过程的基本概念

定义1 设(,,)P ΩF 是概率空间,(, )E E 是可测空间,T 是指标集. 若对任何

t T ∈,有:t X E Ω→,且t X ∈F E ,则称{}(), t X t T ω∈是(, , )P ΩF 上的取值于

(,)E E 中的随机过程,

在无混淆的情况下简称{(), }t X t T ω∈为随机过程,称(,)E E 为状态空间或相空间,称E 中的元素为状态,称T 为时间域. 对每个固定的ω∈Ω,称()t X ω为{}(), t X t T ω∈对应于ω的轨道或现实,对每个固定的t T ∈,称()t X ω为E 值随机元. 有时()t X ω也记为

()()(, )t t X X X t X t ωω===

设T ?R ,{}, t t T ∈F 是F 中的一族单调增的子σ代数(σ代数流)

,即 ① t t T ?∈??F F ,且t F 是σ代数; ② , , s t s t T s t ?∈

X t T ∈?∈F E ,则称{}, t X t T ∈是{}t F 适应的随机过程,或适应于{}

t F 的随机过程. 特别地,若令

1(, , )(())t s s s t

s T

X s t s T X σσ-≤∈≤∈ F E

是由{}, , s X s t s T ≤∈所生成的σ代数,则{}, t X t T ∈是{}t F 适应的随机过程.

当1(, )(, )E =R E B 时,称{}, t X t T ∈为实值随机过程; 当(, )(, )E =C

C E B 时,称{}, t X t T ∈为复值随机过程;

当(, )(, )n n

E =R E B

时,称{}, t X t T ∈为n 维随机过程;

当E 是可列集(有限集)时,称{}, t X t T ∈为可列(有限)随机过程; 当, T =+R R 或[], a b 时,称{}, t X t T ∈为连续参数的随机过程; 当T =Z 或+

Z 时,称{}, t X t T ∈为离散参数的随机过程(随机序列);

当, (), n

n

n

T =+R R Z 或() (2)n

n ≥+Z 时,称{}, t X t T ∈为随机场.

随机过程的四种类型:

(1)指标集T 离散,状态空间E 离散的随机过程; (2)指标集T 离散,状态空间E 连续的随机过程; (3)指标集T 连续,状态空间E 离散的随机过程; (4)指标集T 连续,状态空间E 连续的随机过程.

然而,以上分类是表面的,更深刻的是按随机过程的概率结构而分类.

例如:马尔可夫(Markov )过程、平稳过程、独立增量过程、二阶矩过程、正态过程、泊松(Poisson )过程、生灭过程、分枝过程、更新过程、鞅等.

对于随机过程{}, t X t T ∈而言,可以这样设想,有一个作随机游动的质点M ,以t X 表示在时刻t 质点M 的位置,于是{}, t X t T ∈描绘了质点M 所作的随机运动的变化过程,一般把“t X x =”形象地说成“在时刻t 质点M 处于状态x ”.

定义2 设{}, t X t T ∈是概率空间(,,)P ΩF 上的、以(,)E E 为状态空间的随机过程,T =+R (或R 或直线上的任一区间).如果A ?∈E ,有

{}(,):(,), (,)()t t T X t A T ωωω∈?Ω∈∈?B

F

则称{}, t X t T ∈是可测的.

设{}, t t T ∈F 是F 中的一族单调增的子σ代数. 如果, ,t T A ?∈∈E 有

[]{}[]()(,):(,)0, , (,)0, t

u u t X t A t ωωω∈?Ω∈∈?B F

则称{}, t X t T ∈关于{}, t t T ∈F 循序可测. 命题1 设:t X E Ω→, ()t t

X t T ∈?∈F E ,{}, t t T ∈F 是F 中的一族单调增

的子σ代数. 如果{}, t X t T ∈关于{}, t t T ∈F 循序可测,则{}, t X t T ∈是可测的. 定义3 设{}, t X t T ∈是随机过程,称

{}{}():(), , t t t F x P X x P X x x t T ωω≤=≤∈?∈R

为随机过程{}, t X t T ∈的一维分布函数;称

{}1212,12121212(,),, ,, ,t t t t F x x P X x X x x x t t T ≤≤∈?∈R

为随机过程{}, t X t T ∈的二维分布函数;一般地,称

{}

1212,,,1212(,,,),,,n n t t t n t t t n F x x x P X x X x X x ≤≤≤

1212,,,, ,,,n n x x x t t t T ∈?∈R

为随机过程{}, t X t T ∈的n 维分布函数;而称

{}

12,,,1212(,,,):,,,,1n t t t n n F x x x t t t T n ∈≥ F

为随机过程{}, t X t T ∈的有限维分布函数族.

随机过程{}, t X t T ∈的有限维分布函数族F 具有下列性质:

1. 对1n ?≥,12,,,n t t t T ?∈ ,及12,,,n t t t 的任意排列12,,,n i i i t t t ,有

121212,,,,,,12(,,,)(,,,)i i i n n n t t t i i i t t t n F x x x F x x x = (对称性)

2. 对1m n ?≤≤,有

12121,,,12,,,,,,12(,,,)(,,,,,,)m m m n t t t m t t t t t m F x x x F x x x +=+∞+∞ (相容性)

注 若知道了随机过程{}, t X t T ∈的有限维分布函数族F ,便知道了这一随机过程中任意有限个随机变量的联合分布,也就可以完全确定它们之间的相互关系. 可见,随机

过程的有限维分布函数族能够完整地描述随机过程的统计特征. 但是在实际问题中,要知道随机过程的有限维分布函数族是不可能的,因此,人们想到了用随机过程的某些数字特征来刻画随机过程.

定义4 设{}, t X t T ∈是随机过程,称

()d ()()(d ), t t t m t EX x F x X P t T ωω+∞

-∞

Ω

==∈?

?

为{}, t X t T ∈的均值函数;称

2()(()), t t D t DX E X m t t T =-∈

为{}, t X t T ∈的方差函数;称

(,)Cov(,)(())(()), ,s t s t C s t X X E X m s X m t s t T =--∈

为{}, t X t T ∈的协方差函数;称

(,)(), ,s t R s t E X X s t T =∈

为{}, t X t T ∈的相关函数.

注 若{}, t X t T ∈是复值随机过程,则方差函数的定义为

2

()(), t D t E X m t t T =-∈

协方差函数的定义为

(,)(())(()), ,s t C s t E X m s X m t s t T =--∈

相关函数的定义为

(,)(), ,s t R s t E X X s t T =∈

性质

(1)(,)(), C t t D t t T =∈;

(2)(,)(,)()(), ,C s t R s t m s m t s t T =-∈; (3)若()0m t ≡,则 (,)(,), ,C s t R s t s t T =∈.

§2 马尔可夫链的定义

在实际中有一类很广泛的随机过程,其特点是:过去只影响现在,而不影响将来. 这种随机过程称为马尔可夫过程. 状态离散的马尔可夫过程称为马尔可夫链,本章介绍时间离散的马尔可夫链(简称马尔可夫链).

马尔可夫(Markov )过程的研究始于1906年,是随机过程的一个重要分支,它在近代物理、生物学、管理科学、信息处理、自动控制、金融保险等方面有着许多重要应用. 在本章中,无特别声明我们总是假设: 1. 参数集合{}0,1,2,T = ;

2. 状态空间{}0,1,2,S = 或{},2,1,0,1,2,S =-- 或其子集.

定义5 设{}, 0n X n ≥是定义在概率空间(, , )P ΩF 上的随机过程,状态空间为S . 若对于任意的1n ≥及任意的整数120,n t t t t ≤<<<< 12,,,,n i i i j S ∈ ,有

{}{}1212,,, ()n n t t t t n t t n P X j X i X i X i P X j X i =======* 则称{}, 0n X n ≥为马尔可夫链,简称马氏链. 等式()*称为马氏性或无后效性,且假定

()*式两端的条件概率都有意义(以下涉及到条件概率的式子都作类似的假定).

定理1 随机过程{}, 0n X n ≥是马尔可夫链的充要条件是对任意的1n ≥及任意的

12,,,,n i i i j S ∈ ,有

{}{}111221,,,n n n n n n P X j X i X i X i P X j X i ++=======

§3 转移概率

对于马尔可夫链{}, 0n X n ≥,描述它概率性质最重要的是它在时刻m 的一步转移概率{}

1(), ,ij m m p m P X j X i i j S +===∈.

马尔可夫链是描述某些特定的随机现象的数学模型,而产生这种特定的随机现象的具体模型一般称为系统,因此我们经常把事件{}m X i =说成是在时刻m 时系统处于状态

i ,把{}1m m P X j X i +==说成已知在时刻m 时系统处于状态i ,而在时刻1m +时系统

转移到状态j 的概率等等.

定义6 设{}, 0n X n ≥是状态空间为S 的马尔可夫链,称

{}()(), ,n ij m n m p m P X j X i i j S +===∈

为系统在时刻m 时处于状态i 的条件下, 经n 步转移到状态j 的n 步转移概率,简称时刻m 的n 步转移概率.

显然,()

()n ij p m 具有下列性质: (1)()()0, ,n ij p m i j S ≥∈; (2)

{}()

()1, n ij

m n m j S

j S

p

m P X j X i i S +∈∈====∈∑∑.

上述性质说明了,对于任意给定的i S ∈及0, 1m n ≥≥,{}

()

(), n ij p m j S ∈是一个

概率分布. 规定:

(1)(1)

()()ij ij p m p m =;

(2)(0)

1, ,

()0, .ij

ij i j p m i j δ=?==?≠?

若()

()n ij p m 与m 无关,则称{}, 0n X n ≥是时齐的或齐次的马尔可夫链. 此时,记

()()(), ,, 1n n ij ij p p m i j S n =∈≥;一步转移概率记为(1), ,ij ij p p i j S =∈.

对时齐的马尔可夫链{}, 0n X n ≥,有

{}{}()

0, ,, 0n ij m n m n p P X j X i P X j X i i j S m +======∈?≥

以下恒设马尔可夫链{}, 0n X n ≥是时齐的,并简称为马尔可夫链.

性质 马尔可夫链{}, 0n X n ≥的n 步转移概率()

n ij p 具有下列性质

(1)()

,,0n ij i j S p ?∈≥; (2)()

,

1n ij

j S

i S p

∈?∈=∑.

定理2(Chapman-Kolmogorov ) 设()

n ij p 是马尔可夫链{}, 0n X n ≥的n 步转移概率,

则 ,,,0i j S m n ?∈≥,有

()()()

m n m n ij ik kj k S

p p p +∈=∑ (C-K 方程)

证明 {}{}{}()

00

, m n ij m n m

m n

k S

p P X j X i P

X k X j X i +++∈=======

{}{}{}0

, , m m n m m n k S

k S

P X k X j X i P X k X j X i ++∈∈========∑

{}{}00 ,m

m n m k S P X

k X i P X j X i X k +∈=

=====∑ {}{}0 m

m n m k S

P X

k X i P X j X k +∈=====∑

{}{}()()

00 m n m

n ik kj k S

k S

P X

k X i P X j X k p p ∈∈=

=====∑∑

定理3 马尔可夫链{}, 0n X n ≥的一步转移概率ij p 可以确定所有的n 步转移概率

()

n ij

p . 证明 由C-K 方程,显然. 记()

()

(1),,(), ()n n ij i j S ij i j S p p ∈∈===P

P P . 称()n P 为马尔可夫链{}, 0n X n ≥的n

步转移矩阵,称P 为马尔可夫链{}, 0n X n ≥的(一步)转移矩阵. 此时,C-K 方程可表

示为()

()()m n m n +=P

P P 且()n n =P P .

定义7 设{}, 0n X n ≥是马尔可夫链,对任意的0n ≥,称{}π()i n n P X i ==,

i S ∈为绝对概率,特别地,称{}0π(0), i P X i i S ==∈为初始概率.

显然,绝对概率和初始概率具有下列性质:

π()0, π()1, i i i S

n i S n ∈≥∈??

?=??∑ π(0)0, π(0) 1 i i i S i S

∈≥∈???=??∑

故对任意0n ≥,{}π(), i n i S ∈是概率分布,通常称为绝对(概率)分布;特别,

{}π(0), i i S ∈称为初始(概率)分布. 记{}{}(0)π(0), ()π()i i i S i S n n ∈∈∏=∏=.

定理4 设{}, 0n X n ≥是马尔可夫链, 则它的任意有限维概率分布完全由初始分布和一步转移概率决定.

证明 对任意的1n ≥,任意的整数120n t t t ≤<<< 及任意的12,,,n i i i S ∈ ,有

{}

1212,,,n t t t n P X i X i X i ===

{}{}1

2

1

2

,,,,n

t t t n

i S P X i X i X i X i ∈===== {}{}1

2

1

2

,,,,n

t t t n

i S

P X i X i X i X i ∈===== {}1

2

1

2

,,,,n

t t t n

i S

P X i X i X i X i ∈=====∑ {}{}{}

1

2

10

10

20

1,t t t i S

P X i P X i X i P X i X

i X i ∈=======∑

{}

11011,,,n n t n t t n P X i X i X i X i --==== {}{}{}

12101021t t t i S P X i P X i X i P X i X i ∈======∑

{}

11n n t n t n P X i X i --==

11211121()()()

π(0)n n n n

t t t t t i ii i i i i i S p p p ----∈=∑

§4 若干例子

定义8 设012, , ,ξξξ 是取整数值的独立同分布的随机变量序列,令0

n

n k

k X ξ

==∑,

则称{}, 0n X n ≥为随机游动.

定理5 随机游动{}, 0n X n ≥是时齐的马尔可夫链. (证明略)

例1(无限制的随机游动)若随机游动{}, 0n X n ≥的状态空间为{,2,1,S =--

}0,1,2, ,且转移概率为

, 1,

, 1,

, 0, ij p j i p q j i i j S =+??

==-∈???

其它.

其中01,1p q p <<=-. 求:n 步转移概率()

n ij p .

解 设在n 步转移中,向右移动x 步,向左移动y 步,则经n 步从i 到达j ,x 和y 应满足,x y n x y j i +=-=-. 所以

(())2, (())2x n j i y n j i =+-=--

由于, x y 只能取正整数,故()n j i +-与()n j i --必须是偶数. 又因在n 步转移中

有x 步向右移动,故经n 步转移由i 到j 共有C x

n 种方式,于是

()()()222()C ,()0,

()n j i n j i n j i n n ij

p q n j i p n j i +-+---??+-=??+-?偶奇 特别

222()

C ,0,

n n n

n n ii p q n p n ??=???偶奇

例2(带有一个吸收壁的随机游动) 设{}, 0n X n ≥是随机游动,其状态空间为

{}0,1,2,S = ,若转移概率为

1, 0,, 1, , (01, 1), 10, ij i j p j i p i j S p p q q j i ==??=+?=∈<<+=?=-???,其它.

则称{}, 0n X n ≥为带有吸收壁0的随机游动. 其转移矩阵为

10000000000q p q p q ?? ? ? ?= ? ? ???

P

例3(带有两个吸收壁的随机游动) 设{}, 0n X n ≥是随机游动,其状态空间为

{}0,1,2,,S b = ,其转移概率为

1, 0,1, ,, 1,, (01, 1) , 10, ij i j i j b p p j i i j S

p p q q j i ==??==??

==+∈<<+=??=-???,其它.

其转移矩阵为

100

00000000000000000000000001q p q p q p ??

?

? ?

=

? ? ?

?

??

?

P

例4(带有一个反射壁的随机游动) 设{}, 0n X n ≥是随机游动,其状态空间为

{}0,1,2,S = ,其转移概率为

1,0,1,, 1, , (01, 1), 10, ij i j p j i p i j S p p q q j i ==??=+?=∈<<+=?

=-???,其它.

其转移矩阵为

010********q p q p q ??

? ? ?= ? ? ???

P

例5(带有两个反射壁的随机游动) 设{}, 0n X n ≥是随机游动,其状态空间为

{}, 0n X n ≥,其转移概率为

1, 0,1,1, ,1,, 1,

, (01, 1) , 10, ij i j i b j b p p j i i j S p p q q j i ==??==-??

==+∈<<+=??=-???,其它.

其转移矩阵为

01000000000000000000000000010q p q p q p ??

? ? ?= ? ?

? ? ???

P

例6(带有弹性壁的随机游动) 设{}, 0n X n ≥是随机游动,其状态空间为

{},2,1,0,1,2,S =-- ,其转移概率为

, 1,, ,, (0,,1, 1), 1,0, ij p j i r j i p i j S p q r p q r q j i =+??=?=∈<<++=?

=-???其它.

其转移矩阵为

r p q r p

q r ?? ? ?

?= ?

?

???

P 例7 设{}, 0n X n ≥是只有两个状态({}0, 1S =)的随机游动,其转移矩阵为

1, 01, 011p

p p q q q -??=<<<< ?-??P

这里p q +未必等于1. 求()

n P ,0lim π()n n →∞

,1lim π()n n →∞

.

解 ()()

()

0001()()10

11n n n n n p p p p ??= ???P 由C-K 方程,得

()

(1)(1)(1)000000000110n n n n k k k S

p p p p p p p ---∈=

=+∑

(1)(1)(1)

000000

(1)(1)(1)n n n p p q p q p q p ---=-+-=+-- (1)n q p

p q p q p q =

+--++ 同理可求()()()

011011

, , n n n p p p ,故 ()

(1)(1)(1)(1)n

n n n

n q p

p p p q p q p q p q p q p q

q q

p q p q p q p q p q

p q p q ?

?+----- ?++++ ?= ?---+-- ?++++?

?

P

设初始分布为{}01(0)π(0), π(0)∏=,则

()

()

0000

1

1

0π()π(0)π(0)n n n p p =+ 0π(0)(1)n q p p q p q p q ??=+-- ?++??0(1π(0))(1)n q q p q p q p q ??+---- ?++??

0(1)π(0)n q p p q p q p q ??=

+--- ?++?? 故0lim π()n q n p q →∞=+,同理 1lim π()n p

n p q

→∞=+. 例8 设{}, 0n X n ≥是马尔可夫链,状态空间{}1,2,3S =,转移矩阵为

1434014141203414?? ?= ? ???

P

初始分布为{}{}123π(0), π(0), π(0)16, 12, 13=,求

(1)二步转移矩阵(2)

P

(2){}132,1P X X ==;

(3){}

32102,2, 22P X X X X =≠≠=. 解 (1)(2)

214340143401438314141214141215814034140341431638716??????

??? ?=== ??? ? ??? ???????

P

P .

(2)由全概率公式、乘法公式、马氏性,得 {}

{}{}313013012,12, 1i P X X P X i P X X X i

=====

===∑

{}{}{}3010311

212i P X i P X X i P X X =======∑

3

(2)2211

π(0)i i i p p ==∑

(注:也可直接由定理4得到) 13111311642434816??=?

+?+??= ??

? (3)由加法公式、乘法公式、马氏性,得 {

}

32102,2, 22P X X X X =≠≠=

{}{}3

21032102,1, 122,3, 12P X

X X X P X X X X =====+====

{}{}321032102,1, 322,3, 32P X X X X P X X X X +====+====

211112211332233112233332p p p p p p p p p p p p =+++ 11311390044424464

=??+++??= 例9 设随机游动{}, 0n X n ≥的状态空间{}0,1,2,,S b = ,其中0和b 是吸收壁,

初始状态为a ,转移概率为

1,0 ,

, 1,01, 1, 1,0, ij i j i j b p j i p p q p q j i ====??=+?=<<=-?

=-???或其它.

求质点被0点吸收的概率.

解 设i u 表示质点初始位置为i 而被0点吸收的概率,则

0111, 0, , 1,2,,1b i i i u u u pu qu i b +-===+=-

由于 1 p q +=,故11()()i i i i p u u q u u +--=-.

(1)若12 p q ==,则11i i i i u u u u C +--=-=(常数),故10, .i i i u C u u iC u -=+=+ 因01, 0b u u ==,故1C b =-,从而1i i u b =-+,故1a a

u b

=-. (2)若 p q ≠,则

111()()()(1)i i i i i u u q p u u q p u +--=-==-

()11(1)i

i i u u q p u +=+-

从而

()

1

11(1)b b b u u q p u --=+- ()

2

121(1)b b b u u q p u ---=+-

()0

101(1)u u q p u =+-

相加,得 011()(1)1()

b

b q p u u u q p -=+--,

所以

11()

11()b

q p u q p --=-

-

011()(1)1()

a

a q p u u u q p -=+--

1()1()11()1()a b q p q p q p q p ??--=+- ?--??

1()11()

a

b

q p q p -=--

§5 状态的分类

设{}, 0n X n ≥是马氏链,S 为状态空间,ij p 为转移概率. 定义9 设A S ?,令

{}min :, 1, 1, , 1, n n A n n X A n n X A

T n X A

∈≥?≥∈?=?

+∞?≥??使有 称A T 为{}, 0n X n ≥首达A 的时刻(或A 的首达时),A T 是随机变量,其可能值为

1,2,,.+∞ 当A 为单点集{}j 时,记{}j T 为j T . 令

{}()0, 1n ij j f P T n X i n ===≥

则()n ij f 表示系统从状态i 出发, 经n 步首次到达状态j 的概率. 利用乘法公式、马氏性易知:

{}()1210, ,,, n ij n n f P X j X j X j X j X i -=≠≠≠==

1121121, 1n n i i

i i i j

i j i j

i j

p

p p n --≠≠≠

=

≥∑∑∑ ()n ii f 表示系统从i 出发,

经n 步首次返回i 的概率. 下面的定理是联系转移概率()

n ij p 和首达概率()n ij f 之间关系的常用定理.

定理6 对任意状态 , , 1i j S n ∈≥,有 ()

()()

1.n

n m n m i j

i j j j m p

f p -==∑ 证明 {}{}

()

00,n ij n j n p P X j X i P T n X j X i ====≤==

{}01,n j n m P T m X j X i =??

====????

{}01, n j n m P T m X j X i =??

====????

{}01, n

j

n m P T

m X j X i =====∑

{}{}

001 , n

j

n j m P T

m X i P X j X i T m =======∑

{}()0111

, ,,,n

m ij

n m m m f

P X j X i X j X j X j -====≠≠=∑

{}{}()()011

n

n

m m ij

n m ij n m m m f

P X j X j f P X j X j -========∑∑

()()1

n

m n m ij jj m f

p -==∑

{}{}()001

1

n ij ij

j j n n f f

P T n X i P T X i ∞

=======<∞=∑∑

则ij f 表示系统从状态i 出发,经有限次到达状态j 的概率(也可以说成系统从状态i 出发,迟早要到达状态j 的概率). 显然有

()()

01n n ij ij ij f p f ≤≤≤≤

特别,ii f 表示从状态i 出发,迟早要返回状态i 的概率.

定义10 若1ii f =,则称i 是常返状态;若1ii f <,则称i 是非常返状态或滑过状态. (显然,吸收状态是常返状态)

对于常返状态i ,因{}

01ii i f P T X i =<∞==,这表明,从状态i 出发必定要返回它自身.

记{}{}:, 1, :, 1R ii N ii S i i S f S i i S f =∈==∈<,则,. R N R N S S S S S φ== 例10 设马氏链{}, 0n X n ≥的状态空间为{}1,2,3,4S =,转移矩阵为

140141201212001000001??

? ?= ? ? ???

P 试判别状态的常返性.

解 (1)()

111114,0

(2)n f f n ==≥ 1114.f = (1)(2)()

22222212,12, 0 (3)n f f f n ===≥ 22 1.f = (1)(2)(3)2

3333330,12, (12), f f f === 33 1.f =

(1)()

44441,0(2)n f f n ==≥ 44 1.f =

故 {}{}2,3,4,1.R N S S == 记00(),

()ij j i i E T X i E T X i μμ====,则ij μ表示系统从状态i 出发,首达状态

j 的平均时间;i μ表示系统从状态i 出发,首次返回i 的平均时间. 记j N 表示系统经过

j N 的次数,j N 是随机变量,()j n I X 表示系统在时刻n 经过状态j 的次数,即

1, () (1)0, n j n n X j

I X n X j

=?=≥?≠?

则1

().j j

n n N I

X ∞

==

记00(), ()ij j i i m E N X i m E N X i ====,则ij m 表示系统从状态i 出发,经过状态j 的平均次数;i m 表示系统从状态i 出发,返回状态i 的平均次数. 定理7 设, i j S ∈,则 ()()

1

1

, .n n ij ij i ii n n m p

m p ∞

===

=∑∑

证明 001

()(

())ij j j

n n m E N X i E I

X X i ∞

=====∑

01

(())

j

n n E I

X X i ∞

===∑ ()

01

1()n n

ij

n n P X

j X i p ∞

===

===∑∑

()

01

()n i i ii

n m E N X i p ∞

====∑ 定理8 设 , i j S ∈,则 {}0,

,0,.i j R j N f j S P N X i j S ∈?=∞==?

∈?

证明 {}{}{}00

1l i m j j j n n P N X i P N n X i P N n X

i

∞→∞

=??

=∞==≥==≥=???? 而

{}{}001, j j j m P N n X i P T m N n X i ∞=??

≥===≥=????

{}01, j

j m P T

m N n X i ∞

===≥=∑

{}{}

001, j

j j m P T

m X i P N n X i T m ∞

====≥==∑

{}{}()

001

11m ij

j ij j m f

P N n X j f P N n X j ∞

==

≥-==≥-=∑

{}{}

1

002()0n ij jj j ij jj j f f P N n X j f f P N X j -=≥-==≥=

1

()n ij jj f f -=

{}1

0,lim ()0, ij R

n j ij jj n N

f j S P N X i f f j S -→∞

∈?=∞===?

∈? 推论 设 i S ∈,则 {}01,,

0, .R i N i S P N X i i S ∈?=∞==?∈?

该推论说明了定义6的合理性,即:如果状态i 是常返的,则以概率1,系统无穷次返回状态i ;如果状态i 是非常返的,则以概率1,系统只有有限次返回状态i ,亦即系统无穷次返回状态i 的概率为0. 定理9 设 , i j S ∈,则

00, , 0

(),, 0(1), R ij ij j R ij ij jj N

j S f m E N X i j S f f f j S ?∈=?

===∞∈>??

-∈?

证明 (1)若 , 0R ij j S f ∈=,则 1n ?≥,有()

0n ij f =,由定理6知, 1n ?≥,

有()

0n ij p =,所以 0()0ij j m E N X i ===.

(2)若 , 0R i j j S f ∈>,则由定理8知,{}

00j P N X i =∞=>,所以ij m =

0()j E N X i β==∞.

(3)若 N j S ∈,则由定理8知,{}00j P N X i =∞==,所以

{}0

1

()ij

j

j

n m E N X i nP N n X i ∞

======∑

{}{}00

11j j n n P N n X i P N n X i ∞

=??=≥=-≥+=??

11

1

1

()()(1)()

n n

n ij jj ij jj ij jj

jj

n n n f f f f f f n f ∞

--==??=-=-??∑∑

1(1)()(1)n ij jj jj ij jj n f f f f f ∞

='?

?=-=-??

推论 设 i S ∈,则 0,

,()(1), .

R i i ii ii N i S m E N X i f f i S ∞∈?===?-∈?

定义11 设 , i j S ∈,若 1n ?≥,使()

0n ij p >,

则称状态i 可达状态j ,记为i j →; 若 1n ?≥,有()

0n ij p =,则称状态i 不可达状态j ,记为i j →/;若i j →且j i →,则

称状态i 和状态j 相通(或互通),记为i j ?.

定理10 (1)若i k →且k j →,则i j →;(2)若i k ?且k j ?,则i j ?.

证明 因i k →且k j →,所以 1, 1m n ?≥≥,使得()()0, 0m n ik kj p p >>,由C-K

方程,有()

()()()()

0m n m n m n ij

il lj ik kj l S

p p p p p +∈=≥>∑,故i j →. (2)由(1)易知(2)成立.

定理11 设 , i j S ∈,则(1)0ij i j f →?>;(2)0ij ji i j f f ??>.

证明 只证明(1). ""? 设i j →,则 1n ?≥,使()

0n ij p >,而

()()()1n

n m n m ij ij jj m p f p -==∑

因此,(1)(2)()

, ,,n ij ij ij f f f 中至少有一个大于0,从而()1

0m ij ij m f f ∞

==

>∑

. ""? 设0ij f >,因()1

m ij ij m f f ∞

==

,所以 1n ?≥,使()

0n ij f >,故 ()()()()(0)()

10n

n m n m n n ij ij jj ij jj ij

m p f p f p f -==≥=>∑ 即 i j →. 定理12 (1)()

1

n R ii n i S p ∞

=∈??=∞∑

; (2)()

1

n N ii n i S p ∞

=∈??

<∞∑

,这时必有()lim 0n ii n p →∞

=; (3)若, , , R i j S i S i j ∈∈→,则R j S ∈,且1ii ij ji jj f f f f ====. 证明 由定理7及定理9的推论即知(1)和(2);(3)略. 对于常返状态,即R i S ∈,由于()1

1n ii ii n f f ∞

==

=∑

,因此{}()

, 1n ii f n ≥是一个概率

分布,故

{}()0011()n i i i ii n n E T X i nP T n X i nf μ∞

========∑∑

易知 1.i μ≥

定义12 设状态i 是常返的,即R i S ∈. 若i μ<∞,则称i 是正常返的;若i μ=∞,

则称i 是零常返的. 记 {}{}0

:, , :, R R i R R i S i i S S i i S μμ+=∈<∞=∈=∞.

定义13 (1)设C S ?,若 , i C j S C ?∈∈-,有i j →/,则称C 是闭集;(2)设C S ?,若 , i j C ?∈,都有i j →,则称C 是不可约的;(3)若S 是不可约的,则称马氏链{}, 0n X n ≥是不可约的(即 , i j S i j ?∈?→). 定理13 设C S ?,则下列各条等价 (1)C 是闭集;

(2)() , ,10n ij i C j C n p ?∈?≥?=;

(3) , 0ij i C j C f ?∈??=.

注 (1)若C 是闭集,则 , 1i C n ?∈≥,有

()

1n ij j C

p ∈=∑

.

(2)整个状态空间S 是闭集,是最大的闭集;吸收状态i (即1ii p =)是闭集,是最小的闭集.

定理14 0

,, R R R

S S S +都是闭集. 证明 设, .R R i S j S ∈? 若0ij f >,则i j →,故 R j S ∈. 矛盾,从而0ij f =,由上面的定理13,知R S 是闭集. 同理可证0

, R R S S +

也是闭集.

例11 设{}, 0n X n ≥是马尔可夫链,状态空间为{}1,2,3,4,5S =,转移矩阵为

12

00

12012012 0 00

01 0 0100 0 001

0 0 0??

?

? ?= ?

? ??

?

P

问状态空间{}1,2,3,4,5S =中是否含有真的不可约的闭子集?{}, 0n X n ≥是否是不可约马尔可夫链.

解 易知,状态3是吸收状态,故{}3是闭集;{}{}{}1,4, 1,3,4, 1,2,3,4均是闭集;其中{}{}3, 1,4均是不可约的;因为S 含有真的闭子集,所以{}, 0n X n ≥为非不可约马尔可夫链.

例12 设马尔可夫链{}, 0n X n ≥的状态空间为{}1,2,,9S = ,状态之间的转移概率如图所示. 由图可见:自状态1出发,再返回状态1的可能步数为

T =

{}4, 6, 8, 10,

显然T 的最大公约数为2. 但2T ?,即自状态1出发经过2步不能返回状态1. 但我们仍把2定义为状态2的周期.

定义14 设马尔可夫链{}, 0n X n ≥的状态空间为S ,n 步转移概率为()

n ij p ,对于

i S ∈,若集合{}()1, 0n ii n n p ≥>非空,则称该集合的最大公约数

{}()

()G.C.D 1, 0n ii d d i n n p =≥>

为状态i 的周期. 如果1d >,称状态i 为周期的,如果1d =,称状态i 为非周期的. 若

{}()

:1, 0n ii

n n p

≥>是空集,则不对状态i 定义周期.

由定义知,如果状态i 有周期d ,则对一切[]0mod()n d ≠,都有()

0n ii p =. 但这并不是说对任意的nd ,都有()0nd ii p >. 如在例12中,

状态1的周期2d =,但(2)110p =. 然而有如下结论.

定理15 设状态i 的周期为d ,则存在正整数M ,对一切n M ≥,有()0nd ii p >.

证明 略(证明用到了初等数论的知识).

定理16 设状态i S ∈,如果集合{}

()

1, 0n ii n n f ≥>非空,则

{}{}()

()G.C.D 1, 0G.C.D 1, 0n n ii ii n n p n n f ≥>=≥>

证明 记 {}()G.C.D 1, 0n ii d n n p =≥>,{}

()

G.C.D 1, 0n ii c n n f =≥>,由定理6

()()()()(0)

()1

n

n m n m n n ii

ii ii ii ii ii m p

f p f p f -==≥=∑ 从而{}{

}()()

1, 01, 0

n n ii ii n n p

n n f ≥>?≥>,于是有 1d c ≤≤. 如果1c =,则1d c ==. 如果1c >,只需证明d c ≥,为此只需证明c 是{}()1, 0n ii

n n p ≥>的公约数

即可. 换言之,只需证明对于一切[]0mod()n c ≠,都有()

0n ii p =.

实际上,由c 的定义知,当1r c ≤<时,有()0r ii f =. 于是,对1n c ≤<,有

()()()()

1

1

00n

n

n m n m n m ii

ii

ii

ii m m p

f

p

p --====?=∑∑

现假设当, 0,1,2,,1n kc r k N =+=- 时,有()

0n ii p =. 注意到[]0mod()n c ≠时,

有()0n ii f =. 则由定理6得

()

()()

1

Nc r

Nc r m Nc r m ii

ii ii

m p

f

p +++-==

∑ ()()(2)(2)()()

0c Nc r c c Nc r c Nc r ii ii ii ii ii ii f p f p f p +-+-=+++=

于是,由归纳法知,当[]0mod()n c ≠时,有()

0n ii p =.

定理17 设状态i 常返且有周期d ,则 ()

lim nd ii

i n p d μ→∞

=.(其中i μ为状态i 的平均

返回时间,当i μ=∞时,约定0i d μ=)

证明 见可数状态的马尔可夫过程,胡迪鹤,武汉大学出版社,1983,2526P -. 定义15 非周期的正常返状态称为遍历状态. 定理18 设i 是常返状态,则

① i 是零常返状态??()

lim 0n ii n p →∞

=;

② i 是遍历状态??()

lim 10n ii i n p μ→∞

=>.

证明 ① 设i 是零常返状态,则i μ=∞,由定理17,得 ()

lim 0nd ii

n p →∞

=,而当

[]0mod()n d ≠时,有()0n ii p =,于是得()lim 0n ii n p →∞

=. 反之,设()lim 0n ii n p →∞

=,假设i 是

正常返状态,则i μ<∞,于是由定理17,得()

lim 0nd ii

n p →∞

>,矛盾.

② 设()

lim 10n ii i n p μ→∞

=>,由①知,i 不能是零常返状态,从而只能是正常返状态,

由定理17,得()

lim nd ii

i n p d μ→∞

=,注意到()()

lim lim nd n ii ii n n p p →∞

→∞

=,得1d =,所以,i 是遍历

状态. 反之,由定理17是显然的.

定理19 设状态, i j S ∈,且i j ?,则

(1)i 与j 同为常返状态或同为非常返状态;如果同为常返状态,则它们同为正常返状态或同为零常返状态. (2)i 与j 有相同的周期.

证明 (1)由于i j ?,由可达的定义知存在1l ≥和1n ≥,使得

()

()0, 0l n ij ji p p αβ=>=>

由C-K 方程,对任意的1m ≥,总有

()()()()()l m n l m n m ii ij jj ji jj p p p p p αβ++≥=,()()()()()

l m n n m l m jj ji ii ij ii

p p p p p αβ++≥= 将上式两边从1到∞求和,得

()()

11l m n m ii jj

m m p p αβ∞

++==≥∑

∑, ()()

11

l m n m jj ii m m p p αβ∞

++==≥∑

∑ 可见,

()1k ii k p ∞

=∑与()1k jj k p ∞

=∑相互控制,所以它们同为无穷或同为有限. 由定理12知,

i 与j 同为常返状态或同为非常返状态.

若i 与j 同为常返状态,在以上的不等式中取极限,令m →∞,得

()

()lim lim l m n m ii jj m m p p αβ++→∞→∞

≥, ()()

l i m l i m l m n m j j i i m m p p

αβ++→∞

→∞

因此,()

lim k ii k p →∞

与()

lim k jj k p →∞

同为正或同为零,由定理18知,i 与j 同为正常返状态或同为

零常返状态.

(2)设状态i 的周期为d ,状态j 的周期为c . 因为 ()

()0, 0l n ij ji p p αβ=>=>,所以对任一使()0m jj p >的m ,必有()()()()()0l m n l m n m ii ij jj ji jj p p p p p αβ++≥=>,从而d 可除尽l m n ++. 又因为()()()0l n l n ii ij ji p p p αβ+≥=>,

所以d 也可除尽l n +. 因此,d 可除尽m . 这说明d c ≤. 同理可证d c ≥. 故d c =,即状态i 与j 有相同的周期. 注 状态的常返性与马尔可夫链的初始分布是无关的

定理20 状态空间S 可唯一地分解成有限个或可列无穷个互不相交的子集之和,即

(1)(2)N R R S S S S = ,且使得

(1)每个()k R S 是常返状态组成的不可约闭集;

(2)()k R S 中的状态或全是正常返状态,或全是零常返状态,且有相同的周期; (3)N S 是由全体非常返状态组成,自()k R S 中的状态不能到达N S 中的状态.

定理21 周期为d 的不可约马尔可夫链,其状态空间S 可唯一地分解为d 个互不相交的子集之和,即0121d S S S S S -= ,且使得自r S 中任一状态出发,经一步转移必进入1r S +中,0,1,2,,1r d =- . (约定:0d S S =) 证明 (1)任意取定一状态i S ∈,令

{}()

0, 0nd r r ij S j n p +=?≥>使,0,1,2,,1r d =-

因S 是不可约的,即S 中的状态是互通的,故

10

d r r S S -==

.

(2)若r t j S S ?∈ (r t ≠),由定义知,必存在非负整数n 和m 使得()

0nd r ij

p +>,

()

0md t ij p +>. 又由于j i ?,必存在正整数h ,使()0h ji p >,于是由C-K 方程,有

()()()0nd r h nd r h ii ij ji p p p +++≥>, ()()()0md t h md t h ii ij ji p p p +++≥>

所以,d 能除尽nd r h ++,又能除尽md t h ++,从而d 能除尽

()()()nd r h md t h n m d r t ++-++=-+- 故d 能除尽r t -,注意到 01, 01r d t d ≤≤-≤≤-,故只能0r t -=,这说明,当r t ≠时,r t S S =? .

(3)下面证明对任一r j S ∈,有

1

1r jk k S p +∈=∑

. 事实上

1

1

1

1r r r jk jk jk jk k S k S k S k S p p p p +++∈∈?∈==+=∑∑∑∑

最后一个等式是因:r j S ∈,由定义有()

0nd r ij

p +>,而当1r k S +?时,由定义有

(1)

0nd r ik p ++=,由C-K 方程,有(1)()0nd r nd r ik ij jk p p p +++=≥,故0jk p =.

(4)最后证明分解的唯一性,这只需证明{}0121,,,,d S S S S - 不依赖于最初取定的状态i . 亦即只需证明:对取定的状态i ,如果, r j k S ∈,则对另取的状态i ',仍有

, r j k S '∈(r '可以与r 不同).

设对i 的分解为{}0121,,,,d S S S S - ,对i '的分解为{}0

121,,,,d S S S S -'''' ,又设, r j k S ∈,t i S '∈. 于是,当r t ≥时,自i '出发,以概率1只能在, ,r t r t d --+

2,r t d -+ 等这些步上转移到j 或k ,从而, r t j k S -'∈;当r t <时,自i '出发,以概

率1只能在(), 2, 3,d t r r t d r t d r t d --=-+-+-+ 等这些步上转移到j 或k ,从

而, r

t d j k S -+'∈. 定理22 设{}, 0n X n ≥是周期为d 的不可约的马尔可夫链,则得一新的马尔可夫

链{}, 0nd X n ≥,其一步转移矩阵为()()

()d d ij p =P ,原马尔可夫链{}, 0n X n ≥按照定理

21分解出的每个r S ,是新马尔可夫链{}, 0nd X n ≥的不可约的闭集,且r S 中的状态对新马尔可夫链是非周期的. (证明略)

例13 设不可约马尔可夫链{}, 0n X n ≥的状态空间为{}1,2,3,4,5,6S =,周期为

3d =,转移矩阵为

012012

01001301010000 0 010 0 0 0 100 0 0 0 0

140 340?? ?

? ?

=

? ? ?

?

??

?

P

对取定的状态1,易知

{}{}(3)

010, 01,4,6n j S j n p ?≥>= 使 {}{}(31)110, 03,5n j

S j n p +?≥>= 使 {}{}(32)210, 02n j

S j n p +?≥>= 使 故{}{}{}0121,4,63,52S S S S == . 马尔可夫链{}3, 0n X n ≥的转移矩阵为

(3)

13

0010

13010000007120512013001013007120512013001013?? ?

? ?

= ? ? ?

?

???

P

随机过程 第五章 连续时间的马尔可夫链

第五章 连续时间的马尔可夫链 5.1连续时间的马尔可夫链 考虑取非负整数值的连续时间随机过程}.0),({≥t t X 定义5.1 设随机过程}.0),({≥t t X ,状态空间}0,{≥=n i I n ,若对任意 121...0+<<<≤n t t t 及I i i i n ∈+121,...,,有 })(,...)(,)()({221111n n n n i t X i t X i t X i t X P ====++ =})()({11n n n n i t X i t X P ==++ (5.1) 则称}.0),({≥t t X 为连续时间马尔可夫链. 由定义知,连续时间马尔可夫链是具有马尔可夫性的随机过程,即过程在已知现在时刻n t 及一切过去时刻所处状态的条件下,将来时刻1+n t 的状态只依赖于现在状态而与过去无关. 记(5.1)式条件概率一般形式为 ),(})()({t s p i s X j t s X P ij ===+ (5.2) 它表示系统在s 时刻处于状态i,经过时间t 后转移到状态j 的转移概率. 定义5.2 若(5.2)式的转移概率与s 无关,则称连续时间马尔可夫链具有平稳的或齐次的转移概率,此时转移概率简记为 ),(),(t p t s p ij ij = 其转移概率矩阵简记为).0,,()),(()(≥∈=t I j i t p t P ij 以下的讨论均假定我们所考虑的连续时间马尔可夫链都具有齐次转移概率.简称为齐次马尔可夫过程. 假设在某时刻,比如说时刻0,马尔可夫链进入状态i,而且接下来的s 个单位时间单位中过程未离开状态i,(即未发生转移),问随后的t 个单位时间中过程仍不离开状态i 的概率是多少呢?由马尔可夫我们知道,过程在时刻s 处于状态i 条件下,在区间[s,s+t]中仍然处于i 的概率正是它处于i 至少t 个单位的无条件概率..若记 i h 为记过程在转移到另一个状态之前停留在状态i 的时间,则对一切s,t 0≥有 },{}{t h P s h t s h P i i i >=>+> 可见,随机变量i h 具有无记忆性,因此i h 服从指数分布. 由此可见,一个连续时间马尔可夫链,每当它进入状态i,具有如下性质: (1) 在转移到另一状态之前处于状态i 的时间服从参数为i v 的指数分布;

马尔可夫链蒙特卡罗在实践中的应用

2012年第12期 吉林省教育学院学报 No.12,2012 第28卷JOURNAL OF EDUCATIONAL INSTITUTE OF JILIN PROVINCE Vol .28(总300期) Total No .300 收稿日期:2012—11—14 作者简介:孟庆一(1989—),女,吉林长春人,新加坡籍华人,英国伦敦大学数学系,本科生,研究方向:MCMC 统计学。 浅议马尔可夫链蒙特卡罗在实践中的应用 孟庆一 (英国伦敦大学,英国伦敦) 摘要:本文概括地介绍了马尔可夫链蒙特卡罗(Markov chain Monte Carlo ———MCMC ),一种随机模拟贝叶斯推断的方法。主要的抽样方法包括吉布斯采样(Gibbs Sampling )和Metropolis -Hastings 算法。本文也对MCMC 主题和应用的拓展进行了讨论。 关键词:马尔可夫链;蒙特卡罗;Gibbs 抽样;Metropolis -Hastings 中图分类号:O29 文献标识码:A 文章编号:1671—1580(2012)12—0120—02 统计学中的贝叶斯推理在过去的几十年里有前 所未有的突破,统计学家们发现了一种非常简单,但又非常强大的模拟技术,统称为MCMC 。这种技术可以运用到各种复杂的贝叶斯范例和实际情况。 贝叶斯推理: 贝叶斯方法把所给的模型里所有的未知量的不确定性联系在一起。利用所知的信息,贝叶斯方法用联合概率分布把所有未观察到的数量综合起来,从而得出的推论。在这里,给定已知的未知分布被称为后验分布。有关未知量的推理被称为预测,它们的边缘分布称作为预测分布。 贝叶斯推理根据贝叶斯规则计算后验概率: P (H |E )= P (E |H )·P (H ) P (E )然而,在大多数情况下,所给的模型的复杂性不允许我们运用这个简单的操作。因此,我们需要使用随机模拟, 或蒙地卡罗技术来代替。概述MCMC : MCMC 采用未知量的高维分布,为难度极高的模拟复杂模型的问题提供了一个答案。 一个马尔可夫链是一个序列的随机变量X 1,X 2,X 3,...这个序列有马尔可夫的属性———给予目前的状态,未来和过去的状态是独立的。从数学公 式上看, Pr (X n +1=x |X 1=x 1,X 2=x 2,…,X n =x n )=Pr (X n +1=x |X n =x n )X i 的可能的值可数的集合S 称 为链的状态空间。 幸运的是,在马尔可夫链里,我们也有与大数定律和中心极限定理类似的定理。 另外一个问题存在于如何建立一个马尔可夫链的极限分布与所需的分配一模一样。一种可行的解决方案是Gibbs 抽样。它是基于一个马尔可夫链,其前身的依赖性是由模型中出现的条件分布所决定的。另一种可能性是Metropolis -Hastings 算法。它是基于一个马尔可夫链,其前身的依赖性是分裂成两个部分:一个是建议,另一个是接受这一建议。 Metropolis -Hastings 算法: Metropolis -Hastings 算法,可以从任何概率分布中抽取样品,只要求是可计算函数的密度成正比。在贝叶斯的应用程序中,归一化因子计算往往是非常困难的,所以,和其他常用的抽样算法一样,能够在不知道这个比例常数的情况下产生样本是Metropolis -Hastings 算法的重要特征。 该算法的总体思路是产生一系列在一个马尔可 夫链里的样品。在足够长的时间后,所生成的样品的分布与分布相匹配。 该算法基本上按如下方式工作(这是一个特殊 的例子,其建议密度是对称的情况下):首先,选择一个任意的概率密度Q (x'|x t ),这表明一个新的采样值x'给定样本值x t 。对于简单的Metropolis 算法,这个建议密度必须是对称的Q (x'| 21

第五章 连续时间的Markov链

第五章 连续时间的马尔可夫链 第四章我们讨论了时间和状态都是离散的M arkov 链,本章我们研究的是时间连续、状态离散的M arkov 过程,即连续时间的M arkov 链. 连续时间的M arkov 链可以理解为一个做如下运动的随机过程:它以一个离散时间M arkov 链的方式从一个状态转移到另一状态,在两次转移之间以指数分布在前一状态停留. 这个指数分布只与过程现在的状态有关,与过去的状态无关(具有无记忆性),但与将来转移到的状态独立. 5.1 连续时间马尔可夫链的基本概念 定义 5.1 设随机过程{(),0}X t t ≥,状态空间{,1}n I i n =≥,若对任意的正整数 1210n t t t +≤<<< 及任意的非负整数121,,,n i i i I +∈ ,条件概率满足 {}111122()|(),(),,()n n n n P X t i X t i X t i X t i ++==== {}11()|()n n n n P X t i X t i ++=== (5.1) 则称{(),0}X t t ≥为连续时间的M arkov 链. 由定义知,连续时间的M arkov 链是具有M arkov 性(或称无后效性)的随机过程,它的直观意义是:过程在已知现在时刻n t 及一切过去时刻所处状态的条件下,将来时刻1n t +的状态只依赖于现在的状态而与过去的状态无关. 记(5.1)式条件概率的一般形式为 {()|()}(,)ij P X s t j X s i p s t +=== (5.2) 它表示系统在s 时刻处于状态i ,经过时间t 后在时刻s t +转移到状态j 的转移概率,通常称它为转移概率函数.一般地,它不仅与t 有关,还与s 有关. 定义 5.2 若(5.2)式的转移概率函数与s 无关,则称连续时间M arkov 链具有平稳的转移概率函数,称该M arkov 链为连续时间的齐次(或时齐)M arkov 链. 此时转移概率函数简记为(,)()ij ij p s t p t =.相应地,转移概率矩阵简记为()(()),(,,0)ij P t p t i j I t =∈≥. 若状态空间{0,1,2,}I = ,则有 ()00010210 11 12 012() ()() ...()()()()()... ... .. ....()()( )...... .. .... ij n n n p t p t p t p t p t p t P t p t p t p t p t ?? ? ? ?== ? ? ?? ? (5.3) 假设在某时刻,比如说时刻0,M arkov 链进入状态i ,在接下来的s 个单位时间内过程 未离开状态i (即未发生转移),我们要讨论的问题是在随后的t 个单位时间中过程仍不离开状态i 的概率是多少?由M arkov 性知,过程在时刻s 处于状态i 的条件下,在区间[,] s s t +

马尔可夫链

马尔可夫链 马尔可夫链(Markov chains )是一类重要的随机过程,它的状态空间是有限的或可数无限的。经过一段时间系统从一个状态转到另一个状态这种进程只依赖于当前出发时的状态而与以前的历史无关。马尔可夫链有着广泛的应用,也是研究排队系统的重要工具。 1) 离散时间参数的马尔可夫链 ①基本概念 定义 5.7 设{()0,1,2,}X n n ???=,是一个随机过程,状态空间{0,1,2,}E =,如果对于任意的一组整数 时间120k n n n ???≤<<<,以及任意状态12,, ,k i i i E ∈,都有条件概率 11{()|()}k k k k P X n i X n i --=== (5-17) 即过程{()0,1,2,}X n n ???=,未来所处的状态只与当前的状态有关,而与以前曾处于什么状态无关,则称 {()0,1,2,}X n n ???=,是一个离散时间参数的马尔可夫链。当E 为可列无限集时称其为可列无限状态的马尔可 夫链,否则称其为有限状态的马尔可夫链。 定义5.8 设{()0,1,2,}X n n ???=,是状态空间{0,1,2, }E =上的马尔可夫链,条件概率 (,){()|()}ij p m k P X m k j X m i i j E =+==∈,、 (5-18) 称为马尔可夫链{()0,1,2,}X n n ???=,在m 时刻的k 步转移概率。 k 步转移概率的直观意义是:质点在时刻m 处于状态i 的条件下,再经过k 步(k 个单位时间)转移到状 态j 的条件概率。特别地,当1k =时, (,1){(1)|()}ij p m P X m j X m i =+== (5-19) 称为一步转移概率,简称转移概率。 如果k 步转移概率(,)ij p m k i j E ∈,、,只与k 有关,而与时间起点m 无关,则{()}X n 称为离散时间的齐次马尔可夫链。 定义5.9 设{()0,1,2,}X n n ???=,是状态空间{0,1,2,}E ???=上的马尔可夫链,矩阵 0001010 11101(,)(,)(,)(,)(,)(,)(,)(,)(,) (,) n n j j jn p m k p m k p m k p m k p m k p m k P m k p m k p m k p m k ?? ???? ? ?=? ?????? ? (5-20) 称为{()}X n 在m 时刻的k 步转移概率矩阵。 当1k =时,(,1)P m 称为一步转移概率矩阵。 对于齐次马尔可夫链,容易推得k 步转移概率矩阵与一步转移概率矩阵具有关系 ()(),,1k P m k P m =????,1,2,k ???= (5-21)

课上练习题_离散时间马尔科夫链 423

1、4.23 Trials are performed in sequence. If the last two trials were successes, then the next trial is a success with probability 0.8; otherwise the next trial is a success with probability 0.5. In the long run, what proportion of trials are successes? 2、4.32 Each of two switches is either on or off during a day. On day n, each switch will independently be on with probability [1+#of on switches during day n-1]/4. For instance, if both switches are on during day n-1, then each will independently be on during day n with probability3/4. What fraction of days are both switches on? What fractions are both off?

3、Let ri denote the long-run proportion of time a given irreducible Markov chain is in state i. Explain why ri is also the proportion of transitions that are into state i as well as being the proportion of transition that are from state i. 4、4.44 Suppose that a population consists of a fixed number, say, m, of genes in any generation. Each gene is one of two possible genetic types. If any generation has exactly i (of its m) genes being type 1, then the next generation will have j type 1 genes with probability j m j m i m m i j m- ? ? ? ? ?- ? ? ? ? ? ?? ? ? ? ? . Let Xn denote the number of type 1 genes in the nth generation, and assume that X0 = i. (a) Find E[Xn] (b) What is the probability that eventually all the genes will be type 1?

马尔可夫链预测方法及其一类应用【开题报告】

开题报告 数学与应用数学 马尔可夫链预测方法及其一类应用 一、综述本课题国内外研究动态, 说明选题的依据和意义 概率论自1654年创立以来, 已由最初的博弈分析问题发展成为现今的方法论综合性学科. 而其中随机过程已经是现代概率论发展的必然性. 在这其中, 马尔可夫在1906年的"大数定理关于相依变量的扩展"(Extension de la loi de grands bombers etc)论文中首次创立的马尔可夫链已经成为了概率论的重中之重. 马尔可夫是世界上著名的数学家、社会学家. 他所研究的范围非常的广泛, 涉及到概率论、数论、数的集合、函数逼近论、数理统计、微分方程等方面. 马尔可夫在1906~1912年间, 他提出并研究了一种能用数学分析方法研究自然过程的一般图示, 后人把这种图示以他的姓氏命名为马尔可夫链(Markov Chain). 在当时, 马尔可夫开创性地采用了一种对无后效性的随机过程的研究范式, 即在已知当前状态的情况下, 过程的未来状态与其过去状态无关, 这就是现在大家非常熟悉了解的马尔可夫过程. 在现实生活当中, 有许多过程都能被看作成马尔可夫过程. 如软件可靠性测试、传染病受感染的人数、农村剩余劳动力流动趋势预测、液体中微粒所作的布朗运动、产品市场占有率及利润率的变动等等. 也正是由于马尔可夫链在生活中所具有的普遍存在性, 马尔可夫链理论才被广泛应用于近代的物理学, 生物学, 地质学, 计算机科学, 公共事业, 教育管理、经济管理、以及企业人员管理、桥梁建筑等各个领域. 马尔可夫链运用数学模型对定性问题进行预测提供了一种思路, 丰富了预测的内容. 其大体上可以分为以下几个步骤: 首先, 把现象看作成为一个系统, 并对该系统进行科学的划分. 根据系统的实际和需要划分出多个状态, 系统所划分出来的各个状态就是要预测的内容. 其次, 对现象各种状态的状态概率进行统计测定, 也就是判定出系统当前处于什么状态. 然后, 对各系统未来发展的每次转移概率进行预测, 就是要确定出系统是如何转移的. 最后, 根据系统当前的各种状态和转移概率矩阵, 推测出系统经过若干次转移后, 到达

第章离散时间的马尔可夫链

第1章 离散时间的马尔可夫链 §1 随机过程的基本概念 定义1 设(,,)P ΩF 是概率空间,(, )E E 是可测空间, T 是指标集. 若对任何t T ∈,有 :t X E Ω→,且t X ∈F E ,则称{}(), t X t T ω∈是(, , )P ΩF 上的取值于(,)E E 中的随机过 程,在无混淆的情况下简称{(), }t X t T ω∈为随机过程,称(,)E E 为状态空间或相空间,称E 中的 元素为状态,称T 为时间域. 对每个固定的ω∈Ω,称()t X ω为 {}(), t X t T ω∈对应于ω的轨道或现 实,对每个固定的t T ∈,称()t X ω为E 值随机元. 有时()t X ω也记为 设 T ?R ,{}, t t T ∈F 是F 中的一族单调增的子σ代数(σ代数流),即 ① t t T ?∈??F F ,且t F 是σ代数; ② , , s t s t T s t ?∈

马尔可夫链预测方法及其一类应用【文献综述】

文献综述 数学与应用数学 马尔可夫链预测方法及其一类应用 马尔可夫性是俄国数学家A.A.Mapkov 在1906年最早提出的. 但是, 什么是马尔可夫性呢? 一般来讲,认为它是“相互独立性”的一种自然推广. 设有一串随机事件,...,,...,,121n n A A A A -中(即n A 属于概率空间(P ,,ξΩ)中的σ代数ξ,1≥n ), 如果它们中一个或几个的发生, 对其他事件的发生与否没有影响, 则称这一串事件是相互独立的(用概率空间(P ,,ξΩ)的符号表示, 即))()(11n m n m n n A P A P X I ===, 推广下, 如果在已知,...,1+n n A A 中的某些事件的发生, 与,,...,,121-n A A A 中的事件发生与否无关, 则称这一串事件{1:≥n A n }具有马尔可夫性. 所以说, 马尔可夫性可视为相互独立性的一种自然推广. 从朴素的马尔可夫性, 到抽象出马尔可夫过程的概念, 从最简单的马尔可夫过程到一般的马尔可夫过程, 经历了几十年的发展过程. 它有极其深厚的理论基础, 如拓扑学、函数论、几何学、近世代数、泛函分析. 又有很广泛的应用空间, 如随机分形、近代物理、公共事业中的服务系统、电子信息、计算技术等. 在现实世界中, 有很多过程都是马尔可夫过程, 如软件可靠性测试、传染病受感染的人数、农村剩余劳动力流动趋势预测、液体中微粒所作的布朗运动、产品市场占有率及利润率的变动, 车站排队问题等等, 都可视为马尔可夫过程. 所谓马尔可夫链是指时间连续(或离散)、状态可列、时间齐次的马尔可夫过程. 之所以要研究这种过程, 一方面是由于它的理论比较完整深入, 可以作为一般马尔可夫过程及其他随机过程的借鉴; 二是由于它在自然科学和许多实际问题(如遗传学、教育学、经济学、建筑学、规则论、排队论等)中发挥着越来越大的作用. 自从我国著名数学家、教育家、中科院王梓坤院士在上世纪50年代将马尔可夫理论引入国内以后, 我国数学家对马尔可夫过程的研究也取得了非常好的效果, 在生灭过程的构造和它的积分型泛函的分布、马尔可夫过程的零壹律、Martin 边界与过份函数、马尔可夫过程

课上练习题_连续时间马尔科夫链 619

6.2 Suppose that a one-celled organism can be in one of two states-either A or B. An individual in state A will change to state B at an exponential rate α; an individual in state B divides into two new individuals of type A at an exponential rate β. Define an appropriate continuous-time Markov chain for a population of such organisms and determine the appropriate parameters for this model. 6.3 Consider two machines that are maintained by a single repairman. Machine i functions for an exponential time with rate μbefore breaking down, i = 1,2. The repair times (for either i machine) are exponential with rate μ. Can we analyze this as a birth and death process? If so, what are the parameters? If not, how can we analyze it?

马尔可夫性与马尔可夫链

马尔可夫性与马尔可夫链 【教学目标】 1.掌握马尔可夫性与马尔可夫链。 2.熟练运用马尔可夫性与马尔可夫链解决具体问题。 3.亲历马尔可夫性与马尔可夫链的探索过程,体验分析归纳得出马尔可夫性与马尔可夫链,进一步发展学生的探究、交流能力。 【教学重难点】 重点:掌握马尔可夫性与马尔可夫链。 难点:马尔可夫性与马尔可夫链的实际应用。 【教学过程】 一、直接引入 师:今天这节课我们主要学习马尔可夫性与马尔可夫链,这节课的主要内容有马尔可夫性与马尔可夫链,并且我们要掌握这些知识的具体应用,能熟练解决相关问题。 二、讲授新课 (1)教师引导学生在预习的基础上了解马尔可夫性与马尔可夫链内容,形成初步感知。 (2)首先,我们先来学习马尔可夫性,它的具体内容是: 1n X +的随机变化规律与0X ,1X ,…1n X -的取值都没有关系,随机变量序列{}n X 的所具有的这类性质称为马尔可夫性 它是如何在题目中应用的呢?我们通过一道例题来具体说明。 例: 马尔可夫性描述了一种_____。 解析:状态序列 可以给学生一定的提示。 根据例题的解题方法,让学生自己动手练习。 练习: 序列所有可能取值的集合,被称为_____。 (3)接着,我们再来看下马尔可夫链内容,它的具体内容是:

一般地,我们称具有马尔可夫性的随机变量序列{}n X为马尔可夫链。 它是如何在题目中应用的呢?我们也通过一道例题来具体说明。 例:请同学们查询资料,判断马尔可夫链与布朗运动是否有联系 解析:马尔可夫链与布朗运动以及遍历假说这两个二十世纪初期物理学重要课题是相联系的,但马尔可夫寻求的似乎不仅于数学动机,名义上是对于纵属事件大数法则的扩张。 根据例题的解题方法,让学生自己动手练习。 练习: 请写出马尔科夫链满足的两个假设。 三、课堂总结 (1)这节课我们主要讲了马尔可夫性与马尔可夫链 (2)它们在解题中具体怎么应用? 四、习题检测 1.请同学们写出马尔可夫性的定义。 2.请同学们写出马尔科夫链的定义。 3.请同学们写出马尔科夫性和马尔科夫链之间的联系。

马尔可夫链

马尔可夫过程 编辑词条 一类随机过程。它的原始模型马尔可夫链,由俄国数学家A.A.马尔可夫于1907年提出。该过程具有如下特性:在已知目前状态(现在)的条件下,它未来的演变(将来)不依赖于它以往的演变 ( 过去 ) 。例如森林中动物头数的变化构成——马尔可夫过程。在现实世界中,有很多过程都是马尔可夫过程,如液体中微粒所作的布朗运动、传染病受感染的人数、车站的候车人数等,都可视为马尔可夫过程。关于该过程的研究,1931年A.H.柯尔莫哥洛夫在《概率论的解析方法》一文中首先将微分方程等分析的方法用于这类过程,奠定了马尔可夫过程的理论基础。 目录 马尔可夫过程 离散时间马尔可夫链 连续时间马尔可夫链 生灭过程 一般马尔可夫过程 强马尔可夫过程 扩散过程 编辑本段马尔可夫过程 Markov process 1951年前后,伊藤清建立的随机微分方程的理论,为马尔可夫过程的研究开辟了新的道路。1954年前后,W.费勒将半群方法引入马尔可夫过程的研究。流形上的马尔可夫过程、马尔可夫向量场等都是正待深入研究的领域。 类重要的随机过程,它的原始模型马尔可夫链,由俄国数学家Α.Α.马尔可夫于1907年提出。人们在实际中常遇到具有下述特性的随机过程:在已知它目前的状态(现在)的条件下,它未来的演变(将来)不依赖于它以往的演变(过去)。这种已知“现在”的条件下,“将来”与“过去”独立的特性称为马尔可夫性,具有这种性质的随机过程叫做马尔可夫过程。荷花池中一只青蛙的跳跃是马尔可夫过程的一个形象化的例子。青蛙依照它瞬间或起的念头从一片荷叶上跳到另一片荷叶上,因为青蛙是没有记忆的,当现在所处的位置已知时,它下一步跳往何处和它以往走过的路径无关。如果将荷叶编号并用X0,X1,X2,…分别表示青蛙最初处的荷叶号码及第一次、第二次、……跳跃后所处的荷叶号码,那么{Xn,n≥0} 就是马尔可夫过程。液体中微粒所作的布朗运动,传染病受感染的人数,原子核中一自由电子在电子层中的跳跃,人口增长过程等等都可视为马尔可夫过程。还有些过程(例如某些遗

马尔可夫链模型讲解

马尔可夫链模型(Markov Chain Model) 目录 [隐藏] 1 马尔可夫链模型概述 2 马尔可夫链模型的性质 3 离散状态空间中的马尔可夫链模 型 4 马尔可夫链模型的应用 o 4.1 科学中的应用 o 4.2 人力资源中的应用 5 马尔可夫模型案例分析[1] o 5.1 马尔可夫模型的建立 o 5.2 马尔可夫模型的应用 6 参考文献 [编辑] 马尔可夫链模型概述 马尔可夫链因安德烈·马尔可夫(Andrey Markov,1856-1922)得名,是数学中具有马尔可夫性质的离散时间随机过程。该过程中,在给定当前知识或信息的情况下,过去(即当期以前的历史状态)对于预测将来(即当期以后的未来状态)是无关的。 时间和状态都是离散的马尔可夫过程称为马尔可夫链, 简记为 。 马尔可夫链是随机变量的一个数列。这些变量的范围,即他们所有可能取值的集合,被称为“状态空间”,而Xn的值则是在时间n的状态。如果Xn + 1对于过去状态的条件概率分布仅是Xn的一个函数,则 这里x为过程中的某个状态。上面这个恒等式可以被看作是马尔可夫性质。

马尔可夫在1906年首先做出了这类过程。而将此一般化到可数无限状态空间是由柯尔莫果洛夫在1936年给出的。 马尔可夫链与布朗运动以及遍历假说这两个二十世纪初期物理学重要课题是相联系的,但马尔可夫寻求的似乎不仅于数学动机,名义上是对于纵属事件大数法则的扩张。 马尔可夫链是满足下面两个假设的一种随机过程: 1、t+l时刻系统状态的概率分布只与t时刻的状态有关,与t时刻以前的状态无关; 2、从t时刻到t+l时刻的状态转移与t的值无关。一个马尔可夫链模型可表示为=(S,P,Q),其中各元的含义如下: 1)S是系统所有可能的状态所组成的非空的状态集,有时也称之为系统的状态空间,它可以是有限的、可列的集合或任意非空集。本文中假定S是可数集(即有限或可列)。用小写字母i,j(或S i,S j)等来表示状态。 2)是系统的状态转移概率矩阵,其中P ij表示系统在时刻t处于状态i,在下一时刻t+l处于状态i的概率,N是系统所有可能的状态 的个数。对于任意i∈s,有。 3)是系统的初始概率分布,q i是系统在初始时刻处 于状态i的概率,满足。 [编辑] 马尔可夫链模型的性质 马尔可夫链是由一个条件分布来表示的 P(X | X n) n+ 1 这被称为是随机过程中的“转移概率”。这有时也被称作是“一步转移概率”。二、三,以及更多步的转移概率可以导自一步转移概率和马尔可夫性质:

连续隐马尔科夫链模型简介

4.1 连续隐马尔科夫链模型(CHMM) 在交通规划和决策的角度估计特定出行者的确切的出行目的没有必要,推测出行者在一定条件下会有某种目的的概率就能够满足要求。因此本文提出一种基于无监督机器学习的连续隐马尔科夫链模型(CHMM)来识别公共自行车出行链借还车出行目的,根据个人属性、出行时间和站点土地利用属性数据,得到每次借还车活动属于某种出行目的的概率,进一步识别公共自行车出行链最可能的出行目的活动链。 4.1.1连续隐马尔科夫链模型概述 隐马尔可夫链模型(Hidden Markov Model,HMM)是一种统计模型,它被用来描述一个含有隐含未知状态的马尔可夫链。隐马尔可夫链模型是马尔可夫链的一种,其隐藏状态不能被直接观察到,但能通过观测向量序列推断出来,每个观测向量都是通过状态成员的概率密度分布表现,每一个观测向量是由一个具有相应概率密度分布的状态序列产生。 本文将隐马尔科夫链和混合高斯融合在一起,形成一个连续的隐马尔科夫链模型(CHMM),并应用该模型来识别公共自行车出行链借还车活动目的。连续隐马尔科夫链模型采用无监督的机器学习技术,用于训练的数据无需是标记的数据,该模型既不需要标记训练数据,也没有后续的样本测试,如提示-回忆调查。相反,该模型仅利用智能卡和总的土地利用数据。后者为隐藏活动提供额外的解释变量。出行链内各活动的时间和空间信息是从IC卡数据获得,相关土地利用数据是根据南京土地利用规划图和百度地图POI数据获得。 在本文的研究中,一个马尔可夫链可以解释为出行者在两个连续活动状态之间的状态转换,确定一个状态只取决于它之前的状态,一个状态对应一个出行者未知的借还车活动[48-50]。本研究坚持传统的马尔可夫过程的假设,将它包含进无监督的机器学习模型。“隐藏马尔可夫”源于一个事实,即一系列出行链的活动是不可观察的。 对于CHMM,高斯混合模型负责的是马尔可夫链的输入端,每一个活动模式下的隐藏状态都有属于一个特征空间的集群输出概率,每个集群是观察不到的,隐藏状态集群的数量必须事先给出。一些研究者称这些集群为二级隐状态[51]。

1140503102450451连续时间马尔可夫链

5 连续时间马尔可夫链 5.1引言 本章中我们考虑与离散时间马尔可夫链类似的连续时间马尔可夫链。如离散情形一样,它们由马尔可夫性刻画,即已知现在的状态时将来与过去独立。 在5.2节中。我们定义连续时间马尔可夫链且把它们与第四章的离散时间马尔可夫链相联系。在5.3节中,我们引入一类重要的连续时间马尔可夫链,即所谓生灭过程。这些过程可用作在任何时刻其总量的变化仅为一个单位的群体的模型。在5.4节中,我们导出两组描述系统的概率规律的微分方程——向前与向后方程。5.5节的内容是确定连续时间马尔可夫链的有关的极限(或长时间后的)概率。在5.6节中,我们考虑时间可逆的问题。其中,我们证明一切生灭过程是时间可逆的,而后阐明这事实对于排队系统的重要性。在这一节中也提供了时间可逆性对随机群体模型的应用。在5.7节中,我们阐明逆向链的重要性,即使过程不是时间可逆的。利用它我们研究排队网络模型。导出爱尔朗消失公式,分析共用加工系统。5.8节中我们表面如何“一致化”马尔可夫链——对于数值计算有用的一种技巧。 5.2连续时间马尔可夫链 考虑取非负整数值的连续时间随机过程(){}t ,0X t 3,与第四章中给出的离散时间马尔可夫链的定义类似,过程(){}t ,0X t 3称为连续时间马尔可夫链,如 果对一切,0s t 3及非负整数,i j ,()x u ,0u s # ,有 ()()()(){}|X ,X ,0P X t s j s i u x u u s +===? ()(){}|P X t s j X s i =+== 换言之,连续时间马尔可夫链是具有马尔可夫性的随机过程,即已知现在s 时是 状态及一切过去的状态的套件下在将来时刻t s +的状态的条件分布只依赖现在的状态而与过去独立。若又有()(){}|P X t s j X s i +==与s 无关则称连续时间马尔可夫链具有平稳的或其次的转移概率。将假定我们所考虑的马尔可夫链都有平稳转移概率。 假设在某时刻,比如说时刻0,马尔可夫链进入状态i ,而且假设在接下来的s 个单位时间中过程未离开状态i (即未发生转移)。在随后的t 个单位时间中过程仍不离开状态i 的概率是多少呢?为了回答这个问题。注意到因为在时间s 过程处于状态i ,从马尔可夫性得在区间[],s s t +中它仍然处于状态i 的概率正是他处于状态i 至少t 个单位时间的(无条件)概率。也即若以i t 记过程在转移到另一状态之前停留在状态i 的时间,则对一切,0s t 3有 {}{}|i i i P s t s P t t t t >+>=>

马尔可夫链应用于天气预报

马尔可夫链应用于天气预报 摘要: 在《概率论与随机过程》课中学习了马尔可夫链,马尔可夫过程因其无后效性、遍历性和时齐性,在科学研究、天气预测、农业预测、市场预测等方面应用非常广泛。本文通过对马尔可夫链理论和切普曼-柯尔莫哥洛夫方程的探讨,结合天气因素、降水情况的不确定性和无后效性等诸多特点,构建了基于天气预报的马尔可夫链预测模型,文中给出了马尔可夫链的一步转移概率矩阵和多重转移概率的计算方法,根据此算法可以预报短期天气情况,达到预测天气的目的。 关键字:马尔可夫链 天气预报 转移概率 切普曼-柯尔莫哥洛夫方程 1 引言 天气变化情况与人们的生产、生活息息相关,是人们普遍关注的重点问题之一。所以天气预报的准确性与时效性就显得尤为重要,否则将对人们带来不便,甚至有可能带来重大经济和人员损失。本文借助随机过程中著名的马尔可夫链模型,以某日天气的状态转移数据为例,建立了天气情况预测模型,并借助该模型应用马尔可夫链的遍历性,对未来天气的变化趋势作出了预测分析。由于马尔可夫过程应用广泛,它的重要特征是无后效性和遍历性。因此,运用马尔可夫链,只需要最近或现在的动态资料则可按转移概率可预测将来,这样就可以很方便地达到预测天气变化的目的。 2 马尔可夫链预测模型 2.1 马尔可夫链的概念和特性 马尔可夫过程是指具有以下特性的过程:过程X(t)(或系统)在时刻t 0所处的状态为已知的条件下,过程在时刻t >t 0所处状态的条件分布与过程在时刻t 0之前所处的状态无关,只与时刻t 0所处的状态有关,这种特性称为马尔可夫性或无后效性。则称X(t)为马尔可夫过程。 马尔可夫链实际上就是状态和时间都是离散的马尔可夫过程。这一特性可用分布函数来确切地表出:设随机过程{X(t),t ∈T},状态空间为χ,若对于t 的任意n 个值t 1

连续马尔科夫过程的转移概率及应用

《随机过程》 课程设计(论文) 题目: 连续马尔科夫过程的转移 概率及应用 学院:理学院 专业:应用统计学 班级: 13090501 学生姓名:张志达 学生学号: 1309050131 2015年 12 月 29 日

摘要 选取 1978 ~ 2009 年四川农村居民人均生活消费值的 32 个样本,首先,通过 Markov 预测法预测未来生活消费水平的增长速度以 10% ~ 20% 的概率较大; 然后,为提高预测精度,在传统 ARMA 模型中加入时间变量 t 进行建模并预测,预测结果表明平均相对误差率为 1. 56% ,其中 2006 ~ 2009 年的相对误差的绝对值均小于 0. 5% ; 最后,将 Markov 预测和 ARMA 模型对2010 ~ 2012 年的预测结果对比,发现两者在生活消费增长幅度上吻合,预测结果可靠。结果表明,在与目前相似的政策力度下,短期内四川省农村居民消费需求将持续增长,需进一步扩大消费市场。 关键词农村居民; 生活消费; Markov 预测

目录 一.连续马尔科夫过程的转移概率及其应用 (4) 二.连续时间马尔可夫链基本理论 (5) 2.1定义 (5) 2.2转移概率 (5) 三. 马尔可夫过程研究的问题的分析 (7) 数据来源与研究方法 (7) 2.计算状态转移概率矩阵 (8) 3.结果与分析 (10) 四结论和展望 (11) 五.参考文献 (12) 六计算结果及程序 (12)

一.连续马尔科夫过程的转移概率及其应用 1951年前后,伊藤清建立的随机微分方程的理论,为马尔可夫过程的研究开辟了新的道路。1954年前后, W.费勒将半群方法引入马尔可夫过程的研究。流形上的马尔可夫过程、马尔可夫向量场等都是正待深入研究的领域。 类重要的随机过程,它的原始模型马尔可夫链,由俄国数学家Α.Α.马尔可夫于1907年提出。人们在实际中常遇到具有下述特性的随机过程:在已知它目前的状态(现在)的条件下,它未来的演变(将来)不依赖于它以往的演变(过去)。这种已知“现在”的条件下,“将来”与“过去”独立的特性称为马尔可夫性,具有这种性质的随机过程叫做马尔可夫过程。荷花池中一只青蛙的跳跃是马尔可夫过程的一个形象化的例子。青蛙依照它瞬间或起的念头从一片荷叶上跳到另一片荷叶上,因为青蛙是没有记忆的,当现在所处的位置已知时,它下一步跳往何处和它以往走过的路径无关。如果将荷叶编号并用012,,......x x x 分别表示青蛙最初处的荷叶号码及第一次、第二次、……跳跃后所处的荷叶号码,那么{},0n x n ≥ 就是马尔可夫过程。液体中微粒所作的布朗运动,传染病受感染的人数,原子核中一自由电子在电子层中的跳跃,人口增长过程等等都可视为马尔可夫过程。还有些过程(例如某些遗传过程)在一定条件下可以用马尔可夫过程来近似。 关于马尔可夫过程的理论研究,1931年Α.Η.柯尔莫哥洛夫发表了《概率论的解析方法》,首先将微分方程等分析方法用于这类过程,奠定了它的理论基础。1951年前后,伊藤清在P.莱维和C.H.伯恩斯坦等人工作的基础上,建立了随机微分方程的理论,为研究马尔可夫过程开辟了新的道路。1954年前后,W.弗勒将泛函分析中的半群方法引入马尔可夫过程的研究中,Ε.Б.登金(又译邓肯)等并赋予它概率意义(如特征算子等)。50年代初,角谷静夫和J.L.杜布等发现了布朗运动与偏微分方程论中狄利克雷问题的关系,后来G.A.亨特研究了相当一般的马尔可夫过程(亨特过程)与 位势的关系。目前,流形上的马尔可夫过程、马尔可夫场等都是正待深入研究的领域。

第1章离散时间的马尔可夫链

第1章 离散时间的马尔可夫链 §1 随机过程的基本概念 定义1 设(,,)P ΩF 是概率空间,(, )E E 是可测空间,T 是指标集. 若对任何 t T ∈,有:t X E Ω→,且t X ∈F E ,则称{}(), t X t T ω∈是(, , )P ΩF 上的取值于 (,)E E 中的随机过程, 在无混淆的情况下简称{(), }t X t T ω∈为随机过程,称(,)E E 为状态空间或相空间,称E 中的元素为状态,称T 为时间域. 对每个固定的ω∈Ω,称()t X ω为{}(), t X t T ω∈对应于ω的轨道或现实,对每个固定的t T ∈,称()t X ω为E 值随机元. 有时()t X ω也记为 ()()(, )t t X X X t X t ωω=== 设T ?R ,{}, t t T ∈F 是F 中的一族单调增的子σ代数(σ代数流) ,即 ① t t T ?∈??F F ,且t F 是σ代数; ② , , s t s t T s t ?∈

相关主题
文本预览
相关文档 最新文档