マルコフ決定過程のソースを表示

'''マルコフ決定過程'''（マルコフけっていかてい、{{lang-en-short|Markov decision process; MDP}}）は、状態遷移が確率的に生じる動的システム（確率システム）の確率モデルであり、状態遷移が[[マルコフ性]]を満たすものをいう。
MDP は不確実性を伴う意思決定のモデリングにおける数学的枠組みとして、[[強化学習]]など[[動的計画法]]が適用される幅広い最適化問題の研究に活用されている。
MDP は少なくとも1950年代には知られていた{{sfn|Bellman|1957}}が、研究の中核は1960年に出版された Ronald A. Howard の "Dynamic Programming and Markov Processes" に起因する{{sfn|Howard|1960}}。
MDP は[[ロボット工学]]や自動制御、[[経済学]]、[[製造業]]を含む幅広い分野で用いられている。

== 概要 ==
[[ファイル:Markov_Decision_Process_example.png|右|サムネイル|400x400ピクセル|3つの状態と2つの行動をもつ簡単な MDP の例]]
マルコフ決定過程は離散時間における確率制御過程 (stochastic control process) である。
各時刻において過程 (process) はある'''状態''' (state) を取り、意思決定者 (decision maker) はその状態において利用可能な'''行動''' (action) を任意に選択する。
その後過程はランダムに新しい状態へと遷移し、その際に意思決定者は状態遷移に対応した'''報酬''' (reward) を受けとる。

遷移後の状態 <math>s'</math> 、および得られる報酬の値 <math>r</math> は現在の状態 <math>s</math> と行動 <math>a</math> のみに依存し、<math>s</math> と <math>a</math> が与えられたもとでそれより過去の状態および行動と条件付き独立となる。
言い換えると、マルコフ決定過程の状態遷移は[[マルコフ性]]を満たす。

マルコフ決定過程は[[マルコフ連鎖]]に（選択可能な）行動、および（行動を計画する動機を与える）報酬を追加し拡張したものであると解釈できる。
逆に言えば、各ステップにとる行動がそのステップにおける状態のみ依存するとき、マルコフ決定過程は等価なマルコフ連鎖に置き換えることが出来る。

== 定義 ==
MDP は4つの要素の組 <math display="inline">(S, A, T, R)</math>で表される。ここで各要素はそれぞれ次を意味する。
* <math>S</math> は状態の集合であり、状態空間 (state space) と呼ばれる。
* <math>A</math> は行動の集合であり、行動空間 (action space) と呼ばれる。特に、ある特定の状態 <math>s \in S</math> で取ることの出来る行動の集合を <math>A_s \subset A</math> によって表記する。
* <math>T</math> は遷移関数 (transition function) と呼ばれ、状態 <math>s \in S</math> において行動 <math>a \in A_s</math> を取った場合における状態 <math>s' \in S</math> への状態遷移確率を規定する実数値 <math>T(s, a, s') \in [0, 1]</math> を割り当てる。
* <math>R</math> は報酬関数 (reward function) と呼ばれ、状態 <math>s \in S</math> から <math>s' \in S</math> に行動 <math>a \in A_s</math> を伴い遷移することによって得られる即時報酬 (immediate reward) またはその期待値を表わす実数値 <math>R(s, a, s') \in \mathbb{R}</math> を規定する。

状態空間 <math>S</math> および行動空間 <math>A</math> は、有限集合の場合と（実数値すべての集合のような）無限集合の場合がある。特に <math>S</math> と <math>A</math> がともに有限集合であるような MDP は'''有限マルコフ決定過程''' (finite Markov decision process; finite MDP) と呼ばれる。状態空間が有限である場合においては、時刻 <math>t</math> における状態遷移確率は遷移関数によって
<math display="block">
  \Pr(s_{t+1} = s' \mid s_t = s, a_t = a) = T(s, a, s')
</math>
と記述することが出来る。
（状態空間が有限でない）より一般的なケースにおいては、任意の[[可測集合]] <math>S' \subset S</math> に対して状態遷移確率は
<math display="block">
  \Pr(s_{t+1} \in S' \mid s_t = s, a_t = a) = \int_{s' \in S'} T(s, a, s') \, ds'
