第8章--蒙特卡洛模拟金融衍生产品定价

第8章蒙特卡洛模拟金融衍生产品定价

本章介绍蒙特卡洛模拟期权定价的内容，要求读者掌握随机数生成方式，了解蒙特卡洛定价就是模拟风险中性测度下标的资产的运动过程，学会蒙特卡洛方法模拟欧式期权定价，掌握提高模拟精度的常用方法。

§8.1 随机模拟基本原理

1977年，菲力埔·伯耶勒(Phelim Boyle)提出了模拟方法求解金融资产定价问题，其想法是假设资产价格分布是随机波动，如果知道了这个波动过程，就可以通过随机模拟不同的路径，每做完一次模拟，就产生了一个最终资产价值，再进行若干次这样的过程，那么所得到的结果就是一个最终的资产价值分布，从这个分布中我们可以得到期望的资产价格。

8.1.1 随机数生成函数

1．均匀分布随机数生成函数

MATLAB中的unidrnd函数可以生成1到N的均匀分布随机数。

调用方式

R=unidrnd(N);

R=unidrnd(N,m);

R=unirnd(N,m,n);

其中，N所要生成的随机数个数，m确定输出随机矩阵R的行数,n确定输出随机矩阵R的列数

2．生成服从连续均匀分布的随机数

如果需要生成服从连续分布的随机数，则需调用unifrnd函数，其调用格式为调用方式1

R＝unifrnd(A，B)生成位于A、B之间的一个随机数。

调用方式2

R＝unifrnd(A,B,m)生成位于A、B之间的随机数。m=[m1,m2]表示行数列数。

调用方式3

R＝unifrnd(A,B,m,n),m,n分别表示行数、列数

unifrnd(1,2,[5,6]),unifrnd(1,2,5,6)

8．1．2 生成正态分布随机数

调用方式

R=normrnd(mu,sigma)

R=normrnd(mu,sigma,m)

R=normrnd(mu,sigma,m,n)

8．1．3 特定分布随机数发生器

MATLAB中有统一格式的随机数发生器，函数名称为random，可生成许多服从不同分布的随机数。

y=random('name',A1,A2,A3,m,n)

表8.1 生成特定分布的随机数函数参数表

beta分布：beta,二项分布：bino,卡方:chi2,指数分布：exp,F-分布：f, Gamma:gam Lognormal:logn, uniform:unif;Poisson:poiss,T:t;Normal->norm;

Noncentral F ->ncf, Noncentral->nct

8．1．4 蒙特卡洛模拟方差削减技术

蒙持卡洛模拟精度与模拟次数密切相关，模拟次数越高其精度越高，但是次数增加又会增加计算量。实践证明明减少模拟方差可以提高稳定性，减少模拟次数。有很多种方法可以减小方差，如对偶变量技术、控制变量技术、分层抽样、矩匹配、条件蒙特卡洛模拟等，但最简单并且应用最为广泛的是对偶变量技术与控制变量技术。

对偶变量技术就是先随机抽样得到一组数据，然后以此为基础构造出另一组对偶变量。下面以正态分布为例介绍对偶变量技术。首先从正态分布变量中随机抽取N 个样本值，分别为)N ,...,2,1i (Z i =，由此可以得到N 个模拟值

)N ,...,2,1i (C i =，那么衍生证券蒙特卡洛估计值为

∑=i

i C N 1C ˆ 以)N ,...,2,1i (Z i =为基础，

i Z ~是与i Z 相互对偶的随

机数，由正态分布的性质知，)N ,...,2,1i (Z Z ~i i =-=也是服从正态分布，由对偶随

机数生成的估计值为

∑=i i C ~N 1C ~ 对C ~和C

ˆ取平均得到新的估计 ∑⎪⎪⎭

⎫ ⎝⎛+=+=i i i 2C ˆC ~N 1)C ˆC ~(21C 如果随机抽样的样本)N ,...,3,2,1i (Z i =模拟得到的估计值比较小，那么与之对

偶的随机抽样样本)N ,...,2,1i (Z Z ~i i =-=得到的估计值可能会偏大，二者的平均值

就可能会接近真实值。如果0)C ~,C ˆcov(i i ≤，那么

)C ˆvar(21)C ~,C ˆcov(21)C ˆvar(212C ~C ˆvar i i i i i i ≤+=⎪⎪⎭

⎫ ⎝⎛+ 从上面的不等式可以看出，利用对偶技术可以增加估计稳定性，提高了估计精确度。

8.1.5 随机模拟控制变量技术

摔制变量技术就是将与所估计的未知变量密切相关的另一个已知量的真实值和估计值之间的差异作为控制量，以提高估计精度。在定价实践中，将这两种衍生证券用相同的随机抽样样本和时间间隔，实施同样的蒙特卡洛模拟过程，能

够得到两个模拟估计值，以第二种衍生证券真实值与估计值之间的差异作为控制变量，最后得到第一种衍生证券的蒙特卡洛估计值。

假定1V 是需要估汁的第一种衍生证券的价值，

2V 是价值容易估计的第二种衍生证券的价值，第一种证券与第二种证券相似，而1V ˆ和2

V ˆ分别是第一种衍生证券和第二种衍生证券在同样的随机抽样样本的蒙特卡洛估计值，那么利用控制变量技术得到第一种衍生证券的价格估计值为

)V ˆV (V ˆV ˆ2

21Cl 1-+= 这里2

2V ˆV -就是控制变量，它实际上是第一种衍生证券的蒙特卡洛模拟的估计误差，且上述方程的方差之间的关系为 )V ˆ,V ˆcov(2)V ˆvar()V ˆvar()V ˆvar 2

121Cl 1-+=（ 如果)V ˆ,V ˆcov(2)V ˆvar(2

12<，一定有 )V ˆvar()V ˆvar(1

Cl 1< 因此，当两种衍生证券的协方差很大时，或者当两种衍生证券的价格高度相关时，上述关系是成立的，两种衍生证券的正相关性越强，估计效率越理想。然的从实际应用的角度看，这种控制变量技术的应用十分有限，因此，下面是更一般的控制变量技术，其控制变量的形式为

)V ˆV (V ˆV ˆ2

2t 1-β+=β

方差为

)V ˆ,V ˆcov(2)V ˆvar()V ˆvar()V ˆvar(2

12211β-β+=β 这是关于控制变量系数β的二次三项式，下面的目标是能够找到特殊的β使

方差)V ˆvar(1β最小.这时只要取)V ˆvar()V ˆ,V ˆcov(2

21*=β就可以保证方差)V ˆvar(1β最小，这种控制变量技术的缺点是*β需要提前知道协方差)V ˆ,V ˆcov(2

1的信息，而这一般需