美式期权定价问题的变网格差分方法

本文提出一种求解美式期权定价自由边值问题的变网格差分方法.通过建立一个自由边界所满足的方程,利用变网格技术可同时求出期权的差分解和最佳执行边界.本文分别讨论了显式和隐式变网格差分格式,并给出了差分解的收敛性和稳定性分析.数值实验表明本文算法是一个非常有效的期权定价算法.     

美式利率期权的最佳实施边界的分析

在Vasicek利率模型的假设下,应用变分不等式方法分析了美式利率期权自由边界的性质．首先我们得到美式利率期权自由边界的下界, 然后把自由边界问题化为变分不等式,通过引入惩罚函数证明了该变分不等式解的存在唯一性,最后证明了自由边界的单调性、 有界性和C∞光滑性．     

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

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

双币种博弈期权

博弈期权是一种赋予期权出售方在期权有效期内任意时刻可以赎回合约权利的美式期权.在B-S框架下分析了双币种情形下的博弈期权定价行为,建立了双币种博弈期权的定价模型,分别讨论了敲定价以国内货币计价和国外货币计价下的博弈期权定价问题及其最优赎回策略,通过运用偏微分方程的方法得到了这两种情形下期权价格的表达式及其最优执行边界.最后通过数值模拟,分析了标的资产和汇率的波动水平以及汇率与标的资产的相关系数对期权的最优执行策略和违约金边界的影响.     

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

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

美式期权定价的数值方法

介绍了定价美式期权的几种常见数值方法.对最近几年的主要研究成果做了简单的介绍和比较,并给出了数值算例.特别回顾了美式期权定价的蒙特卡罗模拟加速方法.     

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

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

永久经理期权的最优实施策略

利用偏微分方程数值方法研究金融市场上永久经理期权(ESOs)的最优实施策略问题.讨论了两种实施情况,即通常的整体实施情况以及非限制实施情况.在非限制实施情况下,持有者在任意可实施时刻可以实施其持有的任意份ESOs.两种实施情况下的最优实施策略分别对应着一个抛物型变分不等式定解问题的自由边界(最佳实施边界).通过数值分析的方法分别研究了自由边界的性质,比较了两种情况下自由边界的异同及其所对应的金融意义.     

随机波动率与双指数跳扩散组合模型的美式期权定价

在股价满足Cox-Ingersoll-Ross(CIR)随机波动率与Kou的双指数跳扩散组合模型下,利用随机分析方法讨论了美式看跌期权函数及最佳实施边界的性质.应用一阶线性近似实施边界获得了期权价格的拟解析式和实施边界满足的非线性方程.进一步,应用梯形法离散处理方程式内积分表达式,建立了期权最佳实施边界和价格的数值算法.最后分别给出了常数波动率或CIR随机波动率的数值实例.     

分数Black-Scholes模型下美式期权定价的积分方程式

为得到分数Black-Scholes模型下美式期权价格的公式,文章以看涨期权为例,应用偏微分方程法,推导期权价格的积分方程式.由于美式期权的价格可分解为欧式期权的价格和由于提前实施需要增付的期权金,而提前实施期权金与最佳实施边界的位置有关,所以为导出最佳实施边界所满足的方程,文章首先研究分数Black-Scholes方程的基本解,然后建立美式看涨期权的分解公式,推导最佳实施边界适合的非线性积分方程,从而得到美式看涨期权价格的积分方程式.美式看跌期权价格的积分方程式类似得到.     

Exercisability Randomization of the American Option

Abstract

The valuation of American options is an optimal stopping time problem which typically leads to a free boundary problem. We introduce here the randomization of the exercisability of the option. This method considerably simplifies the problematic by transforming the free boundary problem into an evolution equation. This evolution equation can be transformed in a way that decomposes the value of the randomized option into a European option and the present value of continuously paid benefits. This yields a new binomial approximation for American options. We prove that the method is accurate and numerical results illustrate that it is computationally efficient.     

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.     

A NUMERICAL METHOD FOR DETERMINING THE OPTIMAL EXERCISE PRICE TO AMERICAN OPTIONS

American Options can be exercised prior to the date of expiration,the valuation of American options then constitutes a free boundary value problem.How to determine the free boundary,i.e. the optimal exercise price,is a key problem.In this paper,a nonlinear equation is given.The free boundary can be obtained by solving the nonlinear equation and the numerical results are better.     

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.     

On the asymptotic free boundary for the American put option problem

