二次有限体积法定价美式期权.pdfVIP

2015年 2月 计 算 数 学 第 37卷第 1期 Feb．2015 MATHEMATICA NUMERICA SINICA Vo1．37，No．1 二次有限体积法定价美式期权 ) 甘小艇 (同济大学数学系，上海200092； 楚雄师范学院数学与统计学院 云南楚雄 675000) 殷俊锋 (同济大学数学系，上海200092) 摘 要 本文考虑二次有限体积法定价美式期权．构造了隐式欧拉和 Crank—Nicolson两种全离散二 次有限体积格式，并得到相应的线性互补问题．采用基于超松弛迭代的模方法求解线性互补问题， 并与投影超松弛迭代法作数值比较．数值实验结果表明Crank—Nicolson二次有限体积格式的求 解效率高于隐式欧拉格式，模方法的求解速度较快，二次有限体积法的求解精度较高． 关键词：二次有限体积法；美式期权；模方法；投影超松弛迭代；线性互补问题． MR (2000)主题分类：65M06，65M20，65M50，65G40 1．引 言 期权是金融市场中非常重要的金融衍生工具之一，其理论最初 由Black和 Scholes提出， 经过 Heston等人的改进和完善，成为金融风险管理的有效工具之一，吸引了众多国内外学 者和专家的研究 _1J．期权理论的核心问题是期权定价 问题，对于标准的欧式期权，可以利用 Black—Scholes方程得到解析解的定价公式．然而，对于美式期权，或者考虑跳扩散和随机波动 率等复杂的期权模型，几乎无法找到其解析解，即使运用数值方法求解也是极具挑战力的课 题．因此，研究这些美式期权定价问题的数值方法具有重要的理论意义和实用价值 [213]． 最早出现的期权定价方法是格点法，如二叉树方法、三叉树方法等 _4J．格点法简单易行， 算法也容易实现，但收敛速度较慢，计算量随时间方向节点数的增加而剧烈增长且不够精确． 基于偏微分方程的数值方法近年来受到国内外学者的普遍重视和广泛研究，该方法一般先采 用稳定而精确的数值离散方法，譬如有限差分方法 】，有限元方法 _718】等，结合边界条件、 终止条件及价格约束条件，将美式期权定价问题转化为在每个时间层上的线性互补问题，然 后运用投影三角分解法、投影超松弛方法 (PSOR)或基于超松弛迭代的模方法 (MSOR)等进 行数值求解，最后得到美式期权的价格 [9,10J．另外，对于高维偏微分方程的问题，基于Monte Carlo方法的随机模拟方法吸引了很多研究人员的兴趣，但其缺点在于收敛速度较慢 【¨』． 有限体积法又称广义差分法、有限体积元法，由我国学者李荣华教授最先提出_l2_．该方 法具有格式构造简单，保持局部守恒性和精度高等优点，目前已被广泛应用于许多领域 【]． 近年来，该方法也被国内外学者应用于期权定价问题中 [14,15J．Wang等建立了欧式期权定价 的有限体积格式并进行数值求解 _l6J；孙玉东等对永久美式期权定价模型的三次有限体积法 进行了理论探讨 -17_，此外，李广志等对美式期权定价模型的三次有限体积法进行了较深刻的 2014年 3月 4日收到． )基金项 目：国家 自然科学基金 ，云南省应用基础研究计划青年项 目(2013FD045) 68 计 算 数 学 理论研究 [18]，但他们都没有作进一步数值分析．本文深入研究二次有限体积法定价美式期权， 首先构造了隐式欧拉和Crank—Nicolson(C—N)两种全离散二次有限体积格式，并得到相应的线 性互补问题．其次，采用基于超松弛迭代的模方法 (MSOR)求解线性互补问题，并与投影超 松弛迭代法 (PSOR)作比较．数值实验结果表明两种全离散二次有限体积格式都是有效的，而 C—N格式