自适应隐式有限差分方法求解美式期权定价问题

本文针对美式期权的定价问题设计了基于有限差分方法的预估-校正数值算法.该算法采用显式离散格式先对自由边界条件进行预估,再对经过变量替换后的关于期权价格的偏微分方程采用隐式格式离散,并用Fourier方法分析了此离散格式的稳定性.接下来,引入基于Richardson外推法的后验误差指示子.这个后验误差指示子能够在给定的误差阈值范围内,针对期权价格和自由边界找到合适的网格划分.最后,通过设计多组数值实验并与Fazio^[1]采用显式离散格式算得的数值结果相比较,验证了所提算法的有效性,稳定性和收敛性.     

源于俄式期权定价的自由边界问题

本文主要应用PDE方法对俄式期权定价问题进行理论分析. 类似于美式期权定价问题,俄罗斯期权定价问题可归结为-个-维抛物型变分不等式.我们首先引入惩罚函数证明了该变分不等式的解的存在唯-性,然后研究了自由边界的一些性质,如单调性、光滑性和自由边界的位置.     

连续支付红利的美式障碍期权定价的数值方法

本文利用显式差分格式为连续支付红利的向下触销型美式障碍期权定价.由于障碍的影响,定价模型的边值条件含有间断,故把结点设在障碍水平上,并在障碍附近的区域内运用局部网格加密技术,这样就可以得到较精确的期权价格.本文给出数值算例,验证算法的有效性,并分析障碍对期权价格的影响.     

美式利率期权定价的抛物型变分不等式

应用PDE方法对美式利率期权定价问题进行理论分析.在CIR利率模型下美式利率期权定价问题可归结为一个退化的一维抛物型变分不等式.通过引入惩罚函数证明了该变分不等式的解的存在唯一性,然后研究了自由边界的一些性质,如单调性,光滑性和自由边界在终止期的位置.     

投影三角分解法定价带随机波动率的美式期权

考虑数值求解Heston随机波动率美式期权定价问题,通过在空间方向采用中心差分格式离散二维偏微分算子,在时间方向利用隐式交替方向格式,将美式期权定价问题转化成求解每个时间层上的若干个线性互补问题.针对一般美式期权定价模型离散得到的线性互补问题,构造出投影三角分解法进行求解,并在理论上给出算法的收敛条件.数值实验表明,所构造的数值方法对于求解美式期权定价问题是有效的,并且优于经典的投影超松弛迭代法和算子分裂方法.     

CEV过程下回望期权定价的高精度收敛差分算法

主要研究了CEV过程下一类回望期权的定价的数值解法问题.首先对期权价格所满足的微分方程中的空间变量进行半离散化处理,得到了具体的半离散化差分格式,然后证明了该差分格式具有稳定性和收敛性.数值试验表明本文算法是一个稳定收敛的算法.     

美式期权自由边界的计算方法

美式期权的自由边界问题在金融工程文献中已经引起了广泛的关注,然而,它的数值计算方法一直是一个难点.基于差分技巧,给出了满足具有有限到期日的美式期权自由边界的两种计算方法,即,根据股票期权价格和相应的偏导数来确定自由边界条件.数值结果表明了上述两种方法下自由边界是一致性的.此外研究结果对自由边界的计算提供很好的科学依据.     

基于不光滑边界的变系数抛物型方程的高精度紧格式

本文讨论基于不光滑边界的变系数抛物型方程求解的高精度紧格式.首先构造一般变系数抛物型方程的高精度紧格式,并在理论上证明格式具有空间方向四阶精度.然后针对非光滑边界条件,引入局部网格加密技巧在奇异点附近进行不均匀的网格加密.数值实验以期权定价中Black-Scholes偏微分方程的求解为例,验证高精度紧格式用于光滑边界条件的微分方程离散可以达到四阶精度.对于处理非光滑边界条件,网格局部加密技巧能有效的提高数值解精度,使得高精度紧格式用于定价欧式期权可以接近四阶精度.     

CEV跳-扩散模型中期权定价研究

考虑了CEV与Kou双指数跳-扩散组合模型中的期权定价问题.首先,运用Ito公式和期权定价的无套利原理,得到了模型下期权价格所满足的偏积-微分方程.然后,运用中心差分和Lagrange线性插值,分别对偏积-微分方程中的微分项和积分项进行离散化处理,再由Euler法,最终得了偏积-微分方程的有限差分格式,并且对差分方法的误差和收敛性进行了分析.最后数值实验验证了该算法是一个稳定且收敛的算法.     

