全部版块 我的主页
论坛 经济学人 二区 外文文献专区
542 34
2022-05-19
英文标题:
《Ninomiya-Victoir scheme: strong convergence, antithetic version and
  application to multilevel estimators》
---
作者:
Anis Al Gerbi, Benjamin Jourdain, Emmanuelle Cl\\\'ement
---
最新提交年份:
2015
---
英文摘要:
  In this paper, we are interested in the strong convergence properties of the Ninomiya-Victoir scheme which is known to exhibit weak convergence with order 2. We prove strong convergence with order $1/2$. This study is aimed at analysing the use of this scheme either at each level or only at the finest level of a multilevel Monte Carlo estimator: indeed, the variance of a multilevel Monte Carlo estimator is related to the strong error between the two schemes used on the coarse and fine grids at each level. Recently, Giles and Szpruch proposed a scheme permitting to construct a multilevel Monte Carlo estimator achieving the optimal complexity $O\\left(\\epsilon^{-2}\\right)$ for the precision $\\epsilon$. In the same spirit, we propose a modified Ninomiya-Victoir scheme, which may be strongly coupled with order $1$ to the Giles-Szpruch scheme at the finest level of a multilevel Monte Carlo estimator. Numerical experiments show that this choice improves the efficiency, since the order $2$ of weak convergence of the Ninomiya-Victoir scheme permits to reduce the number of discretization levels.
---
中文摘要:
在本文中,我们对Ninomiya-Victoir格式的强收敛性感兴趣,该格式具有2阶弱收敛性。我们证明了订单为1/2美元时的强收敛性。本研究旨在分析该方案在每个层次上的使用情况,或仅在多层蒙特卡罗估计量的最精细层次上的使用情况:事实上,多层蒙特卡罗估计量的方差与在每个层次的粗网格和细网格上使用的两个方案之间的强误差有关。最近,Giles和Szpruch提出了一个方案,该方案允许构造一个多级蒙特卡罗估计量,以达到精度$\\ε$\\的最优复杂性$\\左(\\ε^{-2}\\右)$。本着同样的精神,我们提出了一种改进的Ninomiya-Victoir方案,该方案可能与Giles-Szpruch方案在多层蒙特卡罗估计的最佳水平上的1美元订单强耦合。数值实验表明,这种选择提高了效率,因为Ninomiya Victoir格式的弱收敛阶数为2$,可以减少离散化层数。
---
分类信息:

一级分类:Quantitative Finance        数量金融学
二级分类:Computational Finance        计算金融学
分类描述:Computational methods, including Monte Carlo, PDE, lattice and other numerical methods with applications to financial modeling
计算方法,包括蒙特卡罗,偏微分方程,格子和其他数值方法,并应用于金融建模
--
一级分类:Mathematics        数学
二级分类:Probability        概率
分类描述:Theory and applications of probability and stochastic processes: e.g. central limit theorems, large deviations, stochastic differential equations, models from statistical mechanics, queuing theory
概率论与随机过程的理论与应用:例如中心极限定理,大偏差,随机微分方程,统计力学模型,排队论
--

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

扫码加我 拉你入群

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

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

全部回复
2022-5-19 19:18:00
Ninomiya-Victoir格式:强收敛,对偶版本以及在多层估计中的应用。Al Gerbi,B.Jourdain*和E.Cl'ement+2018年10月10日在本文中,我们对Ninomiya Victoirscheme的强收敛性质感兴趣,已知Ninomiya Victoirscheme具有弱收敛的2阶。我们证明了1/2阶的强收敛性。本研究旨在分析该方案在每个层次上的使用情况,或仅在多层蒙特卡罗估计量的最高层上的使用情况:事实上,多层蒙特卡罗估计量的方差与在每个层次的粗网格和细网格上使用的两个方案之间的str on g误差有关。最近,Giles和Szpruch在[6]中提出了一个方案,允许构造一个多级蒙特卡罗估值器,以实现-2.为了精度。本着同样的精神,我们提出了一个改进的Ninomiya Victoir方案,该方案可能与Giles Szpr uch sch-eme在多层蒙特卡罗估计量的最低水平上的1阶强耦合。数值实验表明,这种选择提高了效率,因为Ninomiya-Victoir格式的弱收敛阶2会减少离散化层数。1简介本文致力于计算Y=E[f(XT)],其中f:Rn-→ R是一个Payoff函数,Xt是时间T的解∈ R*+, 到形式的多维随机微分方程dXt=b(Xt)dt+dPj=1σj(Xt)dWjt,t∈ [0,T]X=X.(1.1)此处,X∈ Rn是初始条件,W=W西部数据是d-维标准布朗运动,b:Rn-→ Rn是漂移系数,σj:Rn-→ Rn,j∈ {1。
二维码

