×

线性高斯反问题后验协方差矩阵的低秩特征向量压缩。 (英语) Zbl 1398.65057号

摘要:我们考虑线性高斯贝叶斯反问题后验概率密度协方差矩阵的有效计算问题。当噪声和先验概率密度为高斯时,这种统计逆问题的解也是高斯的。因此,基本解的特征是后验概率密度的均值和协方差矩阵。然而,后验概率密度的协方差矩阵是稠密且大的。因此,对于大维参数空间,无法像离散化的偏微分方程那样计算这样的矩阵。最近引入了后验协方差矩阵的低秩近似作为有前途的工具。然而,对于瞬态问题,得到的近似值会受到维数增加的影响。我们在这里利用离散化方程的结构,使得空间和时间分量可以分离,并且显著降低了不断增长的复杂性。特别是,到目前为止,特征向量低阶近似的存储主要由计算和存储复杂性决定,其中,(n_x)是空间域的维数,(n_t)是时间域的维数。我们开发了一种新的方法,将低阶时间算法与低阶Hessian方法结合使用。我们将计算复杂性和存储需求从\(mathcal O(n_xn_t)\)降低到\(mathcal O(n _x+n _t))。我们使用数值实验来说明我们的方法的优点。

MSC公司:

2015财年65 矩阵特征值和特征向量的数值计算
65层10 线性系统的迭代数值方法
65层50 稀疏矩阵的计算方法
93C20美元 偏微分方程控制/观测系统
2015年1月62日 贝叶斯推断
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] R.Andreev和C.Tobler,抛物偏微分方程时空同时离散的多层预处理和低秩张量迭代,数字。线性代数应用。,22(2015),第317–337页·兹比尔1363.65156
[2] P.Benner、M.Koéhler和J.Saak,M.E.S.S.-矩阵方程稀疏解算器库, .
[3] P.Benner和P.Ku¨rschner,用因子ADI方法计算Sylvester方程的实低阶解,计算。数学。申请。,67(2014),第1656–1672页·兹比尔1364.65093
[4] P.Benner和J.Saak,大型稀疏连续时间代数矩阵Riccati和Lyapunov方程的数值解:最新综述、GAMM-Mitt.、。,36(2013),第32–52页·Zbl 1279.65044号
[5] M.Benzi、E.Haber和L.Taralli,一类PDE约束优化问题的预处理技术,高级计算。数学。,35(2011年),第149-173页·Zbl 1293.65041号
[6] A.Borzí,抛物型分布最优控制问题的多重网格方法,J.计算。申请。数学。,157(2003),第365-382页·Zbl 1024.65047号
[7] T.Bui-Thanh、O.Ghattas、J.Martin和G.Stadler,无限维贝叶斯反问题的计算框架我:线性化情况及其在全球地震反演中的应用,SIAM J.科学。计算。,35(2013年),第A2494–A2523页·Zbl 1287.35087号
[8] D.Calvetti和E.Somersalo,贝叶斯科学计算导论,Surv公司。导师。申请。数学。科学。,2,Springer-Verlag,纽约,2007年·Zbl 1137.65010号
[9] J·钟和M·钟,一种计算最优低秩正则逆矩阵的有效方法《反问题》,30(2014),114009·Zbl 1305.65130号
[10] J·钟和M·钟,逆问题的最优正则化逆矩阵,SIAM J.矩阵分析。申请。,38(2017),第458–477页·Zbl 1367.65059号
[11] S.V.多尔戈夫,TT-GMRES:结构张量格式下线性系统的解,俄罗斯J.Numer。分析。数学。《建模》,28(2013),第149-172页·兹比尔1266.65050
[12] S.V.Dolgov和D.V.Savostyanov,高维线性系统的交替最小能量法,SIAM J.科学。计算。,36(2014年),第A2248–A2271页·Zbl 1307.65035号
[13] H.C.Elman、D.J.Silvester和A.J.Wathen,有限元和快速迭代解及其在不可压缩流体动力学中的应用,数字。数学。科学。计算。,牛津大学出版社,纽约,2005年·Zbl 1083.76001号
[14] H.P.Flath、L.C.Wilcox、V.Akçelik、J.Hill、B.van Bloemen Waanders和O.Ghattas,基于低阶部分Hessian逼近的大规模线性反问题Bayes不确定性量化快速算法,SIAM J.科学。计算。,33(2011年),第407-432页·兹比尔1229.65174
[15] G.H.Golub和C.F.van Loan,矩阵计算第三版,约翰霍普金斯大学数学研究生。科学。,约翰·霍普金斯大学出版社,马里兰州巴尔的摩,1996年·Zbl 0865.65009号
[16] L.Grasedyck、D.Kressner和C.Tobler,低阶张量近似技术的文献综述、GAMM-Mitt.、。,36(2013),第53-78页·Zbl 1279.65045号
[17] W.Hackbusch,多重网格方法及应用斯普林格爵士。计算。数学。,4,Springer-Verlag,纽约,1985年·Zbl 0595.65106号
[18] W.Hackbusch,张量空间与数值张量演算Springer Ser.第42卷。计算。数学。,Springer-Verlag,纽约,2012年·Zbl 1244.65061号
[19] R.Herzog和K.Kunisch,PDE约束优化算法、GAMM-Mitt.、。,33(2010年),第163-176页·兹比尔1207.49034
[20] M.Hinze、R.Pinnau、M.Ulbrich和S.Ulbich,基于PDE约束的优化,Springer-Verlag,纽约,2009年·兹比尔1167.49001
[21] K.Ito和K.Kunisch,变分问题的拉格朗日乘子方法及其应用高级设计控制,15,SIAM,费城,2008·兹比尔1156.49002
[22] J.Kaipio和E.Somersalo,统计和计算反问题,申请。数学。科学。,160,Springer-Verlag,纽约,2005年·Zbl 1068.65022号
[23] D.Kressner和C.Tobler,张量积结构线性系统的Krylov子空间方法,SIAM J.矩阵分析。申请。,31(2010),第1688–1714页·Zbl 1208.65044号
[24] D.Kressner和C.Tobler,参数化线性系统的低秩张量Krylov子空间方法,SIAM J.矩阵分析。申请。,32(2011年),第1288–1316页·Zbl 1237.65034号
[25] C.兰索斯,线性微分和积分算子特征值问题的迭代解法《国家标准局研究期刊》,45(1950),第255-282页。
[26] I.V.Oseledets、S.Dolgov、V.Kazeev、D.Savostyanov、O.Lebedeva、P.Zhlobich、T.Mach和L.Song,TT工具箱, .
[27] I.V.Oseledets和S.V.Dolgov,TT格式中线性系统的求解和矩阵反演,SIAM J.科学。计算。,34(2012年),第A2718–A2739页·Zbl 1259.65071号
[28] J.W.Pearson、M.Stoll和A.J.Wathen,含时PDE约束优化问题的正则化预条件,SIAM J.矩阵分析。申请。,33(2012年),第1126–1152页·Zbl 1263.65035号
[29] Y.Saad,稀疏线性系统的迭代方法第二版,SIAM Philadelphia,2003年·Zbl 1031.65046号
[30] Y.Saad,大特征值问题的数值方法,应用经典。数学。66,费城SIAM,2011年·Zbl 1242.65068号
[31] D.Silvester、H.Elman和A.Ramage,不可压缩流和迭代求解器软件(IFISS)版本3.4, (2015).
[32] V.Simoncini,线性矩阵方程的计算方法SIAM Rev.,58(2016),第377–441页·Zbl 1386.65124号
[33] P.Sonneveld和M.B.van Gijzen,IDR(s):求解大型非对称线性方程组的一系列简单快速算法,SIAM J.科学。计算。,31(2008),第1035-1062页·兹比尔1190.65053
[34] A.Spantini、T.Cui、K.Willcox、L.Tenorio和Y.Marzouk,贝叶斯线性反问题的面向目标的最优逼近,SIAM J.科学。计算。,39(2017年),第S167–S196页·Zbl 1373.15027号
[35] A.Spantini、A.Solonen、T.Cui、J.Martin、L.Tenorio和Y.Marzouk,贝叶斯线性反问题的最优低阶逼近,SIAM J.科学。计算。,37(2015),第A2451–A2487页·Zbl 1325.62060号
[36] M.Stoll,一类含时周期PDE约束优化问题的一次全解IMA J.数字。分析。,34(2014),第1554–1577页·Zbl 1303.65050号
[37] M.Stoll和T.Breiten,PDE约束优化的低阶时间方法,SIAM J.科学。计算。,37(2015),第B1–B29页·Zbl 1330.65153号
[38] A.M.Stuart,反问题:贝叶斯观点,实绩数字。,19(2010年),第451-559页·Zbl 1242.65142号
[39] F.Troéltzsch,偏微分方程的最优控制:理论、方法和应用,AMS,普罗维登斯,RI,2010年·Zbl 1195.49001号
[40] P.Wesseling,多重网格方法简介,纯应用。数学。,约翰·威利父子公司,英国奇切斯特,1992年·Zbl 0760.65092号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。