Numerical Methods in Financial and Actuarial Applications: A Stochastic Maximum Principle Approach

Marina Di Giacinto^{1,2}

Show more

1. Introduction

A portfolio simulation approach to the valuation of optimal portfolio policies appears to be the only viable approach in many financial and actuarial applications.

While usually numerical methods are formulated within a dynamic programming approach to the optimization problem, we explore the possibility to state a numerical scheme that goes through the solution of the forward-backward stochastic differential equation (FBSDE, hereafter) which arises from the maximum principle formulation of the optimal control problem. We exemplify situations where the maximum principle approach to portfolio optimization jointly with the simulation approach to the backward stochastic differential equation (BSDE, hereafter) provides an efficient alternative to standard simulation methods based on direct Montecarlo simulation and produces a reliable estimate of the optimal allocation strategies.

Our computation approach relies on the least-squares Monte Carlo approach for estimating the conditional expectations made popular in financial mathematics by of Longstaff and Schwartz [1] and it was first applied to BSDEs and analyzed in this setting by Bouchard and Touzi [2] , Gobet et al. [3] and Lemor et al. [4] .

More recently, Bender and Steiner [5] focus on the least-squares Monte Carlo scheme for approximating the solution of a BSDE by using basis functions which form a system of martingales.

The paper is organized as follows. Section 2 introduces the stochastic maximum principle and the economic implications of the first-order adjoint process. Section 3 illustrates some known results on the FBSDEs. In Section 4 we analyze the simulation approach to the BSDEs. Section 5 is devoted to the discussion of some financial and actuarial applications.

2. The Stochastic Maximum Principle

We consider a finite time horizon $\left[\mathrm{0,}T\right]$ , and a complete probability space $\left(\Omega \mathrm{,}\mathcal{F}\mathrm{,}\mathbb{P}\right)$ with a filtration ${\left\{{\mathcal{F}}_{t}^{B}\right\}}_{t\ge 0}$ , where $t\ge 0$ is the time variable. The filtration ${\left\{{\mathcal{F}}_{t}^{B}\right\}}_{t\ge 0}$ is the completion of the natural filtration generated by a m-dimensional standard Brownian motion $B\left(t\right)$ , $t\ge 0$ . This signifies that the system noise represented by the Brownian motion is the only source of uncertainty in the problem, and the past information about the noise is available to the decision-maker. Sometimes we will use a starting point $s>t$ . In this case ${\left\{{\mathcal{F}}_{t}^{B}\right\}}_{t\ge s}$ will be the complete filtration generated by ${B}^{s}\left(t\right)=B\left(t\right)-B\left(s\right)$ .

Let us consider the time evolution of the system described by a state variable $x(\cdot )$ whose dynamics follow the stochastic differential equation:

