广义维纳过程

Download Report

Transcript 广义维纳过程

投资学专题7:
金融衍生产品定价
复旦大学金融研究院 张宗新
outline
资产价格运动的随机过程
 二叉树模型及其在衍生产品定价中的应用
 Black-Scholes 期权定价模型在衍生产品
定价中的应用
 Monte Carlo模拟在衍生产品定价中的应用


第一节 资产价格运动的随机过程

金融资产价格的运动随时间变化,形成一个随机过程。随
机过程是用来描述随机变量随着时间变化的统计术语。观
测到的价格是随机过程的一个实现,随机过程的理论是对
观测到的价格进行分析和作出统计推断的基础。
资产价格波动的随机过程

一、Wiener过程或Brownian运动

1、维纳过程(Wiener Processes)
股价行为模型通常用维纳过程来表达。理解遵循维纳过程的变量z的
行为,可以考虑在小时间间隔上变量z值的变化。设一个小的时间间
隔长度为△t,定义 z 为在△t时间内z的变化。要使z遵循维纳过程,
△z必须满足两个基本性质:
性质1:△z与△t的关系满足方程式


z   t







其中为从标准正态分布N(0,1)中抽取的一个随机值。
性质2:对于任何两个不同时间间隔△t,△z的值相互独立。
2、广义维纳过程(Generalized Wiener Process)
变量x的广义维纳过程用dz定义如下:
dx=adt+b dz
其中a和b为常数。理解方程较好的方法是分别考虑方程右边的两个组
成部分。adt项说明了x变量单位时间的漂移率期望值为a。如果缺省
bdz项,方程变为dx =adt dx/dt=a
即 x= x0+ at
维纳过程
内幕交易
概率预测
模型 示意
图
《管理世
界》
2008.4

二、Ito(伊藤)引理

一般维纳过程的漂移参数和波动率参数都是不随时间变化的。如果进
一步扩展模型,允许 和  是随机过程 xt的函数,那么我们就可以引
出一个伊藤过程。
伊藤过程,是指如下随机过程:
dx   ( x, t )dt   ( x, t ) dt


其中, dz   dt 是一个标准布朗运动, 、 是变量x和t的函数。






为表述伊藤引理,将资产的随机过程表述为如下方程:
dS   Sdt   Sdz
也就是说,用漂移率  S 和波动率 S 的伊藤过程表示资产
价格的动态。
在时间间隔为 t 后,资产价格的变化比率为:
S
 t   t
S
可见,S / S也具有正态分布特征,其均值为t ,标准差
为  t ,方差为  2 t 。
资产价格运动的随机过程



三、漂移参数  和波动率参数 的估计
上述方程的几何布朗运动中有两个未知参数 和 可以用经验方法进
行估计。假定我们有股价Pt 在等时间间隔 t 上的个观测值,观测到的
股价序列 P1, P2 , , Pn  ,其中 t  1, 2, , n 。
令 rt  ln(Pt )  ln(Pt 1 ),存在 Pt  Pt 1 exp(rt ) ,其中 rt 为第t个时间间隔上的连
续复合收益率。根据Ito引理,并且假定股价 Pt 服从一个几何布朗运动,
2
我们得到 rt 服从均值为(   / 2)t ,方差为 2 t 的正态分布。
第二节二叉树模型及其在衍生产品定价中的应用

一、二叉树模型(Binomial Tree Model)

二叉树期权定价模型假定,在每一期股票价格可以沿两个方向—
—向上或向下——中的任何一个方向变动。因此,可以将将时间T
分为很多小的时间间隔 ,在一个时间间隔内证券价格价格只有两
种运动可能:从开始的S上升到原来的u倍,即Su;或下降到原来
的d倍,其中u>1,d<1(一般假定 u  1 )。也就是说,股价上升或
d
t
下降分别用u和d表示,而在每一个
t ,股票价格变化由S到Su或Sd.
若价格上扬的概率为p, 那么下跌的概率为q=1- p。
p
uS0
S0
1-p
dS0
二叉树模型及其在衍生产品定价中
的应用




