Algebraic solution of the Stein–Stein model for stochastic volatility

We provide an algebraic approach to the solution of the Stein–Stein model for stochastic volatility which arises in the determination of the Radon–Nikodym density of the minimal entropy of the martingale measure. We extend our investigation to the case in which the parameters of the model are time-dependent. Our algorithmic approach obviates the need for Ansätze for the structure of the solution.     

Smooth Solutions to Optimal Investment Models with Stochastic Volatilities and Portfolio Constraints

Abstract. This paper deals with an extension of Merton's optimal investment problem to a multidimensional model with stochastic volatility and portfolio constraints. The classical dynamic programming approach leads to a characterization of the value function as a viscosity solution of the highly nonlinear associated Bellman equation. A logarithmic transformation expresses the value function in terms of the solution to a semilinear parabolic equation with quadratic growth on the derivative term. Using a stochastic control representation and some approximations, we prove the existence of a smooth solution to this semilinear equation. An optimal portfolio is shown to exist, and is expressed in terms of the classical solution to this semilinear equation. This reduction is useful for studying numerical schemes for both the value function and the optimal portfolio. We illustrate our results with several examples of stochastic volatility models popular in the financial literature.     

Smooth Solutions to Optimal Investment Models with Stochastic Volatilities and Portfolio Constraints

   Abstract. This paper deals with an extension of Merton's optimal investment problem to a multidimensional model with stochastic volatility and portfolio constraints. The classical dynamic programming approach leads to a characterization of the value function as a viscosity solution of the highly nonlinear associated Bellman equation. A logarithmic transformation expresses the value function in terms of the solution to a semilinear parabolic equation with quadratic growth on the derivative term. Using a stochastic control representation and some approximations, we prove the existence of a smooth solution to this semilinear equation. An optimal portfolio is shown to exist, and is expressed in terms of the classical solution to this semilinear equation. This reduction is useful for studying numerical schemes for both the value function and the optimal portfolio. We illustrate our results with several examples of stochastic volatility models popular in the financial literature.     

Risk premium and fair option prices under stochastic volatility: the HARA solution

We have solved the problem of finding (HARA) fair option price under a general stochastic volatility model. For a given HARA utility, the ‘risk premium’, i.e., the ‘market price of volatility risk’ is determined via a solution of a certain nonlinear PDE. Equivalently, the fair option price is determined as a solution of an uncoupled system of a non-linear PDE and a Black–Scholes type PDE. To cite this article: S.D. Stojanovic, C. R. Acad. Sci. Paris, Ser. I 340 (2005).     

On the multi-dimensional portfolio optimization with stochastic volatility

In a recent paper by Mnif [18], a solution to the portfolio optimization with stochastic volatility and constraints problem has been proposed, in which most of the model parameters are time-homogeneous. However, there are cases where time-dependent parameters are needed, such as in the calibration of financial models. Therefore, the purpose of this paper is to generalize the work of Mnif [18] to the time-inhomogeneous case. We consider a time-dependent exponential utility function of which the objective is to maximize the expected utility from the investor’s terminal wealth. The derived Hamilton-Jacobi-Bellman(HJB) equation, is highly nonlinear and is reduced to a semilinear partial differential equation (PDE) by a suitable transformation. The existence of a smooth solution is proved and a verification theorem presented. A multi-asset stochastic volatility model with jumps and endowed with time-dependent parameters is illustrated.     

Non-linear filtering and optimal investment under partial information for stochastic volatility models

This paper studies the question of filtering and maximizing terminal wealth from expected utility in partial information stochastic volatility models. The special feature is that the only information available to the investor is the one generated by the asset prices, and the unobservable processes will be modeled by stochastic differential equations. Using the change of measure techniques, the partial observation context can be transformed into a full information context such that coefficients depend only on past history of observed prices (filter processes). Adapting the stochastic non-linear filtering, we show that under some assumptions on the model coefficients, the estimation of the filters depend on a priori models for the trend and the stochastic volatility. Moreover, these filters satisfy a stochastic partial differential equations named “Kushner–Stratonovich equations”. Using the martingale duality approach in this partially observed incomplete model, we can characterize the value function and the optimal portfolio. The main result here is that, for power and logarithmic utility, the dual value function associated to the martingale approach can be expressed, via the dynamic programming approach, in terms of the solution to a semilinear partial differential equation which depends on the filters estimate and the volatility. We illustrate our results with some examples of stochastic volatility models popular in the financial literature.     

模糊集理论在随机波动率期权定价中的应用

目的是对基于随机波动率模型的期权定价问题应用模糊集理论.主要思想是把波动率的概率表示转换为可能性表示,从而把关于股票价格的带随机波动率的随机过程简化为带模糊参数的随机过程.然后建立非线性偏微分方程对欧式期权进行定价.     

Comparison Method of Stability Analysis of Semilinear Time-Delay Stochastic Evolution Equation

Abstract

