期权定价中的蒙特卡洛模拟方法期权作为最基础的金融衍生产品之一，为其定价一直是金融工程的重要研究领域，主要使用的定价方法有偏微分方程法、鞅方法和数值方法。
而数值方法又包括了二叉树方法、有限差分法和蒙特卡洛模拟方法。
蒙特卡洛方法的理论基础是概率论与数理统计，其实质是通过模拟标的资产价格路径预测期权的平均回报并得到期权价格估计值。
蒙特卡洛方法的最大优势是误差收敛率不依赖于问题的维数，从而非常适宜为高维期权定价。
§1. 预备知识◆两个重要的定理：柯尔莫哥洛夫(Kolmogorov)强大数定律和莱维一林德贝格(Levy-Lindeberg)中心极限定理。
大数定律是概率论中用以说明大量随机现象平均结果稳定性的一系列极限定律。
在蒙特卡洛方法中用到的是随机变量序列同分布的Kolmogorov 强大数定律：设12,,ξξ为独立同分布的随机变量序列，若[],1,2,k E k ξμ=<∞=则有11(lim )1nk n k p n ξμ→∞===∑ 显然，若12,,,n ξξξ是由同一总体中得到的抽样，那么由此大数定律可知样本均值11nk k n ξ=∑当n 很大时以概率1收敛于总体均值μ。
中心极限定理是研究随机变量之和的极限分布在何种情形下是正态的，并由此应用正态分布的良好性质解决实际问题。
设12,,ξξ为独立同分布的随机变量序列，若2[],[],1,2,k k E D k ξμξσ=<∞=<∞=(0,1)nkd n N ξμ-−−→∑其等价形式为211lim()exp(),2nxk k n t n P x dt x ξμσ=→∞-∞-≤=--∞<<∞∑⎰。
◆Black-Scholes 期权定价模型 模型的假设条件：1、标的证券的价格遵循几何布朗运动dSdt dW S μσ=+其中，标的资产的价格S 是时间t 的函数，μ为标的资产的瞬时期望收益率，σ为标的资产的波动率，dW 是维纳过程。
2、证券允许卖空、证券交易连续和证券高度可分。
3、不考虑交易费用或税收等交易成本。
4、在衍生证券的存续期内不支付红利。
5、市场上不存在无风险的套利机会。
6、无风险利率r 为一个固定的常数。
下面，通过构造标的资产与期权的资产组合并根据无套利定价原理建立期权定价模型。
首先，为了得到期权的微分形式，先介绍随机微积分中的最重要的伊藤公式。
伊藤Ito 公式：设(,)V V S t =，V 是二元可微函数，若随机过程S 满足如下的随机微分方程(,)(,)dSS t dt S t dW Sμσ=+则有22221((,)(,))(,)2V V V V dV S t S S t S dt S t S dW t S S Sμσσ∂∂∂∂=+++∂∂∂∂根据伊藤公式，当标的资产的运动规律服从假设条件中的几何布朗运动时，期权的价值(,)V V S t =的微分形式为22221()2V V V V dV S S dt S dW t S S Sσμσ∂∂∂∂=+++∂∂∂∂现在构造无风险资产组合V V S S∂∏=-∂，即有d r dt ∏=∏，经整理后得到2222102V V VS rS rV t S S σ∂∂∂++-=∂∂∂这个表达式就是表示期权价格变化的Black-Scholes 偏微分方程。
它同时适合欧式看涨期权、欧式看跌期权、美式看涨期权和美式看跌期权，只是它们的终值条件和边界条件不同，其价值也不相同。
欧式看涨期权的终边值条件分别为{}(,)max 0,T V S T S K =-，00(,)S V S T S S →⎧=⎨→∞⎩ 通过求解带有终边值条件的偏微分方程，得出欧式看涨期权的的解析解：()12(,)()()r T t V S t SN d Ke N d --=-其中，22()x dN d edx--∞=，21d =21d d =-T 为期权的执行日期，K 为期权的执行价格。
欧式看跌期权的终边值条件分别为{}(,)max 0,T V S T K S =-，(,)0K S V S T S →⎧=⎨→∞⎩ 此外，美式看涨期权的终值条件为(,)max{0,}V S t S K ≥-，美式看跌期权的终值条件为(,)max{0,}V S t K S ≥-。
然而，美式期权的价值没有解析解，我们一般可通过数值方法（蒙特卡洛模拟、有限差分法等）求得其近似解。
◆风险中性期权定价模型如果期权的标的资产价格服从几何布朗运动dSrdt dW S σ=+即标的资产的瞬时期望收益率μ取为无风险利率r 。
同理，根据伊藤公式可以得到2ln ()2d S r dt dWσσ=-+222ln ln ()()()~(()(),())22T t T t S S r T t W W N r T t T t σσσσ-=--+----2exp(()()())2T t T t S S r T t W W σσ=--+-对数正态分布的概率密度函数：设2~(,)N ξμσ，e ξη=，则η的密度函数为22(ln ))0()200x x P x x ημσ-- > = ≤ ⎩根据上述公式，得到标的资产T S 的密度函数如下222(ln()())2)0 ()2()00txr T tSxP xT txσσ⎧---⎪- >=-⎪ ≤⎩在风险中性概率测度下，欧式看涨期权定价为：(,)exp(())[max{0,}]QTV S t r T t E S K=---()()222222(ln()())2[max{0,}])2(ln()())2)2QT KKxr T tSE S K dxT txr T tS dxT tσσσσ+∞+∞----=--------⎰⎰接下来，求解以上风险中性期望。
首先，对上式的右边第一个广义积分分别作变量替换2ln()()xr T tyσ---=和u y=-，可以得到()2222222ln()()()()()221 ln()()(ln()())2)2() KSr T tu ur T t r T t r T tKr T txr T tS dxT tSe du Se du Se N dσσσσ+∞++-+∞------+------===⎰⎰再对等式的右边的第二个无穷积分，令2ln ln()()x S r T tuσ----=，可求得2222222ln ln()()2222 ln ln()(2(ln()())2)2()() KS K r T tu uK S r T txr T tS dxT tK du K du KN dσσσσ+∞-+--+∞-----------===⎰⎰将以上的计算结果代入期望等式中，得到欧式看涨期权的价格公式为：()()12(,)[max{0,}]()()r T t Q r T t T V S t e E S K SN d Ke N d ----=-=-其中，21ln ()()S r T t d σ++-=，21d d =-。
可以看出，对于欧式看涨期权的风险中性定价方法的结果与基于资产复制的偏微分方程定价方法的结果是一致的。
基于风险中性的期权定价原理在于：任何资产在风险中性概率测度下，对于持有者来说都是风险偏好中性的，便可用风险中性概率求取期权的期望回报再将其进行无风险折现便是初始时刻的期权价值。
蒙特卡洛模拟方法就是一种基于风险中性原理的期权数值定价方法。
§2. 蒙特卡洛模拟方法及其效率假设所求量θ是随机变量η的数学期望[]E η，那么近似确定θ的蒙特卡洛方法是对η进行n 次重复抽样，产生独立同分布的随机变量序列12,,,n ηηη，并计算样本均值11nn k k n ηη==∑。
那么根据Kolmogorov 强大数定律有(lim )1n n p ηθ→∞==。
因此，当n 充分大时，可用n η作为所求量θ的估计值。
由中心极限定理可得到估计的误差。
设随机变量η的方差2[]D ησ=<∞，对于标准正态分布的上2δ分位数Z δ，有22(exp()12Z n Z t p Z dt δδδηθδ--<≈-=-⎰这表明，置信水平1δ-对应的渐近置信区间是n Z δη±。
实际上，由此可确定蒙特卡洛方法的概率化误差边界，其误差为Z δ，误差收敛速度是12()O n -。
不难看出，蒙特卡洛方法的误差是由σ和决定的。
在对同一个η进行抽样的前提下，若想将精度提高一位数字，要么固定σ，将n 增大100倍；要么固定n 将σ减小10倍。
若两个随机变量12,ηη的数学期望12[][]E E ηηθ==，12σσ≠，那么无论从1η或2η中抽样均可得到θ的蒙特卡洛估计值。
比较其误差，设获得i η的一个抽样所需的机时为i t ，那么在时间T 内生成的抽样数i iT n t =，若使<，则需使1122t t σσ<。
因而，若要提高蒙特卡罗方法的效率，不能单纯考虑增加模拟的次数n 或是减小方差2σ，应当在减小方差的同时兼顾抽取一个样本所耗费的机时，使方差2σ与机时t 的乘积尽量的小。
§3. 蒙特卡洛模拟方法为期权定价的实现步骤期权定价的蒙特卡洛方法的理论依据是风险中性定价原理：在风险中性测度下，期权价格能够表示为其到期回报的贴现的期望值，即12[exp()(,,,)]QT P E rT f S S S =-，其中的QE 表示风险中性期望，r 为无风险利率，T 为期权的到期执行时刻，12(,,,)T f S S S 是关于标的资产价格路径的预期收益。
由此可知，计算期权价格即就是计算一个期望值，蒙特卡洛方法便是用于估计期望值，因此可以得到期权定价的蒙特卡洛方法。
一般地，期权定价的蒙特卡洛模拟方法包含以下几步（以欧式看涨期权为例）:(l)在风险中性测度下模拟标的资产的价格路径 将时间区间[0,]T 分成n 个子区间0120n t t t t T =<<<<=，标的资产价格过程的离散形式是211()()21()()i i ir t t jji i S t S t eσσ+--++=，~(0,1)i z N(2)计算在这条路径下期权的到期回报，并根据无风险利率求得回报的贴现{}exp()max 0,j j T C rT S K =--(3)重复前两步，得到大量期权回报贴现值的抽样样本 (4)求样本均值，得到期权价格的蒙特卡洛模拟值{}11exp()max 0,1exp()mj T mj j MCj rT S K C rT C m m==--=-=∑∑另外，我们还可以得到蒙特卡洛模拟值与真值的概率化误差边界，这也是蒙特卡洛方法为期权定价的优势之一。
由于{}exp()max 0,j j T C rT S K =--，m 条路径的收益均值为11m jmeani C C m ==∑，m 条路径的方差为2var11()1m j mean i C C C m ==--∑，则可得95%的置信区间为[mean mean C C -+。
