Policy Optimization

本文主要对基于策略优化定理的强化学习算法及相关变体进行一个梳理和总结。

Preliminaries

一个马尔可夫决策过程(Markov decision process, MDP)可以由五元组 <S,A,R,T,γ,ρ0> 定义,其中:

  • S 表示状态空间,是由状态构成的集合
  • A 表示动作空间,是由动作构成的集合
  • R:S×AR 表示奖励函数,R(s,a) 表示在状态 s 下执行动作 a 获得的奖励
  • T:S×A×S[0,1] 表示状态转移概率函数,T(s|s,a) 表示在状态 s 下执行动作 a 到达状态s 的概率
  • γ 表示折扣因子
  • ρ0:S[0,1] 表示状态初始分布

Agent 的决策过程由一个随机策略 π:S×A[0,1] 表示,π(a|s) 表示 agent 在状态 s 下选择动作 a 的概率。基于策略 π 可以定义出状态价值函数: (1)Vπ(s)=E[t=0γtR(st,at)|π,s] 状态动作价值函数: (2)Qπ(s,a)=R(s,a)+γEsT(|s,a)[Vπ(s)]

优势函数: (3)Aπ(s,a)=Qπ(s,a)Vπ(s) Agent 的优化目标通常为最大化在初始状态分布上的折扣汇报,即: (4)η(π)=Esρ0[Vπ(s)] 状态的折扣访问频率表示为: (5)ρπ(s)=(P(s0=s)+γP(s1=s)+γ2P(s2=s)+...)

策略梯度定理: (6)η(π)=s,aρπ(s)π(a|s)Qπ(s,a)

Approximately Optimal Approximate Reinforcement Learning

Kakade, S. & Langford, J. Approximately Optimal Approximate Reinforcement Learning. in Proceedings of the Nineteenth International Conference on Machine Learning 267–274 (Morgan Kaufmann Publishers Inc., 2002).

本文提出了三个想要回答的问题:

  1. 是否存在性能度量可以保证每一步更新都有提升?
  2. 验证某个更新提升该性能度量有多么困难?
  3. 在一定合理的次数的策略更新后,策略性能能达到什么样的水平?

考虑如下保守策略更新规则: (7)πnew(a|s)=(1α)π(a|s)+απ(a|s) 其中,α[0,1],如果要在 α=1 时保证策略提升,那么 π 需要在每个状态下都采取比 π 更好的动作。考虑 0<α<1 时,策略提升只需要 π 在大部分而非全部状态选择更好的动作。定义策略优势 Aπ,ρ0(π)(8)Aπ,ρ0(π)=Esρπ[Eaπ(|s)[Aπ(s,a)]] 该策略优势函数衡量了在状态初始分布为 ρ0,行动策略为 π 时,π 选择具有更大优势动作的程度。对于 α=0 时,显然有 ηα|α=0=11γAπ,ρ0

to be continue......

TRPO

Schulman, J., Levine, S., Moritz, P., Jordan, M. & Abbeel, P. Trust region policy optimization. in Proceedings of the 32nd International Conference on International Conference on Machine Learning - Volume 37 1889–1897 (JMLR.org, 2015).

策略优势定理: (9)η(π~)=η(π)+Es0ρ0,aπ~(|s),sT(|s,a)[t=1γtAπ(st,at)] 证明如下(为方便书写,省略期望下标中的 s0ρ0,sT(|s,a)aπ(|s) 简写为 aπ): (10)Eaπ~[t=0γtAπ(st,at)]=Eaπ~[t=0γt[Qπ(st,at)Vπ(st)]]=Eaπ~[t=0γt[rt+γVπ(st+1)Vπ(st)]]=Eaπ~[t=0γtrt]+Eaπ~[t=0γt+1Vπ(st+1)]Eaπ~[t=0γtVπ(st)]=η(π~)Eaπ~[Vπ(s0)]=η(π~)η(π)

根据式 5 ,式 9 表明也可以写为以下形式: (11)η(π~)=η(π)+sρπ~(s)aπ~(a|s)Aπ(s,a) 根据式 11,如果对于每一个状态 s,使 aπ~(a|s)Aπ(s,a)0,那么就保证了新策略是比旧策略好的。然而,式 11 对于 ρπ~ 的依赖导致在实际优化过程中无法直接计算,由此,考虑下面对 η 的局部估计: (12)Lπ(π~)=η(π)+sρπ(s)aπ~(a|s)Aπ(s,a) Lπ(π~)η(π~) 唯一的区别在于把 ρπ~ 换为了 ρπ,接下来对 Lπ(π~) 进行分析。若有一参数化的策略πθ 且对 θ 可微,则有: (13)Lπθ0(πθ0)=η(πθ0)θLπθ0(πθ)|θ=θ0=θη(πθ)|θ=θ0 因此,如果 π 更新的足够小,则提升 Lπ 也会提升 η,考虑如下策略更新: (14)π~(a|s)=(1α)π(a|s)+απ(a|s) 其中,π=argmaxπLπ(π~),从而有(证明见 (Kakade & Langford, 2002)): (15)η(π~)Lπ(π~)2ϵγ(1γ(1α))(1γ)α2 其中 ϵ=maxs|Eaπ~(|s)Aπ(s,a)|

