Ninomiya Victoir方案：强收敛、对偶版本和

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

英文标题：
《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下载：
--> Ninomiya-Victoir_scheme:_strong_convergence,_antithetic_version_and_application_.pdf (514.13 KB)

2022-5-12 14:36:36 上传

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：CTO ICT Applications Quantitative Differential

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=WWdd是a-维标准布朗运动，b:Rn-→ Rn是漂移系数，σj:Rn-→ Rn，j∈ {1, . . .

，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, . . .

使用独立于其他级别的、独立于、同分布的路径模拟。该方法的最优复杂度由方差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方案的相反版本。

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阶强耦合。

这个结果允许我们推导出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, . . .

n}根据以下随机微分方程dxit=bi（Xt）dt+dXj=1σij（Xt）dWjt演变。然后，根据差异系数的规则，人们可以用层叠形式书写（1.1）dXt=σ（Xt）dt+dPj=1σj（Xt）o dWjtX=x（2.2），其中σ=b-民主党=1σjσjandσjis是σjd的雅可比矩阵，定义如下σj=σjiki、 k∈[1；n]]=xkσiji、 k∈[1；n]]。现在，我们介绍[8]中介绍的Ninomiya Victoir方案起点：XNV，ηt=x.o对于k∈ {0…，N- 1} ，如果ηk+1=1:XNV，ηtk+1=exphσ经验Wdtk+1σd. . . 经验Wtk+1σ经验hσXNV，ηtk，（2.3），如果ηk+1=-1:XNV，ηtk+1=exphσ经验Wtk+1σ. . . 经验Wdtk+1σd经验hσXNV，ηtk。（2.4）当我们使用Ninomiya Victoir方案时，首选Stratonovich形式，因为Stratonovich漂移出现在方案的定义中。此外，使用It^o的公式，我们可以t、 s∈ R+，s≤ t、 经验Wjt- Wjs五、y=y+ZtsV经验Wju- Wjs五、Yo dWju。（2.5）然后，重写（2.3）和d（2.4），得到一个XNV，ηtk+1=XNV，ηtk+dXj=1Ztk+1tkσj\'Xj，ηso dWjs+Ztk+1tkσ\'X0，ηs+ σ\'Xd+1，ηsds，（2.6）其中，对于s∈ （tk，tk+1]，\'X0，ηs=expsσXNV，ηtk{ηk+1=1}+\'X1，ηtk+1{ηk+1=-1}, （2.7）对于s∈ （tk，tk+1），j∈ {1，…，d}，\'-Xj，ηs=expWjsσj\'\'Xj-1，ηtk+1{ηk+1=1}+\'Xj+1，ηtk+1{ηk+1=-1}, （2.8）对于s∈ （tk，tk+1]，\'Xd+1，ηs=expsσ\'Xd，ηtk+1{ηk+1=1}+XNV，ηtk{ηk+1=-1}. （2.9）表示“X”-1，ηtk+1=\'Xd+2，ηtk+1=XNV，ηtk，得到一个类似于（2.8）的表达式∈ {0，d+1}和s∈ （tk，tk+1）`Xj，ηs=expsσ\'\'Xj-1，ηtk+1{ηk+1=1}+\'Xj+1，ηtk+1{ηk+1=-1}. （2.10）然后，可以观察到Ninomiya Victoir方案是通过将exactsolution X替换为SDE（1.1）的Stratonovich配方（2.2）中的一个中间过程“Xj，η”而获得的。备注2.1随机过程\'Xj，ηtT∈[0，T]，j∈ {1，…d+1}不适用于自然过滤Ft=σ（Ws，s≤ t） 布朗运动。

为了解决这个问题，我们采用以下过滤方法Fjt=σWjs，s≤ TWk6=jσ工作，工作≤ T, J∈ {1，…，d}。那么，对于j∈ {1，…，d}，根据独立性，Wjis是一个FjBrownian运动，而Xj，η适用于过滤Fj。这确保了每个随机积分都得到了很好的定义。为了研究强收敛性，我们必须建立一个插值格式。允许XNV，ηtT∈[0，T]遵循It^o流程dXNV，ηt=dPj=1σj（\'Xj，ηt）o dWjt+σ\'X0，ηt+ σ\'Xd+1，ηtdtXNV，η=x.（2.11）使用（2.6）和正向感应，可以证明XNV，ηt0≤T≤这是theNinomiya Victoir方案的插值XNV，ηtkK∈[0；N]]。的It^o分解XNV，ηtT∈[0，T]由dXNV，ηt=dPj=1σj（\'Xj，ηt）dWjt+dPj=1σjσj\'Xj，ηtdt+σ\'X0，ηt+ σ\'Xd+1，ηtdtXNV，η=x.（2.12）备注2.2本方案的自然和自适应插值可以是xnv，ηt=hηtTWt，TXNV，η^τt（2.13）其中-1（t，…，td+1；x）=exptσ经验tσ. . . 经验tdσd经验td+1σx（2.14）和h（t，…，td+1；x）=exptσ经验tdσd. . . 经验tσ经验td+1σx、 （2.15）在这两种情况下Wt=Wt，Wdt. 为了得到XNVη的It^o分解，我们必须应用It^o公式。为此，我们必须计算hη的导数。在一般情况下，该函数导数的计算相当复杂。这就是为什么我们不会关注这个插值。2.1强收敛我们记得σj∈ C（注册护士，注册护士），J∈ {1，…，d}，我们假设向量场σj，J∈{0，…，d}，和σjσj，J∈ {1，…，d}是Lipschitz连续函数。显然，b也是solipschitz连续的，因为b=σ+dPj=1σjσj.设L∈ R*+表示它们的共同Lipschitz常数：J∈ {0，…，d}，x、 y∈ Rn，σj（x）- σj（y）≤ LKX- yk，J∈ {1，…，d}，x、 y∈ Rn，σjσj（x）- σjσj（y）≤ LKX- 其中欧几里德范数由k.k.定理2.3表示∈ [1, +∞).

