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

1、期权定价中的蒙特卡洛模拟方法期权作为最基础的金融衍生产品之一，为其定价一直是金融工程的重要研究领域，主要使用的定价方法有偏微分方程法、鞅方法和数值方法。而数值方法又包括了二叉树方法、有限差分法和蒙特卡洛模拟方法。蒙特卡洛方法的理论基础是概率论与数理统计，其实质是通过模拟标的资产价格路径预测期权的平均回报并得到期权价格估计值。蒙特卡洛方法的最大优势是误差收敛率不依赖于问题的维数，从而非常适宜为高维期权定价。1.预备知识两个重要的定理：柯尔莫哥洛夫(Kolmogorov)强大数定律和莱维一林德贝格(Levy-Lindeberg)中心极限定理。大数定律是概率论中用以说明大量随机现象平均结果稳定性的一

2、系列极限定律。在蒙特卡洛方法中用到的是随机变量序列同分布的Kolmogorov强大数定律：设g,g，为独立同分布的随机变量序列，若12Eg=卩8,k=1,2,则有P(liml兰g二卩)=1knsnkk=1显然，若g,g，,g是由同一总体中得到的抽样，那么由12n此大数定律可知样本均值1g当n很大时以概率1收敛于nkk=1总体均值卩。中心极限定理是研究随机变量之和的极限分布在何种情形下是正态的，并由此应用正态分布的良好性质解决实际问题。设g,g，为独立同分布的随机变量序列，若g-npEgk=，冰k=25k=心则有七-N(0,1)1ngpgp其价形式为limP(“k=ix)=1jexp(-)dt，

