第08章期权定价的数值方法

1、第八章 期权定价的数值方法在前面几章中，我们得到了期权价值所满足的偏微分方程，并且解出了一些精确的期权解析定价公式。但是在很多情形中，我们无法得到期权价值的解析解，这时人们经常采用数值方法（Numerical Procedures）为期权定价，其中包括二叉树方法（Binomial Trees）、蒙特卡罗模拟（Monte Carlo Simulation）和有限差分方法（Finite Difference Methods）。当期权收益依赖于标的变量所遵循的历史路径时（如我们将在第九章看到的路径依赖期权），或是期权价值取决于多个标的变量的时候，可以用蒙特卡罗模拟为期权定价。而二叉树图和有限差分方法

2、则比较适用于有提前执行可能性的期权。在这一章里，我们将介绍如何借助上述三种数值方法来为期权定价。为了便于表达，本章中统一假设当前时刻为零时刻，表示为0。第一节 二叉树期权定价模型二叉树期权定价模型是由J. C. Cox、S. A. Ross和M. Rubinstein于1979年首先提出的，已经成为金融界最基本的期权定价方法之一。二叉树模型的优点在于其比较简单直观，不需要太多的数学知识就可以加以应用。一、二叉树模型的基本方法我们从简单的无收益资产期权的定价开始讲解二叉树模型，之后再逐步加以扩展。二叉树模型首先把期权的有效期分为很多很小的时间间隔，并假设在每一个时间间隔内证券价格只有两种运动的可

3、能：从开始的上升到原先的倍，即到达；下降到原先的倍，即。其中，如图8.1所示。价格上升的概率假设为，下降的概率假设为。图8.1 时间内资产价格的变动相应地，期权价值也会有所不同，分别为和。注意，在较大的时间间隔内，这种二值运动的假设当然不符合实际，但是当时间间隔非常小的时候，比如在每个瞬间，资产价格只有这两个运动方向的假设是可以接受的。因此，二叉树模型实际上是在用大量离散的小幅度二值运动来模拟连续的资产价格运动。（一）单步二叉树模型运用单步二叉树为期权定价，可以有两种方法：无套利方法和风险中性定价方法。1.无套利定价法由于期权和标的资产的风险源是相同的，在如图8.1的单步二叉树中，我们可以构造

4、一个证券组合，包括股资产多头和一个看涨期权空头。如果我们取适当的值，使则无论资产价格是上升还是下跌，这个组合的价值都是相等的。也就是说，当时，无论股票价格上升还是下跌，该组合的价值都相等。显然，该组合为无风险组合，因此我们可以用无风险利率对贴现来求该组合的现值。在无套利机会的假设下，该组合的收益现值应等于构造该组合的成本，即将代入上式就可得到：其中2.风险中性定价法在第六章中我们已经探讨过，期权定价可以在风险中性世界中进行，同样，我们也可以在二叉树模型中应用风险中性定价原理，确定参数、和，从而为期权定价。这是二叉树定价的一般方法。在风险中性世界里：（1） 所有可交易证券的期望收益都是无风险利率

5、；（2） 未来现金流可以用其期望值按无风险利率贴现。在风险中性的条件下，标的证券的预期收益率应等于无风险利率，因此若期初的证券价格为，则在很短的时间间隔末的证券价格期望值应为。因此，参数、和的值必须满足这个要求，即： (8.1)二叉树模型也假设证券价格遵循几何布朗运动，根据第六章的讨论，在一个小时间段内证券价格变化的方差是。根据方差的定义，变量的方差等于，因此： (8.2)式（8.1）和(8.2)给出了计算、和的两个条件。第三个条件的设定则可以有所不同， Cox、Ross和Rubinstein所用的条件 这是二叉树模型中最常用的第三个条件，后文我们将会谈到对第三个条件的其他设定方法。是： （8

6、.3）从以上三个条件求得，当很小时： （8.4） （8.5） （8.6）从而比较以上两种方法，我们可以看到，无套利定价法和风险中性定价法具有内在一致性。在风险中性定价过程中，我们无需考虑资产价格上升和下降的概率，也就是说资产预期收益具有无关性，这正好符合风险中性的概念。但是在最后的期权公式中，两种方法都包含了概率，这里的概率是风险中性世界中的概率，参数、和实际上都隐含在给定条件中。一般来说，在运用二叉树方法时，风险中性定价是常用的方法，而无套利定价法则主要是提供了一种定价思想。（二）证券价格的树型结构应用多步二叉树模型来表示证券价格变化的完整树型结构如图8.2所示。图8.2 资产价格的树型结构

7、当时间为0时，证券价格为。时间为时，证券价格要么上涨到，要么下降到；时间为2时，证券价格就有三种可能：、（等于）和，以此类推。一般而言，在时刻，证券价格有种可能，它们可用符号表示为： 其中注意：由于，使得许多结点是重合的，从而大大简化了树图。（三）倒推定价法得到每个结点的资产价格之后，就可以在二叉树模型中采用倒推定价法，从树型结构图的末端T时刻开始往回倒推，为期权定价。由于在到期时刻的预期期权价值是已知的，例如看涨期权价值为，看跌期权价值为，因此在风险中性条件下在求解时刻的每一结点上的期权价值时，都可通过将时刻的期权价值的预期值在时间长度内以无风险利率贴现求出。同理，要求解时的每一结点的期权价

