因此，通过对敲打式和欧佩恩式索赔进行定价，我们也可以对敲打式s型索赔进行定价。在时间t，索赔的价值v与支付f（2.1）≤ T由vt=1{τ>T}u（T，Xt，Yt），u（T，x，y）：=E给出{τ>T}Д（XT）| XT=x∈ 一、 Yt=y, （2.3）在温和条件下，（2.3）中定义的函数u是KolmogorovBackward方程0=(t+A）u，u（t，·）=Д，（2.4），其中A，即（X，Y）的发生器，由A=u（X，Y）明确给出x+σ（x，y）x+c（x，y）y+g（x，y）y+ρσ（x，y）g（x，y）x个y、 定义为作用于两次可微且满足一定边界条件dom（a）：={g的函数∈ C： 林克斯→Ig（x，y）=0}。这里我们使用旋转I表示I的最终终点。例如，如果I=（L，∞), 然后A满足limxLg（x，y）=0的函数。通过这个例子，我们假设（2.4）存在唯一的经典解。我们的目标是找到PDE（2.4）的解决方案u。由于一般系数（u，σ，c，g）不存在（2.4）的显式解，因此我们都会寻求u.3形式渐近展开式的显式近似。在本节中，我们将给出u的形式渐近展开式。为了便于理解，让我们引入一些符号。对于A的任何系数，我们定义χε（x，y）：=χ（\'x+ε（x–\'x），\'y+ε（y–\'y）），χ∈ {u，σ，c，g，σg}，其中（\'x，\'y）是一个固定点，ε∈ [0，1]。接下来，我们引入一个算子Aε，ρ，由Aε显式给出，ρ=με（x，y）x+（σ）ε（x，y）x+cε（x，y）y+（g）ε（x，y）y+ρ（σg）ε（x，y）x个y、 其中d om（Aε，ρ）：=dom（A）。现在，考虑一系列PDE问题，指数为（ε，ρ）0=(t+Aε，ρ）uε，ρ，uε，ρ（t，·，·）=Д。（3.1）注意到Aε，ρ|ε=1=A，根据（2.4）和（3.1）得出，uε，ρ|ε=1=u。

因此，我们不是直接寻求PDE问题（2.4）的近似解，而是通过将uε、ρ展开为ε和ρ的幂来寻求PDE问题（3.1）的近似解，如下所示uε、ρ=∞Xi=0∞Xj=0εiρjui，j，（3.2），其中函数s（ui，j）（目前）未知。一旦我们获得了uερ的近似值，我们将通过设置ε=1来获得u的近似值。假设A中的系数是解析的。我们稍后将看到，我们获得的u doe的近似值不需要这个假设。然而，将此作为假设简化了表示，因此我们将暂时继续。由于A的系数是解析的，因此wehaveAε，ρ=∞Xn=0εnAn，0+ρAn，1, （3.3）An，0=un（x，y）x+（σ）n（x，y）x+cn（x，y）y+（g）n（x，y）y、 （3.4）An，1=（σg）n（x，y）x个y、 （3.5）其中d om（A0,0）：=dom（A），并引入了符号χn（x，y）：=n！dndεnχε（x，y）ε=0=nXi=0i'xn–i'yχ（'x，'y）i！（n–i）！（x–'x）i（y–'y）n–i，（3.6）表示χ∈ {u，σ，c，g，σg}。观察到χ是关于点（\'x，\'y）的χ的泰勒级数展开中的n阶项。将展开式（3.2）和（3.3）插入到偏微分方程问题（3.1）中，并收集ε和ρ具有相同幂的项，我们得到o（1）：0=(t+A0,0）u0,0，u0,0（t，·，·）=Д，（3.7）O（εnρk）：0=(t+A0,0）un，k+nXi=0kXj=0（1–δi+j，0）Ai，jun–i，k–j，un，k（t，·，·）=0。（3.8）为清楚起见，我们明确给出了最低阶项he reO（ε）：0=(t+A0,0）u1,0+A1,0u0,0，u1,0（t，·，·）=0，（3.9）O（ρ）：0=(t+A0,0）u0,1+A0,1u0,0，u0,1（t，·，·）=0，（3.10）O（ε）：0=(t+A0,0）u2,0+A2,0u0,0+A1,0u1,0，u2,0（t，·，·）=0，O（ερ）：0=(t+A0,0）u1,1+A1,1u0,0+A1,0u0,1+A0,1u1,0，u1,1（t，·，·）=0，O（ρ）：0=(t+A0,0）u0,2+A0,1u0,1+A0,2u0,0，u0,2（t，·，·）=0。上述计算激发了以下定义。定义3.1。设u为偏微分方程问题（2.4）的唯一经典解。

我们将uρN定义为uρN（t，x，y）：=NXi=0iXj=0εjρi–juj，i–j（t，x，y）（\'x，\'y，ε）=（x，y，1），（3.11），其中u0,0满足（3.7）和un，k满足（3.8）（n，k）6=（0，0）。备注3.2。观察到，我们在（3.11）中设置了ε=1，因此，参数εp在‘uρN’的定义中不起作用。事实上，ε的引入仅仅是为了在上述形式的渐近展开中解释ol。由于ε未出现在原始PDE问题（2.4）中，因此不应出现在u的近似值中。备注3.3。请注意，我们在（3.11）中设置了（(R)x，y）=（x，y）。这通常是一个混淆点，我们希望澄清应该如何处理。首先，应解决嵌套PDE问题（3.7）–（3.8）的顺序，并固定（\'x，\'y）。为了明确起见，让我们将O（εnρk）PDE的解表示为u'x，'yn，k。如果对点（x，y）处的u的近似值感兴趣，则应计算和（3.11）中的u'x，'yn，k（x，y）|（'x，y）=（x，y）。选择（\'x，y）=（x，y）的原因如下。扩散的小时间行为主要由扩散起点（x，y）附近的扩散系数几何结构决定。反过来，点（x，y）附近任何函数的最精确泰勒级数展开是以（\'x，\'y）=（x，y）为中心的泰勒级数展开。备注3。4、如前所述，不需要分析A的系数。索引d，构造n阶近似值uρNone只需要k的运算符Ak，jk≤ N、 因此，u的N阶近似仅要求A的系数为CN。4显式表达式在本节中，我们提供了计算uρN（u的N阶近似值）所需函数（un，k）的显式表达式。我们首先回顾了Duhamel原理。