3、-gxmax0,S-K,美式看跌期权的终值条件为V(S,t)max0,K-S。然而，美式期权的价值没有解析解，我们一般可通过数值方法(蒙特卡洛模拟、有限差分法等)求得其近似解。风险中性期权定价模型如果期权的标的资产价格服从几何布朗运动=rdt+bdWS即标的资产的瞬时期望收益率卩取为无风险利率r。同理，根据伊藤公式可以得到Q2dInS=(r-)dt+cdWc2Q2InSInS=(r)(T-1)+c(WW)N(r)(T-1),c2(T-1)Tt2Tt2c2S=Sexp(r)(Tt)+c(WW)Tt2Tt对数正态分布的概率密度函数：设gN(卩Q2)，n=eg，则耳的密度函数为1(Inx卩)2、P(

4、x)=x2c2耳0根据上述公式，得到标的资产ST的密度函数如下xry2(InS-(r-)(T-1)2exp(-1一)2兀T-tox2y2(T-1)在风险中性概率测度下，欧式看涨期权定价为：V(S,t)二exp(-r(T-1)Eqmax0,S-KTEQmax0,Sxy2i(Ln-(r-2)(t-t)2-K=J+Wexp(-S一)dxko厉尸t2o2(T-1)丿Kxo2兀、Ttxo2(In-(r-)(T-1)2exp(-S一)dx2o2(T-1)接下来，求解以上风险中性期望。首先，对上式的右边第一个广义积分分别作变量替换In+-(r-罟)(T-1)y=T-7和u二y-yVT-1，可以得到j+8Ke

5、xp(-2o2(T-1)dx=Ser(T-tJ+sJeIn:-(r+0T-t)2兀o；T-tul.2du=Ser(T-t)SIn+(r+)(T-t)J-Ko-T-t-s1-ul.e-2du=Ser(T-t)N(d)v2k1再对等式的右边的第二个无穷积分，令o2，可求得oInx-InS-(r-)(T-1)u=2xo2K5-(r-)(T-t)2+s=exp(-S一)dxkxov2kT-12o2(T-1)KJ+s2du=o2InK-InS-(r-一)(T-t)2oT-tUnSInK+(r-)(T-1)ot-t-se-：du=KN(d)v2k2将以上的计算结果代入期望等式中，得到欧式看涨期权的价格公式

6、为：V(S,t)二e-r(T-t)Eqmax0,SK二SN(d)KemN(d)T12In+(r+兰)(T-1)d=K2其中，1一t,d2J-g、：T-t。可以看出，对于欧式看涨期权的风险中性定价方法的结果与基于资产复制的偏微分方程定价方法的结果是一致的。基于风险中性的期权定价原理在于：任何资产在风险中性概率测度下，对于持有者来说都是风险偏好中性的，便可用风险中性概率求取期权的期望回报再将其进行无风险折现便是初始时刻的期权价值。蒙特卡洛模拟方法就是一种基于风险中性原理的期权数值定价方法。2.蒙特卡洛模拟方法及其效率假设所求量e是随机变量n的数学期望e，那么近似确定e的蒙特卡洛方法是对n进行n次重

7、复抽样，产生独立同分布的随机变量序列q,n2,n，并计算样本均值12n_1ynnnnnk。那么根据Kolmogorov强大数定律有k1p(limne)】。因此，当n充分大时，可用厂作为所求TXnn量e的估计值。由中心极限定理可得到估计的误差。设随机变量n的方差Dn二Q28，对于标准正态分布的上8/2分位数Z&2,112P(代-91暮%而)uJZ82exp（-一）dt=1-8-Z822这表明，置信水平1-8对应的渐近置信区间是实际上，由此可确定蒙特卡洛方法的概率化误差边界，其误差为Z82匚,误差收敛速度是O(n-i2)。不难看出，蒙特卡洛方法的误差是由q和决定的。在对同一个“进行抽样的前提下,若

8、想将精度提高一位数字,要么固定Q，将n增大100倍；要么固定n将Q减小10倍。若两个随机变量3“2的数学期望E叫=En2=，Q12，那么无论从n或H中抽样均可得到。的蒙特卡洛估1212计值。比较其误差，设获得n的一个抽样所需的机时为t，那么在时间T内生成的抽样数niii2，则需使QtQt。因而，若要提高蒙特卡罗方法的效率，不能1122单纯考虑增加模拟的次数n或是减小方差Q2，应当在减小方差的同时兼顾抽取一个样本所耗费的机时，使方差Q2与机时t的乘积尽量的小。3.蒙特卡洛模拟方法为期权定价的实现步骤期权定价的蒙特卡洛方法的理论依据是风险中性定价原理：在风险中性测度下，期权价格能够表示为其到期回报

9、的贴现的期望值，即P=Eqexp(-rT)f(S,S，,S)，其中12T的EQ表示风险中性期望，r为无风险利率，T为期权的到期执行时刻，f(S,S，,S)是关于标的资产价格路径的预期2T收益。由此可知，计算期权价格即就是计算一个期望值，蒙特卡洛方法便是用于估计期望值，因此可以得到期权定价的蒙特卡洛方法。一般地，期权定价的蒙特卡洛模拟方法包含以下几步(以欧式看涨期权为例):在风险中性测度下模拟标的资产的价格路径将时间区间0，T分成n个子区间0=tttt=T，标012n的资产价格过程的离散形式是Sj(+)=Sj(t”)e(r-22)(-0)+。叽-佔，zN(o,i)计算在这条路径下期权的到期回报，

10、并根据无风险利率求得回报的贴现Cj=exp(-rT)maxH0,Sj-K重复前两步，得到大量期权回报贴现值的抽样样本m求样本均值，得到期权价格的蒙特卡洛模拟值C=exp(-rCj=MCmj=1另外，我们还可以得到蒙特卡洛模拟值与真值的概率化误差边界，这也是蒙特卡洛方法为期权定价的优势之一。由于Cj=exp(-rT)maxt,Sj-K，m条路径的收益均值为Tc=Cj，m条路径的方差为c=1(Cjc)，则可i=1i=1meanmvarm-1mean得95%的置信区间为Cmean-1.96mean+196va。m例1：假设无红利的股票A,初始价格为6,价格过程服从几何布朗运动，年预期收益率为10%，

11、收益率的波动率为每年25%，时间步长为0.01年（1年为100时间步），给定数据，S=6,卩=0.1,c=0.25,以及d=100,用蒙特卡洛方0法模拟资产的价格路径如下：5.21JrS(t+5t)=S(t)e(0-1-20-252)0-01+0-2501A01）rcePi6.66.46.265.85.65.4102030405060Period708090100图（1）蒙特卡洛方法模拟股票A价格路径，图（2）蒙特卡洛方法模拟股票B价格路径。若无红利的股票B、C、D,其价格均为6,股票B的期望收益率为0.1,波动率为0.6；股票C的期望收益率为0.5,波动率为0.25；股票D的期望收益率为05

12、,波动率为0.6,分别用蒙特卡洛方法模拟该三种股票在一年内的价格路径如下：S(t)e(0.1-丁0.62)0.01+0.6勺0.01.S(t)e(0.5-0.252)0.01+0.25、0.01i1i8S(t+8t)=S(t)e(0-5-20-62)0-01+.601iDirceP6.16.0565.955.9MonteDarioPricePathSimulationA1020304050Period3）607080901005.850（4）图（3）蒙特卡洛方法模拟股票C价格路径，图（4）蒙特卡洛方法模拟股票D价格路径。从图中可以看出，股票C和股票D的价格上升速度较快，而股票B和股票D的价格波

13、动比较大。这是与股票C和股票D价格的期望收益率较高，股票B和股票D价格的波动率较高相对应的。欧式看涨期权S=6,K=2,r=0.1,b=0.25,T=1，通过0Black-Scholes公式计算得的精确值为c=4.1903，蒙特卡洛模拟的价格为C=4.1787，其蒙特卡洛模拟图如下：5.45.254.8iittsanoEm4.6E4.44.243.800.511.522.533.54log(N)(5)上述同样的条件，路径由100逐渐增加到1000000条，对应地分别得到的期权价值的模拟值和置信区间，结果如下表所示：各种路径下蒙特卡洛方法模拟的95%置信区间N模拟值置信区间1004.31464.

14、0112,4.61805004.22624.0962,4.356310004.22134.1287,4.313920004.16334.0984,4.228150004.16954.1280,4.2111100004.17874.1490,4.2083500004.19604.1826,4.20941000004.18864.1791,4.198010000004.19144.1884,4.19444.蒙特卡洛模拟方法为我国权证定价权证是一种合同，权证投资者在约定时间内有权按约定价格向发行人购入或者出售合同规定的标的证券。权证发行人可以是标的证券的发行人或其之外的第三方。权证主要具有价格发现和

15、风险管理的功能，它是一种有效的风险管理和资源配置工具。现选取我国认股权证中的五粮YGC1、马钢CWB1、伊利CWB1为例，以2006年的价格作为样本区间模拟认股权证的价值，并将这些权证的蒙特卡洛模拟价值和由wind数据库给出的理论值进行比较。本例采用一年期短期利率2.52%作为无风险利率，用这些权证的正股股票价格序列来计算波动率。现实中用等时间间隔观测股票价格序列S(i=0丄2,.n),i股票投资的连续复利收益率u=ln(S/S),(i=12n),则uiii-11的样本标准差b丄兰(U-U)2。如果用日数据计算波动率，n一1：i=1则年度波动率按下式计算：年度波动率=日波动率*(每年的交易日数

16、)12将时间区间取为2006年12月1日2006年12月29日，则由蒙特卡洛方法模拟的认股权证价格与Black-Scholes模型的精确值和市场价格比较的结果如下：蒙特卡洛方法对五粮YGC1认股权证的模拟(b=51.15%)日期实际值蒙特卡洛模拟值理论值日期实际值蒙特卡洛模拟值理论值12-110.16410.0669.82112-1812.10013.52413.35112-410.12010.35710.12112-1912.08013.57413.40112-59.88010.63010.40112-2012.21013.77113.60112-69.39510.38610.15112-2

17、111.90013.37613.20112-79.1479.9989.75112-2211.42012.68712.50112-89.0509.7859.53112-2512.03813.74213.57112-119.8509.2258.95112-2611.97813.40613.23112-129.82510.60010.37112-2713.00114.36414.20112-139.76610.26010.02112-2813.05014.61214.45112-1410.58911.33211.12112-2914.50016.19816.05112-1510.84912.0281

18、1.831一一一一蒙特卡洛方法对马钢CWB1认股权证的模拟（o二53.91%）日期实际值蒙特卡洛模拟值理论值日期实际值蒙特卡洛模拟值理论值12-11.1431.2440.56912-181.7751.7091.05212-41.2091.1880.51712-191.8031.7091.05212-51.2411.2230.54912-201.7301.7561.10312-61.3491.2230.54912-211.6411.7091.05212-71.6331.4160.74312-221.7001.5420.77812-81.7501.6180.95212-251.7071.4530.

19、84812-111.9191.4160.74312-261.8351.5201.05212-121.8741.6180.95212-271.7761.7091.05212-131.7941.7481.09412-281.6441.8111.16312-141.7941.6330.96912-291.7081.7481.09412-151.8301.6330.969蒙特卡洛方法对伊利CWB1认股权证的模拟（o二62.03%）日期实际值蒙特卡洛模拟值理论值日期实际值蒙特卡洛模拟值理论值12-113.32413.53312.62912-1814.76014.81813.98812-413.25013

20、.94713.06912-1915.47915.54114.74812-513.29613.95713.07912-2015.48716.63015.88812-612.91113.95713.07912-2115.59416.44915.69812-712.85313.28812.36912-2215.16816.57315.82812-812.73412.76311.80912-2516.61615.81715.03812-1112.92012.57611.60912-2616.61917.75417.05812-1214.05912.94111.99912-2717.67317.8791

21、7.18812-1313.52814.10813.23912-2817.67319.72619.09812-1414.28113.81512.92912-2917.67319.72619.09812-1514.34914.61913.778从表可看出，由蒙特卡洛方法模拟的认购权证价格的模拟值比由Black-Scholes公式计算的理论值更接近实际值。为了更直观的比较，由蒙特卡洛方法模拟的认股权证价格与Black-Scholes模型的精确值和市场价格比较的结果如下图。其中SJ代表实际值，MC代表蒙特卡洛方法求得的模拟值，BS代表由Black-Scholes公式计算出的理论值。五粮YGC1价格模拟

22、比较图马钢CWB1价格模拟比较图2.51.50.5伊利CWB1价格模拟比较图2120191817161514131210从图中明显看出，五粮YGC1和伊利CWB1的模拟结果比较好，蒙特卡洛模拟值和Black-Scholes模型的理论值均与实际值吻合；而马钢CWB1的实证结果不理想，但是三种结果的走势图有共同的趋势。从比较分析中发现蒙特卡洛方法模拟的价格比Black-Scholes模型更接近实际价格。对于这些认股权证价格的模拟结果的好坏，受诸多因素影响，主要与选取的波动率和中国权证市场的发展特点有关等等。隐含波动率及其数值计算方法隐含波动率是一个在市场上无法观察到的波动率，是通过Black-Sc

23、holes期权定价公式计算出来的波动率。由于我们无法给出它的解析解，因此，只能借助于数值计算给出近似解。下面介绍牛顿迭代法计算隐含波动率。牛顿迭代法是牛顿在17世纪提出的一种在实数域上近似求解方程根的方法。步骤1.将函数f(x)在点x附近展开成泰勒级数f(x)=f(xo)+八x0)(x-x0)+%(x-x0)2+步骤2.取泰勒级数的前两项作为f(x)uf(x)+f(x)(x-x)000假设ff(x)主0，求解方程f(x)=f(x)+ff(x)(x-x)=0，并令其解为0000 x，得x=x-航，这样得到迭代公式x=x-伙，经0f(x)n+1nf(x)0n过n次迭代后，可以求出f(x)=o的近似

24、解。根据牛顿迭代法，隐含波动率的计算步骤如下：假设其他变量保持不变，认为函数f(q)=SN(d)-Ke-r(t-t)N(d)-C是隐含波动率q的一元函12Mar数，其中的c是市场上观察到的期权价格。Mar求函数f(q)的导数广Q)=S、TXe-d22BqV2F由迭代公式q=Q-Qi计算波动率，直至f(Q)8i+1if(Q)八/i(8是期望达到的精度)。此外，为了计算隐含波动率，经济学家和理财专家曾做过种种努力试图寻找一个计算波动率的公式。如Brenner和Subrahmanyam于1988年，Chance于1993年分别提出计算隐含波动率的公式，虽然这些公式对于持有平价期权的波动率的计算还算准

25、确，但是基础资产的价格一旦偏离期权的执行价格的现值，其准确性就会丧失。1996年，Corrado和Miller在前人研究的基础上建立了如下公式，大大提高了隐含波动率的计算的准确性：丄2x(C-頁S+Ke-rTSKe-rTSKe-rT(SKe-rT)2珂(C-)25.服从跳扩散过程的无形资产期权定价问题及其蒙特卡洛模拟分析服从跳扩散过程的期权定价方法正常的波动用几何布朗运动(Brown)来描述一由供需不平衡、利率变动或整个经济的波动等因素引起的。不正常的波动用泊松过程(Poisson)来描述一由未预料到的重要信息的出现引起的。这些信息在不连续的时间点出现，而且出现的时间点不确定，是否会出现也不确

26、定。带跳跃项的伊藤Ito公式：设v=V(S,t)，V是二元可微函数，若随机过程St服从随机微分方程=adt+bdW+dJ其中，dW是标准维纳过程，dJ表示不可预测的跳跃，且EAJ二0。则带跳跃项的伊藤Ito公式为dV.V1d2VdVdV(S,t)=+九乙(V(S+a,t)-V(S,t)p+edt+dS+dJdttii2dS2dSVi=1dJ=V(S,t)-V(S-,t)-九工(V(S+a,t)-V(S,t)pdtVtiii=1其中，S-=limS上式是对跳跃项作如下假定得出的：1、在两个跳跃之间Jt保持不变,而在跳跃时间j（j=1,2,）2、有k种跳跃类型，跳跃尺度为耳i=12,k，跳跃尺度为

27、a的概率为匕，跳跃的发生强度入t依赖于s的最终观测值，跳跃类型和尺度都是独立随机的。则在时间区间（t,t+At内，增量AJ为AJ=AN-XAt正ap）tttii这里N表示的是至时间t发生的跳跃大小的总和，XAt表示跳tt跃发生的概率，唧为跳跃的期望值，则AJ是不可预测的。i=1漂移参数a可看作两个漂移的和a=+X（工ap）tttii这里M”表示S冲连续运动的维纳过程部分，第二项为纯跳跃部分。将Poisson过程引入到期权定价模型中，得到标的资产价格价格的跳扩散方程如下=（卩-Xv）dt+QdW+（Y-1）dNJ1九dt其中，dN1-Xdt，标的资产价格的变化比率为Y,v=E(Y-1)，且y与d

28、N相互独立。令F二InS，根据带跳跃项的伊藤公式可得其微分形式为力F力F1力2FdF=QSdW)2+(F(YS,t)-F(S,t)dN瓦+_(v)SdtSdw)+2茹=(卩-22-Xv)dt+cdW+lnYdN整理上式，得到标的资产价格公式为S=Sexpt0-Xv)t+QW(t)ff()Yii=1在标的资产价格遵循跳扩散过程的假设下，根据上述带跳伊藤公式可得期权价值V=V(S,t)的微分形式如下TOC o 1-5 h zavcqv1a2VavdV=(+(p-Xv)S+Q2S2)dt+QSdW+(V(YS,t)-V(S,t)dNataS2aS2aS构造期权与标的资产的无套利资产组合口，=v-QV

29、s，其微分qs形式为dnr=dv-avds=(+1q2S2asat2竺)dt+(V(YS,t)-V(S,t)dN-(Y-1)S巴dNaS2aS则该无套利资产组合微分形式的期望如下式E(dnr)=(空+1q2S2at2=+Q2S2at2a2VaV)dt+E(V(YS,t)-V(S,t)dN-E(Y-1)SdNaS2aSa2VaV+XEV(YS,t)-V(S,t)-XvSdtaS2aS由于资产组合n=v-空s为无风险组合，因此有如下等式成as立E(d口)=rnrdt=r(V-S)dtaS两式联立并化简得到标的资产价格遵从跳扩散过程的定价公式如下：空+(r-Xv)S空+1Q2S2型-rV+九EV(Y

30、S,t)-V(S,t)=0dtdS2dS2若没有发生跳跃事件，则X=0，将其代入上式所得结果与Black-Scholes微分方程完全一致。当期权分别为欧式看涨、欧式看跌、美式看涨和美式看跌期权时，其边界条件和终值条件与本章第一节的终边值条件相同。TOC o 1-5 h zJli=0Merton假设标的资产价格跳跃高度服从科oi工0，从而推导出欧式看涨期权的定价公式为：V(S,T)=艺e-Xt(XtM)E(W(SHYe-xvtt,K,Q2,r)N(t)!N(t)iN(t)=0i=1其中，W(S,T,K,Q2,r)=S(d)-Ke-rT(,另外，Harworth假设跳跃高度Y服从对数正态分布InY

31、N(匕,Q计，贝0欧式看涨期权的解析解为N(t)!V(S,t)二芸eXt(Xt)w(S,t,K,Q2,r)N(t)!N(t)=0其中，X=Xexp(口+-Q2)r=rX(exp(口+-Q2)一1)+1(口+-q2)Y2Y,Y2YTY2Y,例2.标的资产价格遵从跳扩散过程如下7Q=(y-Xv)dt+QdW+(Y-1)dNSS(0)=20,y=2.5%,Q=20%,X=0.5t1.5v=Y-1,n=500,At=0.004,Y=0.8用蒙特卡洛模拟的资产价格路径如下图所示：无形资产专利池的期权定价模问题专利池的市场价值V依赖于企业使用专利池技术前后生产产品所获得的收益S和成本C及时间t,这三个变量

32、均可用跳扩散模型：dx=(卩一九v)dt+QdW+(Y-1)dNI=J通过构造由V和它所依赖的两个变量S、C组成的资产组合，利用带跳的伊藤引理获得V与S、C所遵循的带跳的随机微分方程，并根据实际情况在一些假设条件下给出该方程的终边值条件，最终获得V的求解公式。构造无风险资产组合n=V-vs-vcSSC一方面n的微分的期望为：E(dn)二r(V-VS-VC)dtVvsc另一方面，E(d口)=(V+1Q2S2V+1Q2C2V+QQSCV)dtst2sss2cccscsc+九E(V(YS,C,t)-V(S,C,t)dt九vVSdtsssssa新产品发明专利池的市场价值V所遵循的方程为V+(r九v)V

33、S+rVC+1a2S2V+1a2C2VtSSSC2SSS2CCC+aaSCV+九E(V(YS,C,t)V(S,C,t)rV=0SCSCSS(S,C,T)=max(S(T)C(T),0)=(S(T)C(T)+(S,C,t)T0asST0(S,C,t)T0asCTg(S,C,t)TSasSTgInYN(p,a2),v二E(Y1)二e%+2aSy)1SSYSYSS期权的价格公式：2SYe-(兀T)Ns(t)N(t)!S-V(S,T,C,a2,r)BSNs(t)=02,a=a2+a2+aaSCS+1a2N(t)2SYSc1r=r尢(epsy+2aS2SY1)+竹(p+1a2)TSY2SY20世纪90年

34、代初,由高分子工程材料的某高校、研究所、设计院和高新技术企业等经过两年的开发研究，研制出新型建材铝塑复合管全套生产工艺，该技术已获多项国家发明专利，且己具备成套设备生产供应能力。当时，该技术在国内只此一项，属新产品发明专利池技术。且其专利技术使用寿命长达50年以上，受专利保护20年。但该技术在国外存在多家供方，不同供方在核心技术内容、原理、流程上基本一致，同时也不排除在一段时间后出现其他更好技术的可能性，一方面时间越长，这种可能性越大。另一方面该技术使用寿命越长，这种可能性越小(l=l(t)。并且,其他同类技术的出现使该专利池技术的收益下降，下降幅度为LnYo因为设备的经济使用寿命是20年，根

35、据市场需求,计划建成一条年生产100吨的生产线，其20年的成本,包括设备的直接制造成本和运营期间的管理费、工资等。若在期初计划投资1000万，以后20年每年的生产量不变，生产成本按每年的通货胀率10%递增。假设在初期预计该项技术20年总收益为4000万，其收益率为25%,方差为20%。九(t)=0.02t1.3,卩=25%,r=卩=10%,Y=0.6SSCSS(0)=4000,C(0)=1000,n=4000,At=0.005新产品发明专利池的市场价值V=8050在一次付清许可费用情况下的价格模型：新产品发明专利池的价格P所遵循的方程为:P+(r一九v)PS+rPC+丄G2S2P+丄G2C2P

36、tsssc2sss2ccc+gGSCP+九E(P(YS,C,t)-P(S,C,t)-rP=0SCSCSSP(S,C,T)=max(a(S(T)-C(T),0)P(S,C,t)T0asST0P(S,C,t)T0asCTgP(S,C,t)taSasStg在一次付清许可费用情况下的新产品发明专利池的价格为:九(t)二0.02ti.3,卩二25%,r二卩二10%,a二0.5,Y二0.6SSCSS(0)二4000,C(0)二1000,n二4000,At二0.005在一次付清许可费用情况下新产品发明专利池的价格P=5450。在首付加每期按收益固定比率支付许可费用情况下的价格模型新产品发明专利池技术产生的收

37、益S遵循模型7Q一=(卩-q九v)dt+cdW+(Y-1)dNSSSSSSSS引进新产品发明专利池技术后的成本C遵循模型dC=ydt+CdWCcc构造无风险资产组合口=P-PSPCPSC一方面口的微分的期望为E(d口)二r(P-PS-PC)dtPPSC新产品发明专利池的价格P所遵循的方程为：另一方面，口的微分及其期望为：PE(d口)=(P-qPS+1Q2s2P+-c2c2P+QQSCP)dtptss2sss2cccscsc+九E(P(Ys,c,t)-P(s,c,t)dt-九vPsdtsssss新产品发明专利池的价格P所遵循的方程为:P+(rq九v)Ps+rPc+-Q2s2P+-c2c2Ptss

38、ssc2sss2ccc+QQscP+九E(P(Ys,c,t)P(s,c,t)rP=0scscssP(s,c,T)=max(a(S(T)c(T),0)P(s,c,t)T0assT0P(s,c,t)T0asCTgP(s,c,t)asasSgInYN(卩,c2),v=E(Y1)=e(b+2q；y)1TOC o 1-5 h zssysyss期权的价格公式：P(s,c,t)“三e-H)：s(t)V(se,T,c,Q2,r)N(t)!bsNs(t)=0sV(se-qsT,T,c,Q2,r)=s(t)e-qj(d)C(t)e-孔(d),T=Ttbs12兀=九esy+2Qsyqs=Q2+Q2+QQscsc+1

39、Q2N(t)2sys.zu+1Q2八N(t)/1、r=rq九(e怙+2Qsy1)+s(|LX+Q2)sstsy2sy九(t)=0.02t1.3,y=25%,r=10%,q=10%,Y=0.6sscsss(0)=4000,c(0)=1000n=4000,At=0.005在首付加每期按收益固定比率支付许可费用情1=1况下新产品发明专利池的价格P=855。6.最小二乘蒙特卡洛模拟与美式期权定价运用最小二乘蒙特卡洛模拟方法为美式期权定价的基本原理与蒙特卡洛模拟方法基本相同，并且用最小二乘回归同时还可解决各样本时点上继续持有期权价值的确定和各样本路径的最优停时的确定。其基本思路是：在期权的有效期内，将其

40、标的资产价格过程离散化，随机模拟出标的资产价格的多条样本路径，从而得到每个时刻资产价格的截面数据。选取以某时刻资产价格为变量的一组基函数作为解释变量，下一时刻期权价值的贴现值作为被解释变量，进行最小二乘法回归求得该时刻期权的持有价值，并与该时刻期权的内在价值作比较，若后者较大，则应该立即执行期权，否则，就应继续持有期权。最小二乘蒙特卡洛模拟方法定价的基本实现步骤：首先，随机生成标的资产价格的多条样本路径；然后，从到期时刻逆向求解，比较期权的内在价值与持有价值，确定出各时刻期权价值和每条样本路径的最优停时；最后，将所有样本的的期权价值求取按无风险利率贴现的算数平均值便是模拟的期权价值。面，我们运

41、用最小二乘蒙特卡洛模拟方法对单个标的资产的美式看跌期权进行定价，其算法实现步骤如下：第一步：随机生成标的资产价格过程的多条样本路径现设一单个标的资产美式看跌期权的持有到期日为T，期权的执行时刻为T*，t*e0,T，标的资产价格为S，期权的执行价格为X。在风险中性条件下，该期权的初始时刻价值为：P=Eqexp(rt*)f(S,S，,S，,S)01t*T其中，S0,S1,S*，,ST为标的资产价格的路径，f(S0,S1,S*，,ST)是在最优执行时刻t*的期权价值。上式定义的P便是将要运用最小二乘蒙特卡洛方法进行模拟的期权价值。将期权的存续区间0,T均分为N个子区间，则每个子区间的长度为At=N，

42、标的资产价格过程的离散形式：S=Sexp(r)At+gAt)ii12i1其中，ie0丄,N，随机变量s服从标准正态分布。因此，利用生成随机数模拟得到标的资产价格S的一条样本路径S0,气，,SN，重复执行M次模拟，我们可得到资产价格的总样本SMN+1)。第二步：计算各个样本的最优停时及各时刻的期权价值对于美式看跌期权，在期权的有效时刻i，样本路径j上的内在价值为I(Sij)=maxXSij,0，持有价值为H(S/)Eexp(rAt)C(Si十)S/。由于美式期权在有效期的任何时候都可行权，所以必须比较该时刻期权的内在价值1(Sj)与持有价值H(叨的大小，以确定该时刻的期权价值以及是否执行期权，即

43、C(Sj)=maxI(Sj),H(Sj)iii由期权的持有价值表达式可知它依赖于下一步期权决策的价值，需通过逆向求解这个期望价值，这正是普通的蒙特卡洛模拟法为美式期权定价的难点所在。小二乘蒙特卡洛模拟方法通过建立一个当前时刻标的资产价格与下一时刻期权价值贴现值的线性回归计量模型：Y=exp(-rAt)C(Sj)=a+aSj+aSj2+卩ji+l12i3ij上述模型以所有样本路径在时刻i的价格S和S2作为解释变量，对应的下一时刻期权价值的现值作为被解释变量。采用普通最小二乘法进行回归，求得回归系数的估计值a2,a3和样本回归方程Y=ai+aS+吧2；再将各个资产价格样本代入到回归方程分别可以得到

44、其期权的持有价值估计值，Y=EYSj=Eexp(-rAt)C(S)|Sj根据计量经济学的理论，这个估计值就是在标的资产价格S/下的期权持有价值的无偏估计值。另外，本例中选取基函数L0(S)=1；L1(S)=S；L2(s)=S2作为解释变量，根据实际情况中也可以选取其他形式的基函数：L(S)=exp(-);七(S)=exp(-)(1-S);S)=exp(-)(1-2S+可)0212222作为解释变量。执行看现在，我们从到期日开始倒推计算求解每条样本路径上的最优停时和每个样本点的期权价值。在到期日T，跌期权的价值为maxX-Sj,0。接着，判断在N-1时刻是否行权。若期权处于实值状态，即SN-1X

45、，则与继续持有期权的价值相比较，若内在价值大于持有价值，则应立即执行期权；否则，继续持有期权。考虑在该时刻期权处于实值的样本子集，近似期权持有价值的回归方程为：yj=exp(rAT)maxXSj,0=a+aSj+aSj2+卩N1N12N13N1j其中，jeJN1，JN1是N-1时刻所有期权处于实值状态的标的资产价格样本集。在时刻N-1的资产价格si,信息下，比较内在价值maxXSj,0与继续持有期权的价值yj1就可做出是否执行期权的决策。同理，我们可倒推继续求得时刻N2,N3,0介值。对于每条样本路径j，期权或是在最优停时t严0丄,N执行，或是永不执行。具体设计程序时，令初值丁N，在时刻j1，

46、如果继续持有期权，则号不变；如果执行期权，则jN-1，依此类推。每个样本上就只有一个最优停时，每次更新：，最后便求得每条样本路径上的最优停时。第三步：对各条样本路径上的期权价值按无风险利率贴现并求其均值经过M次模拟后，得到M条标的资产价格的样本路径，以及每条样本路径上的最优停时；和在该时刻的期权价值：I(Sj)=maxXSj,0j=1,2,Mt*t*由于每条样本路径上的最优执行时间不同，期权价值的贴现因子e-rt；也不同，所以应分别进行贴现求均值，最终得到初始时刻期权价值的最小二乘蒙特卡洛模拟值：弋exp(rt*)I(Sj)TOC o 1-5 h zjt*P二Eqexp(rt*)f(S,S，S

47、，S)二4-01t；TM例3：已知股票价格为50，美式看跌期权执行价为50到期日为5个月，股票年收益率的标准差为0.4，无风险利率为10%，用最小二乘蒙特卡洛模拟其价格。编制最小二乘蒙特卡洛模拟的MATLAB程序如下：functionprice=AmericanOptLSM(S0,K,r,T,sigma,N,M)dt=T/N;R=exp(r-sigmaA2/2)*dt+sigma*sqrt(dt)*randn(N,M);S=cumprod(S0*ones(1,M);R);ExTime=N*ones(M,1);CF=zeros(size(S);CF(end,:)=max(K-S(end,:),0

48、);forii=N:-1:2Idx=find(S(ii,:)C;nIdx=setdiff(1:M),Idx(Jdx);CF(ii,Idx(Jdx)=max(K-X(Jdx),0);ExTime(Idx(Jdx)=ii;CF(ii,nIdx)=exp(-r*dt)*CF(ii+1,nIdx);endPrice=mean(CF(2,:)*exp(-r*dt)%绘制标的股票价格模拟图%x1=0:N;y1=S;y2=mean(S);subplot(2,1,1)plot(x1,y1)subplot(2,1,2)plot(x1,y2)xlabel(期权存续期间)ylabel（股价的模拟路径）%绘制期权价值

49、模拟图%figure;x2=1：N;y3=CF(2:end,:);fori=1:My4(i)=y3(i,ExTime(i);endplot(x2,y3,ExTime,y4,*)xlabel（期权的最优停止时间）ylabel（期权价值的模拟路径）模拟的美式看跌期权的价格路径如下图所示:模拟的期权价值路径及其最优停时如下图:ia1225新巧罚本例中的美式看跌期权价格为：price=AmericanOptLSM(50,50,O1,5/12,O4,50,1OOOOO)Price=426547.改进蒙特卡洛方法计算效率的常用几种方差减少技术方差减少技术的共性是利用模型特点，调整或修正模拟的输出变量，从而

50、降低估计值的方差。在采用方差减少技术时，要具体问题具体分析，针对不同期权类型的特点应用相关的方差减少技术，从而取得效率的最大改进。.对偶变量(Antitheticvariates)技术对偶变量技术是最简单和最常用的方差减少技术。以标准欧式看涨期权为例，其标准蒙特卡洛估计值为exp(-rT)区max0,Sj-KC二exp(-rT)XCj=MCmm标的股票的股价终值抽样为Sj=Se(r-2a2)T+aTzj,j=1,2,mT0由概率论的知识可知-Z也是标准正态分布中相互独立的抽样值,那么用_Z代替Z得到的勺也是股票价格终值的抽样，从而由C=exp(-rT)max也能jT得到期权价格的无偏估计量。因此，由对偶变量技术得到的期权价格蒙特卡洛估计值为C=1Cj+C。AVm2j=1对偶变量技术的有效性：由于VarC=VarC，所以jjVarCj+Cj=-(VarC+CovC,C)；并且，令C=(Z)，对于标22jjjjj准欧式看涨期权，(z)是单调递增函数。由不等式jE(Z)(-Z)E(Z)E(-Z)，可知CovC,C0，从而jjjjjjVarCj+Cj0ng(x)0,VxGRd，则有e=Jh(x)f(x)g(x)dxg(x)将上述积分写成关于密度g的期望形式，可以得到e二Eh(X)二Eg