8、值时，也可以将时的期权价值预期值在时间内以无风险利率r贴现求出。依此类推。采用这种倒推法，最终可以求出零时刻（当前时刻）的期权价值。以上是欧式期权的情况，如果是美式期权，就要在树型结构的每一个结点上，比较在本时刻提前执行期权和继续再持有时间，到下一个时刻再执行期权，选择其中较大者作为本结点的期权价值。例8.1假设标的资产为不付红利股票,其当前市场价为50元,波动率为每年40%，无风险连续复利年利率为10%，该股票5个月期的美式看跌期权协议价格为50元，求该期权的价值。为了构造二叉树，我们把期权有效期分为五段，每段一个月（等于0.0833年）。根据式（8.4）到（8.6），可以算出：据此我们可以

9、画出该股票在期权有效期内的树型图，如图8.3所示。在每个结点处有两个值，上面一个表示股票价格，下面一个表示期权价值。股价上涨概率总是等于0.5076，下降概率总是等于0.4924。在时刻，股票在第个结点（）的价格等于。例如，F结点（）的股价等于。在最后那些结点处，期权价值等于。例如，G结点（）的期权价格等于5035.36=14.64。图8.3 不付红利股票美式看跌期权二叉树从最后一列结点处的期权价值可以计算出倒数第二列结点的期权价值。首先，我们假定在这些结点处期权没被提前执行。这意味着所计算的期权价值是时间内期权价值期望值的现值。例如，E结点（）处的期权价值等于：而F结点处的期权价值等于：然后

10、，我们要检查提前执行期权是否较有利。在E结点，提前执行将使期权价值为0，因为股票市价和协议价格都等于50，显然不应提前执行。因此E结点的期权价值应为2.66元。而在F结点，如果提前执行，期权价值等于50.0039.69元，等于10.31元，大于上述的9.90元。因此，若股价到达F结点，就应提前执行期权，从而F结点上的期权价值应为10.31元，而不是9.90元。用相同的方法我们可以算出各结点处的期权价值，并最终倒推算出初始结点处的期权价值为4.48元。如果我们把期权有效期分成更多小时段，结点数会更多，计算会更复杂，但得出的期权价值会更精确。当非常小时，期权价值将等于4.29元。（四）二叉树方法的

11、一般定价过程下面我们给出用数学符号表示的二叉树期权定价方法，仍然举无收益证券的美式看跌期权为例。假设把该期权有效期划分成N个长度为的小区间，令表示在时间时第j个结点处的美式看跌期权的价值，我们将称为结点的期权价值。同时用表示结点处的证券价格。由于美式看跌期权在到期时的价值是，所以有：，其中当时间从变为时，从结点移动到结点的概率为，移动到的概率为。假定期权不被提前执行，则在风险中性条件下：其中。如果考虑提前执行的可能性的话，式中的必须与期权的内在价值比较，由此可得：按这种倒推法计算，当时间区间的划分趋于无穷大，或者说当每一区间趋于0时，就可以求出美式看跌期权的准确价值。根据实践经验，一般将时间区

12、间分成30步就可得到较为理想的结果。二、基本二叉树方法的扩展（一）有红利资产期权的定价1.支付连续红利率资产的期权定价当标的资产支付连续收益率为的红利时，在风险中性条件下，证券价格的增长率应该为，因此式（8.1）就变为：同时，式（8.4）变为： （8.7）式（8.5）和（8.6）仍然适用。对于股价指数期权来说，为股票组合的红利收益率；对于外汇期来说，为国外无风险利率，因此式（8.5）至（8.7）可用于股价指数和外汇的美式期权定价。对于期货期权，布莱克曾证明，在对期货期权定价时，期货可以看作支付连续红利率的证券参见F.Black, “The Pricing of Commodity Contra

13、cts,” Journal of Financial Economics，3（March 1976），16779）。，即式（8.7）可以进一步修正为： （8.8）这样式（8.5）、（8.6）和（8.8）就可用于美式期货期权的定价。2.支付已知红利率资产的期权定价若标的资产在未来某一确定时间将支付已知红利率（红利与资产价格之比），我们只要调整在各个结点上的证券价格，就可算出期权价格。调整方法如下：如果时刻在除权日之前，则结点处证券价格仍为：如果时刻在除权日之后，则结点处证券价格相应调整为： 对在期权有效期内有多个已知红利率的情况，也可进行同样处理。若为0时刻到时刻之间所有除权日的总红利支付率，则