α 足够小的时候则有: (16)η(π~)Lπ(π~)2ϵγ(1γ)2α2α=DTVmax(π,π~)ϵ=maxsmaxaAπ(s,a),则式 16 依然成立,证明见 (Schulman et al., 2015)。考虑到 DTV(p||q)2DKL(p||q),则有: (17)η(π~)Lπ(π~)CDKLmax(π,π~)Mi(π)=Lπi(π)CDKLmax(π,π~),则有: (18)Mi(πi)=η(πi)Mi(πi+1)η(πi+1)Mi(πi+1)Mi(πi)η(πi+1)η(πi) 因此,如果最大化 Mi(π),则 η(π) 也会随之不断提升,对于参数化的策略 πθ,策略优化问题便可以转换为下面的优化问题: (19)maxθ[Lθold(θ)CDKLmax(θold,θ)] 转换为信赖域约束的形式为: (20)maxθLθold(θ)(21)s.t.DKLmax(θold,θ)δ

上述优化稳定的约束很难进行技术员,可以使用平均 KL 散度作为对该约束启发式的估计: (22)D¯KLρ(θ1,θ2)=Esρ[DKL(πθ1(|s)||πθ2(|s))] 在采样估计时有如下形式: (23)maxθsρθoldaπθ(a|s)Aθold(s,a)(24)s.t.D¯KLρθold(θold,θ)δ 使用重要性采样(Importance Samping)可得: (25)aπθ(a|s)Aθold(s,a)=Eaπθold[πθ(a|s)πθold(a|s)Aθold(s,a)] 从而得到TRPO的优化目标: (26)maxθEsρθold, aπθold[πθ(a|s)πθold(a|s)Aθold(s,a)](27)s.t.Esρθold[DKL(πθold(|s)||πθ(|s))]

PPO

Schulman, J., Wolski, F., Dhariwal, P., Radford, A. & Klimov, O. Proximal Policy Optimization Algorithms. Preprint at http://arxiv.org/abs/1707.06347 (2017).

