强化学习 — 马尔科夫决策过程(MDP)
一、马尔科夫过程(Markov Process)
马尔科夫性 某一状态信息包含了所有相关的历史,只要当前状态可知,所有的历史信息都不再需要,当前状态就可以决定未来,则认为该状态具有马尔科夫性
P ( S t + 1 ∣ S t ) = p ( S t + 1 ∣ S 1 , S 2 , ⋯ , S t )
P(S_{t+1}|S_t) = p(S_{t+1}|S_1, S_2, \cdots , S_t)
P ( S t + 1 ∣ S t ) = p ( S t + 1 ∣ S 1 , S 2 , ⋯ , S t )
马尔科夫过程 又叫做马尔科夫链(Markov Chain)
,它是一个无记忆的随机过程,可以用一个元组<S, P>
表示,其中
S
是有限数量的状态集 S = s 1 , s 2 , s 3 , ⋯ , s t S ={s_1, s_2, s_3, \cdots, s_t} S = s 1 , s 2 , s 3 , ⋯ , s t
P
是状态转移概率矩阵 p ( S t + 1 = s ′ ∣ s t = s ) , 其 中 s ′ 表 示 下 一 时 刻 的 状 态 , s 表 示 当 前 状 态 p(S_{t+1} = s'|s_t=s) , 其中 s' 表示下一时刻的状态,s 表示当前状态 p ( S t + 1 = s ′ ∣ s t = s ) , 其 中 s ′ 表 示 下 一 时 刻 的 状 态 , s 表 示 当 前 状 态
二、马尔科夫奖励过程(Markov Reward Process)
马尔科夫奖励过程 是在马尔科夫过程基础上增加了奖励函数 R R R 和衰减系数 γ \gamma γ , 用 < S , R , P , γ > <S, R,P, \gamma> < S , R , P , γ > 表示
R R R : 表示 S S S 状态下某一时刻的状态S t S_t S t 在下一个时刻 ( t + 1 ) (t + 1) ( t + 1 ) 能获得的奖励的期望
R s = E [ R t + 1 ∣ S t = s ]
R_s = E[R_{t+1}|S_t=s]
R s = E [ R t + 1 ∣ S t = s ]
G t G_t G t : 收获 G t G_t G t 为在一个马尔科夫奖励链上从t时刻开始往后所有的奖励的有衰减的收益总和
G t = R t + 1 + γ R t + 2 + γ 2 R t + 3 + ⋯ + γ T − t − 1 R T
G_t = R_{t+1} + \gamma R_{t+2} + \gamma^2 R_{t+3} + \cdots + \gamma^{T-t-1}R_T
G t = R t + 1 + γ R t + 2 + γ 2 R t + 3 + ⋯ + γ T − t − 1 R T
γ \gamma γ : 折扣因子( D i s c o u n t f a c t o r γ ∈ [ 0 , 1 ] ) (Discount \; factor γ ∈ [0, 1]) ( D i s c o u n t f a c t o r γ ∈ [ 0 , 1 ] )
1、为了避免出现状态循环的情况
2、系统对于将来的预测并不一定都是准确的,所以要打折扣
很显然 越靠近1,考虑的利益越长远。
V ( s ) V(s) V ( s ) : 状态价值函数(state value function) 表示从从该状态开始的马尔科夫链收获G t G_t G t 的期望
v ( s ) = E [ G t ∣ S t = s ]
v(s) = E[G_t|S_t = s]
v ( s ) = E [ G t ∣ S t = s ]
例子:
奖励:S 1 S_1 S 1 +5, S 7 S_7 S 7 +10, 其余奖励为 0,则 R = [ 5 , 0 , 0 , 0 , 0 , 0 , 10 ] R = [5,0,0,0,0,0,10] R = [ 5 , 0 , 0 , 0 , 0 , 0 , 1 0 ]
γ = 0.5 \gamma = 0.5 γ = 0 . 5
$S_4, S_5, S_6, S_7: $ 0 + 0.5*0 + 0.5*0 + 0.125 *10 = 1.25
S 4 , S 3 , S 2 , S 1 S_4,S_3,S_2,S_1 S 4 , S 3 , S 2 , S 1 : 0 + 0.5 ×0 + 0.25×0 + 0.125×5 = 0.625
Bellman Equation 贝尔曼方程
v ( s ) = E [ G t ∣ S t = s ] v(s) = E[G_t|S_t = s] v ( s ) = E [ G t ∣ S t = s ]
= E [ R t + 1 + γ R t + 2 + γ 2 R t + 3 + ⋯ ∣ S t = s ] = E[R_{t+1} + \gamma R_{t+2} + \gamma^2 R_{t+3} + \cdots|S_t = s] = E [ R t + 1 + γ R t + 2 + γ 2 R t + 3 + ⋯ ∣ S t = s ]
= E [ R t + 1 + γ ( R t + 2 + R t + 3 + ⋯ ) ∣ S t = s ] = E[R_{t+1} + \gamma(R_{t+2} + R_{t+3} + \cdots)|S_t=s] = E [ R t + 1 + γ ( R t + 2 + R t + 3 + ⋯ ) ∣ S t = s ]
= E [ R t + 1 + γ v ( S t + 1 ) ∣ S t = s ] =E[R_{t+1} + \gamma v(S_{t+1})|S_t = s] = E [ R t + 1 + γ v ( S t + 1 ) ∣ S t = s ]
= E [ R t + 1 ∣ S t = s ] ⏟ 当 前 的 奖 励 + γ E [ v ( S t + 1 ) ∣ S t = s ] ⏟ 下 一 时 刻 状 态 的 价 值 期 望 =\underbrace{E[R_{t+1}|S_t=s]}_{当前的奖励} + \underbrace{\gamma E[v(S_{t+1})|S_t = s]}_{下一时刻状态的价值期望} = 当 前 的 奖 励 E [ R t + 1 ∣ S t = s ] + 下 一 时 刻 状 态 的 价 值 期 望 γ E [ v ( S t + 1 ) ∣ S t = s ]
使用贝尔曼方程状态价值V V V 可以表示为:
V ( s ) = R ( s ) ⏟ I m m e d i a t e r e w a r d + γ ∑ s ′ ∈ S P ( s ′ ∣ s ) V ( s ′ ) ⏟ D i s c o u n t e d s u m o f f u t u r e r e w a r d
V(s) = \underbrace{R(s)}_{Immediate \; reward} + \underbrace{\gamma \sum_{s' \in S}P(s'|s)V(s')}_{Discounted \; sum \; of \; future \; reward}
V ( s ) = I m m e d i a t e r e w a r d R ( s ) + D i s c o u n t e d s u m o f f u t u r e r e w a r d γ s ′ ∈ S ∑ P ( s ′ ∣ s ) V ( s ′ )
S 表示下一时刻的所有状态
s’ 表示下一时刻可能的状态
三、马尔科夫决策过程(Markov Decision Process)
马尔科夫决策过程 是在马尔科夫奖励过程的基础上加了 Decision
过程,相当于多了一个动作集合,可以用 < S , A , P , R , γ > <S, A, P, R, \gamma> < S , A , P , R , γ >
A A A 表示有限的行为集合
S S S 表示有限的状态集合
P a P^a P a is dynamics / transition model for each action
P ( s t + 1 = s ′ ∣ s t = s , a t = a ) P(s_{t+1} = s'|s_t = s, a_t = a) P ( s t + 1 = s ′ ∣ s t = s , a t = a )
R R R 是奖励函数 R ( s t = s , a t = a ) = E [ R t ∣ s t = s , a t = a ] R(s_t=s, a_t = a) = E[R_t|s_t=s, a_t=a] R ( s t = s , a t = a ) = E [ R t ∣ s t = s , a t = a ]
策略 (Policy)
用 π \pi π 表示策略的集合,其元素 π ( a ∣ s ) \pi(a|s) π ( a ∣ s ) 表示某一状态 s
采取可能的行为 a
的概率
π ( a ∣ s ) = P ( A t = a ∣ S t = s )
\pi(a|s) = P(A_t=a|S_t=s)
π ( a ∣ s ) = P ( A t = a ∣ S t = s )
在马尔科夫奖励过程中 策略 π \pi π 满足以下方程,可以参照下面图来理解
状 态 转 移 概 率 : P π ( s ′ ∣ s ) = ∑ a ∈ A π ( a ∣ s ) P ( s ′ ∣ s , a ) 奖 励 函 数 : R π ( s ) = ∑ a ∈ A π ( a ∣ s ) R ( s , a )
状态转移概率:\quad P^\pi(s'|s) = \sum_{a \in A}\pi(a|s)P(s'|s, a) \\
奖励函数:\quad R^\pi(s) = \sum_{a\in A}\pi(a|s)R(s,a)
状 态 转 移 概 率 : P π ( s ′ ∣ s ) = a ∈ A ∑ π ( a ∣ s ) P ( s ′ ∣ s , a ) 奖 励 函 数 : R π ( s ) = a ∈ A ∑ π ( a ∣ s ) R ( s , a )
状态转移概率可以描述为:在执行策略 π \pi π 时,状态从 s 转移至 s’ 的概率等于执行该状态下所有行为的概率与对应行为能使状态从 s 转移至 s’ 的概率的乘积的和。参考下图
奖励函数可以描述为:在执行策略 π \pi π 时获得的奖励等于执行该状态下所有行为的概率与对应行为产生的即时奖励的乘积的和。
我们引入策略,也可以理解为行动指南,更加规范的描述个体的行为,既然有了行动指南,我们要判断行动指南的价值,我们需要再引入基于策略的价值函数。
基于策略的状态价值函数(state value function)
V ( s ) V(s) V ( s ) 表示从状态 s s s 开始,遵循当前策略时所获得的收获的期望
v π ( s ) = E π [ G t ∣ S t = s ]
v_{\pi}(s) = E_{\pi}[G_t|S_t = s]
v π ( s ) = E π [ G t ∣ S t = s ]
其中 G t G_t G t 可以参照马科夫奖励过程。我们有了价值衡量标准,如果状态 s 是一个好的状态,如何选择动作到达这个状态,这时就需要判断动作的好坏,衡量行为价值。
基于策略的行为价值函数(action value function)
q π ( s , a ) q_{\pi}(s,a) q π ( s , a ) 当前状态s执行某一具体行为a所能的到的收获的期望
q π ( s , a ) = E π [ G t ∣ S t = s , A t = a ]
q_{\pi}(s, a) = E_{\pi}[G_t|S_t=s, A_t=a]
q π ( s , a ) = E π [ G t ∣ S t = s , A t = a ]
根据 Bellman 公式推导可得(参照马尔科夫奖励过程中 V 的推导)
q π ( s , a ) = R ( s , a ) + γ ∑ s ′ ∈ S P ( s ′ ∣ s , a ) ⋅ V π ( s ′ )
q_{\pi}(s,a) = R(s, a) + \gamma \sum_{s' \in S} P_{(s'|s, a)} \cdot V_{\pi}(s')
q π ( s , a ) = R ( s , a ) + γ s ′ ∈ S ∑ P ( s ′ ∣ s , a ) ⋅ V π ( s ′ )
在某一个状态下采取某一个行为的价值,可以分为两部分:其一是离开这个状态的价值,其二是所有进入新的状态的价值于其转移概率乘积的和。参考下图右理解
由状态价值函数和行为价值函数的定义,可得两者的关系
v π ( s ) = ∑ a ∈ A π ( a ∣ s ) ⋅ q π ( s , a )
v_{\pi}(s) = \sum_{a \in A}\pi(a|s) \cdot q_{\pi}(s,a)
v π ( s ) = a ∈ A ∑ π ( a ∣ s ) ⋅ q π ( s , a )
我们知道策略就是用来描述各个不同状态下执行各个不同行为的概率,而状态价值是遵循当前策略时所获得的收获的期望,即状态 s 的价值体现为在该状态下遵循某一策略而采取所有可能行为的价值按行为发生概率的乘积求和。参照下图左理解
上面两个公式组合可得 Bellman Expectation Equation
v π ( s ) = ∑ a ∈ A π ( a ∣ s ) ( R ( s , a ) + γ ∑ s ′ ∈ S P ( s ′ ∣ s , a ) ⋅ V π ( s ′ ) )
v_{\pi}(s) = \sum_{a \in A}\pi(a|s)\left(R(s, a) + \gamma \sum_{s' \in S} P_{(s'|s, a)} \cdot V_{\pi}(s')\right)
v π ( s ) = a ∈ A ∑ π ( a ∣ s ) ( R ( s , a ) + γ s ′ ∈ S ∑ P ( s ′ ∣ s , a ) ⋅ V π ( s ′ ) )
q π ( s , a ) = R ( s , a ) + γ ∑ s ′ ∈ S P ( s ′ ∣ s , a ) ⋅ ∑ a ′ ∈ A π ( a ′ ∣ s ′ ) ⋅ q π ( s ′ , a ′ )
q_{\pi}(s,a) = R(s, a) + \gamma \sum_{s' \in S} P_{(s'|s, a)} \cdot \sum_{a' \in A}\pi(a'|s') \cdot q_{\pi}(s',a')
q π ( s , a ) = R ( s , a ) + γ s ′ ∈ S ∑ P ( s ′ ∣ s , a ) ⋅ a ′ ∈ A ∑ π ( a ′ ∣ s ′ ) ⋅ q π ( s ′ , a ′ )
实例
假设γ = 1 \gamma = 1 γ = 1 上图数字值为对于状态价值v(s),根据公式价值函数需要下一个状态的价值,终点的状态价值为0,利用递归的方式即可得到
五、Decision Making in MDP
可以把 Decision Making 分为两个过程
Prediction (评估给定的策略)
输入: MDP< S , A , P , R , γ > <S, A, P, R, \gamma> < S , A , P , R , γ > 和 策略 $\pi $
输出:value function V π V_\pi V π
Control (寻找最优策略)
输入:MDP$<S, A, P, R, \gamma> $
输出:最优的价值函数 V ∗ V^* V ∗ 和最优的策略 π ∗ \pi^* π ∗
1、Policy evaluation on MDP
Prediction 满足动态规划的条件:
贝尔曼等式可以递归分解
value function 可以被存储和重用
Iteration on Bellman exception backup
我们可以用所有的状态 s 在 t 时刻的价值 v t ( s ′ ) v_t(s') v t ( s ′ ) 来更新 v t + 1 ( s ) v_{t+1}(s) v t + 1 ( s ) 时刻的价值,然后迭代下去
V t + 1 ( s ) = ∑ a ∈ A π ( a ∣ s ) ( R ( s , a ) + γ ∑ s ′ ∈ S P ( s ′ ∣ s , a ) ⋅ V t ( s ′ ) ) V t + 1 ( s ) = R π ( s ) + γ ∑ s ′ ∈ S P π ( s ′ ∣ s ) ⋅ V t ( s ′ )
V_{t+1}(s) = \sum_{a \in A}\pi(a|s)(R(s, a) + \gamma \sum_{s' \in S} P_{(s'|s, a)} \cdot V_t(s')) \\
V_{t+1}(s) = R_\pi(s) + \gamma \sum_{s' \in S} P_{\pi}(s'|s) \cdot V_t(s')
V t + 1 ( s ) = a ∈ A ∑ π ( a ∣ s ) ( R ( s , a ) + γ s ′ ∈ S ∑ P ( s ′ ∣ s , a ) ⋅ V t ( s ′ ) ) V t + 1 ( s ) = R π ( s ) + γ s ′ ∈ S ∑ P π ( s ′ ∣ s ) ⋅ V t ( s ′ )
最优状态价值函数指的是在从所有策略产生的状态价值函数中,选取使状态s价值最大的函数:
V ∗ ( s ) = m a x π V π ( s )
V_*(s) = max_\pi V_\pi(s)
V ∗ ( s ) = m a x π V π ( s )
π ∗ ( s ) = a r g m a x π V ∗ ( s )
\pi^*(s) = argmax_\pi V^*(s)
π ∗ ( s ) = a r g m a x π V ∗ ( s )
对于任何MDP,下面几点成立:
1.存在一个最优策略,比任何其他策略更好或至少相等;
2.所有的最优策略有相同的最优价值函数;
3.所有的最优策略具有相同的行为价值函数。
根据以上几点,我们可以最大化最优行为价值函数 找到最优策略
π ∗ ( a ∣ s ) = { 1 , if a = a r g m a x q ∗ ( s , a ) 0 , otherwise
\pi^*(a|s) =
\begin{cases}
1, & \text{if $a = argmax \; q^*(s,a)$} \\
0, & \text{otherwise}
\end{cases}
π ∗ ( a ∣ s ) = { 1 , 0 , if a = a r g ma x q ∗ ( s , a ) otherwise
策略的个数为 ∣ A ∣ ∣ S ∣ |A|^{|S|} ∣ A ∣ ∣ S ∣ 个,对于状态比较多的环境,计算量太大
更加高效的方法是 policy iteration 和 value iteration
2、MDP Control
(1) Policy Iteration
分为两个部分
Evaluate the policy $ \pi $,(通过给定的 $ \pi $ 计算 V)
Improve the policy $ \pi $,(通过贪心策略)
π ′ = g r e e d y ( v π )
\pi' = greedy(v^{\pi})
π ′ = g r e e d y ( v π )
如果我们有 一个 策略 π \pi π ,我们可以用策略 估计出它的状态价值函数 v π ( s ) v_\pi(s) v π ( s ) , 然后根据策略改进提取出更好的策略 π ′ \pi' π ′ ,接着再计算 v π ′ ( s ) v_{\pi'}(s) v π ′ ( s ) , 然后再获得更好的 策略 π ′ ′ \pi'' π ′ ′ ,直到相关满足相关终止条件。
v i ( s ) = ∑ a ∈ A π ( a ∣ s ) ( R ( s , a ) + γ ∑ s ′ ∈ S P ( s ′ ∣ s , a ) ⋅ v i − 1 ( s ′ ) )
v_i(s) = \sum_{a \in A}\pi(a|s)(R(s, a) + \gamma \sum_{s' \in S} P_{(s'|s, a)} \cdot v_{i-1}(s'))
v i ( s ) = a ∈ A ∑ π ( a ∣ s ) ( R ( s , a ) + γ s ′ ∈ S ∑ P ( s ′ ∣ s , a ) ⋅ v i − 1 ( s ′ ) )
q π i ( s , a ) = R ( s , a ) + γ ∑ s ′ ∈ S P ( s ′ ∣ s , a ) ⋅ v π i ( s ′ ) π i + 1 ( s ) = a r g m a x a q π i ( s , a )
q^{\pi_i}(s,a) = R(s, a) + \gamma \sum_{s' \in S} P_{(s'|s,a)} \cdot v^{\pi_i}(s') \\
\pi_{i+1}(s) = argmax_a \; q^{\pi_i}(s,a)
q π i ( s , a ) = R ( s , a ) + γ s ′ ∈ S ∑ P ( s ′ ∣ s , a ) ⋅ v π i ( s ′ ) π i + 1 ( s ) = a r g m a x a q π i ( s , a )
可以把 q π i ( s , a ) q^{\pi_i}(s,a) q π i ( s , a ) 想象成一个表格,其中横轴代表不同的状态,纵轴代表每种状态下不同的动作产生的价值,然后选择一个最大的行为价值作为当前的状态价值
Bellman Optimality Equation
v ∗ ( s ) = m a x a q ∗ ( s , a )
v_*(s) = max_a q_*(s,a)
v ∗ ( s ) = m a x a q ∗ ( s , a )
q ∗ ( s , a ) = R ( s , a ) + γ ∑ s ′ ∈ S P s ′ s a ⋅ V ∗ ( s ′ )
q_*(s,a) = R(s, a) + \gamma \sum_{s' \in S} P_{s's}^a \cdot V_*(s')
q ∗ ( s , a ) = R ( s , a ) + γ s ′ ∈ S ∑ P s ′ s a ⋅ V ∗ ( s ′ )
把上面两个式子结合起来有Bellman Optimality Equation
v ∗ ( s ) = m a x a ( R ( s , a ) + γ ∑ s ′ ∈ S P s ′ s a ⋅ V ∗ ( s ′ ) )
v_*(s) = max_a (R(s, a) + \gamma \sum_{s' \in S} P_{s's}^a \cdot V_*(s'))
v ∗ ( s ) = m a x a ( R ( s , a ) + γ s ′ ∈ S ∑ P s ′ s a ⋅ V ∗ ( s ′ ) )
q ∗ ( s , a ) = R ( s , a ) + γ ∑ s ′ ∈ S P s ′ s a ⋅ m a x a ′ q ∗ ( s ′ , a ′ )
q_*(s,a) = R(s, a) + \gamma \sum_{s' \in S} P_{s's}^a \cdot max_{a'}q_*(s', a')
q ∗ ( s , a ) = R ( s , a ) + γ s ′ ∈ S ∑ P s ′ s a ⋅ m a x a ′ q ∗ ( s ′ , a ′ )
(2) Value Iteration
如果我们解决了子问题的 V ∗ ( s ′ ) V^*(s') V ∗ ( s ′ )
那么我们就可以通过反复迭代Bellman Optimality Equation 找到 $ V^*(s) $
v ∗ ( s ) = m a x a q ∗ ( s , a ) ⇓ v i + 1 ( s ) ← m a x a ∈ A ( R ( s , a ) + γ ∑ s ′ ∈ S P ( s ′ ∣ s , a ) ⋅ V i ( s ′ ) )
v_*(s) = max_a q_*(s,a) \\
\Downarrow \\
v_{i+1}(s) \leftarrow max_{a \in A} \; \left(R(s, a) + \gamma \sum_{s' \in S} P_{(s'|s,a)} \cdot V_i(s')\right)
v ∗ ( s ) = m a x a q ∗ ( s , a ) ⇓ v i + 1 ( s ) ← m a x a ∈ A ( R ( s , a ) + γ s ′ ∈ S ∑ P ( s ′ ∣ s , a ) ⋅ V i ( s ′ ) )
π ∗ ( s ) = a r g m a x a q π ( s , a ) ⇓ π ∗ ( s ) ← a r g m a x a ( R ( s , a ) + γ ∑ s ′ ∈ S P ( s ′ ∣ s , a ) ⋅ V e n d ( s ′ ) )
\pi^*(s) = argmax_a \; q^{\pi}(s,a) \\
\Downarrow \\
\pi^*(s) \leftarrow argmax_a \; \left(R(s, a) + \gamma \sum_{s' \in S} P_{(s'|s,a)} \cdot V_{end}(s')\right)
π ∗ ( s ) = a r g m a x a q π ( s , a ) ⇓ π ∗ ( s ) ← a r g m a x a ( R ( s , a ) + γ s ′ ∈ S ∑ P ( s ′ ∣ s , a ) ⋅ V e n d ( s ′ ) )
六 后记
关于直观的理解 Policy Iteration 和 Value Iteration 可以访问下面的网址
https://cs.stanford.edu/people/karpathy/reinforcejs/gridworld_dp.html