$\{\begin{array}{l}\text{d}x\left(t\right)=b\left(t,x\left(t\right),\nu \left(t\right)\right)+\sigma \left(t,x\left(t\right),\nu \left(t\right)\right)\text{d}B\left(t\right),\text{\hspace{1em}}t\in \left[0,T\right],\\ x\left(0\right)={x}_{0},\end{array}$ (1)

where $b:\left[\mathrm{0,}T\right]\times {\mathbb{R}}^{n}\times \mathbb{U}\to {\mathbb{R}}^{n}$ , $\sigma :\left[\mathrm{0,}T\right]\times {\mathbb{R}}^{n}\times \mathbb{U}\to {\mathbb{R}}^{d\times m}$ , and $\mathbb{U}$ is a convex body of ${\mathbb{R}}^{k}$ .

The function $\nu (\cdot )$ is the control variable representing the policy of the decision-maker at any time and it is required to be progressively measurable since the controller knows what happened up to the moment, but not able to foretell what is going to happen later on, due to the uncertainty of the system.

The objective is to minimize a given cost functional defined as:

$J\left(\nu (\cdot )\right)=\mathbb{E}\left\{{\displaystyle {\int}_{0}^{T}}f\left(t,x\left(t\right),\nu \left(t\right)\right)dt+h\left(x\left(T\right)\right)\right\}$ (2)

with $f:\left[\mathrm{0,}T\right]\times {\mathbb{R}}^{n}\times \mathbb{U}\to \mathbb{R}$ and $h:{\mathbb{R}}^{n}\to \mathbb{R}$ , over the set of the admissible controls:

${\mathcal{U}}_{ad}\left({x}_{0}\right)\mathrm{:}=\left\{\nu :\left[\mathrm{0,}T\right]\times \Omega \to \mathbb{U}\subseteq {\mathbb{R}}^{k}\mathrm{|}\nu (\cdot )\text{\hspace{0.17em}}\text{is}\text{\hspace{0.17em}}{\left\{{\mathcal{F}}_{t}^{B}\right\}}_{t\in \left[\mathrm{0,}T\right]}\text{-progr}\text{.}\text{\hspace{0.17em}}\text{meas}\text{.}\right\}\mathrm{.}$

Of course, for any $\nu (\cdot )\in {\mathcal{U}}_{ad}\left({x}_{0}\right)$ , by standard assumptions on b, σ, ν that can be found, e.g., in Yong and Zhou ( [6] , Theorem 6.17, p. 50), the above state Equation (1) admits a unique solution $x(\cdot ):=x\left(\cdot \mathrm{;}\nu (\cdot )\right)$ (in the sense of Yong and Zhou [6] , Definition 6.15, pp. 48-49) and $\left(x(\cdot )\mathrm{,}\nu (\cdot )\right)$ is called an admissible pair.

Then, the problem is to find an optimal control ${\nu}^{\star}(\cdot )\in {\mathcal{U}}_{ad}\left({x}_{0}\right)$ such that:

$J\left({\nu}^{\star}(\cdot )\right)=\underset{\nu (\cdot )\in {\mathcal{U}}_{ad}\left({x}_{0}\right)}{\mathrm{inf}}J\left(\nu (\cdot )\right)\mathrm{.}$

The corresponding ${x}^{\star}(\cdot )\mathrm{:}=x\left(\cdot \mathrm{;}{\nu}^{\star}(\cdot )\right)$ and $\left({x}^{\star}(\cdot )\mathrm{,}{\nu}^{\star}(\cdot )\right)$ are called optimal state trajectory and optimal pair, respectively.

Recalling that, for any given function φ,

$\mathrm{min}\phi =-\mathrm{max}\left\{-\phi \right\},$

we may also treat maximization problems with the same steps.

The Hamiltonian function associated to our problem is given by:

$H\left(t\mathrm{,}x\mathrm{,}\nu \mathrm{,}y\mathrm{,}z\right)\mathrm{:}=\langle y\mathrm{,}b\left(t\mathrm{,}x\mathrm{,}\nu \right)\rangle +\text{tr}\left({z}^{{\rm T}}\sigma \left(t\mathrm{,}x\mathrm{,}\nu \right)\right)-f\left(t\mathrm{,}x\mathrm{,}\nu \right)\mathrm{,}$

$\left(t\mathrm{,}x\mathrm{,}\nu \mathrm{,}y\mathrm{,}z\right)\in \left[\mathrm{0,}T\right]\times {\mathbb{R}}^{n}\times \mathbb{U}\times {\mathbb{R}}^{n}\times {\mathbb{R}}^{d\times m}\mathrm{,}$

while the generalized Hamiltonian function is given by the following FBSDE:

$G\left(t\mathrm{,}x\mathrm{,}\nu \mathrm{,}y\mathrm{,}P\right)\mathrm{:}=\langle y\mathrm{,}b\left(t\mathrm{,}x\mathrm{,}\nu \right)\rangle +\frac{1}{2}\text{tr}\left(\sigma {\left(t\mathrm{,}x\mathrm{,}\nu \right)}^{{\rm T}}P\sigma \left(t\mathrm{,}x\mathrm{,}\nu \right)\right)-f\left(t\mathrm{,}x\mathrm{,}\nu \right)\mathrm{,}$

$\left(t\mathrm{,}x\mathrm{,}\nu \mathrm{,}y\mathrm{,}P\right)\in \left[\mathrm{0,}T\right]\times {\mathbb{R}}^{n}\times \mathbb{U}\times {\mathbb{R}}^{n}\times {\mathcal{S}}^{d}\mathrm{.}$

The extended Hamiltonian system corresponding to the pair $\left({x}^{\star}(\cdot )\mathrm{,}{\nu}^{\star}(\cdot )\right)$ is given by the following FBSDE:

$\{\begin{array}{l}\text{d}x\left(t\right)={H}_{y}\left(t,x\left(t\right),\nu \left(t\right),y\left(t\right),z\left(t\right)\right)\text{d}t\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+{H}_{z}\left(t,x\left(t\right),\nu \left(t\right),y\left(t\right),z\left(t\right)\right)\text{d}B\left(t\right),\text{\hspace{0.17em}}t\in \left[0,T\right]\\ \text{d}y\left(t\right)=-{H}_{x}\left(t,x\left(t\right),\nu \left(t\right),y\left(t\right),z\left(t\right)\right)\text{d}t+z\left(t\right)\text{d}B\left(t\right),\text{\hspace{1em}}t\in \left[0,T\right]\\ \text{d}P\left(t\right)=\{-{b}_{x}{\left(t\mathrm{,}{x}^{\star}\left(t\right)\mathrm{,}{\nu}^{\star}\left(t\right)\right)}^{{\rm T}}P\left(t\right)-P\left(t\right){b}_{x}\left(t\mathrm{,}{x}^{\star}\left(t\right)\mathrm{,}{\nu}^{\star}\left(t\right)\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-{\displaystyle \underset{j=1}{\overset{m}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\sigma}_{x}^{j}{\left(t\mathrm{,}{x}^{\star}\left(t\right)\mathrm{,}{\nu}^{\star}\left(t\right)\right)}^{{\rm T}}P\left(t\right){\sigma}_{x}^{j}\left(t\mathrm{,}{x}^{\star}\left(t\right)\mathrm{,}{\nu}^{\star}\left(t\right)\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-{\displaystyle \underset{j=1}{\overset{m}{\sum}}}\left[{\sigma}^{j}{\left(t\mathrm{,}{x}^{\star}\left(t\right)\mathrm{,}{\nu}^{\star}\left(t\right)\right)}^{{\rm T}}{Q}_{j}\mathrm{(}t\mathrm{)}+{Q}_{j}\left(t\right){\sigma}^{j}\left(t\mathrm{,}{x}^{\star}\left(t\right)\mathrm{,}{\nu}^{\star}\left(t\right)\right)\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-{H}_{xx}\left(t\mathrm{,}{x}^{\star}\left(t\right)\mathrm{,}{\nu}^{\star}\left(t\right)\mathrm{,}y\left(t\right)\mathrm{,}z\left(t\right)\right)\}\text{d}t+{\displaystyle \underset{j=1}{\overset{m}{\sum}}}\text{\hspace{0.05em}}{Q}^{j}\left(t\right){\text{d}}^{j}B\left(t\right)\mathrm{,}\text{\hspace{1em}}t\in \left[\mathrm{0,}T\right]\\ x\left(0\right)={x}_{0}\mathrm{,}\text{\hspace{1em}}y\left(T\right)=-{h}_{x}\left({x}^{\star}\left(T\right)\right)\mathrm{,}\text{\hspace{1em}}P\left(T\right)=-{h}_{xx}\left({x}^{\star}\left(T\right)\right)\\ \mathcal{H}\left(t\mathrm{,}{x}^{\star}\left(t\right)\mathrm{,}{\nu}^{\star}\left(t\right)\right)=\underset{\nu \in \mathbb{U}}{max}\text{\hspace{0.05em}}\mathcal{H}\left(t\mathrm{,}{x}^{\star}\left(t\right)\mathrm{,}\nu \right)\text{\hspace{1em}}\text{a}\text{.e}\text{.}\text{\hspace{0.17em}}t\in \left[\mathrm{0,}T\right]\mathrm{,}\text{\hspace{1em}}\mathbb{P}\text{-a}\text{.s}\mathrm{.,}\end{array}$

where:

$\begin{array}{l}\mathcal{H}\left(t\mathrm{,}x\mathrm{,}\nu \right)\mathrm{:}=G\left(t\mathrm{,}x\mathrm{,}\nu \mathrm{,}y\left(t\right)\mathrm{,}P\left(t\right)\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}+\text{tr}\left\{\sigma {\left(t\mathrm{,}x\mathrm{,}\nu \right)}^{{\rm T}}\left[z\left(t\right)-P\left(t\right)\sigma \left(t\mathrm{,}{x}^{\star}\left(t\right)\mathrm{,}{\nu}^{\star}\left(t\right)\right)\right]\right\}\mathrm{,}\end{array}$

$\forall \left(t\mathrm{,}x\mathrm{,}\nu \right)\in \left[\mathrm{0,}T\right]\times {\mathbb{R}}^{n}\times \mathbb{U}\mathrm{,}$

is the H-function.

In the extended Hamiltonian system above, the second and the third relations are clearly BSDEs and called adjoint equation of the first-order and adjoint equation of the second-order respectively, while the last relation is the so-called maximum condition.

The state variable joint with the related first-order adjoint equation are called Hamiltonian system.

Definition 2.1. $\left({x}^{\star}(\cdot )\mathrm{,}{\nu}^{\star}(\cdot )\mathrm{,}y(\cdot )\mathrm{,}z(\cdot )\mathrm{,}P(\cdot )\mathrm{,}Q(\cdot )\right)$ is called an optimal 6-tuple (resp. admissible) if $\left({x}^{\star}(\cdot )\mathrm{,}{\nu}^{\star}(\cdot )\right)$ is an optimal pair, (resp. is an admissible pair), $\left(y(\cdot )\mathrm{,}z(\cdot )\right)$ is an adapted solution of the first-order adjoint equation and $\left(P(\cdot )\mathrm{,}Q(\cdot )\right)$ is an adapted solution of the second-order adjoint equation.

We assume the following hypothesis for $\phi =b\mathrm{,}\sigma \mathrm{,}f\mathrm{,}h$ .

Hypothesis 2.2. The map φ is twice continuously differentiable in $x\in {\mathbb{R}}^{n}$ and there exist a constant $K>0$ and a modulus of continuity $m:\left[\mathrm{0,}+\infty \right)\to \left[\mathrm{0,}+\infty \right)$ such that, for all,

and, , , ,

The following statement is the well-known (Pontryagin) stochastic maximum principle, which gives a set of first-order necessary condition for 6-tuple pairs.

Theorem 2.3 (stochastic maximum principle). Let Hypothesis 2.2 hold and be an optimal pair. Then the optimal 6-tuple solves the above extended stochastic Hamiltonian system.

Proof. See, e.g., Yong and Zhou ( [6] , Chapter 3, pp. 123-137).

Hypothesis 2.4. The maps are locally Lipschitz in and their derivatives in are continuous in.

Theorem 2.5 (sufficient condition of optimality). Let Hypotheses 2.2 and 2.4 hold. Suppose that is an admissible 6-tuple, is concave P-a.s., and the maximum condition is verified, i.e.

then is an optimal pair.

Proof. See, e.g., Yong and Zhou ( [6] , Chapter 3, pp. 139-140).

Economic Implications of the First-Order Adjoint Process

The first-order order adjoint process plays a very important role in economic theory that we aim to introduce.

Assume that the framework of the previous section hold, but for simplicity all the variables are scalar-valued. For any initial data, consider the controlled system:

the set of admissible control:

and the functional:

where is the unique solution to the above state equation.

We have the following definition.

Definition 2.6. The following:

is called value function.

If and is also continuous, then:

At any time t suppose is slightly changed to, then by (3a) we obtain:

Intuitively, the above relation shows that the value of measures the rate of variation of the value function induced by an infinitesimal increase of the level of the state variable. Accordingly, the adjoint variable is called the marginal value of the optimal state.

Suppose the increment of the optimal state value could be purchased in the market. Thus, the adjoint variable represents the maximum price it would be worth paying an additional amount of the state. For this reason, the adjoint variable is also called the shadow price of the optimal state.

Moreover, if is the convex dual value function of V, then the following relationships apply (see, e.g., Di Giacinto et al. [7] ):

It is straightforward from the above relations to derive the following interesting relations:

3. FBSDEs and BSDEs

The conditions of optimality derived from the maximum principle are in general FBSDEs. In the following we review the basic results on the mathematical theory that explores existence uniqueness and solvability of these equations. In particular, we will focus on the interesting examples that will be useful in the applications when the FBSDE is decoupled, i.e. the forward component is independent from the backward one.

While for the forward component the basic properties are known from the theory of general stochastic processes, it is important to review the results about the backward component that is determined by a BSDE. Basic mathematical results on this class of equations are reviewed in the following subsection.

BSDEs

We consider the following controlled (unidimensional) BSDE:

where is the control variable -progressively measurable, the pair is the data, called the terminal condition and the generator (or driver), respectively. Here, x is -measurable and is -measurable for all. A solution to the above BSDE is a pair of processes adapted to the given filtration.

Solving an ordinary differential equation formulated as a terminal value problem or an initial value problem on the same time interval is equivalent under the time-reversing transformation,. On the contrary, when we are looking for a solution to a BSDE that is adapted respect to the given filtration we know only the past information, but are not able to know what is going on to occur in the future. This means that it is not possible simply reverse the time to find a solution as it would remove the adaptiveness.

The theory of BSDEs is an active domain of research due to its connections with stochastic control, mathematical finance, and partial differential equations (see, e.g., the survey elaborated by El Karoui et al. [8] ).

BSDEs were pioneered by Bismut [9] [10] [11] within the stochastic optimal control. Existence and uniqueness of adapted solutions were first proved by Pardoux and Peng [12] in a Brownian setting, with integrable terminal condition and generators that fulfill a canonical Lipschitz assumption in the space variables. These results were extended by Lepeltier and San Martín [13] to continuous drivers with linear growth. Since then, a lot of efforts have been made in order to study the well posedness of these equations.

With regards to the quadratic growth BSDEs, the proof of the existence of solutions are easily obtained under rather general framework. Conversely, the question of the uniqueness results is not trivial and requires stronger assumptions on the solutions.

In the standard framework where the terminal condition is bounded, in Kobylanski [14] obtained for the martingale case the uniqueness of the solution under some Lipschitz conditions. This seminal result has been extended in several directions, among others, by Tevzadze [15] , Morlais [16] , and Briand-Elie [17] . In Tevzadze [15] as well as in Briand-Elie [17] the proof in Kobylanski [14] was simplified in some less general cases. In particular, Tevzadze [15] showed the proof for the existence and uniqueness of a bounded solution in the Lipschitz quadratic, while Briand-Elie [17] incorporated delayed quadratic BSDEs, whose generator depends on the recent past of the component of the solution. Morlais [16] treated quadratic BSDEs driven by a continuous martingale and applied the corresponding results to study some related utility maximization problem.

Concerning the quadratic BSDEs with unbounded terminal value, Briand-Hu [18] first showed the existence of solution, while Barrieu and El Karoui [19] revisited the existence by a direct forward method that does not apply the result of Kobylanski [14] .

Barrieu and El Karoui [19] also showed some uniqueness results that require stronger assumptions comparing with the other existing literature. However, their alternative approach has other potential applications as numerical simulations of quadratic BSDEs and risk measures and their dual representations.

The uniqueness result was tackled for the first time by Briand and Hu [20] under some restrictive conditions. In particular, the authors proved the uniqueness among the solutions of quadratic BSDEs with convex generators and unbounded terminal conditions which admit every exponential moments. Delbaen et al. [21] [22] , strengthened the above result proving that uniqueness holds among solutions reducing the required order of exponential moments.

4. Simulation Approach to BSDE

The simulation approach to BSDE applying Longstaff and Schwartz [1] has been investigated the first time by Bouchard and Touzi [2] . In order to consider the numerical solution of a BSDE consider for simmplicity the simple case where the backward component is a pure martingale, i.e. the drift component is zero. Then the determination of and can be performed observing at time T, and can be obtained simulating forward.

Given the random variable, computation of can be performed computing a conditional expectation using the simulation approach of Longstaff and Schwartz [1] relying on a linear regression and moving backward. The exact choice of the polynomial family is relevant only for the numerical efficiency. For example, using Hermite polynomials we obtain:

where the standard error of the residual is of order:

where is the conditional variance.

Note that the approximation of can be simply derived considering a term by term differentiation:

We will show that in many cases, in particular in the applications that we will discuss later, the approximation to the pair:

provides an approximate solution of the BSDE that can be used to recover the optimal control strategies using the optimality conditions.

5. Applications

In the following applications we will always consider a standard Black & Scholes market, i.e. a risk free bond, whose dynamics are:

where is the risk free rate, and a traded risky asset (stock market index) whose dynamics are:

where is the volatility, is the drift, and is the instantaneous risk premium, i.e. the Sharpe ratio.

5.1. Classical Merton Problem

In this subsection we illustrate the application of the stochastic maximum principle to the solution of the well-known Merton’s portfolio allocation and consumption problem (see, e.g., Merton [23] ) when the utility has a constant relative risk-aversion as in Merton [24] .

The state equation describing the dynamics of wealth is given by:

where is the investment strategy representing the proportion of wealth invested in the risky asset, is the consumption stategy, and are the market parameters of the Black-Scholes model.

The preferences of the agent are given by the following utility function:

while the bequest function is:

The objective functional we aim to optimize is given by:

where the control belongs to the following set of admissible strategies:

and is the subjective discount factor.

The problem we aim to solve is to find an optimal control such that:

The stochastic Hamiltonian system associated to the problem is:

while the second-order adjoint equation is:

It can be easily found that:

with:

is the only adapted solution to the second-order adjoint equation. Observe that for all.

The corresponding H-function is:

and the maximum condition is:

It is easy to show that the function is concave and has a unique maximum point on given by:

(4)

if and only if for all.

Remark 5.1. By the stochastic maximum principle theorem a necessary conditions for to be optimal is:

Clearly the first-order adjoint process y―which represents the risk factor―plays a role in turning the convex function, i.e. the Hamiltonian along the optimal state trajectory, into a concave one, i.e. the H-function along the optimal state trajectory, previously shown.

By plugging (5a)-(5b) into the stochastic Hamiltonian system, we obtain:

Once we find an adapted solution of the above, a candidate for the optimal control can be deduced by (4). To this regard, suppose is an adapted solution and take the following guess function for:

for some deterministic function.

We would like to determine the equation that should satisfy. To this end, we differentiate the above using Itô formula:

i.e.

Comparing the above relation with the first-adjoint equation that appears in the previous stochastic Hamiltonian system, we have:

The above (6b)-(6c) is a Bernoulli equation whose solution is:

(7)

where and subject to the condition, i.e.

(8)

With the position of stated in (7) and the condition on ρ stated in (8), the pair of processes:

is the only square integrable adapted solution to the first-order adjoint equation.

The necessary condition required from the stochastic maximum principle theorem leads to the following policy:

(9)

and state trajectory (simply obtained substituting the expression of the above control policy to the stochastic Hamiltonian system):

(10)

as candidates to be an optimal pair.

Finally, a simple check shows that sufficient conditions of optimality are satisfied. We conclude that pointed out to (10) and (9) is an optimal pair of the classical Merton optimization problem. Therefore, the following result holds.

Proposition 5.2. There exists a unique optimal investment and consumption strategy given by:

where must satisfies the following Bernoulli ordinary differential equation:

whose solution is:

under the condition:

5.2. Merton Model with a Retirement Endowment

In a standard Black & Scholes market, the dynamics of wealth are given by:

where are the usual market parameters.

In addition, the agent is endowed with a nontradable investment I in a pension fund whose risky component is perfectly correlated with the market index. The dynamics of the investement are given by:

where in order to avoid arbitrage opportunities.

The preferences of the agent are described by the following constant relative risk aversion utility function:

while the bequest function is:

where, for simplicity, we set the parameter of the classical Merton problem.

At the terminal date T, the investor receives a liquid amount determined by an initial endowment which is invested in the stock market and cannot be withdrawn until maturity. In her hedging policy, the investor must optimally take into account both the liquid investment which can be continuously traded, and the lump sum amount that she will receive at the final date T and is also exposed to stock market risk fluctuations.

Given the functional:

and the set of admissible control strategies:

our goal is to find an optimal control such that:

The agent keeps into account her investment for retirement in her optimal allocation policy of the disposable wealth x_{0}.

The solution procedure is parallel to the standard Merton problem even though the allocation policy changes and an hedging term takes into account the stochastic opportunity set due to the presence of the investment I for retirement.

By applying the stochastic maximum principle we find the following stochastic Hamiltonian system:

The results can be tested since it is possible to evaluate the optimal proportion. As a matter of fact, due to the market completness it is possible to estimate the exact value of the optimal proportion as:

Following the line applied for solving the classical Merton problem, we can also determine the optimal consumption strategy, and so the optimal control is:

where the deterministic function has been previously defined in (7) in the classical Merton problem.

Now, we show how to derive numerically the allocation exploiting the following observation: due to market completeness the dynamics of the dual variable is known and corresponds to a log-normal process with drift equal to the opposite of the risk free rate r and volatility equal to the Sharpe ratio l.

Given the above observation, the above Hamiltonian system the control variables appear only in the forward equation. By using the observation on duality and Riesz representation (see, e.g., Yong and Zhou [6] , chapter 7, p. 353) it is possible to switch the backward and forward equations. Therefore, in this latter representation the stochastic Hamiltonian system can be rewritten as:

where is determined by the condition

(11)

We need to solve a FBSDE. The forward process is determined by the adjoint process, i.e. deflator whose expression is known. The backward components which must be obtained are and they uniquely determine the optimal pair as functions of the forward process.

The numerical determination of the optimal investment strategy can be obtained considering a least square regression procedure for.

We choose the following parameters:

, , , , , ,

, , , , , ,

where NP stands for number of paths considered and K is the number of discretization steps. In this case, the theoretical and numerical values averaged over 50 evaluations of the extended Merton proportion are:

and the relative error is of order 2.20%.

5.3. Income Drawdown Option with Minimum Guarantee

In this subsection we consider a second interesting application of the same approach that provides a numerical way to extend the analysis carried out in Di Giacinto et al. [7] which is briefly described below. The numerical extension we explore is necessary in order to assess the impact of the running cost as defined in the model.

We consider a defined contribution pension plan for a single representative participant during the decumulation phase, i.e. after retirement. Depending on the law and the rules of the scheme in many countries the retiree is allowed:

1) to defer annuitization at some time after retirement (this option is named, e.g., in UK “income drawdown option”, in US “phased withdrawals”);

2) to withdraw periodic income from the fund;

3) to invest the rest of the fund in the period between retirement and annuitization.

Thus, the pensioner has three degrees of freedom:

1) to decide when to annuitize (if ever);

2) to decide how much of the fund to withdraw at any time between retirement and ultimate annuitization (if any);

3) to decide what investment strategy to adopt in investing the fund at her disposal.