设Γ0,0为p元算子的基本解(t+A0,0）。即0=(t+A0,0）Γ0,0（·，·，·；t，ξ，η），Γ0,0（t，·，·；t，ξ，η）=Δξ，η。（4.1）Duhamel原理指出，唯一经典解为0=(t+A0,0）u+f，u（t，·，·）=h，由u（t，x，y）=P0,0（t，t）h（x，y）+ZTtds P0,0（t，s）f（s，x，y）给出，其中我们引入了P0,0由A0,0生成的半群，定义如下P0,0（t，s）Д（x，y）=ZIdξZRdη0,0（t，x，y；s，ξ，η）Д（ξ，η），（4.2）其中0≤ t型≤ s≤ T、 提案4.1。设函数（un，k）为DE问题嵌套序列（3.7）–（3.8）的唯一类解。然后，省略空间参数（x，y）以简化符号，函数u0,0由u0,0（t）=P0,0（t，t）Д给出，其中P0,0在（4.2）中定义，对于（n，k）6=（0，0），我们有n，k（t）=n+kXj=1XIn，k，jZTtdsZTsds··ZTsj–1dsj（4.3）P0,0（t，s）An，kP0,0（s），An，k··P0,0（sj–1，sj）Anj，kjP0,0（sj，t）Д，带in，k，jgiven byIn，k，j=n、 ····，njk，····，kj！∈ Z2×j+n+····+nj=n，k+···+kj=k，1≤ ni+ki，适用于所有1≤ 我≤ j.证据见附录A。为清楚起见，我们在此给出了最低阶项u1,0（t）=ZTtdsP0,0（t，s）A1,0P0,0（s，t）Д，u0,1（t）=ZTtdsP0,0（t，s）A0,1P0,0（s，t）Д，u2,0（t）=ZTDSZTSDSP0,0（t，s）A1,0P0,0（s，t）Д+ZTTDS0,0（t，s）A2,0P0,0（s，t）Д，u1,1（t）=ZTtdsZTsdsP0,0（t，s）A0,1P0,0（s，s）A1,0P0,0（s，t）Д+ZTtdsZTsdsP0,0（t，s）A1,0P0,0（s，s）A0,1P0,0（s，t）Д+ZTtdsP0,0（t，s）A1,1P0,0（s，t）Д，u0,2（t）=ZTtdsZTsdsP0,0（t，s）A0,1P0,0（s，s）A0,1P0,0（s，t）Д+ZTtdsP0,0（t，s）A0,2P0,0（s，t）Д。为了进一步进行，我们必须明确指定半群P0,0的作用。我们将考虑这些独立案例：欧盟ropean索赔、单一障碍索赔和双重障碍r索赔ms.4.1欧洲索赔在本节中，我们考虑案例I=r。

Asτ=∞ 当I=R时，我们从（2.1）中可以看出，这个案例对应于写在X上的欧洲主张。我们从下面的引理开始。引理4.2。设H为以下线性运算符H=bx+ax、 dom（H）=C（R）。下面的holdsHψω=λωψω，ω∈ R、 ψω（x）=√2πeiωx，λω（b，a）=biω–aω。（4.4）此外，我们有hψω，ψγi=δ（ω–γ），hf，gi：=ZRdxf（x）g（x），（4.5），其中δ（ω–γ）是Diracδ函数，f表示f的复共轭。证据这个引理可以很容易地通过直接计算进行检验。提案4.3。设P0,0是由dom（A0,0）=C（R）的A0,0生成的半群。然后我们有p0,0（t，t）f=ZRdωZRdγe∧ω，γ（t–t）hψω，γ，f iψω，γ，（4.6）hf，gi：=ZRdxZRdyf（x，y）g（x，y），其中ψω，γ和∧ω，γ由ψω给出，γ（x，y）=ψω（x）ψγ（y），λω，γ=λω（b，a）+γ（b，a）=u，（σ）, （b、a）=c、 （g）,（4.4）中定义的ψω和λω（b，a）。证据使用引理4.2和ZRdωψω（x）ψω（y）=δ（x–y），可以通过直接计算检查Γ0,0（t，x，y；t；ξ，η）：=ZRdωZRdγe∧ω，γ（t–t）ψω，γ（x，y）ψω，γ（ξ，η），（4.7）满足（4.1），因此是(t+A0,0）。表达式（4.6）后面直接插入（4.7）到（4.2）。现在，从命题4.1中，我们看到（un，k）a是形式的项之和=jYi=1ZTsi–1dsiP0,0（si–1，si）Ani，kiP0,0（sj，T）Дs=T.（4.8）使用（4.6），我们可以写出这些术语asA=ZRdωj+1ZRdγj+1jYi=1ZTsi–1dsiZRdωiZRdγie∧ωi，γi（si–si–1）hψωi，γi，Ani，kiψωi+1，γi+1ihψωj+1，γj+1，Дie∧ωj+1，γj+1（T–sj）ψω，γ。（4.9）虽然多重积分看起来很难实现，但当我们计算元素SHψωi，γi，Ani，kiψωi+1，γi+1i时，我们会看到除了单个积分以外的所有积分都会崩溃，这些元素出现在（4.9）引理4.4中。假设ψω、γ和h·，·i在Propo sition4.3中定义。定义运算符B：=xiyjkx公司ly，（4.10）其中i，j，k，l∈ Z+。然后我们有hψω′，γ′，Bψω，γi=（iω）k（iγ）l（–iω） i（–iγ） jδ（ω–ω′）δ（γ–γ′）。证据

证明是一个简单的计算。回顾ψω，γ=2πeiωx+iγy，我们有hψω′，γ′，Bψω，γi=ZRdxZRdyψω′，γ′（x，y）xiyjkx公司lyψω，γ（x，y）=（iω）k（iγ）lZRdxZRdyψω′，γ′（x，y）xiyjψω，γ（x，y）=（iω）k（iγ）l（–iω） i（–iγ） jZRdxZRdyψω′，γ′ψω，γ（x，y）=（iω）k（iγ）l（–iω） i（–iγ） jδ（ω–ω′）δ（γ–γ′）。其中，在最后一步中，我们使用了（4.5）。备注4.5。Dirac d e lta函数的导数定义如下：ZRdωf（ω）nωδ（ω–ω′）=ZRdωδ（ω–ω′）（-ω） nf（ω）=（–ω′）nf（ω′），其中我们通过部分进行积分。从（3.4）、（3.5）和（3.6）中注意，运算符（An，k）是形式为（4.10）的运算符之和。因此，根据引理4.4和备注4.5，我们可以看到关于ωi和γi的积分i n（4.9），对于i=1，2。j由于Dirac delta函数而崩溃。由于p ayo fff函数Д不依赖于y，因此关于γj+1的积分也会崩溃。对于i=1，2，…，r e spect to sif的迭代积分。j只涉及指数，且始终可以显式计算。因此，关于ωj+1的积分是什么，通常必须通过数值计算（如果ν（x）=xnepx或一些n∈ Z+和p∈ R、 然后可以解析地计算关于ωj+1的积分）。4.2单屏障索赔在本节中，我们考虑案例I=（L，∞), 对应于写在X上的单个屏障kn锁定请求，屏障L<X。情况I=（–∞, U） 对于U>X，可以进行类似处理。Webegin具有以下引理。引理4.6。设H为以下线性运算符H=bx+ax、 dom（H）={f∈ C（I）：极限→Lf（x）=0}，I=（L，∞).下式为hηω=uωηω，ω∈ R+，ηω（x；b，a）=Rπe–bx/（2a）sinω（x–L）, uω（b，a）=–b4a–aω。（4.11）此外，我们有hηω，ηγi=δ（ω–γ），hf，gi：=ZIdx f（x）g（x）m（x），m（x）=ebx/a。（4.12）这里，δ（ω–γ）是狄拉克δ函数。证据

引理可以通过直接计算进行检验。提案4.7。设P0,0由A0,0withdom（A0,0）={f生成的半群∈ C（E）：极限→Lf（x，y）=0}，E=（L，∞) ×R.那么我们有P0,0（t，t）f=ZR+dωZRdγe∧ω，γ（t–t）hψω，γ，f iψω，γ，（4.13）hf，gi：=ZdxZrdyf（x，y）g（x，y）m（x），i=（L，∞),式中，ψω，γ和∧ω，γ由ψω给出，γ（x，y）=ηω（x；b，a）ψγ（y），λω，γ=uω（b，a）+λγ（b，a），（b，a）=u，（σ）, （b、a）=c、 （g）,其中ψγ和λγ如（4.4）所定义，ηω和γω如（4.11）所定义，m如（4.12）所定义。证据使用Lemma4.6，我们通过直接计算检查Γ0,0（t，x，y；t；ξ，η）：=ZR+dωZRdγe∧ω，γ（t–t）ψω，γ（x，y）ψω，γ（ξ，η）m（ξ），（4.14）满足（4.1），因此是(t+A0,0）。表达式（4.13）后面直接插入（4.14）到（4.2）。再次，为了计算（un，k），我们必须检查形式（4.8）的项。使用（4.13），我们写出了这些项asA=ZRdωj+1ZRdγj+1jYi=1ZTsi–1dsiZRdωiZRdγie∧ωi，γi（si–si–1）hψωi，γi，Ani，kiψωi+1，γi+1ihψωj+1，γj+1，Дie∧ωj+1，γj+1（T–sj）ψω，γ，（4.15），其中ψωi，γi和∧ωi，γi如命题4.7所示。注意到Each An，kc可以表示为（4.10）中定义的B形式的算符之和，我们必须计算Hψω′、γ′、Bψω、γi形式的内积。这激发了以下引理。引理4.8。假设ψω、γ和h·、·i在第4.7条中定义，B在第（4.10）条中定义。然后hψω′，γ′，Bψω，γi=（iγ）l（–iγ） jδ（γ–γ′）Cω′，ω，i，k，Cω′，ω，i，k=iXm=0即时消息Li–mC（1）ω′，ω，m，k，（4.16），其中c（1）ω′，ω，m，k=m！πc（e）ω，k（ω+ω′）–m–1–|ω–ω′）–m–1) 罪mπ+m！πc（o）ω，k（ω+ω′）–m–1–符号（ω–ω′）|ω–ω′）–m–1cos公司mπ, （4.17）c（o）ω，k=k+1Xm=1k2m–1（–1）k–mb2a公司k–2m+1ω2m–1，（4.18）c（e）ω，k=kXm=0k2m（–1）k–mb2a公司k–2mω2m。（4.19）证明。