14、时刻结点的相应的证券价格为：3. 已知红利额若标的资产在未来某一确定日期将支付一个确定数额的红利而不是一个确定的比率，则除权后二叉树的分支将不再重合，这意味着所要估算的结点的数量可能变得很大，特别是如果支付多次已知数额红利的情况将更为复杂（见图8.4）。图8.4 假设红利数额已知且波动率为常数时的二叉树图为了简化这个问题，我们可以把证券价格分为两个部分：一部分是不确定的，而另一部分是期权有效期内所有未来红利的现值。假设在期权有效期内只有一次红利，除息日在到之间，则在时刻不确定部分的价值为： 当时 当时 （8.9）其中表示红利。设为的标准差，假设是常数，用代替式（8.4）到（8.6）中的就可计算

15、出参数、和，这样就可用通常的方法构造出的二叉树了。通过应用式（8.9），把未来收益现值加在每个结点的证券价格上，就会使的二叉树图转化为的二叉树。假设零时刻的值为，则在时刻：当时，这个树上每个结点对应的证券价格为： 当时，这个树上每个结点对应的证券价格为： 这种方法和我们曾经分析过的在已知红利数额的情况下应用布莱克舒尔斯公式中所用的方法一致，通过这种分离，我们可以重新得到重合的分支，减少结点数量，简化了定价过程。同时，这种方法还可以直接推广到处理多个红利的情况。（二）利率是时间依赖的情形在二叉树的风险中性定价过程中，无风险利率是非常重要的参数，通常我们假设其为常数。但是当利率的期限结构呈陡峭的上

16、升或下降趋势的时候，这个假设显然是不合理的。更合理的假设是，即在时刻的结点上，其应用的利率等于到时间内的远期利率（该远期利率是站在零时刻得出的）。由于和并不随着的改变而改变，所以这一假设并不会改变二叉树图的几何形状，改变的是上升和下降的概率，我们仍然可以象以前一样构造出二叉树图：在定价过程中，贴现率也要相应地改为远期利率，其他过程和普通的二叉树方法相同。类似地，在为指数期权、外汇期权和期货期权定价时，可以对红利率和外汇无风险收益率做类似的修正，使其成为时间的函数，进而为这些期权定价。三、构造树图的其他方法和思路（一）的二叉树图在式（8.1）到（8.3）中，前两个式子是确定参数、和的固定条件，而

17、第三个条件是人为给定的，也是最常用的条件，但它并不是唯一的。我们也可以放弃这个假设，转而令，当的高阶小量可以忽略时，我们得到：这种方法的优点在于无论和如何变化，概率总是不变的，缺点在于二叉树图中的中心线上的标的资产价格不会再和初始中心值相等。（二）三叉树图另一种替代二叉树图的方法是三叉树图法，该树图的形状如图8.5所示。在每一个时间间隔内证券价格有三种运动的可能：从开始的上升到原先的倍，即到达；保持不变，仍为；下降到原先的倍，即。、分别为每个结点价格上升、持平和下降的概率。当的高阶小量可以忽略时，满足资产价格变化均值和方差的参数分别为：三叉树图的计算过程与二叉树图的计算过程相似。可以证明：三叉

18、树图的方法与我们第三节将要介绍的显性有限差分方法是一致的。图8.5 资产价格的三叉树图（三）控制方差技术控制方差技术是数值方法的一个辅助技术，可以应用在二叉树模型、蒙特卡罗模拟和有限差分方法上。其基本原理为：期权A和期权B的性质相似（比如其他条件都相同的欧式期权和美式期权，以及在第九章中将谈到的几何平均亚式期权和算术平均亚式期权），我们可以得到期权B的解析定价公式，而只能得到期权A的数值方法解。用代表期权B的真实价值（解析解），表示关于期权A的较优估计值，和表示用同一个二叉树、相同的蒙特卡罗模拟或是同样的有限差分过程得到的估计值。这时，我们假设用数值方法计算出的期权B的误差应等于用数值方法计算

19、出的期权A的误差： 进而得到期权A 的更优估计值为：可以证明，当和之间的协方差较大时，也就是说这个方法减少了对期权A的价值估计的方差，我们利用和的信息改进了对期权A的价值的估计。可以看出，控制方差技术实际上是利用数值方法计算两个类似期权之间的价格差异而不是计算期权价格本身。虽然从计算工作量来看，我们需要计算两个估计值和，但是由于两个期权的性质相似或路径相同，实际增加的工作量并不大。（四）适应性网状模型Figlewski和Gao S.Figlewski and B. Gao, “The Adaptive Mesh Model: A New Approach to Efficient Option

20、 Pricing,” Journal of Financial Economics, Forthcoming.提出了一种适应性网状模型（The Adaptive Mesh Model）来改进数值估计方法的效率。他们的方法是：在使用三叉树图为美式期权定价时，当资产价格接近执行价格时和接近到期时，用高密度的树图来取代原先低密度的树图。即在树图中那些提前执行可能性较大的部分，将一个时间步长进一步细分，如分为，每个小步长仍然采用相同的三叉树定价过程，这样使得树图更好地反映了实际情形，从而大大提高了定价的效率和精确程度。（五）隐含树图在第七章中，我们讨论了隐含波动率问题。隐含波动率是在假设期权市场有效的

