Unit 12, lesson 95

Dynamic allocation and stochastic control

How the Bellman principle leads to the Hamilton-Jacobi-Bellman equation for dynamic allocation, how CRRA utility gives the Merton allocation and its value, and what a verification argument proves about admissible controls.

  • Video 95
  • Growth, dynamic allocation and liabilities

An investor with a three-year horizon may change the fraction of wealth held in a stock at every moment, in response to everything observed so far. Which rule for that fraction maximizes the expected utility of final wealth, and how does the answer differ from a weight chosen once, as in mean-variance optimization? Stochastic control treats the rule itself as the unknown and describes the best achievable continuation by a value function of time and wealth. This lesson derives the Hamilton-Jacobi-Bellman equation for that value function, solves it for power utility in a market with constant opportunities to obtain the Merton allocation together with its value, and proves optimality by a verification argument within a stated class of admissible controls.

Key terms

Stochastic control

Definition Fix a horizon T>0T > 0 and a probability space with the physical measure P\mathbb{P}, a Brownian motion WW under P\mathbb{P} and the filtration (Ft)(\mathcal{F}_t) generated by WW and completed with the null sets. A cash account Bt=ertB_t = e^{rt} earns a constant continuously compounded rate rr, and a stock that pays no dividends follows

dSt=μSt dt+σSt dWtdS_t = \mu S_t\,dt + \sigma S_t\,dW_t

with constants μ\mu and σ>0\sigma > 0. An investor keeps the fraction πt\pi_t of current wealth XtX_t in the stock and the rest in the cash account, trades without costs, may borrow or sell short at the rate rr, and neither consumes nor adds money. The self-financing condition then gives the controlled wealth equation

dXt=Xt[r+πt(μ−r)] dt+πtσXt dWt.\begin{aligned} dX_t &= X_t \bigl[r + \pi_t (\mu - r)\bigr]\,dt \\[0.4em] &\quad + \pi_t \sigma X_t\,dW_t . \end{aligned}

A control on [t,T][t, T] is a real-valued process π\pi that is progressively measurable with respect to (Fs)(\mathcal{F}_s), so that each πs\pi_s uses only the information available at time ss, and that satisfies ∫tTπs2 ds<∞\int_t^T \pi_s^2\,ds < \infty almost surely. For such a control and a starting wealth Xt=x>0X_t = x > 0 the wealth equation has the unique solution

Xs=xexp⁡(∫tsπuσ dWu+∫ts(r+πu(μ−r)−12πu2σ2) du),\begin{aligned} X_s &= x \exp\Bigl(\int_t^s \pi_u \sigma\,dW_u \\[0.4em] &\quad + \int_t^s \bigl(r + \pi_u(\mu - r) \\[0.4em] &\qquad - \tfrac12 \pi_u^2 \sigma^2\bigr)\,du\Bigr), \end{aligned}

which is strictly positive. Because the control is a fraction of current wealth, the positive-wealth constraint holds automatically, while a strategy stated in currency amounts would need it as an explicit condition. The admissible controls form a class A(t)\mathcal{A}(t) that is part of the problem. It contains only controls for which the expected utility below is well defined, it may impose constraints such as a cap on π\pi, and the verification card fixes the class used for power utility. The stochastic control problem is to maximize the expected utility of terminal wealth, and its value function is

J(t,x)=sup⁡π∈A(t)Et,xP[U(XTπ)],\begin{aligned} &J(t, x) \\[0.4em] &= \sup_{\pi \in \mathcal{A}(t)} \mathbb{E}^{\mathbb{P}}_{t,x}\bigl[U(X^{\pi}_T)\bigr], \end{aligned}

where UU is an increasing utility function on (0,∞)(0, \infty), Et,xP\mathbb{E}^{\mathbb{P}}_{t,x} is the expectation under P\mathbb{P} for wealth started at Xt=xX_t = x, and J(T,x)=U(x)J(T, x) = U(x). The expectation uses P\mathbb{P}, because the problem is about preferences over actual outcomes, not about the price of a claim.

In plain words The unknown is a rule rather than a number. At every moment the rule looks at everything observed so far and sets the stock fraction, and the value function records the best expected utility that can still be reached from each time and wealth. Buying some shares once and holding them, with the rest in cash, is one admissible rule among many, and dynamic control searches over all of them.

