Numerical valuation of discrete double barrier options

In the present paper we explore the problem for pricing discrete barrier options utilizing the Black-Scholes model for the random movement of the asset price. We postulate the problem as a path integral calculation by choosing approach that is similar to the quadrature method. Thus, the problem is reduced to the estimation of a multi-dimensional integral whose dimension corresponds to the number of the monitoring dates.We propose a fast and accurate numerical algorithm for its valuation. Our results for pricing discretely monitored one and double barrier options are in agreement with those obtained by other numerical and analytical methods in Finance and literature. A desired level of accuracy is very fast achieved for values of the underlying asset close to the strike price or the barriers.The method has a simple computer implementation and it permits observing the entire life of the option.     

An efficient algorithm for Bermudan barrier option pricing

An efficient option pricing method based on Fourier-cosine expansions was presented by Fang and Oosterlee for European options in 2008,and later,this method was also used by them to price early-exercis...     

A general framework for pricing Asian options under stochastic volatility on parallel architectures

In this paper, we present a transform-based algorithm for pricing discretely monitored arithmetic Asian options with remarkable accuracy in a general stochastic volatility framework, including affine models and time-changed Lévy processes. The accuracy is justified both theoretically and experimentally. In addition, to speed up the valuation process, we employ high-performance computing technologies. More specifically, we develop a parallel option pricing system that can be easily reproduced on parallel computers, also realized as a cluster of personal computers. Numerical results showing the accuracy, speed and efficiency of the procedure are reported in the paper.     

Efficient pricing of discrete Asian options

Asian options are popular path-dependent financial derivatives. This paper uses lattices to price fixed-strike European-style Asian options that are discretely monitored. The algorithm proposed can also be applied to floating-strike Asian options as well because fixed-strike and floating-strike Asian options are related through an equation. The discretely monitored version is usually found in practice instead of the continuously monitored version usually encountered in the literature. This paper presents the first provably quadratic-time convergent lattice algorithm for pricing fixed-strike European-style discretely monitored Asian options. It is the most efficient lattice algorithm with convergence guarantees. The algorithm relies on the Lagrange multipliers to choose the number of states for each node of the lattice. Extensive numerical experiments and comparisons with many existing numerical methods confirm the performance claims and the competitiveness of our algorithm. This result places fixed-strike European-style discretely monitored Asian options in the same complexity class as vanilla options.     

The discontinuous Galerkin method for discretely observed Asian options

Asian options represent an important subclass of the path-dependent contracts that are identified by payoff depending on the average of the underlying asset prices over the prespecified period of option lifetime. Commonly, this average is observed at discrete dates, and also, early exercise features can be admitted. As a result, analytical pricing formulae are not always available. Therefore, some form of a numerical approximation is essential for efficient option valuation. In this paper, we study a PDE model for pricing discretely observed arithmetic Asian options with fixed as well as floating strike for both European and American exercise features. The pricing equation for such options is similar to the Black-Scholes equation with 1 underlying asset, and the corresponding average appears only in the jump conditions across the sampling dates. The objective of the paper is to present the comprehensive methodological concept that forms and improves the valuation process. We employ a robust numerical procedure based on the discontinuous Galerkin approach arising from the piecewise polynomial generally discontinuous approximations. This technique enables a simple treatment of discrete sampling by incorporation of jump conditions at each monitoring date. Moreover, an American early exercise constraint is directly handled as an additional nonlinear source term in the pricing equation. The proposed solving procedure is accompanied by an empirical study with practical results compared to reference values.     

Lie-algebraic approach for pricing moving barrier options with time-dependent parameters

In this paper we apply the Lie-algebraic technique for the valuation of moving barrier options with time-dependent parameters. The value of the underlying asset is assumed to follow the constant elasticity of variance (CEV) process. By exploiting the dynamical symmetry of the pricing partial differential equations, the new approach enables us to derive the analytical kernels of the pricing formulae straightforwardly, and thus provides an efficient way for computing the prices of the moving barrier options. The method is also able to provide tight upper and lower bounds for the exact prices of CEV barrier options with fixed barriers. In view of the CEV model being empirically considered to be a better candidate in equity option pricing than the traditional Black-Scholes model, our new approach could facilitate more efficient comparative pricing and precise risk management in equity derivatives with barriers by incorporating term-structures of interest rates, volatility and dividend into the CEV option valuation model.     