</math>
によって記述される。

== 問題設定 ==
MDP における基本的な問題設定は、現在の状態が与えられたときに意思決定者の取る行動を既定する'''方策''' (policy) を求めることである。
具体的には、方策とは MDP のある状態 <math>s \in S</math> に対して定められる写像であり、次のように意思決定者の行動を規定する。
* '''決定論的方策''' (deterministic policy) : <math>a_t = \pi(s)</math>
* '''確率的方策''' (stochastic policy) : <math>\Pr(a_t = a \mid s_t = s) = \pi(s, a)</math>
 
方策 <math>\pi</math> は、ある目的関数 (objective function) を最大化するようなものが選ばれる。
目的関数には、現在時刻から始まり十分時間が経過した後の時点における報酬の累積値が用いられる。
いま、時刻 <math>\tau</math> において獲得する報酬を <math>r_{\tau} \in \mathbb{R}</math> と表記すると、累積報酬の実現値は <math>r_{t+1} + \cdots + r_{t+H}</math> で求められる。ただし、<math>H</math> は評価区間である。
ここで <math>r_{\tau}</math> は確率変数であり、その実現値は MDP によって規定される状態遷移確率や方策に基づいて確率的に変動する。
したがって、実際の目的関数としては累積報酬の期待値である
<math display="block">
  \mathbb{E}_{\pi} \bigg[ \sum_{\tau = 1}^{H} r_{t+\tau} \mid s_t = s \bigg]
</math>
が用いられる。
ただし、<math>\mathbb{E}_{\pi}[\cdot \mid s_t = s]</math> は方策 <math>\pi</math> のもとでの（条件付き）期待値を意味する。
このように有限の時間区間における累積報酬を評価する問題設定は'''有限時間区間''' (finite horizon) と呼ばれる。

行動が連続的に継続されるなどの場合においては、時間区間 <math>H</math> を適切に設定することが困難となる。
そのような場合、'''無限時間区間''' (infinite horizon) での問題設定が用いられる。
通常、期待値の収束性などの理由から有限時間区間における目的関数の <math>H \to \infty</math> における極限を考えることは不可能であり、典型的には次のように時刻が進むにつれて「割引される」報酬の累積値
<math display="block">
  \mathbb{E}_{\pi} \bigg[ \sum_{\tau=0}^{\infty} \gamma^\tau r_{t+\tau+1} \mid s_t = s \bigg]
</math>
が用いられる。
ここで <math>\gamma \in [0,1]</math> は'''割引率''' (discount rate) と呼ばれる値であり、現在の報酬と未来の報酬との間における重要度の差異を表している。

== 価値関数 ==
モデルのマルコフ性により、ある時刻における目的関数の値は方策 <math>\pi</math> を固定した下で状態 <math>s_t \in S</math> の関数とみなすことが出来る。
このような考えのもと、ある方策 <math>\pi \colon S \times A \to [0, 1]</math> に対して定義される
<math display="block">
  V^\pi(s) = \mathbb{E}_\pi \bigg[ \sum_{\tau=0}^{\infty} \gamma^{\tau} r_{t+\tau+1} \mid s_t = s \bigg]
</math> 
を '''状態価値関数''' (state-value function) と呼ぶ。
状態価値関数は、現在の状態 <math>s</math> からはじまり方策 <math>\pi</math> によって行動を決定することで得られる（割引された）累積報酬の値であると解釈できる。
この状態価値関数 <math>V^{\pi}(s)</math> は、次の整合性条件を満たす。
<math display="block">
  V^\pi(s) = \sum_{a \in A} \pi(s,a) \sum_{s' \in S} T(s, a, s') \Big( R(s, a, s') + \gamma V^\pi(s') \Big)
</math>