Example In a constructed market with μ=0.08\mu = 0.08, r=0.02r = 0.02 and σ=0.2\sigma = 0.2, the constant fraction πt=0.375\pi_t = 0.375 started from x=1x = 1 at time 00 gives the geometric Brownian motion

Xt=exp⁡(0.0396875 t+0.075 Wt),\begin{aligned} X_t &= \exp\bigl(0.0396875\,t \\[0.4em] &\quad + 0.075\,W_t\bigr), \end{aligned}

since 0.02+0.375⋅0.06−12⋅0.3752⋅0.04=0.03968750.02 + 0.375 \cdot 0.06 - \tfrac12 \cdot 0.375^2 \cdot 0.04 = 0.0396875 and 0.375⋅0.2=0.0750.375 \cdot 0.2 = 0.075. After three years the expected wealth is e0.0425⋅3≈1.135985e^{0.0425 \cdot 3} \approx 1.135985 and the median wealth is e0.0396875⋅3≈1.126440e^{0.0396875 \cdot 3} \approx 1.126440. Holding the fraction fixed is itself a trading rule. Because 0.3750.375 lies between 00 and 11, a rise of the stock relative to the cash account would push the stock’s share of wealth above 0.3750.375 if the investor did not trade, so keeping the fraction means selling after rises and buying after falls.

Bellman principle

Definition The Bellman principle, also called the dynamic programming principle, states that for 0≤t<t+h≤T0 \le t < t + h \le T and x>0x > 0

J(t,x)=sup⁡πEt,xP[J(t+h,Xt+hπ)],\begin{aligned} &J(t, x) \\[0.4em] &= \sup_{\pi} \mathbb{E}^{\mathbb{P}}_{t,x}\bigl[J(t + h, X^{\pi}_{t+h})\bigr], \end{aligned}

with the supremum over admissible controls on [t,t+h][t, t + h]. The problem from (t,x)(t, x) splits into a choice over the first interval followed by the best continuation from the state (t+h,Xt+h)(t + h, X_{t+h}) that is reached, and an optimal control remains optimal, with probability one, for the remaining problem from the state it reaches. The split works because the objective is an expectation of a function of terminal wealth and because the pair (t,Xt)(t, X_t) summarizes everything about the past that matters for the future, here since the coefficients are constant. Under standard conditions for controlled diffusions the principle is a theorem, but its proof needs measurability arguments beyond this lesson. Here it serves to derive the Hamilton-Jacobi-Bellman (HJB) equation, and for the explicit solution of the Merton card the verification argument confirms it directly. A criterion that is not of this form, such as a mean-variance criterion containing Var⁡(XT)\operatorname{Var}(X_T), does not satisfy the principle as stated, and its best strategy chosen at time 00 need not remain best at later dates.

In plain words Whatever happens first, the investor then faces the same kind of problem with less time left and a different wealth, and will handle that problem as well as possible. So the best first step is the one that leads to the best expected continuation value, and the optimal rule chosen at the start is still the rule one would choose at any later date.

Example Take the constructed market with γ=4\gamma = 4 and T=3T = 3 and the value function of the Merton card below, for which J(1,x)=−A(1) x−3/3J(1, x) = -A(1)\,x^{-3}/3 with A(1)=e−0.1875≈0.829029A(1) = e^{-0.1875} \approx 0.829029. If the investor holds only cash in the first year and acts optimally afterward, wealth at year 1 is e0.02e^{0.02} and the expected continuation value is −0.829029 e−0.06/3≈−0.260250-0.829029\, e^{-0.06}/3 \approx -0.260250. With the fraction 0.3750.375 in the first year, E0,1P[X1−3]=e−0.09375≈0.910510\mathbb{E}^{\mathbb{P}}_{0,1}[X_1^{-3}] = e^{-0.09375} \approx 0.910510, and the expected continuation value is

−13⋅0.829029⋅0.910510≈−0.251613,\begin{aligned} &-\tfrac13 \cdot 0.829029 \cdot 0.910510 \\[0.4em] &\quad \approx -0.251613 , \end{aligned}