扫码加我 拉你入群

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

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

2022-5-19 19:18:03
,d},是差异系数。标准蒙特卡罗方法包括通过用N离散stoch asticdi微分方程来估计E[f(XT)]∈ N*步骤和使用M逼近期望值∈ N**巴黎理工大学,Cermics(ENPC),INRIA,F-77455,Marne la Vall\'ee,France电子邮件:jourdain@cermics.enpc.fr,anis。al公司-gerbi@cermics.enpc.fr-这项研究得益于“主席式金融家”基金会的支持。+巴黎理工大学,LAMA(UMR 8050),UPEMLV,UPEC,CNRS,F-77454,法国马恩拉瓦利,电子邮件:emmanuelle。clement@u-pem。fr.独立路径模拟。为了清楚起见,crud e蒙特卡罗估计量由^YCMC=MMXk=1f给出XN,kT其中,XN,kar是时间步长为T/N的数值格式XN的独立副本。在一些关于SDE系数的正则性假设和平滑支付的情况下,众所周知,为了确保均方根误差,该方法的计算成本为-(2+α), 当eα是数值格式的弱收敛阶时(见[3]中的定理1)。在[8]中,Ninomiya和Victoir提出了一种数值格式,实现了α=2,与α=1的Euler格式相比,该格式降低了计算复杂度。在时间复杂性方面-(2+α), 项1/α是由于偏差E[f(XT)]-EFXNT公司.为了消除这一项,Giles在[5]中引入了一种多级蒙特卡罗估计器,即介电常数望远镜消除偏差。多层蒙特卡罗估计量的构建如下所示:Y=MMXk=1fX1,0,kT+LXl=1MlMlXk=1FXl、l、kT-FXl码-1,l,kT其中L∈ N*是时间步长为T/2L(Ml)0的最后一级离散化≤L≤L∈(N)*)L+1是各级样本量的向量。此外,对于所有l∈ {1,…,L},两个数值格式Xl,lTand Xl-1、用相同的布朗运动对LTA进行模拟。对于每个隔离级别l∈ {0。
二维码

扫码加我 拉你入群

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

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

2022-5-19 19:18:05
,L},Mlindependent和同分布路径模拟使用独立于其他级别的路径模拟。该方法的最优复杂度由方差V的收敛阶β决定FXl,lT- FXl码-1,lT, 这与格式的强收敛阶γ有关。对于Lip-schitz-Payoff f,利用方案在方差估计中的强收敛性,可以得到β≥ 2γ。对于β>1,最佳复杂度为O-2.. 这种复杂性与具有独立且同分布无偏随机变量的简单MonteCarlo方法中的复杂性相同。Milstein方案满足条件β>1,其中γ=1。不幸的是,要模拟Milstein格式,通常需要模拟当布朗运动d的维数大于2时,没有已知方法的L’evy区域。除非扩散系数σj,j∈ {1,…,d},是常数,Euler格式的强阶为γ=1/2,这导致β=1,并导致最优复杂性为-2.日志.最近,已经开发了两种方法来改进γ=1/2的情况。在[6]中,Giles和Szpruch引入了一个改进的Milstein方案,将L'evy区域设置为零,其对偶版本基于方案中每一对连续的布朗增量的交换。关于多层蒙特卡罗估计,在每个离散化水平l∈ {1,…,L},在最网格上,Giles和Szpruch没有使用简单的方案,而是使用了方案的算术平均值及其对偶版本,如下所示:Y=MMXk=1fX1,0,kT+LXl=1MlMlXk=1FXl、l、kT+ FXl、l、kT- FXl码-1,l,kT其中xl表示Giles-Szpruch方案的对偶版本,时间步长为T/2l。
二维码

