×

一种基于SVD-Krylov迭代的大规模动力系统模型降阶方法。 (英语) Zbl 1137.93006号

摘要:我们提出了一种用于大规模线性时不变动力系统逼近的模型约简算法。该方法是一种双边投影,结合了基于奇异值分解(SVD)和基于Krylov的模型约简技术的特点。虽然投影的SVD侧取决于可观测性gramian,但Krylov侧是通过迭代有理Krylof步骤获得的。简化模型是渐近稳定的,匹配一定的矩,并解决了一个受限的({mathcal H}_2)最小化问题。我们对在约简步骤中使用低秩gramian和约简离散时间系统的方法进行了修改。来自不同学科的几个数值例子验证了该方法的有效性。它的性能明显优于\(q\)-封面[A.优素福R.E.斯克尔顿《协方差等效实现及其在大型系统模型简化中的应用》,载于:C.T.Leondes(编辑),《控制与动态系统》,第1部分,《控制动力学》。系统。,高级理论应用。22, 273–348 (1985;Zbl 0651.93006号);A.尤瑟夫,D.A.瓦吉R.E.斯克尔顿,J.数学。分析。申请。106, 91–115 (1985;Zbl 0573.93073号)]和最小二乘法[S.古格霉素A.C.安托拉斯,线性代数应用。415,第2–3号,290–321(2006年;Zbl 1112.93015号)]与该方法具有类似投影结构的方法。此外,就({mathcal H}_2)和({mathcal H}_2\)误差度量而言,其性能与平衡截断相当,有时甚至优于平衡截断。此外,由于使用了近似的革兰氏算子,该方法对扰动具有鲁棒性。

MSC公司:

93甲15 大型系统
93B11号机组 系统结构简化
93天20分 控制理论中的渐近稳定性
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 安托拉斯,A.C。;鲍尔,J.A。;Kang,J。;Willems,J.C.,关于最小有理插值问题的解,线性代数应用。,137/138, 511-573 (1990) ·Zbl 0715.93020号
[2] A.C.Antoulas,D.C.Sorensen,平衡模型简化的投影方法,《ECE-CAAM部门技术报告》,莱斯大学,1999年(9月)。;A.C.Antoulas,D.C.Sorensen,平衡模型简化的投影方法,《ECE-CAAM部门技术报告》,莱斯大学,1999年(9月)。
[3] A.C.Antoulas,D.C.Sorensen,Sylvester方程和近似平衡归约,见:V.Blondel,D.Hinrichsen,J.Rosenthal,P.M.van Dooren(编辑),线性系统与控制,线性代数及其应用,第351-3522002卷,第671-700页(第四期特刊)。;A.C.Antoulas,D.C.Sorensen,《Sylvester方程和近似平衡约化》,载于:V.Blondel,D.Hinrichsen,J.Rosenthal,P.M.van Dooren(编辑),《线性系统与控制,线性代数及其应用》,第351-3522002卷,第671-700页(第四期特刊)·Zbl 1023.93012号
[4] 安托拉斯,A.C。;索伦森特区。;Gugercin,S.,《大系统模型简化方法调查》,《当代数学》,第280卷(2001年),AMS出版物,第193-219页·Zbl 1048.93014号
[5] 安托拉斯,A.C。;索伦森特区。;周永康,关于Hankel奇异值的衰减率及相关问题,系统控制快报。,第46323-342页(2002年)·Zbl 1003.93024号
[6] Antoulas,A.C.,《大规模动力系统的近似》。大型动力系统的近似,设计与控制进展DC-06(2005),SIAM:SIAM Philadelphia·Zbl 1112.93002号
[7] Arnoldi,W.E.,矩阵特征值问题求解中的最小迭代原则,Quart。申请。数学。,9, 17-29 (1951) ·Zbl 0042.12801号
[8] Bai,Z。;Day,D。;Ye,Q.,ABLE:非厄米特特征值问题的自适应块Lanczos方法,SIAM J.Matrix Ana。申请。,20, 4, 1060-1082 (1999) ·Zbl 0932.65045号
[9] Bai,Z。;Ye,Q.,通过Lanczos程序对传递函数Padé近似的误差估计,Electron。事务处理。数字。分析。,7, 1-17 (1998) ·Zbl 0913.41013号
[10] 鲍尔,J.A。;戈伯格,I。;Rodman,L.,有理矩阵函数的插值(1990),Birkhäuser Verlag:Birkháuser Verlag Basel·Zbl 0835.41005号
[11] Bartels,R.H。;Stewart,G.W.,《矩阵方程的解(AX+XA=C):算法432》,美国通信协会,第15期,第820-826页(1972年)·Zbl 1372.65121号
[12] C.A.Beattie,动力系统模型简化的投影方法,SIAM计算科学与工程会议,美国奥兰多,2005年(2月)。;C.A.Beattie,动力系统模型简化的投影方法,SIAM计算科学与工程会议,美国奥兰多,2005年(2月)·Zbl 1159.93317号
[13] Boley,D.L.,状态空间控制模型的Krylov空间方法,电路系统信号处理。,13, 733-758 (1994) ·Zbl 0833.93024号
[14] C.A.Beattie,S.Gugercin,《Inexact在基于Krylov的模型约简中的求解》,载于:第45届IEEE CDC会议记录,加利福尼亚州圣地亚哥,2006年12月。;C.A.Beattie,S.Gugercin,《Inexact在基于Krylov的模型约简中的求解》,载于:第45届IEEE CDC会议记录,加利福尼亚州圣地亚哥,2006年12月·Zbl 1236.93034号
[15] 本纳,P。;Quintana-Orti,E.S.,用矩阵符号函数求解稳定的广义Lyapunov方程,Numer。算法,2075-100(1999)·兹伯利0940.65035
[16] Y.Chahlaoui,P.Van Dooren,线性时不变动力系统模型简化基准示例集,SLICOT工作说明2002-2:2002年2月<http://www.win.tue.nl/niconet/NIC2/benchmarks.html>; Y.Chahlaoui,P.Van Dooren,线性时不变动力系统模型简化基准示例集,SLICOT工作说明2002-2:2002年2月<http://www.win.tue.nl/nicone/NIC2/benchmarks.html>
[17] J.K.Cullum,W.E.Donath,计算q代数最大特征值和大型稀疏对称矩阵对应特征空间的块Lanczos算法,收录于:IEEE决策与控制会议论文集,1974年,第505-509页。;J.K.Cullum,W.E.Donath,计算q代数最大特征值和大型稀疏对称矩阵对应特征空间的块Lanczos算法,收录于:IEEE决策与控制会议论文集,1974年,第505-509页。
[18] 德维尔马涅,C。;Skelton,R.,使用投影公式进行模型简化,国际。《控制杂志》,40,2141-2169(1987)
[19] 阿利亚加,J.I。;Boley,D.L。;Freund,R.W。;Hernandez,V.,多起始向量的Lanzcos型方法,数学。公司。,69, 1577-1601 (2000) ·Zbl 0953.65018号
[20] Gaier,D.,《复近似讲座》(1987年),Birkhauser·Zbl 0612.30003号
[21] 加里凡,K。;格里姆,E。;Van Dooren,P.,模型简化的理性Lanczos算法,数值。算法,2,1-2,33-63(1996)·Zbl 0870.65053号
[22] 加里凡,K。;范登多普,A。;Van Dooren,P.,Sylvester方程和基于投影的模型简化,J.Compute。申请。数学。,162, 213-229 (2004) ·Zbl 1044.65054号
[23] 加里凡,K。;范登多普,A。;Van Dooren,P.,通过切向插值对mimo系统进行模型简化,SIAM J.矩阵分析。申请。,26, 2, 328-349 (2004) ·Zbl 1078.41016号
[24] Golub,G.H。;Van Loan,C.F.,《矩阵计算》(1996),约翰霍普金斯大学出版社:约翰霍普金大学出版社巴尔的摩·Zbl 0865.65009号
[25] Glover,K.,线性多变量系统的所有最优Hankel-形式近似及其(L^ infty)-误差界,Internat。J.Control,39,1115-1193(1984)·Zbl 0543.93036号
[26] E.J.Grimme,《模型简化的Krylov投影方法》,博士论文,伊利诺伊大学ECE系,厄本那-香槟分校,1997年。;E.J.Grimme,《模型简化的Krylov投影方法》,博士论文,伊利诺伊大学厄本那-香槟分校ECE系,1997年。
[27] 格里姆·E·J。;索伦森特区。;Van Dooren,P.,通过隐式重启Lanczos方法对状态空间系统进行模型简化,Numer。算法,12,1-31(1995)·Zbl 0870.65052号
[28] S.Gugercin,A.C.Antoulas,《7种模型简化算法的比较研究》,载于《第39届IEEE决策与控制会议论文集》,澳大利亚悉尼,2000年12月。;S.Gugercin,A.C.Antoulas,《7种模型简化算法的比较研究》,载于《第39届IEEE决策与控制会议论文集》,澳大利亚悉尼,2000年12月。
[29] S.Gugercin、A.C.Antoulas、N.Bedrossian,国际空间站1R和12A模型的近似值,载于:第40届CDC会议记录,2001年12月。;S.Gugercin,A.C.Antoulas,N.Bedrossian,国际空间站1R和12A模型的近似,载于:第40届美国疾病控制与预防中心会议记录,2001年12月。
[30] S.Gugercin,《大尺度动力系统模型降阶的投影方法》,莱斯大学ECE系博士论文,2002年12月。;S.Gugercin,《大尺度动力系统模型简化的投影方法》,莱斯大学ECE系博士论文,2002年12月。
[31] 古吉丁,S。;索伦森特区。;Antoulas,A.C.,大型Lyapunov方程的修正低阶Smith方法,数值。算法,32,1,27-55(2003)·Zbl 1034.93020号
[32] S.Gugercin,A.C.Antoulas,安\(\数学{H} _2\); S.Gugercin、A.C.Antoulas、An(mathcal){H} _2\)
[33] 古吉丁,S。;Antoulas,A.C.,用最小二乘法对大系统进行模型简化,线性代数应用。,415, 2-3, 290-321 (2006) ·Zbl 1112.93015号
[34] 古吉丁,S。;Li,J.R.,大型稀疏系统平衡截断的Smith型方法,(Benner,P.;Golub,G.;Mehrmann,V.;Sorensen,D.C.,《大尺度系统的降维》,《计算科学与工程讲义》,第45卷(2005年),Springer-Verlag:Springer-Verlag Berlin/Heidelberg),国际标准图书编号3-540-24545-6·Zbl 1215.93026号
[35] S.Gugercin、A.C.Antoulas、C.A.Beattie、{H} _2\); S.Gugercin、A.C.Antoulas、C.A.Beattie、{H} _2\) ·Zbl 1159.93318号
[36] S.Gugercin,C.A.Beattie,Inexact在基于Krylov的大型动力系统模型约简中求解,in:第九届铜山迭代方法会议,科罗拉多州铜山,2006年(4月2日至7日)。;S.Gugercin,C.A.Beattie,Inexact在基于Krylov的大型动力学系统模型约简中求解,见:第九届铜山迭代方法会议,科罗拉多州铜山,2006年(4月2日至7日)。
[37] Hammarling,S.,稳定非负定Lyapunov方程的数值解,IMA J.Numer。分析。,2, 303-323 (1982) ·Zbl 0492.65017号
[38] Hodel,A.S。;Poola,K.P。;Tenison,B.,用近似幂迭代法求解Lyapunov方程,线性代数应用。,236, 205-230 (1996) ·Zbl 0848.65033号
[39] 胡大勇。;Reichel,L.,Sylvester方程的Krylov子空间方法,线性代数应用。,172283-313(1992年)·Zbl 0777.65028号
[40] Jaimoukha,I.M。;Kasenally,E.M.,解大型Lyapunov方程的Krylov子空间方法,SIAM J.Numer。分析。,31, 227-251 (1994) ·Zbl 0798.65060号
[41] Kim,H.M。;Craig,R.R.,使用非对称块Lanczos算法进行结构动力学分析,国际。J.数字。方法工程,26,2305-2318(1988)·兹比尔0661.73062
[42] 李,J.-R。;White,J.,Lyapunov方程的低阶解,SIAM J.矩阵分析。申请。,24, 1 (2002) ·Zbl 1016.65024号
[43] 刘,Y。;安德森,B.D.O.,平衡系统的奇异摄动近似,国际。J.Control,第50期,第1379-1405页(1989年)·Zbl 0688.93008号
[44] Lanczos,C.,求解线性微分和积分算子特征值问题的迭代方法,J.Res.Nat.Bur。标准,45,255-282(1950)
[45] Meier,L。;Luenberger,D.G.,线性常数系统的近似,IEEE Trans。自动化。控制,12585-588(1967)
[46] Moore,B.C.,《线性系统的主成分分析:可控性、可观测性和模型简化》,IEEE Trans。自动化。控制,AC-26,17-32(1981)·Zbl 0464.93022号
[47] Mullis,C.T。;Roberts,R.A.,最小舍入噪声定点数字滤波器的合成,IEEE Trans。电路系统I基金理论应用。,CAS-23,551-562(1976)·Zbl 0342.93066号
[48] Nikishin,A.A。;Yeremin,A.Y.,并行计算机上求解大型稀疏对称正定线性系统的可变块CG算法:I:一般迭代格式,SIAM J.矩阵分析。申请。,16, 1135-1153 (1995) ·Zbl 0837.65029号
[49] O'Leary,D.P.,块共轭梯度算法及相关方法,线性代数应用。,29, 293-322 (1980) ·兹比尔0426.65011
[50] Penzl,T.,大型稀疏Lyapunov方程的循环低秩Smith方法,SIAM J.Sci。计算。,21, 4, 1401-1418 (2000) ·Zbl 0958.65052号
[51] T.Penzl,大型动力系统模型简化算法,技术报告SFB393/99-40,Sonderforschungsbereich 393 Numerische Simulation auf massiv parallelen Rechern,TU Chemnitz,09107,FRG,1999<http://www.tu-chemnitz.de/sfb393/sfb99pr.html>; T.Penzl,大型动力系统模型简化算法,技术报告SFB393/99-40,Sonderforschungsbereich 393 Numerische Simulation auf massiv parallelen Rechern,TU Chemnitz,09107,FRG,1999<http://www.tu-chemnitz.de/sfb393/sfb99pr.html>
[52] Penzl,T.,Lyapunov方程解的特征值衰减界:对称情况,系统。控制信函。,40, 139-144 (2000) ·Zbl 0977.93034号
[53] Ruhe,A.,非对称特征值问题的有理Krylov算法II:矩阵对,线性代数应用。,197, 283-295 (1994) ·Zbl 0810.65031号
[54] Ruhe,A.,计算大型稀疏对称矩阵特征值的带Lanczos算法的实现方面,数学。公司。,33, 680-687 (1979) ·Zbl 0443.65022号
[55] Sorensen,D.C.,多项式滤波器在步长Arnoldi方法中的隐含应用,SIAM J.矩阵分析。申请。,13, 357-385 (1992) ·Zbl 0763.65025号
[56] Varga,A.,SLICOT库中的模型简化软件,(Datta,B.,应用和计算控制,信号和电路(2001),Kluwer学术出版社:Kluwer学术出版社波士顿)
[57] Yousouff,A。;Wagie,D.A。;Skelton,R.E.,通过协方差等效实现的线性系统近似,J.Math。分析。申请。,196, 91-115 (1985) ·Zbl 0573.93073号
[58] 优素福,A。;Skelton,R.E.,《协方差等效实现及其在大系统模型降阶中的应用》(Leondes,C.T.,《控制与动态系统》,第22卷(1985年),学术出版社),273-348·Zbl 0651.93006号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。