In this paper, we set up the comparison theorem between the mild solution of semilinear time-delay stochastic evolution equation with general time-delay variable and the solution of a class (1-dimension) deterministic functional differential equation, by using the Razumikhin–Lyapunov type functional and the theory of functional differential inequalities. By applying this comparison theorem, we give various types of the stability comparison criteria for the semilinear time-delay stochastic evolution equations. With the aid of these comparison criteria, one can reduce the stability analysis of semilinear time-delay stochastic evolution equations in Hilbert space to that of a class (1-dimension) deterministic functional differential equations. Furthermore, these comparison criteria in special case have been applied to derive sufficient conditions for various stability of the mild solution of semilinear time-delay stochastic evolution equations. Finally, the theories are illustrated with some examples.     

Moments and Mellin transform of the asset price in Stein and Stein model and option pricing

In this paper, we derive closed formulas for moments and Mellin transform of the asset price in the stochastic volatility Stein and Stein model. Next, we present applications of our results to pricing power and self-quanto options using numerical methods.     

Efficient L-stable method for parabolic problems with application to pricing American options under stochastic volatility

Efficient L-stable numerical method for semilinear parabolic problems with nonsmooth initial data is proposed and implemented to solve Heston’s stochastic volatility model based PDE for pricing American options under stochastic volatility. The proposed new method is also used to solve two asset American options pricing problem. Cox and Matthews [S.M. Cox, P.C. Matthews, Exponential time differencing for stiff systems, Journal of Computational Physics 176 (2002) 430-455] developed a class of exponential time differencing Runge-Kutta schemes (ETDRK) for nonlinear parabolic problems. Kassam and Trefethen [A.K. Kassam, L.N. Trefethen, Fourth-order time stepping for stiff PDEs, SIAM Journal on Scientific Computing 26 (4) (2005) 1214-1233] showed that while computing certain functions involved in the Cox-Matthews schemes, severe cancelation errors can occur which affect the accuracy and stability of the schemes. Kassam and Trefethen proposed complex contour integration technique to implement these schemes in a way that avoids these cancelation errors. But this approach creates new difficulties in choosing and evaluating the contour integrals for larger problems. We modify the ETDRK schemes using positivity preserving Padé approximations of the matrix exponential functions and construct computationally efficient parallel version using splitting technique. As a result of this approach it is required only to solve several backward Euler linear problems in serial or parallel.     

Ergodicity and Kolmogorov Equations for Dissipative SPDEs with Singular Drift: a Variational Approach

Potential Analysis - We prove existence of invariant measures for the Markovian semigroup generated by the solution to a parabolic semilinear stochastic PDE whose nonlinear drift term satisfies...     

Probabilistic approach to homogenization of a non-divergence form semilinear PDE with non-periodic coefficients

We consider a semilinear partial differential equation (PDE) of non-divergence form perturbed by a small parameter. We then study the asymptotic behavior of Sobolev solutions in the case where the coefficients admit limits in C?esaro sense. Neither periodicity nor ergodicity will be needed for the coefficients. In our situation, the limit (or averaged or effective) coefficients may have discontinuity. Our approach combines both probabilistic and PDEs arguments. The probabilistic one uses the weak convergence of solutions of backward stochastic differential equations (BSDE) in the Jakubowski S-topology, while the PDEs argument consists to built a solution, in a suitable Sobolev space, for the PDE limit. We finally show the existence and uniqueness for the associated averaged BSDE, then we deduce the uniqueness of the limit PDE from the uniqueness of the averaged BSDE.     

Singular perturbations to semilinear stochastic heat equations

We consider a class of singular perturbations to the stochastic heat equation or semilinear variations thereof. The interesting feature of these perturbations is that, as the small parameter ε tends to zero, their solutions converge to the ‘wrong’ limit, i.e. they do not converge to the solution obtained by simply setting ε?=?0. A similar effect is also observed for some (formally) small stochastic perturbations of a deterministic semilinear parabolic PDE. Our proofs are based on a detailed analysis of the spatially rough component of the equations, combined with a judicious use of Gaussian concentration inequalities.     

A class of portfolio selection with a four-factor futures price model

Considering the stochastic exchange rate, a four-factor futures model with the underling asset, convenience yield, instantaneous risk free interest rate and exchange rate, is established. These processes follow jump-diffusion processes (Weiner process and Poisson process). The corresponding partial differential equation (PDE) of the futures price is derived. The general solution of the PDE with parameters is drawn. The weight least squares approach is applied to obtain the parameters of above PDE. Variance is substituted by semi-variance in Markowitzs portfolio selection model. Therefore, a class of multi-period semi-variance model is formulated originally. Then, a continuous-time mean-variance portfolio model is also considered. The corresponding stochastic Hamilton-Jacobi-Bellman (HJB) equation of the problem with nonlinear constraints is derived. A numerical algorithm is proposed for finding the optimal solution in this paper. Finally, in order to demonstrate the effectiveness of the theoretical models and numerical methods, the fuel futures in Shanghai exchange market and the Brent crude oil futures in London exchange market are selected to be examples.     

Pricing Volatility Swaps Under Heston's Stochastic Volatility Model with Regime Switching

