Skip to content
Date: Oct 2, 2026
Builds on Chapter 3.

A machine is run week after week with no planned end. Each week its operator decides whether to keep running it or to pay for an overhaul. There is no last stage from which to start a backward recursion, and the total cost over an unbounded future is infinite for every sensible policy, so “minimize the total cost” does not even pose a question.

Two changes rescue the problem. The first is to discount: a cost incurred \(k\) weeks from now is multiplied by \(\alpha^{k}\) for a fixed \(\alpha\) strictly between \(0\) and \(1\) , so that the infinite sum converges. The second is a change of viewpoint. With a finite horizon the unknowns were \(N+1\) functions \(J_0,\dots,J_N\) linked by a recursion. With an infinite horizon and a model that does not change with time, the problem looks the same from every stage, so there is a single unknown function \(J^\star\) , and the recursion turns into an equation that \(J^\star\) must satisfy. This chapter proves that the equation has exactly one solution, that the solution is the optimal cost, and that an optimal policy can be read from it and uses the same rule at every stage. The whole argument rests on one inequality, the contraction property, which is why the discount factor matters so much.

4.1 A model with no last stage

Definition 4.1 (discounted Markov decision process). A discounted Markov decision process consists of a finite state set \(S=\{1,\dots,n\}\) ; for each state \(i\) a finite nonempty set \(U(i)\) of controls; transition probabilities \(p_{ij}(u)\) , the probability that the next state is \(j\) when control \(u\) is used in state \(i\) , with \(\sum_jp_{ij}(u)=1\) ; expected stage costs \(g(i,u)\) ; and a discount factor \(\alpha\in(0,1)\) . A policy is a sequence \(\pi=(\mu_0,\mu_1,\dots)\) with \(\mu_k(i)\in U(i)\) , and its cost from state \(i\) is

\[ J_\pi(i)=\lim_{N\to\infty}E\Big[\sum_{k=0}^{N-1}\alpha^{k}g\big(x_k,\mu_k(x_k)\big)\ \Big|\ x_0=i\Big]. \]

A policy is stationary if \(\mu_k=\mu\) for all \(k\) ; it is then identified with the single function \(\mu\) , and its cost is written \(J_\mu\) . The optimal cost is \(J^\star(i)=\inf_\pi J_\pi(i)\) .

This is the transition-probability form of the model of Chapter 1, with nothing depending on the stage. A function on \(S\) is a vector in \(\mathbb R^{n}\) , and inequalities between such vectors are meant componentwise. Under a stationary policy the state is a Markov chain with transition matrix \(P_\mu\) , whose entry in row \(i\) and column \(j\) is \(p_{ij}(\mu(i))\) ; the \(k\) -step transition probabilities of that chain are the entries of \(P_\mu^{k}\) , a fact from Stochastic Processes used below without further comment.

The limit in the definition exists. Let \(M=\max_{i,u}\lvert g(i,u)\rvert\) . The \(k\) -th term of the sum is at most \(\alpha^{k}M\) in absolute value, so the series converges absolutely and

\[ \lvert J_\pi(i)\rvert\le\frac{M}{1-\alpha},\qquad \Big\lvert J_\pi(i)-E\Big[\sum_{k=0}^{N-1}\alpha^{k}g\Big]\Big\rvert\le\frac{\alpha^{N}M}{1-\alpha}. \]

The second bound says that the infinite-horizon cost differs from the \(N\) -stage cost by a tail that shrinks geometrically. All of the convergence in this chapter comes from that one fact.

4.2 Two operators

The finite-horizon recursion, written in the discounted form of Exercise 2.13, applies the same transformation at every stage. Give the transformation a name.

Definition 4.2 (Bellman operators). For a function \(J:S\to\mathbb R\) and a stationary policy \(\mu\) , define new functions \(T_\mu J\) and \(TJ\) by

\[ (T_\mu J)(i)=g\big(i,\mu(i)\big)+\alpha\sum_{j}p_{ij}\big(\mu(i)\big)J(j),\qquad (TJ)(i)=\min_{u\in U(i)}\Big[g(i,u)+\alpha\sum_jp_{ij}(u)J(j)\Big]. \]

\((T_\mu J)(i)\) is the cost of a one-stage problem from state \(i\) in which \(\mu\) is used and a terminal cost \(\alpha J\) is charged at the next state; \((TJ)(i)\) is the optimal cost of that one-stage problem. In vector form \(T_\mu J=g_\mu+\alpha P_\mu J\) , an affine map. The operator \(T\) is not affine, because of the minimum, and that is the source of every difficulty that the next chapter’s algorithms have to deal with.

Applying the operators repeatedly builds longer horizons. By Lemma 2.1, \(T_{\mu_0}T_{\mu_1}\cdots T_{\mu_{N-1}}J\) is the cost of using \(\mu_0,\dots,\mu_{N-1}\) over \(N\) stages with terminal cost \(\alpha^{N}J\) ; and by Theorem 2.2, \(T^{N}J\) is the optimal cost of the \(N\) -stage problem with that terminal cost. In particular \(T^{N}0\) is the optimal \(N\) -stage discounted cost.

Example 4.1: Two policies for a wearing machine

A machine is good (state \(1\) ) or worn (state \(2\) ). A good machine runs at no cost and becomes worn with probability \(0.3\) each week. A worn machine can be run at a cost of \(2\) per week, staying worn, or overhauled at a cost of \(5\) , which makes it good at the start of next week. Take \(\alpha=0.9\) . In state \(1\) there is one control; in state \(2\) there are two, so there are two stationary policies.

