马尔可夫链建模

离散时间马尔可夫链对象框架概述

DTMC对象框架提供了建模和分析离散时间马尔可夫链的基本工具。该对象支持链具有有限数万博1manbetx量的,在离散时间具有时间均相过渡结构演变状态。

DTMC标识每个马尔可夫链与NumStates-通过-NumStates转移矩阵P初始状态的独立X0或状态的初始分布π0。您可以指定P无论是作为右随机矩阵或经验计数的矩阵。

  • 作为右随机矩阵:

    • PIJ是从状态的转变的非负概率一世国家Ĵ

    • 每行P总和为1。

    • π Ť + 1 = π Ť P 描述了从时间状态分布的演变Ť时间Ť+ 1。

    在时间的分布状态ŤπŤ是长度的行向量NumStates

  • 作为经验计数的矩阵,PIJ是次状态观察到的数一世转换到状态Ĵ。该DTMC物体规格化的行P所以,这是一个正确的,随机矩阵。

mcmix函数是一个替代的马尔可夫链对象创建者;它产生与指定的零图案和随机转变概率的链。mcmix非常适合用于创建具有不同的混合时间,用于测试的目的链。

为了显现有向图,有向图或者,具有链相关联,可以使用graphplot目标函数。graphplot类似于情节一个MATLAB的目标函数®有向图对象,但它包括用于分析马尔可夫链结构的附加功能。参数设置亮点通信类(即,强连通有向图的组件)以及影响收敛具体特点,如复发,短暂,和周期性。您可以突出显示转移概率在P通过使用着色强度热图的曲线边缘。

为了显现大型结构中的链,graphplot可冷凝通信类来代表节点。此选项是基于缩合一个的目标函数有向图目的。

分类目标函数是在该图类高亮的数值模拟。分类决定限制行为通信类别的返回特性。状态分类结合图论的算法,如bfsearch(广度优先搜索)一个MATLAB的目标函数图形对象,而是用更直接的矩阵计算特定于马尔可夫链的理论。该子链方法,可以从链作进一步的分析中提取特定通信类。

isreducibleisergodic对象函数给链结构的简明摘要。总之,他们提供了一个独特的极限分布存在的充分必要条件 π * ,其中 π * = π * P π 0 π * 每一个初次分配π0。该渐近目标函数计算 π * ,如果它存在,并使用特征值分析估计混合时间。该eigplot目标函数绘制的本征值P。在此图中的特征值,图表的示例通过返回eigplot

一个障碍收敛周期性。该目标函数通过调整状态惯性(即,通过加权的对角元素消除周期性P),以产生在链“懒惰”的特定量。极限分布是通过这些变化的影响。

模拟重新分配目标函数提供过程的实现方式,因为它从一个指定的初始状态或分布的演变。该辛普劳distplot目标函数提供各种可视化。该图是表示从均匀初始状态分布在开始状态分布的演变的一个分布图的一个例子。

马尔可夫链分析流程

你可以开始在两个方面构建马尔可夫链模型对象:

  • 确定相关的离散状态的过程,然后估计它们之间的转移概率。在最简单的情况下,理论表明链结构和转移矩阵P。在这种情况下,你的理论在实践中,一些并不总是显而易见的从理论如何发挥出感兴趣。一旦你知道P通过传递创建一个Markov链对象PDTMC,它实现了理论链。

  • 如果你有一个过程不太具体的信息,那么你就必须与各国和可行的过渡模式的不同数量的实验来重现实证结果。该mcmix功能提供的洞察链,其可以捕获在数据基本特征的骨架结构。通过迭代过程中,可以调整随机生成的转换矩阵P以满足建模目标。

对于一个经济计量模型构建,选用的最重要的后果P是连锁的渐近行为。要理解这种行为,识别和区分来自经常性状态的过渡状态(一些国家对恢复时的概率趋于零渐近)(那些其恢复时的概率去一个渐近)。短暂和复发是由所有状态以连通的类共享属性。以确定状态是否视觉上是瞬时的或复发,通过马尔可夫链对象到graphplot对象功能,并指定'ColorNodes',真。可替代地,的输出分类目标函数为评价提供的数字工具。该图是与被分类节点有向图的一个例子。