The first choice (i.e. the optimal time of annuitization) can be tackled by defining an optimal stopping time problem. The last two choices (i.e. the optimal consumption and investment selection) represent a classical investment and consumption problem, which can be solved using, e.g., stochastic optimal control techniques.

The contribution by Di Giacinto et al. [7] is the solution in closed form of the optimal control problem with constraints regarding the third choice (i.e. the investment choice) using a quadratic loss function and by applying the dynamic programming approach.

In a standard Black & Scholes market, the state equation is given by:

where is the process representing the fund wealth of the pensioner, x_{0} is the fund wealth at the retirement date, is the dollar amount invested in the risky asset, b_{0} is the consumption rate of the pensioner, and are the usual market parameters.

The objective functional is given by:

where and x is the sum of the subjective discount factor ρ and of the force of mortality δ, i.e..

The loss function has the following quadratic shape:

while is the target function. We choose:

where.

Constraints on the strategies and on the wealth are considered in the model. More precisely:

i) “no short selling” on the strategies:

ii) “no ruin” on the wealth:

with the “safety level” set as:

where.

Given the set of admissible controls:

the optimal control problem to be solved is to find an optimal control such that:

We are able to express the related stochastic Hamiltonian system as:

where y_{0} is determined following the same procedure applied in (11).