二、二叉树期权定价模型
二项式期权定价模型(Binomial Option Pricing Model,简
称BOPM)是对期权进行估价方法,它是通过统计中的二
项分布,假定只有两种可能结果而推算出来的。
下面,我们可以分六步骤对看涨期权的二项式期权定价模
型进行分析:
第一步:分析股价的未来可能运动形态;第二步:列出期
权的价格分布;第三步:构建对冲投资组合;第四步:对保
值比率进行求解;第五步:用净现值法(NPV)解出买入
期权的价格;第六步,将单期扩展之多期。
假定不支付红利股票的的3个月期的美式看跌期权,股票价格15元,
执行价格15元,无风险利率为3%,波动率为50%。即:S=15,
X=15,r=0.03, =0.5,T=0.25。
为构造二叉树,假定到期期限分为4个阶段,每段长度
=0.25/4=0.0625。
二叉树模型及其在衍生产品定价中
的应用



三、二叉树模型在可转债定价中的应用
可转债的二叉树定价步骤如下:
第一步,先计算出对应股票的二叉树节点上的数值。利用股价的历史
数据(一般利用过去3个月或者半年的股价数据)估计出股票的波动
率  ,然后计算出二叉树的几个重要参数。
r t
u  e t , d  1/ u, p  (e
 d ) (u  d )
rf

t

(
T

t
)
m
其中
,t,T分别指的是可转债的初始和期末时刻, 为无风
险利率,使用这些参数就可以推出股票的价格树图。
第二步,通过可转债的相关条件来递推价格树中各个节点的可转债的
价格。
f


二叉树模型及其在衍生产品定价中
的应用
105
实际价格
理论价格
104
103
102
101
100
99
98
-1
-8
13
20
27
10
15
22
29
-6
-3
9
9
0
1
10
10
10
11
-9
-9
-9
-1
-1
05
05
5
5
5
5
5
0
0
0
0
0
05
05
05
05
应用二叉树模型,Matlab程序对西钢转债进行拟合。
第三节 Black-Scholes期权定价模型
在衍生产品



一、 Black-Scholes 期权定价模型
1973年,美国芝加哥大学教授费希尔·布莱克和迈伦·斯科
尔斯提出了有史以来的第一个期权定价模型,即布莱克斯科尔期权定价模型(Black-Scholes Options Pricing
Model,BSOPM)。
布莱克-斯科尔斯推导出了一个确定期权价格的明确公式,
 r (T t )
即:
c  SN(d1 )  Xe

其中
N (d 2 )
1
ln( S / X )  (r   2 )(T  t )
2
d1 
 T t
1
ln( S / X )  (r   2 )(T  t )
2
d2 
 d1   T  t
 T t

二、B-S期权定价求解

由于BS公式是关于期权定价的连续时间公式,因此容易分析期权价
格的敏感性,即可以利用BS公式求出的看涨期权的价格同看涨期权
的内在价值进行比较分析,分析两者随着股票价格变化的差异。
12
期权价值
10
8
6
4
2
0
1
2
3
4
5
6
7
8
9
10
11
看涨期权的价格和内在价值
资产价格
BS公式求出的看涨-看跌期权价格
BS公式的EXCEL求解过程。在此,股票价格S为25元,执行价格X为25
元,无风险利率为8%,股票的波动率为30%,期权的到期年限为0.5年
。计算相应的看涨期权的价格。
运用B-S公式进行期权定价求解
当前股价
25
执行价格
25
无风险利率
8%
到期时间(年)T
0.5
股价波动率
30%
股票价格
看涨期权价格
内在价值
2.597032043
0
10
1.29408E-05
0
12.5
0.001053997
0
15
0.018528597
0
d1
0.294627825
17.5
0.129037806
0
d2
0.082495791
20
0.49746469
0
0.615860834
22.5
1.297035624
0
N(d1)
N(d2)
0.532873834
25
2.597032043
0
27.5
4.341892958
2.5
30
6.411299266
5
32.5
8.684918903
7.5
35
11.07341867
10
看涨期权价格
看跌期权价格(利用平
权)
看跌期权价格(利用
BS公式)
2.597032043 S*N(-d1)-X*exp(-r*T)*N(d2)
1.616768021 C-S+X*exp(-r*T)
1.616768021 X*exp(-r*T)*N(-d2)-S*N(d1)
20
05
-1
20 205 05
-1
20 206 21
20 0106 10
-0
20 206 06
-0
20 206 22
20 0306 13
-0
20 306 29
20 0406 14
-0
20 506 09
-0
20 506 25
20 0606 12
-0
20 606 28
-0
20 706 14
20 0806 01
-0
20 806 17
-0
20 906 04
20 0906 20
-1
20 006 13
-1
031
Black-Scholes期权定价模型在衍生
产品
6
鞍钢JTC1 - 实际价格
5
鞍钢JTC1 - 理论价值
4
3
2
1
0
鞍钢权证的理论价格和实际价格
三、波动率与波动率微笑
1、历史波动率


