×

用于近似Robin边值问题和反射扩散遍历极限的最简单随机游走。 (英语) Zbl 1525.60086号

摘要:提出并分析了一种简单易行的弱意义数值方法来逼近反射随机微分方程。证明了该方法具有一阶弱收敛性。结合蒙特卡罗技术,它可以用来数值求解具有Robin边界条件的线性抛物型和椭圆型偏微分方程。本文的一个关键结果是使用所提出的方法计算遍历极限,即RSDE在(mathbb{R}^d)域内及其边界上的不变律的期望。这允许在紧凑的支持下高效地从分布中采样。考虑并分析了时间平均和集合平均估计量。考虑了一些扩展,包括二阶弱近似、任意斜向反射的情况,以及一种新的自适应弱格式来求解带有Neumann边界条件的泊松偏微分方程。给出的理论结果得到了几个数值实验的支持。

MSC公司:

60华氏35 随机方程的计算方法(随机分析方面)
65立方米 随机微分和积分方程的数值解
60 H10型 随机常微分方程(随机分析方面)
37甲10 生成、随机和随机差分及微分方程
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] AFSHAR,H.和DOMKE,J.(2015)。反射、折射和哈密顿蒙特卡罗。会议记录28第十届神经信息处理系统国际会议2 3007-3015.
[2] Agmon,S.、Douglis,A.和Nirenberg,L.(1959年)。满足一般边界条件的椭圆型偏微分方程解的边界附近估计。一、。普通纯应用程序。数学。12 623-727. ·Zbl 0093.10401号 ·doi:10.1002/cpa.3160120405
[3] AHN,S.、KORATTIKARA,A.和WELLING,M.(2012)。基于随机梯度Fisher评分的贝叶斯后验抽样程序。这个29th实习生。机器学习会议1591-1598年。
[4] 安德森·D.F.、HIGAM·D.J.、莱特·S.C.和威廉斯·R.J.(2019)。关于约束Langevin方程和(生物)化学反应网络。多尺度模型。模拟。17 1-30. ·Zbl 1481.60146号 ·doi:10.1137/18M1190999
[5] ASMUSSEN,S.(1992)。繁忙交通中的排队模拟。数学。操作。物件。17 84-111. ·Zbl 0754.60101号 ·doi:10.1287/门17.1.84
[6] ASMUSSEN,S.、GLYNN,P.和PITMAN,J.(1995)。一维反射布朗运动模拟中的离散化误差。附录申请。普罗巴伯。5 875-896. ·Zbl 0853.65147号
[7] BAYER,C.、SZEPESSY,A.和TEMPONE,R.(2010)。反射扩散和停止扩散的自适应弱近似。蒙特卡罗方法应用。16 1-67. ·Zbl 1196.65029号 ·doi:10.1515/MCMA.2010.001
[8] BENCHÉRIF-MADANI,A.和PARDOUX,É。(2009). 具有Neumann边界条件的泊松方程的概率公式。斯托克。分析。申请。27 739-746. ·Zbl 1171.35493号 ·网址:10.1080/07362990902976520
[9] BENSOUSSAN,A.和LIONS,J.L.(1987)。脉冲控制与拟变分不等式威利。
[10] BERNAL,F.(2019年)。一般有界扩散的Milstein方法的实现。科学杂志。计算。79 867-890. ·Zbl 1417.60062号 ·doi:10.1007/s10915-018-0884-6
[11] BLEI,D.M.、NG,A.Y.、JORDAN,M.I.和LAFFERTY,J.(2003)。潜在Dirichlet分配。J.马赫。学习。物件。3 993-1022. ·Zbl 1112.68379号
[12] BOSSY,M.、GOBET,E.和TALAY,D.(2004年)。用于有效近似反射扩散的对称欧拉格式。J.应用。普罗巴伯。41 877-889. ·兹比尔1076.65009 ·doi:10.1017/s002190020002060x
[13] BROSSE,N.、DURMUS,A.、MOULINES,E.和PEREYRA,M.(2017年)。使用近端Langevin Monte Carlo从具有紧支撑的对数曲线分布中采样。程序。马赫。学习。物件。65 319-342.
[14] BUBECK,S.、ELDAN,R.和LEHEC,J.(2018)。使用投影朗之万蒙特卡罗从对数曲线分布中采样。离散计算。地理。59 757-783. ·Zbl 1397.65010号 ·doi:10.1007/s00454-018-9992-1
[15] Byrne,S.和Girolma,M.(2013)。嵌入流形上的测地蒙特卡罗。扫描。J.统计。40 825-845. ·Zbl 1349.62186号 ·doi:10.1111/sjos.12036
[16] CATTIAUX,P.、LE?N,J.R.和PRIEUR,C.(2017年)。使用欧拉格式对反射扩散进行不变密度估计。蒙特卡罗方法应用。23 71-88. ·Zbl 1370.60132号 ·doi:10.1515/mcma-2017-0104
[17] CELEUX,G.、EL ANBARI,M.、MARIN,J.-M.和ROBERT,C.P.(2012)。回归中的正则化:在信息不足的情况下比较贝叶斯方法和频率论方法。贝叶斯分析。7 477-502·Zbl 1330.62284号 ·doi:10.1214/12-BA716
[18] COSTANTINI,C.、PACCHIAROTTI,B.和SARTORETTO,F.(1998年)。反射扩散过程泛函的数值逼近。SIAM J.应用。数学。58 73-102. ·Zbl 0913.60031号 ·doi:10.1137/S0036139995291040
[19] DA LIO,F.(2008)。具有Neumann边界条件的抛物型方程解的大时间行为。数学杂志。分析。申请。339 384-398. ·Zbl 1136.35011号 ·doi:10.1016/j.jmaa.2007.06.052
[20] DA PRATO,G.和ZABCZYK,J.(1996年)。无限维系统的遍历性.伦敦数学学会讲座笔记系列229.剑桥大学出版社,剑桥·Zbl 0849.60052号 ·doi:10.1017/CBO9780511662829
[21] 丹格菲尔德,C.E.,KAY,D.和BURRAGE,K.(2012年)。通过反射随机微分方程模拟离子通道动力学。物理学。版本E85 051907.
[22] DAVIS,M.H.A.和NORMAN,A.R.(1990)。具有交易成本的投资组合选择。数学。操作。物件。15 676-713. ·Zbl 0717.90007号 ·doi:10.1287/门15.4.676
[23] DELIGIANNIDIS,G.、BOUCHARD-COTE,A.和DOUCET,A.(2019年)。弹性粒子采样器的指数遍历性。安。统计师。47 1268-1287. ·Zbl 1467.60057号 ·doi:10.1214/18-AOS1714
[24] ELLIOTT,R.J.(2001)。连续时间Kronecker引理与鞅收敛性。斯托克。分析。申请。19 433-437. ·Zbl 0990.60045号 ·doi:10.1081/SAP-100002020
[25] FAUGERAS,O.,CLéMENT,F.,DERICHE,R.,KERIVEN,R.、PAPADOPOULO,T.、ROBERTS,J.、VIéVILLE,T.,DEVERNAY,F.、GOMES,J等人(1999年)。脑电和脑电逆问题:伴随状态方法I:连续情况。研究报告编号RR-3673,INRIA。
[26] 费舍尔,右(1953年)。球体上的色散。程序。R.Soc.伦敦。序列号。A、 数学。物理学。科学。217 295-305. ·Zbl 0051.37105号 ·doi:10.1098/rspa.1953.0064
[27] FRAGA FILHO,C.A.D.(2019年)。关于拉格朗日粒子方法中的边界条件和连续介质力学的物理基础。Contin公司。机械。Thermodyn公司。31 475-489. ·doi:10.1007/s00161-018-0702-2
[28] Freidlin,M.(1985)。泛函积分与偏微分方程.数学研究年鉴109.新泽西州普林斯顿,普林斯顿大学出版社·Zbl 0568.60057号 ·doi:10.1515/9781400881598
[29] 弗里德曼(1961)。任意阶抛物型方程解的渐近性。数学学报。106 1-43. ·Zbl 0102.09001号 ·doi:10.1007/BF02545812
[30] GARRONI,M.G.和MENALDI,J.-L.(1992年)。二阶抛物型积分微分问题的格林函数.皮特曼数学系列研究笔记275.哈洛朗曼科技公司·Zbl 0806.45007号
[31] GIHMAN,I.I.和SKOROHOD,A.V.(1972年)。随机微分方程.Ergebnisse der Mathematik und Ihrer Grenzgebiete,乐队72.纽约斯普林格·Zbl 0242.60003号
[32] GILBARG,D.和TRUDINGER,N.S.(1983年)。二阶椭圆偏微分方程,第2版。德国数学研究所[数学科学基本原理] 224. 柏林施普林格·Zbl 0562.35001号 ·doi:10.1007/978-3-642-61798-0
[33] Girolma,M.和Calderhead,B.(2011年)。黎曼流形朗之万和哈密顿蒙特卡罗方法。J.R.统计社会服务。B.统计方法。73 123-214. ·Zbl 1411.62071号 ·文件编号:10.1111/j.1467-9868.2010.00765.x
[34] GOBET,E.(2001年)。反射扩散弱近似的有效方案。蒙特卡罗方法应用。193-202年7月·兹比尔0986.65002
[35] GOBET,E.(2001年)。模拟区域内扩散的欧拉格式和半空间近似。ESAIM概率。斯达。5 261-297. ·Zbl 0998.60081号 ·doi:10.1051/ps:2001112
[36] HANKS,E.M.、JOHNSON,D.S.和HOOTEN,M.B.(2017年)。约束动物运动的反射随机微分方程模型。《农业杂志》。生物与环境。斯达。22 353-372. ·Zbl 1388.62342号 ·doi:10.1007/s13253-017-0291-8
[37] HERNáNDEZ-LERMA,O.和LASSERRE,J.B.(2003)。马尔可夫链与不变概率.数学进步211.巴塞尔Birkhäuser·Zbl 1036.60003号 ·doi:10.1007/978-3-0348-8024-4
[38] 池田,N.和渡边,S.(1989)。随机微分方程与扩散过程,第2版。北荷兰数学图书馆24.荷兰北部,阿姆斯特丹·Zbl 0684.60040号
[39] KABANOV,Y.和SAFARIAN,M.(2009年)。有交易成本的市场:数学理论.斯普林格金融柏林施普林格·Zbl 1186.91006号 ·数字对象标识代码:10.1007/978-3-540-68121-2
[40] KLEIN,J.P.和MOESCHBERGER,M.L.(2005)。生存分析:截尾数据技术柏林施普林格。
[41] KURITZ,K.、ST-HR,D.、POLLAK,N.和ALLG-WER,F.(2017年)。用遍历原理进行细胞周期分析与年龄结构细胞种群模型之间的关系。J.理论。生物。414 91-102. ·Zbl 1369.92030号 ·doi:10.1016/j.jtbi.2016.11.024
[42] O.A.女士、V.A.SOLONNIKOV和N.N.乌拉尔·塞瓦(1968)。抛物型线性和拟线性方程.数学专著翻译,第卷第二十三条。阿米尔。数学。佛罗里达州普罗维登斯Soc·Zbl 0174.15403号
[43] LEIMKUHLER,B.和MATTHEWS,C.(2015)。分子动力学:用确定性和随机数值方法.跨学科应用数学39.查姆施普林格。
[44] LEIMKUHLER,B.和SHANG,X.(2016)。噪音梯度系统的自适应恒温器。SIAM J.科学。计算。38 A712-A736·Zbl 1382.65019号 ·数字对象标识码:10.1137/15M102318X
[45] LEITE,S.C.和WILLIAMS,R.J.(2019年)。化学反应网络的约束Langevin近似。附录申请。普罗巴伯。29 1541-1608. ·Zbl 1467.60059号 ·doi:10.1214/18-AAP1421
[46] LIEBERMAN,G.M.(1986年)。斜导数问题的中间Schauder估计。架构(architecture)。定额。机械。分析。93 129-134·Zbl 0603.35025号 ·doi:10.1007/BF00279956
[47] 狮子,P.-L.和SZNITMAN,A.-S.(1984)。具有反射边界条件的随机微分方程。普通纯应用程序。数学。37 511-537. ·Zbl 0598.60060号 ·doi:10.1002/cpa.3160370408
[48] 刘毅(1993)。具有边界条件的随机微分方程的数值方法。普渡大学博士论文。
[49] Lunardi,A.(1995)。抛物问题中的解析半群与最优正则性.现代Birkhä用户经典Birkhäuser/Springer Basel AG,巴塞尔·兹伯利0816.35001
[50] Mardia,K.V.和Jupp,P.E.(2000年)。方向统计.概率统计中的威利级数奇切斯特·威利·Zbl 0935.62065号
[51] Mardia,K.V.、Taylor,C.C.和Subramaniam,G.K.(2007年)。角度数据的蛋白质生物信息学和双变量von Mises分布的混合物。生物识别63 505-512. ·兹比尔1136.62082 ·doi:10.1111/j.1541-0420.2006.00682.x
[52] J.C.Mattingly、A.M.Stuart和M.V.Tretyakov(2010年)。通过泊松方程数值时间平均和平稳测度的收敛性。SIAM J.数字。分析。48 552-577. ·Zbl 1217.65014号 ·doi:10.1137/090770527
[53] MIL'SHTEIN,G.N.(1996)。随机方程数值积分在求解具有Neumann边界条件的边值问题中的应用。理论问题。申请。41 170-177. ·Zbl 0888.60050号
[54] MIL’SHTEIN,G.N.和TRET'YAKOV,M.V.(2002)。Dirichlet问题的最简单随机游动。理论问题。申请。47 53-68. ·Zbl 1038.60066号
[55] MILSTEIN,G.N.和TRETYAKOV,M.V.(1999)。时空有界扩散的模拟。附录申请。普罗巴伯。9 732-779. ·兹标0964.60065 ·doi:10.1214/aoap/1029962812
[56] MILSTEIN,G.N.和TRETYAKOV,M.V.(2002年)。非线性抛物方程Neumann问题的概率解法。IMA J.数字。分析。22 599-622. ·Zbl 1018.65006号 ·doi:10.1093/imanum/224.599
[57] MILSTEIN,G.N.和TRETYAKOV,M.V.(2004)。数学物理中的随机数值.科学计算柏林施普林格·Zbl 1085.60004号 ·doi:10.1007/978-3-662-10063-9
[58] MILSTEIN,G.N.和TRETYAKOV,M.V.(2007年)。计算Langevin方程的遍历极限。物理学。D类229 81-95. ·Zbl 1117.65012号 ·doi:10.1016/j.physd.2007.03.011
[59] 米兰达·C(1970)。椭圆型偏微分方程.Ergebnisse der Mathematik und Ihrer Grenzgebiete,乐队2.施普林格,纽约-柏林。
[60] MOLČANOV,S.A.(1964年)。关于扩散过程理论中的一个问题。理论问题。申请。9 472-477. ·Zbl 0161.37101号
[61] NAGASAWA,M.和SATO,K.(1963年)。关于马尔可夫过程的时间变化和消亡的一些定理。Kodai数学。学期代表。15 195-219. ·Zbl 0123.35202号
[62] PAYNE,L.E.和PHILIPPIN,G.A.(1995年)。二阶抛物问题及其导数解的衰减界。数学。模型方法应用。科学。5 95-110. ·Zbl 0832.35019号 ·doi:10.1142/S0218202595000061
[63] 彼得森R.(1997)。反映随机微分方程的惩罚方案。伯努利3 403-414. ·兹比尔0899.60053 ·doi:10.2307/3318456
[64] PINEDA,M.和STAMATAKIS,M.(2018年)。通过反射化学朗之万方程对表面反应进行随机建模。计算。化学。工程师。117 145-158.
[65] Roberts,G.O.和Tweedie,R.L.(1996年)。Langevin分布及其离散近似的指数收敛性。伯努利2 341-363. ·Zbl 0870.60027号
[66] SAISHO,Y.(1987)。具有反射边界的多维区域的随机微分方程。可能性。理论相关领域74 455-477. ·Zbl 0591.60049号 ·doi:10.1007/BF00699100
[67] SATO,K.和UENO,T.(1964/65)。边界上的多维扩散和马尔可夫过程。数学杂志。京都大学。4 529-605. ·Zbl 0219.60057号 ·doi:10.1215/kjm/1250524605
[68] SŁOMIŃSKI,L.(2001)。具有反射边界的SDE解的欧拉近似。随机过程。申请。94 317-337. ·Zbl 1053.60062号 ·doi:10.1016/S0304-4149(01)00087-4
[69] TALAY,D.(1990年)。用于计算不变律的随机微分系统的二阶离散化方案。斯托克。斯托克。代表。29 13-36. ·Zbl 0697.60066号
[70] TANAKA,H.(1979年)。凸区域中具有反射边界条件的随机微分方程。广岛数学。J。9 163-177. ·Zbl 0423.60055号
[71] VOLKONSKI,V.A.(1958)。强马尔可夫过程中时间的随机替换。理论问题。申请。3 310-326. ·Zbl 0102.14004号
[72] WELLING,M.和TEH,Y.W.(2011年)。基于随机梯度Langevin动力学的贝叶斯学习。会议记录28第十届国际机器学习会议681-688.
[73] 温策尔,A.D.(1981)。随机过程理论课程.麦格劳-希尔国际图书公司,纽约·兹比尔0502.60001
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。