您如何看待马尔可夫链是不可约的?


12

我有一些很难理解的马尔可夫链财产束缚

所谓不可简化是指随机过程可以“从任何状态进入任何状态”。

但是,什么定义了它可以从状态进入状态还是不能进入状态?Ĵij


维基百科页面给出了形式化:

状态j可访问的(写入ij)从状态i,如果存在整数nij>0 ST

P(Xnij=j | X0=i)=pij(nij)>0

那么交流是如果ijji

从这些不可还原性以某种方式得出。


关于“可访问性”的直觉是什么?我不明白为什么有条件概率会使某些东西“可访问”?
mavavilj

您可能会从难以接近的角度看待。国家据说是从人迹罕至有没有机会从那里,这是对任何数量的步骤这一事件的概率仍然。为了定义可访问性,应该将量化器(即切换为,将切换到(与相同,因为概率为正)。Ñ 0 ∀ = 0 0 > 0jiin0=00>0
nmerci 2015年

Answers:


12

这是转移矩阵的三个示例,前两个用于可约的情况,最后一个用于不可约的情况。

P1

P1=(0.50.5000.90.100000.20.8000.70.3)P2=(0.10.10.40.40.50.10.10.30.20.40.20.20001)
对于,当您处于状态3或4时,您将停留在那里,并且例如,状态1和状态2相同。无法从状态1进入状态3或4。P1

对于,您可以从状态1到3进入任何状态,但是一旦进入状态4,您将停留在那里。 为此例如,您可以从任何状态开始,但仍然可以达到任何其他状态,尽管不一定要一步一步。P 3 = 0.5 0.5 0 0 0 0 0.9 0 0 0 0 0.1 0 0 0 0.8 0 0.2 0.7 0 0.1 0 0.2 0 0 0 0 0.1 0.9 0 0.9 0 0 0 0.1 0P2

P3=(0.50.500000.900000.10000.800.20.700.100.200000.10.900.90000.10)

5

状态被说成是从一个状态访问(通常表示为),如果存在一些,使得: 即,一个可以从状态得到到状态在与概率的步骤。Ĵ Ñ 0 p Ñ Ĵ = PX Ñ = Ĵ | X 0 = > 0 Ĵ Ñ p Ñ Ĵjiijn0

pijn=P(Xn=jX0=i)>0
ijnpijn

如果和成立,则状态和通信(通常由)。因此,如果两个国家进行通信,则马尔可夫链是不可约的Ĵ Ĵ Ĵijjiijij


是在电源或索引?p Ñ Ĵñpijñ
mavavilj 2015年

这是一个索引。但是,它有一个解释:如果是转移概率矩阵,则是第个元素(此处是幂) 。p n i j i j P n nP=p一世Ĵp一世Ĵñ一世ĴPññ
nmerci 2015年

2

令和为马氏链的两个不同状态。如果进程从状态到状态概率为正,无论步数是多少(例如1、2、3),那么我们就说状态可从状态访问。一世Ĵ一世ĴĴ一世

从概念上讲,我们将其表示为。用概率表示如下:如果存在整数使得,则状态可从状态访问。一世ĴĴ一世>0p一世Ĵ>0

类似地,我们说,如果存在整数,使得。Ĵ一世ñ>0pĴ一世ñ>0

现在,如果和 都为真,那么我们说状态和相互通信,并且用符号表示为。就概率而言,这意味着存在两个整数使得和。 一世ĴĴ一世一世Ĵ一世Ĵp m i j > 0 p n j i >>0ñ>0p一世Ĵ>0pĴ一世ñ>0

如果马尔可夫链中的所有状态都属于一个封闭的沟通类,则该链称为不可约马尔可夫链。不可约性是链的属性。

在不可约的马尔可夫链中,过程可以从任何状态进入任何状态,无论它需要多少步骤。


1

现有的一些答案对我来说似乎是不正确的。

正如J. Medhi的“ 随机过程”(第79页,第4版)中所引用的那样,如果马尔可夫链除了状态空间之外不包含任何适当的“闭合”子集,则它是不可约的。

因此,如果在您的转移概率矩阵中有状态的子集,使得您无法“到达”(或访问)除这些状态之外的任何其他状态,则马尔可夫链是可约的。否则,马尔可夫链是不可约的。


-1

首先要提个警告:除非有充分的理由,否则请不要看矩阵:我唯一想到的就是检查错误键入的数字或阅读教科书。

如果是您的转换矩阵,请计算。如果所有条目都不为零,则矩阵是不可约的。否则,它是可还原的。如果过大,计算与一样大就可以了。相同的测试,准确性稍差。exp P P P n nP经验值PPPññ

不可约性意味着:您可以以有限的步骤从任何状态进入任何其他状态。

在Christoph Hanck的示例,您不能直接从状态1转到状态6,但是可以进入1-> 2-> 6P3


1
您如何定义“可以从状态到状态 ”?Ĵ一世Ĵ
mavavilj 2015年

1
您真的需要问老师。你知道他不会吃你的。
titus 2015年

当您使用exp(P)时,您是在指矩阵指数吗?或,其中i,j是矩阵P的ij项?ËP一世Ĵ
Hunle

我指的是矩阵指数
titus
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.