A model is developed for pricing volatility derivatives, such as variance swaps and volatility swaps under a continuous‐time Markov‐modulated version of the stochastic volatility (SV) model developed by Heston. In particular, it is supposed that the parameters of this version of Heston's SV model depend on the states of a continuous‐time observable Markov chain process, which can be interpreted as the states of an observable macroeconomic factor. The market considered is incomplete in general, and hence, there is more than one equivalent martingale pricing measure. The regime switching Esscher transform used by Elliott et al. is adopted to determine a martingale pricing measure for the valuation of variance and volatility swaps in this incomplete market. Both probabilistic and partial differential equation (PDE) approaches are considered for the valuation of volatility derivatives.     

The Minimal Entropy Martingale Measure and Numerical Option Pricing for the Barndorff–Nielsen–Shephard Stochastic Volatility Model

Abstract

We develop and apply a numerical scheme for pricing options in the stochastic volatility model proposed by Barndorff–Nielsen and Shephard. This non-Gaussian Ornstein–Uhlenbeck type of volatility model gives rise to an incomplete market, and we consider the option prices under the minimal entropy martingale measure. To numerically price options with respect to this risk neutral measure, one needs to consider a Black and Scholes type of partial differential equation, with an integro-term arising from the volatility process. We suggest finite difference schemes to solve this parabolic integro-partial differential equation, and derive appropriate boundary conditions for the finite difference method. As an application of our algorithm, we consider price deviations from the Black and Scholes formula for call options, and the implications of the stochastic volatility on the shape of the volatility smile.     

Robust investment–reinsurance optimization with multiscale stochastic volatility

This paper investigates the investment and reinsurance problem in the presence of stochastic volatility for an ambiguity-averse insurer (AAI) with a general concave utility function. The AAI concerns about model uncertainty and seeks for an optimal robust decision. We consider a Brownian motion with drift for the surplus of the AAI who invests in a risky asset following a multiscale stochastic volatility (SV) model. We formulate the robust optimal investment and reinsurance problem for a general class of utility functions under a general SV model. Applying perturbation techniques to the Hamilton–Jacobi–Bellman–Isaacs (HJBI) equation associated with our problem, we derive an investment–reinsurance strategy that well approximates the optimal strategy of the robust optimization problem under a multiscale SV model. We also provide a practical strategy that requires no tracking of volatility factors. Numerical study is conducted to demonstrate the practical use of theoretical results and to draw economic interpretations from the robust decision rules.     

Stochastic viscosity solutions for nonlinear stochastic partial differential equations. Part I

This paper, together with the accompanying work (Part II, Stochastic Process. Appl. 93 (2001) 205–228) is an attempt to extend the notion of viscosity solution to nonlinear stochastic partial differential equations. We introduce a definition of stochastic viscosity solution in the spirit of its deterministic counterpart, with special consideration given to the stochastic integrals. We show that a stochastic PDE can be converted to a PDE with random coefficients via a Doss–Sussmann-type transformation, so that a stochastic viscosity solution can be defined in a “point-wise” manner. Using the recently developed theory on backward/backward doubly stochastic differential equations, we prove the existence of the stochastic viscosity solution, and further extend the nonlinear Feynman–Kac formula. Some properties of the stochastic viscosity solution will also be studied in this paper. The uniqueness of the stochastic viscosity solution will be addressed separately in Part II where the relation between the stochastic viscosity solution and the ω-wise, “deterministic” viscosity solution to the PDE with random coefficients will be established.     

American Options in the Heston Model with Stochastic Interest Rate and Its Generalizations

Abstract

We consider the Heston model with the stochastic interest rate of Cox–Ingersoll–Ross (CIR) type and more general models with stochastic volatility and interest rates depending on two CIR-factors; the price, volatility and interest rate may correlate. Time-derivative and infinitesimal generator of the process for factors that determine the dynamics of the interest rate and/or volatility are discretized. The result is a sequence of embedded perpetual options arising in the time discretization of a Markov-modulated Lévy model. Options in this sequence are solved using an iteration method based on the Wiener–Hopf factorization. Typical shapes of the early exercise boundary are shown, and good agreement of option prices with prices calculated with the Longstaff–Schwartz method and Medvedev–Scaillet asymptotic method is demonstrated.     

Energy methods for stochastic differential equations

This article is devoted to the existence of strong solutions to stochastic differential equations (SDEs). Compared with Ito's theory, we relax the assumptions on the volatility term and replace the global Lipschitz continuity condition with a local Lipschitz continuity condition and a Hoelder continuity condition. In particular, our general SDE covers the Cox–Ingersoll–Ross SDE as a special case. We note that the general weak existence theory presumably extends to our general SDE (although the explicit time dependence of the drift term and the volatility term might require some extra considerations). However, avoiding weak existence theory we prove the existence of a strong solution directly using a priori estimates (the so-called energy estimates) derived from the SDE. The benefit of this approach is that the argument only requires some basic knowledge about stochastic and functional analysis. Moreover, the underlying principle has developed to become one of the cornerstones of the modern theory of partial differential equations (PDEs). In this sense, the general goal of this article is not just to establish the existence of a strong solution to the SDE under consideration but rather to introduce a new principle in the context of SDEs that has already proven to be successful in the context of PDEs.