×

用因子ADI方法计算Sylvester方程的实低阶解。 (英语) Zbl 1364.65093号

摘要:我们研究大型稀疏Sylvester方程的因子交替方向隐式(ADI)迭代。建立了相关Sylvester残差的一个新的低秩表达式,它可以在迭代过程中廉价地计算残差范数,并产生重新计算的因子ADI迭代。还考虑了广义Sylvester方程的应用。
我们还讨论了复杂移位参数的有效处理,并揭示了ADI迭代与这些复杂移位之间的互连。这将进一步修改因子化ADI迭代,该迭代仅使用绝对必要的复杂算术运算和存储,并产生完全由真实数据组成的低阶解因子。
某些线性矩阵方程,例如交叉Gramian Sylvester和Stein方程,实际上是广义Sylvestr方程的特例,我们展示了如何从广义因子ADI迭代中推导出专门定制的低阶ADI迭代。

MSC公司:

65楼30 其他矩阵算法(MSC2010)
15A24号 矩阵方程和恒等式
65层10 线性系统的迭代数值方法
PDF格式BibTeX公司 XML格式引用
全文: DOI程序

参考文献:

[1] 兰卡斯特,P。;Tismenetsky,M.,《矩阵理论》(1985),学术出版社:奥兰多学术出版社·Zbl 0516.15018号
[2] 霍恩,R。;Johnson,C.,《矩阵分析主题》(1991),剑桥大学出版社:剑桥大学出版社·Zbl 0729.15001号
[3] 巴特尔斯,R。;Stewart,G.,矩阵方程的解(A X+X B=C):算法432,Commun。美国医学会,15,9,820-826(1972)·Zbl 1372.65121号
[4] Sorensen,D。;周,Y.,矩阵Sylvester和Lyapunov方程的直接方法,J.Appl。数学。,2003, 277-303 (2003) ·Zbl 1028.65039号
[5] Hammarling,S.,稳定非负定Lyapunov方程的数值解,IMA J.Numer。分析。,2, 303-323 (1982) ·Zbl 0492.65017号
[6] Golub,G.H。;纳什,S。;Van Loan,C.F.,问题的Hessenberg-Schur方法(A X+X B=C),IEEE Trans。自动垫。控制,24,6909-913(1979)·Zbl 0421.65022号
[7] Epton,M.,解(A X D-B X C=E)的方法及其在隐式常微分方程数值解中的应用,BIT,20,341-345(1980)·Zbl 0452.65015号
[8] 加德纳,J。;Laub,A.J。;阿马托,J。;Moler,C.,Sylvester矩阵方程的解(A X B+C X D=E),ACM Trans。数学。软件,18,223-231(1992)·兹伯利0893.65026
[9] Roberts,J.,线性模型简化和使用符号函数求解代数Riccati方程,Internat。《控制杂志》,32,4,677-687(1980)·Zbl 0463.93050号
[10] 本纳,P。;昆塔纳-奥尔蒂,E。;Quintana-Ortyí,G.,通过有理迭代方案求解稳定Sylvester方程,J.Sci。计算。,28, 1, 51-83 (2005) ·Zbl 1098.65041号
[11] Sorensen,D。;安托拉斯,A.,《西尔维斯特方程和近似平衡约简》,《线性代数应用》。,351-352, 0, 671-700 (2002) ·Zbl 1023.93012号
[12] 本纳,P。;科勒,M。;Saak,J.,《(H_2)模型降阶中的稀疏Sylvester方程》,技术代表MPIMD/11-11,马克斯·普朗克研究所马格德堡预印本,(2011年12月),可从http://www.mpi-magdeburg.mpg.de/prints/2011/11/
[13] El Guennouni,A。;Jbilou,K。;Riquet,A.,求解大型Sylvester方程的Block Krylov子空间方法,Numer。算法,29,75-96(2002)·Zbl 0992.65040号
[14] Bao,L。;Lin,Y。;Wei,Y.,广义Sylvester方程的Krylov子空间方法,应用。数学。计算。,175, 1, 557-573 (2006) ·Zbl 1094.65036号
[15] Bao,L。;Lin,Y。;Wei,Y.,求解大型Sylvester方程的一种新投影方法,应用。数字。数学。,57, 5-7, 521-532 (2007) ·Zbl 1118.65028号
[16] Jbilou,K.,大型Sylvester矩阵方程的低阶近似解,应用。数学。计算。,177, 365-376 (2006) ·Zbl 1095.65041号
[17] 胡,D。;Reichel,L.,Sylvester方程的Krylov-子空间方法,线性代数应用。,172, 283-313 (1992) ·Zbl 0777.65028号
[18] Bouhamidi,A。;哈奇,M。;Heyouni,M。;Jbilou,K.,大型Sylvester矩阵方程的预处理块Arnoldi方法,数值。线性代数应用。,20, 2, 208-219 (2013) ·Zbl 1289.65099号
[19] 本纳,P。;Breiten,T.,关于大规模矩阵方程近似低阶解的最优性,SCL,67,55-64(2014)·Zbl 1292.65040号
[20] Wachspress,E.,Lyapunov矩阵方程的迭代解,应用。数学。莱特。,10787-90(1988年)·Zbl 0631.65037号
[21] Calvetti,D。;Reichel,L.,ADI迭代方法在噪声图像恢复中的应用,SIAM J.矩阵分析。申请。,17, 165-186 (1996) ·Zbl 0849.65101号
[22] Sorensen,D。;Zhou,Y.,李亚普诺夫方程解的特征值衰减率和灵敏度的界限,技术代表TR02-07,公司部。申请。数学。,德克萨斯州休斯顿莱斯大学,(2002年6月),可在线访问http://www.caam.rice.edu/caam/trs/tr02.html
[23] Grasedyck,L.,Sylvester方程解的低秩或(H)-矩阵近似的存在性,Numer。线性代数应用。,11, 371-389 (2004) ·Zbl 1164.65381号
[24] Penzl,T.,大型稀疏Lyapunov方程的循环低秩Smith方法,SIAM J.Sci。计算。,21, 4, 1401-1418 (2000) ·Zbl 0958.65052号
[25] Benner,P.,《解决大规模控制问题》,IEEE控制系统。Mag.,14,1,44-59(2004)
[27] 本纳,P。;李,J.-R。;Penzl,T.,大型Lyapunov方程、Riccati方程和线性二次控制问题的数值解,Numer。线性代数应用。,15, 9, 755-777 (2008) ·兹比尔1212.65245
[29] 本纳,P。;李,R.-C。;Truhar,N.,《关于Sylvester方程的ADI方法》,J.Compute。申请。数学。,233, 4, 1035-1045 (2009) ·Zbl 1176.65050号
[30] 李,R.-C。;Truhar,N.,《关于Sylvester方程的ADI方法》,《2008-2年技术报告》,德克萨斯大学阿灵顿分校数学系(2008年),可在http://www.uta.edu/math/print/rep2008_02.pdf
[31] 本纳,P。;Kürschner,P。;Saak,J.,低阶Cholesky因子adi方法中复杂移位参数的有效处理,Numer。算法,62,2,225-251(2013)·Zbl 1267.65047号
[33] 李,J.-R。;White,J.,Lyapunov方程的低阶解,SIAM J.矩阵分析。申请。,24260-280(2002年)·Zbl 1016.65024号
[35] Feitzinger,F。;Hylla,T。;Sachs,E.W.,Riccati方程的非精确Kleinman-Newton方法,SIAM J.矩阵分析。申请。,31, 2, 272-288 (2009) ·Zbl 1191.49033号
[36] 本纳,P。;Kürschner,P。;Saak,J.,对称二阶系统平衡截断的改进数值方法,数学。计算。模型。动态。系统。,19, 6, 593-615 (2013) ·Zbl 1305.93043号
[38] 弗拉格,G.M。;Gugercin,S.,关于Sylvester方程的ADI方法和最佳(H_2)点,Appl。数字。数学。,64,0,50-58(2013)·Zbl 1255.65087号
[39] Day,D。;Heroux,M.A.,通过等效实公式求解复值线性系统,SIAM J.Sci。计算。,23480-498(2001年)·Zbl 0992.65020号
[40] 本纳,P。;Kürschner,P。;Saak,J.,有效避免低阶ADI方法中的复杂算法,Proc。申请。数学。机械。,12, 1, 639-640 (2012)
[41] 本纳,P。;Kürschner,P。;Saak,J.,《具有复杂位移的低阶ADI方法的真实版本》,技术代表MPIMD/12-11,马克斯·普朗克研究所马格德堡预印本(2012),可从http://www.mpi-magdeburg.mpg.de/prints网站/
[42] Sadkane,M.,大型离散Lyapunov方程的低秩Krylov平方Smith方法,线性代数应用。,436, 8, 2807-2827 (2012) ·Zbl 1260.65040号
[43] 本纳,P。;Khoury,G.E。;Sadkane,M.,《关于大规模Stein方程的平方史密斯方法》,Numer。线性代数应用。(2013)
[44] Li,T.等人。;翁,P.C.-Y。;wah Chu,E.K。;Lin,W.-W.,大尺度Stein和Lyapunov方程,Smith方法和应用,数值。算法,63,4,727-752(2013)·Zbl 1271.65079号
[45] 本纳,P。;Faßbender,H.,关于大规模稀疏离散时间Riccati方程的数值解,高级计算。数学。,35, 119-147 (2011) ·Zbl 1230.65070号
[47] 本纳,P。;Kürschner,P。;Saak,J.,《大型Lyapunov和Sylvester方程ADI方法中的自生成和有效移位参数》,技术代表MPIMD/13-18,马克斯·普朗克研究所马格德堡预印本(2013),可从http://www.mpi-magdeburg.mpg.de/prints网站/
[48] 比尼,D.A。;伊安娜佐,B。;Meini,B.,代数Riccati方程的数值解(2011),工业和应用数学学会
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。