21、条件下，从期权的市场价格倒推出其中隐含的波动率信息。类似地，隐含树图方法 参见E. Derman and I. Kani, “The Volatility Smile and Its Implied Tree,” Quantitative Strategies Publications, Goldman Sachs, (January 1994); E. Derman and I. Kani, “Riding on a Smile,” RISK, (February 1994), 32-39; M. Rubinstein, “Implied Binomial Trees,” Journal o

22、f Finance, 49, 3 (July 1994), 771-818.（Implied Trees）通过构建一个与目前市场上的期权价格信息相一致的资产价格树图，从而得到市场对标的资产价格未来概率分布的看法。其具体方法是在二叉树图中，通过前一时刻每个结点的期权价格向前推出（注意不是倒推）下一时刻每个结点的资产价格和相应概率：假设我们已知时刻各结点的值。在时刻应有个结点。这时，我们要求个变量：1. 时刻个结点上的资产价格2. 在和之间个资产价格向上运动的概率（向下运动的概率可以用1减去向上运动的概率得到）。为了求解这个变量，树图的构造应满足：1. 在和之间使用远期利率而非一个常数利率，在时刻

23、各结点资产价格的预期收益率应等于以上这个远期利率，这可以写出个方程；2. 树图的构造应使得市场中在时刻到期的个欧式期权准确定价，这些期权的执行价格等于时刻各结点上的资产价格，这可以再构造出个方程；3. 应保证树图的中心结点等于零时刻的初始资产价格。这样，我们得到个方程，可以解出以上的个变量，从而把树图的构造向前推进了一个时间步长。以此类推，最后我们可以得到整个树图，了解市场对未来资产价格分布和波动率的看法。值得注意的是：由于这是从市场期权价值信息倒推出来的隐含树图，因此不存在统一的上升下降概率和上升下降幅度，即、和是各不相同的。这也导致了有时会出现负的概率，需要对那个导致负概率的期权价格进行替

24、换。隐含树图的主要作用在于从交易活跃的常规期权中得到的关于波动率微笑和期限结构的信息，来为奇异期权定价。该树图可以和当天市场上的波动率矩阵相吻合，但是同时也隐含了未来的波动率微笑和期限结构，它们可能与当天市场上观测到的这些值很不相同。当一个期权的价值需要依赖于未来的某个波动率时，应用隐含树图的时候需要非常谨慎。四、二叉树定价模型的深入理解由上可见，二叉树图模型的基本出发点在于：假设资产价格的运动是由大量的小幅度二值运动构成，用离散的随机游走模型模拟资产价格的连续运动可能遵循的路径。同时二叉树模型与风险中性定价原理相一致，即模型中的收益率和贴现率均为无风险收益率，资产价格向上运动和向下运动的实际

25、概率并没有进入二叉树模型，模型中隐含导出的概率是风险中性世界中的概率，从而为期权定价。实际上，当二叉树模型相继两步之间的时间长度趋于零的时候，该模型将会收敛到连续的对数正态分布模型，即布莱克舒尔斯偏微分方程。取当前时刻为（这是为了后面计算的方便，并不影响结论），在给定参数、和的条件下（注意这里并未限定求、和的第三个条件，而是一般适用的），当时，二叉树公式：可以在进行泰勒展开，最终可以化简为：的高阶小量可以忽略，从而说明离散二叉树模型和连续布莱克舒尔斯模型是十分相似的，在时，二叉树模型收敛于布莱克舒尔斯偏微分方程。最后，二叉树模型和布莱克舒尔斯模型的另一个相似点在于：它们都可以通过选取适当的值，

26、构造一个由份的标的资产多头和一份期权空头组成的无套利组合。二叉树模型中的值满足；布莱克舒尔斯模型中的则满足，之后两者都可以利用这个无套利组合为期权定价。这里我们可以看到的极限就是，又一次验证了二叉树模型和布莱克舒尔斯模型的一致性。但是，三叉树模型则无法实现这样一个无套利组合，需要运用别的方法来构造。第二节蒙特卡罗模拟蒙特卡罗模拟是一种通过模拟标的资产价格的随机运动路径得到期权价值期望值的数值方法，也是一种应用十分广泛的期权定价方法。一、蒙特卡罗模拟的基本过程蒙特卡罗模拟要用到风险中性定价原理，其基本思路是：由于大部分期权价值实际上都可以归结为期权到期回报(payoff)的期望值的贴现，因此，尽