Note that the above Hamiltonian system is linear, thus the Riesz duality can be applied to switch the forward and the backward components. While for the Hamiltonian system is represented by a decoupled FBSDE, so the previous numerical approach can be applied, for the FBSDE is coupled but still linear and therefore numerically solvable by applying a different numerical approach, e.g., the four step scheme illustrated by Ma et al. [25] .

In the following we illustrate the simulation paths of the optimal investment strategy for and some realistic parametrization. The parameter chosen are:

where j is the retiree age, b_{0} is computed by applying the projected Italian males life table and NP stands for performed number of paths. The final target and the final guarantee are subjective and determined by the pensioner’s risk aversion.

6. Conclusions

In the above analysis, we showed how to use the stochastic maximum principle to reformulate conventional portfolio consumption models in order to explore their solution through the numerical solution of a linear stochastic hamiltonian system. We showed in a number of applications the flexibility of this methodology.

In the next future we aim to explore the performance of this method in two important extensions: a multidimensional state variable problems and in presence of optimal stopping control strategies. This last extension is particularly relevant for actuarial applications.

References

[1] Longstaff, F.A. and Schwartz, E.S. (2001) Valuing American Options by Simulation: A Simple Least Squares Approach. The Review of Financial Studies, 14, 113-147.

https://doi.org/10.1093/rfs/14.1.113