一种亚式风格可重置执行价格期权设计

本文设计了一种亚式风格的可重置执行价格期权;严格证明了可重置执行边界的存在性,以及连续区域与重置区域的单连通性;利用Hartman-Watson分布,写出了可重置期权的定价公式,并利用此公式给出了可重置执行边界的一种新的数值算法.     

Laplace transforms and American options

Laplace transform methods are used to study the valuation of American call and put options with constant dividend yield, and to derive integral equations giving the location of the optimal exercise boundary. In each case studied, the main result of this paper is a nonlinear Fredholm-type integral equation for the location of the free boundary. The equations differ depending on whether the dividend yield is less than or exceeds the risk-free rate. These integral equations contain a transform variable, so the solution of the equations would involve finding the free boundary that satisfies the equations for all values of this transform variable. Expressions are also given for the transform of the value of the option in terms of this free boundary.     

Differential quadrature method for pricing American options

In this article, differential quadrature method (DQM), a highly accurate and efficient numerical method for solving nonlinear problems, is used to overcome the difficulty in determining the optimal exercise boundary of American option. The following three parts of the problem in pricing American options are solved. The first part is how to treat the uncertainty of the early exercise boundary, or free boundary in the language of the PDE treatment of the American option, because American options can be exercised before the date of expiration. The second part is how to solve the nonlinear problem, because the problem of pricing American options is nonlinear. And the third part is how to treat the initial value condition with the singularity and the boundary conditions in the DQM. Numerical results for the free boundary of American option obtained by both DQM and finite difference method (FDM) are given and from which it can be seen the computational efficiency is greatly improved by DQM. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 711–725, 2002; Published online in Wiley InterScience (www.interscience.wiley.com); DOI 10.1002/num.10028.     

Penalty methods for the numerical solution of American multi-asset option problems

We derive and analyze a penalty method for solving American multi-asset option problems. A small, non-linear penalty term is added to the Black–Scholes equation. This approach gives a fixed solution domain, removing the free and moving boundary imposed by the early exercise feature of the contract. Explicit, implicit and semi-implicit finite difference schemes are derived, and in the case of independent assets, we prove that the approximate option prices satisfy some basic properties of the American option problem. Several numerical experiments are carried out in order to investigate the performance of the schemes. We give examples indicating that our results are sharp. Finally, the experiments indicate that in the case of correlated underlying assets, the same properties are valid as in the independent case.     

A parabolic variational inequality arising from the valuation of strike reset options

A strike reset option is an option that allows its holder to reset the strike price to the prevailing underlying asset price at a moment chosen by the holder. The pricing model of the option can be formulated as a one-dimensional parabolic variational inequality, or equivalently, a free boundary problem, where the free boundary just corresponds to the optimal reset strategy adopted by the holder of the option. This paper is concerned with the theoretical analysis of the model. The existence and uniqueness of the solution are established. Furthermore, we study properties of the free boundary. The monotonicity and C^∞ smoothness of the free boundary are proven in some situations.     

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.     

Approximations for the values of american options

The solution of the American option valuation problem is the solution of a parabolic partial differential equation satisfying free boundary conditions. The free boundary represents the critical price, at which the option should be exercised. In this paper the free boundary is determined by an algebraic relation and an approximate solution derived. A suitable modification of the approximate solution gives the exact solution. The uniqueness of the free boundary implies the expression determined by the algebraic relation is the true critical price     

Nonlocal boundary value problems in differential and difference settings

We study nonlocal boundary value problems of the first and second kind for the heat equation with variable coefficients in the differential and difference settings. By the method of energy inequalities, we find a priori estimates for the differential and difference problems.     

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.     

The regularity of the free boundary for strike reset option

A strike reset option is an option that allows its holder to reset the strike price to the prevailing underlying asset price at a moment chosen by the holder. The pricing model of the option can be formulated as a parabolic variational inequality and the optimal reset strategy is the free boundary. The smoothness of the free boundary in some cases was showed in our article published in JDE. We would prove its smoothness in the other case in this paper by a generalized comparison principle for the variational inequality.     

THE REGULARITY OF THE FREE BOUNDARY FOR STRIKE RESET OPTION

A strike reset option is an option that allows its holder to reset the strike price to the prevailing underlying asset price at a moment chosen by the holder. The pricing model of the option can be formulated as a parabolic variational inequality and the optimal reset strategy is the free boundary. The smoothness of the free boundary in some cases was showed in our article published in JDE. We would prove its smoothness in the other case in this paper by a generalized comparison principle for the variational inequality.