图解马尔可夫链：从“无记忆”到“预测未来”

1. 当天气预报遇上马尔可夫链

想象一下你正在计划周末的野餐，最关心的就是天气会不会下雨。气象台给出的预报说："如果今天晴天，明天有70%概率继续晴天；如果今天下雨，明天有60%概率继续下雨。"这种预测方式背后，其实就用到了马尔可夫链的核心思想——无记忆性。

我刚开始接触这个概念时，总觉得"无记忆"这个说法太抽象。直到有次观察家里的智能空调才恍然大悟：它根据当前室温自动调节风速，完全不会纠结昨天这个时候的温度是多少。这种"活在当下"的特性，正是马尔可夫链的精髓所在。

1.1 状态转换就像翻牌子

让我们用更形象的例子来说明。假设你面前有两个牌子，A面写着"晴天"，B面写着"雨天"。每天早晨你根据特定规则翻动牌子：

• 如果今天显示A（晴天）：

  • 明天保持A的概率：70%
  • 翻到B的概率：30%

• 如果今天显示B（雨天）：

  • 明天保持B的概率：60%
  • 翻到A的概率：40%

这个翻牌游戏完美展示了状态转移概率。我们可以用矩阵表示这个规律：

实测下来，这种表示方法特别适合计算多步预测。比如已知今天晴天，想知道后天雨天的概率，只需要把矩阵自乘一次：

import numpy as np transition = np.array([[0.7, 0.4], [0.3, 0.6]]) two_day_transition = np.linalg.matrix_power(transition, 2) # 输出 [[0.61 0.52], [0.39 0.48]]

计算结果显示：今天晴天时，后天雨天的概率是39%；今天雨天时，后天雨天的概率是48%。这种预测方式在金融、物流等需要短期预测的领域特别实用。

2. 游戏闯关中的马尔可夫智慧

最近帮我侄子设计了一个简单的棋盘游戏，意外发现马尔可夫链可以用来建模游戏进程。假设玩家每回合可能处于三个状态："起点"、"中途"、"终点"，转移规则如下：

• 起点有80%概率留在原地，20%进入中途
• 中途有30%退回起点，50%保持原位，20%到达终点
• 终点是吸收状态（游戏结束）

2.1 用Python模拟游戏过程

为了验证这个设计是否合理，我写了个模拟程序：

import random def simulate_game(): position = '起点' steps = 0 while position != '终点': steps += 1 if position == '起点': position = random.choices(['起点', '中途'], weights=[0.8, 0.2])[0] elif position == '中途': position = random.choices(['起点', '中途', '终点'], weights=[0.3, 0.5, 0.2])[0] return steps # 模拟10000次游戏 results = [simulate_game() for _ in range(10000)] print(f"平均需要{sum(results)/len(results):.1f}步通关")

跑出来的结果大约是15.3步通关，这个数据帮助我调整了游戏难度。在实际项目中，这种建模方法常被用于用户行为分析，比如预测用户从浏览商品到下单的转化概率。

2.2 稳态分布的秘密

反复运行这个游戏会发现一个有趣现象：无论从哪个状态开始，长期来看停留在各状态的概率会趋于稳定。通过计算转移矩阵的n次方，可以找到这个稳态分布：

transition = np.array([[0.8, 0.3, 0], [0.2, 0.5, 0], [0, 0.2, 1]]) # 计算100次转移后的矩阵 steady_state = np.linalg.matrix_power(transition, 100) print(steady_state[0]) # 输出：[0. 0. 1.]

结果显示最终所有概率都集中在终点状态，这符合我们的游戏设计。在电商场景中，类似的稳态分析可以帮助优化用户转化路径。

3. 文本生成中的马尔可夫魔法

去年做智能客服项目时，我发现一个惊艳的应用：用马尔可夫链生成自然语言。比如分析大量客服对话后，可以建立词语间的转移概率：

• "您好"后面有60%概率接"请问"，30%接"欢迎"，10%接"很高兴"
• "问题"后面有70%概率接"已经"，20%接"正在"，10%接"需要"

3.1 构建词级马尔可夫模型

下面这段代码展示了如何用《红楼梦》文本训练简单生成器：

from collections import defaultdict import random text = "黛玉葬花宝玉哭宝钗劝..." # 实际使用更长文本 model = defaultdict(lambda: defaultdict(int)) # 训练模型 words = list(text) for i in range(len(words)-1): model[words[i]][words[i+1]] += 1 # 生成文本 current = random.choice(words) output = [current] for _ in range(100): next_words = list(model[current].keys()) weights = list(model[current].values()) current = random.choices(next_words, weights=weights)[0] output.append(current) print(''.join(output))

虽然生成的文本缺乏深层语义，但在数据增强方面很实用。我在处理医疗问答短缺问题时，就用这种方法扩展了训练数据。

4. 股市预测中的马尔可夫局限

曾经有段时间我痴迷于用马尔可夫链预测股价，结果踩了个大坑。假设股票有三种状态：涨（+1%）、平（-1%~+1%）、跌（-1%），历史数据统计的转移矩阵如下：

4.1 为什么金融预测会失效

编写回测代码后发现问题：

def simulate_market(initial_state, days): current = initial_state portfolio = 100 for _ in range(days): if current == '涨': portfolio *= 1.01 elif current == '跌': portfolio *= 0.99 current = random.choices( ['涨', '平', '跌'], weights=transition[current] )[0] return portfolio

实际运行发现预测结果与真实市场偏差很大，因为股价变化还受外部政策、市场情绪等因素影响，违背了马尔可夫链的"无记忆性"假设。这个教训让我明白：技术方案必须匹配问题特性。后来改用在高频交易中的极短期预测（秒级），效果才有所改善。