[2] Bouchard, B. and Touzi, N. (2004) Discrete-Time Approximation and Monte-Carlo Simulation of Backward Stochastic Differential Equations. Stochastic Processes and their Applications, 111, 175-206.

https://doi.org/10.1016/j.spa.2004.01.001

[3] Gobet, E., Lemor, J.-P. and Warin, X. (2005) A Regression-Based Monte Carlo Method to Solve Backward Stochastic Differential Equations. The Annals of Applied Probability, 15, 2172-2202.

https://doi.org/10.1214/105051605000000412

[4] Lemor, J.-P., Gobet, E. and Warin, X. (2006) Rate of Convergence of an Empirical Regression Method for Solving Generalized Backward Stochastic Differential Equations. The Annals of Applied Probability, 12, 889-916.

[5] Bender, C. and Steiner, J. (2012) Least-Squares Monte Carlo for backward SDEs. Carmona, R.A., Del Moral, P., Hu, P. and Oudjane, N., Eds., Numerical Methods in Finance, Vol. 12, Springer Proceedings in Mathematics, Springer, Berlin Heidelberg, 257-289.

https://doi.org/10.1007/978-3-642-25746-9_8

[6] Yong, J. and Zhou, X.Y. (1999) Stochastic Controls: Hamiltonian Systems and HJB Equations. Springer Verlag, New York.