任意の <math>\pi'</math> および <math>s \in S</math> に対し <math>V^{\pi^*}(s) \geq V^{\pi'}(s)</math> を満たす方策 <math>\pi^*</math> を'''最適方策''' (optimal policy) と呼ぶ。
また、最適方策のもとでの状態価値関数 <math>V^*(s) = V^{\pi^{*}}(s)</math> は特に'''最適状態価値関数''' (optimal state-value function)と呼ばれる。
最適状態価値関数は、次の[[ベルマン方程式]]を満たす{{sfn|Suton|1998}}。
<math display="block">
  V^*(s) = \max_{a \in A} \sum_{s' \in S} T(s,a,s') \Big( R(s, a, s') + \gamma V^*(s') \Big)
</math>

== アルゴリズム ==
MDP は[[線型計画法|線形計画法]]または[[動的計画法]]で解くことができる。ここでは後者によるアプローチを示す．

=== 価値反復法 ===
価値反復法 (value iteration){{sfn|Bellman|1957}}は後ろ向き帰納法 (backward induction) とも呼ばれ、ベルマン方程式を満たす価値関数を繰り返し計算により求める。 [[ロイド・シャープレー]] が1953年に発表した{{仮リンク|確率ゲーム|en|stochastic games}}に関する論文{{sfn|Shapley|1953}}には価値反復法の特殊な場合が含まれるが、このことが認知されたのは後になってからである{{sfn|Kallenberg|2002}}．

ステップ <math>i</math>における価値関数の計算結果を <math>V_i(s)</math> と表記すると、価値反復法における更新式はつぎのように記述される:
<math display="block">
  V_{i+1}(s) \leftarrow \max_{a \in A_s} \sum_{s' \in S} T(s, a, s') \Big( R(s, a, s') + \gamma V_i(s') \Big) \quad \forall s \in S
</math>

上式をすべての状態において値が収束するまで繰り返したときの値を <math>V^\infty(s)</math> とし、最適方策  <math>\pi^*</math> を次式で求める。

<math display="block">
  \pi^*(s) \leftarrow \arg\max_{a \in A_s} \sum_{s' \in S} T(s, a, s') \Big( R(s, a, s') + \gamma V^\infty(s') \Big) \quad \forall s \in S
</math>

=== 方策反復法 ===
方策反復法 (policy iteration){{sfn|Howard|1960}}では、方策固定の下で行われる価値関数の更新 (policy evaluation) と、価値関数固定のもとで行われる方策の更新 (policy improvement) を交互に行うことで最適方策を求める。
# 次の線形方程式を解き、価値関数を更新する
#:<math>
  V^\pi(s) = \sum_{a \in A} \pi(s,a) \sum_{s' \in S} T(s, a, s') \Big( R(s, a, s') + \gamma V^\pi(s') \Big)
</math>
# 方策を次式で更新する
#:<math>
  \pi(s) \leftarrow \arg\max_{a \in A_s} \sum_{s' \in S} T(s, a, s') \Big( R(s, a, s') + \gamma V^\pi(s') \Big) \quad \forall s \in S
</math>

これらの操作を <math>\pi</math> がすべての状態に対し変化しなくなるまで繰り返すことで、最適方策を得る。
方策反復法は離散値を取る方策の値が変化しなくなるという明確な終了条件を持つため有限時間でアルゴリズムが終了するという利点を持つ。

== 拡張と一般化 ==
=== 部分観測マルコフ決定過程 ===
{{seealso|部分観測マルコフ決定過程}}

MDP では方策 <math>\pi(s)</math> を計算する際に現在の状態 <math>s</math> が既知であることを仮定している。
実際には状態観測に不確実性が伴う場合などこの仮定が成り立たない場合が多く、このような場合の一般化として'''[[部分観測マルコフ決定過程]]''' (Partially Observable Markov Decision Process; POMDP) が用いられる。

=== 強化学習 ===
{{seealso|強化学習|Q学習}}
状態遷移確率 <math>T(s,a,s')</math> や報酬関数 <math>R(s,a,s')</math> が未知の場合，環境との相互作用を通じてこれらの情報を得ながら行動を決定する必要がしばしば生じる．このような問題は[[強化学習]]の枠組みで議論される{{sfn|Sutton|Barto|1998}}．

