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 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…     

A Modified Binomial Tree Method for Currency Lookback Options

Abstract The binomial tree method is the most popular numerical approach to pricing options. However, for currency lookback options, this method is not consistent with the corresponding continuous models, which leads to slow speed of convergence. On the basis of the PDE approach, we develop a consistent numerical scheme called the modified binomial tree method. It possesses one order of accuracy and its efficiency is demonstrated by numerical experiments. The convergence proofs are also produced in terms of numerical analysis and the notion of viscosity solution. Supported by National Science Foundation of China (No. 19871062)     

Pricing of perpetual American and Bermudan options by binomial tree method

In this paper, we consider the binomial tree method for pricing perpetual American and perpetual Bermudan options. The closed form solutions of these discrete models are solved. Explicit formulas for the optimal exercise boundary of the perpetual American option is obtained. A nonlinear equation that is satisfied by the optimal exercise boundaries of the perpetual Bermudan option is found.      

求解美式回望期权的有限元方法

期权作为一种金融衍生产品,在欧美国家一直很受欢迎.由于其规避风险的特性,期权也吸引了中国投资者的兴趣.基于市场的需求,2015年初,上海证券交易所推出了中国首批期权产品,期权定价问题的研究热潮正席卷全球.本文研究的美式回望期权,是一种路径相关的期权,其支付函数不仅依赖于标的资产的现值,也依赖其历史最值.分析回望期权的特点,不难发现:1)这类期权空间变量的变化范围为二维无界不规则区域,难以应用数值方法直接求解;2)最佳实施边界未知,使得该问题变得高度非线性.本文的主要工作就是解决这两个困难,得到回望期权和最佳实施边界的数值逼近结果.现有的处理问题1)的有效方法是采用标准变量替换、计价单位变换以及Landau变换将定价模型化为一个[0,1]区间上的非线性抛物问题,本文也将沿用这些技巧处理问题1).进一步,采用有限元方法离散简化后的定价模型,并论证了数值解的非负性,提出了利用Newton法求解离散化的非线性系统.最后,通过数值模拟,验证了本文所提算法的高效性和准确性.     

CEV下有交易费用的回望期权的定价研究

本文在研究服从CEV过程且无交易费用的回望期权定价模型的基础上,推导出CEV下有交易费用的回望期权定价模型,并利用变量转换和二叉树方法求解,最终给出了CEV下有交易费用的回望期权的近似解。     

Pricing and Hedging of Lookback Options in Hyper-exponential Jump Diffusion Models

ABSTRACT

In this article, we consider the problem of pricing lookback options in certain exponential Lévy market models. While in the classic Black-Scholes models the price of such options can be calculated in closed form, for more general asset price model, one typically has to rely on (rather time-intense) Monte-Carlo or partial (integro)-differential equation (P(I)DE) methods. However, for Lévy processes with double exponentially distributed jumps, the lookback option price can be expressed as one-dimensional Laplace transform (cf. Kou, S. G., Petrella, G., & Wang, H. (2005). Pricing path-dependent options with jump risk via Laplace transforms. The Kyoto Economic Review, 74(9), 1–23.). The key ingredient to derive this representation is the explicit availability of the first passage time distribution for this particular Lévy process, which is well-known also for the more general class of hyper-exponential jump diffusions (HEJDs). In fact, Jeannin and Pistorius (Jeannin, M., & Pistorius, M. (2010). A transform approach to calculate prices and Greeks of barrier options driven by a class of Lévy processes. Quntitative Finance, 10(6), 629–644.) were able to derive formulae for the Laplace transformed price of certain barrier options in market models described by HEJD processes. Here, we similarly derive the Laplace transforms of floating and fixed strike lookback option prices and propose a numerical inversion scheme, which allows, like Fourier inversion methods for European vanilla options, the calculation of lookback options with different strikes in one shot. Additionally, we give semi-analytical formulae for several Greeks of the option price and discuss a method of extending the proposed method to generalized hyper-exponential (as e.g. NIG or CGMY) models by fitting a suitable HEJD process. Finally, we illustrate the theoretical findings by some numerical experiments.     

A fast high-order finite difference algorithm for pricing American options

