×

高斯随机场的多级近似:快速模拟。 (英语) Zbl 1448.60110号

摘要:我们提出并分析了几种用于快速模拟可能非平稳高斯随机场(GRF)的多级算法,例如,通过有界域(mathcal{D}\subset\mathbb{R}^n)的闭包,或者更一般地,通过紧度量空间(mathca{X}),例如紧流形(mathcal{M})。我们算法允许的彩色GRF(mathcal{Z})解一些(beta>n/4)的随机分数阶微分方程(mathcal{A}^ beta\mathcal}Z}=mathcal{W}),其中(mathcali{A})是发散形式的线性、局部、二阶椭圆自共轭微分算子是\(\mathcal{X}\)上的白噪声。因此,我们考虑了具有形式为(mathcal{C}=mathcal{A}^{-2\beta})的协方差算子的\(mathcal{X}\)上的GRF。所提出的算法对嵌套序列上的(mathcal{Z})样本进行了数值逼近{T}_\规则的简单分区(mathcal){T}_\(\mathcal{D}\)和(\mathcal{M}\)的ell\)。计算GRF(mathcal{Z})在三角剖分上的一个近似实现的工作和内存{T}_\具有一致性的\(\mathcal{X}\)的\(\mathcal{O}\left(N_\ell^{-\rho}\right)\),对于某些一致性阶\(\rho>0\),在\(N_\ell=\#(\mathcal{T}_\ell)\),独立于GRF可能的低正则性。这些算法基于负分数阶椭圆“着色”算子(mathcal{a}^{-\beta})到白噪声(mathcal{W})的积分表示(应用)的正弦求积。对于所提出的数值逼近,我们证明了各种范数下计算量和一致性误差的界。

MSC公司:

60G60型 随机字段
35兰特 分数阶偏微分方程
35卢比60 随机偏微分方程的偏微分方程
65立方厘米 马尔可夫链的数值分析或方法
65N30型 含偏微分方程边值问题的有限元、Rayleigh-Ritz和Galerkin方法
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Abramowitz,M.和Stegun,I.A.,《公式、图形和数学表数学函数手册》,第55卷(美国政府印刷局,华盛顿特区,1964年)·Zbl 0171.38503号
[2] Balakrishnan,A.V.,闭算子的分数幂及其生成的半群,Pac。《数学杂志》10(1960)419-437·Zbl 0103.33502号
[3] Balakrishnan,A.V.,应用功能分析,第二版。,第3卷(Springer-Verlag,纽约,1981年)·Zbl 0459.46014号
[4] Beaky,M.M.,Scherrer,R.J.和Villumsen,J.V.,种子热暗物质模型中大尺度结构的拓扑,天体物理学。J.387(1992)443-448。
[5] Bertozzi,A.L.,Luo,X.,Stuart,A.M.和Zygalakis,K.C.,高维数据基于图形分类的不确定性量化,SIAM/ASA J.不确定性。Quantif.6(2018)568-595·Zbl 1394.62083号
[6] Bolin,D.和Kirchner,K.,高斯随机场的有理SPDE方法与一般平滑度,J.Compute。图表。统计数据(2019年),https://doi.org/101080/10618600.2019.1665537。 ·兹比尔07499255
[7] Bolin,D.,Kirchner,K.和Kovács,M.,具有空间白噪声的分数阶椭圆随机偏微分方程的数值解,IMA J.Numer。分析。(2018), https://doi.org/10.1093/imanum/dry091。 ·Zbl 1405.65147号
[8] Bolin,D.,Kirchner,K.和Kovács,M.,具有空间白噪声的分数阶椭圆随机偏微分方程的Galerkin近似的弱收敛性,BIT-Numer。数学58(2018)881-906·Zbl 1405.65147号
[9] Bolin,D.和Lindgren,F.,《嵌套随机偏微分方程生成的空间模型及其在全球臭氧测绘中的应用》,Ann.Appl。统计数字5(2011)523-550·Zbl 1235.60075号
[10] Bonito,A.和Pasciak,J.E.,椭圆算子分数幂的数值逼近,数学。计算84(2015)2083-2110·Zbl 1331.65159号
[11] Bramble,J.H.和Pasciak,J.E.,多层算法的新估计,包括V循环,数学。计算60(1993)447-471·Zbl 0783.65081号
[12] Bramble,J.H.,Pasciak,J.E.和Steinbach,O.,关于\(H^1(\ operatorname{\Omega})\)中\(L^2)投影的稳定性,数学。计算71(2002)147-156·Zbl 0989.65122号
[13] Cameletti,M.、Lindgren,F.、Simpson,D.和Rue,H.,通过SPDE方法对颗粒物浓度进行时空建模,AStA Adv.Stat.Anal.97(2013)109-131·Zbl 1443.62401号
[14] Chan,G.和Wood,A.T.A.,《算法AS 312:模拟平稳高斯随机场的算法》,J.R.Stat.Soc.C,Appl。《美国联邦法律大全》第46卷(1997年)第171-181页·Zbl 0913.65142号
[15] J.Chen和M.L.Stein,高斯随机场的线性成本协方差函数,预印本(2017),arXiv:1711.05895[stat.ME]。
[16] Coifman,R.R.、Kevrekidis,I.G.、Lafon,S.、Maggioni,M.和Nadler,B.,《扩散图、还原坐标和随机系统的低维表示》,多尺度模型。模拟7(2008)842-864·Zbl 1175.60058号
[17] Coifman,R.R.和Maggioni,M.,扩散小波,应用。计算。哈蒙。分析21(2006)53-94·邮编1095.94007
[18] S.G.Cox和K.Kirchner,广义Whittle-Matérn域在Sobolev和Hölder空间中的正则性和收敛性分析,2019-23号技术报告,瑞士苏黎世联邦理工学院数学系应用数学研讨会(2019年)·兹比尔1455.65210
[19] Creasey,P.E.和Lang,A.,球面上各向同性高斯随机场的快速生成,蒙特卡罗方法应用24(2018)1-11·Zbl 1386.60180号
[20] Croci,M.,Giles,M.B.,Rognes,M.E.和Farrell,P.E.,《非嵌套网格多层Monte Carlo的高效白噪声采样和耦合》,SIAM/ASA J.不确定性。Quantif.6(2018)1630-1655·Zbl 07003649号
[21] Dahmen,W.,算子方程的多尺度和小波方法,《数值模拟中的多尺度问题和方法》,第1825卷(Springer,柏林,2003),第31-96页·Zbl 1036.65095号
[22] Dahman,W.,Harbrecht,H.和Schneider,R.,边界积分方程的压缩技术——渐近最优复杂度估计,SIAM J.Numer。分析43(2006)2251-2271·Zbl 1113.65114号
[23] Davies,E.B.,《谱、理论和微分算子》,《剑桥高等数学研究》第42卷(剑桥大学出版社,剑桥,1995年)·Zbl 0893.47004号
[24] Demlow,A.,曲面上椭圆问题的高阶有限元方法和逐点误差估计,SIAM J.Numer。分析47(2009)805-827·Zbl 1195.65168号
[25] Demmel,J.W.,Eisenstat,S.C.,Gilbert,J.R.,Li,X.S.和Liu,J.W.H.,稀疏部分旋转的超节点方法,SIAM J.Matrix Anal。申请20(1999)720-755·Zbl 0931.65022号
[26] Dietrich,C.R.和Newsam,G.N.,通过协方差矩阵的循环嵌入快速准确地模拟平稳高斯过程,SIAM J.Sci。计算18(1997)1088-1107·Zbl 0890.65149号
[27] Dodwell,T.J.、Ketelsen,C.、Scheichl,R.和Teckentrup,A.L.,一种多层马尔可夫链蒙特卡罗算法及其在地下水流不确定性量化中的应用,SIAM/ASA J.不确定性。Quantif.3(2015)1075-1108·Zbl 1330.65007号
[28] Dölz,J.,Harbrecht,H.和Schwab,Ch.,粗糙随机场的协方差正则性和(mathcal{H})-矩阵近似,Numer。数学135(2017)1045-1071·Zbl 1362.65124号
[29] Feischl,M.,Kuo,F.Y.和Sloan,I.H.,用H矩阵快速生成随机场,数值。数学140(2018)639-676·Zbl 1404.65033号
[30] Gantner,R.N.,Herrmann,L.和Schwab,Ch.,仿射参数椭圆偏微分方程的拟蒙特卡罗积分:局部支持和乘积权重,SIAM J.Numer。分析56(2018)111-135·Zbl 1378.65072号
[31] D.Gilbarg和N.S.Trudinger,二阶椭圆型偏微分方程,1998年再版。,数学经典(Springer-Verlag,柏林,2001)·Zbl 1042.35002号
[32] Golub,G.H.和Van Loan,C.F.,《矩阵计算》,第四版。,(约翰·霍普金斯大学出版社,巴尔的摩,2013年)·Zbl 1268.65037号
[33] Graham,I.G.,Kuo,F.Y.,Nuyens,D.,Scheichl,R.和Sloan,I.H.,用QMC嵌入循环:对数正态系数椭圆偏微分方程的分析,Numer。数学140(2018)479-511·Zbl 1417.65203号
[34] Grisvard,P.,《非光滑域中的椭圆问题》,第69卷(工业和应用数学学会(SIAM),费城,2011年)·兹比尔1231.35002
[35] Hale,N.,Higham,N.J.和Trefethen,L.N.,《计算》(\mathbf{A}^ alpha,\log(\mathbf{A{)),以及通过轮廓积分实现的相关矩阵函数,SIAM J.Numer。分析46(2008)2505-2523·Zbl 1176.65053号
[36] Hammond,D.K.,Vandergheynst,P.和Gribonval,R.,通过谱图理论研究图上的小波,应用。计算。哈蒙。分析30(2011)129-150·Zbl 1213.42091号
[37] Herrmann,L.,Lang,A.和Schwab,Ch.,球面上对数正态扩散的数值分析,Stoch。部分差异。Equ.、。,分析。计算6(2018)1-44·兹比尔1402.60062
[38] Herrmann,L.和Schwab,Ch.,对数正态系数椭圆偏微分方程的带乘积权重的多层拟蒙特卡罗积分,ESAIM Math。模型。数字。分析53(2019)1507-1552·Zbl 07135561号
[39] L.Herrmann和Ch.Schwab,对流-扩散反应的多级准蒙特卡罗不确定度量化,第2019-06号技术报告,ETH Zürich数学系应用数学研讨会(2019年)[发表在第13届国际蒙特卡罗与准蒙特卡洛方法科学计算会议(2018年)上]·Zbl 07240088号
[40] Herrmann,L.和Schwab,Ch.,对数正态参数椭圆偏微分方程的QMC集成:局部支持和产品权重,数值。数学141(2019)63-102·Zbl 07006664号
[41] R.Hiptmair和L.Kielhorn,《BETL:通用边界元模板库》,第2012-36号技术报告,瑞士苏黎世联邦理工学院数学系应用数学研讨会(2012年)。
[42] Kolev,T.V.,Xu,J.和Zhu,Y.,不连续系数反应扩散问题的多层预条件,科学杂志。计算67(2016)324-350·Zbl 1353.65131号
[43] Lang,A.和Schwab,Ch.,《球面上的各向同性高斯随机场:正则性、快速模拟和随机偏微分方程》,Ann.Appl。Probab.25(2015)3047-3094·Zbl 1328.60126号
[44] Latz,J.,Eisenberger,M.和Ullmann,E.,参数化高斯随机场的快速采样,计算机。方法应用。机械。工程348(2019)978-1012·Zbl 1441.65001号
[45] M.Ledoux和M.Talagrand,《巴拿赫空间中的概率:等高线和过程》,1991年再版。,数学经典(Springer-Verlag,柏林,2011)·Zbl 1226.60003号
[46] Lindgren,F.,Rue,H.和Lindström,J.,高斯场和高斯-马尔可夫随机场之间的显式联系:随机偏微分方程方法,J.R.Stat.Soc.B,Stat.Methodol.73(2011)423-498·Zbl 1274.62360号
[47] J.-L.Lions和E.Magenes,非均匀边值问题和应用,第一卷,由P.Kenneth、Grundlehren der Mathematischen Wissenschaften译自法语,第181卷(Springer-Verlag,纽约,1972)·Zbl 0223.35039号
[48] Lunardi,A.,《插值理论》,第16卷(Edizioni della Normale,比萨,2018)·Zbl 1394.46001号
[49] Lund,J.和Bowers,K.L.,《求积和微分方程的Sinc方法》(工业和应用数学学会(SIAM),费城,1992年)·Zbl 0753.65081号
[50] Marinucci,D.和Peccati,G.,《球面上的随机场:表示、极限定理和宇宙学应用》,第389卷(剑桥大学出版社,剑桥,2011年)·Zbl 1260.60004号
[51] Y.Meyer,《小波与算子》,D.H.Salinger译自1990年法文原著,《剑桥高等数学研究》,第37卷(剑桥大学出版社,剑桥,1992年)·Zbl 0776.42019号
[52] Nguyen,H.和Stevenson,R.,《具有改进定量特性的有限元小波》,J.Compute。申请。数学230(2009)706-727·兹比尔1173.65077
[53] Osborn,S.、Vassilevski,P.S.和Villa,U.,《空间相关随机场的多级分层抽样技术》,SIAM J.Sci。计算39(2017)S543-S662·Zbl 1422.65406号
[54] Osborn,S.、Zulian,P.、Benson,T.、Villa,U.、Krause,R.和Vassilevski,P.S.,使用非匹配网格生成空间相关随机场的可缩放分层PDE采样器Numer。线性代数应用25(2018)e2146·Zbl 1499.65006号
[55] Penny,W.D.、Trujillo-Barreto,N.J.和Friston,K.J.,贝叶斯功能磁共振成像时间序列分析与空间先验,神经影像24(2005)350-362。
[56] Rekatsinas,N.和Stevenson,R.,二次有限元小波Riesz基,Int.J.Wavelets Multiresolut。信息处理16(2018)1850033·Zbl 1507.65306号
[57] Sain,S.R.,Furrer,R.和Cressie,N.,《区域气候模型多元输出的空间分析》,Ann.Appl。统计数据5(2011)150-175·Zbl 1220.62152号
[58] Scheichl,R.,Stuart,A.M.和Teckentrup,A.L.,椭圆反问题中计算后验期望的拟蒙特卡罗和多级蒙特卡罗方法,SIAM/ASA J.不确定性。Quantif.5(2017)493-518·兹比尔1516.65118
[59] O.Steinbach,《椭圆边值问题的数值逼近方法:有限元和边界元》,译自2003年德国原文(Springer,纽约,2008)·Zbl 1130.65001号
[60] Stenger,F.,《基于Sinc和解析函数的数值方法》,第20卷(Springer-Verlag,纽约,1993年)·Zbl 0803.65141号
[61] Stevenson,R.,《求解算子方程的自适应小波方法:概述》,载于《多尺度、非线性和自适应逼近》(Springer,Berlin,2009),第543-597页·Zbl 1192.65063号
[62] Taylor,M.E.,《伪微分算子》,第34卷(普林斯顿大学出版社,普林斯顿,1981年)·Zbl 0453.47026号
[63] Torquato,S.,《随机异质材料:微观结构和宏观性能》,第16卷(Springer-Verlag,纽约,2002年)·Zbl 0988.74001号
[64] Whittle,P.,《多维随机过程》,公牛。Inst.Int.Stat.40(1963)974-994·Zbl 0129.10603号
[65] Xu,J.,通过空间分解和子空间校正的迭代方法,SIAM Rev.34(1992)581-613·Zbl 0788.65037号
[66] Xu,J.和Zikatanov,L.,希尔伯特空间中的交替投影方法和子空间校正方法,J.Amer。数学。Soc.15(2002)573-597·Zbl 0999.47015号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。