期权定价中的蒙特卡洛模拟方法

  • 格式:doc
  • 大小:1.61 MB
  • 文档页数:55

实用标准文档 文案大全 期权定价中的蒙特卡洛模拟方法

期权作为最基础的金融衍生产品之一,为其定价一直是金融工程的重要研究领域,主要使用的定价方法有偏微分方程法、鞅方法和数值方法。而数值方法又包括了二叉树方法、有限差分法和蒙特卡洛模拟方法。 蒙特卡洛方法的理论基础是概率论与数理统计,其实质是通过模拟标的资产价格路径预测期权的平均回报并得到期权价格估计值。蒙特卡洛方法的最大优势是误差收敛率不依赖于问题的维数,从而非常适宜为高维期权定价。

§1. 预备知识 ◆两个重要的定理:柯尔莫哥洛夫(Kolmogorov)强大数定律和莱维一林德贝格(Levy-Lindeberg)中心极限定理。

大数定律是概率论中用以说明大量随机现象平均结果实用标准文档

文案大全 稳定性的一系列极限定律。在蒙特卡洛方法中用到的是随机变量序列同分布的Kolmogorov强大数定律: 设12,,L为独立同分布的随机变量序列,若

[],1,2,kEkL则有11(lim)1nknkpn 显然,若12,,,nL

是由同一总体中得到的抽样,那么由

此大数定律可知样本均值11nkkn当n很大时以概率1收敛于

总体均值。

中心极限定理是研究随机变量之和的极限分布在何种情形下是正态的,并由此应用正态分布的良好性质解决实际问题。 设12,,L为独立同分布的随机变量序列,若

2[],[],1,2,kkEDkL

则有

1(0,1)nkdknNn



其等价形式为2111lim()exp(),22nxkkntnPxdtxn。

◆Black-Scholes期权定价模型 模型的假设条件: 1、标的证券的价格遵循几何布朗运动 实用标准文档 文案大全 dSdtdWS

其中,标的资产的价格S是时间t的函数,为标的资产的瞬时期望收益率,为标的资产的波动率,dW是维纳过程。 2、证券允许卖空、证券交易连续和证券高度可分。 3、不考虑交易费用或税收等交易成本。 4、在衍生证券的存续期不支付红利。 5、市场上不存在无风险的套利机会。 6、无风险利率r为一个固定的常数。

下面,通过构造标的资产与期权的资产组合并根据无套利定价原理建立期权定价模型。首先,为了得到期权的微分形式,先介绍随机微积分中的最重要的伊藤公式。 伊藤Ito公式:设(,)VVSt,V是二元可微函数,若随机过程S满足如下的随机微分方程 (,)(,)dSStdtStdWS 则有 222

2

1((,)(,))(,)2VVVVdVStSStSdtStSdWtSSS



根据伊藤公式,当标的资产的运动规律服从假设条件中的几何布朗运动时,期权的价值(,)VVSt的微分形式为 实用标准文档 文案大全 222

2

1()2VVVVdVSSdtSdWtSSS



现在构造无风险资产组合VVSS,即有drdt,

经整理后得到 222

2

102VVVSrSrVtSS



这个表达式就是表示期权价格变化的Black-Scholes

偏微分方程。它同时适合欧式看涨期权、欧式看跌期权、美式看涨期权和美式看跌期权,只是它们的终值条件和边界条件不同,其价值也不相同。 欧式看涨期权的终边值条件分别为 (,)max0,TVSTSK

,00(,)SVSTSS

通过求解带有终边值条件的偏微分方程,得出欧式看涨期权的的解析解: ()12(,)()()rTtVStSNdKeNd

其中,221()2xdNdedx,21ln(/)(/2)()SKrTtdTt,21ddTt,T为期权的执行日期,K为期权的执行价格。

欧式看跌期权的终边值条件分别为 (,)max0,TVSTKS

,0(,)0KSVSTS 实用标准文档 文案大全 此外,美式看涨期权的终值条件为(,)max{0,}VStSK,美式看跌期权的终值条件为(,)max{0,}VStKS。然而,美式期权的价值没有解析解,我们一般可通过数值方法(蒙特卡洛模拟、有限差分法等)求得其近似解。

◆风险中性期权定价模型 如果期权的标的资产价格服从几何布朗运动 dSrdtdWS

即标的资产的瞬时期望收益率取为无风险利率r。同理,根据伊藤公式可以得到 2ln()2dSrdtdW

222lnln()()()~(()(),())22TtTtSSrTtWWNrTtTt

2exp(()()())2TtTtSSrTtWW

对数正态分布的概率密度函数:设2~(,)N,e,则的密度函数为 221(ln)exp()0()2200xxPxxx







根据上述公式,得到标的资产TS的密度函数如下 实用标准文档 文案大全 22

2(ln()())21exp()0()2()200txrTtSxPxTtTtxx







在风险中性概率测度下,欧式看涨期权定价为: (,)exp(())[max{0,}]QTVStrTtESK



22

222

2

(ln()())12[max{0,}]exp()22(ln()())2exp()22QTKKxrTtSESKdxTtTtxrTtKSdxTtxTt





接下来,求解以上风险中性期望。首先,对上式的右边第一个广义积分分别作变量替换 2ln()()2xrTtSyTt

和uyTt,可以得到

222

2

22

2ln()()2()()()221ln()()2(ln()())12exp()2211()22KSrTtuuKrTtrTtrTtTtKrTtSTtxrTtSdxTtTtSeeduSeeduSeNd









 再对等式的右边的第二个无穷积分,令 2lnln()()2xSrTtuTt

,可求得 实用标准文档 文案大全 222

2

22

2lnln()()2222lnln()()2(ln()())2exp()2()211()22KSKrTtuuTtKSrTtTtxrTtKSdxTtxTtKeduKeduKNd









将以上的计算结果代入期望等式中,得到欧式看涨期权的价格公式为:

()()12(,)[max{0,}]()()rTtQrTtTVSteESKSNdKeNd

其中,21ln()()2SrTtKdTt,21ddTt。

可以看出,对于欧式看涨期权的风险中性定价方法的结果与基于资产复制的偏微分方程定价方法的结果是一致的。基于风险中性的期权定价原理在于:任何资产在风险中性概率测度下,对于持有者来说都是风险偏好中性的,便可用风险中性概率求取期权的期望回报再将其进行无风险折现便是初始时刻的期权价值。蒙特卡洛模拟方法就是一种基于风险中性原理的期权数值定价方法。

§2. 蒙特卡洛模拟方法及其效率