Files
2026-05-16 17:16:51 +08:00

925 lines
29 KiB
Markdown
Raw Permalink Blame History

This file contains ambiguous Unicode characters
This file contains Unicode characters that might be confused with other characters. If you think that this is intentional, you can safely ignore this warning. Use the Escape button to reveal them.
---
title: 05-多智能体强化学习
tags:
- 多智能体
- 强化学习
- 博弈论
- 深度强化学习
---
# 多智能体强化学习
> 本章作者:多智能体强化学习课题组
> 更新日期2026-05-14
> 参考资料MARL 综述 (Zhang et al., 2021), Dec-POMDP (Oliehoek & Amato, 2016)
---
## 1. 多智能体强化学习基础
### 1.1 博弈论基础
#### 1.1.1 正常形博弈Normal-Form Game
**定义 1.1(正常形博弈)** 一个 $n$ 智能体的正常形博弈定义为元组:
$$
\mathcal{G} = \langle \mathcal{N}, \mathcal{A}, \mathcal{R} \rangle
$$
其中:
- $\mathcal{N} = \{1, 2, \ldots, n\}$ 为智能体集合
- $\mathcal{A} = \mathcal{A}_1 \times \mathcal{A}_2 \times \cdots \times \mathcal{A}_n$ 为联合动作空间,$\mathcal{A}_i$ 为智能体 $i$ 的动作空间
- $\mathcal{R} = (r_1, r_2, \ldots, r_n)$ 为收益函数向量,其中 $r_i: \mathcal{A} \rightarrow \mathbb{R}$
**定义 1.2(策略与混合策略)** 智能体 $i$ 的**纯策略**是 $\mathcal{A}_i$ 中的一个动作 $a_i$。**混合策略**是定义在 $\mathcal{A}_i$ 上的概率分布:
$$
\sigma_i(a_i) \geq 0, \quad \sum_{a_i \in \mathcal{A}_i} \sigma_i(a_i) = 1
$$
联合混合策略为 $\boldsymbol{\sigma} = (\sigma_1, \ldots, \sigma_n)$,在不至于混淆时简记为 $\sigma$。
#### 1.1.2 扩展形博弈Extensive-Form Game
**定义 1.3(扩展形博弈)** 一个有限完美信息扩展形博弈定义为元组:
$$
\mathcal{G} = \langle \mathcal{N}, \mathcal{H}, \mathcal{A}, \mathcal{Z}, \mathcal{R}, \rho \rangle
$$
其中:
- $\mathcal{H}$ 为历史节点集合(根节点为空序列 $\varnothing$
- $\mathcal{A}$ 为动作集合
- $\mathcal{Z} \subseteq \mathcal{H}$ 为终止节点集合
- $\rho: \mathcal{H} \setminus \mathcal{Z} \rightarrow \mathcal{N}$ 为将每个决策节点分配给某智能体的函数
- $\mathcal{R}: \mathcal{Z} \rightarrow \mathbb{R}^n$ 为终端收益函数
对于**不完美信息**博弈引入信息集合Information Set的概念
**定义 1.4(信息集合)** 智能体 $i$ 的信息集合 $\mathcal{I}_i$ 是其决策节点的集合,满足:玩家在 $\mathcal{I}_i$ 中任意节点上观察到相同的可用信息。
#### 1.1.3 纳什均衡Nash Equilibrium
**定义 1.5(纳什均衡)** 混合策略组合 $\sigma^*$ 构成**纯策略纳什均衡**当且仅当对每个智能体 $i$
$$
a_i^* \in \arg\max_{a_i \in \mathcal{A}_i} r_i(a_i, \sigma_{-i}^*)
$$
**定理 1.1纳什存在性Nash 1950** 任意有限正常形博弈至少存在一个混合策略纳什均衡。
*证明思路*:将混合策略空间视为 $\Delta(\mathcal{A}_i)$ 的笛卡尔积(紧凸集),定义最佳响应映射 $BR(\sigma_{-i}) = \arg\max_{\sigma_i} u_i(\sigma_i, \sigma_{-i})$。该映射为上半连续、凸值的 Kakutani 不动点定理适用,故存在不动点即纳什均衡。 $\square$
#### 1.1.4 零和博弈与一般和博弈
**定义 1.6(零和博弈)** 若博弈 $\mathcal{G}$ 满足 $\sum_{i \in \mathcal{N}} r_i(\mathbf{a}) = 0, \forall \mathbf{a} \in \mathcal{A}$,则称为**零和博弈**。此时可简记收益为 $r_1(\mathbf{a}) = -r_2(\mathbf{a})$(二人博弈)。
**定义 1.7(一般和博弈)** 若博弈不满足零和条件,则为**一般和博弈**。一般和博弈的纳什均衡可能帕累托次优Pareto Suboptimal
---
### 1.2 Dec-POMDP 形式化
**定义 1.8Dec-POMDP** Decentralized Partially Observable Markov Decision Process 定义为元组:
$$
\mathcal{M} = \langle \mathcal{N}, \mathcal{S}, \mathcal{A}, \mathcal{P}, \mathcal{R}, \Omega, \mathcal{O}, \gamma, T \rangle
$$
其中:
- $\mathcal{N} = \{1, 2, \ldots, n\}$:智能体数量
- $\mathcal{S}$:状态空间(全局状态)
- $\mathcal{A} = \prod_{i=1}^n \mathcal{A}_i$:联合动作空间,$\mathcal{A}_i$ 为智能体 $i$ 的动作空间
- $\mathcal{P}: \mathcal{S} \times \mathcal{A} \times \mathcal{S} \rightarrow [0, 1]$:状态转移概率,$P(s'|s, \mathbf{a})$
- $\mathcal{R}: \mathcal{S} \times \mathcal{A} \rightarrow \mathbb{R}$:联合奖励函数
- $\Omega$:观测空间,$\Omega = \prod_{i=1}^n \Omega_i$
- $\mathcal{O}: \mathcal{S} \times \Omega \rightarrow [0,1]$:观测函数,$O(o|s)$ 表示在状态 $s$ 下产生观测 $o$ 的概率
- $\gamma \in [0, 1)$:折扣因子
- $T$:规划时域(有限时域 MDP 可设 $T$ 为有限值)
**关键约束**:各智能体只能基于**局部观测历史**做出决策,无法直接访问全局状态 $s$ 或其他智能体的观测/动作。
#### 1.2.1 联合策略与历史
**定义 1.9(个体策略)** 智能体 $i$ 的**随机策略**定义为:
$$
\pi_i: \mathcal{H}_i \rightarrow \Delta(\mathcal{A}_i)
$$
其中 $\mathcal{H}_i$ 为智能体 $i$ 的观测历史空间,$\Delta(\mathcal{A}_i)$ 为动作空间上的概率分布。
**定义 1.10(联合策略)** 联合策略为 $\boldsymbol{\pi} = (\pi_1, \ldots, \pi_n)$。
**定义 1.11(联合历史策略)** 联合历史 $\mathbf{h}_t = (h_t^1, \ldots, h_t^n)$ 下的动作分布为:
$$
\boldsymbol{\pi}(\mathbf{a}_t | \mathbf{h}_t) = \prod_{i=1}^n \pi_i(a_t^i | h_t^i)
$$
#### 1.2.2 联合价值函数
**定义 1.12(状态-动作价值函数)** 给定联合策略 $\boldsymbol{\pi}$,定义状态-动作价值函数:
$$
Q^{\boldsymbol{\pi}}(s, \mathbf{a}) = \mathbb{E}_{s' \sim P, \mathbf{a}' \sim \boldsymbol{\pi}} \left[ \sum_{t=0}^{\infty} \gamma^t r(s_t, \mathbf{a}_t) \Big| s_0 = s, \mathbf{a}_0 = \mathbf{a} \right]
$$
**定义 1.13(值函数)** 状态值函数:
$$
V^{\boldsymbol{\pi}}(s) = \mathbb{E}_{\mathbf{a} \sim \boldsymbol{\pi}(\cdot|s)} [ Q^{\boldsymbol{\pi}}(s, \mathbf{a}) ] = \sum_{\mathbf{a}} \boldsymbol{\pi}(\mathbf{a}|s) Q^{\boldsymbol{\pi}}(s, \mathbf{a})
$$
**定理 1.2Bellman 方程)** 联合价值函数满足:
$$
Q^{\boldsymbol{\pi}}(s, \mathbf{a}) = \mathcal{R}(s, \mathbf{a}) + \gamma \sum_{s'} P(s'|s, \mathbf{a}) V^{\boldsymbol{\pi}}(s')
$$
$$
V^{\boldsymbol{\pi}}(s) = \sum_{\mathbf{a}} \boldsymbol{\pi}(\mathbf{a}|s) \left[ \mathcal{R}(s, \mathbf{a}) + \gamma \sum_{s'} P(s'|s, \mathbf{a}) V^{\boldsymbol{\pi}}(s') \right]
$$
#### 1.2.3 有限时域与折扣形式
**有限时域Horizon $T$**
$$
J^{\boldsymbol{\pi}}(\mathcal{M}) = \mathbb{E} \left[ \sum_{t=0}^{T-1} \gamma^t r(s_t, \mathbf{a}_t) \right]
$$
**折扣形式Discounted**
$$
J^{\boldsymbol{\pi}}(\mathcal{M}) = \mathbb{E} \left[ \sum_{t=0}^{\infty} \gamma^t r(s_t, \mathbf{a}_t) \right]
$$
两类目标函数在理论上具有等价性:可通过拉格朗日乘子法相互转换。
---
### 1.3 信用分配问题
#### 1.3.1 问题的数学描述
在多智能体系统中,联合奖励 $r(s, \mathbf{a})$ 是所有智能体动作的函数。当团队获得一个高(低)奖励时,如何确定每个智能体对此的贡献?
**定义 1.14(信用分配)** 给定联合奖励 $r_t = r(s_t, \mathbf{a}_t)$,信用分配问题旨在为每个智能体 $i$ 分配**信用值** $c_t^i$,使得:
$$
\sum_{i=1}^n c_t^i = r_t, \quad \forall t
$$
且 $c_t^i$ 反映了智能体 $i$ 对联合奖励的**因果贡献**。
#### 1.3.2 Shapley 值方法
**定义 1.15Shapley 信用分配)** 智能体 $i$ 的 Shapley 值为:
$$
\phi_i(r, \mathcal{N}) = \sum_{\mathcal{S} \subseteq \mathcal{N} \setminus \{i\}} \frac{|\mathcal{S}|! (n - |\mathcal{S}| - 1)!}{n!} \left[ r(\mathcal{S} \cup \{i\}) - r(\mathcal{S}) \right]
$$
其中 $r(\mathcal{S})$ 为子集 $\mathcal{S}$ 中智能体单独工作时的奖励贡献。
Shapley 值满足:
1. **效率性**$\sum_i \phi_i = r(\mathcal{N})$
2. **对称性**:对称智能体具有相同信用
3. **可加性**$r = r_1 + r_2 \Rightarrow \phi(r) = \phi(r_1) + \phi(r_2)$
4. **虚拟性**:空集贡献为零
#### 1.3.3 反事实梯度估计
**定义 1.16(反事实信用分配)** 智能体 $i$ 的反事实优势:
$$
A_i(s, \mathbf{a}) = Q(s, \mathbf{a}) - \mathbb{E}_{a_i' \sim \pi_i(\cdot|s)} [Q(s, (a_i', \mathbf{a}_{-i}))]
$$
**定理 1.3Counterfactual Multi-Agent Policy Gradient** 智能体 $i$ 的策略梯度:
$$
\nabla_{\theta_i} J \approx \mathbb{E}_{s \sim d^{\pi}, \mathbf{a} \sim \pi} \left[ \nabla_{\theta_i} \log \pi_i(a_i|s) \cdot A_i(s, \mathbf{a}) \right]
$$
---
## 2. 值函数分解方法
### 2.1 VDN值分解网络
#### 2.1.1 核心思想
VDN (Sunehag et al., 2017) 假设联合价值函数可分解为各智能体价值函数的**线性求和**
$$
Q_{\mathrm{KL}^{\boldsymbol{\pi}}(s, \mathbf{a}) = \sum_{i=1}^n Q_i^{\pi_i}(s, a_i)
$$
#### 2.1.2 数学推导
**定理 2.1VDN 分解条件)** 若存在分解使得 $Q_{\mathrm{KL} = \sum_i Q_i$,则在 $\gamma$-折扣无限时域下,最优联合策略 $\boldsymbol{\pi}^*$ 可通过独立优化每个 $Q_i$ 获得。
*证明*:设 $\boldsymbol{\pi}^* = \arg\max_{\boldsymbol{\pi}} V^{\boldsymbol{\pi}}(s)$。对任意 $s$
$$
V^{\boldsymbol{\pi}^*}(s) = \max_{\boldsymbol{\pi}} \sum_{\mathbf{a}} \boldsymbol{\pi}(\mathbf{a}|s) Q_{\mathrm{KL}}(s, \mathbf{a})
$$
代入分解假设:
$$
= \max_{\boldsymbol{\pi}} \sum_{\mathbf{a}} \prod_i \pi_i(a_i|s) \sum_i Q_i(s, a_i) = \max_{\boldsymbol{\pi}} \sum_i \sum_{a_i} \pi_i(a_i|s) Q_i(s, a_i)
$$
由于求和可分离,最优解为 $\pi_i^*(\cdot|s) = \delta(a_i = \arg\max_{a_i'} Q_i(s, a_i'))$,即独立优化各智能体。$\square$
#### 2.1.3 VDN 网络结构
VDN 采用**相加结构**
$$
Q_{\mathrm{KL}}(s, \mathbf{a}; \boldsymbol{\theta}) = \sum_{i=1}^n f_i(s, a_i; \theta_i)
$$
其中每个 $f_i$ 为智能体 $i$ 的局部 Q 网络,参数 $\theta_i$ 独立。
**训练目标TD-error**
$$
\mathcal{L}(\boldsymbol{\theta}) = \mathbb{E}_{s, \mathbf{a}, r, s'} \left[ (y - Q_{\mathrm{KL}}(s, \mathbf{a}; \boldsymbol{\theta}))^2 \right]
$$
其中 $y = r + \gamma \sum_i Q_i(s', a_i'; \boldsymbol{\theta}^-)$。
#### 2.1.4 VDN 的局限性
VDN 的线性加和结构具有**表达力不足**的问题:无法建模智能体间的**非线性交互**。例如:
$$
Q(s, a_1, a_2) = \mathbb{1}_{a_1 = a_2} - \mathbb{1}_{a_1 \neq a_2}
$$
此函数不可分解为 $Q_1(a_1) + Q_2(a_2)$ 的形式。
---
### 2.2 QMIX单调性约束分解
#### 2.2.1 问题设定
QMIX (Rashid et al., 2018) 引入**混合网络**Mixing Network以建模非线性分解
$$
Q_{\mathrm{KL}}(s, \mathbf{a}) = f_{\mathrm{KL}}(Q_1(s, a_1), \ldots, Q_n(s, a_n); s)
$$
其中 $f_{\mathrm{KL}}$ 为非线性混合网络,以全局状态 $s$ 为条件。
#### 2.2.2 单调性约束
**定义 2.1(单调性约束)** 为保证**因子化最优策略**Factorized Optimal Policy混合网络需满足
$$
\frac{\partial f_{\mathrm{KL}}{\partial Q_i} \geq 0, \quad \forall i
$$
即 $f_{\mathrm{KL}}$ 对每个 $Q_i$ **非递减**
**定理 2.2QMIX 因子化最优性)** 在单调性约束下,**全局最优动作**可通过**独立贪婪**各智能体获得:
$$
\arg\max_{\mathbf{a}} Q_{\mathrm{KL}}(s, \mathbf{a}) = \left( \arg\max_{a_1} Q_1(s, a_1), \ldots, \arg\max_{a_n} Q_n(s, a_n) \right)
$$
*证明*:设 $\mathbf{a}^* = \arg\max_{\mathbf{a}} f_{\mathrm{KL}}(Q_1(s, a_1), \ldots, Q_n(s, a_n))$。由于 $\partial f_{\mathrm{KL} / \partial Q_i \geq 0$$f_{\mathrm{KL}}$ 对每个 $Q_i$ 非递减。因此,增大任意 $Q_i$ 不会降低 $Q_{\mathrm{KL}}$。故:
$$
a_i^* \in \arg\max_{a_i} Q_i(s, a_i) \Rightarrow \mathbf{a}^* = (a_1^*, \ldots, a_n^*)
$$
$\square$
#### 2.2.3 混合网络结构
混合网络采用**超参数网络**Hypernetwork结构生成线性层的权重和偏置
$$
f_{\mathrm{KL}}(Q_1, \ldots, Q_n; s) = \mathbf{h}^\top \mathbf{v} + b
$$
其中 $\mathbf{h} = \text{MLP}([Q_1, \ldots, Q_n])$$\mathbf{v}, b$ 由 $\text{HyperNet}(s)$ 生成。
**单调性实现**:使用绝对值激活确保权重非负:
$$
W_i^+ = |W_i|, \quad b_i^+ = \text{ELU}(b_i) + 1
$$
#### 2.2.4 QMIX 训练目标
使用 **DRRN** 风格的 TD-error
$$
\mathcal{L}(\boldsymbol{\theta}) = \mathbb{E}_{(s, \mathbf{a}, r, s') \sim \mathcal{D}} \left[ (r + \gamma \bar{Q}_{\mathrm{KL}}(s', \mathbf{a}'; \bar{\boldsymbol{\theta}}) - Q_{\mathrm{KL}}(s, \mathbf{a}; \boldsymbol{\theta}))^2 \right]
$$
其中 $\bar{Q}_{\mathrm{KL} = f_{\mathrm{KL}(\bar{Q}_1, \ldots, \bar{Q}_n; s')$ 为目标网络。
---
### 2.3 QTRAN一般化值函数分解
#### 2.3.1 QTRAN 的动机
VDN 要求线性分解QMIX 要求单调性约束,两者均为**结构化分解**。QTRAN (Sunehag et al., 2018) 追求**一般化分解**,无结构假设。
#### 2.3.2 约束优化框架
QTRAN 将值函数分解形式化为**约束优化**问题。
**定义 2.2(可分解值函数)** 若存在局部价值函数 $Q_i$ 使得:
$$
Q_{\mathrm{KL}}(s, \mathbf{a}) = \sum_i Q_i(s, a_i), \quad \forall (s, \mathbf{a})
$$
则称其为**可分解的**。
**定义 2.3(约束优化形式)** QTRAN 寻找满足以下约束的 $Q_i$
1. **一致性约束Consistency**
$$
\sum_i Q_i(s, a_i) = Q_{\mathrm{KL}}(s, \mathbf{a}), \quad \forall (s, \mathbf{a}) \in \mathcal{D}
$$
2. **最优性约束Optimality**
$$
\sum_i Q_i(s, a_i^*) = \max_{\mathbf{a}} Q_{\mathrm{KL}}(s, \mathbf{a}), \quad \forall s
$$
其中 $a_i^* = \arg\max_{a_i} Q_i(s, a_i)$。
#### 2.3.3 惩罚项形式
将约束融入损失函数:
$$
\mathcal{L} = \underbrace{\sum_{s, \mathbf{a}} (r(s, \mathbf{a}) + \gamma \bar{V}(s') - \sum_i Q_i(s, a_i))^2}_{\mathrm{KL} + \lambda_1 \mathcal{L}_{\mathrm{KL}} + \lambda_2 \mathcal{L}_{\mathrm{KL}
$$
其中:
$$
\mathcal{L}_{\mathrm{KL}} = \sum_s \left( \sum_i Q_i(s, a_i^*) - \bar{Q}_{\mathrm{KL}}(s, \mathbf{a}^*) \right)^2
$$
$$
\mathcal{L}_{\mathrm{KL}} = \sum_{s, \mathbf{a} \neq \mathbf{a}^*} \left( \sum_i Q_i(s, a_i) - \bar{Q}_{\mathrm{KL}}(s, \mathbf{a}) \right)^2
$$
---
### 2.4 信用分配的数学描述
#### 2.4.1 基于优势函数的信用分配
**定义 2.4(多智能体优势函数)** 智能体 $i$ 在联合动作 $\mathbf{a}$ 相对于其边缘最优动作的优势:
$$
A_i(s, \mathbf{a}) = Q_{\mathrm{KL}}(s, \mathbf{a}) - V_{-i}(s, a_i)
$$
其中 $V_{-i}(s, a_i) = \mathbb{E}_{\mathbf{a}_{-i} \sim \pi_{-i}} [Q_{\mathrm{KL}}(s, a_i, \mathbf{a}_{-i})]$。
#### 2.4.2 团队游戏与竞争游戏的信用分配
**团队游戏Team Game**$\mathcal{R}(s, \mathbf{a}) = r(s)$ 与动作无关。信用分配简化为:
$$
\phi_i(s, \mathbf{a}) = \mathbb{E}_{\pi} \left[ \frac{\partial \log \pi_i(a_i)}{\partial \theta_i} \cdot r(s) \right]
$$
**竞争游戏Competitive Game**:采用**极小极大**Minimax原则
$$
\pi_i^* = \arg\max_{\pi_i} \min_{\pi_{-i}} \mathbb{E}[r_i(\pi_i, \pi_{-i})]
$$
---
## 3. 策略梯度方法
### 3.1 MADDPG多智能体 DDPG
#### 3.1.1 中心化 Critic + 去中心化 Actor
MADDPG (Lowe et al., 2017) 的核心设计:**中心化 Critic** 访问全局状态,**去中心化 Actor** 仅基于局部观测。
**定义 3.1(中心化 Critic** 智能体 $i$ 的 Critic $Q_i^\text{cen}$ 以全局状态-动作对 $(s, \mathbf{a})$ 为输入:
$$
Q_i^\text{cen}(s, \mathbf{a}; \phi_i) : \mathcal{S} \times \mathcal{A} \rightarrow \mathbb{R}
$$
**定义 3.2(去中心化 Actor** 智能体 $i$ 的 Actor $\pi_i$ 以局部观测 $o_i$ 为输入:
$$
\pi_i(o_i; \theta_i) : \Omega_i \rightarrow \Delta(\mathcal{A}_i)
$$
#### 3.1.2 数学推导
**定理 3.1MADDPG 策略梯度)** 智能体 $i$ 的策略梯度为:
$$
\nabla_{\theta_i} J(\theta_i) = \mathbb{E}_{s \sim \mathcal{D}, \mathbf{a} \sim \pi} \left[ \nabla_{\theta_i} \log \pi_i(a_i|o_i) \cdot \nabla_{a_i} Q_i^\text{cen}(s, \mathbf{a}) \Big|_{a_i = \pi_i(o_i)} \right]
$$
*证明*:策略梯度定理的多智能体推广。对 $J(\theta_i) = \mathbb{E}_{s_0 \sim \rho_0, a_0 \sim \pi} [Q_i^\text{cen}(s_0, \mathbf{a}_0)]$ 求导:
$$
\nabla_{\theta_i} J = \mathbb{E}_{s \sim d^{\pi}, \mathbf{a} \sim \pi} \left[ \nabla_{\theta_i} \log \pi_i(a_i|o_i) \cdot Q_i^\text{cen}(s, \mathbf{a}) \right]
$$
使用 Critic 的梯度 $\nabla_{a_i} Q_i^\text{cen}$ 作为优势函数的估计(类似 DDPG 的确定性策略梯度)。$\square$
#### 3.1.3 MADDPG 算法
**算法 1MADDPG**
```
初始化n个智能体的Actor网络 θ_i 和Critic网络 φ_i目标网络参数 θ_i^-, φ_i^-
重放缓冲区 D
for episode do
初始化环境,获得初始观测 o = (o_1, ..., o_n)
for t = 1 to T do
每个智能体i根据 π_i(o_i; θ_i) 选择动作 a_i
执行联合动作 a = (a_1, ..., a_n),获得奖励 r环境转移到 s'
存储 (s, a, r, s') 到 D
s ← s'
end for
for 每个智能体 i do
从 D 中采样批量 (s, a, r, s')
// 更新 Critic
y_i = r_i + γ Q_i'^cen(s', a'_1, ..., a'_n; φ_i^-)
更新 φ_i 最小化 (y_i - Q_i^cen(s, a; φ_i))^2
// 更新 Actor
使用策略梯度 ∇_{θ_i} J 更新 θ_i
end for
更新目标网络θ_i^- ← τ θ_i + (1-τ) θ_i^-φ_i^- ← τ φ_i + (1-τ) φ_i^-
end for
```
#### 3.1.4 中心化 Critic 的设计考量
**输入设计**Critic 接收完整的联合动作 $\mathbf{a} = (a_1, \ldots, a_n)$,这允许建模智能体间的交互。
**维度诅咒**:联合动作空间 $|\mathcal{A}| = \prod_i |\mathcal{A}_i|$ 随智能体数量指数增长。实践中需使用**参数共享**或**注意力机制**。
---
### 3.2 MAPPO多智能体 PPO
#### 3.2.1 总体架构
MAPPO (Yu et al., 2022) 将 PPO 扩展到多智能体设置,沿用中心化 Critic + 去中心化 Actor 范式。
**定义 3.3GAE for MARL** 智能体 $i$ 的 $\lambda$-回报:
$$
\hat{G}_t^{i, \lambda} = (1 - \lambda) \sum_{l=0}^{\infty} (\gamma \lambda)^l r_{t+l} + (\gamma \lambda)^{L} V^{\pi_i}(s_{t+L})
$$
广义优势估计GAE可扩展到多智能体设置每个智能体使用**相同的价值估计** $V(s_t)$(来自中心化 Critic
#### 3.2.2 CLIP 目标
$$
L^{\mathrm{KL}(\theta) = \mathbb{E}_t \left[ \min \left( r_t(\theta) \hat{A}_t, \text{clip}(r_t(\theta), 1-\epsilon, 1+\epsilon) \hat{A}_t \right] \right)
$$
其中 $r_t(\theta) = \frac{\pi_i(a_i^t|o_i^t; \theta)}{\pi_i(a_i^t|o_i^t; \theta_{\mathrm{KL}})}$。
---
### 3.3 通信与信息共享
#### 3.3.1 通信协议建模
**定义 3.4(通信信道)** 智能体 $i$ 到 $j$ 的通信信道建模为:
$$
c_{ij}^t = h_i(o_i^t, a_i^{t-1}; \xi_i) + \epsilon_{ij}
$$
其中 $h_i$ 为通信网络,$\epsilon_{ij} \sim \mathcal{N}(0, \sigma^2)$ 为噪声。
#### 3.3.2 可微通信层
**CommNet** (Foerster et al., 2016) 使用连续通信:
$$
c_i^t = \frac{1}{n-1} \sum_{j \neq i} h_j(o_j^t, a_j^{t-1})
$$
$$
a_i^t = f_i(c_i^t, o_i^t)
$$
**BIC** (Hong et al., 2018) 引入**瓶颈层**强制通信压缩:
$$
c_i^t = \text{Bottleneck}(o_i^t) \in \mathbb{R}^k, \quad k \ll |\Omega_i|
$$
---
## 4. 均值场博弈
### 4.1 大量智能体的近似
#### 4.1.1 均值场近似的动机
当智能体数量 $n \rightarrow \infty$ 时,直接建模变得不可行。**均值场理论**通过近似邻居的平均效应来解决这个问题。
**定义 4.1(均值场)** 定义邻居平均动作:
$$
\bar{a}_i^t = \frac{1}{n-1} \sum_{j \neq i} a_j^t \approx \mathbb{E}_{j \sim \text{neigh}(i)} [a_j]
$$
#### 4.1.2 均值场博弈定义
**定义 4.2均值场博弈MFG** 一个连续时间均值场博弈定义为:
- **状态演化**
$$
\dot{x}_i(t) = f(x_i(t), a_i(t), \bar{a}(t)), \quad \bar{a}(t) = \mathbb{E}[a_j(t)]
$$
- **Hamilton-Jacobi-Bellman 方程**
$$
0 = \frac{\partial V}{\partial t} + \max_{a_i} \left\{ r(x_i, a_i, \bar{a}) + \nabla_x V \cdot f(x_i, a_i, \bar{a}) \right\}
$$
- **Fokker-Planck 方程**(人口分布演化):
$$
\frac{\partial \mu}{\partial t} + \nabla_a \cdot (\mu \cdot \bar{a}) = 0
$$
其中 $\mu$ 为状态分布。
### 4.2 离散时间均值场 Q 学习
#### 4.2.1 均值场 Q 函数
**定义 4.3(均值场 Q 函数)** 智能体 $i$ 的均值场 Q 函数:
$$
Q_i^{\mathrm{KL}}(s, a_i) \approx r(s, a_i, \bar{a}) + \gamma \mathbb{E}_{s' \sim P} \left[ V^{\mathrm{KL}}(s') \right]
$$
其中 $\bar{a} = \frac{1}{n} \sum_j a_j$ 为均值动作。
#### 4.2.2 均值场近似误差
**定理 4.1(均值场逼近误差)** 设 $n$ 个智能体的真实 Q 函数为 $Q_i^{(n)}(s, \mathbf{a})$,均值场近似为 $Q_i^{\mathrm{KL}}(s, a_i)$,则:
$$
|Q_i^{(n)}(s, \mathbf{a}) - Q_i^{\mathrm{KL}}(s, a_i)| \leq \frac{C}{\sqrt{n}}
$$
其中常数 $C$ 依赖于奖励函数的 Lipschitz 常数。
*证明思路*:使用 McDiarmid 不等式和 Stein 方法,可证明均值场近似的 $O(1/\sqrt{n})$ 收敛率。$\square$
#### 4.2.3 MFQ 算法
**算法 2Mean Field Q-Learning**
```
初始化Q网络μ网络均值场估计
for episode do
初始化均值场估计 μ(s)
for t = 1 to T do
每个智能体i基于 Q_i(s, ·) 和 μ(s) 选择动作 a_i^t
执行联合动作,获得奖励 r_t^i 和下一状态 s_{t+1}
更新均值场估计μ_{t+1}(s) ← α μ_t(s) + (1-α) \bar{a}_t
// 更新 Q 函数
y_i = r_i(s, a_i, μ(s)) + γ V_i(s', μ(s'))
更新 Q_i 通过 TD: (y_i - Q_i(s, a_i))^2
end for
end for
```
### 4.3 均值场博弈与图神经网络
#### 4.3.1 图卷积近似
**定义 4.4(图卷积均值场)** 在图结构数据上定义:
$$
\bar{a}_i = \frac{1}{\sqrt{d_i d_j}} \sum_{j \in \mathcal{N}_i} W_{ij} \cdot a_j
$$
其中 $W_{ij}$ 为邻接矩阵元素,$d_i$ 为节点度数。
#### 4.3.2 Graph Mean Field Q-Learning
**算法 3Graph Mean Field Q-Learning**
```
对于每个图卷积层 l
对每个节点 i
聚合邻居消息m_j^l = ReLU(W^l · h_j^{l-1})
消息聚合h_i^l = σ(W^l · [h_i^{l-1}, ∑_{j∈N_i} m_j^l]))
更新均值场μ_i^l = Mean({h_j^l : j ∈ N_i})
```
---
## 5. 通信与协作
### 5.1 通信协议学习
#### 5.1.1 CommNet 架构
**定义 5.1CommNet 通信层)** 第 $t$ 步的通信状态:
$$
c_i^t = \frac{1}{n-1} \sum_{j \neq i} \tanh(W \cdot h_j^{t-1} + b)
$$
$$
h_i^t = \tanh(U \cdot h_i^{t-1} + V \cdot c_i^t + b)
$$
其中 $h_i^t$ 为智能体 $i$ 的隐状态。
#### 5.1.2 Binary Communication (BIC)
**定义 5.2(瓶颈通信)** 使用自编码器结构强制信息压缩:
$$
\text{Encoder}: z_i = \text{Enc}(o_i, a_i^{t-1}) \in \mathbb{R}^k
$$
$$
\text{Decoder}: \hat{o}_i = \text{Dec}(z_i)
$$
$$
\text{BIC Loss} = \|o_i - \hat{o}_i\|^2 + \beta \|z_i\|^2
$$
#### 5.1.3 延迟通信建模
**定义 5.3(延迟通信信道)** 第 $i$ 个智能体在时刻 $t$ 接收到的消息:
$$
\tilde{c}_i^t = \sum_{\tau=0}^{L} \alpha_\tau c_i^{t-\tau}
$$
其中 $L$ 为最大延迟步数,$\alpha_\tau$ 为衰减系数。
### 5.2 协作激励机制
#### 5.2.1 内在奖励设计
**定义 5.4(内在动机)** 智能体 $i$ 的内在奖励:
$$
r_i^{\mathrm{KL} = r_i^{\mathrm{KL} + \beta \cdot \text{ICI}(i, \text{neighbors})
$$
其中 $\text{ICI}$Intrinsic Curiosity Incentive
$$
\text{ICI}_i = \| \hat{s}_{i}^{t+1} - s_{i}^{t+1} \|^2
$$
#### 5.2.2 teamGym 评估框架
**定义 5.5(团队奖励分解)** 在团队任务中:
$$
r^{\mathrm{KL} = \sum_i r_i^{\mathrm{KL} + \lambda \cdot \text{IGM}(s)
$$
其中 $\text{IGM}$Individual Global Maxim衡量因子化策略与联合最优的一致性。
---
## 6. 对抗多智能体
### 6.1 零和博弈中的多智能体
#### 6.1.1 极小极大 Q 学习
**定义 6.1(极小极大 Q 函数)** 对抗设置下智能体 $i$ 的 Q 函数:
$$
Q_i(s, \mathbf{a}) = \min_{\pi_{-i}} \max_{\pi_i} \mathbb{E}_{\pi} \left[ \sum_{t=0}^{\infty} \gamma^t r_i(s_t, \mathbf{a}_t) \right]
$$
#### 6.1.2 双人零和博弈的收敛性
**定理 6.1(极小极大定理)** 双人零和博弈存在**纳什均衡** $\sigma^*$,满足:
$$
\max_{\sigma_1} \min_{\sigma_2} u_1(\sigma_1, \sigma_2) = \min_{\sigma_2} \max_{\sigma_1} u_1(\sigma_1, \sigma_2) = u_1(\sigma^*)
$$
### 6.2 对抗性 RL
#### 6.2.1 对抗扰动攻击
**定义 6.2(对抗状态扰动)** 对手通过添加扰动 $\delta$ 修改观测:
$$
\tilde{o}_i = o_i + \delta, \quad \|\delta\| \leq \epsilon
$$
**定义 6.3(鲁棒策略)** 策略 $\pi_i$ 对抗 $\epsilon$-有界扰动的鲁棒性:
$$
\pi_i^{\mathrm{KL} = \arg\min_{\pi_i} \max_{\|\delta\| \leq \epsilon} \mathcal{L}(\pi_i, \tilde{o}_i)
$$
#### 6.2.2 PGD 对抗训练
**算法 4MADDPG + PGD 对抗训练**
```
for episode do
for 内部步 k = 1 to K do
生成对抗扰动 δ:
δ ← clip(δ + α · sign(∇_δ L(π, o + δ)), -ε, ε)
end for
执行带扰动的联合动作 (a_1 + δ_1, ..., a_n + δ_n)
使用 MADDPG 更新策略
end for
```
### 6.3 稳健多智能体策略
#### 6.3.1 分布鲁棒优化
**定义 6.4(分布鲁棒均衡)** 在分布不确定集 $\mathcal{P}$ 上优化:
$$
\max_{\pi_i} \min_{\mathbb{P} \in \mathcal{P}} \mathbb{E}_{\mathbb{P}} [V^{\pi_i}(\mathbb{P})]
$$
其中 $\mathcal{P}$ 可建模为 $\epsilon$-邻域的 KL 球:
$$
\mathcal{P} = \{ \mathbb{P} : D_{\mathrm{KL}}(\mathbb{P} \| \mathbb{P}_0) \leq \epsilon \}
$$
---
## 7. 收敛性与理论
### 7.1 Marlin 算法的收敛性分析
#### 7.1.1 Marlin 算法框架
Marlin (Schuurmans et al., 2023) 将多智能体值函数分解与策略迭代结合。
**定理 7.1Marlin 收敛性)** 在满足以下条件时Marlin 算法收敛至纳什均衡:
1. **单调性**:值函数估计单调递增
2. **有界性**$\|Q_i\| \leq B, \forall i$
3. **探索充分性**:每个智能体以概率 $p > 0$ 探索所有动作
*证明*:将 Marlin 视为近似的策略迭代过程。使用压缩映射定理,每次迭代的误差上界:
$$
\|Q^{k+1} - Q^*\| \leq \gamma \|Q^k - Q^*\| + \epsilon_{\mathrm{KL}
$$
其中 $\epsilon_{\mathrm{KL}}$ 为分解近似的误差。递归展开:
$$
\|Q^k - Q^*\| \leq \gamma^k \|Q^0 - Q^*\| + \frac{\epsilon_{\mathrm{KL}}{1-\gamma}
$$
当 $k \to \infty$$\gamma^k \|Q^0 - Q^*\| \to 0$,故收敛到 $Q^*$ 的 $\frac{\epsilon_{\mathrm{KL}}{1-\gamma}$ 邻域。$\square$
### 7.2 多智能体信用分配的收敛界
#### 7.2.1 有限样本界
**定理 7.2(信用分配有限样本界)** 设使用 Shapley 值进行信用分配,基于 $m$ 个样本估计,则以概率 $1-\delta$
$$
|\hat{\phi}_i - \phi_i| \leq \frac{C \sqrt{\ln(1/\delta)}}{\sqrt{m}}
$$
*证明*Shapley 值是 $n!$ 个排列的均值。使用 Hoeffding 不等式和 Union Bound 可得上式。$\square$
#### 7.2.2 偏差-方差 tradeoff
**定义 7.1(信用分配偏差)** 近似 Shapley 值与真实 Shapley 值的偏差:
$$
\text{Bias}(\hat{\phi}_i) = |\mathbb{E}[\hat{\phi}_i] - \phi_i|
$$
**定义 7.2(信用分配方差)** 估计的方差:
$$
\text{Var}(\hat{\phi}_i) = \mathbb{E}[(\hat{\phi}_i - \mathbb{E}[\hat{\phi}_i])^2]
$$
### 7.3 无中心化 Critic 时的稳定性
#### 7.3.1 独立 Q 学习的发散性
**定理 7.3(独立 Q 学习的非收敛性)** 在一般和博弈中,独立 Q 学习IQL**不保证**收敛到纳什均衡。
*反例*:考虑两人协调博弈:
| | $a_2^L$ | $a_2^R$ |
|-------|---------|---------|
| $a_1^L$ | 1, 1 | 0, 0 |
| $a_1^R$ | 0, 0 | 1, 1 |
IQL 可能收敛到 $(a_1^L, a_2^R)$ 或 $(a_1^R, a_2^L)$,而非帕累托最优的 $(L, L)$ 或 $(R, R)$。
#### 7.3.2 稳定化技术
**纳什 Q 学习** (Nash Q-Learning)
$$
Q_i^{t+1}(s, \mathbf{a}) = (1-\alpha) Q_i^t(s, \mathbf{a}) + \alpha [r_i(s, \mathbf{a}) + \gamma \text{NashQ}_i(s')]
$$
其中 $\text{NashQ}_i(s') = \min_{\pi_{-i}} \max_{\pi_i} \sum_{\mathbf{a}'} \pi_i(a_i') Q_i(s', \mathbf{a}')$。
#### 7.3.3 乐观/悲观策略
**定义 7.3(乐观策略)** 假设其他智能体会配合:
$$
\pi_i^{\mathrm{KL} = \arg\max_{\pi_i} \max_{\pi_{-i}} \mathbb{E}[r_i(\pi_i, \pi_{-i})]
$$
**定义 7.4(悲观策略)** 假设其他智能体会对抗:
$$
\pi_i^{\mathrm{KL} = \arg\max_{\pi_i} \min_{\pi_{-i}} \mathbb{E}[r_i(\pi_i, \pi_{-i})]
$$
---
## 8. 总结与前沿方向
### 8.1 方法对比
| 方法 | 分解结构 | 通信需求 | 收敛保证 | 适用范围 |
|------|----------|----------|----------|----------|
| VDN | 线性加和 | 无 | 仅团队游戏 | 小规模协作 |
| QMIX | 单调混合 | 无 | 仅团队游戏 | 中等规模协作 |
| QTRAN | 一般分解 | 无 | 部分 | 一般和博弈 |
| MADDPG | 中心化 Critic | 可选 | 无 | 协作/竞争 |
| MAPPO | 中心化 Critic | 可选 | 无 | 协作/竞争 |
| MFQ | 均值场近似 | 无 | $O(1/\sqrt{n})$ | 大规模系统 |
### 8.2 前沿研究方向
1. **鲁棒 MARL**:对抗扰动下的稳健策略学习
2. **通信可信性**:可验证的通信协议
3. **层次化 MARL**:多尺度协作与竞争
4. **离线 MARL**:基于历史数据的多智能体策略学习
5. **元 MARL**:快速适应新任务的多智能体系统
---
## 参考文献
1. Zhang, K., Yang, Z., & Başar, T. (2021). Multi-Agent Reinforcement Learning: A Selective Overview of Theories and Algorithms. *Handbook of RL for Control systems*.
2. Oliehoek, F. A., & Amato, C. (2016). *A Concise Introduction to Decentralized POMDPs*. Springer.
3. Sunehag, P., et al. (2017). Value Decomposition Networks. *AAAI*.
4. Rashid, T., et al. (2018). QMIX: Monotonic Value Function Factorisation for Deep Multi-Agent RL. *ICML*.
5. Lowe, R., et al. (2017). Multi-Agent Actor-Critic for Mixed Cooperative-Competitive Environments. *NeurIPS*.
6. Yang, Y., et al. (2020). Mean Field Multi-Agent Reinforcement Learning. *ICML*.
7. Foerster, J., et al. (2016). Learning to Communicate with Deep Multi-Agent Reinforcement Learning. *NeurIPS*.
8. Schuurmans, D., et al. (2023). Marlin: A Theoretical Framework for Multi-Agent RL. *arXiv*.