Continuity correction for discrete barrier options with two barriers

Discrete barrier options are the options whose payoffs are determined by underlying prices at a finite set of times. We consider the discrete barrier option with two barriers. Broadie et al. (1997) [16] proposed a continuity correction for the discretely monitored barrier option. We extend this idea to barrier option with two barriers. The proof for discrete chained barrier option is provided and numerical results show the continuity correction approximation is remarkably accurate.     

Robust barrier option pricing by frame projection under exponential Lévy dynamics

We present an efficient method for robustly pricing discretely monitored barrier and occupation time derivatives under exponential Lévy models. This includes ordinary barrier options, as well as (resetting) Parisian options, delayed barrier options (also known as cumulative Parisian or Parasian options), fader options and step options (soft-barriers), all with single and double barriers, which have yet to be priced with more general Lévy processes, including KoBoL (CGMY), Merton’s jump diffusion and NIG. The method’s efficiency is derived in part from the use of frame-projected transition densities, which transform the problem into the Fourier domain and accelerate the convergence of intermediate expectations. Moreover, these expectations are approximated by Toeplitz matrix-vector multiplications, resulting in a fast implementation. We devise an augmentation approach that contributes to the method’s robustness, adding protection against mis-specifying a proper truncation support of the transition density. Theoretical convergence is verified by a series of numerical experiments which demonstrate the method’s efficiency and accuracy.     

Comment on ‘Correcting for Simulation Bias in Monte Carlo Methods to Value Exotic Options in Models Driven by Lévy Processes’ by C. Ribeiro and N. Webber

Abstract

This paper proposes a pricing method for path-dependent derivatives with discrete monitoring when an underlying asset price is driven by a time-changed Lévy process. The key to our method is to derive a backward recurrence relation for computing the multivariate characteristic function of the intertemporal joint distribution of the time-changed Lévy process. Using the derived representation of the characteristic function, we obtain semi-analytical pricing formulas for geometric Asian, forward start, barrier, fader and lookback options, all of which are discretely monitored.     

Option pricing with mean reversion and stochastic volatility

Many underlying assets of option contracts, such as currencies, commodities, energy, temperature and even some stocks, exhibit both mean reversion and stochastic volatility. This paper investigates the valuation of options when the underlying asset follows a mean-reverting lognormal process with stochastic volatility. A closed-form solution is derived for European options by means of Fourier transform. The proposed model allows the option pricing formula to capture both the term structure of futures prices and the market implied volatility smile within a unified framework. A bivariate trinomial lattice approach is introduced to value path-dependent options with the proposed model. Numerical examples using European options, American options and barrier options demonstrate the use of the model and the quality of the numerical scheme.     

NUMERICAL ANALYSIS ON BINOMIAL TREE METHODS FOR AMERICAN LOOKBACK OPTIONS

1　IntroductionLookback options are path-dependent options whose payoffs depend on the maximumor the minimum of the underlying asset price during the life of the options( see[6] [1 0 ][1 4] ) .Here the maximum or minimum realized asset price may be monitored either con-tinuously or discretely.An American lookback call( put) option allows to be exercised atany time prior to expiry and gives the holder the rightto buy( sell) atthe historical mini-mum( maximum) of the underlying asset price on ex…     

几何平均亚式期权定价方法的探析

本文对几何平均亚式期权不同的定价方法进行了详细的论述，从随机偏微分方程途径与概率论途径两个角度仔细描述了亚式期权定价的过程中，每个具体的主要演算步骤．本文采用几何平均法计算资产价格的平均值，并以连续时间的情形为例，用两种不同的方法得到几何平均亚式期权的解析定价公式，并通过比较得出两种结论是完全一致的．     

Convergence of the binomial tree method for Asian options in jump-diffusion models

The binomial tree methods (BTM), first proposed by Cox, Ross and Rubinstein [J. Cox, S. Ross, M. Rubinstein, Option pricing: A simplified approach, J. Finan. Econ. 7 (1979) 229-264] in diffusion models and extended by Amin [K.I. Amin, Jump diffusion option valuation in discrete time, J. Finance 48 (1993) 1833-1863] to jump-diffusion models, is one of the most popular approaches to pricing options. In this paper, we present a binomial tree method for Asian options in jump-diffusion models and show its equivalence to certain explicit difference scheme. Employing numerical analysis and the notion of viscosity solution, we prove the uniform convergence of the binomial tree method for European-style and American-style Asian options.     