強化学習における代表的な学習アルゴリズムは[[Q学習]]と呼ばれるものである。
Q学習では、'''行動価値関数''' (action-value function) と呼ばれる関数 <math>Q^\pi(s,a)</math> に着目する。ここで <math>Q^\pi(s,a)</math> は次のように定義される:
<math display="block">
  Q^\pi(s,a) = \mathbb{E}_\pi[\sum_{t=0}^\infty \gamma^t r_{t+1} | s_0=s, a_0=a]
</math>

いま，最適方策のもとでの行動価値関数 <math>Q^*(s,a) = \max_\pi Q^\pi(s,a) </math> は <math>V^*(s) = \max_{a} Q^*(s,a)</math> を満たす。
すなわち、<math>Q^*</math> を学習することができれば（モデルのパラメータを直接求めることなく）最適方策を獲得することができる。
Q学習では、各試行における遷移前後の状態と入力、および試行で得られる即時報酬の実現値をもとに <math>Q(s,a)</math> の値を逐次更新する。
実際の学習プロセスでは、すべての状態を十分サンプリングするため確率的なゆらぎを含むよう学習時の行動が選択される。

強化学習では最適化に必要なパラメータの学習を状態遷移確率・報酬関数を介することなくおこなうことが出来る（価値反復法や方策反復法ではそれらの明示的な仕様（各状態間の遷移可能性，報酬関数の関数形など）を与える必要がある）。
状態数（および行動の選択肢）が膨大な場合、強化学習はしばしばニューラルネットワークなどの関数近似と組み合わせられる。

=== 学習オートマトン  ===
機械学習理論における MDP のもう一つの応用は学習オートマトン (Learning Automata) と呼ばれる。
これは環境が確率的な挙動を示す場合における強化学習の一つでもある。
学習オートマトンに関する最初の詳細な論文は 1974 年に Narendra と Thathachar によりまとめられた{{sfn|Narendra|Thathachar|1974}}（そこでは有限状態オートマトンと明示的に記載されている）。
強化学習と同様，学習オートマトンのアルゴリズムも確率や報酬が未知の場合の問題を解くことができる。
Q学習の違いは，価値関数ではく学習の結果を探すために行動の確率を直接求めることである。
学習オートマトンは収束性が解析学の要領で厳密に証明されている{{sfn|Narendra|Thathachar|1989}}．

=== 制約付きマルコフ決定過程 ===
制約付きマルコフ決定過程 (Constrained Markov Decision Process; CMDP) はマルコフ決定過程の拡張である。
MDP と CMDP には3つの基本的な違いがある{{sfn|Altman|1999}}:
* ある行動をほかのものの代わりに適用した後で（複数の）コストが発生する
* CMDP は線形計画法のみで解くことが出来る（動的計画法を用いることはできない）
* 終端時刻における方策が初期状態に依存する
CMDP の応用例は数多く存在し、最近ではロボット工学におけるモーションプランニングに用いられている{{sfn|Feyzabadi|Carpin|2014}}。

== 参考文献 ==
{{reflist hide}}
{{refbegin|2}}

* {{cite journal
 |last=Bellman
 |first=R.
 |year=1957
 |title=A Markovian Decision Process
 |url=http://www.iumj.indiana.edu/IUMJ/FULLTEXT/1957/6/56038
 |journal=Journal of Mathematics and Mechanics
 |volume=6
 |ref=harv}}
* {{cite book
 |last=Howard
 |first=Ronald. A.
 |title=Dynamic Programming and Markov Processes
 |year=1960
 |publisher=The M.I.T. Press
 |ref=harv}}
* {{cite journal
 |last=Shapley
 |first=Lloyd.
 |year=1953
 |title=Stochastic Games
 |journal=Proceedings of National Academy of Science
 |volume=39
 |pages=1095–1100
 |ref=harv}}
* {{cite book
 |last=Kallenberg
 |first=Lodewijk.
 |title=Handbook of Markov decision processes: methods and applications
 |year=2002
 |publisher=Springer
 |ISBN=0-7923-7459-2
 |chapter=Finite state and action MDPs
 |editor-last1=Feinberg
 |editor-first1=Eugene A.
 |editor-last2=Shwartz
 |editor-first2=Adam
 |ref=harv}}