扫码加我 拉你入群

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

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

2022-5-19 19:18:08
Gilesand Szpruch在一些关于SDE系数和平滑支付的正则性假设下表明,Sγ等于1/2,β等于2,这导致了-2.. 第二种方法称为多级Richardson-Romber g方法,由Lemaire和Pag\'es在[7]中研究,它充分利用了弱误差展开的存在,同时保持了多级Monte Carlo估计的性质。多级Rich ardson-Romber估计量是多级蒙特卡罗方法的加权版本,该方法集成了Pag\'es在[10]中开发的多级Richardson-Romber g外推。Lemaire和Pag\'es得到了一个最优复杂性O-2日志当β=1时,改进了标准多级蒙特卡罗方法。当β>1时,最优复杂度为-2.ISP保留。在这篇论文中,我们建议使用Ninomiya Victoir方案,该方案已知在多级蒙特卡罗估值器最后一级L的最细网格上表现出2阶弱收敛。这一想法受到Debrabant和R¨ossler[2]的启发,他们建议在多层MonteCarlo方法的最底层L的最底层网格上使用高阶弱收敛的方案。通过这种方式,Debrabant和R¨ossler通过减少离散化级别的数量来降低计算复杂性中的常数。在第二节中,为了保证Ninomiya-Victoir格式的强收敛阶,我们提出了在时间网格点之间进行适当插值的思想。然后,在一些关于SDE系数的正则性假设下,我们证明了阶γ=1/2的stron g收敛性。在第3节中,我们提出了一个改进的Ninomiya Victoir方案,该方案可能与Giles Szpruch方案的1阶强耦合。
二维码

扫码加我 拉你入群

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

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

2022-5-19 19:18:11
这一结果允许我们推导出NinomiyaVictoir方案的对立版本,并将Giles Szpru-ch和Debrabant-R¨ossler的思想结合起来,通过将Giles Szpruch方案从0级构建到L级的多层蒙特卡罗估计量- 1以及Ninomiya Victoir方案和Giles Szpruch方案在最后一级L之间的耦合。第4节证实了该估计器的效率,我们在其中详细介绍和评论,如【6】所示,在Clark Cameron SDE和HestonSDE上进行的数值实验。2 Ninomiya-Victoir模式的强收敛性。本节首先介绍本文将使用的一些符号。为了离散化(1.1),我们考虑一个时间步长为h=T/N的非if-orm网格,其中N∈ N*和wedenote:o(tk)k∈[[0;N]]=kh具有相等时间步长h的[0,T]细分,oτs s之前的最后一次离散化∈ [0,T],即^τs=tkif s∈ (tk,tk+1),d表示s=t=0,我们设置^τ=t=0,oeτs在s之后的第一次离散化∈ [0,T],即ˇτs=tk+1if s∈ (tk,tk+1),对于t=0,我们设置ˇτ=t=0,oJ∈ {1,…,d},s∈ [0,T],使得tk<s≤ tk+1,Wjs=Wjs- Wjtk,os∈ [0,T],使得tk<s≤ tk+1,s=s- tk,oη=(η,…,ηN)一系列独立的、同分布的Rademacher随机变量,与W无关,o通过稍微滥用N旋转,我们设置ηs=ηk+1if s∈ (传统知识,传统知识+1),o十、∈ R+,十、 表示唯一的n∈ N*满足n- 1<x≤ n、 o十、∈ R+,十、 表示唯一的n∈ N*满足n≤ x<n+1。让V:Rn-→ RnLipschitz连续,并考虑Rn中的普通微分方程:dx(t)dt=V(x(t))x(0)=x.(2.1)时间t时(2.1)的解,t∈ R表示为x(t)=exp(tV)x,(2.1)的积分形式为x(t)=exp(tV)x=x+ZtV(x(s))ds=x+ZtV(exp(sV)x)ds。我们回忆起(1.1)中的th,每个坐标i∈ {1。
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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