×

线性离散不适定问题的Arnoldi分解、GMRES和预处理。 (英语) Zbl 1417.65119号

概述:GMRES是求解大型线性方程组的最常用迭代方法之一,这些方程组是由线性适定问题的离散化引起的,例如椭圆偏微分方程的边值问题。该方法也适用于线性方程组的迭代解,这些方程组是通过离散线性不适定问题(例如许多反问题)而获得的。然而,当应用于后一类问题时,GMRES并不总是表现良好。本文试图阐明GMRES在某些情况下性能不佳的原因,并讨论基于特定类型预处理的一些补救措施。GMRES的标准实现基于Arnoldi过程,该过程也可用于定义Tikhonov或TSVD正则化的解子空间,分别产生Arnoldi Tikhonov和Arnoldi TSVD方法。讨论了GMRES、Arnoldi-Tikhonov和Arnoldi-TSVD方法的性能。数值例子说明了这些方法的特性。

MSC公司:

65层22 数值线性代数中的不适定性和正则化问题
65层10 线性系统的迭代数值方法
65层35 矩阵范数、条件、缩放的数值计算
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Baart,M.L.,《在噪声条件下的最小二乘问题中使用自相关确定伪库》,IMA J.Numer。分析。,2, 241-247 (1982) ·Zbl 0484.65021号
[2] 澳大利亚比约克。,矩阵计算中的数值方法(2015),Springer:Springer New York·Zbl 1322.65047号
[3] Brown,P.N。;Walker,H.F.,(几乎)奇异系统的GMRES,SIAM J.矩阵分析。申请。,18, 37-51 (1997) ·Zbl 0876.65019号
[4] Calvetti,D。;刘易斯,B。;Reichel,L.,《利用GMRES方法恢复具有空间变化模糊的图像》,(Luk,F.T.,《光电仪器工程师学会(SPIE)学报》,X.《光电仪器工程学会(SPEE)学报》,第4116卷(2000年),国际光学工程学会:国际光学工程协会(华盛顿州贝灵汉),364-374
[5] Calvetti,D。;刘易斯,B。;Reichel,L.,《关于线性离散不定问题迭代方法的子空间选择》,Int.J.Appl。数学。计算。科学。,11, 1069-1092 (2001) ·Zbl 0994.65043号
[6] Calvetti,D。;刘易斯,B。;Reichel,L.,Krylov子空间迭代方法,用于图像恢复中的非对称离散不适定问题,(Luk,F。高级信号处理算法、架构和实现XI。高级信号处理算法、架构和实现XI,《光电仪器工程师学会学报》,第4474卷(2001年),国际光学工程学会:国际光学工程协会,华盛顿州贝灵汉,224-233
[7] Calvetti,D。;刘易斯,B。;Reichel,L.,关于GMRES方法的正则化性质,Numer。数学。,91, 605-625 (2002) ·兹比尔1022.65044
[8] Calvetti,D。;刘易斯,B。;Reichel GMRES,L.,L曲线和离散不定问题,BIT,42,44-65(2002)·Zbl 1003.65040号
[9] Calvetti,D。;Morigi,S。;赖切尔,L。;Sgallari,F.,Tikhonov正则化和大型离散不定问题的L曲线,J.Compute。申请。数学。,123, 423-446 (2000) ·Zbl 0977.65030号
[10] Chan,R.H.F。;Jin,X.-Q.,迭代Toeplitz解算器简介(2007),SIAM:SIAM Philadelphia·Zbl 1146.65028号
[11] Chan,T.F.,Toeplitz系统的最优循环预条件,SIAM J.Sci。统计计算。,9, 766-771 (1988) ·Zbl 0646.65042号
[12] Chan,T.F。;Jackson,K.R.,离散牛顿算法的非线性预处理Krylov子空间方法,SIAM J.Sci。统计计算。,5, 533-542 (1984) ·Zbl 0574.65043号
[13] Davis,P.J.,循环矩阵(1994),切尔西:切尔西纽约·Zbl 0898.15021号
[14] 多纳泰利,M。;马丁·D·。;Reichel,L.,具有抗反射边界条件的图像去模糊的Arnoldi方法,Appl。数学。计算。,253, 135-150 (2015) ·Zbl 1338.94006号
[15] 杜,K。;Duintjer Tebbens,J。;Meurant,G.,对于GMRES,Electron,任何允许的谐波Ritz值集都是可能的。事务处理。数字。分析。,47, 37-56 (2017) ·Zbl 1372.65107号
[16] 戴克斯,L。;Noschese,S.公司。;Reichel,L.,离散不定Toeplitz系统的循环预条件,数值。算法,75,477-490(2017)·Zbl 1368.65039号
[17] 戴克斯,L。;Reichel,L.,线性离散不适定问题的一系列范围受限迭代方法,(Dolomites Research Notes on Approximation,vol.6(2013)),27-36·Zbl 1291.65117号
[18] Eisenstat,S.C.,偏对称线性系统的Krylov子空间方法的等效性(2015)
[19] 英国,H.W。;汉克,M。;Neubauer,A.,反问题的正则化(1996),Kluwer:Kluwer-Dordrecht·Zbl 0859.65054号
[20] Gazzola,S。;Hansen,P.C。;Nagy,J.,IR工具:迭代正则化方法和大规模测试问题的MATLAB包,Numer。算法(2018)
[21] Gazzola,S。;诺瓦蒂,P。;Russo,M.R.,Tikhonov正则化中选择参数的嵌入式技术,Numer。线性代数应用。,21, 796-812 (2014) ·Zbl 1340.65068号
[22] Gazzola,S。;诺瓦蒂,P。;Russo,M.R.,《关于Krylov投影方法和Tikhonov正则化》,电子。事务处理。数字。分析。,44, 83-123 (2015) ·Zbl 1312.65065号
[23] 格林鲍姆,A。;Ptak,V。;Strakoš,Z.,对于GMRES,任何非增量收敛曲线都是可能的,SIAM J.Matrix Anal。申请。,17, 465-469 (1996) ·兹比尔0857.65029
[24] Greengard,L。;Rokhlin,V.,三维拉普拉斯方程快速多极子方法的新版本,《数值学报》。,6, 229-269 (1997) ·Zbl 0889.65115号
[25] Hanke,M.,《病态问题的共轭梯度型方法》(1995),Longman:Longman-Harlow·Zbl 0830.65043号
[26] 汉克,M。;Nagy,J。;Plemmons,R.,不适定问题的预条件迭代正则化,(Reichel,L.;Ruttan,A.;Varga,R.S.,《数值线性代数》(1993),德格鲁特:德格鲁伊特-柏林),141-163·Zbl 0794.65039号
[27] Hansen,P.C.,秩亏和离散病态问题(1998),SIAM:SIAM Philadelphia
[28] Hansen,P.C.,针对MATLAB 7.3的正则化工具版本4.0,数值。算法,46,189-194(2007)·Zbl 1128.65029号
[29] Hansen,P.C。;Jensen,T.K.,正则化最小残差方法的平滑范数预处理,SIAM J.矩阵分析。,29, 1-14 (2006) ·Zbl 1154.65028号
[30] Hansen,P.C。;Jensen,T.K.,图像去模糊正则化迭代中的噪声传播,电子。事务处理。数字。分析。,31, 204-220 (2008) ·Zbl 1171.65032号
[31] 海厄姆,N.J.,计算最近对称半正定矩阵,线性代数应用。,103, 103-118 (1988) ·Zbl 0649.65026号
[32] Huckle,T.,正规矩阵的Arnoldi方法,SIAM J.矩阵分析。申请。,15, 479-489 (1994) ·兹伯利0796.65037
[33] Jensen,T.K。;Hansen,P.C.,最小残差法迭代正则化,BIT,47,103-120(2007)·Zbl 1113.65037号
[34] Kindermann,S.,线性不适定问题基于最小化的无噪声级参数选择规则的收敛性分析,Electron。事务处理。数字。分析。,38, 233-257 (2011) ·Zbl 1287.65043号
[35] 刘易斯,B。;Reichel,L.,Arnoldi-Tikhonov正则化方法,J.Compute。申请。数学。,226, 92-102 (2009) ·Zbl 1166.65016号
[36] 洛根,D。;鲁伊斯,D。;Touhami,A.,非线性方程组的自适应预条件,J.Compute。申请。数学。,189, 362-374 (2006) ·Zbl 1088.65048号
[37] Meurant,G.,《大型线性系统的计算机解决方案》(1999),Elsevier:Elsevier Amsterdam·Zbl 0934.65032号
[38] Moret,I.,关于GMRES超线性收敛性的注记,SIAM J.Numer。分析。,34, 513-516 (1997) ·Zbl 0873.65054号
[39] 新墨西哥州纳希蒂加尔。;南卡罗来纳州雷迪。;Trefethen,L.N.,非对称矩阵迭代有多快?,SIAM J.矩阵分析。申请。,13, 778-795 (1992) ·Zbl 0754.65036号
[40] 新墨西哥州纳希蒂加尔。;赖切尔,L。;Trefethen,L.N.,非对称线性系统的混合GMRES算法,SIAM J.矩阵分析。申请。,13, 796-825 (1992) ·Zbl 0757.65035号
[41] Neuman,A。;赖切尔,L。;Sadok,H.,线性离散不适定问题的范围受限迭代方法的实现,线性代数应用。,436, 3974-3990 (2012) ·Zbl 1241.65045号
[42] Ng,M.K.,《Toeplitz系统的迭代方法》(2004),牛津大学出版社:牛津大学出版社·Zbl 1059.65031号
[43] Noschese,S.公司。;帕斯奎尼,L。;Reichel,L.,不可约实三对角矩阵到正规的结构化距离,电子。事务处理。数字。分析。,28, 65-77 (2007) ·Zbl 1171.65037号
[44] Noschese,S.公司。;帕斯奎尼,L。;Reichel,L.,《三对角Toeplitz矩阵:性质和新应用》,数值。线性代数应用。,20, 302-326 (2013) ·兹比尔1289.65082
[45] Noschese,S.公司。;Reichel,L.,Toeplitz矩阵的正规结构距离及其在预处理中的应用,Numer。线性代数应用。,18, 429-447 (2011) ·Zbl 1245.65036号
[46] Noschese,S.公司。;Reichel,L.,离散不定问题的修正TSVD方法,数值。线性代数应用。,21813-822(2014)·Zbl 1340.65070号
[47] Noschese,S.公司。;Reichel,L.,关于超优广义循环预条件的注记,应用。数字。数学。,75, 188-195 (2014) ·Zbl 1302.65072号
[48] Novati,P.,基于Arnoldi的线性不适定问题方法的一些性质,SIAM J.Numer。分析。,55, 1437-1455 (2017) ·Zbl 1367.65084号
[49] 诺瓦蒂,P。;Russo,M.R.,基于GCV的Arnoldi-Tikhonov正则化方法,BIT,54,501-521(2014)·Zbl 1317.65104号
[50] 佩奇,C.C。;桑德斯,M.A.,稀疏不定线性方程组的求解,SIAM J.Numer。分析。,12, 617-629 (1975) ·Zbl 0319.65025号
[51] 佩奇,C.C。;Saunders,M.A.,LSQR:稀疏线性方程组和稀疏最小二乘算法,ACM Trans。数学。软质。,8, 43-71 (1982) ·Zbl 0478.65016号
[52] 赖切尔,L。;Rodriguez,G.,离散不定问题的新旧参数选择规则,数值。算法,63,65-87(2013)·兹比尔1267.65045
[53] 赖切尔,L。;Ye,Q.,奇异系统的无故障GMRES,SIAM J.Matrix Ana。申请。,26, 1001-1021 (2005) ·Zbl 1086.65030号
[54] 赖切尔,L。;Ye,Q.,简单平方平滑正则化算子,Electron。事务处理。数字。分析。,33, 63-83 (2009) ·Zbl 1171.65033号
[55] Ringrose,J.R.,《紧凑非自伴算子》(1971),Van Nostrand Reinhold:Van Nostrand Reinho尔德London·Zbl 0223.47012号
[56] Saad,Y.,《稀疏线性系统的迭代方法》(2003),SIAM:SIAM Philadelphia·Zbl 1002.65042号
[57] 萨阿德,Y。;Schultz,M.H.,GMRES:求解非对称线性系统的广义最小残差法,SIAM J.Sci。统计计算。,7, 856-869 (1986) ·Zbl 0599.65018号
[58] 斯特雷拉,V.V。;Trytyshnikov,E.E.,哪个循环预处理器更好?,数学。竞争。,65, 137-150 (1996) ·Zbl 0847.65029号
[59] Tyrtyshnikov,E.E.,最优和超最优循环预处理,SIAM J.矩阵分析。申请。,13, 459-473 (1992) ·Zbl 0774.65024号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。