https://doi.org/10.1007/978-1-4612-1466-3

[7] Di Giacinto, M., Federico, S., Gozzi, F. and Vigna, E. (2014) Income Drawdown Option with Minimum Guarantee. European Journal of Operational Research, 234, 610-624.

https://doi.org/10.1016/j.ejor.2013.10.026

[8] El Karoui, N., Peng, S. and Quenez, M.-C. (1997) Backward Stochastic Differential Equations in Finance. Mathematical Finance, 7, 1-71.

https://doi.org/10.1111/1467-9965.00022

[9] Bismut, J.-M. (1973) Conjugate Convex Function in Optimal Stochastic Control. Journal of Mathematical Analysis and Applications, 44, 384-404.

https://doi.org/10.1016/0022-247X(73)90066-8

[10] Bismut, J.-M. (1978) An Introductory Approach to Duality in Optimal Stochastic Control. SIAM Review, 20, 62-78.

https://doi.org/10.1137/1020004

[11] Bismut, J.-M. (1978b) Contrôle des systèmes linaires quadratiques: Applications de l’intégrale stochastique. Séminaire de probabilités de Strasbourg, 12, 180-264.

[12] Pardoux, E. and Peng, S.G. (1990) Adapted Solution of a Backward Stochastic Differential Equation. Systems & Control Letters, 14, 55-61.

