用拟蒙特卡罗模拟的欧式看涨期权的定价.pdfVIP

使用拟蒙特卡罗模拟的欧式看涨期权的定价 汪 东 ，张为黎2 (1．上海交通大学安泰管理学院，上海200052；2．上海市发展和改革委员会，上海200003) 摘【 要】传统泉特卡罗模拟使用伪随机数序列，而拟泉特卡罗模拟采用的是完全确定的拟随机数序列(又称低 差异数序列)进行模拟。本文对比了低差异序列与伪随机数序列的统计特点，应用两种模拟方法对欧式期权的价格进行 了模拟计算。实验结果显示拟泉特卡罗模拟在计算精度高于传统泉特卡罗模拟，并且计算速度也更快。 关【 键 词】拟泉特卡罗；泉特卡罗模拟；低差异数序列；期权 【中图分类号】F224．0 文【献标识码】A 文【章编号】1O04—2768(2O04)07—01O9_—02 一 、 介绍 随机数 (低差异数序列)的质量对于模拟的结果有非常重要的影 拟蒙特卡罗模拟(Quasi—MonteCarlosimulation)是采用拟 响。对于模拟人员来说，了解其伪随机数发生器的好坏是很重 随机数序列代替伪随机数的蒙特卡罗模拟。这些随机数是实际 要的。 问题中需要模拟的概率分布的代表样本。 伪随机数序列看起来是不可预测的，就像随机数一样。然 拟随机序列也被称为低差异序列 (1ow—discrepancyse． 而，伪随机数是由确定的算法生成的，因此前面用了一个形容词 quenc~)，实际上，低差异数是完全确定的，因此 “拟随机数”这种 “伪”。Excel软件中的函数rnad()就是一个常用的随机数发生 叫法是不严格的，不过本文仍使用这种叫法。在一些应用中采 器，它使用的算法是线性同余法，很多模拟中都会用到它。 用低差异数进行模拟，可以提高蒙特卡罗模拟的效果，使计算精 更可靠的随机数发生器可以产生与随机数序列几乎一样的 度更高，计算时间更少。本文将通过一个欧式期权的定价的例 伪随机数序列，这里说的一样是指相关统计性质的一样。但是 子来验证这一点。 在许多情况下，用它们来做蒙特卡罗模拟却在主要问题上收敛 进行蒙特卡罗模拟的主要步骤是：第一步在 [0，1]区间上生 得太慢。金融、实物期权以及其它领域中的主要问题是求一个 成均匀分布的伪随机数，拟蒙特卡罗模拟要生成低差异数序列； 函数的(高维)积分。例如，期权得益函数中自变量的取值服从 第二步通常是用逆变换将均匀分布的随机数转换成正态分布的 于一些概率密度函数。在这类应用中，可以牺牲一些实际问题 随机数；第三步是计算出一个模拟出的期权价格；最后是将 以上 中不重要的统计性质(比如独立性)，换取能更快速地收敛到正 的模拟过程反复进行n次，计算出n次模拟期权价格的平均价 确值的性质。 值。根据大数定律，当n趋于无穷大时，平均价格等于期权的价 差异性是对点的分布均匀性的一种测度，拟随机数具有低 值。 差异性的特点。要代表均匀分布很需要低差异性。在多维情况 本文在第二部分分析低差异数序列的性质，并将它与传统 下，它测度点与点之间有没有大的差距以及大的聚集现象。 的伪随机数序列的统计性质进行对比，另外也介绍了如何产生 表 1解释了拟随机数序列的良好统计性质，比较了低差异 这些低差异数序列。第三部分介绍产生正态分布的Moro逆变 数序列与理论的[0，1]均匀分布以及两个典型的由EXCEL生 换。第四部分是以一个欧式看涨期权为例，实证检验拟蒙特卡 成的伪随机数序列的几个常用的统计指标。 罗模拟方法的计算精度和计算速度，结论部分指出实际应用中 裹 1 低差异数序列、均匀分布、伪随机数库列的统计性质的比较 拟蒙特卡罗模拟方法面临的问题。 低差异数序列 f0．1]均匀分布 伪随机数序列 1伪随机数序列2