matlab对国内生产总值（GDP）建立马尔可夫链模型（MC）并可视化|附代码数据

原文链接：http://tecdat.cn/?p=17549

最近我们被客户要求撰写关于马尔可夫链模型的研究报告，包括一些图形和统计输出。

本示例说明如何创建并可视化Markov链模型的结构和演化 。考虑从随机转移矩阵中创建马尔可夫链的四状态马尔可夫链，该模型模拟了国内生产总值（GDP）的动态 

创建实际GDP的马尔可夫链模型。指定状态名称。

P = [0.5 0.5 0.0 0.0;
     0.5 0.0 0.5 0.0;
     0.0 0.0 0.0 1.0;
     0.0 0.0 1.0 0.0];
stateNames = ["Regime 1" "Regime 2" "Regime 3" "Regime 4"];
mc = mc(P,'StateNames'

可视化马尔可夫链的一种方法是绘制转移矩阵的热图。

figure;
imagesc(P);
colormap(jet);
colorbar;
axis squar

图片

点击标题查阅往期内容

图片

R语言连续时间马尔科夫链模拟案例 Markov Chains

图片

左右滑动查看更多

图片

01

图片

02

图片

03

图片

04

图片

有向图

有向图将链中的状态显示为节点，并将状态之间的可行转换显示为有向边。

绘制马尔可夫链的默认有向图。

figure;
plot(mc);

图片

通过基于转移概率指定边缘颜色来比较转移概率。

figure;
plot(mc,'ColorEdges'

图片

通过根据状态类型指定节点颜色和标记来识别循环状态和瞬态状态。

figure;
h = hplot(mc,'ColorEdges',true,

图片

低均值状态是瞬态的，最终转移到递归的高均值状态。

节点标签的默认字体大小为8。将字体大小减少到7。

FontSize = 7;

图片

混合图

可以通过在有向图中绘制目标概率和预期的第一次命中时间来可视化。

 从马尔可夫链中的每个状态开始计算命中目标状态的指定子集的概率。其中节点颜色表示命中概率。

绘制马尔可夫链的有向图，其中节点颜色表示命中方案1的概率。

htp(mc,"Regime 1",'Graph

图片

 从马尔可夫链中的每个状态开始，计算目标状态的指定子集的预期首次命中时间。其中节点颜色表示命中时间。

绘制马尔可夫链的有向图，其节点颜色表示包含状态3和4的目标子类的预期首次命中时间。

target = ["Regime 3" "Regime 4"];
htime(mc,target

图片

从方案1开始，该子类的预期首次命中时间为6个时间步长。

特征值图

特征值图显示了复平面上的特征值。特征值图并标识：

• Perron-Frobenius特征值，使用粗体星号为非负矩阵。
• 光谱间隙，即长度等于第二大特征值幅度（SLEM）的半径和长度为1的半径之间的区域。光谱间隙决定了马尔可夫链的混合时间。大间隙表示混合较快，而细间隙表示混合较慢。

在复杂平面上绘制并返回转换矩阵的特征值。

figure;
eigplt(mc)

图片

eVals = 4×1

    0.8090
   -0.3090
    1.0000
   -1.0000

两个特征值的模量为1，表明马尔可夫链的周期为2。

重新分配图

重新分布图从初始分布绘制了状态重新分布。 使用马尔可夫链对象生成的数据来绘制重新分布 。可以将重新分布绘制为静态热图或动画直方图或有向图。

从初始分布生成10步重新分布。

redis(mc,numSteps,'X0',x0);

将重新分布绘制为热图。

图片

由于状态1和状态2是瞬态的，因此马尔可夫链最终将概率集中在状态3和状态4。此外，如特征值图所示，状态3和状态4的周期为2。

绘制动画直方图。将帧速率设置为一秒。

图片

仿真图

仿真图绘制了从特定初始状态开始的马尔可夫链的随机游动图。 

生成100个十步随机游走，其中每个状态都会初始化游走25次。

simu(mc,numSteps,

将模拟绘制为热图，以显示每个步骤达到的状态比例。

图片

绘制已实现转移矩阵的热图。

figure;
simp(mc,X,'Type','transition');

图片

所实现的转移矩阵看起来类似于理论转移矩阵。

图片

点击文末 “阅读原文”

获取全文完整代码数据资料。

本文选自《matlab对国内生产总值（GDP）建立马尔可夫链模型（MC）并可视化》。

点击标题查阅往期内容

PYTHON用时变马尔可夫区制转换（MARKOV REGIME SWITCHING）自回归模型分析经济时间序列 R语言中实现马尔可夫链蒙特卡罗MCMC模型matlab贝叶斯隐马尔可夫hmm模型实现 R语言如何做马尔科夫转换模型markov switching model matlab中的隐马尔可夫模型(HMM)实现 R语言马尔可夫体制转换模型Markov regime switching R语言马尔可夫转换模型研究交通伤亡人数事故预测 R语言中的马尔科夫机制转换(Markov regime switching)模型 用机器学习识别不断变化的股市状况—隐马尔科夫模型(HMM)股票指数预测实战 用机器学习识别不断变化的股市状况—隐马尔科夫模型(HMM)的应用 R语言中实现马尔可夫链蒙特卡罗MCMC模型用R语言模拟M / M / 1随机服务排队系统 用R语言模拟混合制排队随机服务排队系统 用R语言模拟M/M/1随机服务排队系统