Policy Gradient,策略梯度法是强化学习中的一种常用方法。
比较详细的推导可以看:https://datawhalechina.github.io/easy-rl/#/chapter4/chapter4?id=tips
强化学习由三个组件构成:Actor、Environment、Reward Function。
Actor 决定如何执行下一个动作;
Environment 就是电脑主机或者对手会根据Actor的动作产生变化;
Reward Function 会根据我们做的动作进行打分。
在强化学习里面 Environment 和 Reward Function 都是不可操作的,我们的目标就是调整Actor的策略(Policy)获得最大的奖励期望。
一场游戏叫做一个 回合(episode) 或者 试验(trial);
过程就是Env先初始化一个状态 s 1 s_1 s1,而后我们的Actor根据该初始状态做出动作;
产生一个Trajectory,我们用 τ \tau τ表示,如下:
τ = { s 1 , a 1 , s 2 , a 2 , . . . , s t , a t } \tau = \{s_1, a_1, s_2, a_2, ..., s_t, a_t \} τ={s1,a1,s2,a2,...,st,at}
我们使用神经网络来完成动作的选择,以打电动游戏为例子:
我们有策略网络来实现策略 π \pi π,网络的参数为 θ \theta θ,输入为游戏的图像,输出的我们的操作。
在一场游戏中(一个episode),我们的Actor可以与环境交互产生一个回合的记录序列 Trajectory:
τ = { s 1 , a 1 , s 2 , a 2 , . . . , s t , a t } \tau = \{s_1, a_1, s_2, a_2, ..., s_t, a_t \} τ={s1,a1,s2,a2,...,st,at}
每一个 τ \tau τ产生的概率为:
p θ ( τ ) = p ( s 1 ) p θ ( a 1 ∣ s 1 ) p ( s 2 ∣ s 1 , a 1 ) p θ ( a 2 ∣ s 2 ) p ( s 3 ∣ s 2 , a 2 ) . . . = p ( s 1 ) ∏ t = 1 T p θ ( a t ∣ s t ) p ( s t + 1 ∣ s t , a t ) \begin{aligned} p_\theta(\tau) &= p(s_1)p_\theta(a_1|s_1)p(s_2|s_1, a_1)p_\theta(a_2|s_2)p(s_3|s_2,a_2)... \\ &=p(s_1)\prod_{t=1}^Tp_\theta(a_t|s_t)p(s_{t+1}|s_t, a_t) \end{aligned} pθ(τ)=p(s1)pθ(a1∣s1)p(s2∣s1,a1)pθ(a2∣s2)p(s3∣s2,a2)...=p(s1)t=1∏Tpθ(at∣st)p(st+1∣st,at)
每一个 τ \tau τ的总回报为:
R ( τ ) = ∑ t = 1 T r t R(\tau) = \sum_{t=1}^Tr_t R(τ)=t=1∑Trt
我们的 τ \tau τ是利用Actor与环境互动产生的,在动作选择过程中存在很多随机性,环境本身也存在很多随机性,所以 R ( τ ) R(\tau) R(τ)是一个随机变量(random variable)。
所以我们不能用 R ( τ ) R(\tau) R(τ)评价一个策略网络 θ \theta θ的好坏,但是我们可以使用 R ( τ ) R(\tau) R(τ)的期望评价:
E [ R ( τ ) ] = R ˉ ( τ ) = ∑ τ R ( τ ) p θ ( τ ) = E τ ~ p θ ( τ ) [ R ( τ ) ] \begin{aligned} E[R(\tau)] = \bar{R}(\tau) = \sum_\tau R(\tau)p_\theta(\tau) = E_{\tau~p_\theta(\tau) }[R(\tau)] \end{aligned} E[R(τ)]=Rˉ(τ)=τ∑R(τ)pθ(τ)=Eτ~pθ(τ)[R(τ)]
因此我们的目标就是最大化(maximize) R ˉ ( τ ) \bar{R}(\tau) Rˉ(τ),那么我们只需要朝着 R ˉ ( τ ) \bar{R}(\tau) Rˉ(τ)梯度上升的方向调整参数 θ \theta θ即可,
∇ R ˉ ( τ ) = ∑ τ R ( τ ) ∇ p θ ( τ ) = ∑ τ R ( τ ) p θ ( τ ) ∇ p θ ( τ ) p θ ( τ ) = ∑ τ R ( τ ) p θ ( τ ) ∇ l o g p θ ( τ ) = E τ ~ p θ ( τ ) [ R ( τ ) ∇ l o g p θ ( τ ) ] ≈ 1 N ∑ n = 1 N R ( τ ) ∇ l o g p θ ( τ ) = 1 N ∑ n = 1 N ∑ t = 1 T n R ( τ n ) ∇ l o g p θ ( a t n ∣ s t n ) \begin{aligned} \nabla\bar{R}(\tau) &= \sum_\tau R(\tau)\nabla p_\theta(\tau) \\ &= \sum_\tau R(\tau) p_\theta(\tau) \frac {\nabla p_\theta(\tau)}{p_\theta(\tau)} \\ &= \sum_\tau R(\tau) p_\theta(\tau) \nabla log p_\theta(\tau) \\ &= E_{\tau~p_\theta(\tau) }[R(\tau) \nabla log p_\theta(\tau) ] \\ & \approx \frac1N \sum_{n=1}^N R(\tau) \nabla log p_\theta(\tau) \\ &= \frac1N \sum_{n=1}^N \sum_{t=1}^{T_n} R(\tau^n) \nabla log p_\theta(a_t^n|s_t^n) \end{aligned} ∇Rˉ(τ)=τ∑R(τ)∇pθ(τ)=τ∑R(τ)pθ(τ)pθ(τ)∇pθ(τ)=τ∑R(τ)pθ(τ)∇logpθ(τ)=Eτ~pθ(τ)[R(τ)∇logpθ(τ)]≈N1n=1∑NR(τ)∇logpθ(τ)=N1n=1∑Nt=1∑TnR(τn)∇logpθ(atn∣stn)
直观理解这个梯度公式:我们知道在 s t s_t st状态下执行了 a t a_t at,如果最终的回报 R ( τ n ) R(\tau^n) R(τn)是好的,那么就增加 ( s t ∣ a t ) ( s_t| a_t) (st∣at)的选择概率,否则减少。
回顾一下我们的参数更新方法:
θ ← θ + η ∇ R ˉ ( τ ) \theta \larr \theta + \eta\nabla\bar{R}(\tau) θ←θ+η∇Rˉ(τ)
上面我们将梯度的计算方法调整为:
∇ R ˉ ( τ ) = 1 N ∑ n = 1 N ∑ t = 1 T n ( R ( τ n ) − b ) ∇ l o g p θ ( a t n ∣ s t n ) \nabla\bar{R}(\tau) = \frac1N \sum_{n=1}^N \sum_{t=1}^{T_n} (R(\tau^n) - b) \nabla log p_\theta(a_t^n|s_t^n) \\ ∇Rˉ(τ)=N1n=1∑Nt=1∑Tn(R(τn)−b)∇logpθ(atn∣stn)
从公式中我们可以看出,在同一个episode里面,所有的 ( s t n ∣ a t n ) ( s_t^n| a_t^n) (stn∣atn)都具有相同的调整权重 ( R ( τ n ) − b ) (R(\tau^n) - b) (R(τn)−b),这显然是不合理的:
一种解决的思路是:
我们继续观察现在的梯度公式,这里还存在一个问题,就是未来十分久远的决策的反馈并不一定受当前决策的影响:
所以最终的梯度公式如下:
∇ R ˉ ( τ ) = 1 N ∑ n = 1 N ∑ t = 1 T n ( ∑ t ′ = t T r t ′ n γ t ′ − t − b ) ∇ l o g p θ ( a t n ∣ s t n ) , γ ∈ [ 0 , 1 ] \nabla\bar{R}(\tau) = \frac1N \sum_{n=1}^N \sum_{t=1}^{T_n} (\sum_{t'=t}^T r_{t'}^n \gamma^{t'-t} - b) \nabla log p_\theta(a_t^n|s_t^n) , \gamma \in[0,1] ∇Rˉ(τ)=N1n=1∑Nt=1∑Tn(t′=t∑Trt′nγt′−t−b)∇logpθ(atn∣stn),γ∈[0,1]
我们将 R − b R-b R−b这一项定义为 Advantage function(比较优势) 符号表示为 A θ ( s t , a t ) A^{\theta}(s_t, a_t) Aθ(st,at):
A θ ( s t , a t ) = ∑ t ′ = t T r t ′ n γ t ′ − t − b A^{\theta}(s_t, a_t) = \sum_{t'=t}^T r_{t'}^n \gamma^{t'-t} - b Aθ(st,at)=t′=t∑Trt′nγt′−t−b