《期权定价中的蒙特卡洛模拟方法new.doc

期权定价中的蒙特卡洛模拟方法 期权作为最基础的金融衍生产品之一，为其定价一直是金融工程的重要研究领域，主要使用的定价方法有偏微分方程法、鞅方法和数值方法。而数值方法又包括了二叉树方法、有限差分法和蒙特卡洛模拟方法。 蒙特卡洛方法的理论基础是概率论与数理统计，其实质是通过模拟标的资产价格路径预测期权的平均回报并得到期权价格估计值。蒙特卡洛方法的最大优势是误差收敛率不依赖于问题的维数，从而非常适宜为高维期权定价。 §1. 预备知识 ◆两个重要的定理：柯尔莫哥洛夫(Kolmogorov)强大数定律和莱维一林德贝格(Levy-Lindeberg)中心极限定理。 大数定律是概率论中用以说明大量随机现象平均结果稳定性的一系列极限定律。在蒙特卡洛方法中用到的是随机变量序列同分布的Kolmogorov强大数定律： 设为独立同分布的随机变量序列，若 则有 显然，若是由同一总体中得到的抽样，那么由此大数定律可知样本均值当n很大时以概率1收敛于总体均值。 中心极限定理是研究随机变量之和的极限分布在何种情形下是正态的，并由此应用正态分布的良好性质解决实际问题。 设为独立同分布的随机变量序列，若 则有 其等价形式为。 ◆Black-Scholes期权定价模型 模型的假设条件： 1、价格遵循几何布朗运动 其中，标的资产的价格是的函数，为标的资产的期望率，为标的资产的波动率，是2、证券允许卖空、证券交易连续和证券高度可分。 3交易费用或税收。 4在衍生证券的支付红利。 5市场上不存在套利机会。 6无风险利率为一个固定的常数。，是二元可微函数，若随机过程满足如下的随机微分方程 则有 根据伊藤公式，当标的资产的运动规律时，期权的价值 现在构造无风险资产组合，，经整理后得到 这个表达式就是表示期权价格变化的Black-Scholes微分方程。， 通过求解带有终边值条件的偏微分方程，得出欧式看涨期权的的解析解： 其中，，为期权的执行日期，为期权的执行价格。， 此外，美式看涨期权的终值条件为，美式看跌期权的终值条件为。然而，美式期权的价值没有解析解，我们一般可通过数值方法（蒙特卡洛模拟、有限差分法等）求得其近似解。 ◆风险中性期权定价模型 如果期权的标的资产价格服从几何布朗运动 即标的资产的瞬时期望率。同理，根据伊藤公式可以得到 对数正态分布的概率密度函数：设，，则的密度函数为 根据上述公式，得到标的资产的密度函数如下 在风险中性概率测度下，欧式看涨期权定价为： 接下来，求解以上风险中性期望。首先，对上式的右边第一个广义积分分别作变量替换 和，可以得到 再对等式的右边的第二个无穷积分，令 ，可求得 将以上的计算结果代入期望等式中，得到欧式看涨期权的价格公式为： 其中，，。 可以看出，对于欧式看涨期权的风险中性定价方法的结果与基于资产复制的偏微分方程定价方法的结果是一致的。基于风险中性的期权定价原理在于：任何资产在风险中性概率测度下，对于持有者来说都是风险偏好中性的，便可用风险中性概率求取期权的期望回报再将其进行无风险折现便是初始时刻的期权价值。蒙特卡洛模拟方法就是一种基于风险中性原理的期权数值定价方法。 §2. 蒙特卡洛模拟方法及其效率 假设所求量是随机变量的数学期望，那么近似确定的蒙特卡洛方法是对进行n次重复抽样，产生独立同分布的随机变量序列，并计算样本均值 。那么根据Kolmogorov强大数定律有 。因此，当n充分大时，可用作为所求量的估计值。 由中心极限定理可得到估计的误差。设随机变量的方差，对于标准正态分布的上分位数，有 这表明，置信水平对应的渐近置信区间是 。实际上，由此可确定蒙特卡洛方法的概率化误差边界，其误差为，误差收敛速度是。 不难看出，蒙特卡洛方法的误差是由和决定的。在对同一个进行抽样的前提下，若想将精度提高一位数字，要么固定，将n增大100倍；要么固定n将减小10倍。若两个随机变量的数学期望，，那么无论从或中抽样均可得到的蒙特卡洛估计值。比较其误差，设获得的一个抽样所需的机时为，那么在时间T内生成的抽样数，若使，则需使。因而，若要提高蒙特卡罗方法的效率，不能单纯考虑增加模拟的次数n或是减小方差，应当在减小方差的同时兼顾抽取一个样本所耗费的机时，使方差与机时t的乘积尽量的小。 §3. 蒙特卡洛模拟方法为期权定价的实现步骤 期权定价的蒙特卡洛方法的理论依据是风险中性定价原理：在风险中性测度下，期权价格能够表示为其到期回报的贴现的期望值，即，其中的表示风险中性期望，r为无风险利率，T为期权的到期执行时刻，是关于标的资产价格路径的预期收益。 由此可知，计算期权价格即就是计算一个期望值，蒙特卡洛方法便是用于估计期望值，因此可以得到期