Search

In this paper, we propose a new approach for stochastic control problems arising from utility maximization. The main idea is to directly start from the dynamical programming equation and compute the conditional expectation using a novel representation of the conditional density function through the Dirac Delta function and the corresponding series representation. We obtain an explicit series representation of the value function, whose coefficients are expressed through integration of the value function at a later time point against a chosen basis function. Thus we are able to set up a recursive integration time-stepping scheme to compute the optimal value function given the known terminal condition, e.g. utility function. Due to tensor decomposition property of the Dirac Delta function in high dimensions, it is straightforward to extend our approach to solving high-dimensional stochastic control problems. The backward recursive nature of the method also allows for solving stochastic control and stopping problems, i.e. mixed control problems. We illustrate the method through solving some two-dimensional stochastic control (and stopping) problems, including the case under the classical and rough Heston stochastic volatility models, and stochastic local volatility models such as the stochastic alpha beta rho (SABR) model.

Rough volatility models have recently been empirically shown to provide a good fit to historical volatility time series and implied volatility smiles of SPX options. They are continuous-time stochastic volatility models, whose volatility process is driven by a fractional Brownian motion with Hurst parameter less than half. Due to the challenge that it is neither a semimartingale nor a Markov process, there is no unified method that not only applies to all rough volatility models, but also is computationally efficient. This paper proposes a semimartingale and continuous-time Markov chain (CTMC) approximation approach for the general class of rough stochastic local volatility (RSLV) models. In particular, we introduce the perturbed stochastic local volatility (PSLV) model as the semimartingale approximation for the RSLV model and establish its existence , uniqueness and Markovian representation. We propose a fast CTMC algorithm and prove its weak convergence. Numerical experiments demonstrate the accuracy and high efficiency of the method in pricing European, barrier and American options. Comparing with existing literature, a significant reduction in the CPU time to arrive at the same level of accuracy is observed., Comment: 29 pages

In this paper we study a utility maximization problem with both optimal control and optimal stopping in a finite time horizon. The value function can be characterized by a variational equation that involves a free boundary problem of a fully nonlinear partial differential equation. Using the dual control method, we derive the asymptotic properties of the dual value function and the associated dual free boundary for a class of utility functions, including power and non-HARA utilities. We construct a global closed-form approximation to the dual free boundary, which greatly reduces the computational cost. Using the duality relation, we find the approximate formulas for the optimal value function, trading strategy, and exercise boundary for the optimal investment stopping problem. Numerical examples show the approximation is robust, accurate and fast., Comment: 30 pages, 10 figures

The aim of this paper is to study the fast computation of the lower and upper bounds on the value function for utility maximization under the Heston stochastic volatility model with general utility functions. It is well known there is a closed form solution of the HJB equation for power utility due to its homothetic property. It is not possible to get closed form solution for general utilities and there is little literature on the numerical scheme to solve the HJB equation for the Heston model. In this paper we propose an efficient dual control Monte Carlo method for computing tight lower and upper bounds of the value function. We identify a particular form of the dual control which leads to the closed form upper bound for a class of utility functions, including power, non-HARA and Yarri utilities. Finally, we perform some numerical tests to see the efficiency, accuracy, and robustness of the method. The numerical results support strongly our proposed scheme., Comment: 25 pages, 2 figures

In this paper we propose a new robust algorithm to find the optimal static replicating portfolios for general nonlinear payoff functions and give the estimate of the rate of convergence that is absent in the literature. We choose the static replication by minimizing the error bound between the nonlinear payoff function and the linear spline approximation and derive the equidistribution equation for selecting the optimal strike prices. The numerical tests for variance swaps and swaptions and also for the static quadratic replication and the model with counterparty risk show that the proposed algorithm is simple, fast and accurate. The paper has generalized and improved the results of the static replication and approximation in the literature., Comment: 15 pages