×

关于大规模微分矩阵方程求解器的(L D L ^ T)因式分解的好处。 (英语) Zbl 1320.65110号

摘要:我们提出了求解大规模矩阵微分方程的有效算法。特别地,我们讨论了微分Riccati方程(DRE),并说明了微分Lyapunov方程(DLE)的适用性。我们将重点放在矩阵设置中基于标准版本的常微分方程的方法上。这些方法的应用产生了代数李亚普诺夫方程(ALE),每个步骤都有一定的结构要求解。交替方向隐式(ADI)算法和基于Krylov子空间的方法可以利用这种特殊结构。然而,直接应用经典低阶公式需要使用复杂算法。使用右侧和方程解的(LDL^T)型分解,我们避免了这个问题。因此,对于应用中出现的大规模问题,所提出的方法是一种更实用的替代方法。此外,它们使高阶方法的应用变得可行。数值结果表明,与早期的公式相比,所提出的方法具有更好的性能。

MSC公司:

65升06 常微分方程的多步、Runge-Kutta和外推方法
15A23型 矩阵的因式分解
93甲15 大型系统
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 本纳,P。;Tomljanović,Z。;Truhar,N.,使用降维法优化选定本征频率的阻尼,数值。线性代数应用。,20, 1-17 (2013) ·Zbl 1289.70028号
[2] 布兰奇尼,F。;卡萨格伦德,D。;加多尼奥,P。;Miani,S.,《振动结构半主动阻尼的恒定增益和开关增益》,国际。J.Control,85,1886-1897(2012)·Zbl 1253.93137号
[3] Hoepffner,J.,《随机激励下剪切流的稳定性和控制》(2006年),KTH皇家理工学院:KTH瑞典皇家理工大学博士论文
[4] 宾夕法尼亚州。;弗鲁盖尔,M。;Chang,P.,随机强迫在调节ENSO可预测性中的作用,《气候杂志》,17,3125-3140(2004)
[5] 宾夕法尼亚州。;Sardeshmukh,P.D.,热带海面温度异常的最佳增长,《气候杂志》,1999-224(1995)
[6] Sandberg,H.,线性时变系统的模型降阶(2004),隆德理工学院:瑞典隆德技术学院,博士论文
[7] Sandberg,H.,线性时变系统模型约简的案例研究,Automatica,43,34646-472(2006)·兹比尔1121.93331
[8] Abou-Kandil,H。;Freiling,G。;伊奥内斯库,V。;Jank,G.,《控制与系统理论中的矩阵Riccati方程》(2003),Birkhäuser:Birkháuser Basel,Switzerland·Zbl 1027.93001号
[9] 一川,A。;Katayama,H.,关于时变Riccati方程的备注,系统控制快报。,37, 5, 335-345 (1999) ·Zbl 1050.93507号
[10] Jacobs,O.L.R.,《控制理论导论》(1993),牛津科学出版社:牛津科学出版社·Zbl 0818.93001号
[11] 彼得森,I.R。;乌格里诺夫斯基,V.A。;Savkin,A.V.,《使用(H^\infty)方法的鲁棒控制设计》(2000),Springer-Verlag:Springer-Verlag London,UK·Zbl 0963.93003号
[12] Choi,C。;Laub,A.J.,求解刚性Riccati微分方程的高效矩阵值算法,IEEE Trans。自动化。控制,35,770-776(1990)·Zbl 0714.93011号
[13] 戴维森·E·J。;Maki,M.C.,矩阵Riccati微分方程的数值解,IEEE Trans。自动化。控制,18,71-73(1973)·Zbl 0263.93031号
[14] Dieci,L.,微分Riccati方程的数值积分和一些相关问题,SIAM J.Numer。分析。,2981-815(1992年)·Zbl 0768.65037号
[15] 肯尼,C。;Leipnik,R.B.,微分矩阵Riccati方程的数值积分,IEEE Trans。自动化。控制,30962-970(1985)·Zbl 0594.65054号
[16] 本纳,P。;Mena,H.,无限维LQR问题和相关微分Riccati方程的数值解(2012),Max Planck Iinstitute Magdeburg,可从
[17] 本纳,P。;Mena,H.,Rosenbrock求解微分Riccati方程的方法,IEEE Trans。自动化。控制,58,11,2950-2957(2013)·Zbl 1369.65088号
[18] Hansen,E。;Stillfjord,T.,抽象微分Riccati方程分裂的收敛性分析,SIAM J.Numer。分析。,523128-3139(2014年)·Zbl 1405.65082号
[19] Mena,H.,抛物型偏微分方程最优控制问题中微分Riccati方程的数值解(2007),Escuela Politecnica Nacional博士论文
[20] Antoulas,A.,《大尺度动力系统的近似》,高级设计。Control(2005),工业和应用数学学会:美国宾夕法尼亚州费城工业和应用算术学会·Zbl 1112.93002号
[21] Datta,B.,《线性控制系统设计和分析的数值方法》(2003),爱思唯尔学术出版社
[22] Kjelgard-Mikkelsen,C.,大型Lyapunov方程的数值方法(2009),普渡大学博士论文
[23] Li,J.R。;White,J.,Lyapunov方程的低阶解,SIAM J.矩阵分析。申请。,24, 1, 260-280 (2002) ·Zbl 1016.65024号
[24] Penzl,T.,大型稀疏Lyapunov方程的循环低秩Smith方法,SIAM J.Sci。计算。,21, 4, 1401-1418 (2000) ·Zbl 0958.65052号
[25] Simoncini,V.,求解大规模Lyapunov矩阵方程的新迭代方法,SIAM J.Sci。计算。,29, 3, 1268-1288 (2007) ·Zbl 1146.65038号
[27] 西蒙西尼,V。;德鲁斯金,V。;Knizhnerman,L.,有理Krylov子空间的分析和求解Lyapunov方程的ADI方法,SIAM J.Numer。分析。,49, 5, 1875-1898 (2011) ·Zbl 1244.65060号
[28] Vandereycken,B。;Vandewalle,S.,计算Lyapunov方程低阶解的黎曼优化方法,SIAM J.矩阵分析。申请。,31, 5, 2553-2579 (2010) ·Zbl 1221.65108号
[29] 阿舍尔,U.M。;Petzold,L.R.,《常微分方程和微分代数方程的计算机方法》(1998),SIAM:SIAM Philadelphia,PA·Zbl 0908.65055号
[30] 本纳,P。;Saak,J.,大型稀疏连续时间代数矩阵Riccati和Lyapunov方程的数值解:最新综述,GAMM-Mitt。,36, 1, 32-52 (2013) ·Zbl 1279.65044号
[31] 海尔,E。;Wanner,G.,《求解常微分方程II-Stiff和微分代数问题》,Springer Ser。计算。数学。(2000),Springer-Verlag:纽约Springer-Verlag
[32] 布隆,J.G。;亨茨多夫,W。;斯佩,E.J。;Verwer,J.G.,应用于光化学弥散问题的二阶Rosenbrock方法,SIAM J.Sci。计算。,20, 4, 1456-1480 (1999) ·Zbl 0928.65116号
[33] 霍恩,R.A。;Johnson,C.R.,矩阵分析(1985),剑桥大学出版社:剑桥大学出版社·Zbl 0576.15001号
[34] 安德森,E。;Bai,Z。;比肖夫,C。;德梅尔,J。;Dongarra,J。;杜克罗兹,J。;格林鲍姆,A。;Hammarling,S。;麦肯尼,A。;Sorensen,D.,《LAPACK用户指南》(1999),SIAM:SIAM Philadelphia,PA·Zbl 0934.65030号
[35] 本纳,P。;李,R.-C。;Truhar,N.,《关于Sylvester方程的ADI方法》,J.Compute。申请。数学。,233, 4, 1035-1045 (2009) ·Zbl 1176.65050号
[36] 本纳,P。;Kürschner,P。;Saak,J.,低阶Cholesky因子ADI方法中复杂移位参数的有效处理,Numer。算法,62,2,225-251(2013)·Zbl 1267.65047号
[37] 本纳,P。;Kürschner,P。;Saak,J.,大型Lyapunov和Sylvester方程ADI方法中的自生成和有效移位参数,电子。事务处理。数字。分析。,43, 142-162 (2014) ·Zbl 1312.65068号
[38] 本纳,P。;Kürschner,P。;Saak,J.,对称二阶系统平衡截断的改进数值方法,数学。计算。模型。动态。系统。,19, 6, 593-615 (2013) ·兹比尔1305.93043
[39] Druskin,V.公司。;Simoncini,V.,大型动力系统的自适应有理Krylov子空间,系统控制快报。,60, 8, 546-560 (2011) ·Zbl 1236.93035号
[40] 博尔霍夫,M。;Eppler,A.,广义投影Lyapunov方程的低秩Cholesky因子Krylov子空间方法
[41] Locatelli,A.,最优控制(2001),Birkhäuser:Birkháuser巴塞尔,波士顿,柏林·1096.49500兹罗提
[42] 本纳,P。;Saak,J.,《钢型材最佳冷却的半离散传热模型》(Benner,P.;Mehrmann,V.;Sorensen,D.,《大尺度系统的降维》,Lect.Notes计算科学工程(2005),Springer-Verlag:Springer-Verlag Berlin/Heidelberg,德国),353-356·Zbl 1170.80341号
[43] Shampine,L.F.,《Rosenbrock方法的实施》,ACM Trans。数学。软件,8,2,93-103(1982)·Zbl 0483.65041号
[44] Abels,J。;Benner,P.,CAREX-连续时间代数Riccati方程的基准示例集合(2.0版),工作注释1999-14,SLICOT(1999年11月),可从
[45] 彭兹尔,T。,利亚帕克《用户指南》(2000),TU Chemnitz:TU Chennitz 09107 Chemnitz,Germany,网址:
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。