https://doi.org/10.1016/0167-6911(90)90082-6

[13] Lepeltier, J.-P. and San Martín, J. (1997) Backward Stochastic Differential Equations with Continuous Coeffcient. Statistics & Probability Letters, 32, 425-430.

https://doi.org/10.1016/S0167-7152(96)00103-4

[14] Kobylanski, M. (2000) Backward Stochastic Differential Equations and Partial Differential Equations with Quadratic Growth. The Annals of Probability, 28, 558-602.

https://doi.org/10.1214/aop/1019160253

[15] Tevzadze, R. (2008) Solvability of Backward Stochastic Differential Equations with Quadratic Growth. Stochastic Processes and their Applications, 118, 503-515.

https://doi.org/10.1016/j.spa.2007.05.009

[16] Morlais, M.-A. (2009) Quadratic BSDEs Driven by a Continuous Martingale and Applications to the Utility Maximization Problem. Finance and Stochastics, 13, 121-150.

https://doi.org/10.1007/s00780-008-0079-3

[17] Briand, P. and Elie, R. (2013) A Simple Constructive Approach to Quadratic BSDEs with or without Delay. Stochastic Processes and Their Applications, 123, 2921-2939.

https://doi.org/10.1016/j.spa.2013.02.013

[18] Briand, P. and Hu, Y. (2006) BSDE with Quadratic Growth and Unbounded Terminal Value. Probability Theory and Related Fields, 136, 604-618.

