全部版块 我的主页
论坛 经济学人 二区 外文文献专区
772 34
2022-05-13
英文标题:
《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提出了一个方案,允许构造一个多级蒙特卡罗估值器,以达到最优复杂度$O\\left(\\epsilon^{-2}\\right)$来获得精度$\\epsilon$。本着同样的精神,我们提出了一个改进的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-13 10:46:35
Ninomiya Victoir格式:强收敛,对偶版本和在多级估计中的应用。阿尔·格比,B·乔丹*在本文中,我们对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会减少离散级别的数量。本文致力于计算Y=E[f(XT)],其中f:Rn-→ R是一个支付函数,Xt是时间T的解∈ R*+, 到形式的多维随机微分方程dXt=b(Xt)dt+dPj=1σj(Xt)dWjt,t∈ [0,T]X=X(1.1)这里,X∈ Rn是初始条件,W=WWd是d-维标准布朗运动,b:Rn-→ Rn是漂移系数,σj:Rn-→ Rn,j∈ {1, . . .
二维码

扫码加我 拉你入群

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

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

2022-5-13 10:46:38
,d}是差异因素。标准的蒙特卡罗方法包括通过用N离散stoch散光微分方程来估计E[f(XT)]∈ N*步骤和使用M逼近期望值∈ N**巴黎东部大学,欧洲核子研究中心(ENPC),因里亚,F-77455,法国马恩拉瓦利,电子邮件:jourdain@cermics.enpc.fr,安妮丝。艾尔-gerbi@cermics.enpc.fr-这项研究得益于“主席式金融家”Risque基金会的支持。+巴黎东部大学,喇嘛(UMR 8050),UPEMLV,UPEC,CNRS,F-77454,法国马恩拉瓦利,电子邮件:emmanuelle。clement@u-佩姆。fr.独立路径模拟。明确地说,积垢蒙特卡罗估计量由^YCMC=MMXk=1f给出XN,kT其中XN,kar是数值格式XN的独立副本,时间为s步T/N。在一些关于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我在哪里∈ N*是时间步长为T/2L(Ml)0的最后一级离散化≤L≤L∈(N)*)L+1是各层级样本量的向量。而且,我∈ {1,…,L},两个数值格式Xl,LTX和Xl-1.用相同的布朗运动模拟了LTA。对于每个隔离级别l∈ {0, . . .
二维码

扫码加我 拉你入群

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

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

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

扫码加我 拉你入群

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

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

2022-5-13 10:46:46
Gilesand Szpruch在一些关于SDE系数的正则性假设下,为了平滑支付,证明了Sγ等于1/2,β等于2,这导致了一个最优复杂性-2.. Lemaire和Pag`es在[7]中研究的第二种方法称为多级Richardson-Romber g方法,它充分利用了弱误差展开的存在,同时保持了多级Monte Carlo估计的性质。多级Rich ardson-Romber估计量是多级蒙特卡罗方法的加权版本,该方法集成了Pag`es在[10]中开发的多级Richardson-Romber g外推。Lemaire和Pag`es得到了一个最优复杂度O-2log当β=1时,改进了标准多级蒙特卡罗方法。当β>1时,最优复杂度为0-2.我保留了。在本文中,我们建议使用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-13 10:46:49
这个结果允许我们推导出NinomiyaVictoir方案的一个相反版本,并将Giles Szpru-ch和Debrabant-R¨ossler的思想结合起来,从0级到L级,通过Giles Szpruch方案建立多层蒙特卡罗估计- 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[0,T]的细分,具有相等的时间步长h,o^τs s之前的最后一次离散化∈ [0,T],即τs=tkif s∈ (tk,tk+1),d表示s=t=0,我们设置^τ=t=0,oGτs在s之后的首次离散化∈ [0,T],即ˇτs=tk+1if s∈ (tk,tk+1]和fors=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*令人满意的- 1<x≤ n、 o十、∈ R+,十、 表示唯一的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)中的每个坐标i∈ {1, . . .
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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