27、可能地模拟风险中性世界中标的资产价格的多种运动路径，计算每种路径结果下的期权回报均值，之后贴现就可以得到期权价值。以一个简单的欧式期权（即只有两个状态变量资产价格和时间，且利率为常数）为例，我们可以说明蒙特卡罗模拟的基本方法： 1. 从初始时刻的标的资产价格开始，直到到期为止，为取在风险中性世界中跨越整个有效期的一条随机路径。这就给出了标的资产价格路径的一个实现。2. 计算出这条路径下期权的回报。3. 重复第一步和第二步，得到许多样本结果，即风险中性世界中的期权回报的值。4. 计算这些样本回报的均值，得到风险中性世界中预期的期权回报值。5. 用无风险利率贴现，得到这个期权的估计价值。二、蒙特卡

28、罗模拟的技术实现（一）随机路径我们已经知道，风险中性世界中，标的资产价格变量所遵循的过程可以写作 注意，在风险中性世界中和在现实世界中，资产的波动率是一样的，改变的只是资产的预期收益率。 （8.10）也可以写作 （8.11）为了模拟的路径，我们把期权的有效期分为个长度为的时间段，则式（8.10）和（8.11）的近似方程分别为(8.12)（8.13）式（8.13）还可以写为 （8.14）其中代表时刻的价值，是从标准正态分布中抽取的一个随机样本。因此，蒙特卡罗模拟就是离散地模拟资产价格的时间序列，这样，只要得知初始时刻的值，随机抽取一个，就能算出时刻的值；接着时的值又能从时的值计算得到。因此，通过

29、个正态分布的随机取样就可以组建一个资产价格路径的蒙特卡罗模拟样本，并得到相应的回报值。重复以上的模拟至足够大的次数，计算回报值的平均值，折现后就得到了期权的期望值，另外也可以顺便得到所估计的期权价值的标准差。在以上两中模拟路径中，用比用本身更准确。使用式（8.12）模拟会存在的高阶小项的误差，仅仅在时是完全正确的。但是式子（8.13）或（8.14）却是精确的，因而对于所有的都是正确的。例8.1 假设无红利的股票价格运动服从式（8.12），年预期收益率为14，收益波动率为每年20，时间步长为0.01年，则根据式（8.12）有通过不断从标准正态分布样本中抽取的值，代入上式，我们可以得到股票价格运动

30、的一条路径。表8.1显示了模拟的一条特殊路径。假设股票价格的初始值为20，的第一个样本值为0.52，则第一个时间步长结束后，因此，第二步开始时的股票价格上升为20.236，这次抽到的为1.44，因此因此第三步开始时的股票价格变为20.847，以此类推，最终可以得到股票价格的一条模拟路径，其最后的价格21.124可以看成是10个时间步长或是年末股票价格分布中的一个随机抽样值。应当注意的是：表8.1仅仅表示了股票价格运动的一种可能方式，不同的随机取样将会导致不同的结果。进行这样的模拟达到足够多次，就可以得到年末股票价格的一个完整的概率分布。当然，由于我们采用的是式（8.12），因此只有在时，这种模

31、拟过程才是完全正确的。（该表格中假设股票价格精确到0.001）表8.1 当，时的股票价格模拟每步开始时的股票价格的随机抽样值该时间步长中的股票价值变化20.00020.23620.84720.51821.14620.88320.60320.71920.29220.61721.1240.521.440.861.460.690.740.211.100.731.162.560.2360.6110.3290.6280.2620.2800.1150.4270.3250.5071.111（二）单个变量和多个变量的蒙特卡罗模拟蒙特卡罗模拟的优点之一在于无论回报结果依赖于标的变量S所遵循的路径还是仅仅取决于的

32、最终价值，都可以使用这一方法。同时，这个过程也可以扩展到那些回报取决于多个标的市场变量的情况。1.当回报仅仅取决于到期时的最终价值时我们前面给出的蒙特卡罗方法可以模拟资产价格运动的整个路径，从而为那些价值可能依赖于路径发展过程的期权定价提供了便利，如我们将在后面看到的亚式期权（其期权价值依赖于有效期内的资产价格平均值）。如果期权价值仅仅取决于到期时的最终价值，而与其中间发展路径无关，比如最简单的欧式期权，蒙特卡罗模拟则可以非常方便迅速地进行处理。由于式（8.13）或（8.14）对所有的都是精确的，我们就可以直接用一个大步（）（假设初始时刻为零时刻）来多次模拟最终的资产价格，得到期权价值：（8.

33、15）2.当回报依赖于多个市场变量时当存在多个标的变量时，每次模拟运算中对每个变量的路径都必须进行抽样，从样本路径进行的每次模拟运算可以得出期权的终值。假设期权依赖于个变量，。定义为的波动率，为在风险中性世界中的期望增长率，为和之间的瞬间相关系数 变量、不一定是常数，可能取决于。与在单变量情况下一样，期权的有效期要分成个长度为的时间段，如果采用(9.12)的形式，的离散过程可以写为：（8.16）其中是从多元标准正态分布中抽取的一组随机样本，和之间的相关系数是，一次模拟运算包括从多元标准正态分布中获得个的样本。把这些样本值代入式（8.16）可以产生每个的模拟路径，并由此计算出期权的一个样本值。显

