[论文解读] High-order accurate implicit methods for the pricing of barrier options

J. C. Ndogmo, Davis Bundi Ntwiga|ArXiv.org|Sep 29, 2007

Stochastic processes and financial applications参考文献 16被引用 8

一句话总结

本文提出了一种用于定价障碍期权的高阶精度隐式有限差分格式，结合了四阶空间精度与基于概率的最优边界条件确定方法。该方法在速度和精度方面表现卓越——尤其适用于离散监控的障碍期权——其计算结果最接近半解析基准值，大多数情况下计算时间低于1.2秒，即使在最复杂的配置下也低于5秒。

ABSTRACT

This paper deals with a high-order accurate implicit finite-difference approach to the pricing of barrier options. In this way various types of barrier options are priced, including barrier options paying rebates, and options on dividend-paying-stocks. Moreover, the barriers may be monitored either continuously or discretely. In addition to the high-order accuracy of the scheme, and the stretching effect of the coordinate transformation, the main feature of this approach lies on a probability-based optimal determination of boundary conditions. This leads to much faster and accurate results when compared with similar pricing approaches. The strength of the present scheme is particularly demonstrated in the valuation of discretely monitored barrier options where it yields values closest to those obtained from the only semi-analytical valuation method available.

研究动机与目标

• 开发一种高阶精度的数值方法，用于定价各类障碍期权，包括含回收费用、股息及双障碍期权。
• 解决现有有限差分格式在计算效率与精度方面的不足，特别是针对离散监控障碍期权的场景。
• 通过基于概率的最优边界条件确定方法改进边界条件处理，提升精度与收敛性。
• 在连续与离散监控下，通过实证证明该方法在速度与精度方面优于标准格式（如Crank-Nicolson、显式、全隐式格式）。
• 通过与唯一可用的离散监控障碍期权半解析解进行对比，验证方法的有效性，结果显示结果近乎完全一致。

提出的方法

• 采用θ-格式，通过将θ与网格尺寸关联，实现空间四阶精度与时间二阶精度。
• 应用坐标变换，将资产价格域映射到计算域，通过网格拉伸提升关键区域（如障碍区域）的分辨率。
• 边界条件采用基于概率的最优确定方法，确保边界处期权价值精确，显著提升精度。
• 在变换后的x-τ坐标系中采用非均匀网格实现该格式，但x方向的均匀网格取得了最佳结果。
• 该方法扩展用于处理多种障碍类型：单重与双重敲出、敲入期权、支付股息的股票，以及回收费用。
• 通过连续与离散监控的数值实验对该方法进行验证，网格尺寸最高达4000×6000，以实现高精度。

实验结果

研究问题

• RQ1 高阶精度隐式有限差分格式是否能在障碍期权定价中优于标准格式（如Crank-Nicolson）？
• RQ2 基于概率的最优边界条件确定方法如何提升障碍期权估值的精度，特别是在离散监控情形下？
• RQ3 所提方法在离散监控障碍期权上能达到多接近半解析解的精度？此类解析解极为稀少。
• RQ4 坐标变换与高阶空间离散化是否能显著降低计算成本，同时保持高精度？
• RQ5 该方法能否有效应用于具有多种特征的复杂障碍期权，如回收费用、股息与双障碍？

主要发现

• 该方法在所有连续监控情形下计算时间均低于1.2秒，在最复杂的离散监控期权（如q=0.045的双敲出期权）下也低于5秒，其中最长达4.20秒。
• 对于离散监控障碍期权，该方法的计算结果最接近唯一的半解析基准解，所有测试案例的绝对差异均小于0.01。
• 即使仅使用20×20的最小网格，该方法在连续监控下仍表现出极高的精度，展现出卓越的收敛性与稳定性。
• 基于概率的边界条件确定方法显著提升了精度，尤其在单障碍与离散监控双障碍期权中表现突出。
• 在变换坐标系中使用非均匀网格并未提升性能；x坐标方向的均匀网格表现最佳，表明该格式在变换空间中保持网格均匀性为最优策略。
• 该方法在速度与精度方面均优于标准格式（Crank-Nicolson、显式、全隐式），尤其在具有挑战性的离散监控情形下表现更优。