×

\计算线性微分代数方程的Lyapunov和Sacker–Sell谱区间的(QR)方法和误差分析。 (英语) Zbl 1242.65150号

作者提出了新的离散和连续(QR)方法,用于计算线性DAE的全部或几个Lyapunov和Sacker-Shell谱区间。特别是,他们考虑了无stangenes-free微分代数方程(DAEs),它们具有积分分离特性。文中提供了详细的算法,并讨论了与文献中已提出的其他方法的比较。详细的扰动和误差分析使本文引人注目。在得出结论之前,最后一节将介绍重要的数值示例,这些示例表明,连续(QR)算法通常比离散算法更精确,计算成本相当。

MSC公司:

65升80 微分代数方程的数值方法
65升07 常微分方程解稳定性的数值研究
34D08型 常微分方程的特征和Lyapunov指数
34A09号 隐式常微分方程,微分代数方程
65升70 常微分方程数值方法的误差界
34D09型 常微分方程解的二分法、三分法
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 亚得里亚诺娃:微分方程线性系统导论。输入:事务处理。数学。专著,第146卷。AMS,普罗维登斯(1995)·Zbl 0844.34001号
[2] Ascher,U.M.,Petzold,L.R.:常微分和微分代数方程的计算机方法。SIAM,费城(1998)·Zbl 0908.65055号
[3] Balla,K.,Linh,V.H.:微分代数方程和哈密顿系统的伴随对。申请。数字。数学。53, 131–148 (2005) ·Zbl 1081.34004号 ·doi:10.1016/j.apnum.2004.08.015
[4] Benettin,G.、Galgani,G.,Giorgilli,L.、Strelcyn,J.M.:光滑动力系统和哈密顿系统的Lyapunov指数;一种计算所有这些参数的方法。第一部分:理论和第二部分:数值应用。麦加尼卡15、9–20、21–30(1980)·Zbl 0488.70015号
[5] Benner,P.,Byers,R.:矩阵铅笔的算法:理论和新算法。数字。数学。103539–573(2006年)·Zbl 1101.65046号 ·数字对象标识代码:10.1007/s00211-006-0001-x
[6] Bohl,P.:《差异》。J.f.d.Reine和Angew。数学。144, 284–313 (1913) ·JFM 45.0498.01标准
[7] Brenan,K.E.,Campbell,S.L.,Petzold,L.R.:微分代数方程初值问题的数值解,第2版。SIAM,费城(1996)·兹比尔0844.65058
[8] Bridges,T.J.,Reich,S.:计算Stiefel流形上的Lyapunov指数。《物理学》,D 156,219–238(2001)·Zbl 0979.37044号 ·doi:10.1016/S0167-2789(01)00283-4
[9] Campbell,S.L.:DAE沿轨迹的线性化。Z.安圭。数学。物理学。46, 70–84 (1995) ·Zbl 0837.34007号 ·doi:10.1007/BF00952257
[10] Chern,J.-L.,Dieci,L.:一些矩阵分解的光滑性和周期性。SIAM J.矩阵分析。申请。22, 772–792 (2000) ·Zbl 0981.65048号 ·doi:10.1137/S0895479899353622
[11] Chyan,C.J.,Du,N.H.,Linh,V.H.:关于线性时变微分代数系统指数稳定性和稳定性半径的数据依赖性。J.差异Equ。245, 2078–2102 (2008) ·Zbl 1162.34004号 ·doi:10.1016/j.jde.2008.07.016
[12] Cong,N.D.,Nam,H.:线性微分代数方程的Lyapunov不等式。数学学报。越南。28, 73–88 (2003) ·Zbl 1048.34089号
[13] Cong,N.D.,Nam,H.:指数为1的线性微分代数方程的Lyapunov正则性。数学学报。越南。29, 1–21 (2004) ·Zbl 1080.34039号
[14] Daleckii,J.L.,Krein,M.G.:Banach空间中微分方程解的稳定性。美国数学学会,普罗维登斯(1974)
[15] Dieci,L.,Eirola,T.:关于矩阵的光滑分解。SIAM J.矩阵分析。申请。20, 800–819 (1999) ·Zbl 0930.15014号 ·doi:10.1137/S0895479897330182
[16] Dieci,L.,Russell,R.D.,Van Vleck,E.S.:幺正积分器及其在连续正交归一化技术中的应用。SIAM J.数字。分析。31, 261–281 (1994) ·Zbl 0815.65096号 ·doi:10.1137/0731014
[17] Dieci,L.,Russell,R.D.,Van Vleck,E.S.:关于连续动力系统Lyapunov指数的计算。SIAM J.数字。分析。34, 402–423 (1997) ·Zbl 0891.65090号 ·doi:10.1137/S0036142993247311
[18] Dieci,L.,Van Vleck,E.S.:连续和离散动力系统的几个Lyapunov指数的计算。申请。数字。数学。17, 275–291 (1995) ·Zbl 0834.65069号 ·doi:10.1016/0168-9274(95)00033-Q
[19] Dieci,L.,Van Vleck,E.S.:Lyapunov和其他光谱:一项调查。收录于:《在离散化条件下保持稳定性的讲座汇编》(科罗拉多州柯林斯堡,2001年),第197-218页。SIAM,费城(2002)
[20] Dieci,L.,Van Vleck,E.S.:李亚普诺夫谱间隔:理论与计算。SIAM J.数字。分析。40, 516–542 (2002) ·Zbl 1021.65067号 ·doi:10.1137/S0036142901392304
[21] Dieci,L.,Van Vleck,E.S.:关于用QR方法计算Lyapunov指数的错误。数字。数学。101, 619–642 (2005) ·Zbl 1083.65069号 ·doi:10.1007/s00211-005-0644-z
[22] Dieci,L.,Van Vleck,E.S.:用QR方法近似Lyapunov指数的微扰理论。J.戴恩。差异Equ。18, 815–842 (2006) ·兹比尔1166.34319 ·doi:10.1007/s10884-006-9024-3
[23] Dieci,L.,Van Vleck,E.S.:Lyapunov和Sacker–Sell光谱区间。J.戴恩。差值等于。19, 265–293 (2007) ·Zbl 1130.34027号 ·doi:10.1007/s10884-006-9030-5
[24] Dieci,L.,Van Vleck,E.S.:关于QR积分中的错误。SIAM J.数字。分析。46, 1166–1189 (2008) ·Zbl 1170.65068号 ·数字对象标识码:10.1137/06067818X
[25] Diehl,M.、Leineweber,D.B.、Schäfer,A.、Bock,H.G.、Schlöder,J.P.:利用回收的废料切割优化多段分批精馏。AIChE J.48(12),2869–2874(2002)·数字对象标识代码:10.1002/aic.690481214
[26] Diehl,M.,Uslu,I.,Findeisen,R.,Schwarzkopf,S.,Allgöwer,F.,Bock,H.G.,Bürner,T.,Gilles,E.D.,Kienle,A.,Schlöder,J.P.,Stein,E.:大型过程的实时优化:高纯度精馏塔的非线性模型预测控制。收录于:Grötschel,M.,Krumke,S.O.,Rambau,J.(编辑)《大型系统的在线优化:最新技术》,第363–384页。施普林格,海德堡(2001)·兹比尔0999.93056
[27] Eich-Soellner,E.,Führer,C.:多体系统中的数值方法。图布纳,斯图加特(1998)·兹比尔0899.70001
[28] Geist,K.,Parlitz,U.,Lauterborn,W.:计算Lyapunov指数的不同方法的比较。掠夺。定理。物理学。83, 875–893 (1990) ·Zbl 1058.65512号 ·doi:10.1143/PTP.83.875
[29] Golub,G.H.,Van Loan,C.F.:《矩阵计算》,第三版。约翰·霍普金斯大学出版社,巴尔的摩(1996)·Zbl 0865.65009号
[30] Greene,J.M.,Kim,J.-S.:李亚普诺夫谱的计算。《物理学D》第24卷,第213-225页(1983年)·Zbl 0619.34047号 ·doi:10.1016/0167-2789(87)90076-5
[31] Griepentlog,E.,März,R.:微分代数方程及其数值处理。Teubner,Leipzig(1986年)·Zbl 0629.65080号
[32] Günther,M.,Feldmann,U.:基于CAD的工业电路建模I.网络方程的数学结构和指数。Surv公司。数学。Ind.8,97–129(1999)·Zbl 0923.65039号
[33] Günther,M.,Feldmann,U.:工业II中基于CAD的电路建模。电路配置和参数的影响。Surv公司。数学。Ind.8131-157(1999)·Zbl 0923.65040号
[34] Hairer,E.,Wanner,G.:求解常微分方程II:刚性和微分-代数问题,第2版。施普林格,柏林(1996)·Zbl 0859.65067号
[35] Kunkel,P.,Mehrmann,V.:奇异线性微分代数方程类的特征。选举人。J.Lin.阿尔及利亚。13, 359–386 (2005) ·Zbl 1102.34003号
[36] Kunkel,P.,Mehrmann,V.:微分代数方程。分析和数值求解。EMS,苏黎世(2006)·Zbl 1095.34004号
[37] Kunkel,P.,Mehrmann,V.:微分代数方程和自旋稳定离散化的稳定性。选举人。事务处理。数字。分析。26, 385–420 (2007) ·Zbl 1171.65418号
[38] Kunkel,P.,Mehrmann,V.,Rath,W.,Weickert,J.:线性微分代数方程的新软件包。SIAM J.科学。计算。18, 115–138 (1997) ·Zbl 0868.65041号 ·doi:10.1137/S1064827595286347
[39] Kunkel,P.,Mehrmann,V.,Seidel,S.:一般非线性微分代数方程数值解的MATLAB包。《技术报告16/2005》,德国柏林TU Berlin Mathematik研究所(2005)。http://www.math.tu-berlin.de/prints网站/
[40] Linh,V.H.,Mehrmann,V.:Lyapunov,Bohl和Sacker–微分代数方程的Sell谱区间。J.戴恩。差异Equ。21, 153–194 (2009) ·Zbl 1165.65050号 ·doi:10.1007/s10884-009-9128-7
[41] Lyapunov,A.M.:运动稳定性的一般问题。A·T·富勒(A.T.Fuller)翻译自爱德华·达沃(Edouard Davoux)于1907年翻译的1892年俄语原文的法语译文。《国际控制杂志》,55,521–790(1992)
[42] Otter,M.、Elmqvist,H.、Mattson,S.E.:使用modelica进行多域建模。收录于:Fishwick,P.(ed.)《CRC动态系统建模手册》。CRC,博卡拉顿(2006)
[43] Rabier,P.J.,Rheinboldt,W.C.:微分代数方程的理论和数值分析。IN:《数值分析手册》,第八卷,第183-540页。荷兰北部,阿姆斯特丹(2002)·兹比尔1007.65057
[44] Riaza,R.:微分代数系统。分析方面和电路应用。世界科学,哈肯萨克(2008)·Zbl 1184.34004号
[45] Riaza,R.:铅笔特征值散度导致拟线性DAE中的稳定性损失。SIAM J.数学。分析。41, 2226–2245 (2010) ·Zbl 1217.34012号 ·数字对象标识码:10.1137/08072557X
[46] Sacker,R.J.,Sell,G.R.:线性微分系统的谱理论。J.差异Equ。27, 320–358 (1978) ·Zbl 0372.34027号 ·doi:10.1016/0022-0396(78)90057-8
[47] Stewart,G.W.:矩阵的$\(\backslash\)mathit{QR}$因式分解的扰动界。SIAM J.Numer。分析。14, 509–518 (1977) ·兹比尔0358.65038 ·doi:10.1137/0714030
[48] Van Vleck,E.S.:关于乘积$\(\backslash\)mathit{QR}$分解中的错误。SIAM J.矩阵分析。申请。31, 1775–1791 (2010) ·Zbl 1206.65147号 ·doi:10.1137/090761562
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。