34、然，如果各个变量在期权有效期内的不同时刻到期，蒙特卡罗模拟也可以加以反映和处理，而且可以应用于任何一种随机过程（但是所有的随机过程都必须是各变量在风险中性世界中所遵循的过程），这正是蒙特卡罗模拟强大和适应性广泛的体现。（三）常数利率和随机利率的蒙特卡罗模拟对于一般的股票期权、股票指数期权、外汇期权和商品期权等，我们大多假设利率是常数。在风险中性世界中，标的资产的预期收益率是，之后折现时使用的也是无风险利率。假设到期时期权的回报为，则期权价值可以用公式表示为（初始时刻设为0）：. 其中，表示风险中性世界中的期望。以上是利率为常数的情况，如果期权模型中的变量之一本身就是短期无风险利率或是其他与有关

35、的变量，例如利率衍生产品，则蒙特卡罗模拟方法与前类似，只是要模拟风险中性世界中的路径，每次模拟时既要计算到期时终值相应带来的期权回报，又要计算期权有效期内的平均值。最后折现的时候使用的贴现率是这个平均值，用数学符号表示为：为有效期内瞬间无风险利率的平均值。（四）随机样本的产生和模拟运算次数的确定在介绍了蒙特卡罗模拟的主要过程和方法之后，下面我们来讨论两个细节问题：1. 的产生如前所述，是服从标准正态分布的一个随机数。大多数程序语言都为抽取0到1之间的随机数编制了程序。如果只有一个单变量，则可以通过下式获得：其中是0到1的相互独立的随机数。如果需要从二元标准正态分布中抽取样本，则可以用如下的方法

36、：和是用上述方法从单变量标准正态分布中抽取的独立样本，则其中是相关系数。如果从元标准正态分布中取样，同样先从单变量标准正态分布中抽取个独立变量，则为了使有正确的方差，并使与之间有正确的相关系数，必须满足且。令第一个样本等于，就可以解出这些的方程，通过和计算出，再通过、和计算出，依次类推。2. 模拟运算次数的确定蒙特卡罗过程是用随机数序列实现有限次数的模拟，进行模拟运算的次数取决于所要求的精度。在模拟运算得到期权价值均值的时候，我们也可以得到其标准差。设是如上所述进行运算的个数，为均值，是标准差，则期权估计值的标准误差为。如果对估计值要求95的置信度，则期权价值应满足可见，我们对期权价值的不确定

37、性与模拟运算次数的平方根成反比，如果要将精确度提高为原来的10倍，则模拟运算次数应为原来的100倍。实际中常用的一个值是10万次。三、减少方差的技巧如果蒙特卡罗模拟过程是按照我们上述的方法进行，往往需要极大，才能使期权价值的估计值较为精确合理，从而使得计算效率很低，因此人们运用了多种方法来降低估计的方差，以大大减少实验次数，加快收敛。其中最常用的是对偶变量技术和控制方差技术。（一）对偶变量技术对偶变量技术（Antithetic Variable Technique）是指在一次模拟运算中，计算两个期权价值：用通常方法计算得到，改变计算中所有的的符号得到，这条模拟路径得到的期权价值是和的平均值，期

38、权的最终估计值则是有限个的均值。通过这种方法，当一次模拟中的一个值高于真实值时，则另一个值必然偏低，反之亦然。人们发现，这种方法相当有效。如果是的标准差，是模拟运算的次数（每一次中包含计算两个值），估计值的标准差就是，这通常大大小于运算次得到的标准差。（二）控制方差技术蒙特卡罗模拟和二叉树模型中的控制方差技术是相同的，只是在蒙特卡罗中，对于两个类似的期权，必须使用相同的随机数流和相同的平行地进行两次模拟。（三）重点抽样法重点抽样法（Importance Sampling）适合于那些大部分路径对定价意义不大的情形，比如一个深度虚值的看涨期权，大部分路径上的终值为零。这时我们只选取那些标的资产的到

39、期日价值大于其执行价格的路径，即重要路径为期权定价。这样等于缩小了样本空间，从而加速了收敛。运用重点抽样法进行模拟时要注意，最后从重要路径中获得的贴现平均值还要再乘上（重要路径出现的概率），才能获得期权价值的最终估计值。（四）间隔抽样法间隔抽样法（Stratified Sampling）将市场变量在未来时刻的基本概率分布分为多个区间，并根据它的概率从每个间隔中抽样。假设存在10个可能性基本相等的区间，那么我们的抽样方案就可以设计为每个间隔中抽取的样本各占10。如果样本间隔的数量很多，我们就可以取每个区间内的均值或是中位数作为该区间的代表样本值。这样也可以提高模拟运算的效率。Curran和Mor

40、o的文章中介绍了间隔抽样的两种具体方法 参见M. Curran, “Strata Gems,” RISK, (March 1994), 70-71; B. Moro, “The Full Monte,” RISK, (February 1985), 57-58.。（五）样本矩匹配法样本矩匹配法（Moment Matching）是指对从标准正态分布中抽取的样本进行调整，使其一阶矩、二阶矩甚至高阶矩都匹配。假设我们为了计算特定时间内某个特定变量的变化值而使用的正态分布样本为，为了使前两阶矩相匹配，我们分别计算样本的均值和标准差，之后定义调整后的样本为经调整后的样本均值为零，标准差为1.0，其后所有