Numerical pricing of options using high-order compact finite difference schemes

We consider high-order compact (HOC) schemes for quasilinear parabolic partial differential equations to discretise the Black–Scholes PDE for the numerical pricing of European and American options. We show that for the heat equation with smooth initial conditions, the HOC schemes attain clear fourth-order convergence but fail if non-smooth payoff conditions are used. To restore the fourth-order convergence, we use a grid stretching that concentrates grid nodes at the strike price for European options. For an American option, an efficient procedure is also described to compute the option price, Greeks and the optimal exercise curve. Comparisons with a fourth-order non-compact scheme are also done. However, fourth-order convergence is not experienced with this strategy. To improve the convergence rate for American options, we discuss the use of a front-fixing transformation with the HOC scheme. We also show that the HOC scheme with grid stretching along the asset price dimension gives accurate numerical solutions for European options under stochastic volatility.     

Pricing annuity guarantees under a double regime-switching model

This paper is concerned with the valuation of equity-linked annuities with mortality risk under a double regime-switching model, which provides a way to endogenously determine the regime-switching risk. The model parameters and the reference investment fund price level are modulated by a continuous-time, finite-time, observable Markov chain. In particular, the risk-free interest rate, the appreciation rate, the volatility and the martingale describing the jump component of the reference investment fund are related to the modulating Markov chain. Two approaches, namely, the regime-switching Esscher transform and the minimal martingale measure, are used to select pricing kernels for the fair valuation. Analytical pricing formulas for the embedded options underlying these products are derived using the inverse Fourier transform. The fast Fourier transform approach is then used to numerically evaluate the embedded options. Numerical examples are provided to illustrate our approach.     

A systematic approach to pricing and hedging international derivatives with interest rate risk: analysis of international derivatives under stochastic interest rates

This paper deals with the valuation and the hedging of non-path-dependent European options on one or several underlying assets in a model of an international economy allowing for both, interest rate risk and exchange rate risk. Using martingale theory and, in particular, the change of numeraire technique we provide a unified and easily applicable approach to pricing and hedging exchange options on stocks, bonds, futures, interest rates and exchange rates. We also cover the pricing and hedging of compound exchange options.     

Analytic solution for American barrier options with two barriers

This paper concerns American barrier options with two barriers. Standard American Options are difficult to price but there exist good numerical or analytical approximation methods. The situation is different for American barrier options. These options cease to exist or come into being if some price barrier is hit during the option's life. The paper studies analytic valuation of American barrier options with two barriers where the barriers become active by turns. In this paper, analytic valuation formulas for these options are derived by using both constant and exponential barriers for optimal early exercise policies.     

Asymptotic expansion for pricing options for a mean-reverting asset with multiscale stochastic volatility

This work investigates the valuation of options when the underlying asset follows a mean-reverting log-normal process with a stochastic volatility that is driven by two stochastic processes with one persistent factor and one fast mean-reverting factor. Semi-analytical pricing formulas for European options are derived by means of multiscale asymptotic techniques. Numerical examples demonstrate the use of the model and the quality of the numerical scheme.     

Convergence analysis of a monotonic penalty method for American option pricing

This paper is devoted to study the convergence analysis of a monotonic penalty method for pricing American options. A monotonic penalty method is first proposed to solve the complementarity problem arising from the valuation of American options, which produces a nonlinear degenerated parabolic PDE with Black-Scholes operator. Based on the variational theory, the solvability and convergence properties of this penalty approach are established in a proper infinite dimensional space. Moreover, the convergence rate of the combination of two power penalty functions is obtained.     

On modified Mellin transforms, Gauss-Laguerre quadrature, and the valuation of American call options

We extend a framework based on Mellin transforms and show how to modify the approach to value American call options on dividend-paying stocks. We present a new integral equation to determine the price of an American call option and its free boundary using modified Mellin transforms. We also show how to derive the pricing formula for perpetual American call options using the new framework. A result due to Kim (1990) [24] regarding the optimal exercise price at expiry is also recovered. Finally, we apply Gauss-Laguerre quadrature for the purpose of an efficient and accurate numerical valuation.