在之前的Lipschitz假设下，存在一个确定性常数CNV∈ R*+以至于N∈ N*, E“supt≤TXt- XNV，ηt2pη#≤ CNV1+kxk2p惠普。当然，这个结果意味着N∈ N*, E“supt≤TXt- XNV，ηt2p#≤ CNV1+kxk2p惠普。明显地XNV，ηt0≤T≤和h依赖于N，但为了保持符号的简单性，对N的依赖性没有明确表示。下面的命题将被用来证明理论。命题2.4让p≥ 1，Y=（Yt）0≤T≤hbe由d维布朗运动驱动的下列n维SDE的解，直到t=h（dYs=α（Ys）ds+β（Ys）dwsy独立于（Wt）t∈[0，h]使得EhkYk2pi<+∞.假设α和β是Lipschitz连续函数，那么：C∈ R*+, t、 s∈ [0，h]，s≤ t、 （一）Eh1+kZtk2pi≤ Eh1+kZk2piexp（Ch）。（2.16）（ii）EhkZt- Zsk2pi≤ C1+EhkZk2pi（t）- s） p.（2.17）如果β=0，我们有一个更好的结果：EhkZt- Zsk2pi≤ C1+EhkZk2pi（t）- s） 2便士。（2.18）常数Conly取决于kα（0）k，kβ（0）k，T，p以及函数α和β的Lipschitz常数。所有这些结果都是众所周知的（例如参见[11]）。2.2中间结果通过使用前面的命题，可以证明该方案具有一致的负弯矩。引理2.5P≥ 1.C∈ R*+, T∈ [0，T]，N∈ N*, J∈ {0，…，d+1}，E“1+\'Xj，ηt2pη#≤ exp（Cˇτt）1+kxk2p.证据：L et p≥ 1和t∈ [0，T]，然后K∈ {0，…，N- 1} 使tk<t≤ tk+1。

对于j=0，\'X0，ηstk<s≤tk+1是以下ODE的解决方案dZs=σ（Zs）dsZtk=XNV，ηtk{ηk+1=1}+\'X1，ηtk+1{ηk+1=-1}.η和W之间的独立性与（2.16）相结合，确保：E“1+\'X0，ηt2pη#≤ E“1+XNV，ηtk{ηk+1=1}+\'X1，ηtk+1{ηk+1=-1}2pη#exp中国={ηk+1=1}E“1+XNV，ηtk2pη#+1{ηk+1=-1} E“1+\'X1，ηtk+12pη#!经验中国.（2.19）同样，对于1≤ J≤ d:\'Xj，ηstk<s≤tk+1是以下SDE的解决方案：=σjσj（Zs）ds+σj（Zs）dWjsZtk=`Xj-1，ηtk+1。使用相同的参数，可以得到：E“1+\'Xj，ηt2pη#≤{ηk+1=1}E“1+\'\'Xj-1，ηtk+12pη#+1{ηk+1=-1} E“1+\'Xj+1，ηtk+12pη#!exp（Ch）。（2.20）显然，对于j=d+1，有一个类似的结果：E“1+\'Xj，ηt2pη#≤{ηk+1=1}E“1+\'Xd，ηtk+12pη#+1{ηk+1=-1} E“1+XNV，ηtk2pη#!经验中国.（2.21）对于所有向量场，全局Lipschitz常数L都是相同的，因此，这三个不等式涉及相同的常数。在这两种常微分方程中，向量场σ乘以1/2，相当于将方程积分到h/2，只需去除乘法因子1/2。这就是为什么在不等式（2.19）和（2.21）中都有一个因子1/2。自从福尔克∈ {0，…，N}，XNV，ηtk=1{ηk=1}Xd+1，ηtk+1{ηk=-1} \'-X0，ηtk，可以使用正向感应onk和正向感应（分别向后）结合使用j∈ {0，…，d+1}如果ηk+1=1（分别为ηk+1=-1） 获得：E“1+\'Xj，ηt2pη#≤ 实验（Ctk+1）1+kxk2p其中C=（d+1）C。下面的引理是位置2.4和引理2.5的直接应用。引理2.6P≥ 1.C∈ R*+, T∈ [0，T]，N∈ N*, J∈ {1, . . .