which is J(0,1)J(0, 1). The fractions 0.250.25 and 0.50.5 in the first year give about −0.252559-0.252559 each. Among constant first-year fractions the best is 0.3750.375, and its continuation value reproduces the value of the whole problem.

Hamilton-Jacobi-Bellman equation

Definition Suppose JJ is continuous on [0,T]×(0,∞)[0, T] \times (0, \infty) and has continuous partial derivatives JtJ_t, JxJ_x and JxxJ_{xx} on [0,T)×(0,∞)[0, T) \times (0, \infty). For a fraction π∈R\pi \in \mathbb{R} the controlled generator is

LπJ=x[r+π(μ−r)]Jx+12π2σ2x2Jxx.\begin{aligned} \mathcal{L}^{\pi} J &= x \bigl[r + \pi(\mu - r)\bigr] J_x \\[0.4em] &\quad + \tfrac12 \pi^2 \sigma^2 x^2 J_{xx} . \end{aligned}

By Itô’s formula, the drift of J(s,Xs)J(s, X_s) under a constant fraction π\pi is Jt+LπJJ_t + \mathcal{L}^{\pi} J evaluated at (s,Xs)(s, X_s), where the second-order term comes from the quadratic variation d⟨X⟩s=π2σ2Xs2 dsd\langle X \rangle_s = \pi^2 \sigma^2 X_s^2\,ds. If the stochastic integral in Itô’s formula has mean zero, the Bellman principle over [t,t+h][t, t + h], divided by hh with h↓0h \downarrow 0, leads formally to the Hamilton-Jacobi-Bellman equation

Jt+sup⁡π∈RLπJ=0,J(T,x)=U(x),\begin{gathered} J_t + \sup_{\pi \in \mathbb{R}} \mathcal{L}^{\pi} J = 0, \\[0.4em] J(T, x) = U(x), \end{gathered}

for 0≤t<T0 \le t < T and x>0x > 0. The derivation is heuristic. It assumes the smoothness of JJ and the Bellman principle, and it exchanges a limit with an expectation. When the value function is not smooth, the equation is interpreted in the weaker viscosity sense, which lies beyond this lesson. When x>0x > 0, σ>0\sigma > 0, Jx>0J_x > 0 and Jxx<0J_{xx} < 0, the expression LπJ\mathcal{L}^{\pi} J is a strictly concave quadratic in π\pi. Setting its derivative x(μ−r)Jx+πσ2x2Jxxx(\mu - r) J_x + \pi \sigma^2 x^2 J_{xx} to zero gives the unique maximizer

π∗=−(μ−r)Jxσ2xJxx,\pi^{*} = -\frac{(\mu - r) J_x}{\sigma^2 x J_{xx}} ,

and substituting it turns the HJB equation into the nonlinear PDE

Jt+rxJx−(μ−r)2Jx22σ2Jxx=0.\begin{aligned} &J_t + r x J_x \\[0.4em] &\quad - \frac{(\mu - r)^2 J_x^2}{2 \sigma^2 J_{xx}} = 0 . \end{aligned}

If Jxx>0J_{xx} > 0, or if Jxx=0J_{xx} = 0 while (μ−r)Jx≠0(\mu - r) J_x \ne 0, the supremum is +∞+\infty and no maximizer exists. If the fraction is restricted to an interval [a,b][a, b], the supremum is taken over that interval, and strict concavity puts the maximizer at min⁡(max⁡(π∗,a),b)\min(\max(\pi^{*}, a), b), the point of the interval nearest to π∗\pi^{*}.

In plain words Along any admissible path the value of the remaining problem has drift at most zero, and along an optimal path its drift is exactly zero. The equation balances the passage of time in JtJ_t, the expected growth of wealth priced by the marginal value JxJ_x, and the variance of wealth penalized by the curvature JxxJ_{xx}. Expected growth rises linearly in the risky fraction while variance rises quadratically, so negative curvature produces an interior best fraction.