对于理想的欧式期权而言,BS期权定价模型仅依赖于五
个参数:股票价格、期权的执行价格、期权的到期时间、
无风险利率和股票的价格波动率。在这些参数中,和由发
行的金融合约的条款所定,和可从市场得到。唯一需要确
定的参数就是波动率。
请注意,BS模型中波动率是指在t  0 到 t  T 的未来时期
内的标的资产的波动率。由于在现实金融市场上,证券价
格的波动是一个随机过程,估计波动率并不是一件简单的
事情。通常,有两种方法可以对波动率进行估计,即历史
波动率(historical volatility)与隐含波动(imp volatility)。


(1)方差估计法
计算方式如下:先计算出标的资产价格S第i天的报酬u ,
即ui =ln(Si/Si-1),利用此前一段时间(可选择3个月、半年)
资产报酬数据,估计日报酬的标准差。即:
t

1 m
 
(uni  u )2

m  1 i 1
2
u

m
1
 un  i
m i 1

1 m
(u n  i  u ) 2

m  1 i 1
这里,
为 u i的算术平均。u i 的标准差相当于  的
估计值,其中  为时间间隔长度(以年为计算单位)。
鞍钢股份的历史波动率
2006-11-23
2006-11-08
2006-10-24
2006-10-09
2006-09-15
2006-08-30
2006-07-11
2006-06-01
2006-05-11
2006-04-19
2006-04-04
2006-03-14
2006-02-22
2006-02-07
2006-01-11
2005-12-20
2005-12-05
2005-11-18
2005-11-03
2005-09-07
2005-08-23
2005-08-08
2005-06-24
2005-06-09
2005-02-25
2005-02-01
0.14
Volatility
0.12
0.1
0.08
0.06
0.04
0.02
0
鞍钢股份的波动性(GARCH估计)
2006-11-23
2006-11-08
2006-10-24
2006-10-09
2006-09-15
2006-08-30
2006-07-11
2006-06-01
2006-05-11
2006-04-19
2006-04-04
2006-03-14
2006-02-22
2006-02-07
2006-01-11
2005-12-20
2005-12-05
2005-11-18
2005-11-03
2005-09-07
2005-08-23
2005-08-08
2005-06-24
2005-06-09
2005-02-25
2005-02-01

(2)GARCH(1,1)模型估计
volatility
0.003
0.0025
0.002
0.0015
0.001
0.0005
0

2、隐含波动率
确定波动率的第二种方法是估计隐含波动率。隐含波动率是另外
一种定义,假定:
 S0 为当前股票价格;K为执行价格;T为到期时间;r为无风险利率;
V为期权当前的市场价格。
 利用上述参数,通过数值方法求解下式,可以得到隐含波动率的
V (S , T  t )  SN (d1 )  Ker (T t ) N (d2 )
值:
 其中,时间从到期日起以天计,且:

ln(S / K )  (r   I 2 / 2)(T  t )
d1 
I T t
d2  d1   I T  t
华菱权证的历史波动和隐含波动
2日
9日
0.5
11
月
12
日
14
日
9月
17
日
8月
20
日
7月
22
日
6月
25
日
5月
27
日
4月
30
日
3月
10
月
20
06
年
20
06
年
20
06
年
20
06
年
20
06
年
20
06
年
20
06
年
20
06
年
20
06
年
3月
20
06
年
0.6
隐含波动率
波动率
0.4
0.3
0.2
0.1
0