Never overhaul. The cost vector satisfies \(J=T_\mu J\) :

\[ J(1)=0.9\big(0.7J(1)+0.3J(2)\big),\qquad J(2)=2+0.9J(2). \]

The second equation gives \(J(2)=20\) , and the first then gives \(0.37J(1)=0.27\cdot20\) , so \(J(1)=\tfrac{540}{37}\approx14.595\) .

Overhaul when worn. Now

\[ J(1)=0.9\big(0.7J(1)+0.3J(2)\big),\qquad J(2)=5+0.9J(1). \]

Substituting the second into the first, \(0.37J(1)=0.27\,(5+0.9J(1))\) , so \(0.127J(1)=1.35\) and

\[ J(1)=\tfrac{1350}{127}\approx10.630,\qquad J(2)=\tfrac{1850}{127}\approx14.567 . \]

Overhauling is better from both states, by about \(4\) from the good state and \(5.4\) from the worn one. That an identity of the form \(J=T_\mu J\) determines the cost of a policy, as these calculations assumed, is part of Theorem 4.5.

Two properties of the operators carry the whole theory.

Lemma 4.3 (monotonicity and constant shift). Let \(J\) and \(J'\) be functions on \(S\) and \(c\) a real number, and let \(\mathbf 1\) denote the function equal to \(1\) everywhere. If \(J\le J'\) then \(TJ\le TJ'\) and \(T_\mu J\le T_\mu J'\) . Moreover \(T(J+c\mathbf 1)=TJ+\alpha c\mathbf 1\) and \(T_\mu(J+c\mathbf 1)=T_\mu J+\alpha c\mathbf 1\) .

Proof. For fixed \(i\) and \(u\) the bracket \(g(i,u)+\alpha\sum_jp_{ij}(u)J(j)\) is nondecreasing in each \(J(j)\) , since the weights \(\alpha p_{ij}(u)\) are nonnegative; and replacing \(J\) by \(J+c\mathbf 1\) adds \(\alpha c\sum_jp_{ij}(u)=\alpha c\) to it. Both properties pass to \(T_\mu\) , which is one such bracket, and to \(T\) , because a minimum of brackets inherits an inequality that holds for each bracket and is shifted by a constant that is added to every bracket.

Monotonicity is the property already singled out in Chapter 2: lowering the cost of the future lowers the cost of the present. The shift property is new, and it is where \(\alpha\) enters: a uniform error of size \(c\) in the terminal cost becomes an error of size \(\alpha c\) after one application.

4.3 The contraction property

Measure the distance between two functions on \(S\) by the maximum norm \(\lVert J\rVert=\max_i\lvert J(i)\rvert\) .

