LINEAR-ALGEBRA Β· Interactive Practice | Unit 23 Β· Video 2

Why Every Markov Chain Forgets Where It Started: The Ξ»=1 Attractor

IKey Formulas

Formula Name Description
uk+1=A uku_{k+1} = A\, u_k Markov iteration One time step of the chain
uk=c1Ξ»1kx1+c2Ξ»2kx2+β‹―+cnΞ»nkxnu_k = c_1 \lambda_1^k x_1 + c_2 \lambda_2^k x_2 + \cdots + c_n \lambda_n^k x_n Spectral master formula Powers of AA in the eigenvector basis
Ξ»1=1,  ∣λi∣<1\lambda_1 = 1,\ \ \lvert\lambda_i\rvert < 1 for iβ‰₯2i \ge 2 Markov spectrum Always true for column-stochastic AA
ukβ†’c1x1u_k \to c_1 x_1 as kβ†’βˆžk \to \infty Steady state Only the Ξ»=1\lambda=1 direction survives

IIVisualization 1 β€” Forgetting the Starting Split

Each step fraction pp leaves City A and fraction qq leaves City B; run the chain forward and watch the split settle.

πŸ’‘ For A=[1βˆ’pqp1βˆ’q]A = \begin{bmatrix} 1-p & q \\ p & 1-q \end{bmatrix} the steady state is 1p+q(q, p)\tfrac{1}{p+q}(q,\,p) β€” fixed by the rates alone, never by the start.

IIIVisualization 2 β€” How Fast a Chain Forgets

Every non-dominant term decays like ∣λ2∣k|\lambda_2|^k; the closer ∣λ2∣|\lambda_2| sits to 1, the slower the chain forgets.

πŸ’‘ The steps needed to shrink to 1% is the chain's mixing time; the ratio ∣λ2∣|\lambda_2| that sets it is the mixing rate.

IVVisualization 3 β€” All Roads to One Steady State

Six different starting splits, all on the line a+b=1a+b=1 β€” every trajectory slides to the same steady-state point.

VQuiz Questions

Question 1

For any column-stochastic Markov matrix AA, what is the largest eigenvalue Ξ»1\lambda_1?

βœ… Correct! Every column-stochastic matrix has Ξ»1=1\lambda_1 = 1 β€” guaranteed by the columns summing to 1.

❌ Not quite. The columns-sum-to-1 property forces a specific eigenvalue. Think about what 1TA\mathbf{1}^T A equals.

Show solution

Solution:

Every column-stochastic matrix has columns summing to 1. This means the all-ones row vector 1T\mathbf{1}^T satisfies 1TA=1T\mathbf{1}^T A = \mathbf{1}^T, so 1T\mathbf{1}^T is a left eigenvector with eigenvalue 11. Since AA and ATA^T share eigenvalues, Ξ»=1\lambda = 1 is also an eigenvalue of AA.

Combined with the Perron–Frobenius bound ∣λiβˆ£β‰€1|\lambda_i| \le 1 for all eigenvalues of a stochastic matrix, we conclude:

Ξ»1=1Β exactly,Β always.\lambda_1 = 1 \text{ exactly, always.}

Question 2

A Markov chain has eigenvalues Ξ»1=1\lambda_1 = 1, Ξ»2=0.95\lambda_2 = 0.95, and Ξ»3=0.2\lambda_3 = 0.2. Approximately how many steps kk are needed for the non-dominant components to decay by a factor of about 100 (i.e., ∣λ2∣kβ‰ˆ0.01|\lambda_2|^k \approx 0.01)?

βœ… Correct! ln⁑(0.01)/ln⁑(0.95)β‰ˆ90\ln(0.01)/\ln(0.95) \approx 90. The slowest-decaying eigenvalue Ξ»2=0.95\lambda_2 = 0.95 sets the mixing time.

❌ Not quite. Convergence is controlled by the second-largest eigenvalue. Solve 0.95kβ‰ˆ0.010.95^k \approx 0.01.

Show solution

Solution:

The convergence rate is governed by the second-largest eigenvalue magnitude, ∣λ2∣=0.95|\lambda_2| = 0.95. We want:

0.95kβ‰ˆ0.010.95^k \approx 0.01

Taking logs: kβ‰ˆln⁑(0.01)ln⁑(0.95)=βˆ’4.605βˆ’0.0513β‰ˆ89.8k \approx \dfrac{\ln(0.01)}{\ln(0.95)} = \dfrac{-4.605}{-0.0513} \approx 89.8.

So about 90 steps. Note that Ξ»3=0.2\lambda_3 = 0.2 decays much faster (0.2900.2^{90} is astronomically tiny), so it's irrelevant β€” the slowest non-dominant eigenvalue dominates the mixing time.

Question 3

True or False: For a (well-behaved) Markov chain, the steady state u∞u_\infty depends on the starting distribution u0u_0.

βœ… Correct! The steady state c1x1c_1 x_1 is determined by the matrix's Ξ»=1\lambda=1 eigenvector, and the probability constraint pins c1c_1 β€” independent of u0u_0.

❌ Not quite. Look back at Visualization 3: six different starting splits all flowed to the same red star.

Show solution

Solution:

The master formula gives uk=c1β‹…1kβ‹…x1+c2Ξ»2kx2+β‹―u_k = c_1 \cdot 1^k \cdot x_1 + c_2 \lambda_2^k x_2 + \cdots. As kβ†’βˆžk \to \infty, only the first term survives, so:

u∞=c1x1.u_\infty = c_1 x_1.

The direction x1x_1 is fixed by the matrix AA β€” it does not depend on u0u_0 at all. The scalar c1c_1 does depend on u0u_0 in principle, but for any probability vector u0u_0 (entries non-negative, summing to 1), the constraint that u∞u_\infty also sums to 1 pins down c1c_1 uniquely.

Result: every probability starting distribution converges to the same steady state. The chain forgets where it began.

(Caveat: this assumes Ξ»=1\lambda = 1 has multiplicity one and the chain is irreducible/aperiodic β€” the well-behaved case.)

Question 4

Consider the column-stochastic matrix A=[0.70.40.30.6].A = \begin{bmatrix} 0.7 & 0.4 \\ 0.3 & 0.6 \end{bmatrix}.

What is the steady-state distribution (the eigenvector for Ξ»=1\lambda = 1, normalized so its entries sum to 1)?

βœ… Correct! Solving (Aβˆ’I)x=0(A-I)x = 0 and normalizing gives (4/7,3/7)(4/7, 3/7).

❌ Not quite. (0.5,0.5)(0.5, 0.5) would only be the steady state if the matrix were symmetric (e.g., p=qp = q). Here pβ‰ qp \ne q.

❌ Not quite. Reading the diagonal isn't the same as solving Ax=xAx = x. Set up the eigenvalue equation.

❌ Not quite. You may have swapped the components. Solve (Aβˆ’I)x=0(A-I)x = 0 carefully.

Show solution

Solution:

We need Ax=xA x = x, i.e., (Aβˆ’I)x=0(A - I)x = 0:

Aβˆ’I=[βˆ’0.30.40.3βˆ’0.4].A - I = \begin{bmatrix} -0.3 & 0.4 \\ 0.3 & -0.4 \end{bmatrix}.

The first row gives βˆ’0.3x1+0.4x2=0-0.3 x_1 + 0.4 x_2 = 0, so x2=34x1x_2 = \tfrac{3}{4} x_1, or (x1,x2)∝(4,3)(x_1, x_2) \propto (4, 3).

Normalizing so the entries sum to 1:

x1=(47,Β 37)β‰ˆ(0.571,Β 0.429).x_1 = \left(\frac{4}{7},\ \frac{3}{7}\right) \approx (0.571,\ 0.429).

Shortcut for 2Γ—2 column-stochastic [1βˆ’pqp1βˆ’q]\begin{bmatrix} 1-p & q \\ p & 1-q \end{bmatrix}: steady state is 1p+q(q,p)\dfrac{1}{p+q}(q, p). Here p=0.3p = 0.3, q=0.4q = 0.4, giving (0.4,0.3)/0.7=(4/7,3/7)(0.4, 0.3)/0.7 = (4/7, 3/7).

Solved: 0 / 4