Optimal stopping, free boundary, and American option in a jump-diffusion model

This paper considers the American put option valuation in a jump-diffusion model and relates this optimal-stopping problem to a parabolic integro-differential free-boundary problem, with special attention to the behavior of the optimal-stopping boundary. We study the regularity of the American option value and obtain in particular a decomposition of the American put option price as the sum of its counterpart European price and the early exercise premium. Compared with the Black-Scholes (BS) [5] model, this premium has an additional term due to the presence of jumps. We prove the continuity of the free boundary and also give one estimate near maturity, generalizing a recent result of Barleset al. [3] for the BS model. Finally, we study the effect of the market price of jump risk and the intensity of jumps on the American put option price and its critical stock price.     

Stochastic volatility asymptotics of stock loans: Valuation and optimal stopping

A stock loan, or equity security lending service, is a loan which uses stocks as collateral. The borrower has the right to repay the principal with interest and regain the stock, or make no repayment and surrender the stock. Therefore, the valuation of stock loan is an optimal stopping problem related to a perpetual American option with a negative effective interest rate. The negative effective interest rate makes standard techniques for perpetual American option pricing failure. Using a fast mean-reverting stochastic volatility model, we applied a perturbation technique to the free-boundary value problem for the stock loan price. An analytical pricing formula and optimal exercise boundary are derived by means of asymptotic expansion.     

Valuing American options under the CEV model by Laplace-Carson transforms

This paper investigates American option pricing under the constant elasticity of variance (CEV) model. Taking the Laplace-Carson transform (LCT) to the corresponding free-boundary problem enables the determination of the optimal early exercise boundary to be separated from the valuation procedure. Specifically, a functional equation for the LCT of the early exercise boundary is obtained. By applying Gaussian quadrature formulas, an efficient method is developed to compute the early exercise boundary, American option price and Greeks under the CEV model.     

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

An extension of the Euler Laplace transform inversion algorithm with applications in option pricing

We show that the Euler algorithm for Laplace transform inversion can be extended to functions defined on the entire real line, if they have specific decay features. Our objective is to apply the method to option pricing problems, specifically when inverting Laplace transforms of the option price in the logarithm of the strike.     

Pricing American options when asset prices jump

We present a transformation that helps price American options on assets that are modeled by a diffusion as well as a jump component. The presence of a jump component circumvents some shortcomings of the Black-Scholes model. The proposed transformation essentially transforms the arising free-boundary partial integro-differential equation (PIDE) into a sequence of fixed-boundary PIDEs which are much easier to solve. Finally, we provide numerical results illustrating convergence of the scheme and comparisons to other methods.     

基于LSMC模型的煤炭资源价值定价研究

煤炭资源价值定价可以抽象为一种美式期权定价问题.最小二乘蒙特卡洛模拟(LSMC)方法是解决美式期权定价问题的一个有效途径.详尽地分析了Cortazar等人的基于资源价格、利率和便利收益随机变动的三因素定价模型,利用向量Ito定理提出了三因素模型中价格、利率和便利收益变量的递推公式.对LSMC方法原理进行了细致的阐述,总结出实现LSMC方法的完整过程,并在Matlab环境下编制了LSMC算法实现程序,进行算例计算.算例结果表明,LSMC方法用于资源定价是有效可靠的.研究为煤炭资源价值定价提供了一个完整具有可操作性的工具.     

美式期权定价问题的数值方法

本文研究美式股票看跌期权定价问题的数值方法。通过将问题转化为等价的变分不等式方程，分别建立了半离散和全离散有限元逼近格式。并给出了有限元解的收敛性和稳定性分析。数值实验表明本文算法是一个高效和收敛的算法。     

Radial basis function partition of unity operator splitting method for pricing multi-asset American options

The operator splitting method in combination with finite differences has been shown to be an efficient approach for pricing American options numerically. Here, the operator splitting formulation is extended to the radial basis function partition of unity method. An approach that has previously often been used together with radial basis function methods to deal with the free boundary arising in American option pricing is to solve a penalised version of the Black–Scholes equation. It is shown that the operator splitting technique outperforms the penalty approach when used with the radial basis function partition of unity method. Numerical experiments are performed for one, two and three underlying assets. The advantage of the operator splitting technique grows with the number of dimensions.     

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

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

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.     

Compact finite difference method for American option pricing

A compact finite difference method is designed to obtain quick and accurate solutions to partial differential equation problems. The problem of pricing an American option can be cast as a partial differential equation. Using the compact finite difference method this problem can be recast as an ordinary differential equation initial value problem. The complicating factor for American options is the existence of an optimal exercise boundary which is jointly determined with the value of the option. In this article we develop three ways of combining compact finite difference methods for American option price on a single asset with methods for dealing with this optimal exercise boundary. Compact finite difference method one uses the implicit condition that solutions of the transformed partial differential equation be nonnegative to detect the optimal exercise value. This method is very fast and accurate even when the spatial step size h   is large (h?0.1)$(h?0.1)$. Compact difference method two must solve an algebraic nonlinear equation obtained by Pantazopoulos (1998) at every time step. This method can obtain second order accuracy for space x and requires a moderate amount of time comparable with that required by the Crank Nicolson projected successive over relaxation method. Compact finite difference method three refines the free boundary value by a method developed by Barone-Adesi and Lugano [The saga of the American put, 2003], and this method can obtain high accuracy for space x. The last two of these three methods are convergent, moreover all the three methods work for both short term and long term options. Through comparison with existing popular methods by numerical experiments, our work shows that compact finite difference methods provide an exciting new tool for American option pricing.     