有向图的冷凝视图通过合并每个通信类成简化了这一评价“超级节点”。在冷凝图,可以很容易地通过超级节点的出度识别短暂与复发(一个出度大于0意味着短暂)。束缚链由单一的,必然反复发作,通信类的。Unichains由单个复发类和任意数量的卫星瞬时类的。Unichains保持不可约链的期望限制行为。考虑冷凝图的经常是前体修整不相干的瞬时状态的链。该子链功能修剪瞬态类链。这个数字是在上图的图的压缩视图。

统一的两个主要障碍,限制行为是:

  • 归约中,多于一个的通信类的存在

  • 周期性中,单个类中的倾向周期的子类中

的组合graphplot分类目标函数可以识别这些问题。如果一个链是可还原的,而不是一个unichain,是很常见的分析拆分独立复发类之间或完全重新制定链。如果链是周期性的(也就是说,它包含了周期性的反复类),但整体结构捕获应用程序的基本细节,目标函数提供一种补救措施。懒链扰动的对角线元素P消除周期性,留下渐近性不受影响。

isreducibleisergodic对象的功能总结状态分类。每一个链条具有平稳分布 π * ,其中 π * = π * P , 后果P是随机的,具有一个的特征值。如果链是不可约的,平稳分布是独一无二的。然而,不可还原性,而足够的,不是唯一的必要条件。Unichains也导致其在过渡状态零个概率大众独特的平稳分布。在这方面,国家分类分析是必不可少的,因为isreducible回报真正仅当链作为一个整体由单个连通的类。isreducible回报任意unichains,在这种情况下,你必须决定瞬态类是否是模型的相关部分。

遍历, 要么原始性是不可约和非周期性的组合。遍历链具有独特的极限分布,也就是说,π0要收敛 π * 每一个初次分配π0。你可以决定是否链条,作为一个整体,是遍历使用isergodic。该功能通过评估的唯一复发类标识遍历unichains。A链是周期性的,如果它是束缚而不是遍历,也就是说,如果〜tfirreduc+〜tfergo=,其中tfirreductfergo通过返回isreducibleisergodic, 分别。

一旦您已确认链是遍历的,你可以通过使用确定唯一极限分布渐近目标函数。渐近返回的极限分布 π * 和的混合时间,这对于瞬态行为的衰减的时间常数的估计。该门阶 - 弗罗贝纽斯定理不可约非负矩阵(见[1])是用于解释这些结果是有用的。任何随机矩阵具有一个的谱半径。周期性矩阵期间,ķ,有ķ特征值均匀地分布绕单位圆的ķ团结的根源。在单位圆内的最大特征值的大小决定过渡状态的衰减率。该eigplot目标函数提供的此信息的快速可视化。这个数字是马尔可夫链的一个周期为3的特征值曲线。

不管链的渐近性的,可以研究通过施加有限步骤分析其混合率。该hitprobhittime函数返回命中概率和预期的第一击球次为目标状态的一个子集,从链中的每个状态开始。这两个函数任选绘制与节点颜色指定击球概率或次有向图。此图显示了使用节点的颜色指定预期的第一击球次有向图的一个例子。有向图也表示开始状态是否是远程的目标。

模拟和再分配让你产生上是难以直接从理论上推导链的统计信息。该模拟辛普劳对象的功能,和重新分配distplot对象的功能,这样的分析提供计算和图形工具。模拟中,例如,产生通过链独立随机游动。与模拟在计量经济学工具箱别处对象功能™,相关统计数据的整体平均预报中发挥了重要作用。相应的辛普劳目标函数提供了几种方法来可视化。该图显示后长度的100个随机游动10个步骤,通过在前面的图中的周期马尔可夫链访问的状态的比例。

参考

[1]霍恩,R.,和C·R·约翰逊。矩阵分析。英国剑桥:剑桥大学出版社,1985年。

也可以看看

对象

功能

相关话题