41、的运算中都使用这个经调整的样本。样本矩匹配法可以节省计算时间，但是由于所有的抽样值都要保存到模拟运算结束，存储负担很重。样本矩匹配法有时也被叫做二次抽样法（Quadratic Resampling），常常与对偶变量技术结合使用，因为后者可以自动地匹配奇数阶矩，这样样本矩匹配法就只需要匹配二阶矩，或者再加上四阶矩。（六）准随机序列抽样法准随机序列（Quasi-Random Sequences）是指从一个概率分布中抽取的代表样本组成的序列（它们在本质上是确定性的），准随机序列抽样法类似于间隔抽样，目的都是为了得到标的变量的代表值。只是在间隔抽样中我们假设已知需要抽取的样本数，而在准随机序列抽样法中

42、，每次抽样都试图填补之前已存在的样本之间的空缺，使得抽取的样本值总是能大致均匀地分布在整个概率空间中，这使得模拟的收敛速度得到了改进。由于数论中有一个差异（Discrepancy）的概念衡量的是点集分布的均匀性，因此准随机序列法又被称为低差异序列法。通过这种方法，标准差能从降低为。（七）树图取样法我们还可以通过一个树图来为资产价格路径取样，进行蒙特卡罗模拟。这样，在间隔抽样中，我们就可以应用树图而不是随机地抽取代表路径。如果一条路径经过任何给定结点的比例等于（或接近于）该结点被到达的概率，就被称为代表路径。Mintz(1997)最早介绍了这种方法。 具体方法参见D. Mintz, “Less

43、is more,” RISK, 10, 7 (July 1997), 42-45四、蒙特卡罗模拟的理解和应用蒙特卡罗方法的实质是模拟标的资产价格的随机运动，预测期权的平均回报，并由此得到期权价格的一个概率解。蒙特卡罗模拟的主要优点包括：1. 在大多数情况下，人们可以很直接地应用蒙特卡罗模拟方法，而无需对期权定价模型有深刻的理解，所用的数学知识也很基本；为了获得更精确的答案，只需要进行更多的模拟；无需太多工作就可以转换模型。以上这些优点使得蒙特卡罗方法成为一个相当广泛和强大的期权定价技术。2. 蒙特卡罗模拟的适用情形相当广泛，其中包括：（1）期权的回报仅仅取决于标的变量的最终价值的情况；（2）期

44、权的回报依赖于标的变量所遵循的路径，即路径依赖的情形；（3）期权的回报取决于多个标的变量的情况，尤其当随机变量的数量增加时，蒙特卡罗模拟的运算时间近似为线性增长而不象其他方法那样以指数增长，因此该方法对依赖三种以上风险资产的多变量期权模型很有竞争力。因此，蒙特卡罗模拟可以适用于复杂随机过程和复杂终值的计算，比如将在第九章介绍的奇异期权和路径依赖期权，同时，在运算过程中蒙特卡罗模拟还能给出估计值的标准误差，这也是该方法的优点之一。另一方面，蒙特卡罗模拟的缺点主要是：1. 只能为欧式期权定价，难以处理提前执行的情形。尝试使用蒙特卡罗模拟技巧来为美式期权定价，成为近年来这个领域的发展方向之一。2.

45、为了达到一定的精确度，一般需要大量的模拟运算。尤其在处理三个以下的变量时，蒙特卡罗模拟相对于其他方法来说偏慢，例如在第九章中处理一些路径依赖期权时，人们常常用二叉树模型等来取代蒙特卡罗模拟，就是因为其耗费的计算时间太多。而本节中那些减少方差的技巧，都是为了改进计算的效率而提出的。第三节有限差分方法在应用微分方程建模的学科和工程当中，有限差分方法是最常用的数值求解方法。在金融界，这个方法也日益受到人们的重视，越来越多地用在期权定价当中。其主要思想是：应用有限差分方法将衍生证券所满足的偏微分方程转化为一系列近似的差分方程，即用离散算子逼近、和各项，之后用迭代法求解，得到期权价值。在坐标图上，有限差

46、分方法则体现为格点（Grids），如图8.6所示。图8.6 有限差分方法的格点图具体地说，有限差分方法就是用有限的离散区域来替代连续的时间和资产价格：首先，把从零时刻（初始时刻设为零时刻）到到期日时刻之间的时间分为有限个等间隔的小时间段，设，就有共个时间段；其次，把资产价格的变化也分成个等间隔的小价格段，定义，就得到个资产价格 从数学意义上说，的范围是无限的，但是在经济意义上，资产价格一般会有一个特定的合理变化边界，超出这一边界的价格是没有意义的，因此只需考虑从0到价格上限之间的变化即可。这个上限也不需要太大，一般在执行价格的三到四倍就足够了。如果划分合理，初始的资产价格会落在零时刻的一个格点

