跳跃—扩散模型资产定价公式的数值计算方法

龙源期刊网 http://www.qikan.com.cn

跳跃—扩散模型资产定价公式的数值计算方法

作者:张鸿雁,李 强,张 志 来源:《经济数学》2010年第02期

摘 要 假定资产价格变化过程服从跳跃扩散过程,那么基于它的欧式期权就满足一个偏积分微分方程(PIDE),本文利用差分法来离散这个PIDE方程,用两种迭代方法得到方程的数值解:基于雅可比正则分裂法和预条件共轭梯度法.

关键词 跳跃扩散模型;差分法;FFT算法;欧式看涨期权 中图分类号 O241.82文献标识码:A

引 言

美国芝加哥大学教授Black和Scholes[1]在1973年发表了“The Pricing of Options and Corporate Liabilities”一文,提出了著名的

期权定价公式,在

公式中,假设股票

的价格过程是连续的几何布朗运动,但是,在现实市场中,一些突发情况会引起股票价格发生跳跃.基于上述考虑,Merton在1976年首先提出了跳跃扩散模型,在Merton模型中,资产价格在没有受到外界重大影响时服从布朗运动,当资产价格受到突发事件的影响而发生跳跃时,就用跳跃过程来描述.

本文首先介绍PIDE[2]的具体形式,在Merton模型下对其差分离散,得到一个Toeplitz矩阵方程,用两种方法解这个矩阵方程,一是基于雅可比正则分裂的迭代方法[3-4],二是预条件共轭梯度方法.考虑到Toeplitz[5]矩阵的特殊性,在迭代的过程中,将其植入到一个循环矩阵中,利用循环矩阵和向量的乘积来计算Toeplitz和向量的乘积,而循环矩阵向量乘积可以通过快速富里叶变换(FFT)快速计算,这样就加快了迭代速度.共轭梯度法是解决Toeplitz线性方程组的主要方法之一,在利用共轭梯度法的情况下,快速傅里叶变换的作用是加快共轭梯度法的迭代速度,但不改变其收敛速度,共轭梯度法的收敛速度取决于线性方程组系数矩阵的条件数,基于此考虑,本文采用预条件共轭梯度算法,选用R.Chan优化循环预条件器[6],预条件器的使用是为了改善系数矩阵的条件数,以便提高收敛速度. 2 跳跃扩散模型

假设市场是完备无套利的市场,在跳跃扩散模型下,资产的价格变化过程服从随机微分方程

-

龙源期刊网 http://www.qikan.com.cn

其中,υ(t)是漂移率,σ(t)是波动率,ω(t)标准布朗运动是1-独立的.

由文献[7]可知,在上述假设下,基于资产价格S与时间τ的未定权益V(S,τ)满足

--

的概率是

函数,是一个随机变量,用ζ表示平均跳跃幅度E(η-1),泊松过程

是泊松过程

的概率

是泊松到达强度,η-1是由S跳跃到Sη的跳跃幅度

与布朗运动ω(t)是相互

式中,t=T-τ是到期时间为T的时间,r是无风险利率,g(η)是跳跃幅度η的密度函数. 对式(2)的积分部分进行指数变量变换,令

再对其余部分进行变换,令f(y)是跳跃幅度

(t,x)∈[0,T)×R,边界条件 (τ,x)∈

3 Merton模型下的有限差分离散 在

中,由于称为截断点

通常限定x的取值范围是

---

-(r---

令u(τ,?)=υ(T-τ,?),则式(4)变为-

的密度函数,则式(3)变为---

--

函数

则式(2)变为

称为截断区域,式(5)中的积分部分可以分割为两部分

龙源期刊网 http://www.qikan.com.cn

\\

在计算欧式看涨期权的情况下,在R\\

上的积分u(τ,z)

可以进行近似代替

--当x→+

时.(6)

u(τ,x)→0, 当x→-

经 济 数 学第 27 卷

第2期张鸿雁等:跳跃扩散模型资产定价公式的数值计算方法 对于式(5)中的积分部分,进行变量变换z=x+y,则 -

定义函数

--- 在模型中

--

在的情况下,有表达式

--

其中是标准分布函数

-

- 考虑时间和空间节点,使--

-(---1),和

-1)k,m=1,…,q,

k=T/q.记-

在[-上用复合梯形原则,有积分近似

--

联系客服:779662525#qq.com(#替换为@) 苏ICP备20003344号-4