二叉树和三叉树的期权定价方法

第七章期权定价的二叉树和三叉树方法在这一章中，我们利用二叉树和三叉树方法为期权定价。在第2.1节中我们已经介绍了利用基础途径的二叉树方法解决期权价格不确定性的模型。二叉树方法依赖于对相关随机过程的离散化并利用计算和内存的结合以满足易于管理的要求。我们也在，我们必须把原来的单步格方法扩展到多步格方法，但是我们必须校对格使它能够反映出相关模型，且这个模型是连续时间、连续状态的随机微分方程。然后我们就可以推广到多步的二叉树格和三叉树格。

在7.1节中，我们从如何利用在离散概率分布的时刻下随机价格波动校准简单的二叉树格。从这点来看，弄清楚网格技术和蒙特卡洛模拟之间的联系是非常重要的，而利用时刻匹配技术缩减方差可以看作一种快捷的抽样排序。然后我们讨论内存效率的实现是如何设计的，美式期权定价是7.2节的主题。同时，还是要注重它和其他技术方法的联系。现在我们要做的本质上是一个非常简单满足动态规划原则的程序，我们将在第10章程序中进一步拓展。在7.3节中，我们把上述方法推广到双标的资产的情形，虽然这是一个最简单的情形，但是我们可以从这个情形中看出内存控制是这一情形的基础。另一种一般化的代表是三叉树格方法，三叉树格方法可以作为一种更普遍的有限差分方法（具体将在，最后，我们在7.5节中具体讨论网格化方法的优势和劣势。

期权定价的二叉树和三叉树格方法

图7.1单时期二叉树格

7.1二叉树定价方法

在，我们已经考虑过单步二叉树方法在无套利情况下的期权定价，这里我们为了方便直接利用图7.1。其主要思想是复制两个资产，

一个是无风险资产，另一个是相关股票。利用这两项资产，我们可

以通过它们的组合塑造任何收益率的资产。如果我们令u和d为任意两个价格的角标，我们可以看到期权的价格应该为f

则，

f0=e-rδt[pf u+(1-p)f d]（7.1）

在公式7.1中f

u 和f

d

是标的资产在涨跌两种情况的期权价格，p是

风险中性前提下相关资产升值的概率。

为了寻找一个更好的不确定性模型，我们可以增加分类的情况，复制期权收益，甚至我们可以使用更多的资产，或允许中间日期交易。第二种可能性更为实际，并且也是必不可少的，例如，对于在期权的存续期内可以随时执行的美式期权来说。对其求极限，就会得到连续时间模型，并且其最后收敛于Black—sholes方程。当Black—sholes方程没有解析解的时候，我们必须采取一些离散化的途径，比如说可以通过蒙特卡洛模拟从而估计出风险中性条件下预期收益，或者建立一个自适应网格的有限差分方法去解决相应的PDE模型。就像我们在图7.2中展示的一样，多级二叉树格方法就是一种可以选择的离散化方法。我们也可以考虑利用树图，但是要注意使计算方法易于控制。

二叉树格定价

图7.2新生成的二叉树图

这里我们为了方便令u=1/d。虽然这个不是必须的，但是在后面我们可以看到，这个假设令模型简化了很多即每上一步紧接着下

δ t 一步都会得到相同的初始价格。

正如我们从图中看到的一样，我们仅用了有限个价格步。这个

有可能就是实施该方法的优势。但是，我们该怎么恰当的确定 u 和

d 的值呢，我们应该利用近似相关的连续过程去校对网格。

二叉树格方法应该是风险中性过程一个良好的相似。

因此，我们应以这样的方式参数设置晶格，即保持着连续时间模型

的一些基本属性，这一过程就叫做校准。从 S t 开始，经过一个小的

时间步 δt ，从 2.5 节我们可以看到新价格是一个随机变量 S t +δt ，且

利用对数正态对数分布的特性，我们得到

E [S t +δt / S t ] = e r δt

(7.2) 和

Var [S t +δt / S t ] = e 2r δt (e σ

2 - 1) （7.3）

一个合理的要求是这些离散的动态点必须和它们的时刻相匹配。要

注意的是，我们只有两个个等式，却有 3 个参数，p,u 和 d,所以三

个变量有一个为自由变量，我们令 u = 1/ d ，这样做是为了计算简便，

但不是必须的。在网格点上，我们有：

E [S t +δt ] = pu .S t + (1 - p )d .S t ，

和（7.2）联立得

注意， p 是风险中性条件下的概率，它不依赖于真实浮动，为了和

方差匹配，在晶格上我们看到

从（7.3）中我们也可以看到

把最后两个等式联立可得

最终得到

将 p 带入最后一个等式的右侧，化简得

最后我们得到这样的等式

其中，利用 u = 1/ d ，可以转化为二次方程：

方程的一个跟为

利用一阶条件拓展，只受 δt 的影响，我们可以简化表达式，对平方

根近似化简可得

因此

但是对于二阶条件，我们对 e σ δt 拓展，最终获得参数

u = e σ δt ， d = e -σ δt ， p = e r δt - d u - d ， （7.4）

这就是著名的 CRR 公示

这里强调一下：这个方法以及文献中所用的参数都不是唯一的，

例如我们可以取 p = 0.5 ，经计算可得：

p = 0.5 ， u = e (r - σ 2 2 )δt +σ δt ，

d =

e (r - σ 2 2 )δt -σ δt 这就是杰诺-拉德参数，此外，我们一直在努力结束涉及计算以

及线性方程组的计算，通过对数转换的方法，我们尽量的避免这些

困难。在以后，我们都将采用这个方法。

假设无风险利率和波动是时间常数，我们所得的结果适用于整

个晶格参数，为一个期权定价，我们需要对标的资产制定一个网格，