Pricing perpetual American options under a stochastic-volatility model with fast mean reversion

In this paper, we present a “correction” to Merton’s (1973) well-known classical case of pricing perpetual American puts by considering the same pricing problem under a general fast mean-reverting SV (stochastic-volatility) model. By using the perturbation method, two analytic formulae are derived for the option price and the optimal exercise price, respectively. Based on the newly obtained formulae, we conduct a quantitative analysis of the impact of the SV term on the price of a perpetual American put option as well as its early exercise strategies. It shows that the presence of a fast mean-reverting SV tends to universally increase the put option price and to defer the optimal time to exercise the option contract, had the underlying been assumed to be falling. It is also noted that such an effect could be quite significant when the option is near the money.     

分数随机利率模型中的期权定价

本文研究分数随机利率模型中的期权定价问题.通过选取不同的资产作为计价单位及相应的测度交换,将经典模型中的测度变换方法推广到分数布朗运动市场环境,既丰富了分数期权定价的拟鞅方法,也得到了股票价格与利率分别服从几何分数布朗运动时的期权定价公式.     

American Call Options Under Jump‐Diffusion Processes – A Fourier Transform Approach

We consider the American option pricing problem in the case where the underlying asset follows a jump‐diffusion process. We apply the method of Jamshidian to transform the problem of solving a homogeneous integro‐partial differential equation (IPDE) on a region restricted by the early exercise (free) boundary to that of solving an inhomogeneous IPDE on an unrestricted region. We apply the Fourier transform technique to this inhomogeneous IPDE in the case of a call option on a dividend paying underlying to obtain the solution in the form of a pair of linked integral equations for the free boundary and the option price. We also derive new results concerning the limit for the free boundary at expiry. Finally, we present a numerical algorithm for the solution of the linked integral equation system for the American call price, its delta and the early exercise boundary. We use the numerical results to quantify the impact of jumps on American call prices and the early exercise boundary.     

Numerical valuation of options under Kou's model

Numerical methods are developed for pricing European and American options under Kou's jump-diffusion model which assumes the price of the underlying asset to behave like a geometrical Brownian motion with a drift and jumps whose size is log-double-exponentially distributed. The price of a European option is given by a partial integro-differential equation (PIDE) while American options lead to a linear complementarity problem (LCP) with the same operator. Spatial differential operators are discretized using finite differences on nonuniform grids and time stepping is performed using the implicit Rannacher scheme. For the evaluation of the integral term easy to implement recursion formulas are derived which have optimal computational cost. When pricing European options the resulting dense linear systems are solved using a stationary iteration. Also for pricing American options similar iterations can be employed. A numerical experiment demonstrates that the described method is very efficient as accurate option prices can be computed in a few milliseconds on a PC. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A new predictor-corrector scheme for valuing American puts

In this paper, we present a new numerical scheme, based on the finite difference method, to solve American put option pricing problems. Upon applying a Landau transform or the so-called front-fixing technique [19] to the Black-Scholes partial differential equation, a predictor-corrector finite difference scheme is proposed to numerically solve the nonlinear differential system. Through the comparison with Zhu’s analytical solution [35], we shall demonstrate that the numerical results obtained from the new scheme converge well to the exact optimal exercise boundary and option values. The results of our numerical examples suggest that this approach can be used as an accurate and efficient method even for pricing other types of financial derivative with American-style exercise.     

Equilibrium pricing bounds on option prices

We consider the problem of valuing European options in a complete market but with incomplete data. Typically, when the underlying asset dynamics is not specified, the martingale probability measure is unknown. Given a consensus on the actual distribution of the underlying price at maturity, we derive an upper bound on the call option price by putting two kinds of restrictions on the pricing probability measure. First, we put a restriction on the second risk-neutral moment of the underlying asset terminal value. Second, from equilibrium pricing arguments one can put a monotonicity restriction on the Radon-Nikodym density of the pricing probability with respect to the true probability measure. This density is restricted to be a nonincreasing function of the underlying price at maturity. The bound appears then as the solution of a constrained optimization problem and we adopt a duality approach to solve it. Explicit bounds are provided for the call option. Finally, we provide a numerical example.     

美式看跌期权定价的一种混合数值方法

本文对美式看跌期权的定价提供了一种新的混合数值方法,即快速傅里叶变换法加龙格-库塔法.首先将美式看跌期权价格所满足的Black-Scholes微分方程定解问题转化为一个标准的抛物型初、边值问题,然后通过傅里叶变换,使之转换为一个不带股价变量的常微分方程初值问题,再利用龙格-库塔法对其进行数值求解.数值实验表明,本文算法是一种快速的高精度的算法.     

随机利率下期权定价

本文在随机利率的基础上,考虑股票价格过程和利率过程分别为扩散过程和Ito过程,并且在相关的假设下,运用鞅方法推导出欧式期权价值过程所满足的微分方程;以及利率满足一种特殊方程时,运用最优停止的鞅方法,得到了随机利率下美式期权的价格和最优停时.