全部版块 我的主页
论坛 经济学人 二区 外文文献专区
843 26
2022-06-01
英文标题:
《Multilevel estimation of expected exit times and other functionals of
  stopped diffusions》
---
作者:
Michael B. Giles, Francisco Bernal
---
最新提交年份:
2018
---
英文摘要:
  This paper proposes and analyses a new multilevel Monte Carlo method for the estimation of mean exit times for multi-dimensional Brownian diffusions, and associated functionals which correspond to solutions to high-dimensional parabolic PDEs through the Feynman-Kac formula. In particular, it is proved that the complexity to achieve an $\\varepsilon$ root-mean-square error is $O(\\varepsilon^{-2}\\, |\\!\\log \\varepsilon|^3)$.
---
中文摘要:
本文提出并分析了一种新的多层蒙特卡罗方法,用于估计多维布朗扩散的平均退出时间,以及通过费曼-卡克公式对应于高维抛物型偏微分方程解的相关泛函。特别地,证明了实现$\\varepsilon$均方根误差的复杂度为$$O(\\varepsilon ^{-2}、\\124;\\!\\ log\\varepsilon ^ 3)$。
---
分类信息:

一级分类:Mathematics        数学
二级分类:Numerical Analysis        数值分析
分类描述:Numerical algorithms for problems in analysis and algebra, scientific computation
分析和代数问题的数值算法,科学计算
--
一级分类:Quantitative Finance        数量金融学
二级分类:Computational Finance        计算金融学
分类描述:Computational methods, including Monte Carlo, PDE, lattice and other numerical methods with applications to financial modeling
计算方法,包括蒙特卡罗,偏微分方程,格子和其他数值方法,并应用于金融建模
--