47、上。其中，和是相互独立的。这样，我们就构造了一个共有个格点的图，时间、资产价格和期权价值都仅仅在相应的格点处离散计算。点对应时刻和资产价格，则表示处的期权价值。应用这些格点之间的关系和已知的边界条件，我们可以把连续偏微分方程转化为一系列的差分方程，逐次求解，就可以得到零时刻初始资产价格所对应格点的期权价值。下面我们将具体介绍如何使用这些格点逼近微分，求出期权价值。这可以用多种方法实现，其中包括隐性有限差分法、显性有限差分法和其他的一些方法。为了方便说明，我们使用一个无红利股票的美式看跌期权作为例子。实际上有限差分方法不仅适用于布莱克舒尔斯公式，也适用于其他类似的问题，如随机利率的情形等。一、隐

48、性有限差分法为了将偏微分方程化为差分方程，需要应用离散算子分别逼近、和。其中的一种近似方法是隐性有限差分法（Implicit Finite Difference Method）。（一）、和的差分近似1. 的近似对于坐标方格内部的点，期权价值对资产价格的一阶导数可以用三种差分来表示：、和这三种逼近方法分别称为前向差分近似（Forward Difference Approximation）、后向差分近似（Backward Difference Approximation）和中心差分近似（Central Difference Approximation）。可以看到，这三种方法是针对时刻进行的差分近似

49、，只是的取值方向不同，中心差分实际上是前两者的平均值。应用泰勒展开式考察这三种近似方法的精确度，可以发现前向和后向差分近似的误差均为的高阶小项，而中心差分的误差则为，精确度更高。这是因为这一近似方法关于对称，使得一些误差项可以相互抵消。因此，大多数时候人们采用中心差分法来逼近，但有时也根据需要使用单向差分方法。2. 的近似对于点处的，我们则采取前向差分近似以使时刻的值和时刻的值相关联：这一近似的误差是，可以进一步改进。但是在这里，这样的精确度已经足够了。3. 的近似点处的的后向差分近似为，因此点处期权价值对标的资产价格的二阶差分为这个二阶差分也是中心差分，其误差为。从以上三个近似我们可以发现，

50、除了对时间的差分涉及时刻和时刻的期权值，对资产价格S的一阶和二阶差分都只使用了时刻的不同期权值。（二）差分方程把以上三个近似代入布莱克舒尔斯偏微分方程，整理得到： （8.17）其中， , 由于和使用中心差分，整个方程的误差为。这个差分方程的思想可以在图8.7中得到体现。也就是说，根据这个方程，可以从时刻的三个相邻格点的期权价值（其对应的资产价格分别为、和）求出时刻资产价格为时的期权价值。因此隐性有限差分法可以理解为从格点图内部向外推知外部格点的期权价值。图8.7隐性有限差分方法 图8.8显性有限差分方法（三）边界条件1.从到期时刻的期权回报值可以得到边界上所有格点的期权价值。时刻看跌期权的价值

51、为其中 2. 当股票价格为零时，看跌期权的价值为，由此可以得到下方边界上所有格点的期权价值： 3. 当股票价格趋于无穷时，看跌期权的价值趋于零。可以近似认为上方边界上 （四）求解期权价值已知差分方程（8.17）和三个边界条件之后，我们的目的是求出格点图左边界的价值，其中的一个格点就是我们要求的期权价值。利用方程（8.17）和边界条件，我们可以写出时刻的个联立方程： 和时，时，由此，可以解出每个的期权价值，然后再与每个格点的期权内在价值进行比较，判断是否要提前执行，从而得到时刻每个格点的期权价值。依此类推，最后可以计算出，当等于初始资产价格时，该格点对应的就是我们要求的期权价值。例8.2 用隐性

52、有限差分方法解例8.1中的期权价值。其中,。可以从表格8.2中看到，无红利股票美式看跌期权的价格为4.07元。股票价格到期时间（月）54.543.532.521.510.501000.000.000.000.000.000.000.000.000.000.000.00950.020.020.010.010.000.000.000.000.000.000.00900.050.040.030.020.010.010.000.000.000.000.00850.090.070.050.030.020.010.010.000.000.000.00800.160.120.090.070.040.030.

53、020.010.000.000.00750.270.220.170.130.090.060.030.020.010.000.00700.470.390.320.250.180.130.080.040.020.000.00650.820.710.600.490.380.280.190.110.050.020.00601.421.271.110.950.780.620.450.300.160.050.00552.432.242.051.831.611.361.090.810.510.220.00504.073.883.673.453.192.912.572.171.660.990.00456.586.446.296.135.965.775.575.365.175.025.004010.1510.1010.0510.0110.0010.0010.0010.0010.0010.0010.003515.0015.0015.0015.0015.0015.0015.0015.0015.0015.0015.003020.0020.0020.0020.0020.0020.002