We describe an improvement of Han and Wu’s algorithm [H. Han, X.Wu, A fast numerical method for the Black–Scholes equation of American options, SIAM J. Numer. Anal. 41 (6) (2003) 2081–2095] for American options. A high-order optimal compact scheme is used to discretise the transformed Black–Scholes PDE under a singularity separating framework. A more accurate free boundary location based on the smooth pasting condition and the use of a non-uniform grid with a modified tridiagonal solver lead to an efficient implementation of the free boundary value problem. Extensive numerical experiments show that the new finite difference algorithm converges rapidly and numerical solutions with good accuracy are obtained. Comparisons with some recently proposed methods for the American options problem are carried out to show the advantage of our numerical method.     

Uniqueness of positive radial solutions for quasilinear elliptic equations in an annulus

Extending a previous result of Tang [1] we prove the uniqueness of positive radial solutions of Δ_pu+f(u)=0, subject to Dirichlet boundary conditions on an annulus in Rⁿ with 2<p≤n, under suitable hypotheses on the nonlinearity f. This argument also provides an alternative proof for the uniqueness of positive solutions of the same problem in a finite ball (see [9]), in the complement of a ball or in the whole space Rⁿ (see [10], [3] and [11]).     

基于跳扩散模型欧式期权定价的条件二叉树方法

对股票价格的跳扩散模型进行了分析,在CRR二叉树期权定价模型的基础上考虑标的股票价格发生跳跃的情况,得出基于跳扩散过程的股票期权的条件二叉树定价模型,并且证明在极限情况下,该条件二叉树模型的期权定价公式趋于Merton的解析定价公式,数值试验证实该条件二叉树模型的有效性。     

Pricing American Asian options with higher moments in the underlying distribution

We develop a modified Edgeworth binomial model with higher moment consideration for pricing American Asian options. With lognormal underlying distribution for benchmark comparison, our algorithm is as precise as that of Chalasani et al. [P. Chalasani, S. Jha, F. Egriboyun, A. Varikooty, A refined binomial lattice for pricing American Asian options, Rev. Derivatives Res. 3 (1) (1999) 85–105] if the number of the time steps increases. If the underlying distribution displays negative skewness and leptokurtosis as often observed for stock index returns, our estimates can work better than those in Chalasani et al. [P. Chalasani, S. Jha, F. Egriboyun, A. Varikooty, A refined binomial lattice for pricing American Asian options, Rev. Derivatives Res. 3 (1) (1999) 85–105] and are very similar to the benchmarks in Hull and White [J. Hull, A. White, Efficient procedures for valuing European and American path-dependent options, J. Derivatives 1 (Fall) (1993) 21–31]. The numerical analysis shows that our modified Edgeworth binomial model can value American Asian options with greater accuracy and speed given higher moments in their underlying distribution.     

The Calderón-Zygmund property for quasilinear divergence form equations over Reifenberg flat domains

The results by Palagachev (2009) [3] regarding global Hölder continuity for the weak solutions to quasilinear divergence form elliptic equations are generalized to the case of nonlinear terms with optimal growths with respect to the unknown function and its gradient. Moreover, the principal coefficients are discontinuous with discontinuity measured in terms of small BMO norms and the underlying domain is supposed to have fractal boundary satisfying a condition of Reifenberg flatness. The results are extended to the case of parabolic operators as well.     

On the global regularity of generalized Leray-alpha type models

We generalize Leray-alpha type models studied in Cheskidov et al. (2005) [1] and Linshiz and Titi (2007) [4] via fractional Laplacians and employ Besov space techniques to obtain global regularity results with the logarithmically supercritical dissipation.     

Lattice Boltzmann methods for solving partial differential equations of exotic option pricing

This paper establishes a lattice Boltzmann method (LBM) with two amending functions for solving partial differential equations (PDEs) arising in Asian and lookback options pricing. The time evolution of stock prices can be regarded as the movement of randomizing particles in different directions, and the discrete scheme of LBM can be interpreted as the binomial models. With the Chapman-Enskog multi-scale expansion, the PDEs are recovered correctly from the continuous Boltzmann equation and the computational complexity is O(N), where N is the number of space nodes. Compared to the traditional LBM, the coefficients of equilibrium distribution and amending functions are taken as polynomials instead of constants. The stability of LBM is studied via numerical examples and numerical comparisons show that the LBM is as accurate as the existing numerical methods for pricing the exotic options and takes much less CPU time.     

跳扩散模型下美式期权的对称性

利用分析方法得到了跳扩散模型下美式看涨、看跌期权的价格和最佳实施边界间的对称性公式．美式看涨和看跌期权价格问的对称关系通常是利用概率理论得到,这里给出了这些结果在跳扩散模型下的另一种证明．此外,由本文所得结果和偏微分方程理论,可以得到跳扩散模型下美式看涨期权的最佳实施边界以及永久美式期权的若干性质．     