Proposition 4.4 (contraction). For all functions \(J\) , \(J'\) on \(S\) and every stationary policy \(\mu\) ,

\[ \lVert TJ-TJ'\rVert\le\alpha\,\lVert J-J'\rVert,\qquad \lVert T_\mu J-T_\mu J'\rVert\le\alpha\,\lVert J-J'\rVert . \]

Proof. Let \(c=\lVert J-J'\rVert\) . Then \(J\le J'+c\mathbf 1\) , so by Lemma 4.3, \(TJ\le T(J'+c\mathbf 1)=TJ'+\alpha c\mathbf 1\) . Exchanging the roles of \(J\) and \(J'\) gives \(TJ'\le TJ+\alpha c\mathbf 1\) . Together these say \(\lvert(TJ)(i)-(TJ')(i)\rvert\le\alpha c\) for every \(i\) . The same argument applies to \(T_\mu\) .

Example 4.2: The contraction at work

For the machine of Example 4.1, \((TJ)(1)=0.9\,(0.7J(1)+0.3J(2))\) and \((TJ)(2)=\min\{2+0.9J(2),\ 5+0.9J(1)\}\) .

Take \(J=(0,0)\) and \(J'=(10,20)\) , at distance \(\lVert J-J'\rVert=20\) . Then \(TJ=(0,\ \min\{2,5\})=(0,2)\) and

\[ TJ'=\big(0.9\cdot13,\ \min\{20,\,14\}\big)=(11.7,\ 14), \]

so \(\lVert TJ-TJ'\rVert=\max\{11.7,\,12\}=12\) , comfortably below the bound \(0.9\cdot20=18\) .

Iterating from \(J_0=(0,0)\) gives \(J_1=(0,2)\) , \(J_2=(0.54,\,3.8)\) , \(J_3=(1.3662,\,5.42)\) , \(J_4=(2.3241,\,6.2296)\) . The minimizing control in state \(2\) is “run” in computing \(J_1\) , \(J_2\) , and \(J_3\) , and “overhaul” from \(J_4\) on. The distance from \(J_k\) to the vector \((10.630,14.567)\) of Example 4.1 is \(12.567\) , \(10.767\) , \(9.264\) , \(8.337\) for \(k=1,2,3,4\) : successive ratios \(0.857\) , \(0.860\) , \(0.900\) , never above \(\alpha\) . Reaching a distance of \(10^{-6}\) takes \(156\) iterations.

A contraction on \(\mathbb R^{n}\) has exactly one fixed point, and iteration finds it.

Theorem 4.5 (fixed points of the Bellman operators). The operator \(T\) has a unique fixed point \(\bar J\) , and for every function \(J\) on \(S\) ,

\[ \lVert T^{k}J-\bar J\rVert\le\alpha^{k}\lVert J-\bar J\rVert\qquad(k=0,1,2,\dots). \]

The same holds for each \(T_\mu\) , whose unique fixed point is the cost \(J_\mu\) of the stationary policy \(\mu\) ; equivalently \(J_\mu\) is the unique solution of the linear system \((I-\alpha P_\mu)J=g_\mu\) .

Proof. Let \(J_k=T^{k}J\) . By Proposition 4.4, \(\lVert J_{k+1}-J_k\rVert\le\alpha^{k}\lVert J_1-J_0\rVert\) , so for \(m\) greater than \(k\) ,

\[ \lVert J_m-J_k\rVert\le\big(\alpha^{k}+\alpha^{k+1}+\dots+\alpha^{m-1}\big)\lVert J_1-J_0\rVert\le\frac{\alpha^{k}}{1-\alpha}\lVert J_1-J_0\rVert . \]

Each coordinate sequence \(J_k(i)\) is therefore a Cauchy sequence of real numbers and converges; call the limit \(\bar J(i)\) . Since \(\lVert TJ_k-T\bar J\rVert\le\alpha\lVert J_k-\bar J\rVert\to0\) , the sequence \(J_{k+1}=TJ_k\) converges to \(T\bar J\) as well as to \(\bar J\) , so \(T\bar J=\bar J\) . If \(J'\) is another fixed point, then \(\lVert\bar J-J'\rVert=\lVert T\bar J-TJ'\rVert\le\alpha\lVert\bar J-J'\rVert\) , which forces \(\lVert\bar J-J'\rVert=0\) . The rate follows from \(T^{k}J-\bar J=T^{k}J-T^{k}\bar J\) and \(k\) applications of the proposition.

The argument used only the contraction property, so it applies to \(T_\mu\) . Its fixed point is \(J_\mu\) because \(T_\mu^{N}0\) is the \(N\) -stage cost of \(\mu\) , which converges to \(J_\mu\) by the tail bound of Section 4.1 and to the fixed point by what was just proved. Finally \(J=T_\mu J\) is the system \(J=g_\mu+\alpha P_\mu J\) .

The step from “Cauchy” to “convergent” is the completeness of the real numbers, from Real Analysis (math.real-analysis.cauchy-completeness). Nothing else from outside is needed.

4.4 Bellman’s equation

Theorem 4.6 (optimality in discounted problems). In a discounted Markov decision process: the optimal cost \(J^\star\) is the unique solution of Bellman’s equation \(J=TJ\) , that is,

\[ J^\star(i)=\min_{u\in U(i)}\Big[g(i,u)+\alpha\sum_jp_{ij}(u)J^\star(j)\Big]\qquad\text{for all } i\in S; \]

a stationary policy \(\mu\) is optimal if and only if \(\mu(i)\) attains the minimum in this equation for every \(i\) , that is, \(T_\mu J^\star=TJ^\star\) ; and consequently an optimal stationary policy exists.

Proof. Let \(\bar J\) be the fixed point of \(T\) from Theorem 4.5. First, no policy does better than \(\bar J\) . For any policy \(\pi=(\mu_0,\mu_1,\dots)\) the \(N\) -stage cost is \(T_{\mu_0}\cdots T_{\mu_{N-1}}0\) . Since \(T_\mu J\ge TJ\) for every \(J\) , and both operators are monotone, an induction on \(N\) gives \(T_{\mu_0}\cdots T_{\mu_{N-1}}0\ge T^{N}0\) . Letting \(N\to\infty\) , the left side tends to \(J_\pi\) and the right side to \(\bar J\) , so \(J_\pi\ge\bar J\) . Hence \(J^\star\ge\bar J\) .

Second, some stationary policy attains \(\bar J\) . Since each \(U(i)\) is finite, there is a \(\mu\) with \(\mu(i)\) attaining the minimum in \((T\bar J)(i)\) for every \(i\) , that is, \(T_\mu\bar J=T\bar J=\bar J\) . So \(\bar J\) is a fixed point of \(T_\mu\) , and by Theorem 4.5 it equals \(J_\mu\) . Thus \(J_\mu=\bar J\le J^\star\le J_\mu\) , which proves \(J^\star=\bar J\) and that \(\mu\) is optimal.

For the characterization: if \(T_\mu J^\star=TJ^\star=J^\star\) then \(J_\mu=J^\star\) as just shown. Conversely, if \(J_\mu=J^\star\) then \(T_\mu J^\star=T_\mu J_\mu=J_\mu=J^\star=TJ^\star\) .

The theorem turns an optimization over infinitely many policies into \(n\) equations in \(n\) unknowns. It also gives a test that needs no algorithm: to prove a guessed policy optimal, compute its cost by solving a linear system and check whether that cost satisfies Bellman’s equation.

Example 4.3: Bellman’s equation as a test

Test the two policies of Example 4.1. The first component of Bellman’s equation involves no choice and holds for both cost vectors by construction, so only state \(2\) needs checking.

For “overhaul when worn”, with \(J=(10.630,\,14.567)\) : running costs \(2+0.9\cdot14.567=15.110\) and overhauling costs \(5+0.9\cdot10.630=14.567\) . The minimum is \(14.567=J(2)\) , attained by overhauling. Bellman’s equation holds, so this vector is \(J^\star\) and the policy is optimal.

For “never overhaul”, with \(J=(14.595,\,20)\) : running costs \(2+0.9\cdot20=20\) and overhauling costs \(5+0.9\cdot14.595=18.135\) . The minimum is \(18.135\) , which is not \(J(2)=20\) . Bellman’s equation fails, so this policy is not optimal. The test even says how to improve it: overhaul.

Common mistake. Confusing the two equations \(J=T_\mu J\) and \(J=TJ\) . Every stationary policy satisfies the first with its own cost vector; that is how \(J_\mu\) was computed. Only an optimal policy’s cost satisfies the second. An argument that runs “my policy’s cost satisfies its own equation, so it satisfies Bellman’s equation, so it is optimal” has replaced the minimum over all controls by the one control the policy happens to use. The check in Example 4.3 must compare every control in every state.

4.5 What the discount factor does

The discount factor is a statement about patience, and the optimal policy depends on it.

Example 4.4: Patience changes the policy

Solve the machine of Example 4.1 for a general \(\alpha\) . Under “never overhaul”, \(J(2)=\dfrac{2}{1-\alpha}\) . Under “overhaul when worn”, eliminating \(J(1)=\dfrac{0.3\alpha J(2)}{1-0.7\alpha}\) from \(J(2)=5+\alpha J(1)\) gives

\[ J(2)=\frac{5\,(10-7\alpha)}{(1-\alpha)(10+3\alpha)} . \]

Overhauling is the better policy when this is the smaller of the two, that is, when \(5(10-7\alpha)\le2(10+3\alpha)\) , or \(\alpha\ge\tfrac{30}{41}\approx0.732\) . A decision maker with \(\alpha\) below \(\tfrac{30}{41}\) prefers to pay \(2\) every week for ever; one with \(\alpha\) above it pays \(5\) now for cheaper weeks later. At \(\alpha=0.5\) the costs from the worn state are \(4\) against \(\tfrac{130}{23}\approx5.652\) ; at \(\alpha=0.9\) they are \(20\) against \(14.567\) .

Two curves of normalized cost from the worn state against the discount factor: a horizontal line at 2 for the policy that never overhauls, and a decreasing curve, from about 3.6 at a discount factor of 0.3 down to about 1.15 near 1, for the policy that overhauls when worn; the curves cross at a discount factor of about 0.73, and the lower envelope is emphasized.
Figure 4.1. The quantity plotted is the cost from the worn state multiplied by one minus the discount factor, a cost per week. The optimal cost is the lower envelope, and the optimal policy changes where the curves cross: patience makes the overhaul worth its price.

Multiplying by \(1-\alpha\) , as in the figure, converts a discounted total into a weighted average cost per stage, which makes different discount factors comparable. As \(\alpha\to1\) the normalized cost of the overhaul policy tends to \(\tfrac{15}{13}\) , its long-run average cost per week.

There is a second reading of \(\alpha\) that involves no impatience at all.

Proposition 4.7 (discounting as random termination). Suppose that after each stage the process is terminated with probability \(1-\alpha\) , independently of the states and controls, and that no cost is incurred after termination. Then for every policy the expected total undiscounted cost up to termination equals the discounted cost \(J_\pi(i)\) .

Proof. Stage \(k\) is reached if and only if the process survives \(k\) times, which has probability \(\alpha^{k}\) and is independent of the state and control at stage \(k\) . The expected cost contributed by stage \(k\) is therefore \(\alpha^{k}E[g(x_k,\mu_k(x_k))]\) , and summing over \(k\) gives \(J_\pi(i)\) .

So a machine that is scrapped each week with probability \(0.1\) for outside reasons, and whose owner does not discount at all, is governed by the model with \(\alpha=0.9\) . The expected remaining lifetime is \(\dfrac{1}{1-\alpha}\) stages, which gives a rule of thumb: a discount factor \(\alpha\) means an effective horizon of about \(\dfrac{1}{1-\alpha}\) stages.

Example 4.5: The shop with no last day

Return to the shop of Example 1.1 and let it trade for ever with discount factor \(\alpha\) . In Chapter 2 the three-day policy filled the shelf on the first two days and held back on the last. With no last day, is “order up to two” optimal?

Its cost can be found by reasoning. Under this policy the shelf holds \(2\) units after every order, the holding cost is \(1\) in expectation, and there are no shortages. From stock \(x\) the first day costs \((2-x)+1\) . On every later day the order replaces the previous day’s demand, whose mean is \(1\) , so the expected cost is \(1+1=2\) . Hence

\[ J(x)=3-x+\frac{2\alpha}{1-\alpha},\qquad\text{and for }\alpha=0.9:\quad J=(21,\,20,\,19). \]

Now test Bellman’s equation in state \(0\) . Ordering two units costs \(J(0)\) by construction. Ordering one unit costs

\[ \tfrac52+\alpha\big(\tfrac14J(1)+\tfrac34J(0)\big)=\tfrac52+\alpha J(0)-\tfrac{\alpha}{4}, \]

using \(J(1)=J(0)-1\) , and since \((1-\alpha)J(0)=3-\alpha\) this exceeds \(J(0)\) by \(\tfrac34\alpha-\tfrac12\) . Ordering nothing is worse still. The same comparison arises in state \(1\) . So “order up to two” satisfies Bellman’s equation, and is optimal, exactly when \(\alpha\ge\tfrac23\) . For \(\alpha=0.9\) the alternatives in state \(0\) cost \(21.175\) for one unit and \(23.9\) for none, against \(21\) .

For smaller \(\alpha\) the shop behaves as it did on its last day in Chapter 2. A discount factor below \(\tfrac23\) corresponds to an effective horizon shorter than three days, and Example 2.1 found that with one day to go the second unit is not worth buying.

4.6 The edge of the theory

Everything above used \(\alpha\) strictly less than \(1\) , a finite state set, and bounded costs. It is worth seeing what each protects against.

If \(\alpha=1\) the operator \(T\) still makes sense, but Lemma 4.3 now says \(T(J+c\mathbf 1)=TJ+c\mathbf 1\) , so \(T\) does not shrink distances and Bellman’s equation can lose uniqueness. In a model with all costs equal to zero, every constant function solves \(J=TJ\) . In the machine example with \(\alpha=1\) every policy has infinite cost from every state, and the criterion no longer separates good policies from bad ones. Undiscounted problems need different tools: an absorbing cost-free state that every policy eventually reaches, which is the setting of the shortest-path problems of Chapter 3, or a criterion of average cost per stage, hinted at by the normalization in Figure 4.1. Neither is developed in this book.

With infinitely many states the same theorems hold for bounded costs, with suprema in place of maxima; with unbounded costs, as in the linear-quadratic problem, the contraction argument is not available and each problem needs its own analysis, as in Theorem 3.6.

Theorem 4.6 says what the optimal cost is. It does not say how to compute it, since \(J=TJ\) is a nonlinear system. Example 4.2 already contains one method, iterating \(T\) , and Example 4.3 contains the germ of another, improving a policy by testing it. Chapter 5 develops both into algorithms and compares them.

Exercises

  1. Write the definitions of \(T_\mu\) and \(T\) , state which of the two is an affine map of \(J\) , and say what the fixed point of each one is.

    Answer

    \((T_\mu J)(i)=g(i,\mu(i))+\alpha\sum_jp_{ij}(\mu(i))J(j)\) and \((TJ)(i)=\min_{u\in U(i)}[g(i,u)+\alpha\sum_jp_{ij}(u)J(j)]\) . \(T_\mu\) is affine, \(T_\mu J=g_\mu+\alpha P_\mu J\) ; \(T\) is not, because of the minimum. The fixed point of \(T_\mu\) is the cost \(J_\mu\) of the stationary policy \(\mu\) ; the fixed point of \(T\) is the optimal cost \(J^\star\) .

  2. For the machine of Example 4.1 with \(\alpha=0.9\) , compute \(TJ\) for \(J=(5,10)\) , and name the control that attains the minimum in state \(2\) .

    Answer

    \((TJ)(1)=0.9\,(0.7\cdot5+0.3\cdot10)=5.85\) and \((TJ)(2)=\min\{2+0.9\cdot10,\ 5+0.9\cdot5\}=\min\{11,\ 9.5\}=9.5\) , attained by “overhaul”. So \(TJ=(5.85,\ 9.5)\) .

  3. Evaluate both stationary policies of Example 4.1 with \(\alpha=0.5\) by solving the two linear systems, and use the test of Example 4.3 to decide which is optimal.

    Solution

    Never overhaul. \(J(2)=2+\tfrac12J(2)\) gives \(J(2)=4\) . Then \(J(1)=\tfrac12\big(0.7J(1)+0.3\cdot4\big)\) gives \(0.65J(1)=0.6\) , so \(J(1)=\tfrac{12}{13}\) .

    Overhaul when worn. \(J(2)=5+\tfrac12J(1)\) and \(J(1)=\tfrac12\big(0.7J(1)+0.3J(2)\big)\) . Substituting, \(0.65J(1)=0.15\,(5+\tfrac12J(1))=0.75+0.075J(1)\) , so \(J(1)=\tfrac{0.75}{0.575}=\tfrac{30}{23}\) and \(J(2)=5+\tfrac{15}{23}=\tfrac{130}{23}\) .

    Test. For “never overhaul”, in state \(2\) : running costs \(2+\tfrac12\cdot4=4\) and overhauling costs \(5+\tfrac12\cdot\tfrac{12}{13}=\tfrac{71}{13}\approx5.46\) . The minimum is \(4=J(2)\) , attained by running. Bellman’s equation holds in state \(2\) , and it holds in state \(1\) because there is only one control. So \(\big(\tfrac{12}{13},4\big)\) is \(J^\star\) and “never overhaul” is optimal at \(\alpha=\tfrac12\) .

    This agrees with Example 4.4, since \(\tfrac12\) is below the threshold \(\tfrac{30}{41}\) . A decision maker whose effective horizon is two weeks will not pay \(5\) to avoid a weekly cost of \(2\) .

  4. Prove the two bounds of Section 4.1: that \(\lvert J_\pi(i)\rvert\le M/(1-\alpha)\) , and that the \(N\) -stage cost of a policy differs from \(J_\pi(i)\) by at most \(\alpha^{N}M/(1-\alpha)\) .

    Solution

    Let \(c_k=E[g(x_k,\mu_k(x_k))]\) be the expected cost at stage \(k\) . Since \(\lvert g(i,u)\rvert\le M\) for all \(i\) and \(u\) , the random variable \(g(x_k,\mu_k(x_k))\) lies in \([-M,M]\) and so does its expectation: \(\lvert c_k\rvert\le M\) .

    The \(N\) -stage cost is \(\sum_{k=0}^{N-1}\alpha^{k}c_k\) . The series \(\sum_k\alpha^{k}c_k\) converges absolutely by comparison with the geometric series \(\sum_k\alpha^{k}M\) , so the limit defining \(J_\pi(i)\) exists, and

    \[ \lvert J_\pi(i)\rvert\le\sum_{k=0}^{\infty}\alpha^{k}M=\frac{M}{1-\alpha}. \]

    The difference between \(J_\pi(i)\) and the \(N\) -stage cost is the tail \(\sum_{k=N}^{\infty}\alpha^{k}c_k\) , whose absolute value is at most \(\sum_{k=N}^{\infty}\alpha^{k}M=\dfrac{\alpha^{N}M}{1-\alpha}\) . Both bounds hold for every policy and every initial state, with the same constant, which is what makes the convergence uniform.

  5. Let \(P\) be a transition matrix and \(\alpha\in(0,1)\) . Prove that \(I-\alpha P\) is invertible and that its inverse is \(\sum_{k=0}^{\infty}\alpha^{k}P^{k}\) , and interpret the \((i,j)\) entry of the inverse in terms of the Markov chain.

    Solution

    Every entry of \(P^{k}\) lies in \([0,1]\) , because \(P^{k}\) is a transition matrix. So each entry of the partial sums \(S_m=\sum_{k=0}^{m}\alpha^{k}P^{k}\) is a partial sum of a series dominated by \(\sum_k\alpha^{k}\) , and \(S_m\) converges entrywise to a matrix \(S\) . Multiplying out,

    \[ (I-\alpha P)S_m=I-\alpha^{m+1}P^{m+1}\to I\qquad(m\to\infty), \]

    because the entries of \(\alpha^{m+1}P^{m+1}\) are at most \(\alpha^{m+1}\) . Hence \((I-\alpha P)S=I\) , and similarly \(S(I-\alpha P)=I\) ; so \(I-\alpha P\) is invertible with inverse \(S\) .

    The \((i,j)\) entry of \(P^{k}\) is the probability that the chain started at \(i\) is at \(j\) at time \(k\) . So the \((i,j)\) entry of \(S\) is \(\sum_k\alpha^{k}P(x_k=j\mid x_0=i)\) , the expected discounted number of visits to \(j\) starting from \(i\) . The formula \(J_\mu=(I-\alpha P_\mu)^{-1}g_\mu\) then reads: the cost of a policy is the sum over states of the expected discounted number of visits times the cost per visit.

  6. Using the machine of Example 4.1, find two vectors \(J\) and \(J'\) for which \(T(J+J')\ne TJ+TJ'\) , and explain which part of the definition of \(T\) is responsible.

    Hint

    Choose \(J\) and \(J'\) so that different controls are minimizing in state \(2\) . For instance \(J=(0,10)\) makes overhauling attractive and \(J'=(10,0)\) makes running attractive; compute \((TJ)(2)\) , \((TJ')(2)\) , and \((T(J+J'))(2)\) . The minimum of a sum is not the sum of the minima.

  7. A student evaluates the policy “never overhaul” in Example 4.1, obtains \(J=(14.595,20)\) , checks that \(J(2)=2+0.9J(2)\) and \(J(1)=0.9(0.7J(1)+0.3J(2))\) , and concludes that “\(J\) satisfies Bellman’s equation, so the policy is optimal”. State precisely which equation the student verified, which one was required, and carry out the required check.

    Solution

    The student verified \(J=T_\mu J\) for the policy \(\mu\) “never overhaul”: in each state the equation uses only the control that \(\mu\) prescribes. By Theorem 4.5 every stationary policy’s cost satisfies its own such equation, so the check carries no information about optimality.

    What was required is \(J=TJ\) , in which the right-hand side is a minimum over all controls available in the state. In state \(1\) there is one control and nothing more to check. In state \(2\) the two controls give \(2+0.9\cdot20=20\) for running and \(5+0.9\cdot14.595=18.135\) for overhauling. So \((TJ)(2)=18.135\) , which is less than \(J(2)=20\) : Bellman’s equation fails, and the policy is not optimal.

    The failed check is informative. It shows that using “overhaul” once in state \(2\) and then following “never overhaul” already beats the policy, and Chapter 5 turns that observation into the policy improvement step.

  8. Consider a model with two states, one control in each, all costs zero, and \(\alpha=1\) . Show that every constant function is a fixed point of \(T\) , and explain which step of the proof of Theorem 4.5 fails.

    Answer

    With zero costs and \(\alpha=1\) , \((TJ)(i)=\sum_jp_{ij}J(j)\) , and for a constant function \(J=c\mathbf 1\) this is \(c\) . So every constant is a fixed point. The failing step is the uniqueness argument: \(\lVert\bar J-J'\rVert\le\alpha\lVert\bar J-J'\rVert\) forces the distance to be zero only when \(\alpha\) is less than \(1\) . The Cauchy estimate also fails, since \(\sum_k\alpha^{k}\) diverges.

  9. A contract is cancelled at the end of each month with probability \(0.2\) , independently of everything else, and costs \(3\) in every month in which it is active. Compute the expected total cost in two ways: by summing a series directly, and by Proposition 4.7.

    Answer

    The contract is active in month \(k\) with probability \(0.8^{k}\) , so the expected total cost is \(\sum_{k=0}^{\infty}0.8^{k}\cdot3=\dfrac{3}{1-0.8}=15\) . By Proposition 4.7 this is the discounted cost of a constant stream of \(3\) with \(\alpha=0.8\) , again \(\dfrac{3}{1-\alpha}=15\) .

  10. In Example 4.1 let the overhaul cost be \(c\) in place of \(5\) , with \(\alpha=0.9\) . Find all values of \(c\) for which “overhaul when worn” is optimal.

    Solution

    Evaluate “overhaul when worn” with overhaul cost \(c\) . From \(J(1)=0.9\,(0.7J(1)+0.3J(2))\) , \(J(1)=\tfrac{0.27}{0.37}J(2)\) . Then \(J(2)=c+0.9J(1)=c+\tfrac{0.243}{0.37}J(2)\) , so \(J(2)=\tfrac{0.37}{0.127}\,c=\tfrac{370}{127}c\) .

    By Theorem 4.6 this policy is optimal exactly when its cost satisfies Bellman’s equation, and the only state with a choice is state \(2\) . There the policy’s own control gives \(J(2)\) , and the alternative, running, gives \(2+0.9J(2)\) . The condition \(J(2)\le2+0.9J(2)\) is \(J(2)\le20\) , that is,

    \[ \tfrac{370}{127}c\le20,\qquad c\le\tfrac{254}{37}\approx6.865 . \]

    For \(c\) above this value “never overhaul” is optimal, with \(J(2)=20\) ; at the threshold both policies have cost \(20\) from the worn state and both are optimal. With \(c=5\) the condition holds with room to spare, as Example 4.3 found.

  11. For a model with \(M=5\) and \(\alpha=0.9\) , how many stages \(N\) guarantee that the optimal \(N\) -stage cost \(T^{N}0\) is within \(0.01\) of \(J^\star\) in every state? Use Theorem 4.5 with the bound \(\lVert J^\star\rVert\le M/(1-\alpha)\) .

    Answer

    \(\lVert T^{N}0-J^\star\rVert\le\alpha^{N}\lVert J^\star\rVert\le0.9^{N}\cdot50\) . This is at most \(0.01\) when \(0.9^{N}\le2\times10^{-4}\) , that is, \(N\ge\dfrac{\log(2\times10^{-4})}{\log 0.9}\approx80.8\) . So \(N=81\) stages suffice.

  12. For the shop of Example 4.5 with \(\alpha=\tfrac12\) , compute the cost vector of the policy “order up to one” by solving \(J=T_\mu J\) , and verify that it satisfies Bellman’s equation.

    Solution

    Under “order up to one” the controls are \(1,0,0\) in states \(0,1,2\) , with expected stage costs \(\tfrac52,\tfrac32,1\) and the transition matrix of Exercise 1.4. With \(\alpha=\tfrac12\) the equations \(J=T_\mu J\) are

    \[ J(0)=\tfrac52+\tfrac12\big(\tfrac34J(0)+\tfrac14J(1)\big),\quad J(1)=\tfrac32+\tfrac12\big(\tfrac34J(0)+\tfrac14J(1)\big),\quad J(2)=1+\tfrac12\big(\tfrac14J(0)+\tfrac12J(1)+\tfrac14J(2)\big). \]

    The first two differ only in their constants, so \(J(1)=J(0)-1\) . Then \(J(0)=\tfrac52+\tfrac12\big(J(0)-\tfrac14\big)\) , giving \(J(0)=\tfrac{19}{4}\) and \(J(1)=\tfrac{15}{4}\) . The third gives \(\tfrac78J(2)=1+\tfrac{19}{32}+\tfrac{15}{16}=\tfrac{81}{32}\) , so \(J(2)=\tfrac{81}{28}\) . Thus \(J=\big(\tfrac{19}{4},\tfrac{15}{4},\tfrac{81}{28}\big)\) .

    Now Bellman’s equation. In state \(0\) : ordering nothing costs \(5+\tfrac12J(0)=\tfrac{59}{8}\) ; ordering one costs \(J(0)=\tfrac{19}{4}=4.75\) ; ordering two costs \(3+\tfrac12\big(\tfrac14J(2)+\tfrac12J(1)+\tfrac14J(0)\big)=\tfrac{137}{28}\approx4.893\) . The minimum is \(\tfrac{19}{4}\) , the policy’s own value. In state \(1\) : ordering nothing costs \(J(1)=\tfrac{15}{4}=3.75\) and ordering one costs \(\tfrac{137}{28}-1\approx3.893\) . In state \(2\) there is no choice. Bellman’s equation holds, so “order up to one” is optimal at \(\alpha=\tfrac12\) , consistently with the threshold \(\tfrac23\) of Example 4.5.

  13. Prove that if a function \(J\) satisfies \(J\le TJ\) then \(J\le J^\star\) , and that if \(J\ge TJ\) then \(J\ge J^\star\) .

    Solution

    Suppose \(J\le TJ\) . Applying the monotone operator \(T\) to both sides gives \(TJ\le T^{2}J\) , and by induction \(T^{k}J\le T^{k+1}J\) for every \(k\) . So \(J\le TJ\le T^{2}J\le\dots\le T^{k}J\) for every \(k\) . By Theorem 4.5 and Theorem 4.6, \(T^{k}J\to J^\star\) as \(k\to\infty\) , and a limit of vectors each of which is at least \(J\) is at least \(J\) . Hence \(J\le J^\star\) .

    The second statement is proved in the same way with every inequality reversed: \(J\ge TJ\) implies \(T^{k}J\ge T^{k+1}J\) , so \(J\ge T^{k}J\to J^\star\) .

    The result says that \(J^\star\) is the largest function satisfying \(J\le TJ\) and the smallest satisfying \(J\ge TJ\) . The first description is the basis of the linear programming method of Chapter 5; the second gives upper bounds on the optimal cost from any function that the Bellman operator does not increase.

  14. Suppose the stage quantities are rewards \(r(i,u)\) to be maximized. Define the corresponding operator with a maximum, prove that it is a contraction with modulus \(\alpha\) in the maximum norm, and relate its fixed point to the \(J^\star\) of the cost problem with \(g=-r\) .

    Hint

    Define \((T^{r}J)(i)=\max_u[r(i,u)+\alpha\sum_jp_{ij}(u)J(j)]\) . Lemma 4.3 holds for it with the same proof, since a maximum of nondecreasing brackets is nondecreasing and is shifted by a common constant; then repeat the proof of Proposition 4.4. For the relation, show that \(T^{r}(-J)=-TJ\) when \(g=-r\) , so the fixed point of \(T^{r}\) is \(-J^\star\) .

  15. Let \(P=\begin{pmatrix}0&1\\0&1\end{pmatrix}\) , \(g_\mu=0\) , and \(\alpha=0.9\) . Show that \(T_\mu\) is not a contraction in the Euclidean norm by computing \(\lVert T_\mu J-T_\mu J'\rVert_2\) for \(J=(0,1)\) and \(J'=(0,0)\) , and explain why this does not contradict Proposition 4.4.

    Answer

    \(T_\mu J-T_\mu J'=\alpha P(J-J')=0.9\,(1,1)\) , whose Euclidean norm is \(0.9\sqrt2\approx1.273\) , while \(\lVert J-J'\rVert_2=1\) . Distances grew. Proposition 4.4 is a statement about the maximum norm, in which \(\lVert T_\mu J-T_\mu J'\rVert=0.9\) and \(\lVert J-J'\rVert=1\) . The contraction property depends on the norm; the maximum norm fits because each row of \(P\) is a probability distribution.

  16. Suppose the cost of a transition depends also on the next state, as \(g(i,u,j)\) . Show that the theory of this chapter applies with the expected stage cost \(\bar g(i,u)=\sum_jp_{ij}(u)g(i,u,j)\) , and compute \(\bar g(0,2)\) for the shop from the stage cost of Example 1.1.

    Hint

    Write the expected cost of stage \(k\) given \(x_k=i\) and \(u_k=u\) by conditioning on the next state. Only this expectation enters \(J_\pi\) , \(T_\mu\) , and \(T\) . For the shop, the stage cost of Example 1.1 depends on demand, and so on the next state; its expectation for \(x=0\) , \(u=2\) is the entry \(\bar g(0,2)=3\) of the table there.

  17. For the policy “overhaul when worn” of Example 4.1, find the stationary distribution of its Markov chain and the long-run average cost per week, and compare the latter with \((1-\alpha)J_\mu(1)\) at \(\alpha=0.9\) and \(\alpha=0.99\) .

    Solution

    Under “overhaul when worn” the transition matrix is \(P=\begin{pmatrix}0.7&0.3\\1&0\end{pmatrix}\) . A stationary distribution \(\pi=(\pi_1,\pi_2)\) satisfies \(\pi P=\pi\) : \(\pi_2=0.3\pi_1\) , and with \(\pi_1+\pi_2=1\) this gives \(\pi=\big(\tfrac{10}{13},\tfrac{3}{13}\big)\) . In the long run the machine is worn, and an overhaul is paid for, in a fraction \(\tfrac{3}{13}\) of the weeks. The long-run average cost per week is \(0\cdot\tfrac{10}{13}+5\cdot\tfrac{3}{13}=\tfrac{15}{13}\approx1.1538\) .

    From Example 4.4, \(J_\mu(1)=\dfrac{15\alpha}{(1-\alpha)(10+3\alpha)}\) , so \((1-\alpha)J_\mu(1)=\dfrac{15\alpha}{10+3\alpha}\) . At \(\alpha=0.9\) this is \(\tfrac{13.5}{12.7}\approx1.0630\) , and at \(\alpha=0.99\) it is \(\tfrac{14.85}{12.97}\approx1.1449\) . As \(\alpha\to1\) it tends to \(\tfrac{15}{13}\) .

    The normalized discounted cost is a weighted average of the expected stage costs with weights \((1-\alpha)\alpha^{k}\) . As \(\alpha\to1\) the weights spread over an ever longer horizon, the influence of the initial state fades, and the average approaches the stationary average. That the chain forgets its starting state is a property of Markov chains taken from Stochastic Processes (math.stochastic-processes.stationary-distribution).

  18. Suppose a policy \(\mu\) and a function \(J\) satisfy \(T_\mu J=TJ\) and \(\lVert J-J^\star\rVert\le\varepsilon\) , but \(J\) is not exactly \(J^\star\) . Give an example, using Example 4.2, in which \(\mu\) is not optimal, and explain why this does not contradict Theorem 4.6.

    Hint

    Take \(J=J_2=(0.54,\,3.8)\) from Example 4.2. Which control attains the minimum in \((TJ)(2)\) ? Is that the optimal control? Theorem 4.6 characterizes optimal policies by \(T_\mu J^\star=TJ^\star\) , with the exact optimal cost; a policy that is greedy for an approximation of \(J^\star\) need not be optimal, and Chapter 10 bounds how much it can lose.

Search knowledge base

Explore

On this page

Note