27차시 · 정의와 표기

확률벡터와 전이행렬

유한한 상태를 \(1,\ldots,d\)라 하자. 상태분포 \(x_n\)은 \(n\)단계에 각 상태에 있을 확률을 담은 열벡터이며 \(x_{n,i}\ge0\), \(\sum_i x_{n,i}=1\)이다.

\[P_{ij}=\Pr(\text{다음 상태 }i\mid\text{현재 상태 }j),\quad P_{ij}\ge0,\quad\sum_iP_{ij}=1\]

다음 상태의 조건부 확률이 과거 전체가 아니라 현재 상태만으로 정해지는 모형을 마르코프 연쇄라고 한다. 여기서는 전이행렬 \(P\)가 매 단계 같다고 가정하여 \(x_{n+1}=Px_n\)으로 계산한다. 열은 출발, 행은 도착이다.

27차시 · 계산과 논증

열의 합이 1이라는 조건은 왜 필요한가?

① \(P\)가 전이행렬이고 \(x\)가 확률벡터일 때 \(Px\)도 확률벡터임을 증명하라. ② 행의 합만 1인 다음 행렬에 \(x=(1,0)^T\)를 넣고 문제를 설명하라.

\[R=\begin{bmatrix}.8&.2\\.3&.7\end{bmatrix}\]
풀이와 판단 근거

모든 항이 음이 아니므로 \((Px)_i\ge0\)이다. 또한

\[\sum_i(Px)_i=\sum_j\left(\sum_iP_{ij}\right)x_j=\sum_jx_j=1\]

\(Rx=(.8,.3)^T\)는 합이 1.1이므로 확률벡터가 아니다. 행벡터로 곱하는 약속에서는 행의 합을 1로 두지만, 이 수업의 열벡터 약속에 그대로 적용할 수 없다.

27차시 · 마르코프 연쇄

같은 전이 규칙을 반복하면 상태분포는 수렴할까?

다음 상태는 현재 상태에만 좌우되고 전이확률 \(P\)는 매 단계 같다고 가정합니다.

\(\mathbf{x}_{n+1}=P\mathbf{x}_n\)을 반복하면서 ‘변하지 않는 분포’와 ‘그 분포로 가까워지는 현상’을 구분해 봅시다.

A 마을

100%

B 마을

0%

27차시 · 정의와 표기

정지분포와 수렴의 구별

전이행렬 \(P\)에 대해 \(P\pi=\pi\)를 만족하는 확률벡터 \(\pi\)를 정지분포라고 한다. 고윳값 1의 고유벡터 중 성분이 음이 아니고 합이 1인 벡터이다. \(\mathbf1\)은 성분이 모두 1인 열벡터이므로 \(\mathbf1^T\pi\)는 성분의 합이다.

\[P\pi=\pi,\qquad \pi_i\ge0,\qquad \mathbf1^T\pi=1\]

정지분포에서 출발하면 분포가 그대로이다. 모든 초기분포 \(x_0\)에 대해 \(P^nx_0\to\pi\)가 되는지는 다른 질문이다. 다음 계산에서 정지분포를 먼저 찾고 반복의 수렴을 따로 확인한다.

27차시 · 유도와 계산

정지분포를 방정식으로 계산

\(P\pi=\pi\)와 확률의 합 1을 함께 사용합니다.

\[.8a+.3b=a,\quad a+b=1\Rightarrow .2a=.3b\Rightarrow(a,b)=(.6,.4)\]

정지분포는 움직이지 않는 분포입니다. 모든 초기분포가 여기에 수렴하는지는 주기성과 다른 고윳값을 별도로 살펴야 합니다.

직접 반복

한 단계씩 눌러 비율의 변화가 어떻게 줄어드는지 살펴봅시다

열 = 출발 · 행 = 도착
AB
A0.80.3
B0.20.7
0일
A · 100.0%B · 0.0%

초기분포 $\mathbf{x}_0=(1,0)^\mathsf T$

매일 이동률: A→B 20% · B→A 30%

$A_n$$B_n$

27차시 · 정의와 표기

기약성·주기·수렴 조건

어느 상태에서 출발해도 다른 모든 상태에 유한 단계 안에 양의 확률로 갈 수 있으면 기약이라고 한다. 상태 \(i\)의 주기는 돌아올 수 있는 단계 수들의 최대공약수이다.

\[d(i)=\gcd\{n\ge1:(P^n)_{ii}>0\}\]

기약 연쇄에서 이 주기가 1이면 비주기라고 한다. 유한 상태의 기약·비주기 연쇄는 유일한 정지분포를 가지며 모든 초기분포가 그곳으로 수렴한다.

모든 \(P_{ij}>0\)이면 두 조건이 성립한다. 앞의 2상태 예에서는 \(x_{n,A}-.6=.5^n(x_{0,A}-.6)\)이므로 오차가 매번 절반이 된다.

27차시 · 정의와 표기

페이지랭크의 전이 규칙

링크를 따라 이동하는 행렬을 \(M\)이라 한다. 나가는 링크가 없는 페이지의 열은 미리 확률벡터 \(v\)로 채워 모든 열의 합을 1로 만든다. \(\mathbf1\)은 성분이 모두 1인 열벡터이다.

\[G=\alpha M+(1-\alpha)v\mathbf1^T,\quad 0<\alpha<1,\quad v_i>0,\quad\mathbf1^Tv=1\]

확률 \(\alpha\)로 링크를 따르고 확률 \(1-\alpha\)로 \(v\)에 따라 새 페이지를 선택한다. \(v\mathbf1^T\)의 모든 열은 \(v\)이다. 페이지랭크는 \(G\pi=\pi\)인 정지분포 \(\pi\)의 각 성분이다.

\(G_{ij}\ge(1-\alpha)v_i>0\)이므로 유일한 정지분포로 수렴한다. \(v_i=0\)을 허용하면 이 양수성 근거를 그대로 쓸 수 없다.

27차시 · 계산과 논증

주기적인 왕복에 임의 이동을 섞기

\[M=\begin{bmatrix}0&1\\1&0\end{bmatrix},\quad v=\binom{1/2}{1/2},\quad x_0=\binom10\]

① \(M\)의 정지분포와 \(M^nx_0\)의 거동을 구하라. ② \(\alpha=.8\)인 \(G\)를 구하고, 차이 \(d_n=x_{n,1}-x_{n,2}\)의 점화식으로 수렴을 설명하라.

풀이와 판단 근거

정지분포는 \((1/2,1/2)^T\)이지만 \(M^nx_0\)는 두 상태를 번갈아 지나 수렴하지 않는다.

\[G=\begin{bmatrix}.1&.9\\.9&.1\end{bmatrix},\qquad d_{n+1}=-.8d_n,\quad d_n=(-.8)^n\]

확률의 합은 1이고 차이는 0으로 가므로 두 성분 모두 1/2로 간다. 부호가 교대하여 정지분포 양쪽을 번갈아 지나지만 진폭은 줄어든다. \(\alpha=1\)이면 이 감쇠가 사라진다.

핵심 정리

확률 보존과 상태분포의 수렴은 별개다

정지분포를 찾은 뒤, 기약성과 비주기성까지 살펴봅시다.

전이행렬정지분포기약성·비주기성