 Markov decision process

Markov decision processes (MDPs), named after Andrey Markov, provide a mathematical framework for modeling decisionmaking in situations where outcomes are partly random and partly under the control of a decision maker. MDPs are useful for studying a wide range of optimization problems solved via dynamic programming and reinforcement learning. MDPs were known at least as early as the 1950s (cf. Bellman 1957). Much research in the area was spawned due to Ronald A. Howard's book, Dynamic Programming and Markov Processes, in 1960. Today they are used in a variety of areas, including robotics, automated control, economics and manufacturing.
More precisely, a Markov Decision Process is a discrete time stochastic control process. At each time step, the process is in some state s, and the decision maker may choose any action a that is available in state s. The process responds at the next time step by randomly moving into a new state s', and giving the decision maker a corresponding reward R_{a}(s,s').
The probability that the process moves into its new state s' is influenced by the chosen action. Specifically, it is given by the state transition function P_{a}(s,s'). Thus, the next state s' depends on the current state s and the decision maker's action a. But given s and a, it is conditionally independent of all previous states and actions; in other words, the state transitions of an MDP possess the Markov property.
Markov decision processes are an extension of Markov chains; the difference is the addition of actions (allowing choice) and rewards (giving motivation). Conversely, if only one action exists for each state and all rewards are zero, a Markov decision process reduces to a Markov chain.
Contents
Definition
A Markov decision process is a 4tuple , where
 S is a finite set of states,
 A is a finite set of actions (alternatively, A_{s} is the finite set of actions available from state s),
 is the probability that action a in state s at time t will lead to state s' at time t + 1,
 R_{a}(s,s') is the immediate reward (or expected immediate reward) received after transition to state s' from state s with transition probability P_{a}(s,s').
(The theory of Markov decision processes does not actually require S or A to be finite,^{[citation needed]} but the basic algorithms below assume that they are finite.)
Problem
The core problem of MDPs is to find a policy for the decision maker: a function π that specifies the action π(s) that the decision maker will choose when in state s. Note that once a Markov decision process is combined with a policy in this way, this fixes the action for each state and the resulting combination behaves like a Markov chain.
The goal is to choose a policy π that will maximize some cumulative function of the random rewards, typically the expected discounted sum over a potentially infinite horizon:
 (where we choose a_{t} = π(s_{t}))
where is the discount factor and satisfies . (For example, γ = 1 / (1 + r) when the discount rate is r.) γ is typically close to 1.
Because of the Markov property, the optimal policy for this particular problem can indeed be written as a function of s only, as assumed above.
Algorithms
MDPs can be solved by linear programming or dynamic programming. In what follows we present the latter approach.
Suppose we know the state transition function P and the reward function R, and we wish to calculate the policy that maximizes the expected discounted reward.
The standard family of algorithms to calculate this optimal policy requires storage for two arrays indexed by state: value V, which contains real values, and policy π which contains actions. At the end of the algorithm, π will contain the solution and V(s) will contain the discounted sum of the rewards to be earned (on average) by following that solution from state s.
The algorithm has the following two kinds of steps, which are repeated in some order for all the states until no further changes take place. They are
Their order depends on the variant of the algorithm; one can also do them for all states at once or state by state, and more often to some states than others. As long as no state is permanently excluded from either of the steps, the algorithm will eventually arrive at the correct solution.
Notable variants
Value iteration
In value iteration (Bellman 1957), which is also called backward induction, the π array is not used; instead, the value of π(s) is calculated whenever it is needed. Shapley's 1953 paper on stochastic games included as a special case the value iteration method for MDPs, but this was recognized only later on.^{[1]}
Substituting the calculation of π(s) into the calculation of V(s) gives the combined step:
This update rule is iterated for all states s until it converges with the lefthand side equal to the righthand side (which is the Bellman equation for this problem).
Policy iteration
In policy iteration (Howard 1960), step one is performed once, and then step two is repeated until it converges. Then step one is again performed once and so on.
Instead of repeating step two to convergence, it may be formulated and solved as a set of linear equations.
This variant has the advantage that there is a definite stopping condition: when the array π does not change in the course of applying step 1 to all states, the algorithm is completed.
Modified policy iteration
In modified policy iteration (van Nunen, 1976; Puterman and Shin 1978), step one is performed once, and then step two is repeated several times. Then step one is again performed once and so on.
Prioritized sweeping
In this variant, the steps are preferentially applied to states which are in some way important  whether based on the algorithm (there were large changes in V or π around those states recently) or based on use (those states are near the starting state, or otherwise of interest to the person or program using the algorithm).
Extensions and generalizations
A Markov decision process is a stochastic game with only one player.
Partial observability
Main article: partially observable Markov decision processThe solution above assumes that the state s is known when action is to be taken; otherwise π(s) cannot be calculated. When this assumption is not true, the problem is called a partially observable Markov decision process or POMDP.
Reinforcement Learning
If the probabilities or rewards are unknown, the problem is one of reinforcement learning (Sutton and Barto, 1998).
For this purpose it is useful to define a further function, which corresponds to taking the action a and then continuing optimally (or according to whatever policy one currently has):
While this function is also unknown, experience during learning is based on (s,a) pairs (together with the outcome s'); that is, "I was in state s and I tried doing a and s' happened"). Thus, one has an array Q and uses experience to update it directly. This is known as Qlearning.
The power of reinforcement learning lies in its ability to solve the Markov decision process without computing the transition probabilities; note that transition probabilities are needed in value and policy iteration. Also, reinforcement learning can be combined with function approximation, and thereby one can solve problems with a very large number of states. Reinforcement Learning can also be handily performed within Monte Carlo simulators of systems.
Continuoustime Markov Decision Process
In Discretetime Markov Decision Process, decisions are made at discrete time epoch. However, for Continuoustime Markov Decision Process, decisions can be made at any time when decision maker wants. Different than Discretetime Markov Decision Process, Continuoustime Markov Decision Process could better model the decision making process when the interested system has continuous dynamics, i.e., the system dynamics is defined by Partial Differential Equations(PDEs).
Definition
In order to discuss the Continuoustime Markov Decision Process, we introduce two sets of notations:
If the state space and action space are finite,
 : State space;
 : Action space;
 q(i  j,a): , transition rate function;
 R(i,a): , a reward function.
If the state space and action space are continuous,
 : State space.;
 : Space of possible control;
 f(x,u): , a transition rate function;
 r(x,u): , a reward rate function,((dR(x(t),u(t)) = r(x(t),u(t))dt and R(x,u) is the reward function we discussed in previous case.)
Problem
Like the Discretetime Markov Decision Processes, in Continuoustime Markov Decision Process we want to find the optimal policy or control which could give us the optimal expected integrated reward:
Where
Linear programming formulation
If the state space and action space are finite, we could use linear programming formulation to find the optimal policy, which was one of the earliest solution approaches. Here we only consider the ergodic model, which means our continuoustime MDP become ergodic Continuoustime Markov Chain under stationary policy. Under this assumption, although the decision maker could make decision at any time, on the current state, he could not get more benefit to make more than one actions. It is better for him to take action only at the time when system transit from current state to another state. Under some conditions,(for detail check Corollary 3.14 of ContinuousTime Markov Decision Processes), if our optimal value function V ^{*} is independent of state i, we will have a following equation:
If there exists a function h, then will be the smallest g could satisfied the above equation. In order to find the , we could have the following linear programming model:
 Primal linear program(PLP)
 Dual linear program(DLP)
y(i,a) is a feasible solution to the DLP if y(i,a) is nonnative and satisfied the constraints in the DLP problem. A feasible solution y * (i,a) to the DLP is said to be an optimal solution if
for all feasible solution y(i,a) to the DLP. Once we found the optimal solution y * (i,a), we could use those optimal solution to establish the optimal policies.
HamiltonJacobiBellman equation
In Continuoustime MDP, if the state space and action space are continuous, the optimal criterion could be found by solving HamiltonJacobiBellman partial differential equation. In order to discuss the HJB equation, we need to reformulate our problem
D() is the terminal reward function, x(t) is the system state vector, u(t) is the system control vector we try to find. f() shows how the state vector change over time. HamiltonJacobiBellman equation is as follows:
We could solve the equation to find the optimal control u(t), which could give us the optimal value V ^{*}
Application
Queueing system, epidemic processes, Population process.
Alternative notations
The terminology and notation for MDPs are not entirely settled. There are two main streams — one focuses on maximization problems from contexts like economics, using the terms action, reward, value and calling the discount factor β or γ, while the other focuses on minimization problems from engineering and navigation, using the terms control, cost, costtogo and calling the discount factor α. In addition, the notation for the transition probability varies.
in this article alternative comment action a control u reward R cost g g is the negative of R value V costtogo J J is the negative of V policy π policy μ discounting factor discounting factor α transition probability P_{a}(s,s') transition probability p_{ss'}(a) In addition, transition probability is sometimes written Pr(s,a,s'), Pr(s'  s,a) or, rarely, p_{s's}(a).
See also
 Partially observable Markov decision process
 Dynamic programming
 Bellman equation for applications to economics.
 Hamilton–Jacobi–Bellman equation
 Optimal control
Notes
 ^ Lodewijk Kallenberg, Finite state and action MDPs, in Eugene A. Feinberg, Adam Shwartz (eds.) Handbook of Markov decision processes: methods and applications, Springer, 2002, ISBN 0792374592
References
 R. Bellman. A Markovian Decision Process. Journal of Mathematics and Mechanics 6, 1957.
 R. E. Bellman. Dynamic Programming. Princeton University Press, Princeton, NJ, 1957. Dover paperback edition (2003), ISBN 0486428095.
 Ronald A. Howard Dynamic Programming and Markov Processes, The M.I.T. Press, 1960.
 D. Bertsekas. Dynamic Programming and Optimal Control. Volume 2, Athena, MA, 1995.
 M. L. Puterman. Markov Decision Processes. Wiley, 1994.
 H.C. Tijms. A First Course in Stochastic Models. Wiley, 2003.
 Sutton, R. S. and Barto A. G. Reinforcement Learning: An Introduction. The MIT Press, Cambridge, MA, 1998.
 J.A. E. E van Nunen. A set of successive approximation methods for discounted Markovian decision problems. Z. Operations Research, 20:203208, 1976.
 S. P. Meyn, 2007. Control Techniques for Complex Networks, Cambridge University Press, 2007. ISBN 9780521884419. Appendix contains abridged Meyn & Tweedie.
 S. M. Ross. 1983. Introduction to stochastic dynamic programming. Academic press
 X. Guo and O. HernándezLerma. ContinuousTime Markov Decision Processes, Springer, 2009.
 M. L. Puterman and Shin M. C. Modified Policy Iteration Algorithms for Discounted Markov Decision Problems, Management Science 24, 1978.
External links
 MDP Toolbox for Matlab  An excellent tutorial and Matlab toolbox for working with MDPs.
 Reinforcement Learning An Introduction by Richard S. Sutton and Andrew G. Barto
 SPUDD A structured MDP solver for download by Jesse Hoey
 Learning to Solve Markovian Decision Processes by Satinder P. Singh
Categories: Optimal decisions
 Dynamic programming
 Markov processes
 Stochastic control
Wikimedia Foundation. 2010.