---
PDF下载:
-->
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2022-6-1 14:18:31
s-Toping DiffusionSM的预期退出时间和其他函数的多级估计。B、 Giles F.Bernalstember 2018年9月5日摘要本文提出并分析了一种新的多层蒙特卡罗方法,用于估计多维布朗扩散的平均退出时间,以及通过Feynman-Kac公式对应于高维抛物型偏微分方程解的相关泛函。特别地,证明了实现ε均方根误差的复杂性为O(ε-2 | logε|)。(将发表在《SIAM/ASA不确定性量化杂志》上)1简介在本文中,我们关注随机微分方程(SDE)dXt=a(Xt,t)dt+b(Xt,t)dWt的解,0<t≤ T、 (1)具有确定性初始条件X=X。我们假设Xt∈ RdandW={Wt:t≥0}是标准的d′维布朗运动。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2022-6-1 14:18:34
此外,我们假设a:Rd×[0,T]→ Rd和b:Rd×[0,T]→ Rd×Rd′是Lipschitz连续的,Lipschitz常数为La,Lb,例如Ka(x,t)- a(y,s)k≤ La(kx-yk+| t-s |),kb(x,t)- b(y,s)k≤ 磅(kx-yk+| t-s |)),(x,t),(y,s)∈ Rd×[0,T],并且b bt处处都是非奇异的。给定一个单连通有界开域D Rd和固定点X∈D、 我们感兴趣的是估计预期的退出时间Ex,0[τ],其中Ex上的su ffex,t对于时间t:0≤t型≤T表示这是路径Xs上的期望条件,s>T,从Xt=x开始∈ D、 相应的退出时间τ定义为τ=min(T,inf{s>T:Xs/∈D} )。此外,我们有兴趣估算u(x,0),它由u(x,t)=Ex,t定义ZτtE(t,s)f(Xs,s)ds+E(t,τ)g(Xτ,τ), (2) 式中,e(t,t)=exp-ZttV(Xt,t)dt.我们假设函数f:Rd×[0,T]→ R、 g:Rd×[0,T]→ R、 V:Rd×[0,T]→ R、 所有Lipschitz c是否分别与Lipschitz constantsLf、Lg、lv以及g连续∈C2,1(D×[0,T]),及其空间HessianHgand时间导数˙g≡g级/t都有界。Feynman-Kac公式[1 5,16]确定,如果D发行充分平滑,u(x,t)满足PDEut+Xjaj(x,t)uxj+xj,k,lbj,k(x,t)bl,k(x,t)uxj公司xl码-V(x,t)u(x,t)+f(x,t)=0,(3)在域D×[0,t]内,当x∈D、 t=t或x∈D、 0<t<t。当对SDE使用Euler-Maruyama离散化时,对于大小为h的统一时间步长,标准分析给出了一个O(h1/2)强误差,并对路径近似进行了修正,Gobet和Menozzi[11]证明了预期停止时间的弱误差也是O(h1/2)。要获得ε的均方根误差,需要O(ε-2) 路径样本,h=O(ε),因此每个路径样本的平均成本为O(ε-2).
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2022-6-1 14:18:44
因此,总计算成本为O(ε-4).O(h1/2)弱错误是由于SDE溶液中O(h1/2)在每个时间步内的移动;Euler-Maruyama离散化的标准分段常数插值无法解释这一点。一个可能的改进是使用布朗桥插值函数[9],它允许在每个时间步内对平滑边界的最小距离进行采样,从而将弱阶改进为O(h)。Gobet和Menozzi【11】开发的另一种方法基于Glasserman、Broadie和Kou【3】的原始想法,引入了边界校正,将边界设置在法线方向上,距离为o(h1/2)。这也改善了弱err或too(h),因此总体计算成本降低到O(ε-3).Primozic在Giles的障碍选项多层蒙特卡罗(MLMC)方法[6]的基础上,开发了一种MLMC算法,用于估算具有平面边界的一维差异的平均退出时间。这使用了Milstein离散化,该离散化给出了O(h)强误差,并结合了每个时间步内跨越边界概率的Brownian Bridge估计。这给出了一个O(ε-2) 计算复杂性,但其对多维应用的一般化仅限于以下情况:基本SDE满足交换性条件,这意味着在实施Milstein离散化时不需要每个时间步的L'evy区域[9]。Higham等人[13]提出了一种基于EulerMaruyama离散化的MLMC方法,并证明其复杂性为O(ε-3 | logε|)。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2022-6-1 14:18:47
如果他们使用前面提到的两种方法中的任何一种来制作最弱的一阶误差,这将提高到O(ε-2.5 | logε|)。本文的目的是改进Higham等人分析的数值算法,将计算复杂度降低到O(ε-2 | logε|)。这是通过解决【13】分析中确定的问题来实现的,即多级校正方差的衰减很差,如h→ 这大约是O(h1/2),因为在边界O(h1/2)范围内开始的路径有O(h1/2)的概率,在退出do main之前将持续isO(1)的时间。我们通过在另一个退出域后将细路径或粗路径模拟分割成多个独立副本来解决这个问题。这些子路径上的平均值近似于条件期望值,并将多级方差减少到大约O(h)。然而,每个层级上的期望值没有改变,因此位于MLMC中心的伸缩求和仍然有效。本文首先用一个引理来限定FeynmanKac泛函的方差,然后回顾了SDE路径和路径泛函的Euler-Maruyama近似、改进弱收敛性的Gob et&Menozzi技术以及Higham等人使用的MLMC算法。然后提出并分析了新的alg算法。数值分析在很大程度上依赖于Gobet&Menozzi[11]的理论结果,并且与Higham et al[13]的分析结构相似。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2022-6-1 14:18:51
通过数值实验证明了新算法的有效性,并在最后的结论中讨论了未来的研究方向。2退出时间和Feynman-Kac功能如果我们定义uexit(x,t)对应于fexit(x,t)的特定Feynman-Kac功能≡ 0,gexit(x,t)≡ t、 Vexit(x,t)≡ 0,则uexit(x,0)等于之前定义的预期退出时间Ex,0[τ]。我们现在做出一个关键假设。假设1:存在Lipschitz常数Lu,Lexitsuch,即| u(x,t)- u(y,s)|≤ Lu(kx-yk+| t-s |),uexit(x,t)- uexit(y,s)|≤ Lexit(kx-yk+| t-s |)),(x,t),(y,s)∈ D×[0,T]。注释:如果边界D是非常光滑的,假设1中u(x,t)和uexit(x,t)的Lips-chitz条件遵循假设a、b、f、g、V的Lipschitz性质。然而,如果该领域内存在重新进入的角落,则这一假设可能无法满足。下面的引理限制了SDE泛函的方差lpt=ZτtE(t,s)f(Xs,s)ds+E(t,τ)g(Xτ,τ)。Euler-Maruyama离散化的等价理论将在以后的数值分析中发挥重要作用。引理2.1。存在常数c,对于任何x∈ D、 0个≤t<TVx,t[磅]≤ c Ex,t[τ-t] 。注:Vx,tre表示以Xt=x为条件的方差。通过I t^o演算,dE(t,s)g(Xs,s)=E(t,s)-V g+˙g+(g) Ta+跟踪(bTHgb)ds+(g) Tb数据仓库,带a、b、g、˙g≡ g级/t,g、 g的空间梯度和Hg,g的黑森值,都在(Xs,s)处计算。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群