跳到主要内容

稳态分布

一、直观理解：链条转久了，各状态“占比”趋于稳定​

对不可约非周期的有限马尔可夫链，存在唯一稳态分布 π\pi，使 π=πP\pi=\pi P，且从任意初分布出发 π(n)→π\pi^{(n)}\to\pi。长期来看，处于各状态的时间比例收敛到 πi\pi_i。

这回答：“长期有多少比例时间是雨天？”“顾客长期处于哪一服务阶段？”等。

二、适用场景与优缺点​

类型说明
适用场景长期占比、均衡结构；品牌份额稳态；排队网络节点占用的马氏描述。
核心优点线性方程可解；解释力强；可做参数敏感性。
局限与缺点需正则条件；吸收链稳态退化在吸收态；现实可能未达稳态。

三、核心原理​

求 π\pi 满足

πP=π,∑iπi=1,πi≥0\pi P=\pi,\quad \sum_i\pi_i=1,\quad \pi_i\ge0

等价于解 (PT−I)πT=0(P^\mathsf{T}-I)\pi^\mathsf{T}=0 并加归一化。也可用特征值 1 的左特征向量。

混合时间与谱隙有关：∣λ2∣|\lambda_2| 越接近 1，收敛越慢。

四、实践要点​

  • 先检查不可约、非周期；
  • 吸收链改研究吸收概率/时间；
  • 比较 π\pi 与长时间模拟频率；
  • 报告收敛速度直觉。

五、MATLAB 示例​

%% 求解稳态分布
P = [
    0.7 0.2 0.1
    0.3 0.4 0.3
    0.2 0.3 0.5
];

% 方法1: 线性方程
A = (P' - eye(3));
A(end,:) = 1;                 % 替换最后一行为归一化
b = [0;0;1];
pi = A \ b;
pi = pi';                     % 行向量
fprintf("稳态 pi = [%s]\n", num2str(pi, '%.4f '));
fprintf("校验 ||pi P - pi|| = %.3e\n", norm(pi*P - pi));

% 方法2: 迭代
pi_k = [1 0 0];
hist = pi_k;
for k=1:50
    pi_k = pi_k * P;
    hist(end+1,:) = pi_k; %#ok<AGROW>
end
figure;
plot(hist, "LineWidth", 1.5);
legend("\\pi1","\\pi2","\\pi3");
xlabel("步数"); ylabel("概率");
title("分布收敛到稳态"); grid on;

% 长时间频率对照
rng(5);
s=1; cnt=zeros(1,3);
T=20000;
for t=1:T
    s = randsample(1:3,1,true,P(s,:));
    cnt(s)=cnt(s)+1;
end
fprintf("模拟频率 = [%s]\n", num2str(cnt/T, '%.4f '));

六、如何解释结果​

给出 π\pi、含义（长期占比）、收敛图与模拟验证。若某状态 πi\pi_i 很小，说明长期很少停留。政策若改变 PP，应比较新的 π\pi。

七、常见误区​

1. 有吸收态仍谈唯一正稳态​

2. 把瞬时分布当稳态​

3. 不解归一化导致任意倍数​

4. 收敛很慢却只迭代几步​

附：建模流程中的落地清单（稳态分布）​

要把方法真正写进论文，建议补齐以下模块，而不是只贴一个最终数字。

1. 问题对齐与假设​

说明题目中的输入/输出如何映射到本方法，并显式列出关键假设。假设越清楚，后面的敏感性分析越有对象。

2. 关键关系​

π=πP\pi=\pi P 且归一；正则链下长期占比收敛到 π\pi。

3. 补充实现片段​

P=[0.7 0.3;0.4 0.6]; A=[P'-eye(2)]; A(end,:)=1; pi=A\[0;1]; disp(pi');

4. 报告要素​

π\pi、收敛图、模拟频率对照、政策含义。

5. 不适用情形​

吸收链仍谈唯一内点稳态；迭代过少。

6. 与验证章节的衔接​

至少做一项：参数扰动、与基线方法对照、极限情形检查。把“方法正文”和“全局方法论（第15门类）”连起来，结论才站得住。

附：结果解读示例与常见审稿问题​

结果解读示例​

假设已经得到主要数值输出。写作时应把数字翻译回题目语言，例如：

  • “在基准参数下指标为 A；当关键参数变化 ±10% 时，指标波动约 B%，方案排序未变。”
  • “算法在 R 次独立重复中有 S 次收敛到同一结构解，说明结论对初始化不敏感。”
  • “与简化基线相比，本方法在验证集误差降低 Δ，但计算时间增加，权衡后仍采用。”

如果只有表格没有句子，评委很难判断你是否真正理解结果。

参数与数据的最小敏感性​

即使篇幅有限，也建议给出一个小表：选取 1–3 个最关键参数，各取偏低/基准/偏高三档，报告输出变化。这比大段空泛的“模型具有较强鲁棒性”更有说服力。

可复现信息​

在附录或正文脚注中记录：软件版本、随机种子、关键容差、数据截取时间范围。另一人按文档应能得到同量级结果。

审稿式自问​

  1. 方法假设在题目数据上是否可检查？
  2. 若最敏感参数估计错误 20%，结论是否翻盘？
  3. 是否存在更简单方法达到相近效果？
  4. 图表是否都能用一句话概括？
  5. 局限是否写得具体，而不是套话？

把上述自问的答案写进论文，方法章节会明显更“像完成了建模”，而不是“像完成了公式摘录”。

八、小结​

稳态分布描述马尔可夫链的长期占有规律，由 π=πP\pi=\pi P 与归一化唯一确定（在正则条件下）。它是转移矩阵分析的自然下一章，也是排队、评级迁移与份额演化的常用终点指标。