资源描述
股票价格运行的幂律特征及幂律跳跃扩散模型* 本文受到中国国家自然基金项目(项目编号:70801066;71071168)、教育部人文社科规划基金、中山大学高校科研基本业务费支持
作者简介:曹宏铎(1972—),男,河北秦皇岛人,博士,副教授,研究方向:金融工程.Email:caohan@
曹宏铎 何智 李旲
(中山大学管理学院,广州 510275)
摘要:在Merton提出的跳跃扩散模型的逻辑框架之下,进行了两个方面修正:将计数过程由过程修正为带有幂律性质的更新过程,同时对跳跃的幅度也进行修正,赋予股票价格运动过程发生跳跃的时间和幅度以幂律分布特征;通过实证研究,本文所做出的两个修正可以更加准确的描述股票价格的运动过程,可以同时得到具有尖峰胖尾的收益率分布和波动聚集性;以此为基础可以更加准确地为期权等金融衍生品进行定价,同时也对金融风险管理提供有效的工具。
关键词:风险;期权定价模型;跳跃扩散模型;更新过程;幂律分布
0 引言
早在20世纪初,法国数学家Bachelier开始利用运动来研究股票价格的运动[1].Fama提出了著名的有效市场假说[2],基于有效市场假说,Black和Scholes提出了期权定价模型,可以得到欧式看涨期权和欧式看跌期权的价格解析解[3],对金融理论和实务的发展产生了重大影响.
大量的实证研究对传统金融经济学的研究假设和框架提出了质疑与批评.Osborne发现股票价格收益率分布具有胖尾的特征[4],同时很多学者也通过实证研究发现股票收益率分布的尖峰胖尾特性.因此,假设股票收益率服从正态分布的模型也需要进行一定的修正.对模型的修正主要从两个方面出发:引入波动率的随机性和股票收益的非连续性.
Engle提出模型[5],Bollerslev将其进行推广得到模型[6].Cox和Ross提出价格依赖的波动率扩散模型[7], Hull和White提出随机波动率模型[8],Heston建立仿射随机波动率模型[9],而模型则可以被看做随机波动率模型的特殊形式.Peters提出分数布朗运动[10],但Rogers指出该模型存在套利机会[11].Merton提出要在股票价格的运动中加入属于非系统风险跳跃的过程,认为股票价格的运动过程不是连续的[12],Kou提出期权定价的双指数跳跃扩散模型[13].Bates认为将随机波动率模型和跳跃扩散模型相互结合,则可以更好的刻画出股票价格运动方式[14],Scott[15]、Duffie 等[16]、Bakshi[17]提出了各自的结合模型.Barndorf-Nielsen提出广义双曲分布[18],Eberlein和Keller将广义双曲分布的子族双曲分布用于金融模型分析[19],Blattberg和Gonedes[20]、Bibby和Sørensen[21]也进行了相关的研究.
叶中行等[22]、杨宏林和陈收[23]、钟春仿和佟孟华[24]通过实证研究发现,股票价格收益率的尾部分布呈现出幂律特征,而对尾部分布的研究则具有很强的现实意义,尤其体现在的应用领域[25].的使用需要对股票收益率损失部分的尾部分布做出规定,而目前普遍的应用方法是让损失的尾部分布服从正态分布或历史统计数据分布.幂律分布为股票收益率的尾部分布赋予了较高的概率,这样就可以使可以更好的应用于金融风险管理.
同时,Barabási[26]、Vázquez等[27]发现个体的行动也具有幂律特征,即事件的发生具有在长时间的“了无一事”中存在高频度阵发活动的特点。个体行动的幂律性体现在两次事件发生之间的等待时间,等待时间的长短则具有幂律分布的特点.股票价格的运动也符合这样的的结论,即允许股票价格可以在较长的时间里不发生跳跃,也可以在短时间内密集的发生跳跃.相比之下,广泛的运用在对事件发生的描述的过程则无法刻画出这样的规律.本文将人类行为动力学引入到金融学理论中,更确切的说是对股票价格的跳跃扩散模型进行修正.
1 模型介绍
实证研究表明股票价格的变化并不服从对数正态分布,而这正是模型的重要假设,这一假设隐含认为股票价格的变化路径是连续的.Merton将路径连续的扩散过程和随机发生的跳跃过程相互结合,提出跳跃扩散模型,如式(1)所示。
, (1)
其中, 和为上可积函数,分别表示风险资产收益率和风险资产收益率波动率;是标准运动;是一个计数过程,表示发生跳跃的次数;表示股票价格发生跳跃的幅度,是具有独立同分布的非负随机变量序列。
即是不确定性定义在完备概率空间上的市场中时刻的股票价格,该完备概率空间由计数过程、标准运动和独立同分布的随机变量列组成.简言之,Merton模型由线性漂移、布朗运动和复合泊松过程三部分构成。
Merton(1976)提出的模型中,是一个Poisson过程,服从对数正态分布。该模型形式简洁、逻辑清晰,因此,自发表后便成为分析标的资产价格不连续变化的标准模式。在此框架下的后续研究一般接受的假设,而针对的分布进行研讨。其中,Kou(1999,2002)提出的双指数跳跃扩散模型影响最大。
大量实证研究结果表明资产价格行为中幂律的存在,同时,在其他领域中的研究也表明幂律是人类行为规律的一种普适的现象。本文认为是一个更新过程,更新过程的事件间隔时间服从含有幂律性质的分布,其概率密度函数为
, (2)
即方程[28].
过程假设股票市场历史上的波动情况与现在的波动情况无关,关注于等待时间的长短,而不是起始时刻的初始状态.但是,市场中的交易主体会对市场信息反映过度或者反映不足,他们对股票市场走势的预期往往基于历史交易资料和当前市场的宏观环境以及所有能够掌握的信息,而任何信息的不完备或是具有偏差的解读都会造成事件发生的时间记忆特点.从这个角度来讲,本文所定义的更新过程比过程可以更加准确的模拟出股票价格的跳跃情况.
跳跃幅度也具有幂律性质。令,则是非对称的分布,其概率密度函数为
, (3)
其中,,,,,,,,和分别为股票价格上跳和下跳的概率,价格向上和向下发生跳跃幅度的概率密度函数分别为和,、为常数.
经过以上定义后的式(1)称为幂律跳跃扩散模型,其所描述的股票价格运动过程可以得到收益率的胖尾特征,也可以得到波动聚集性,从而该模型对金融资产定价理论和实践具有重要的意义.
式(1)的解为:
. (4)
利用式(4),可以得到收益率的表达式:
(5)
为了分析的简便,假设资产期望收益率和波动率为常数,将式(5)进行幂级数展开,利用,可以得到:
, (6)
其中,为服从标准正态分布的随机数,为服从伯努利分布的随机数,即有和,的概率密度函数由式(3)定义.式(6)中的概率密度函数可以进行如下表示:
(7)
可以看出,式(7)所展示的股票收益率概率密度函数实际上是对正态分布的一种修正,修正的幅度由式(8)所示.
(8)
将修正函数加入到正态分布中,即可得到股票收益率的概率密度函数,修正函数的图形通常呈现出形.图1所展示的股票日收益率的概率密度函数呈现出明显的胖尾特征.
2 模型参数估计及实证研究
股票价格跳跃扩散模型认为跳跃行为是带有幂律特征的更新过程,跳跃幅度的分布也具有幂律特征。本文以上证综合指数从2006年1月4日到2009年7月2日共849个交易日的收盘点数为样本,对模型进行实证检验并进行参数估计。
2.1 :具有幂律特征的更新过程
本文对的定义,蕴含股票价格收益率在经过极短或极长的等待时间后发生跳跃的概率较小,而每两次跳跃的间隔时间则更加倾向于某一段时间的假设.我们以上证综合指数为对象,将样本按照周分为173组。
1.根据式(2)构造待检验的模型:
. (9)
为了分析的简便,式(9)中的为每两次跳跃之间的等待天数,为等待天后发生跳跃这种情形的频率.
2.跳跃识别:
按照式(10)和(11)可以得到每周的指数波动率,其中为每周第天的收盘点数,为每周第天的收益率。
(10)
(11)
通过构造模型来得到预期的指数波动率.在式(12)中,为第周的收益率,和分别为第周和第周的指数波动率,为代表长期波动因素的常数项,为残差项.
(12)
利用第3周到第173周共171组数据对式(12)进行拟合,回归结果见表1.
表1 历史收益率和波动率对当期波动率的回归结果
系数估计
检验值
值
0.000276
5.365538
0.0000
0.062437
4.122887
0.0001
0.229602
2.981498
0.0033
0.222697
0.213443
24.06595
根据回归结果,第周的指数波动率可如下表示:
(13)
利用式(13)可以计算出对第周的预期指数波动率,然后将每周第天的与该周的相比较,如果,则说明第天的收盘指数发生了跳跃,而如果,则说明第天的收盘指数没有发生跳跃.我们根据这一标准对从2006年1月16日到2009年7月2日共841个交易日进行收盘点数是否发生跳跃的判断,并统计出两次跳跃之间的等待天数及其频数,统计结果见表3.
3. 等待时间分布的参数估计
由此,我们得到18组数据,根据式(9)进行回归分析,回归结果见表2.
表2 等待时间概率分布的回归分析
系数估计
0.839324
2.064301
1.983040
-1.944025
-0.176437
检验值
1.422957
9.730175
2.483533
-6.931492
-8.221085
值
0.1752
0.0000
0.0253
0.0000
0.0000
0.949222
0.808581
0.942451
0.796617
140.2004
67.58624
根据回归结果,我们可以得到等待时间分布的概率密度函数:
. (14)
该概率密度函数的图形如图1所示.
在此,我们注意到,为了保证定积分结果为1,满足概率密度函数的性质,我们将式(14)修正为式(15).
(15)
由计算在经过天的等待后指数发生跳跃的概率,并由计算累积概率,计算结果见表3.
a b
c d
e f
图1 等待时间分布的概率密度函数
注:图中虚线代表本文所定义的分布,实线代表指数分布.(a)虚线为分布的概率密度函数,由式(15)给出其具体表达式,实线为指数分布的概率密度函数,由式给出其具体表达式;(b)显示出分布和指数分布概率密度函数在尾部的图形,可以看到分布较指数分布在尾部赋予更大的权重;(c)给出在双对数坐标下,分布概率密度函数的图形,可以看出其为一条斜率为-2的直线,显示出幂律分布的特点;(d)在半对数坐标下,指数分布的概率密度函数为一条斜率为负的直线;(e)将回归分析所使用的18个散点加入到(a)中;(f)将(e)放置在双对数坐标系下所显示出的图形.
表3 等待天数的频数及概率
天数
频数
频率
分布概率 累积概率
指数分布概率 累积概率
1
65
0.296804
0.1523
0.1523
0.1617
0.1617
2
44
0.200913
0.2465
0.3988
0.1356
0.2973
3
27
0.123288
0.1483
0.5471
0.1137
0.4110
4
25
0.114155
0.0926
0.6397
0.0953
0.5063
5
18
0.082192
0.0623
0.7021
0.0799
0.5861
6
7
0.031963
0.0445
0.7466
0.0669
0.6531
7
8
0.03653
0.0333
0.7799
0.0561
0.7092
8
7
0.031963
0.0258
0.8057
0.0470
0.7562
9
5
0.022831
0.0205
0.8262
0.0394
0.7957
10
3
0.013699
0.0167
0.8429
0.0331
0.8287
11
2
0.009132
0.0139
0.8568
0.0277
0.8564
12
0.0117
0.8685
0.0232
0.8796
13
2
0.009132
0.0009
0.8785
0.0195
0.8991
14
1
0.004566
0.0086
0.8871
0.0163
0.9154
15
1
0.004566
0.0075
0.8945
0.0137
0.9291
16
0.0066
0.9011
0.0115
0.9406
17
0.0058
0.9069
0.0096
0.9502
18
1
0.004566
0.0052
0.9122
0.0081
0.9582
19
0.0047
0.9168
0.0068
0.9650
20
0.0042
0.9211
0.0057
0.9707
21
1
0.004566
0.0038
0.9249
0.0047
0.9754
22
1
0.004566
0.0035
0.9283
0.0040
0.9794
23
0.0032
0.9315
0.0033
0.9827
24
0.0031
0.9344
0.0030
0.9855
25
0.0027
0.9371
0.0023
0.9879
26
0.0025
0.9396
0.0020
0.9898
27
1
0.004566
0.0023
0.9419
0.0016
0.9915
合计
219
1
0.9420
------
0.9915
------
4. 等待时间不同分布比较
再来考虑等待时间服从指数分布的情况,此时股票价格的跳跃是过程.等待时间分布的概率密度函数为:
. (16)
实证分析的过程如前文所述,回归结果见表2.同样考虑到概率密度函数的性质,我们可以得到式(17)所示的指数分布概率密度函数,并计算在经过天的等待后收盘点数发生跳跃的概率及累积概率,见表3.
(17)
由表3发现,本文所定义的分布而得到的概率存在一个缺点,即等待天数为1天的概率小于等待天数为2天的概率,这与频率统计结果不符,但等待天数在2天以内的累积概率达到0.3988,较指数分布更加符合频率统计结果.同样,由分布所得到的概率比指数分布得到的结果更加符合频率统计结果.而且,分布赋予尾部事件更高的概率:等待天数在27天以内的累积概率为94%,而指数分布却达到了99%.
指数收盘点数发生一次跳跃称为发生一次事件,将事件描绘在时间轴上可以得到一个事件序列.比较图2中真实发生的事件序列、过程和前文得到的更新过程所模拟出的事件序列,可以发现:本文所定义的更新过程与真实事件序列更为接近,都呈现出在较长时间的等待中出现密集的阵发事件的特点.我们从图1中也可以发现, 的概率密度函数较指数分布而言,为密集的阵发事件和不活跃的偶发事件都赋予了较高的概率,这一点更接近真实世界的情况.
a
b
c
图2 过程、真实过程、更新过程的事件序列
注:图中的横轴代表时间,每一个点代表在时刻发生的指数跳跃,每两个点之间的线段长度则代表等待时间的长短.过程所模拟出的100个连续事件的事件序列,本图使用的等待时间分布的概率密度函数为实证所得到的式(17);描绘出上证综合指数从2006年1月16日到2007年11月29日共451个交易日所发生的100次跳跃事件的序列;由本文所定义的更新过程产生的模拟100个连续事件的事件序列,每两个点之间的等待时间的长度服从本文所定义的分布,其概率密度函数为实证所得到的式(15).
当等待时间分布服从分布时,期望值的存在取决于参数的取值,即仅时存在期望值,否则不收敛。根据实证分析结果,等待天数的期望值大约为30.84天,因而允许股票价格跳跃过程出现较长的等待时间.同时,分布是具有时间记忆性的,在时刻等待天后股票价格出现跳跃的概率依赖于已经等待的天数.由图3,我们可以发现,当等待天数足够大的时候,函数呈现出幂律分布的特点,如式(18)所示.
(18)
由式(18)所存在的关系,我们可以得到如式(19)所示的等待天数的条件期望.
(19)
其中,为自上一次跳跃发生之后已经等待的天数,为继续等待的天数,即在一次跳跃发生天之后股票价格会再次发生跳跃.
同理,我们也可以得到股票价格发生波动的条件概率,如式(20)所示.
(20)
其中,为预期下一次跳跃距离此次跳跃的等待天数,即自一次跳跃发生天之后,再等待天股票价格将会发生新的跳跃.
由式(20),可以发现已知在第期股票价格发生跳跃之后,随着时间的推移再发生第二次跳跃的概率在逐渐减少,即在一次跳跃发生之后往往会紧跟着另一次跳跃.等待天数的概率存在幂率性是由等待时间概率密度函数的幂律性造成的,图4显示出等待天数的概率呈现出幂律性.
a b
图3 等待时间期望在双对数坐标系中的函数关系
注:图中实线代表函数,将其在上积分即可得到等待时间的数学期望,虚线代表函数.(a)在双对数坐标系中,函数在等待天数较小的阶段与函数有明显的不同,其斜率(绝对值)在为0之后随着等待天数的增加而逐渐增加;(b)当等待天数足够大的时候,函数在双对数坐标系中成为斜率为-1.064301的直线,即体现出幂律分布的特点.
a b
图4 分布下跳跃等待天数概率的分布散点图
注:(a)展示了由等待时间概率密度函数(式(15))计算得到的跳跃发生等待天数的概率分布,而(b)是在双对数坐标系中的情形.
2.2 跳跃幅度的概率分布
由于股票价格向上和向下发生跳跃的幅度是不对称的,因此我们需要对向上跳跃和向下跳跃进行分别讨论.从2006年1月16日到2009年6月30日上证综合指数收盘点数所发生的220次跳跃中,向上跳跃发生了111次,而向下跳跃发生了109次.
考虑式(10),我们将定义为,即对第天的收益率取自然对数,将定义为股票价格发生幅度为的跳跃的概率.将上证综合指数收盘点数的跳跃幅度∈[0.015,0.095]以0.0025为步长划分为28个区间,并按照跳跃方向的不同分别统计各个区间内跳跃发生次数的频率.
由统计结果,可以发现随着跳跃幅度的增加,跳跃发生的频率呈现出由小到大、再到小的变化过程.通过前文的实证分析,分布的概率密度函数也较好的拟合了真实世界所出现的情形.因此,识别跳跃的判断过程使得前半段呈现出上升趋势的分布概率密度函数在描述跳跃幅度中显得具有一定的合理性.
对式(9)进行回归分析,其中和分别为区间中位数和由统计得出的各个区间内跳跃发生次数的频率。回归结果由表4所示.由此,我们可以分别得到收盘点数向上跳跃和向下跳跃幅度的概率密度函数,其图形由图5所示.同时,考虑到概率密度函数的性质,我们可以得到收盘点数跳跃幅度的如下表达式:
(21)
表4 上证综合指数跳跃幅度的回归分析
向上跳跃
向下跳跃
系数估计
-12.03822
3.230466
0.060289
-15.35520
5.057079
0.155018
检验值
-5.164447
2.004270
1.414449
-6.593187
4.888499
4.181217
值
0.0001
0.0078
0.1743
0.0000
0.0001
0.0004
0.749303
0.581960
修正
0.721448
0.542147
统计值
26.89997
14.61721
a b
c d
图5 上证综合指数跳跃幅度的概率密度函数
注:图中虚线代表指数收盘点数向下跳跃的情形,实线代表指数收盘点数向上跳跃的情形.(a)收盘点数向上跳跃和向下跳跃幅度的概率密度函数;(b)将概率密度函数代表向下跳跃的左支和代表向上跳跃的右支重叠放置,可以更加清晰的分辨出两种情形的差异;(c)在双对数坐标中概率密度函数左支的图形,可以看出其为一条斜率为-5的直线,显示出幂律分布的特点;(d)在双对数坐标中概率密度函数右支的图形,可以看出其为一条斜率为-3的直线,也显示出幂律分布的特点,同时也可以看出左右两支的差异.
3 幂律跳跃扩散模型特征分析
对式(7)中的各个参数进行赋值,以得到股票日收益率的概率密度函数.首先,结合前文的实证结果,股票收益率每日发生跳跃的概率为,时间为1日(即0.004年),股票市场向盈利方向发生跳跃的概率为0.5,6个方程参数赋值为、、、、、,然后再假设无风险年收益率为5%,股票市场年波动率为25%,由此可以画出概率密度函数的图形(如图6所示).
a b
c d
图6 股票日收益率的概率密度函数
注:图中虚线代表股票收益率的概率密度函数,实线代表具有相同均值与方差的正态分布概率密度函数.(a)同时描绘了股票收益率的概率密度函数和正态分布的图形;(b)修正函数的图形,显示出对正态分布进行修正的幅度与方向;(c)显示出图形左尾的情形;(d)显示出图形右尾的情形,由此可以看到股票收益率所具有的胖尾特征.
下面进行250天(一个交易年)股票价格走势的数值模拟:
产生服从的随机数序列,根据表3确定服从分布的等待天数序列,由此确定250个交易日中股票价格发生跳跃的日期,并以相同的概率确定股票价格向上跳跃或者是向下跳跃.向上发生跳跃的幅度,由计算得出;向下发生跳跃的幅度,由计算得出.其中,为服从的随机数.考虑到涨跌停板制度的影响,假设股票价格每日涨跌幅度不能超过,因此的最大取值为0.09531,的最小取值为-0.10536,即赋予股票价格发生跳跃以至于涨跌停板的概率分别为0.094和0.058.没有发生跳跃过程的股票价格则由模型来生成,假设无风险年收益率始终为5%,第一周的波动率标准差为25%,之后每周的波动率方差则由上一周的收益率和方差共同确定,每日股票收益率就可以由服从正态分布的随机数模拟产生.
表5对模拟产生的股票收益率时间序列进行了描述性统计,可以看出股票收益率不服从正态分布,并具有左偏和尖峰胖尾的特征.图7给出了模拟得到的股票价格走势图(假设股票初始价格为100元)、收益率统计直方图和波动率图.
表5 模拟股票收益率时间序列的描述性统计
均 值
-0.000895916
中位数
0.000150
标准差
0.019627
方 差
0.00385205
最大值
0.095310
最小值
-0.086250
偏 度
-0.26238
峰 度
10.8003
Jarque-Bera统计量
636.6695
Jarque-Bera概率
0.000000
在数值模拟过程中,本文实际使用的股票价格收益率并不满足式(6),而是满足:
. (22)
这种处理方式与前文的实证过程相互一致,并且在实际应用中更加简便,不需要从股票价格跳跃幅度中剔除出波动率的影响.
a
b
c d
图7 股票价格走势数值模拟
注:(a)通过数值模拟得到的股票价格250个交易日的走势图;(b)在时间轴上描绘出模拟产生的在250个交易日中所发生的42次跳跃事件序列,时间轴上每个圆点代表发生一次跳跃;(c)通过数值模拟得到的股票收益率统计直方图,同时将其与正态分布进行比较,可以发现模拟得到的股票收益率具有尖峰胖尾的特征;(d)通过数值模拟得到的250个交易日股票价格时间序列的波动图.
当股票价格运动不考虑波动率的变化时,即股票价格变化如式(23)所示.
(23)
通过进行1000天(4个交易年)股票收益率的数值模拟可以发现股票价格数列存在波动聚集性.在理论上,由式(6)可以得到其方差,如式(24)所示.
(24)
可以看出,第期的波动率与第期的波动率相关性与两次跳跃之间的等待时间分布密切相关(见式(19)和(20)),即的计算受到的影响,波动聚集性的存在是状态依赖的.由此,更新过程的引入可以带来状态依赖的波动聚集性.
至此,本文所给出的幂律跳跃扩散模型可以同时得到股票价格收益率的尖峰胖尾特征和波动聚集性特征.
4 幂律跳跃扩散模型的应用
在金融衍生品定价方面,由于幂律跳跃扩散模型全面地反映了基础标的物价格运行特征,因此,该模型在无套利假设下可以更准确地为衍生品定价。
对于给定的概率,定义
, (25)
其中,的概率密度函数如式(3)所示,为本文所定义的更新过程.欧式看涨期权的定价公式可以表示为
, (26)
其中,为标的资产股票的价格,为欧式看涨期权执行价格,为无风险收益率,.
漂移幂律(Shifted Power Law)是一种介于幂律与指数函数之间的分布[29].对于,当,时,带有幂律性的函数就会转化为典型的指数函数.因此,本文所定义的更新过程可以退化为过程,跳跃幅度的幂律性分布也可以退化为指数分布,由此本文所定义的模型可以退化到Kou所提出的指数跳跃扩散模型,因此幂律跳跃扩散模型是指数跳跃扩散模型的更一般形式.当股票价格跳跃幅度变得很小或更新过程的强度为0时,本文所提出的幂律跳跃扩散模型则退化为模型.
在风险管理方面,更新过程的引入除了可以更加准确地描述股票价格的运动过程,而且在实际应用领域对金融风险管理也是具有一定参考意义的.如果将过程作为描述股票价格的跳跃过程,实际上隐含着默认在任何一个时刻股票价格发生跳跃的概率都是相同的,显然这与实际情况并不相符.而本文所定义的更新过程则可以模拟得到股票价格在较长时间不发生跳跃过程的情形,也可以得到跳跃过程在短时间内密集发生的状况,这更加符合真实的情况.更新过程的引入使得股票收益率动力学方程可以描绘股票价格在任何一个时点发生跳跃的概率并不相同,而这种不相同的概率赋予是依赖于当时市场所处的状态的——这表明市场存在一定的记忆性,即今天市场的现状并非凭空而来,而是基于其历史的,而这与现实世界的情形更为贴切.引入更新过程以后获得的跳跃发生的时间期望以及某一时间跳跃发生的概率可以成为一个运用于实际市场操作领域的指标,而且该指标来源于市场本身,获取和实际应用起来简便易行。从实际操作角度来说,基于股票价格收益率服从正态分布的显得对极端事件无可奈何,因而需要引入其他指标对进行修正,而跳跃扩散模型则可以较好地解决这个问题,同时更新过程可以得到量化的结果,因此对金融风险管理的实际操作具有一定的现实指导作用.
5 结论
本文在Merton的跳跃扩散模型的思想基础上进行了两个方面的拓展,提出了幂律跳跃扩散模型.首先,将计数过程由过程修正为带有幂律性质的更新过程,可以描述股票价格运动过程中密集发生的跳跃,同时也允许在较长时间内不出现跳跃的情形;然后,在考虑收益损失状态下跳跃幅度不同的同时,也赋予跳跃幅度以幂律分布的特征.通过实证研究发现幂律跳跃扩散模型所作出的两个修正可以更加准确的描述股票价格的运动过程,同时得到具有胖尾特征的收益率分布和波动聚集性,可以更细致地刻画股票价格运动,为金融衍生资产定价奠定基础;同时,引入更新过程后也为金融风险管理提供了两个可操作的指标。
参考文献:
[1] Bachelier L. Théorie de la Spéculation, Annales de I'Ecole Normale Supérieure,3, Paris:Gauthier-Villars,1900.
[2] Fama, E. F. Efficient Capital Markets: A Review of the Theory and Empirical Work,” Journal of Finance, 1971, 25(2):383-417.
[3] Black F.,M.Scholes . The Pricing of Options and Corporate Liabilities. Journal of Political Economy, 1973,81(3):637-654.
[4] Osborne, M. F. M. (1964), “Brownian Motion in the Stock Market,” in P. Coonter, ed., The Random Character of Stock Market Prices, Cambridge, MA: M.I.T. Press.
[5] Engle R.F.. Autoregressive Conditional Heteroskedasticity with Estimates of the Variance of U.K. Inflation. Econometrica, 1982, 50:987-1008.
[6] Bollerslev, Tim. Generalized Autoregressive Conditional Heteroskedasticity. Journal of Econometrics,1986,31:307-327.
[7] Cox J.C.,S.A.Ross.The Valuation of Options for Alternative Stochastic Processes. Journal of Financial Economics,1976,3:145-166.
[8] Hull J.C.,A. White. The Pricing of Options on Assets with Stochastic Volatilities. Journal of Financial Studies, 1987,4(4):727-757.
[9] Heston S.L. A Closed Form Solution for Option with Stochastic Volatility with Application to Bond and Currencies. Review of Financial Studies, 1993,6(3):327-343.
[10] Peters, E.E. Fractal Market Analysis: Applying Chaos Theory to Investment and Economics. New York: John Wiley and Sons, 1994.
[11] Rogers, L.C.G. Arbitrage from Fractional Brownian Motion. Math Finance, 1997, 7:95-105.
[12] Merton R.C. Option Pricing when Underlying Stock Returns are Discontinuous. Journal of Financial Economics, 1976,3(1-2): 125-144.
[13] Kou S.G. A Jump-Diffusion Model for Option Pricing. Management Science,2002, Vol.48,N.O. 8:1086-1101
[14] Bates D. Jumps and Stochastic Volatility: Exchange Rates Processes Implicit in Deutsche Mark Options. Review of Financial Studies, 1996,9(1):69-107.
[15] Scott L.O. Pricing Stock Option in a Jump-diffusion Model with Stochastic Volatility and Interest Rates: Application of Fourier Inversion Methods. Mathemtical Finance, 1997,7(4):413-423.
[16] Duffie D., J.Pan,
展开阅读全文