回顾命题4.3中ψ的定义，我们共同计算ψω′，γ′，Bψω，γi=ZIdxZRdy m（x）ψω′，γ′（x，y）xiyjkx公司lyψω，γ（x，y）=（iγ）l（–iγ） jδ（γ–γ′）ZIdx m（x）xiηω′（x）kxηω（x）。（4.20）通过直接计算，我们发现kxηω（x）=rπe–b2axc（o）ω，kcos（ω（x–L））+c（e）ω，ksin（ω（x–L））,其中c（o）ω，kand c（e）ω，kare分别由（4.18）和（4.19）给出。因此，ZIdx m（x）xiηω′（x）kxηω（x）=πZIdx xisinω′（x–L）c（o）ω，kcos（ω（x–L））+c（e）ω，ksin（ω（x–L））=iXm=0即时消息Li–mπZ∞dx xmsin（ω′x）c（o）ω，kcos（ωx）+c（e）ω，ksin（ωx）=iXm=0即时消息Li–mC（1）ω′，ω，m，k，（4.21），其中C（1）ω′，ω，m，kis由（4.17）给出。将（4.21）插入（4.20）给出（4.16）。（3.4）、（3.5）和（3.6）中的注释，运算符An，kar是形式（4.10）的运算符之和。我们从（4.16）中可以看出，（4.15）中对γi，i=1，2，3，······，j的积分将due压缩为Diracdelta函数。由于ν与y无关，关于γj+1的积分也会崩溃。此外，在（4.15）中，关于si，i=1，2，3，····j的迭代积分只涉及sia中的指数，因此可以显式求值。我们只剩下关于ωi的积分，i=1，2，3，·····，j+1，可以用数值计算。4.3双势垒索赔在本节中，我们考虑情况I=（L，U），它对应于写在X上的双势垒敲除索赔，其中势垒L和U满足L<X<U。我们从以下引理开始。引理4.9。设H为以下线性运算符H=bx+ax、 dom（H）={f∈ C（I）：极限→Lf（x）=0，limx→Uf（x）=0}，I=（L，U）。以下为holdsHφl= νlφl, l ∈ Nφl（x；b，a）=rU–Le–bx2asinπl（x–L）U–L, νl（b，a）=–b4a–aπl（U–L）。（4.22）此外，我们还有hφl, φki=δl,k、 hf，gi：=ZIdx m（x）f（x）g（x）。这里，δl,kis是Kronecker del ta函数，m由（4.12）给出。证据