* {{cite book
 |last1=Sutton |first1=R. S.
 |last2=Barto |first2=A. G.
 |title=Reinforcement Learning: An Introduction
 |year=1998
 |url=http://webdocs.cs.ualberta.ca/~sutton/book/the-book.html
 |publisher=The MIT Press
 |location=Cambridge, MA
 |ref=harv}}
* {{cite journal
 |last1=Narendra |first1=K. S.
 |last2=Thathachar |first2=M. A. L.
 |year=1974
 |title=Learning Automata - A Survey
 |url=http://ieeexplore.ieee.org/stamp/stamp.jsp?arnumber=5408453
 |journal=IEEE Transactions on Systems, Man, and Cybernetics|volume=SMC-4
 |issue=4
 |pages=323–334
 |doi=10.1109/TSMC.1974.5408453
 |issn=0018-9472
 |ref=harv}}
* {{cite book
 |last1=Narendra |first1=Kumpati S.
 |last2=Thathachar |first2=Mandayam A. L.
 |title=Learning automata: An introduction
 |url=https://books.google.com/books?id=hHVQAAAAMAAJ
 |year=1989
 |publisher=Prentice Hall
 |isbn=9780134855585
 |ref=harv}}
* {{cite book
 |last=Altman |first=Eitan
 |title=Constrained Markov decision processes
 |year=1999
 |publisher=CRC Press
 |volume=7
 |ref=harv}}
* {{cite conference
 |last1=Feyzabadi |first1=S.
 |last2=Carpin |first2=S.
 |year=2014
 |pages=297,303
 |title=Risk-aware path planning using hierarchical constrained Markov Decision Processes
 |conference=IEEE International Conference
 |book-title=Automation Science and Engineering (CASE)
 |url=https://ieeexplore.ieee.org/abstract/document/6899341/
 |doi=10.1109/CoASE.2014.6899341}}
* {{cite journal
 | lang=ja
 | last=木村 |first=元
 | title=《第1回》強化学習の基礎
 | journal=計測と制御
 | publisher=計測自動制御学会
 | year=2013
 | volume=52
 | number=1
 | pages=72-77
 | url=https://doi.org/10.11499/sicejl.52.72
 | naid=10031140795
 |ref=harv}}
{{refend}}

==外部リンク==
* [http://www.cs.ualberta.ca/~sutton/book/ebook Reinforcement Learning] An Introduction by Richard S. Sutton and Andrew G. Barto
* [http://www.eecs.umich.edu/~baveja/Papers/Thesis.ps.gz Learning to Solve Markovian Decision Processes] by [http://www.eecs.umich.edu/~baveja/ Satinder P. Singh]
* [http://www.jstor.org/stable/3690147 Optimal Adaptive Policies for Markov Decision Processes] by Burnetas and Katehakis (1997)
* ソフトウェアパッケージ
** [http://www7.inra.fr/mia/T/MDPtoolbox/ MDP Toolbox for MATLAB, GNU Octave, Scilab and R] The Markov Decision Processes (MDP) Toolbox.
** [http://www.ai.mit.edu/~murphyk/Software/MDP/mdp.html MDP Toolbox for Matlab] - An excellent tutorial and Matlab toolbox for working with MDPs.
** [https://pypi.python.org/pypi/pymdptoolbox MDP Toolbox for Python] A package for solving MDPs
** [http://www.cs.uwaterloo.ca/~jhoey/research/spudd/index.php SPUDD] A structured MDP solver for download by Jesse Hoey

== 関連項目 ==
* [[動的計画法]]
* [[ベルマン方程式]]
* [[強化学習]]
* [[Q学習]]

{{確率論}}
{{Computer-stub}}

{{デフォルトソート:まるこふけつていかてい}}
[[Category:理論計算機科学]]
[[Category:論理学]]
[[Category:アンドレイ・マルコフ]]
[[Category:数学に関する記事]]
[[Category:エポニム]]