Example For the constructed values μ=0.08\mu = 0.08, r=0.02r = 0.02, σ=0.2\sigma = 0.2, γ=4\gamma = 4 and T=3T = 3, the value function of the Merton card has at (t,x)=(0,1)(t, x) = (0, 1) the derivatives Jx≈0.754840J_x \approx 0.754840, Jxx≈−3.019358J_{xx} \approx -3.019358 and Jt≈−0.023589J_t \approx -0.023589. There

LπJ(0,1)≈0.015097+0.045290 π−0.060387 π2,\begin{aligned} &\mathcal{L}^{\pi} J(0, 1) \approx 0.015097 \\[0.4em] &\quad + 0.045290\,\pi - 0.060387\,\pi^2 , \end{aligned}

a concave quadratic with maximizer 0.045290/(2⋅0.060387)=0.3750.045290/(2 \cdot 0.060387) = 0.375 and maximum about 0.0235890.023589, which cancels JtJ_t as the equation requires. If the fraction were capped at 0.250.25, the maximum over the allowed interval would sit at the endpoint 0.250.25, with value about 0.0226450.022645.

CRRA utility

Definition For a relative risk aversion γ>0\gamma > 0, the CRRA (constant relative risk aversion) or power utility on (0,∞)(0, \infty) is

U(x)=x1−γ1−γ,γ≠1,U(x)=ln⁡x,γ=1.\begin{gathered} U(x) = \frac{x^{1 - \gamma}}{1 - \gamma}, \quad \gamma \ne 1, \\[0.4em] U(x) = \ln x, \quad \gamma = 1 . \end{gathered}

Its derivatives are U′(x)=x−γ>0U'(x) = x^{-\gamma} > 0 and U′′(x)=−γx−γ−1<0U''(x) = -\gamma x^{-\gamma - 1} < 0, so UU is strictly increasing and strictly concave, and the coefficient of relative risk aversion is