TRPO的优化目标为: (28)maxθEsρθold, aπθold[πθ(a|s)πθold(a|s)Aθold(s,a)](29)s.t.Esρθold[DKL(πθold(|s)||πθ(|s))] 将约束转换为一个惩罚项,即可将该优化问题转换为无约束优化 (30)maxθEsρθold, aπθold[πθ(a|s)πθold(a|s)Aθold(s,a)βDKL(πθold(|s)||πθ(|s))]rt(θ)=πθ(at|st)πθold(at|st),则rt(θold)=1,TRPO的近似优化目标为: (31)LCPI(θ)=Et[rt(θ)A^t] 在无约束的情况下,LCPI 可能会导致一次优化产生很大的策略更新,为此,约束 rt(θ) 在 1 附近来限制优化幅度: (32)LCLIP(θ)=Et[min(rt(θ)A^t,clip(rt(θ),1ϵ,1+ϵ)A^t]

IMPALA

Espeholt, L. et al. IMPALA: Scalable Distributed Deep-RL with Importance Weighted Actor-Learner Architectures. in Proceedings of the 35th International Conference on Machine Learning 1407–1416 (PMLR, 2018).

IMPALA 采用异步收集经验的方式来提高采样效率,使用 V-trace 解决采样策略和更新策略不一致的问题。

V-trace Target

(33)vt=V(st)+τ=tt+n1γτt(i=tτ1ci)δτV

其中,δτV=ρτ(rτ+γV(sτ+1)V(sτ))ρτ=min(ρ¯,π(aτ|sτ)μ(aτ|sτ))ci=min(c¯,π(ai|si)μ(ai|si)),规定当t=τ时,i=tτ1ci=1,此外ρ¯c¯,当π=μ时,ρτ=1ci=1,则上式转变为 on-policy n-step TD的形式: vt=V(st)+τ=tt+n1γτt(rτ+γV(sτ+1)V(sτ))(34)=τ=tt+n1γτtrτ+γnV(st+n)

V-trace 可以由如下递归形式计算: (35)vt=V(st)+δtV+γct(vt+1V(st+1)) 可以考虑给 ci 的计算中加入类似 Retrace(λ) 的折扣系数 λ[0,1],得到 ci=λmin(c¯,π(ai|si)μ(ai|si)),在 on-policy 场景下,当 n= 时,V-trace 变为 TD(λ)

Off-Policy TRPO

Meng, W., Zheng, Q., Shi, Y. & Pan, G. An Off-Policy Trust Region Policy Optimization Method With Monotonic Improvement Guarantee for Deep Reinforcement Learning. IEEE Transactions on Neural Networks and Learning Systems 33, 2223–2235 (2022).

策略优势定理为9,在 TRPO 中使用其 on-policy 的近似形式 12,两者的区别在于 TRPO 使用了近似状态分布 ρπ 替代策略分布 ρπ~。Off-policy TRPO 更进一步,采用行为策略 μ 来进行经验收集更新策略 π,从而提出新的近似形式: (36)Lπ,μ(π~)=η(π)+sρμ(s)aπ~(a|s)Aπ(s,a) 由该近似推出 off-policy 的优化目标: maxθEsρμ,aμ[πθ(a|s)μ(a|s)Aπθold(s,a)](37)s.t.D¯KLρμ,sqrt(μ,θold)DKLρμ,sqrt(θold,θ)+DKLρμ(θold,θ)δ 详细证明见(Meng et al., 2022)

Off-Policy PPO

Meng, W., Zheng, Q., Pan, G. & Yin, Y. Off-Policy Proximal Policy Optimization. Proceedings of the AAAI Conference on Artificial Intelligence 37, 9162–9170 (2023).

off-policy trpo 的 clip 近似形式,使用如下 clipped surrogate object: (38)LOff-PolicyCLIP=Esρμ,aμ[min[r(s,a)Aπold(s,a),clip(r(s,a),l(s,a),h(s,a)),A^πold(s,a)]] 其中 r(s,a)=π(a|s)μ(a|s)l(s,a)=πold(a|s)μ(a|s)(1ϵ)h(s,a)=πold(a|s)μ(a|s)(1+ϵ)

详见(Meng et al., 2023)

Behavior PPO

Zhuang, Z., Lei, K., Liu, J., Wang, D. & Guo, Y. Behavior Proximal Policy Optimization. in Proceedings of the eleventh International Conference on Learning Representation (2023).

在线同策略算法天然可以解决离线强化学习问题

该工作的架构为 BC(行为克隆)+RL(强化学习),RL部分通过保守更新和 ϵ 来实现在 offline dataset 上的多步更新。

根据 9 策略优势可写为: (39)JΔ(π,π^β)=Esρπ,aπ(|s)[Aπ^β(s,a)] 其中,π^β 为行为策略(behavior policy,可以从 dataset 通过行为克隆学习)

提出其在 offline dataset D 上的近似形式为: (40)J^Δ(π,π^β)=EsρD,aπ(|s)[Aπ^β(s,a)] 两者之间的差距为: JΔ(π,π^β)J^Δ(π,π^β)4γAπ^βmaxsDTV(π||π^β)[s]Esρπ^β()[DTV(π||π^β)[s]](41)2γAπ^βmaxsDTV(π||π^β)[s]EsρD()[1π^β(a|s)] 其中,Aπ^β=maxs,a|Aπ^β(s,a)|,详细证明见(Zhuang et al., 2023)。由式 41 可知,要保证优化目标 JΔ(π,π^β) 不下降,应该在最大化 EsρD,aπ(|s)[Aπ^β(s,a)] 的同时最小化 maxsDTV(π||π^β)[s],由此便得到了从 π^β 出发的单步更新。

考虑在离线数据集 D 上进行多步更新,即在 πk 的基础上优化得到 πk+1 ,根据策略优势定理,优化目标为: (42)JΔ(π,πk)=Esρπ,aπ(|s)[Aπk(s,a)]

考虑状态从离线数据集 D 进行采样的如下近似形式: (43)J^Δ(π,πk)=EsρD,aπ(|s)[Aπk(s,a)] 两者之间的差距为: JΔ(π,πk)J^Δ(π,πk)4γAkmaxsDTV(π||πk)[s]Esρπk()[DTV(π||πk)[s]]4γAkmaxsDTV(π||πk)[s]Esρπ^β()[DTV(πk||π^β)[s]](44)2γAπkmaxsDTV(π||πk)[s]EsρD()[1π^β(a|s)] 其中,Aπk=maxs,a|Aπk(s,a)|,详细证明见(Zhuang et al., 2023)。由式 44 可知,若保证 JΔ(π,πk) 非降,则需要在最大化 EsρD,aπ(|s)[Aπk(s,a)] 的同时最小化 AkmaxsDTV(π||πk)[s],从而实现在离线数据集 D 上的多步更新。

上述结论可写为: maxπEsρD,aπ(|s)[Aπk(s,a)]s.t.maxsDTV(π||πk)ϵ 经过一系列推导之后得到 clip 形式的无约束优化目标为: (45)Lk=EsρD,aπk(|s)[min(π(a|s)πk(a|s)Aπk(s,a),clip(π(a|s)πk(a|s),12ϵ,1+2ϵ)Aπk(s,a))]