引理可以通过直接计算进行检验。提案4.10。设P0,0由A0,0withdom（A0,0）={f生成的半群∈ C（E）：极限→Lf（x，y）=0，limx→Uf（x，y）=0}，E=（L，R）×R。那么我们有p0,0（t，t）f=∞十、l=1ZRdγe∧l,γ（T–T）hψl,γ、 f iψl,γ、 （4.23）hf，gi：=ZIdxZRdyf（x，y）g（x，y）m（x），I=（L，U），其中ψl,γ和∧l,γ由ψ给出l,γ（x，y）=φl（x；b，a）ψγ（y），λl,γ=νl（b，a）+λγ（b，a），（b，a）=u，（σ）, （b、a）=c、 （g）,（4.4）中给出了ψγ和λγ（b，a），φl和νl在（4.2.2）中给出，m在（4.12）中给出。证据使用引理4.9，我们通过直接计算检查Γ0,0（t，x，y；t；ξ，η）：=∞十、l=1ZRdγe∧l,γ（T–T）ψl,γ（x，y）ψl,γ（ξ，η）m（ξ），（4.24）满足（4.1），因此是(t+A0,0）。表达式（4.23）后面直接插入（4.24）到（4.2）。与欧洲和单屏障情况一样，为了计算函数（un，k），我们必须评估表（4.8）中的术语。使用g（4.23），我们将这些条款写为=∞十、lj+1=1ZRdγj+1jYi=1ZTsi–1dsi∞十、li=1ZRdγie∧li、 γi（si–si–1）hψli、 γi，Ani，kiψli+1，γi+1ihψlj+1，γj+1，Дie∧lj+1，γj+1（T–sj）ψω，γ。（4.25）由于每个An，kis是（4.10）中定义的B形式的运算符之和，我们必须计算hψ形式的项li、 γi，Bψli+1，γi+1i。引理4.11。Letψl,γ和h·，·i在命题4.10中定义，B在（4.10）中定义。

THENψl′,γ′，Bψl,γi=（iγ）l（–iγ） jδ（γ–γ′）Cl′,l,i、 k，（4.26），其中Cl′,l,i、 k=iXm=0即时消息Li–mU–Lπm+1C（1）l′,l,m、 k级l′6=l+ C（2）l,m、 kδl′,l,andC（1）l′,l,m、 k=πm+ΓEm+1c（e）l,keF公司m+1；，m+3；–(l – l′)π（4.27）–πm+ΓEm+1c（e）l,keF公司m+1；，m+3；–(l + l′)π+πm+22（m+2）l′c（o）l,k–lc（o）l,kFm+1；，m+2；–(l – l′)π+πm+22（m+2）l′c（o）l,k+lc（o）l,kFm+1；，m+2；–(l + l′)π,C（2）l,m、 k=πm+2lc（o）l，km+2Fm+1；，m+2；–πl（4.28）–2πm+3lc（e）l，km+4m+3Fm+；，m+；–πlc（o）l,k=rU–Lk+1Xj=1k2j–1（–1）k–jb2a公司k–2j+1πlU–L2j–1，c（e）l,k=rU–LkXj=0k2j（–1）k–jb2a公司k–2jπlU–L2j。这里，Γ分别是Euler-ga-mma函数和Fandefare超几何函数以及正则化超几何函数。证据回顾命题4.10中ψ的定义，我们计算ψl′,γ′，Bψl,γi=ZdxZrdy m（x）ψl′,γ′（x，y）xiyjkx公司jyψl,γ（x，y）=（iγ）l（–iγ） jδ（γ–γ′）ZIdx m（x）xiφl′（十）kxφl（十）。（4.29）我们通过直接计算得出：kxφl（x；b，a）=e–b2axc（o）l,KCO公司lπ（x–L）U–L+ c（e）l,ksin公司lπ（x–L）U–L.因此，我们有zidx m（x）xiφl′（十）kxφl（x）=ZIdx xisinl′π（x–L）U–Lc（o）l,KCO公司lπ（x–L）U–L+ c（e）l,ksin公司lπ（x–L）U–L=iXm=0即时消息Li–mU–Lπm+1Zπdx xmsin(l′x）c（o）l,KCO公司(lx）+c（e）l,ksin公司(lx）=iXm=0即时消息Li–mU–Lπm+1C（1）l′,l,m、 k级l′6=l+ C（2）l,m、 kδl′,l, （4.30）其中C（1）的公式l′,l,m、 kand C（2）l′,l,m、 kare分别在（4.27）和（4.28）中给出。将（4.30）插入（4.29）收益率（4.26）。备注4。12、函数sfandef，出现在C的表达式中l′,l,i、 k，由计算πdx xmsin形式的积分得出(l′x）cos(lx）andRπdx xmsin(l′x）sin(lx）。对于任何一个m，l, l′∈ N、 这些积分等于包含x、正弦和余弦幂的项的有限和（可以通过部分积分看到）。

因此，可以用最小的计算工作量来评估函数Fandefc。（3.4）、（3.5）和（3.6）中的注释，运算符An，kar是形式（4.10）的运算符之和。我们从（4.26）中看到，（4.25）中关于γi，i=1，2，3，·····，j的积分由于Dirac Deltafunction而崩溃。由于ν与y无关，因此关于γj+1的积分也会崩溃。此外，在（4.15）中，关于si，i=1，2，3，·························································································································································································································································································。因此，（4.25）是一个明确的SUM，不需要任何数值积分。备注4.13。抛物算子的基本解(t+A0,0+ρA0,1）可以在我们考虑的所有三种情况（欧洲、单势垒和双势垒）中明确获得。因此，有人可能会奇怪，为什么我们要将操作符A扩展为（x–'x）和（y–'y）的幂，以及sρ的幂（而不是仅扩展为（x–'x）和（y–'y）的幂）。我们展开ρ的i npowers的原因是，在没有展开的情况下，（4.3）中关于s，s，…，的积分，SJ不能在单势垒或双势垒情况下显式计算。因此，通过在ρ中展开，可以避免计算多维数值积分。5准确性结果在本节中，我们确定了欧洲期权正式定价近似值的准确性。为了保证精度，我们先介绍一些加法符号。对于集合E Rd，用Cn表示，1b（E）E上具有全局Lipschitz连续导数且阶数小于或等于n的有界函数的类别。设kf kCn，1b表示∞-n阶f导数的范数。我们也用C–1,1b（E）=L表示∞（E） 我们设置k·kC–1,1b=k·kL∞.

下面的定理证明了在局部随机波动率动力学描述的集合上书写的欧式期权价格的一阶近似值的准确性。定理5.1。考虑情况I=R。支持σ，u，c，g的非负整数N∈CN，1b（R），存在一个正常数M，使得≤ kσkCN，1b，kukCN，1b，kckCN，1b，kgkCN，1b≤ M、 此外，假设∈ Ch–1,1b（R），约0≤ h类≤ 2、那么我们有|（u–'uρ）（t，x，y）|≤ C（T–T）h+1，0≤ t<t，x∈ 一、 y型∈ R、 （5.1）对于N≥ 1，我们有|（u–'uρN）（t，x，y）|≤ C（（T–T）+ρ|）NXi=0 |ρ| i（T–T）N–i+h0≤ t<t，x∈ 一、 y型∈ R、 （5.2）（5.1）和（5.2）中的正常数C仅取决于M、N和kДkCh–1,1b。证据见附录B。建立屏障式索赔的渐近精度仍然是一个悬而未决的问题，原因如下。定理5.1的证明利用了zerothorder欧洲问题的冰核中存在的高斯对称性。这种对称性在单势垒和双势垒情况下都不存在，因此无法应用相同的技术来证明准确性。在下一节中，我们将在几个数值示例中探讨我们对barri-er式索赔的应用程序模拟的准确性。6数值示例6.1 Heston模型在本节中，我们对具有Heston dynamicsHeston（1993）的基础S=Ext进行定价近似。具体而言，我们假设（X，Y）s atis dxt=–Ytdt+pYtdWt，dYt=κ（θ–Yt）dt+δpYtdBt，dhW，Bit=ρdt，其中2κθ≥ δ使得Y过程保持严格的正。在我们的数值实验中，我们考虑了具有以下参数的双势垒淘汰调用s和put fixedxyk Tρκθδ0.62 0.0 4。62 0.083-0.4 1.15 0.04 0.2其中EK代表罢工，T代表到期日。

我们首先考虑看涨期权支付（x）=（ex–eK）+，下限L=0固定，上限U>K变化。我们计算了一小时的零阶和二阶价格近似值uρ和uρ以及“精确”价格u，我们通过蒙特卡罗模拟得到了这些价格。在图1中，我们绘制了零阶和二阶近似值的误差u–‘uρ和u–‘uρ，作为上势垒u的函数。为了了解误差的大小，我们还绘制了图2中的精确价格u作为u的函数。在图3和图4中，我们为看跌期权支付（x）=（eK–ex）+提供了类似的p批次，上势垒U=1固定，同时改变下势垒L<K。我们从图s1和图3中可以看出，对于几乎所有级别的看跌期权和看涨期权，U的近似值比Uρ的近似值更精确。备注6.1。出于以下原因，我们省略了图1和图3中的一阶近似值“uρ”。与| uρ–| uρ–| uρ|相比，差值| uρ–| uρ|很小，因为当支付函数Д仅依赖于x时（在调用和放置支付函数的情况下），我们有u0,1=0。因此，我们近似中的第一个相关校正项出现在u1,1的第二阶。由于A.6.2 CEV模型系数的y依赖性，与第一次校正相比，包括第一次相关校正的影响更大。在本节中，我们对具有方差恒定性（CEV）DynamicBox（1975）的S=Ext进行了定价近似。具体而言，我们假设X满足dxt=–σe2Xt（γ–1）dt+σeXt（γ–1）dWt。其中σ>0，γ>0。我们考虑使用以下参数fixedxk Tσγ0.62 0.6 2 0.083 0.32 0.019的双屏障敲出期权和认沽期权，其中ek表示行使，T表示到期日。我们首先考虑看涨期权支付（x）=（ex–eK）+，下限L=0固定，上限U>K变化。

我们分别计算了零和二阶价格近似值uA和u，以及通过蒙特卡罗模拟得到的“精确”价格u。注意，我们省略了u和u中的su perc criptρ，因为相关性在局部波动设置中不起作用。在图5中，我们将第零阶和第二阶近似值的误差u–‘u’和u–‘uO’绘制为上势垒u的函数。为了了解误差的大小，我们还在图6中绘制了精确的价格u作为u的函数。在图7和图8中，我们为看跌期权提供了类似的p批次，其中上势垒U=1固定，而下势垒L<K。我们从图s5和图7中可以看出，在看涨期权和看跌期权的情况下，二阶近似优于二阶近似。7结论在本文中，我们提出了一种在局部随机波动率环境下欧洲和壁垒式索赔的形式化定价近似。我们为欧式索赔提供了严格的准确度结果。我们还提供了几个数值例子，说明了我们对屏障式索赔的近似方法的准确性和多功能性。未来的研究将集中于将我们的技术扩展到其他路径依赖的派生，如回溯和方差风格的声明。致谢作者感谢Stefano Pagliarani和Andrea Pascucci对本手册的有益反馈。命题证明4.1在本节中，我们提供命题证明4.1。命题证明4.1。我们首先注意到，通过将Duhamel原理应用于（3.9）和（3.10），公式（4.3）适用于（n，k）=（1，0）和（n，k）=（0，1）。接下来，假设一个n归纳假设，对于非负整数n和k，n+k≥ 1公式（4.3）适用于非负整数对（m，j），使得m+j≤ n+k。

定义bn，k：={（i，j）| 0≤ 我≤ n、 0个≤ j≤ k，1≤ i+j≤ b} 。将D uhamel原理应用于（3.8），我们可以看到UN+1，k=ZTtdsP0,0（t，s）n+1Xi=0kXj=0（1–δi+j，0）Ai，jun–i+1，k–j=ZTtdsP0,0（t，s）X（i，j）∈An+k+1n+1，kAi，jun–i+1，k–j=ZTtdsP0,0（t，s）An+1，kP0,0（s，t）Д+X（i，j）∈An+k+1n+1，kn+k-i-j+1Xl=1XIn-i+1，k-j，lZTtdsZTsds··ZTsl-1dslP0,0（t，s）Ai，jP0,0（s，s）A（n-i+1），（k-j）··P0,0（sl-1，sl）A（n-i+1）l，（k-j）lP0,0（sl，t）Д，（A.1），其中（A.1）遵循我们的归纳假设。对（A.1）中的和进行重新排序，我们得到UN+1，k=ZTtdsP0,0（t，s）An+1，kP0,0（s，t）Д+n+kXl=1X（i，j）∈An+k–l+1n+1，kXIn–i+1，k–j，lZTtdsZTsds··ZTsl–1dslP0,0（t，s）Ai，jP0,0（s，s）A（n–i+1），（k–j）··P0,0（sl–1，sl）A（n–i+1）l，（k–j）lP0,0（sl，t）Д。（A.2）接下来，注意in+1，k，l=[（i，j）∈An+k–l+1n+1，k（i nn···nl–1j kk···kl–1.nn···nl–1kk···kl–1.∈ In–i+1，k–j，l–1） 。（A.3）因此，将（A.3）与（A.2）相结合，我们得到UN+1，k=ZTtdsP0,0（t，s）An+1，kP0,0（s，t）Д+n+kXl=1XIn+1，k，l+1ZTtdsZTsds··ZTsl–1dslP0,0（t，s）A（n+1），kP0,0（s，s）A（n+1），k··P0,0（sl–1，sl）A（n+1）l+1，kl+1P0,0（sl，t）^1。重新标记（s、s、s、···、sl) → （s，s，···，s）l+1） 重新编制索引并给定Sun+1，k=ZTtdsP0,0（t，s）An+1，kP0,0（s，t）Д+n+k+1Xl=2XIn+1，k，lZTtdsZTsds···ZTsl–1dslP0,0（t，s）a（n+1），kP0,0（s，s）a（n+1），k··P0,0（sl–1，sl）a（n+1）l，klP0,0（sl，t）Д=n+k+1Xl=1XIn+1，k，LZTL TDSZTSDS···ZTsldslP0,0（t，s）A（n+1），kP0,0（s，s）A（n+1），k···P0,0（sl-1，sl）A（n+1）l，klP0,0（sl，t）~n，对于（n+1，k）情况为（4.3）。（n，k+1）的证明是类似的。B定理5.1的证明在本节中，我们证明定理5.1。我们的策略是根据我们目前的情况调整inLorig et al.（2015a）的渐近精度证明。因此，许多定理5.1所需的命题和引理源自Lorig等人的类似命题和引理。

（2015a）。在本节中，我们将t z=（x，y），\'z=（\'x，y）和ζ=（ξ，η）作为R的元素。还可以为操作符A和An，k引入多索引符号。我们有A=x |α|≤2aα（z）Dαz，An，k=Xα∈Akaα，n（z）Dαz，A={（1，0），（0，1），（2，0），（0，2）}，A={（1，1）}，其中α=（α，α）∈ A.∪ A、 |α|=α+α，Dαz=αz在证明定理5.1之前，我们需要一些初步结果。下面，我们用Γ表示与抛物算子对应的基本解(t+A）。引理B.1。对于任何δ>0和α，β∈ N带β≤ N+2，我们有（z–ζ）αDβzΓ（t，z；t，ζ）≤ C·（T–T）|α|–|β| bΓ（T，z；T，ζ），0≤ t<t，z，ζ∈ I×R，（B.1）和（z–ζ）αDβζ（t，z；t，ζ）≤ C·（T–T）|α|–|β| bΓ（T，z；T，ζ），0≤ t<t，z，ζ∈ I×R，（B.2），其中BΓ是算子的基本解(t+（M+δ）(z+z） ，而C是仅依赖于M、N、δ和|β|的正常数。证据结果（B.1）是（Lorig等人，2015a，引理6.21）。通过检查Kolmogorov正演方程，可以看到线性（B.2）。以下因素也会有所帮助。设a和b为常数，如a，b≥ 1/2。然后，对于0≤ t<TZTtds（t–s）a（s–t）b=ΓE（a+1）ΓE（b+1）ΓE（a+b+2）（t–t）a+b+1，（b.3），其中Γ是Euler gamma函数。提案B.2。在定理5.1的假设下，对于任何多指数β∈ 我们有Dβzu0,0（t，z）≤ C·（T–T）min{h–|β|，0}，0≤ t<t，z，z∈ R、 （B.4）如果N≥ 1，则对于任何n，k∈ N、 1个≤ n+k≤ N、 我们有Dβzun，k（t，z）≤ C·（T–T）n+h–|β|1+| z–| z | n（T–T）–n, 0≤ t<t，z，z∈ R、 （B.5）（B.4）和（B.5）中的常数仅取决于M、N、|β|和kΝkCh–1,1b。证据该证明与（Lorig等人，2015a，引理6.24）的证明是一致的。提案B.3。i的定义≥ 0Aρi：=Ai，0+ρAi–1,1，\'Aρn：=nXi=0Aρi，（B.6）uρn：=nXi=0εiρn–iui，n–iε=1，（B.7），约定A–1,1=0。

然后对于N≥ 0，我们有（u–‘uρN）（t，z）=ZTtdsZRdζ（t，z；s，ζ）NXi=0（A–‘Aρi）uρN–i（s，ζ）。（B.8）证明。我们会证明(t+A）（u–'uρN）+NXi=0A–AρiuρN–i=0，（u–'uρN）（T，·）=0，（B.9），其中（B.8）后面是应用Duhamel的pr incipal。请注意，如果我们显示(t+A）(R)uρN=NXi=0A–AρiuρN–i，（B.10），因为(t+A）u=0且u（t，·）=uρN（t，·）=Д。根据方程式（3.7）、（3.8）、（B.6）和（B.7），我们可以推断(t+A0,0）uρn+nXi=1Aρiuρn–i=0。（B.11）我们现在通过诱导来展示（B.10）。当N=0时，由于uρ=(R)uρ，我们有(t+A）’uρ=A–Aρuρ。现在假设（B.10）适用于N≥ 那么我们在（B.11）中得到(t+A）(R)uρN+1=(t+A）(R)uρN+(t+A）uρN+1=NXi=0A–AρiuN–i+A–AρuρN+1–N+1Xi=1AρiuρN–i+1=N+1Xi=1A–Aρi–1uρN–i+1+A–AρuρN+1–N+1Xi=1AρiuρN–i+1=N+1Xi=0A–AρiuρN–i+1。因此，（B.10）适用于所有N。我们现在可以证明T heorem 5.1。定理5.1的证明。从（B.6）和（B.8）中，我们得到（u–’uρN）（t，z）=NXk=0ZTtdsZRdζ（t，z；s，ζ）A–kXj=0Aj，0+ρAj–1,1N–kXi=0ρiuN–k–i，i（s，ζ）。（B.12）设Taαk（z）是aα（z）的第k个泰勒多项式近似，其约定为Taα–1（z）=0。我们将（B.12）改写为（u–'uρN）（t，z）=NXk=0N–kXi=0ρi+1ZTtdsZRdζa（1,1）–Ta（1,1）k–1（ζ） Γ（t，z；s，ζ）D（1,1）ζuN–k–i，i（s，ζ）+NXk=0N–kXi=0Xα∈AρiZTtdsZRdζaα–Taαk（ζ） Γ（t，z；s，ζ）DαζuN–k–i，i（s，ζ）=NXk=0N–kXi=0ρi+1J（1）i，k+NXk=0N–kXi=0ρiJ（2）i，k，1+J（2）i，k，2, （B.13）其中j（1）i，k：=ZTtdsZRdζa（1,1）–Ta（1,1）k–1（ζ） Γ（t，z；s，ζ）D（1,1）ζuN–k–i，i（s，ζ），J（2）i，k，1：=X |α|≤1ZTtdsZRdζaα–Taαk（ζ） Γ（t，z；s，ζ）DαζuN–k–i，i（s，ζ），J（2）i，k，2：=X |α|=2α6=（1,1）ZTtdsZRdζaα–Taαk（ζ） Γ（t，z；s，ζ）DαζuN–k–i，i（s，ζ）。我们首先考虑J（1）i，k。我们不认为J，0=0，因为根据惯例A–1，1=0。

对于k≥ 1，我们按零件进行整合，以获得|α|=|α|=1，J（1）i，k=–ZTtdsZRdζhDαζa（1,1）–Ta（1,1）k–1（ζ） Γ（t，z；s，ζ）ihDαζuN–k–i，i（s，ζ）i。根据乘积规则a和（B.2），在z=’z处进行评估得出| J（1）i，k |≤ CZTtdsZRdζ| z–ζ| k–1Γ（t，z；s，ζ）DαζuN–k–i，i（s，ζ）+ CZTtdsZRdζ| z–ζ| k | DαzΓ（t，z；s，ζ）|DαζuN–k–i，i（s，ζ）.应用（B.1）和（B.5）得出| J（1）i，k |≤ CZTtdsZRdζbΓ（t，z；s，ζ）（s–t）k–1（t–s）N+h–k–i–11+| z–ζ| N–k–i（T–s）N–k–i≤ CZTtdsh（s–t）k–1（t–s）N+h–k–i–1+（s–t）N–i–1（t–s）h–1iZRdζbΓ（t，z；s，ζ）≤ C（T–T）N+h–i.（by（B.3））类似的参数显示| J（2）i，k，1 |≤ C（T–T）N+h–i+2，| J（2）i，k，2 |≤ C（T–T）N+h–i+1，对于0≤ k≤ N、 当N=0时，J（1）i，k=0，那么由（B.13）我们得到|（u–'uρ）（t，z）|≤ C（T–T）h+1。当N≥ 1，通过（B.13），我们得到了|（u–'uρN）（t，z）|≤ CNXk=0N–kXi=0|ρ| i+1（T–T）N–i+h+|ρ| i（T–T）N+h–i+1≤ C（（T–T）+ρ|）NXk=0N–kXi=0 |ρ| i（T–T）N–i+h=C（（T–T）+ρ|）NXk=0（N–k+1）|ρ| k（T–T）N–k+h≤ C（（T–T）+ρ|）NXk=0 |ρ| k（T–T）N–k+h，这证明了定理5.1。参考Antonov，A.和M.Spector（2012年3月）。sabr模型的高级分析。SSRN。Bates，D.（1988年）。崩溃溢价：不对称过程下的期权定价，应用于德意志马克期货期权。宾夕法尼亚大学工作文件。Bowie，J.和P.Carr（1994，8）。静态简单性。危险Carr，P.、K.Ellis和V.Gupta（1998年）。奇异选项的静态hedg。《金融杂志》53（3），1165–1190。Carr，P.和R.Lee（2009年）。Put调用对称：扩展和应用程序。数学金融1 9（4），523–560。Carr，P.和M.Lorig（2015年8月2日）。针对价格和波动性的壁垒式索赔进行了稳健的复制。ArXiv电子打印。Carr，P.和S.Nadtochiy（2011年）。时间同质差异下的统计套期保值。暹罗金融数学杂志2（1），794–838。Cox，J.（1975）。

期权定价注释I：差异的恒定弹性。未出版的草稿，斯坦福大学。该论文的修订版于1996年由《投资组合管理杂志》出版。Davydov，D.和V.Linetsky（2001年）。cev流程下的定价和决策路径依赖操作。《管理科学》47（7），949–96 5。Faulhaber，O.（2002年2月）。存在随机波动性的双障碍期权定价分析方法。Kaiserslautern技术大学博士论文。Fouque，J.-P.，S.Jai munga l和M.Lorig（2011年）。快速均值回复随机波动率模型中期权价格的谱分解。暹罗金融数学杂志2（1）。Fouque，J.-P.，G.Papanicolaou，R.Sircar和K.Solna（2011年）。股票、利率和信贷衍生品的多尺度随机波动性。剑桥：剑桥大学出版社。Hagan，P.、D.Kumar、A.Lesniewski和D.Woodward（2002年）。管理微笑风险。Wilmott杂志1000，84–108。Heston，S.（1993年）。随机波动性期权的封闭形式解，应用于债券和货币期权。修订版。财务部。螺柱。6（2），327–343。Lipton，A.、A.Gal和A.Lasis（20 14）。本地随机波动率框架下普通和第一代奇异期权的定价：研究和新结果。《定量金融》第14（11）章，1899-1922年。Lorig，M.（2014）。基于多尺度差异的衍生品定价：特征函数展开法。数学金融24（2），331–363。Lorig，M.、S.Pagliarani和A.Pascucci（2015a）。抛物方程的解析展开式。出版：暹罗应用数学杂志。Lorig，M.、S.Pagliarani和A.Pascucci（2015b）。多因素局部随机波动率模型的显式隐含波动率。出现：数学金融。默顿，R.（1973）。理性期权定价理论。

《贝尔经济学和管理科学杂志》4（1），141–183。Pagliarani，S.和A.Pascucci（2012年）。局部挥发模型中过渡密度的解析近似。分欧元。J、 数学。10（1），250–270.0.65 0.70 0.75 0.80 0.85 0.900.00050.00100.00150.00200.0025图e 1：对于第6.1节中考虑的赫斯顿模型，我们将u–uρ（蓝色虚线）和u–uρ（橙色虚线）作为看涨期权上限r u的函数。0.65 0.70 0.75 0.80 0.85 0.900.010.020.030.04图e 2：对于第6.1节中考虑的赫斯顿模型，我们将u绘制为看涨期权上屏障u的函数。0.30 0.35 0.40 0.45 0.50 0.55 0.60-0.0015-0.0010-0.00050.00050.0010图e 3：对于第6.1节中考虑的赫斯顿模型，我们绘制u–“uρ”（蓝色虚线）和u–“uρ”（橙色虚线）作为看跌期权下屏障L的函数。0.30 0.35 0.40 0.45 0.50 0.55 0.600.010.020.030.04图e 4：对于第6.1节中考虑的赫斯顿模型，我们将u绘制为看跌期权下载波L的函数。0.65 0.70 0.75 0.800.00050.00100.00150.0020图e 5：对于第6.2节中考虑的CEV模型，我们绘制了u–'u（蓝色虚线）和u–'u（橙色虚线）作为看涨期权上屏障u的函数。图e 6：对于第6.2节中考虑的CEV模型，我们将u绘制为看涨期权的上屏障u的函数。图e 7：对于第6.2节中考虑的CEV模型，我们将u–‘u（蓝色虚线）和u–‘u（橙色虚线）绘制为看跌期权的下载波L的函数。图e 8：对于第6.2节中考虑的CEV模型，我们将u绘制为看跌期权下载波L的函数。