A renewal jump-diffusion process with threshold dividend strategy

In this paper, we consider a jump-diffusion risk process with the threshold dividend strategy. Both the distributions of the inter-arrival times and the claims are assumed to be in the class of phase-type distributions. The expected discounted dividend function and the Laplace transform of the ruin time are discussed. Motivated by Asmussen [S. Asmussen, Stationary distributions for fluid flow models with or without Brownian noise, Stochastic Models 11 (1) (1995) 21–49], instead of studying the original process, we study the constructed fluid flow process and their closed-form formulas are obtained in terms of matrix expression. Finally, numerical results are provided to illustrate the computation.     

Supremum principles for the higher-order derivatives of solutions to the Dirichlet problem for parabolic equations without compatibility conditions between the data

The Dirichlet problem is considered for the heat equation ut=auxx, a>0 a constant, for (x,t)∈[0,1]×[0,T], without assuming any compatibility condition between initial and boundary data at the corner points (0,0) and (1,0). Under some smoothness restrictions on the data (stricter than those required by the classical maximum principle), weak and strong supremum and infimum principles are established for the higher-order derivatives, ut and uxx, of the bounded classical solutions. When compatibility conditions of zero order are satisfied (i.e., initial and boundary data coincide at the corner points), these principles allow to estimate the higher-order derivatives of classical solutions uniformly from below and above on the entire domain, except that at the two corner points. When compatibility conditions of the second order are satisfied (i.e., classical solutions belong to on the closed domain), the results of the paper are a direct consequence of the classical maximum and minimum principles applied to the higher-order derivatives. The classical principles for the solutions to the Dirichlet problem with compatibility conditions are generalized to the case of the same problem without any compatibility condition. The Dirichlet problem without compatibility conditions is then considered for general linear one-dimensional parabolic equations. The previous results as well as some new properties of the corresponding Green functions derived here allow to establish uniformL¹-estimates for the higher-order derivatives of the bounded classical solutions to the general problem.     

Pricing catastrophe options in discrete operational time

We employ a doubly-binomial process as in Gerber [Gerber, H.U., 1988. Mathematical fun with the compound binomial process. ASTIN Bull. 18, 161-168] to discretize and generalize the continuous “randomized operational time” model of Chang et al. ([Chang, C.W., Chang, J.S.K., Yu, M.T., 1996. Pricing catastrophe insurance futures call spreads: A randomized operational time approach. J. Risk Insurance 63, 599-616] and CCY hereafter) from a complete-market continuous-time setting to an incomplete-market discrete-time setting, so as to price a richer set of catastrophe (CAT) options. For futures options, we derive the equivalent martingale probability measures by benchmarking to the shadow price of a bond to span arrival uncertainty, and the underlying futures price to span price uncertainty. With a time change from calendar time to the operational transaction-time dimension, we derive CCY as a limiting case under risk-neutrality when both calendar-time and transaction-time intervals shrink to zero. For a cash option with non-traded underlying loss index, we benchmark to the market reinsurance premiums to span claim uncertainty, and with a time change to claim time, we derive the cash option price as a binomial sum of claim-time binomial Asian option prices under the martingale measures.     

The term structure of interest rates in the economic and monetary union

The main purpose of this article is to present a new numerical procedure that can be used to implement a variety of different interest rate models. The new approach allows to construct no-arbitrage models for the term structure, where the stochastic process driving the rates is infinitely divisible, as in the cases of pure-diffusion and jump-diffusion mean reverting models. The new method determines a unique fully specified hexanomial tree, consistent with risk neutral probabilities. A simple forward recursive procedure solves for the entire tree. The proposed lattice model, which generalized the Hull and White [37] single-factor model, is relatively simple, computational efficient and can fit any initial term structure observed in the market. Numerical experiments demonstrate how the jump-diffusion mean reverting model is particularly suited to describe the European money market rates behavior. Interest rates controlled by the monetary authorities behave as if they are jump processes and the term structure, at short maturity, is contingent upon the levels of these official rates.     

An analytic pricing formula for lookback options under stochastic volatility

In this work, an analytic pricing formula for floating strike lookback options under Heston’s stochastic volatility model is derived by means of the homotopy analysis method. The fixed strike lookback options can then be priced on the basis of the results of floating strike and the put–call parity relation for lookback options.