−x U′′(x)U′(x)=γ-\frac{x\,U''(x)}{U'(x)} = \gamma

at every wealth. For γ≠1\gamma \ne 1 the utility is homogeneous of degree 1−γ1 - \gamma, with U(cx)=c1−γU(x)U(cx) = c^{1 - \gamma} U(x) for c>0c > 0, and for log utility ln⁡(cx)=ln⁡c+ln⁡x\ln(cx) = \ln c + \ln x. For γ>1\gamma > 1 the utility is negative, and larger wealth still gives a larger, less negative value. The log case is the limit of the shifted utility (x1−γ−1)/(1−γ)(x^{1 - \gamma} - 1)/(1 - \gamma) as γ→1\gamma \to 1, and the shift is needed because the constant 1/(1−γ)1/(1 - \gamma) that separates it from UU diverges. Adding a constant to a utility or multiplying it by a positive constant changes no decision.

In plain words Risk aversion measured relative to wealth is the same at every wealth level, so a gamble that scales with wealth is judged the same way by a rich investor and a poor one. This is why the best fraction of wealth in the stock will turn out not to depend on wealth. A larger γ\gamma means a stronger dislike of proportional risk.

Example With γ=4\gamma = 4 the utility is U(x)=−x−3/3U(x) = -x^{-3}/3, so U(1)=−1/3U(1) = -1/3 and U(2)=−1/24≈−0.041667U(2) = -1/24 \approx -0.041667. A gamble that leaves wealth 0.50.5 or 1.51.5 with probability one half each has expected utility 12⋅(−8/3)+12⋅(−0.098765)≈−1.382716\tfrac12 \cdot (-8/3) + \tfrac12 \cdot (-0.098765) \approx -1.382716. Its certainty equivalent, the sure wealth with the same utility, is

(3⋅1.382716)−1/3≈0.6224,(3 \cdot 1.382716)^{-1/3} \approx 0.6224 ,

far below the mean wealth 11. Doubling both outcomes doubles the certainty equivalent to about 1.24471.2447, as homogeneity requires. For log utility the certainty equivalent of the original gamble is 0.75≈0.8660\sqrt{0.75} \approx 0.8660, and for γ=2\gamma = 2 it is 0.750.75, so it falls as γ\gamma grows.

Merton allocation

Definition In the market of the stochastic control card, with σ>0\sigma > 0 and CRRA utility with γ>0\gamma > 0 and γ≠1\gamma \ne 1, look for a solution J(t,x)=A(t) U(x)J(t, x) = A(t)\,U(x) of the HJB equation with A(T)=1A(T) = 1 and A>0A > 0. Then Jx=Ax−γ>0J_x = A x^{-\gamma} > 0 and Jxx=−γAx−γ−1<0J_{xx} = -\gamma A x^{-\gamma - 1} < 0, and

LπJ=A(t) x1−γ g(π),g(π)=r+π(μ−r)−12γπ2σ2.\begin{aligned} &\mathcal{L}^{\pi} J = A(t)\, x^{1 - \gamma}\, g(\pi), \\[0.4em] &g(\pi) = r + \pi(\mu - r) \\[0.4em] &\qquad - \tfrac12 \gamma \pi^2 \sigma^2 . \end{aligned}

The factor A(t) x1−γA(t)\, x^{1 - \gamma} is positive, so the supremum over π\pi is attained where gg is largest, at the Merton allocation

π∗=μ−rγσ2,\pi^{*} = \frac{\mu - r}{\gamma \sigma^2} ,

which agrees with the general formula of the HJB card, since −Jx/(xJxx)=1/γ-J_x/(x J_{xx}) = 1/\gamma. The largest value of gg is

b=g(π∗)=r+(μ−r)22γσ2.b = g(\pi^{*}) = r + \frac{(\mu - r)^2}{2 \gamma \sigma^2} .

The HJB equation becomes A′(t) x1−γ/(1−γ)+A(t) b x1−γ=0A'(t)\, x^{1 - \gamma}/(1 - \gamma) + A(t)\, b\, x^{1 - \gamma} = 0. Dividing by x1−γ>0x^{1 - \gamma} > 0 and multiplying by the nonzero number 1−γ1 - \gamma, which is negative when γ>1\gamma > 1, gives A′+(1−γ) b A=0A' + (1 - \gamma)\, b\, A = 0 with A(T)=1A(T) = 1. Hence A(t)=e(1−γ)b(T−t)>0A(t) = e^{(1 - \gamma) b (T - t)} > 0 and

J(t,x)=x1−γ1−γ×e(1−γ) b (T−t).\begin{aligned} J(t, x) &= \frac{x^{1 - \gamma}}{1 - \gamma} \\[0.4em] &\quad \times e^{(1 - \gamma)\, b\, (T - t)} . \end{aligned}

For γ=1\gamma = 1 the same steps with J(t,x)=ln⁡x+a(t)J(t, x) = \ln x + a(t) give π∗=(μ−r)/σ2\pi^{*} = (\mu - r)/\sigma^2, the fraction that maximizes the expected log growth rate and the continuous-time counterpart of the fraction in Expected utility and Kelly growth, and J(t,x)=ln⁡x+b (T−t)J(t, x) = \ln x + b\,(T - t) with b=r+(μ−r)2/(2σ2)b = r + (\mu - r)^2/(2\sigma^2). In every case the certainty equivalent of J(t,x)J(t, x), the sure terminal wealth with that utility, is x eb(T−t)x\, e^{b (T - t)}, so bb is the certainty-equivalent growth rate. The verification card proves that π∗\pi^{*} is optimal within a stated class of admissible controls. The fraction is constant because rr, μ\mu and σ\sigma are constant and relative risk aversion is constant. The fraction exceeds 11 when μ−r>γσ2\mu - r > \gamma \sigma^2, and the investor then borrows. It is negative, a short position, when μ<r\mu < r, and it is 00 when μ=r\mu = r. With random investment opportunities the value function depends on further state variables and the optimal fraction acquires hedging terms, and consumption or transaction costs change the problem. These extensions are beyond this lesson.

In plain words Keep a fixed proportion of wealth in the stock, equal to the excess return divided by risk aversion times variance, and trade continuously to hold it there. The formula has the same form as the one-period weight that maximizes excess return minus γ/2\gamma/2 times variance, but here it is the best rule among all rules that react to information, and the value function also says what that rule is worth.

Example For the constructed values μ=0.08\mu = 0.08, r=0.02r = 0.02, σ=0.2\sigma = 0.2, γ=4\gamma = 4, T=3T = 3 and x=1x = 1, the excess return is 0.060.06 and γσ2=0.16\gamma \sigma^2 = 0.16, so π∗=0.375\pi^{*} = 0.375. Then b=0.02+0.0036/0.32=0.03125b = 0.02 + 0.0036/0.32 = 0.03125 and (1−γ) b T=−0.28125(1 - \gamma)\, b\, T = -0.28125, so

J(0,1)=−13 e−0.28125≈−13⋅0.754840≈−0.251613.\begin{aligned} J(0, 1) &= -\tfrac13\, e^{-0.28125} \\[0.4em] &\approx -\tfrac13 \cdot 0.754840 \\[0.4em] &\approx -0.251613 . \end{aligned}

The value is negative because power utility with γ>1\gamma > 1 is negative, not because wealth is lost. Its certainty equivalent is e0.09375≈1.098285e^{0.09375} \approx 1.098285, compared with e0.06≈1.061837e^{0.06} \approx 1.061837 from the cash account alone. The investor holds 37.5 percent of wealth in the stock and 62.5 percent in cash at all times. With log utility the same market gives π∗=1.5\pi^{*} = 1.5, a leveraged position that borrows half of wealth.

Verification argument

Definition Let vv be a function on [0,T]×(0,∞)[0, T] \times (0, \infty) with continuous partial derivatives vtv_t, vxv_x and vxxv_{xx} that satisfies the HJB equation on [0,T)×(0,∞)[0, T) \times (0, \infty) with v(T,x)=U(x)v(T, x) = U(x). For an admissible control π\pi and Xt=xX_t = x, Itô’s formula gives

U(XT)=v(t,x)+MT+∫tT(vt+Lπsv)(s,Xs) ds,\begin{aligned} U(X_T) &= v(t, x) + M_T \\[0.4em] &{}+ \int_t^T \bigl(v_t + \mathcal{L}^{\pi_s} v\bigr)(s, X_s)\,ds , \end{aligned}

where

Ms=∫tsvx(u,Xu)×πuσXu dWu.\begin{aligned} M_s &= \int_t^s v_x(u, X_u) \\[0.4em] &\quad \times \pi_u \sigma X_u\,dW_u . \end{aligned}

The HJB equation makes the integrand of the dsds integral at most zero for every value of πs\pi_s. First, if E[MT]=0\mathbb{E}[M_T] = 0, for example because MM is a square-integrable martingale, then Et,xP[U(XT)]≤v(t,x)\mathbb{E}^{\mathbb{P}}_{t,x}[U(X_T)] \le v(t, x). Second, if a control π^\hat\pi attains the supremum in the HJB equation at (s,Xs)(s, X_s) for almost every ss, its dsds integral vanishes, and if also E[MT]=0\mathbb{E}[M_T] = 0 for it, then Et,xP[U(XTπ^)]=v(t,x)\mathbb{E}^{\mathbb{P}}_{t,x}[U(X^{\hat\pi}_T)] = v(t, x). On a class of admissible controls for which the first step applies and that contains π^\hat\pi, this proves v=Jv = J and the optimality of π^\hat\pi. For a general control MM is only a local martingale. One then stops it at times τn\tau_n increasing to TT, obtains E[v(τn,Xτn)]≤v(t,x)\mathbb{E}[v(\tau_n, X_{\tau_n})] \le v(t, x), and needs uniform integrability of the stopped values to pass to the limit, so these integrability conditions are part of the optimality statement. For the Merton candidate the integrand of MM is A(s) σπsXs1−γA(s)\, \sigma \pi_s X_s^{1 - \gamma}, and the lesson takes as admissible the controls with

Et,xP[∫tTπs2Xs2(1−γ) ds]<∞,\begin{aligned} &\mathbb{E}^{\mathbb{P}}_{t,x}\Bigl[\int_t^T \pi_s^2 X_s^{2(1 - \gamma)}\,ds\Bigr] \\[0.4em] &\quad < \infty , \end{aligned}

for which MM is a square-integrable martingale. This class contains every bounded progressively measurable fraction, because such fractions give wealth with finite moments of every order. Within it the Merton fraction is optimal and the candidate is the value function. For 0<γ<10 < \gamma < 1 the candidate is nonnegative, so MM is a local martingale bounded below by −v(t,x)-v(t, x), hence a supermartingale, and the upper bound extends to every control with ∫tTπs2 ds<∞\int_t^T \pi_s^2\,ds < \infty.

In plain words A smooth solution of the HJB equation is only a candidate. The verification argument compares it with every admissible strategy, shows that none does better on average, and shows that the candidate’s own control reaches the candidate value. The integrability conditions are what allow the expectation of the stochastic integral to be set to zero, so they belong to the result.

Example For a constant fraction π\pi and τ=T−t\tau = T - t, the logarithm ln⁡XT\ln X_T is normal with mean ln⁡x+(r+π(μ−r)−12π2σ2)τ\ln x + \bigl(r + \pi(\mu - r) - \tfrac12 \pi^2 \sigma^2\bigr)\tau and variance π2σ2τ\pi^2 \sigma^2 \tau, and the normal moment generating function gives

Et,xP[XT1−γ]=x1−γ e(1−γ) g(π) τ\begin{aligned} &\mathbb{E}^{\mathbb{P}}_{t,x}\bigl[X_T^{1 - \gamma}\bigr] \\[0.4em] &\quad = x^{1 - \gamma}\, e^{(1 - \gamma)\, g(\pi)\, \tau} \end{aligned}

with the function gg of the Merton card. The expected utility x1−γe(1−γ)g(π)τ/(1−γ)x^{1 - \gamma} e^{(1 - \gamma) g(\pi) \tau}/(1 - \gamma) increases with g(π)g(\pi) for every γ≠1\gamma \ne 1, so constant fractions are ranked by gg, their certainty-equivalent growth rate. For the constructed values with γ=4\gamma = 4, T=3T = 3 and x=1x = 1 the expected utility is −e−9g(π)/3-e^{-9 g(\pi)}/3.

Fraction π\pig(π)g(\pi)Expected utility
000.020.02−0.278423-0.278423
0.250.250.030.03−0.254460-0.254460
0.3750.3750.031250.03125−0.251613-0.251613
0.50.50.030.03−0.254460-0.254460
1100−0.333333-0.333333

The fraction 0.3750.375 reproduces J(0,1)J(0, 1) exactly, which confirms the value as well as the control. Full investment has g(1)=0g(1) = 0 and certainty equivalent 11, the initial wealth, although its expected wealth is e0.24≈1.271249e^{0.24} \approx 1.271249. With a cap πt≤0.25\pi_t \le 0.25 the same ansatz solves the constrained HJB equation with g(0.25)=0.03g(0.25) = 0.03 in place of bb, and the verification argument applies to the capped admissible controls. The capped optimum is therefore the constant 0.250.25, with value −e−0.27/3≈−0.254460-e^{-0.27}/3 \approx -0.254460. This is below −0.251613-0.251613, as it must be, since the cap removes strategies. The less negative number is the larger one.

Check yourself

  1. In a constructed market with μ=0.07\mu = 0.07, r=0.03r = 0.03 and σ=0.25\sigma = 0.25, an investor has CRRA utility with γ=2\gamma = 2 and horizon T=5T = 5. Find the Merton fraction, the values J(0,1)J(0, 1) and J(0,2)J(0, 2), and the certainty equivalent of J(0,1)J(0, 1). Does the fraction change when initial wealth doubles?
Show the answer

The Merton fraction is π∗=0.04/(2⋅0.0625)=0.32\pi^{*} = 0.04/(2 \cdot 0.0625) = 0.32. The certainty-equivalent growth rate is

b=0.03+0.00160.25=0.0364.\begin{aligned} b &= 0.03 + \frac{0.0016}{0.25} \\[0.4em] &= 0.0364 . \end{aligned}

With γ=2\gamma = 2 the utility is U(x)=−1/xU(x) = -1/x, so J(0,x)=−e−0.0364⋅5/x=−e−0.182/xJ(0, x) = -e^{-0.0364 \cdot 5}/x = -e^{-0.182}/x. This gives J(0,1)≈−0.833601J(0, 1) \approx -0.833601 and J(0,2)≈−0.416801J(0, 2) \approx -0.416801, half the size, because JJ is homogeneous of degree 1−γ=−11 - \gamma = -1 in wealth. The certainty equivalent of J(0,1)J(0, 1) is e0.182≈1.199614e^{0.182} \approx 1.199614, above the value e0.15≈1.161834e^{0.15} \approx 1.161834 of the cash account alone. The fraction stays 0.320.32 at every wealth and every date, since the formula does not contain xx or tt.

  1. A classmate wants a risk-seeking investor and takes γ=−1\gamma = -1 in the power utility, so that U(x)=x2/2U(x) = x^2/2. With μ=0.08\mu = 0.08, r=0.02r = 0.02, σ=0.2\sigma = 0.2 and T=3T = 3, the formula (μ−r)/(γσ2)(\mu - r)/(\gamma \sigma^2) gives −1.5-1.5. Is this fraction optimal?
Show the answer

No. The formula is the stationary point of the Hamiltonian, and it is a maximizer only when Jxx<0J_{xx} < 0, which the Merton card obtains from γ>0\gamma > 0. With γ=−1\gamma = -1 the same ansatz gives J=Ax2/2J = A x^2/2 with Jxx=A>0J_{xx} = A > 0, and LπJ=Ax2g(π)\mathcal{L}^{\pi} J = A x^2 g(\pi) with

g(π)=0.02+0.06 π+0.02 π2,g(\pi) = 0.02 + 0.06\,\pi + 0.02\,\pi^2 ,

a convex quadratic. So −1.5-1.5 minimizes the Hamiltonian. Among constant fractions the expected utility is 12e6g(π)\tfrac12 e^{6 g(\pi)}, which equals about 0.4303540.430354 at π=−1.5\pi = -1.5, the smallest value, 0.5637480.563748 at π=0\pi = 0 and 68.50130768.501307 at π=5\pi = 5, and it grows without bound as ∣π∣\lvert \pi \rvert grows. The supremum is +∞+\infty, so there is no optimal control and the value function is infinite. The hypothesis that fails is γ>0\gamma > 0, which makes the utility concave and the curvature JxxJ_{xx} negative.

  1. Under log utility with μ=0.08\mu = 0.08, r=0.02r = 0.02, σ=0.2\sigma = 0.2 and T=3T = 3, the investor may neither borrow nor sell short, so πt\pi_t must lie in [0,1][0, 1]. Find the optimal fraction, the value J(0,1)J(0, 1) and its certainty equivalent, and compare them with the unconstrained problem.
Show the answer

With J(t,x)=ln⁡x+a(t)J(t, x) = \ln x + a(t) the derivatives are Jx=1/xJ_x = 1/x and Jxx=−1/x2J_{xx} = -1/x^2, so the Hamiltonian does not depend on xx and equals

g(π)=0.02+0.06 π−0.02 π2.g(\pi) = 0.02 + 0.06\,\pi - 0.02\,\pi^2 .

This concave quadratic has its unconstrained maximum at π=1.5\pi = 1.5, outside [0,1][0, 1], so over the allowed interval the maximum is at the nearest point π=1\pi = 1, where g(1)=0.06g(1) = 0.06. Solving a′+0.06=0a' + 0.06 = 0 with a(3)=0a(3) = 0 gives J(0,x)=ln⁡x+0.18J(0, x) = \ln x + 0.18, so J(0,1)=0.18J(0, 1) = 0.18 with certainty equivalent e0.18≈1.197217e^{0.18} \approx 1.197217. Verification is direct here. For every control with values in [0,1][0, 1] the stochastic integral in ln⁡XT\ln X_T has mean zero, so E0,1P[ln⁡X3]\mathbb{E}^{\mathbb{P}}_{0,1}[\ln X_3] is the expected integral of g(πs)g(\pi_s) over three years, at most 3⋅0.06=0.183 \cdot 0.06 = 0.18, with equality for πt=1\pi_t = 1. Without the constraint the fraction is 1.51.5, the value is 3⋅0.065=0.1953 \cdot 0.065 = 0.195 and the certainty equivalent is e0.195≈1.215311e^{0.195} \approx 1.215311. The constraint lowers the value, as it must when the set of strategies shrinks, and the constrained optimum holds all wealth in the stock, so it needs no rebalancing.

How this lesson connects