https://doi.org/10.1007/s00440-006-0497-0

[19] Barrieu, P. and El Karoui, N. (2013) Monotone Stability of Quadratic Semimartingales with Applications to Unbounded General Quadratic BSDEs. Annals of Probability, 41, 1831-1863.

https://doi.org/10.1214/12-AOP743

[20] Briand, P. and Hu, Y. (2008) Quadratic BSDEs with Convex Generator and Unbounded Terminal Conditions. Probability Theory and Related Fields, 141, 543-567.

https://doi.org/10.1007/s00440-007-0093-y

[21] Delbaen, F., Hu, Y., and Richou, A. (2011) On the Uniqueness of Solutions to Quadratic BSDEs with Convex Generators and Unbounded Terminal Conditions. Annales de l’Institut Henri Poincaré—Probabilités et Statistiques, 47, 559-574.

https://doi.org/10.1214/10-AIHP372

[22] Delbaen, F., Hu, Y., and Richou, A. (2015) On the Uniqueness of Solutions to Quadratic BSDEs with Convex Generators and Unbounded Terminal conditions: The Critical Case. Discrete and Continuous Dynamical Systems, 35, 5273-5283.

https://doi.org/10.3934/dcds.2015.35.5273

[23] Merton, R. (1992) Continuous-Time Finance. United States of America, Blackwell.

[24] Merton, R.C. (1969) Lifetime Portfolio Selection under Uncertainty: The Continuous-Time Case. The Review of Economics and Statistics, 51, Article ID: 24757.

https://doi.org/10.2307/1926560

[25] Ma, J., Protter, P., and Yong, J. (1994) Solving Forward-Backward Stochastic Differential Equations Explicitly—A Four Step Scheme. Probability Theory and Related Fields, 98, 339-359.

https://doi.org/10.1007/BF01192258