▶《强化学习的数学原理》(2024春)_西湖大学赵世钰 Ch4 值迭代 与 策略迭代 【动态规划 算法】

news2024/12/23 10:53:00

PPT 截取必要信息。 课程网站做习题。总体 MOOC 过一遍

  • 1、视频 + 学堂在线 习题
  • 2、过 电子书 补充 【下载: 本章 PDF 电子书 GitHub】 [又看了一遍视频。原来第一次跳过了好多内容。。。]
  • 3、总体 MOOC 过一遍 习题

学堂在线 课程页面链接
中国大学MOOC 课程页面链接
B 站 视频链接

PPT和书籍下载网址: 【GitHub 链接】


在这里插入图片描述

总述:

开始介绍第一个可以找到最优策略的算法。 ——> 动态规划算法

介绍 3 种 迭代算法:
1、值迭代算法: 上一章讨论的求解 Bellman 最优方程的压缩映射定理 所提出的算法。
2、策略迭代算法
3、截断策略迭代算法: 值迭代 和 策略迭代 算法是该算法的极端情况。

动态规划 算法,需要系统模型。
本章介绍的策略迭代算法 扩展得到 第 5 章介绍的蒙特卡洛算法。
——————————————
model-based 算法

值迭代 上一章 的延伸
策略迭代 下一章 蒙特卡洛学习的基础

在这里插入图片描述

值迭代 和 策略迭代 是 截断策略迭代 的两个极端情况

4.1 值迭代

贝尔曼最优公式 的 矩阵向量形式:

v = f ( v ) = max ⁡ π ( r π + γ P π v ) \bm v=f(\bm v) =\max\limits_\pi({\bm r}_\pi+\gamma {\bm P}_\pi {\bm v}) v=f(v)=πmax(rπ+γPπv)

求解方法: 上一章 的 压缩映射定理 建议的迭代算法 【值迭代】

v k + 1 = f ( v k ) = max ⁡ π ( r π + γ P π v k ) ,     k = 1 , 2 , 3... {\bm v}_{k+1} = f({\bm v}_k)=\max\limits_\pi({\bm r}_\pi+\gamma {\bm P}_\pi {\bm v}_k), ~~~k=1, 2, 3... vk+1=f(vk)=πmax(rπ+γPπvk),   k=1,2,3...

其中 v 0 {\bm v}_0 v0 可为任意值。

两步:
1、策略 更新 (policy update)

  • v k {\bm v}_k vk 给定, 求解 π k + 1 = arg ⁡ max ⁡ π ( r π + γ P π v k ) \pi_{k+1} = \arg \max\limits_{\pi}({\bm r}_\pi+\gamma {\bm P}_\pi {\bm v}_k) πk+1=argπmax(rπ+γPπvk)

2、值 更新 (value update)

  • 上一步得到的策略 π k + 1 \pi_{k+1} πk+1, 更新 v k + 1 = r π k + 1 + γ P π k + 1 v k {\bm v}_{k +1}={\bm r}_{\pi_{k+1}}+\gamma {\bm P}_{\pi_{k+1}}{\bm v}_k vk+1=rπk+1+γPπk+1vk

在这里插入图片描述

v k v_k vk 是否是一个状态值?
答案是否定的。虽然 v k v_k vk 最终收敛于最优状态值,但不能保证满足任何策略的 Bellman方程。例如,它一般不满足 v k = r π k + γ P π k v k v_k=r_{\pi_k}+\gamma P_{\pi_k}v_k vk=rπk+γPπkvk v k = r π k + 1 + γ P π k + 1 v k v_k=r_{\pi_{k+1}}+\gamma P_{\pi_{k+1}}v_k vk=rπk+1+γPπk+1vk。它只是算法生成的一个中间值。另外,由于 v k v_k vk 不是状态值,所以 q k q_k qk 不是动作值。

编程实现 需要知道 更具体的形式 elementwise form
更新策略 更具体的形式为:
π k + 1 ( s ) = arg ⁡ max ⁡ π ∑ a π ( a ∣ s ) ( ∑ r p ( r ∣ s , a ) r + γ ∑ s ′ p ( s ′ ∣ s , a ) v k ( s ′ ) ) ⏟ q k ( s , a ) , s ∈ S \pi_{k+1}(s) = \arg \max\limits_{\pi}\sum_a\pi(a|s)\underbrace{\Big(\sum_rp(r|s, a)r+\gamma \sum_{s^{\prime}}p(s^{\prime}|s, a)v_k(s^{\prime})\Big)}_{q_k(s, a)},s \in S πk+1(s)=argπmaxaπ(as)qk(s,a) (rp(rs,a)r+γsp(ss,a)vk(s)),sS

