分数阶BS方程中自由边界的数值计算.pdf

Pure Mathematics 理论数学理论数学,2023,13(9),2666-2677 Published Online September 2023 in Hans.https:/www.hanspub.org/journal/pm https:/doi.org/10.12677/pm.2023.139273 文章引用文章引用:钱玥,黄小涛.分数阶 BS 方程中自由边界的数值计算J.理论数学,2023,13(9):2666-2677.DOI:10.12677/pm.2023.139273 分数阶分数阶BS方程中自由边界的数值计算方程中自由边界的数值计算 钱钱 玥玥，黄小涛，黄小涛 南京航空航天大学数学学院，江苏 南京 收稿日期：2023年8月11日；录用日期：2023年9月12日；发布日期：2023年9月19日 摘摘 要要 本文研究一类分数阶美式期权定价模型的数值计算问题。由于分数阶本文研究一类分数阶美式期权定价模型的数值计算问题。由于分数阶Brown运动能够更好地反映风险资运动能够更好地反映风险资产价格的自相关性和长期记忆性，本文将研究由分数阶产价格的自相关性和长期记忆性，本文将研究由分数阶Brown运动推导得到的期权价格模型，即分数阶运动推导得到的期权价格模型，即分数阶Black-Scholes方程。由于美式期权可以提前执行，该模型的边界条件可归结为分数阶抛物型方程的自由方程。由于美式期权可以提前执行，该模型的边界条件可归结为分数阶抛物型方程的自由边值问题。自由边界问题的困难之处在于需要同时求出期权的最佳执行边界及分数阶边值问题。自由边界问题的困难之处在于需要同时求出期权的最佳执行边界及分数阶BS方程的数值解，方程的数值解，因此本文在美式期权的持有区域使用因此本文在美式期权的持有区域使用-差分格式预估数值解，再利用边界条件及差分格式预估数值解，再利用边界条件及Taylor公式校正得到收公式校正得到收敛的方程数值解和自由边界。敛的方程数值解和自由边界。关键词关键词 分数阶美式期权，自由边值问题，数值计算，分数阶美式期权，自由边值问题，数值计算，-差分格式差分格式 Numerical Calculation of Free Boundaries in the Fractional BS Equations Yue Qian,Xiaotao Huang College of Mathematics,Nanjing University of Aeronautics&Astronautics,Nanjing Jiangsu Received:Aug.11th,2023;accepted:Sep.12th,2023;published:Sep.19th,2023 Abstract This paper studies the numerical computation of a class of fractional-order American option pricing models.Since the fractional-order Brownian motion can better reflect the autocorrela-tion and long-term memory of risky asset prices,this paper will study the option price model derived from the fractional-order Brownian motion,i.e.,the fractional-order Black-Scholes equ-ation.And since American options can be executed in advance,the boundary conditions of this model can be reduced to the free boundary problem of the fractional order parabolic equation.钱玥，黄小涛 DOI:10.12677/pm.2023.139273 2667 理论数学 The difficulty of the free boundary problem is that the optimal execution boundary of the option and the numerical solution of the fractional-order BS equation need to be found simultaneously,so in this paper,the numerical solution is predicted in the holding region of the American option using the-difference formats,and then the converged numerical solution of the equation and the free boundary are obtained by using the boundary conditions and the correction of Taylors formula.Keywords Fractional-Order American Options,Free Edge Value Problems,Numerical Computation,-Difference Scheme Copyright 2023 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 Black-Scholes(BS)期权定价模型为期权定价、风险管理以及其他衍生金融产品的定价提供了一个重要的基础框架。1900 年，Louis Bachelier 的投机理论1率先将金融理论和数学思想融合到一起，提出利用布朗运动刻画股票价格变化的随机模型。Sprekle 2指出了股票价格遵循对数正态分布的假设，极大地推进 Louis Bachelier 1的模型，避免了股票价格为负数的情形。1963 年，Boness 3提出投资者具有风险偏好，并在 Sprekle 模型的基础之上添加了风险因素。1964 年，Paul Samnelson 4按照期权相比原生标的资产具有更高风险等性质证实了更具实用性的欧式看涨期权模型。Black 和 Scholes(1973)在有效市场假说和股票价格遵循几何布朗运动等一系列的假设下，推导出了著名的 BS 期权定价模型5。但是由于资本市场的本质特征和状态都是随机波动的，与传统 BS 期权定价模型的假设并不完全吻合，研究人员凭借对股票市场的观察和探索发现，该模型定价与实际市场价格有着较大的差别。Bojdecki 6和 Tudor 7先后研究了修正后的次分数 Brown 运动，提出了基于次分数 Brown 运动的 BS 方程能够更好地刻画原生标的资产价格的变动趋势。伴随微分方程的分形结构在金融领域内被发现，大部分研究人员开始关注金融领域中的分数阶布朗运动，这种模型能更好地反映金融市场中存在的长记忆效应和非正态分布特征。Wyss 8(2000)最早得出时间分数 BS 方程，2006 年，Cartea 等9给出分数阶跳跃扩散期权定价模型和分数阶障碍期权定价模型，2008 年，Jumarie 10基于分数阶 Taylor 公式以及整数阶 BS 方程的推导过程，得出改进的时间分数阶 BS 方程。与欧式期权不同，美式期权的持有人可以在期权到期之前任何时候行使期权，因此期权持有人在整个期权的有效期内必须时刻注意是继续持有还是立即实施。根据期权的内在价值是大于还是等于实施期权的收益，美式期权定价分成两个区域，继续持有区域和终止持有区域，它们之间的交界面称为最佳实施边界。在美式期权中，自由边界是一个重要的概念，因为它可以影响期权的价格和交易策略。本文针对分数阶美式期权定价方程自由边值问题，构造出一种自由边界的预估校正方法。这种方法的基本思想是：对每一行，综合利用分数阶方程及自由边界满足的条件，对期权价格在自由边界附近做泰勒展开，从而得到自由边界点的近似估计，然后调整网格进行下一行计算。与现有的数值方法相比，本文方法有两个优点：一是求解区域限定在自由边界内，从而减少了计算量；二是在求出期权值的同时Open AccessOpen Access钱玥，黄小涛 DOI:10.12677/pm.2023.139273 2668 理论数学 也求出了期权的最佳执行边界。2.时间分数阶期权定价模型时间分数阶期权定价模型 期权定价理论最早由 Black 和 Scholes 提出，在假设原生资产价格服从几何 Brown 运动，基于随机微分方程理论推导出了整数阶期权定价模型5：222210.2CCCSrSrCtSS+=分数阶布朗运动可以被用于更准确地描述市场中的非局部和非马尔可夫性质，例如长记忆、长程相关性和非对称性等，其模型更应用于现实世界。由此得到期权价格所满足的欧式期权定价模型，即时间分数阶欧式 BS 方程11 12：()()()()()()1221,0.22tSSStVS trD SVrVS V+=其中相关符号定义如下(后续章节继续沿用)：S：风险资产的价格(以简单为计，本文都认为是股票价格)；D：期权在有效期内标的资产的红利(或回报率)；t：时间变量，且0t；r：投资的无风险利率；：标的资产的价格波动率；V：欧式看涨期权的价格；()(),tVS t：(),V S t对 t 的 Riemann-Liouville 分数阶导数，且01。定理定理(时间分数阶Black-Scholes方程)假设股票价格S的运动过程为白噪声下的分数阶布朗运动，则：()()2ddd.SrS tSw tt=+(1)式中，r 为无风险利率(令股票价格的漂移率r=以便于计算)，为期权价格波动率，()w t是与分数阶布朗运动()(),2,01HBtB t=相关的 Guassian 白噪声。设(),V S t是关于股票价格 S 和时间 t 的欧式看涨期权价格函数，且(),V S t关于 t 有 阶导数，关于 S 可导。则有：()()()()()()1221,0.22tSSStVS trD SVrVS V+=证明证明(),V S t作泰勒展开，有()()()()()22211dddddd.12VVVVtSSotSStS=+(2)式中，()()()()()()()()()222212222222dddd2dd.SrS tSw ttr Str S w ttS wtt+=+=+(3)当d0t 时，忽略等式右边阶数高于dt的项，即有()()()2222dd.SS wtt=(4)钱玥，黄小涛 DOI:10.12677/pm.2023.139273 2669 理论数学 将(1)和(4)代入(2)，则()()()()()()()22222211ddddd.12VVVVtrS tSw ttS wttStS=+(5)对上式取均值，得到()()222211ddd.12VVVVStrStStS=+(6)利用对冲建立相应的期权定价连续模型及计算公式。形成投资组合 VS=选取，使得在(),dt tt+内是无风险的。设在时刻 t 形成投资组合，并在时段(),dt tt+内不改变份额。那么由于是无风险的，因此在时刻dtt+，投资组合的回报是 ddttttr t+=由于考虑到支付股息，因此ddddttttttttVS D tS+=，代入上式，有 ddddtttttttVSrtS D t=+利用伊藤公式，同时取tVS=，得到()dddVVVrVrD StSSS=+(7)联立(6)和(7)整理可得考虑红利情况下，时间分数阶 Black-Scholes 微分方程：()()()()()()12210,0,0.22tSSStVS VrD SVrVSSttT+=美式期权定价问题的数学模型一般可归结为抛物型方程自由边值问题。自由边界问题13是指微分方程求解区域的某些边界是未知的(甚至是随时间变化的)，需要和微分方程的解共同确定。Holmes 14和Yang 15提出有限元方法，利用 Front-fixing 变换，将自由边界消除，得到一个矩形求解域，由此化简计算。但一旦网格剖分极其细密，也会致使运算量增加，进而出现计算数据溢出计算机内存的情况，降低计算效率。Janson 和 Tyskc 16提出神经网络法，这种方法是将回归方法和神经网络相结合，在某些情况下会出现较大误差或不符合实际情况的结果。Muthuraman 17提出移动边界法，只要求解一个通常的线性初边值问题，用新的初始边界代替自由边界，重复上述步骤即可求出自由边界并解决期权问题。张铁和祝丹梅18提出变网格差分方法，首先建立一个自由边界所满足的方程，再利用变网格技术求解此问题。Khaliq，Voss 和 Kazm 19 20提出惩罚函数法，但这种方法的缺点是无法求出自由边界。此外还有蒙特卡罗模拟方法21，紧致差分算法22 23等许多算法24-29。至今对时间分数阶 BS 方程数值解的探讨不多，SongLina 等30(2013)结合欧式看跌期权提出时间分数阶 BS 方程的隐式差分算法，2017 年，张瑜等31提出时间分数阶 BS 方程的奇数层显性差分与偶数层隐性差分交替的并行差分的方法。假设市场是风险中性的，在期权有效期内有红利支付。此时，美式看涨期权的分数阶 BS 方程为：钱玥，黄小涛 DOI:10.12677/pm.2023.139273 2670 理论数学 ()()()()()()12210,0,0.22tSSStVS VrD SVrVSSttT+=(8)其边界条件为：()()(),max,0,0,V S TSESST=(9)()0,0,0,VttT=(10)()()()()(),1,0,SV SttStE VStttT=(11)()(),.V S tSE StS=+(12)其中 T 为期权到期日，E 为期权的执行价格，()St是期权提前执行的自由边界，在期权执行日有()*max,rSTEED=。分数阶美式期权边界条件的物理意义体现在投资者在期权价格、时间和市场价格等方面做出行权或继续持有的决策，反映了投资者对期权价值和市场情况的判断，对于期权定价和投资决策具有重要的参考价值。为了便于构造数值方法，首先对问题进行变换，目的是将变系数方程化为常系数方程，将反向时间问题化为正向时间问题，引进变量变换()()()()()()2*21,2ln,ln,21,.fTtStSxkxkkrDEEV S tEu xSE=+=+=+在此变换下，自由边界问题(8)(12)问题等同于如下问题：()()()()()1,0,xxxfuWuhufug xTxx+=(13)()()(),0max 1e,0,0,xfu xxx=(14)()1lim,1 e,0,x kxu xT=(15)()()()()1,0,0,0,fxfu xuxT=(16)()(),0,.fu xxx=+(17)其中，()()()()()()()()()()()11112221212122,222224e,.22x krDrhfWhETrrDg xTT+=+=此时()fx为新的自由边界，并且()()()0max ln,0fxr D=。另外，为方便计算将公式(13)、(14)中的左边界条件限定为xL=。钱玥，黄小涛 DOI:10.12677/pm.2023.139273 2671 理论数学 3.分数阶期权定价方程的数值算法分数阶期权定价方程的数值算法 3.1.方程解的差分近似方程解的差分近似 由于 B-S 方程在时间上是分数阶求导而在空间上是整数阶求导，在此需要用不同的离散化方法。3.1.1.时间与空间上的差分近似时间与空间上的差分近似 方程(1317)数值求解的困难在于必须同时确定持有区域的函数值以及自由边界的位置，而具有特定性质和边界条件的分数阶抛物型偏微分方程，可以通过适当的离散化和迭代求解，由此得出分数阶抛物型偏微分方程的数值解。同时，由于 BS 方程关于空间仅有一维，有限元和差分法等效。差分方法是通过差商代替微商对方程以及定解问题离散化。通过建立与偏微分方程相应的差分方程有多种方式，从求解的方式可以分为两大类：一类是显式差分格式，求解的过程是显式，通过直接运算求出它的值；另一种是隐式差分格式，求解过程必须通过一个代数方程组才能得到它的值32。众所周知，如果()yf x=充分光滑，那么()fx和()fx可以用以下差商形式代替33：()()()()()()()2222,.2f xxf xxf xxf xf xxfffxfxxxxx+=而对于时间分数阶期权定价模型，其中分数阶导数使用 Riemann-Liouville 定义33，时间分数阶导数的离散格式为34 35：()()()()()()()111211,1,2mimimjimjju xdu xu xjjO k+=+忽略误差，可以得到()()()()()()111111,.2mimimjimjmiju xdu xd u xl u x+=其中，()()()111111,211.jjljjdjjj=+将求解区域划分为均匀网格，将时间区间10,T划分成 M 个间隔相等的小区间，记()1,1,2,1.mmdmM M=+将区间(),0fL x划分成 N 个间隔相等的区间段，()1,1,2,1.ixLidxiN N=+=+不妨记(),miimuu x=，边界条件转化为：()()()()100max 1 e,0,1 exp1,0.LidxmmiNuuLk mdu+=建立问题(6)(10)的古典显式格式(见图 1)：()()()()()111112,2.2mmmmmimmmiiiiimmmiiu xuuuuuWhfugdxdx+=+古典隐式格式(见图 2)：()()()()()1111111111111112,2.2mmmmmimmmiiiiimmmiiu xuuuuuWhfugdxdx+=+钱玥，黄小涛 DOI:10.12677/pm.2023.139273 2672 理论数学 Figure 1.Fractional order explicit differ-ence scheme 图图 1.分数阶显式差分格式 Figure 2.Fractional order implicit differ-ence scheme 图图 2.分数阶隐式差分格式 3.1.2.-差分格式差分格式 由于显示差分的条件稳定性和隐式差分的绝对稳定性，在此取加权差分格式36 37：()()()()()()()()()()()()()()1121111112111111111111112,122122,2nimjimjjmmmmmmmiiiiimmmiimmmmmmmiiiiimmmiidu xu xjjuuuuuWhfugdxdxuuuuuWhfugdxdx+=+=+其中01，当1 2=时，上式即为分数阶 BS 方程的 C-N 格式。整理得：()()()()()()()()()()1111111111111111112111211,mmmmmimmimmimmmmmimmimmimmmjmijimiijBAuBCuABuBAuBCuABuRgd ul uRg+=+=+钱玥，黄小涛 DOI:10.12677/pm.2023.139273 2673 理论数学 即()111111111111111.1mmmmimimimmmmmmmjiimimimijii mjcuaubuRRc ua ub ud uu lgg+=+=+其矩阵形式的具体表达式 11122111331111111111111100mmmmmmmmmmmmmmmmmmmmmmNNmmmmmmNNmmmmbcbcUUababUUccabcabcUUababUUa UaU+=+()11222113331111111111111!,1mmmmmmjmjjmmNNNmmmmmNmNNNNUggUggld URRUggc UcUUgg+=+1122222221122333331123111122111111122mjmmmmjmmmmjmjmjmmmNNNNNmjmmmNNNNNUUUUUUUUUUdddddUUUUUUUUUU+=+=+.其中，()()()()()()()()()()()()()()()()()()()2111111111111112,.2,21.1,1,12.1,211.mmmmmmmiimmmmmmmmmmmmmmmmmmmjjRWRhRdABCR fgg xdxdxcBAaABbCBcBAaABbCBljjdjjj+=+=+=+=+=+3.2.自由边界的差分近似自由边界的差分近似 自由边界的确定可以看作是分数阶期权定价模型在边界条件下的最优决策问题。在分数阶抛物型方程中，自由边界的确定涉及到期权价格在边界上的数值解。通过求解分数阶抛物型方程，可以得到期权价格在不同时间和价格水平下的数值解，从而确定自由边界。关于自由边界的数值求解产生了许多数值方法，包括焓方法(E.L.Albasiny,1956)，前沿固定法，有限差分法(J.Crank,1957)和变网格法等。本文主要使用前沿追踪法，建立了变网格，在每一个时间步长，离边界远处用正常的有限差分，在自由边界附近将差分公式进行修正，以适应不均匀的空变化的自由边界。上一节给出了期权定价方程在持有区域内部的离散算法，本节用“预估校正”的方法同时确定自由边界及边界附近的函数值。3.2.1.预估预估 由于函数值和自由边界的位置必须同时确定，因此利用 u 在时刻各点的值来确定自由边界的位置()fx。钱玥，黄小涛 DOI:10.12677/pm.2023.139273 2674 理论数学 设第 m 行值已知，即自由边界()fmx的值已知，现确定 m+1 行函数的值及自由边界的位置。首先假设()()1fmfmxx+=由上节的-差分格式或显隐差分格式可以求出执行区域上的函数值，在此建立方程估计()fx的值。已知()1Nfxx是邻近()fx的一个网格节点，设()1fNxxqdx=+，其中 dx 为网格步长。我们要确定()fx的值，只需要求出 q。由泰勒公式得，()()()223121.2fffNxxxuuuuqdxqdxOqdxxx=+利用边界条件(16)，且舍去高阶小项()()3Odx，得到近似方程()221210.2fNxuuqdxx=(18)再次利用(16)式得到 dd0.ddfffffxxxxxuuuux=+=故由方程(13)得()()(),fxxfxhug x=代入(18)式有()()()21,10.2fNg xuqdxh+=因此可得 q 所满足的非线性方程，()()()()112222122224e0.212Nxqdx kNqdxrrDTu+=+(19)在方程(19)中解 q 存在唯一的基础上，我们可利用牛顿法计算 q：()()()()()()()11112112211,224e.448224eNmmmNmmmmmmxq dx kmmmNmxq dx kmmmmH qqqHqq rrDquH qHqHqrrDqdxDq dxqdx+=+=+其中，()()()()()()()()()()()11122221222211222224e,212448224e.2212NNxqh kNxqh kqdxrrDTH quq dxrrDqdxDqdxqdxHqHqT+=+=+由于执行区域内10Nu，由牛顿迭代可以由方程(19)求出 q，即求出()1fmx+。钱玥，黄小涛 DOI:10.12677/pm.2023.139273 2675 理论数学 如果01q的值时，可根据边界条件设 0,0,1,.mimUiNmM=下图(见图 3)是从到m到1m+步网格扩充示意图，其中黑色小方框点表示扩充的节点。Figure 3.Grid expansion diagram 图图 3.网格扩充示意图 3.2.2.校正校正 当计算第 m+1 行()fx的值时，上节首先假设()()11fmfmxx+=，再由方程得到边界的预估值()211mfmNxxqdx+=+。对新的预估值()21fmx+重复预估过程，得到数列()1,kfmxkZ+，其极限即为自由边界的精确值。给定精度，当()()111kkfmfmxx+时，停止计算，以()1kfmx+表示精确值。重复以上过程即可得到执行区域上分数阶 BS 方程的数值解及自由边界。4.总结总结 由于资本市场的本质特征和状态都是随机波动的，与传统 BS 期权定价模型的假设并不完全吻合，而基于次分数 Brown 运动的 BS 方程能够更好地刻画原生标的资产价格的变动趋势，基于修正后的布朗运动能够更多地反映自相关性、长期记忆性和增量相关性等众多性质，最终得出时间分数阶 BS 方程。本文研究了一种包含分数阶微分方程的美式期权定价方程。该模型更适用于实际金融市场是相对易于求解的期权定价模型。本文给出了该模型在持有区域上的-差分格式，其中离散化的精度和稳定性可以通过选择合适的网格大小和差分近似方法来提高。此外，根据该格式得到了分数阶 BS 微分方程数值解的算法。由于自由边界的问题是一个非线性问题，同样在一般情况下也不容易得不到其显式表达式。本文根据自由边界所满足的条件，通过构建求解美式期权定价自由边值问题的预估校正网格差分方法，同时求出期权的最佳执行边界及分数阶 BS 方程的数值解。基金项目基金项目 国家自然科学基金面上项目支持：NSFC12071215。参考文献参考文献 1 Bachelier,L.(1900)Theory de la speculation.Annals Scientifiques de Leola Normal Superior,17,21-86.钱玥，黄小涛 DOI:10.12677/pm.2023.139273 2676 理论数学 https:/doi.org/10.24033/asens.476 2 Sprenkle,C.M.(1961)Warrant Prices as Indicators of Expectations and Preferences.Yale Economic Essays,1,178-231.3 Boness,A.J.(1964)Elements of a Theory of Stock Option Value.The Journal of Political Economy,72,163-175.https:/doi.org/10.1086/258885 4 Samuelson,P.A.(1965)Rational Theory of Warrant Pricing.Industrial Management Review,6,441-454.5 Black,F.and Scholes,M.(1976)The Pricing of Options and Corporate Liabilities.Science,4,141-183.6 Bojdecki,T.,Gorostiza,L.and Talarczyk,A.(2007)Some Extensions of Fractional Brownian Motion and Sub-Fractional Brownian Motion Related to Particle Systems.Electronic Communications in Probability,12,161-172.https:/doi.org/10.1214/ECP.v12-1272 7 Tudor,C.(2008)Inner Product Spaces of Integrands Associated to Subfractional Brownian Motion.Statistics and Probability Letters,78,2201-2209.https:/doi.org/10.1016/j.spl.2008.01.087 8 Wyss,W.(2000)The Fractional Black-Scholes Equation.Fractional Calculus and Applied Analysis,3,51-61.9 Cartea,A.and Del-Castillo-Negrete,D.(2006)Fractional Diffusion Models of Option Prices in Markets with Jumps.Physica A:Statistical Mechanics and its Applications,374,749-763.https:/doi.org/10.1016/j.physa.2006.08.071 10 Jumarie G.(2008)Stock Exchange Fractional Dynamics Defined as Fractional Exponential Growth Driven by(Usual)Gaussian White Noise.Application to Fractional Black Scholes Equations.Insurance Mathematics and Economics,42,271-287.https:/doi.org/10.1016/j.insmatheco.2007.03.001 11 张秀芝.美式期权的定价方法介绍与比较D:硕士学位论文.济南:山东大学,2012.12 Holmes,A.D.and Yang,H.(2008)A Front-Fixing Finite Element Method for the Valuation of American Options.SIAM Journal on Scientific Computing,30,2158-2180.https:/doi.org/10.1137/070694442 13 Nielsen,B.F.,Shavhaug,O.and Tveito,A.(2002)Penalty and Front-Fixing Methods for the Numerical Solution of American Option Problems.The Journal of Computational Finance,5,69-97.https:/doi.org/10.21314/JCF.2002.084 14 Janson,S.and Tyskc,J.(2003)Volatility Time and Properties of Option Prices.Annals of Applied Probability,13,890-913.https:/doi.org/10.1214/aoap/1060202830 15 Muthuraman,K.(2008)A Moving Boundary Approach to American Option Pricing.Journal of Economic Dynamics and Control,32,3520-3537.https:/doi.org/10.1016/j.jedc.2007.12.007 16 张铁,祝丹梅.美式期权定价问题的变网格差分方法J.计算数学,2008,30(4):379-387.17 Khaliq,A.Q.M.,Voss,D.A.and Kazmi,S.H.K.(2006)A Linearly Implicit Predictor-Corrector Scheme for Pricing American Options Using a Penalty Method Approach.Journal of Banking and Finance,30,489-502.https:/doi.org/10.1016/j.jbankfin.2005.04.017 18 Khaliq,A.Q.M.,Voss,D.A.and Kazmi,K.(2008)Adaptive-Methods for Pricing American Options.Journal of Computational and Applied Mathematics,222,210-227.https:/doi.org/10.1016/j.cam.2007.10.035 19 Zhao,J.C.,Davison,M.and Corless,R.M.(2007)Compact Finite Difference Method for American Option Pricing.Journal of Computational and Applied Mathematics,206,306-321.https:/doi.org/10.1016/j.cam.2006.07.006 20 Tangman,D.Y.,Gopaul,A.and Bhuruth,M.(2008)Numerical Pricing of Options Using High-Order Compact Finite Difference Schemes.Journal of Computational and Applied Mathematics,218,270-280.https:/doi.org/10.1016/j.cam.2007.01.035 21 Han,H.D.and Wu,X.N.(2003)A Fast Numerical Method for the BlackScholes Equation of American Options.SIAM Journal on Numerical Analysis,41,2081-2095.https:/doi.org/10.1137/S0036142901390238 22 Jaillet,P.,Lamberton,D.and Lapeyere,B.(1990)Variational Inequalities and the Pricing of American Options.Acta Applicandae Mathematica,21,263-289.https:/doi.org/10.1007/BF00047211 23 Wilmott,P.,Dewynne,J.and Howison,S.(1993)Option Pricing:Mathematical Models and Computation.Oxford Fi-nancial Press,London.24 Jacka,S.D.(1991)Optimal Stopping and the American Put.Mathematical Finance,1,1-14.https:/doi.org/10.1111/j.1467-9965.1991.tb00007.x 25 Evans,J.D.,Kuske,R.and Keller,J.R.(2002)American Options on Assets with Dividends Near Expiry.Mathematical Finance,12,219-237.https:/doi.org/10.1111/1467-9965.02008 26 Tangman,D.Y.,Gopaul,A.and Bhuruth,M.(2