In practical work with American put options, it is important to be able to know when to exercise the option, and when not to do so. In computer simulation based on the standard theory of geometric Brownian motion for simulating stock price movements, this problem is fairly easy to handle for options with a short lifespan, by analyzing binomial trees. It is considerably more challenging to make the decision for American put options with long lifespan. In order to provide a satisfactory analysis, we look at the corresponding free boundary problem, and show that the free boundary—which is the curve that separates the two decisions, to exercise or not to—has an asymptotic expansion, where the coefficient of the main term is expressed as an integral in terms of the free boundary. This raises the perspective that one could use numerical simulation to approximate the integral and thus get an effective way to make correct decisions for long life options.     

Front-fixing FEMs for the pricing of American options based on a PML technique

In this paper, efficient numerical methods are developed for the pricing of American options governed by the Black–Scholes equation. The front-fixing technique is first employed to transform the free boundary of optimal exercise prices to some a priori known temporal line for a one-dimensional parabolic problem via the change of variables. The perfectly matched layer (PML) technique is then applied to the pricing problem for the effective truncation of the semi-infinite domain. Finite element methods using the respective continuous and discontinuous Galerkin discretization are proposed for the resulting truncated PML problems related to the options and Greeks. The free boundary is determined by Newton’s method coupled with the discrete truncated PML problem. Stability and nonnegativeness are established for the approximate solution to the truncated PML problem. Under some weak assumptions on the PML medium parameters, it is also proved that the solution of the truncated PML problem converges to that of the unbounded Black–Scholes equation in the computational domain and decays exponentially in the perfectly matched layer. Numerical experiments are conducted to test the performance of the proposed methods and to compare them with some existing methods.     

Weak Galerkin finite element method for valuation of American options

We introduce a weak Galerkin finite element method for the valuation of American options governed by the Black-Scholes equation. In order to implement, we need to solve the optimal exercise boundary and then introduce an artificial boundary to make the computational domain bounded. For the optimal exercise boundary, which satisfies a nonlinear Volterra integral equation, it is resolved by a higher-order collocation method based on graded meshes. With the computed optimal exercise boundary, the front-fixing technique is employed to transform the free boundary problem to a one- dimensional parabolic problem in a half infinite area. For the other spatial domain boundary, a perfectly matched layer is used to truncate the unbounded domain and carry out the computation. Finally, the resulting initial-boundary value problems are solved by weak Galerkin finite element method, and numerical examples are provided to illustrate the efficiency of the method.     

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

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

Spectral Method for the Black-Scholes Model of American Options Valuation

In this paper, we devote ourselves to the research of numerical methods for American option pricing problems under the Black-Scholes model. The optimal exercise boundary which satisfies a nonlinear Volterra integral equation is resolved by a high-order collocation method based on graded meshes. For the other spatial domain boundary, an artificial boundary condition is applied to the pricing problem for the effective truncation of the semi-infinite domain. Then, the front-fixing and stretching transformations are employed to change the truncated problem in an irregular domain into a one-dimensional parabolic problem in [−1,1]. The Chebyshev spectral method coupled with fourth-order Runge-Kutta method is proposed for the resulting parabolic problem related to the options. The stability of the semi-discrete numerical method is established for the parabolic problem transformed from the original model. Numerical experiments are conducted to verify the performance of the proposed methods and compare them with some existing methods.     

Solving complex PDE systems for pricing American options with regime‐switching by efficient exponential time differencing schemes

In this article, we study the numerical solutions of a class of complex partial differential equation (PDE) systems with free boundary conditions. This problem arises naturally in pricing American options with regime‐switching, which adds significant complexity in the PDE systems due to regime coupling. Developing efficient numerical schemes will have important applications in computational finance. We propose a new method to solve the PDE systems by using a penalty method approach and an exponential time differencing scheme. First, the penalty method approach is applied to convert the free boundary value PDE system to a system of PDEs over a fixed rectangular region for the time and spatial variables. Then, a new exponential time differncing Crank–Nicolson (ETD‐CN) method is used to solve the resulting PDE system. This ETD‐CN scheme is shown to be second order convergent. We establish an upper bound condition for the time step size and prove that this ETD‐CN scheme satisfies a discrete version of the positivity constraint for American option values. The ETD‐CN scheme is compared numerically with a linearly implicit penalty method scheme and with a tree method. Numerical results are reported to illustrate the convergence of the new scheme. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013