根据上一章的分析,上述优化问题的最优策略解为:
π k + 1 ( a ∣ s ) = { 1 a = a k ∗ ( s ) 0 a ≠ a k ∗ ( s ) \pi_{k+1}(a|s)=\begin{cases}1\quad a=a_k^*(s) \\ 0\quad a\neq a_k^*(s) \end{cases} πk+1(as)={1a=ak(s)0a=ak(s)
其中 a k ∗ ( s ) = arg ⁡ max ⁡ a q k ( a , s ) a_k^*(s)=\arg\max\limits_aq_k(a, s) ak(s)=argamaxqk(a,s)

  • 如果 a k ∗ ( s ) = arg ⁡ max ⁡ a q k ( a , s ) a_k^*(s)=\arg\max\limits_aq_k(a, s) ak(s)=argamaxqk(a,s) 有多个解,我们可以选择任意一个解而不影响算法的收敛性。 只是若是没选中最优策略该选的动作,需要多迭代几次后才能获得最终的最优策略

贪心策略 π k + 1 \pi_{k+1} πk+1:贪心地选择 q k ( a , s ) q_k(a, s) qk(a,s) 最大的 action     ~~~     确定的策略
————————————
更新值 更具体的形式为:
v k + 1 ( s ) = ∑ a π k + 1 ( a ∣ s ) ( ∑ r p ( r ∣ s , a ) r + γ ∑ s ′ p ( s ′ ∣ s , a ) v k ( s ′ ) ) ⏟ q k ( s , a ) , s ∈ S = max ⁡ a q k ( a , s ) \begin{align*}v_{k +1}(s)&=\sum_a\pi_{k+1}(a|s)\underbrace{\Big(\sum_rp(r|s, a)r+\gamma \sum_{s^{\prime}}p(s^{\prime}|s, a)v_k(s^{\prime})\Big)}_{q_k(s, a)},s \in S\\ &=\max_aq_k(a, s)\end{align*} vk+1(s)=aπk+1(as)qk(s,a) (rp(rs,a)r+γsp(ss,a)vk(s)),sS=amaxqk(a,s)

迭代流程

v k ( s ) → q k ( s , a ) → v_k(s)\to q_k(s, a)\to vk(s)qk(s,a) 贪心策略   π k + 1 ( a ∣ s ) → ~\pi_{k+1}(a|s)\to  πk+1(as) 新的值   v k + 1 = max ⁡ a q k ( s , a ) ~v_{k+1}=\max\limits_{a}q_k(s, a)  vk+1=amaxqk(s,a)

伪代码: 值迭代算法
目标: 搜索 求解 贝尔曼最优公式的 最优状态值 和 最优策略。

在这里插入图片描述

遍历 每个状态 中的 每个动作, 计算 q k q_k qk

  • 策略 更新: 选择 q k q_k qk 最大的 action
  • 值 更新: 将 v k + 1 ( s ) v_{k+1}(s) vk+1(s) 更新为 计算得到的最大 q k q_k qk

4.1.2 例子

在这里插入图片描述

对 每个状态 的 每个动作 ,初始化 q q q 值表

在这里插入图片描述

在这里插入图片描述
按照这里
策略更新 是将 每个状态 的 q q q 值最大的动作 的选取概率 π ( a ∣ s ) \pi(a|s) π(as) 置为 1。     ~~~    等效于 让策略在这一步 做这个 q q q 值最大的动作
值更新 是将 每个状态 的 值更新为 相应状态的最大 q q q 值。

v 0 v_0 v0 可以任意选取,这里选择为 0。 不同的初值选取对迭代过程影响多大?如何根据具体情况选取合适的初值?
——> 比较直觉的是若是初始值选得离最优状态值较远, 需要的迭代次数会多些。

对于 状态 s 1 s_1 s1,动作 a 3 a_3 a3 a 5 a_5 a5 对应的 q q q 都是最大的, 这里直接选了 a 5 a_5 a5, 有没有可能在这里选 a 3 a_3 a3 得到的才是最优策略呢?
——> 确实有可能, 所以要多次迭代,收敛后迭代结束获得的就是 最优策略。

在这里插入图片描述

第一次 迭代, s1 没有达到 最优。

在这里插入图片描述

在这里插入图片描述

这里 迭代两次 就获得了 最优策略。

其它 更复杂情况 的迭代停止条件为:

在这里插入图片描述

迭代停止 则认为 获得了 最优策略。

4.2 策略迭代

主要内容: 是什么?——> 性质 ——> 如何 编程实现

任意给定的初始策略 π 0 \pi_0 π0

两步:
1、策略 评估 (policy evaluation, PE)

  • 计算 π k \pi_k πk 的状态值: v π k = r π k + γ P π k v π k         {\bm v}_{\pi_k}={\bm r}_{\pi_k}+\gamma {\bm P}_{\pi_k}{\bm v}_{\pi_k}~~~~~~~ vπk=rπk+γPπkvπk       求解 贝尔曼方程

策略评估做的事:通过计算相应的 状态值 来评估给定策略。

2、策略 优化 (policy improvement,PI)

  • π k + 1 = arg ⁡ max ⁡ π ( r π + γ P π v π k ) \pi_{k+1}=\arg\max\limits_\pi({\bm r}_\pi+\gamma {\bm P}_\pi {\bm v}_{\pi_k}) πk+1=argπmax(rπ+γPπvπk)

迭代流程

π 0 → P E v π 0 → P I π 1 → P E v π 1 → P I π 2 → P E v π 2 → P I . . . \pi_0\xrightarrow{PE}v_{\pi_0}\xrightarrow{PI}\pi_1\xrightarrow{PE}v_{\pi_1}\xrightarrow{PI}\pi_2\xrightarrow{PE}v_{\pi_2}\xrightarrow{PI}... π0PE vπ0PI π1PE vπ1PI π2PE vπ2PI ...

PE: 策略 评估
PI:策略 优化

现在处理以下几个问题:

Q1:在策略评估步骤中,如何通过求解 Bellman 方程得到状态值?
Q2:在策略优化步骤中,为什么新策略 π k + 1 \pi_{k+1} πk+1 优于 π k π_k πk?
Q3:为什么这样的迭代算法最终可以达到最优策略?
Q4:这个策略迭代算法和之前的值迭代算法是什么关系?

Q1:在策略评估步骤中,如何通过求解 Bellman 方程得到状态值?

如何 获取 v π k v_{\pi_k} vπk

已知: v π k = r π k + γ P π k v π k {\bm v}_{\pi_k}={\bm r}_{\pi_k}+\gamma {\bm P}_{\pi_k}{\bm v}_{\pi_k} vπk=rπk+γPπkvπk

方法一: 矩阵求逆

v π k = ( I − γ P π k ) − 1 r π k {\bm v}_{\pi_k}=({\bm I}-\gamma {\bm P}_{\pi_k})^{-1}{\bm r}_{\pi_k} vπk=(IγPπk)1rπk

方法二: 迭代 ✔

v π k ( j + 1 ) = r π k + γ P π k v π k ( j ) ,     j = 0 , 1 , 2 , . . . {\bm v}_{\pi_k}^{(j+1)}={\bm r}_{\pi_k}+\gamma {\bm P}_{\pi_k}{\bm v}_{\pi_k}^{(j)}, ~~~j=0,1,2,... vπk(j+1)=rπk+γPπkvπk(j),   j=0,1,2,...

策略迭代 是在策略评估步骤中嵌入另一个迭代算法的迭代算法!

Q2:在策略优化步骤中,为什么新策略 π k + 1 \pi_{k+1} πk+1 优于 π k π_k πk?

在这里插入图片描述

  • 证明 1: 在策略优化步骤中,为什么新策略 π k + 1 \pi_{k+1} πk+1 优于 π k π_k πk?    ~~    P73-

Q3:为什么策略迭代算法最终可以找到最优策略?

由于每次迭代都会改进策略, 即

v π 0 ≤ v π 1 ≤ v π 2 ≤ ⋯ ≤ v π k ≤ ⋯ ≤ v ∗ \bm v_{\pi_0}\leq\bm v_{\pi_1}\leq\bm v_{\pi_2}\leq\cdots\leq\bm v_{\pi_k}\leq\cdots\leq\bm v^* vπ0vπ1vπ2vπkv

v π k \bm v_{\pi_k} vπk 不断减小并最终收敛。仍需证明 将收敛到 v ∗ \bm v^* v

在这里插入图片描述

定理 4.1 (策略迭代的收敛性)。策略迭代算法生成的状态值序列 { v π k } k = 0 ∞ \{v_{\pi_k}\}_{k=0}^\infty {vπk}k=0 收敛到最优状态值 v ∗ v^* v。因此,策略序列 { π k } k = 0 ∞ \{\pi_k\}_{k=0}^\infty {πk}k=0 收敛到最优策略。

  • 证明 2: 证明策略迭代会收敛到 最优策略 P75

证明的思路是证明 策略迭代算法 比 值迭代算法 收敛得更快。

如果 策略迭代 和 值迭代 从相同的初始猜测开始,由于 策略迭代 算法的收敛性,策略迭代 将比 值迭代 收敛得更快

Q4:这个策略迭代算法和之前的值迭代算法是什么关系?

值迭代 和 策略迭代 是 截断策略迭代 的两个极端, 后续将进一步说明。

——————————————————
如何 实现 策略迭代算法?

在这里插入图片描述
在这里插入图片描述

策略迭代 算法:
目标: 搜索 最优状态值 和 最优策略

在这里插入图片描述

策略迭代算法 生成的中间值是是 状态值。 因为这些值是当前策略的 Bellman 方程的解。

4.2.3 例子

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
一个 示例 P79
发现一个有趣的现象:接近目标的状态 的策略 先变好, 远离目标的状态的策略会后变好。

在某一个状态, 选择 greedy action 时, 严重依赖于 其它状态的策略。
若其它状态的策略是不好的, 此时虽然选一个 动作值 ( q q q) 最大的 动作, 可能意义不大;
如果 其它状态 有能够到达目标区域 的策略, 选择变到那个状态,也能到达目标区域, 得到正的 reward。

当某个状态周围 没有状态 能够到达 目标区域 的时候, 这个状态无法到达目标区域。
当周围有状态能够到达目标区域策略时, 新的策略也能到达目标区域。

1、观察策略是如何演变的,一个有趣的模式是靠近目标区域的状态比远离目标区域的状态更早找到最优策略。只有较近的状态能先找到到达目标的轨迹,较远的状态才能找到经过较近状态到达目标的轨迹
2、状态值的空间分布呈现出一种有趣的模式:靠近目标的状态具有更大的状态值。这种模式的原因是,一个 agent 从更远的状态出发,必须走很多步才能获得正的奖励。这样的奖励将严重打折扣,因此相对较小。

4.3 截断策略迭代

值迭代 算法和 策略迭代 算法是截断策略迭代算法的两种特殊情况。

策略迭代:初始策略 为 π 0    \pi_0~~ π0   【任意猜测的】

  • 策略评估 (PE): v π k = r π k + γ P π k v π k {\bm v}_{\pi_k}={\bm r}_{\pi_k}+\gamma {\bm P}_{\pi_k}{\bm v}_{\pi_k} vπk=rπk+γPπkvπk
  • 策略优化 (PI): π k + 1 = arg ⁡ max ⁡ π ( r π + γ P π v π k ) \pi_{k+1}=\arg\max\limits_\pi({\bm r}_\pi+\gamma {\bm P}_\pi {\bm v}_{\pi_k}) πk+1=argπmax(rπ+γPπvπk)

值迭代:初始值 为 v 0 {\bm v}_0 v0

  • 策略更新 (PU): π k + 1 = arg ⁡ max ⁡ π ( r π + γ P π v k ) \pi_{k+1}=\arg\max\limits_\pi({\bm r}_\pi+\gamma {\bm P}_\pi {\bm v}_k) πk+1=argπmax(rπ+γPπvk)
  • 值更新 (VU): v k + 1 = r π k + 1 + γ P π k + 1 v k {\bm v}_{k+1}={\bm r}_{\pi_{k+1}}+\gamma {\bm P}_{\pi_{k+1}}{\bm v}_k vk+1=rπk+1+γPπk+1vk

!!每一步的等号右侧 都有的: r + γ P v {\bm r} +\gamma {\bm P} {\bm v} r+γPv

在这里插入图片描述

在这里插入图片描述

从相同的初始条件开始。
前三个步骤是相同的。
第四步就不一样了:

  • 策略迭代,求解 v π 1 = r π 1 + γ P π 1 v π 1 v_{π_1} = r_{π_1} + γP_{\pi_1}v_{\pi_1} vπ1=rπ1+γPπ1vπ1 需要一个迭代算法 ( 迭代无数次 )
  • 值迭代, v 1 = r π 1 + γ P π 1 v 0 v_1 = r_{π_1} + \gamma P_{π_1}v_0 v1=rπ1+γPπ1v0一步迭代

在这里插入图片描述

每步求解 v \bm v v 值时, 值迭代 需要一步, 策略迭代需要无穷步,迭代次数取中间值如何呢?

值迭代算法:计算一次。
策略迭代算法:计算无限次迭代。
截断策略迭代算法:计算一个有限次迭代(例如 j j j )。从 j j j ∞ \infty 的其余迭代被截断。

在这里插入图片描述

算法中的 v k v_k vk v k ( j ) v_k^{(j)} vk(j) 不是状态值,是真实状态值的近似值,因为在策略评估步骤中只执行有限次迭代。

只有当我们在 策略评估 步骤中运行无限次迭代时,才能获得真实的状态值。

截断策略迭代 会不会 结束迭代时是一个 发散的结果?

在这里插入图片描述

  • 证明。参考 电子书 PDF P83

——————————————

证明: 截断策略迭代算法 的收敛性。

因为

v π k ( j ) = r π k + γ P π k v π k ( j − 1 ) v_{\pi_k}^{(j)}=r_{\pi_k}+\gamma P_{\pi_k}v_{\pi_k}^{(j-1)} vπk(j)=rπk+γPπkvπk(j1)

v π k ( j + 1 ) = r π k + γ P π k v π k ( j ) v_{\pi_k}^{(j+1)}=r_{\pi_k}+\gamma P_{\pi_k}v_{\pi_k}^{(j)} vπk(j+1)=rπk+γPπkvπk(j)

v π k ( j + 1 ) − v π k ( j ) = γ P π k ( v π k ( j ) − v π k ( j − 1 ) ) = ⋯ = γ j P π k j ( v π k ( 1 ) − v π k ( 0 ) ) v_{\pi_k}^{(j+1)}-v_{\pi_k}^{(j)}=\gamma P_{\pi_k}(v_{\pi_k}^{(j)}-v_{\pi_k}^{(j-1)})=\cdots=\gamma^j P^j_{\pi_k}(v_{\pi_k}^{(1)}-v_{\pi_k}^{(0)}) vπk(j+1)vπk(j)=γPπk(vπk(j)vπk(j1))==γjPπkj(vπk(1)vπk(0))

v π k ( 0 ) = v π k − 1      v_{\pi_k}^{(0)}=v_{\pi_{k-1}}~~~~ vπk(0)=vπk1     上一轮迭代的结果

v π k ( 1 ) = r π k + γ P π k v π k ( 0 ) = r π k + γ P π k v π k − 1 ≥ r π k − 1 + γ P π k − 1 v π k − 1          ① = v π k − 1 = v π k ( 0 ) \begin{aligned}v_{\pi_k}^{(1)}&=r_{\pi_k}+\gamma P_{\pi_k}v_{\pi_k}^{(0)}\\ &=r_{\pi_k}+\gamma P_{\pi_k}\textcolor{blue}{v_{\pi_{k-1}}}\\ &\geq r_{\pi_{\textcolor{blue}{{k-1}}}}+\gamma P_{\pi_{\textcolor{blue}{{k-1}}}}\textcolor{blue}{v_{\pi_{k-1}}}~~~~~~~~\textcolor{blue}{①}\\ &=v_{\pi_{k-1}}\\ &=v_{\pi_k}^{(0)}\end{aligned} vπk(1)=rπk+γPπkvπk(0)=rπk+γPπkvπk1rπk1+γPπk1vπk1        =vπk1=vπk(0)

v π k ( j + 1 ) ≥ v π k ( j ) v_{\pi_k}^{(j+1)}\geq v_{\pi_k}^{(j)} vπk(j+1)vπk(j)

π k = arg ⁡ max ⁡ π ( r π + γ P π v π k − 1 ) \pi_k=\arg\max\limits_\pi(r_\pi+\gamma P_\pi v_{\pi_{k-1}}) πk=argπmax(rπ+γPπvπk1)

——————————————
在这里插入图片描述
相比于策略迭代算法,截断的策略迭代算法在策略评估步骤中只需要有限次数的迭代,因此计算效率更高。与值迭代相比,截断策略迭代算法可以在策略评估步骤中多运行几次迭代,从而加快收敛速度

Pl 【策略迭代】 的收敛性证明是基于 VI 【值迭代】 的收敛性证明。由于 VI 收敛,得到 PI 收敛。

小结:

在这里插入图片描述
4.5
Q:值迭代算法一定能找到最优策略吗?
是的。值迭代正是上一章求解 Bellman 最优性方程的 压缩映射定理 所提出的算法。利用 压缩映射定理 保证了算法的收敛性。

model-based VS model-free
虽然本章介绍的算法可以找到最优策略,但由于它们需要系统模型,通常被称为动态规划算法而不是强化学习算法。
强化学习算法可以分为两类:基于模型的和免模型的。
这里,“基于模型的”并不是指系统模型的需求。相反,基于模型的强化学习使用数据来估计系统模型,并在学习过程中使用该模型。相比之下,免模型强化学习在学习过程中不涉及模型估计

——————
习题

值迭代、策略迭代、截断策略迭代

值迭代算法中间产生的值不一定对应某些策略的状态值,这些只是产生的一些中间过程的数值,没有特别的含义。

压缩映射定理给出的算法 实际是 值迭代算法。

策略迭代算法 同时获得 最优状态值 和 最优策略。 【策略评估 需要计算状态值】

补充

证明 1: 在策略优化步骤中,为什么新策略 π k + 1 \pi_{k+1} πk+1 优于 π k π_k πk?    ~~    P73-

证明:
状态值 v π k + 1 v_{\pi_{k+1}} vπk+1 v π k v_{\pi_k} vπk 满足贝尔曼公式:
v π k + 1 = r π k + 1 + γ P π k + 1 v π k + 1 v_{\pi_{k+1}}=r_{\pi_{k+1}}+\gamma P_{\pi_{k+1}}v_{\pi_{k+1}} vπk+1=rπk+1+γPπk+1vπk+1
v π k = r π k + γ P π k v π k v_{\pi_k}=r_{\pi_k}+\gamma P_{\pi_k}v_{\pi_k} vπk=rπk+γPπkvπk
由于 π k + 1 = arg ⁡ max ⁡ π ( r π + γ P π v π k ) \pi_{k+1} = \arg\max\limits_\pi(r_\pi+\gamma P_\pi v_{\pi_k}) πk+1=argπmax(rπ+γPπvπk)
r π k + 1 + γ P π k + 1 v π k + 1 ≥ r π k + γ P π k v π k r_{\pi_{k+1}}+\gamma P_{\pi_{k+1}}v_{\pi_{k+1}}\geq r_{\pi_k}+\gamma P_{\pi_k}v_{\pi_k} rπk+1+γPπk+1vπk+1rπk+γPπkvπk
v π k − v π k + 1 = ( r π k + γ P π k v π k ) − ( r π k + 1 + γ P π k + 1 v π k + 1 ) ≤ ( r π k + 1 + γ P π k + 1 v π k ) − ( r π k + 1 + γ P π k + 1 v π k + 1 ) ≤ γ P π k + 1 ( v π k − v π k + 1 ) ≤ γ 2 P π k + 1 2 ( v π k − v π k + 1 ) ≤ . . . ≤ γ n P π k + 1 n ( v π k − v π k + 1 ) ≤ lim ⁡ n → ∞ γ n P π k + 1 n ( v π k − v π k + 1 ) = 0 \begin{align*}v_{\pi_k}-v_{\pi_{k+1}}&= (r_{\pi_k}+\gamma P_{\pi_k}v_{\pi_k}) - (r_{\pi_{k+1}}+\gamma P_{\pi_{k+1}}v_{\pi_{k+1}})\\ &\leq(r_{\pi_{k+1}}+\gamma P_{\pi_{k+1}}v_{\pi_k}) - (r_{\pi_{k+1}}+\gamma P_{\pi_{k+1}}v_{\pi_{k+1}})\\ &\leq \gamma P_{\pi_{k+1}}(v_{\pi_k} -v_{\pi_{k+1}})\\ &\leq \gamma^2 P_{\pi_{k+1}}^2(v_{\pi_k} -v_{\pi_{k+1}})\\ &\leq ...\\ &\leq \gamma^n P_{\pi_{k+1}}^n(v_{\pi_k} -v_{\pi_{k+1}})\\ &\leq \lim\limits_{n\to\infty}\gamma^n P_{\pi_{k+1}}^n(v_{\pi_k} -v_{\pi_{k+1}})\\ &=0\end{align*} vπkvπk+1=(rπk+γPπkvπk)(rπk+1+γPπk+1vπk+1)(rπk+1+γPπk+1vπk)(rπk+1+γPπk+1vπk+1)γPπk+1(vπkvπk+1)γ2Pπk+12(vπkvπk+1)...γnPπk+1n(vπkvπk+1)nlimγnPπk+1n(vπkvπk+1)=0
在这里插入图片描述

证明 2: 证明策略迭代会收敛到 最优策略 P75

在这里插入图片描述

定理 4.1 (策略迭代的收敛性)。策略迭代算法生成的状态值序列 { v π k } k = 0 ∞ \{v_{\pi_k}\}_{k=0}^\infty {vπk}k=0 收敛到最优状态值 v ∗ v^* v。因此,策略序列 { π k } k = 0 ∞ \{\pi_k\}_{k=0}^\infty {πk}k=0 收敛到最优策略。

证明的思路是证明 策略迭代算法 比 值迭代算法 收敛得更快。

——————————
证明:
为了 证明 { v π k } k = 0 ∞ \{v_{\pi_k}\}_{k=0}^\infty {vπk}k=0 的收敛性, 引入由以下式子生成的 另一个序列 { v k } k = 0 ∞ \{v_k\}_{k=0}^\infty {vk}k=0

v k + 1 = f ( v k ) = max ⁡ π ( r π + γ P π v k ) v_{k+1}=f(v_k)=\max\limits_\pi(r_\pi+\gamma P_\pi v_k) vk+1=f(vk)=πmax(rπ+γPπvk)

这个迭代算法 正是 值迭代算法,则给定任意初始值 v 0 v_0 v0, v k v_k vk 收敛到 v ∗ v^* v

k = 1 k=1 k=1, 对任意 π 0 \pi_0 π0, 有 v π 0 ≥ v 0 v_{\pi_0}\geq v_0 vπ0v0

通过 归纳法 证明 对任意 k k k, 有 v k ≤ v π k ≤ v ∗ v_k\leq v_{\pi_k}\leq v^* vkvπkv

k ≥ 0 k\geq0 k0, 假设 v π k ≥ v k v_{\pi_k}\geq v_k vπkvk

用到的一些中间式:

在这里插入图片描述
v π k + 1 ≥ v π k    v_{\pi_{k+1}}\geq v_{\pi_k}~~ vπk+1vπk   【上面的 证明 1 已证。即 策略优化后的策略的状态值 比之前的大】 , P π k + 1 ≥ 0 P_{\pi_{k+1}}\geq0 Pπk+10
② 令 π k ′ = arg ⁡ max ⁡ π ( r π + γ P π v k ) {\textcolor{blue}{{\pi_k^\prime}}}=\arg \max\limits_\pi(r_\pi+\gamma P_\pi v_k) πk=argπmax(rπ+γPπvk)
π k + 1 = arg ⁡ max ⁡ π ( r π + γ P π v π k ) \pi_{k+1}=\arg \max\limits_\pi(r_\pi+\gamma P_\pi v_{\pi_k}) πk+1=argπmax(rπ+γPπvπk)

对于 k + 1 k + 1 k+1 有:

v π k + 1 − v k + 1 = ( r π k + 1 + γ P π k + 1 v π k + 1 ) − max ⁡ π ( r π + γ P π v k ) ≥ ( r π k + 1 + γ P π k + 1 v π k ) − max ⁡ π ( r π + γ P π v k )            ① = ( r π k + 1 + γ P π k + 1 v π k ) − ( r π k ′ + γ P π k ′ v k )            ② ≥ ( r π k ′ + γ P π k ′ v π k ) − ( r π k ′ + γ P π k ′ v k )            ③ = γ P π k ′ ( v π k − v k ) \begin{aligned}v_{\pi_{k+1}}-v_{k+1}&=(r_{\pi_{k+1}}+\gamma P_{\pi_{k+1}}v_{\pi_{k+1}})-\max\limits_\pi(r_\pi+\gamma P_\pi v_k)\\ &\geq(r_{\pi_{k+1}}+\gamma P_{\pi_{k+1}}v_{\pi_{\textcolor{blue}{k}} })-\max\limits_\pi(r_\pi+\gamma P_\pi v_k)~~~~~~~~~~\textcolor{blue}{①}\\ &=(r_{\pi_{k+1}}+\gamma P_{\pi_{k+1}}v_{\pi_k })-(r_{\textcolor{blue}{{\pi_k^\prime}}}+\gamma P_{\textcolor{blue}{{\pi_k^\prime}}}v_k)~~~~~~~~~~\textcolor{blue}{②}\\ &\geq(r_{\textcolor{blue}{{\pi_k^\prime}}}+\gamma P_{\textcolor{blue}{{\pi_k^\prime}}}v_{\pi_k })-(r_{\textcolor{blue}{{\pi_k^\prime}}}+\gamma P_{\textcolor{blue}{{\pi_k^\prime}}}v_k)~~~~~~~~~~\textcolor{blue}{③}\\ &=\gamma P_{\pi_k^\prime}(v_{\pi_k}-v_k)\end{aligned} vπk+1vk+1=(rπk+1+γPπk+1vπk+1)πmax(rπ+γPπvk)(rπk+1+γPπk+1vπk)πmax(rπ+γPπvk)          =(rπk+1+γPπk+1vπk)(rπk+γPπkvk)          (rπk+γPπkvπk)(rπk+γPπkvk)          =γPπk(vπkvk)

因为 v π k − v k ≥ 0 v_{\pi_k}-v_k\geq0 vπkvk0 P π k ′ P_{\pi_k^\prime} Pπk 非负。

γ P π k ′ ( v π k − v k ) ≥ 0 \gamma P_{\pi_k^\prime}(v_{\pi_k}-v_k)\geq0 γPπk(vπkvk)0

v π k + 1 − v k + 1 ≥ 0 v_{\pi_{k+1}}-v_{k+1}\geq0 vπk+1vk+10

归纳得到, 对任意 k > 0 k > 0 k>0 v k ≤ v π k ≤ v ∗ v_k\leq v_{\pi_k}\leq v^* vkvπkv
v k v_k vk 收敛到 v ∗ v^* v, 由夹逼准则可得, v π k v_{\pi_k} vπk 也收敛到 v ∗ v^* v

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1840726.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

【宠粉赠书】科技图表绘制:R语言数据可视化

为了回馈粉丝们的厚爱,今天小智给大家送上一套科研绘图的必备书籍——《R语言数据可视化:科技图表绘制》。下面我会详细给大家介绍这套图书,文末留有领取方式。 图书介绍 《R语言数据可视化:科技图表绘制》结合编者多年的数据分析…

jenkins安装和使用 (二)

参考视频资料 https://www.bilibili.com/video/BV1bS4y1471A?p10&vd_sourcee0dcd147bd5d730317de804d788cd6f9 安装maven插件 新建item 配置构建信息 项目地址替换为自己的实际地址 其余保持先保持默认 先然后在主页就看到了这个项目 查看控制台输出 稍等一…

Java 17的新特性

Java 17引入了多项新特性,以下是一些重要的更新: 增强的伪随机数生成器(JEP 356) Java 17为伪随机数生成器(PRNG)提供了新的接口类型和实现,包括可跳转的PRNG和另一类可拆分的PRNG算法&#xf…

聊聊 Mybatis 动态 SQL

这篇文章,我们聊聊 Mybatis 动态 SQL ,以及我对于编程技巧的几点思考 ,希望对大家有所启发。 1 什么是 Mybatis 动态SQL 如果你使用过 JDBC 或其它类似的框架,你应该能理解根据不同条件拼接 SQL 语句有多痛苦,例如拼…

得物面试:什么是零复制?说说 零复制 底层原理?(吊打面试官)

尼恩说在前面 在40岁老架构师 尼恩的读者交流群(50)中,最近有小伙伴拿到了一线互联网企业如得物、阿里、滴滴、极兔、有赞、希音、百度、网易、美团的面试资格,遇到很多很重要的零复制的问题: 说一说Rocketmq、是如何实现每秒上百万数据的超…

ACIS采用自底向上的方式创建box

在学习任何CAD几何内核的时候,首先需要学习的是这个几何内核的几何拓扑数据结构,学习完毕后,一个很好的练习是自己采用底层的几何、拓扑构建API,创建一个box。通过这个练习,可以加深初学者对所学几何内核数据结构的理解…

AXI三板斧之Outstanding、Out-of-order、interleaving

1、AXI三板斧之Outstanding 可以不用等单个命令的响应,直接连续发送N个命令(N>1),假设Slave端的Outstanding能力为N时(N>1),那么Master端可以在Slave不返回读数据的情况下,连…

从小白到大神:算法工程师的核心竞争力养成计划!

从小白到大神:算法工程师的核心竞争力养成计划! 我们会发现,本科及研究生所学的知识,比如高等数学、线性代数、传统机器学习方法及深度学习理论等这些,都只是作为算法学习的基础,并不能成为算法工程师的核…

React@16.x(34)动画(中)

目录 3,SwitchTransition3.1,原理3.1.2,key3.1.2,mode 3.2,举例3.3,结合 animate.css 4,TransitionGroup4.1,其他属性4.1.2,appear4.1.2,component4.1.3&…

【Python】对应接口url 被编码后的处理

Python 系列 文章目录 Python 系列前言一、网页链接是什么?二、使用步骤1.解码 总结 前言 提示:这里可以添加本文要记录的大概内容: 我们在查找网页的开发代码的时候,可能经常查看到接口的链接是: %7B%22funName%22%…

昇思25天学习打卡营第1天|快速入门

一、简介: 本节通过MindSpore已经封装好的API,快速实现一个深度学习模型的数据集准备、训练评估,模型参数保存和加载,对新手朋友十分友好。这里非常感谢华为昇思团队在算力和代码方面的指导。 二、环境准备: 在开始…

YOLOv10改进 | Conv篇 |YOLOv10引入DBB卷积(助力涨点)

1. DBB介绍 1.1 摘要:我们提出了卷积神经网络(ConvNet)的通用构建块来提高性能,而无需任何推理时间成本。 该块被称为多样化分支块(DBB),它通过组合不同尺度和复杂度的不同分支来丰富特征空间,从而增强单个卷积的表示能力,包括卷积序列、多尺度卷积和平均池化。 训练后…

参花实业增收不增利:曾被执行近3000万,资产负债比率高达166%

《港湾商业观察》廖紫雯 日前,参花实业控股有限公司(以下简称:参花实业)再度递表港交所,保荐机构为中泰国际,参花实业国内运营主体为固安县参花面粉有限公司、廊坊固安县参花企业管理有限公司。 作为一家…

【linux】Valgrind工具集详解(十六):交叉编译、移植到arm(失败)

1、源码下载 官网:https://valgrind.org/ 源码:https://valgrind.org/downloads/current.html 2、配置 ./configure CC=arm-linux-gnueabihf-gcc \CXX=arm-linux-gnueabihf-g++ \AR=arm-linux-gnueabihf-ar \--host=arm-linux-gnueabihf \--pr

Mkdocs中文系列教程补充(1)

什么是requirements.txt 我的理解是mkdocs依赖的py库 第一次建立MKdocs文档使用 mkdocs new . 完后,比较建议执行一下: pip install -r requirements.txt 不然mkdocs serve后会出现什么 xxx not found ,比如下面这位老哥 示例 mkdocs …

CleanShot X for Mac v4.7 屏幕截图录像工具(保姆级教程,小白轻松上手,简单易学)

Mac分享吧 文章目录 一、准备工作二、部分特有功能效果1、截图软件的普遍常用功能(画框、箭头、加文字等)都具备,不再详细介绍2、ABCD、1234等信息标注(每按一下鼠标,即各是A、B、C、D...等)3、截图更换背…

MySQL日志(三):数据安全

先来看一个结论:只要redo log和binlog保证持久化到磁盘, 就能确保MySQL异常重启后, 数据可以恢复。 binlog写入逻辑 binlog的写入逻辑比较简单: 事务执行过程中, 先把日志写到binlog cache, 事务提交的时候…

【论文阅读】-- MIRIA:用于时空交互数据原位可视化和分析的混合现实工具包

MIRIA: A Mixed Reality Toolkit for the In-Situ Visualization and Analysis of Spatio-Temporal Interaction Data 摘要1 引言2 背景及相关工作2.1 用户交互和运动数据分析2.2 沉浸式分析 3 时空用户交互分析3.1 现有系统的用例和分析3.2 要求 4 MIRIA 工具包4.1 一般概念4.…

基于matlab的RRT算法路径规划(附带案例源码)

文章中的所有案例均为博主手动复现,用于记录博主学习路径规划的过程,如有不妥,欢迎在评论区交流 目录 1 标准RRT1.1 算法原理1.2 演示 2 GBRRT2.1 算法原理2.2 算法演示 3 RRT-STAR3.1 算法原理3.2 算法演示 4 RRT-CONNECT4.1 算法原理4.2 算…

[NewStarCTF 2023 公开赛道]R!C!E!

好久没写了。今天儿弄一个rce 题很直接&#xff0c;好久没这么直白的题了&#xff0c;看源码 <?php highlight_file(__FILE__); if(isset($_POST[password])&&isset($_POST[e_v.a.l])){$passwordmd5($_POST[password]);$code$_POST[e_v.a.l];if(substr($password…