3、波动率微笑(Volatility Smiles)

应用期权市场价格和BS公式推算出来的隐含波动率具有以下两个方
向的变动规律:
一是“波动率微笑”,即隐含波动率会随着期权执行价格不同而不同。
由于隐含波动率是执行价格和到期日的函数,特别地,当执行价格等
于股票最初价格S0时,隐含波动率最小,当执行价格偏离0时,隐含
波动率会增加,这种现象通常称为“波动率微笑”。
二是波动率期限结构(Volatility Term Structure),即隐含波动率会
随期权到期时间不同而变化。从长期来看,波动率大多表现出均值回
归,即到期日接近时,隐含波动率的变化较剧烈,随着到期时间的延
长,隐含波动率将逐渐向历史波动率的平均值靠近。波动率微笑的形
状也受到期权到期时间的影响。一般而言,期权到期日越近,波动率
“微笑”就越显著,到期日越长,不同价格的隐含波动率差异越小,
接近于常数。



四、期权的衍生物及其风险对冲

1.德尔塔( )

在任何确定的时间内,衍生证券的价值是标的资产价格的函数。这个
函数对标的资产价格变化的敏感度用希腊字母德尔塔(Delta,)来
描述。德尔塔( )是Black-Schols期权定价模型的一个重要衍生概
f
念,在证券组合中对投资者具有重要意义。其公式表达为:

14

S
12
10
C
8
6
4
2
0
10
C:\C_doc\教学\Matlab_金融试验\衍生
市场
12
14
16
18
20
S
22
24
26
28
30
看涨期权对股票价格的敏感性和Delta策略




在Black-Schols期权定价模型中,德尔塔( )的决定十
分简单:它就等于N (d1 ) 。德尔塔特性如下:
(1)认购权证的Delta一定为正值,认沽权证的Delta一定
为负值。这正负号表示期权价格和标的资产价格之间的变
动关系。正号表示同向变动,负号表示异向变动;
(2)Delta数值的范围介于-1和+1之间。
(3)平价期权的Delta数值约为0.5。
Black-Scholes期权定价模型在衍生
产品



Delta中性组合
对于价格低于理论价值的权证,还可以按比例购买
股票+认沽权证,构造Delta 中性组合,也就是买入
波动率,在股价向任何一方向变动时,组合价值都
将上升。
但是,如果波动率下降,则会影响到套期保值效
果,隐含波动率变动会使得组合价值曲线发生位
移,一旦权证的隐含波动率下降,即使股价发生
了较大变动,组合价值仍会受到损失。
Black-Scholes期权定价模型在衍生
产品




2.伽马()
仅在标的股票的价格只发生微小变动时,德尔塔对冲才是
有效的,因为它只考虑了一阶导数。如果标的股票价格可
能发生较大的变化,那么,对冲组合就要考虑二阶导数。
于是,引入伽马(Gamma, )的概念。
伽马度量的是衍生资产的凸性,伽马度量的是期权价格曲
线上该点的二阶导数。对于不支付红利的欧式期权来说,
存在:
2

c 
N ( d1 )
  C
 2 
S S
S T  t
X=30;sigma=0.3; r_f=0.05
价格=15-50
周期:12M
0.14
0.12
Gamma
0.1
0.08
0.06
0.04
0.02
12
10
50
8
40
6
30
4
2
20
2010
0.1
0.2
Delta
0.3
0.4
0.5
0.6
0.7
0.8
0.9
C:\Program Files\MATLAB71\work\Delta_Gamma_option01
Black-Scholes期权定价模型在衍生
产品




3.西塔( )
西塔( ,Theta)是期权定价中的另一个重要参数。西塔( )
被定义为:

f
1
 re rT XN (d 2 )   2 S 2
t
2
西塔度量的是衍生证券价值的变动方向。如果时间增加,
期权曲线将向右移动。西塔正是度量的曲线的这种移动。
Black-Scholes期权定价模型在衍生
产品



4、维加(vega,)
当波动率变化一个单位时(通常为1%),衍生证券的价
值变化称为维加(vega, )。用公式表达为:

f

反映的是证券价格本身波动对衍生证券价格的影响。若
构造的组合使  值等于零,则该组合的价值不受波动率变
化的影响。按照BS期权定价公式,可以得到不支付红利
股票的欧式看涨期权和看跌期权的  表达式:
 
C 
d
C
 SN (d1 ) 1  Xe  r (T t ) N (d 2 )  S T  tz (d1 )


P 
P S



( Xe  r (T t )  S )   C
  
Black-Scholes期权定价模型在衍生
产品




5、罗( )
当利率变化一个单位时(通常为1%),衍生证券的价值
变化称为罗( )。用公式表达为:

f
r
可见,反映的是衍生产品价格对利率变化的比率。按照
BS期权定价公式,可以得到不支付红利股票的欧式看涨
期权和看跌期权的  表达式:
C 
C
 Xe  r (T t ) N (d 2 )
r
P 
P
  X (T  t )e  r (T t ) N (d 2 )
r
蒙特卡罗模拟在衍生产品定价中的
应用

第四节 蒙特卡罗模拟在衍生产品定价中的应用

一、蒙特卡罗模拟方法介绍

蒙特卡罗模拟(Monte Carlo Simulation)是一种通过模拟标的资产
价格随机运动路径得到权证价值期望值的数值方法,是一种应用十分
广泛的金融衍生产品定价方法。
如果股价运动服从伊藤过程,则当然股价如果服从其他分布,只要给
出具体的表达式,我们就可利用蒙特卡罗模拟法进行模拟。蒙特卡罗
模拟进行期权定价的核心在于生成股价价格的随机过程。

蒙特卡罗模拟在衍生产品定价中的
应用







蒙特卡罗模拟的实质是模拟标的变量的随机运动,预测其衍生产品的
平均回报,并由此得到衍生品价格的一个概率解。
蒙特卡罗模拟的优点:
(1)提供一个相当广泛和强大的期权定价技术:
(2)得到广泛应用。正如Campell所指出的:“虽然在依赖于路径的
的衍生证券定价的确方面存在一些近似的解析解,但最为有效的方法
还是蒙特卡罗模拟。”
蒙特卡罗模拟的缺点:
(1)只能用于欧式期权,期权不能提前执行。美式期权的提前执行
性增加了期权定价的复杂性,而这必须用一个动态程序分析进行递归。
因而尝试应用蒙特卡罗模拟技术来对美式期权定价,成为近年来者来
领域的发展方向之一。
(2)为了达到一定的精度,必须进行大量的模拟运算。
蒙特卡罗模拟在衍生产品定价中的
应用

二、金融衍生产品定价的Monte Carlo模拟





Boyle最早利用Monte Carlo方法进行期权定价。Monte Carlo模拟
相关步骤为:
第1步:设定基础资产数据的生成过程。通常使用带漂移的随机游
走模型,同时指定漂移项的大小和波动参数的大小,设定执行加
价格K和期限T。
第2步:从正态分布中抽取长度为T的序列,该序列为误差项序列,
 t ~ N。
(0,1)
即
第3步:构建观测值长度为T的基础资产序列。
第4步:分析基础资产在到期日T时的价格。对看涨期权而言,若
到期日基础资产价格 PT  K ,则此次重复试验下,期权价值为0并
到期作废。
蒙特卡洛模拟计算欧式期权价格
function eucall=blsmc(s0,K,r,T,sigma,Nu)
 s0=50;K=52;r=0.05;T=5/12;sigma=0.4;Nu
=1000;


MoterCarlo_01.M
鞍钢权证市场价格、BS理论价格、蒙特卡罗拟合图
2006-10-23
2006-10-9
2006-9-25
5
2006-9-11
2006-8-28
2006-8-14
2006-7-31
2006-7-17
2006-7-3
2006-6-19
2006-6-5
2006-5-22
2006-5-8
2006-4-24
2006-4-10
2006-3-27
2006-3-13
2006-2-27
2006-2-13
2006-1-30
2006-1-16
2006-1-2
2005-12-19
2005-12-5
6
BS理论价格
实际价格
蒙特卡洛模拟价格
4
3
2
1
0