×

计算时滞微分方程Lyapunov指数的伪谱约简。 (英语) Zbl 1401.37091号

将延迟微分方程(DDE)简化为常微分方程(ODE)系统是基于[D.布雷达等,SIAM J.Appl。动态。系统。第15期,第1期,第1-23页(2016年;兹比尔1352.34101)].
在这里,作者考虑非线性自治DDE\[x^{\素数}(t)=f(x{t}),\]其中,对于\(tau>0)和\(X:=C([-\tau,0];\mathbb{R})\),\(f:X\rightarrow\mathbb{R}\)是光滑的,并且\(X{t}\ in X\)表示\(X_{t}(θ):=X(t+\theta),\ theta\ in[-\tao,0]\)。通过使用伪谱离散化,他们得到了简化的(M+1)ODE系统\[u^{\prime}_{0}(t)=f(u(t))\]\[u^{\prime}_{i}(t)=\sum_{j=0}^{M}d_{i,j}u{j}(t),\]\(i=1,\dots,M\),初始条件为\(u{i}(0)=\varphi(\theta{i})\),\(i=0,1,\ dots,M \)。作者提到,有几篇论文从数值或计算的角度解决了这个问题,但提出的方法通常很难被非专业人士复制,相关算法和代码很少可用或用户友好。作者的意图是填补那些对应用程序感兴趣的人对离线例程的需求与严格方法之间的差距。
作者考虑了DDE的Lyapunov指数的计算问题。为了计算常微分方程的Lyapunov指数,作者选择了离散QR技术。他们在论文末尾提供了Matlab代码,如(i)所示数据qr,实现了计算线性非自治常微分方程Lyapunov指数的离散QR方法(ii)解决DE_MG实现非线性DDE的时间积分和(iii)lrhs_MG(lrhs_MG),实现了线性(化)非自治矩阵的构造。他们提到,这是第一个用于计算DDE的Lypaunov指数的公开可用代码。作者写道,他们提出的方法具有潜在的广泛应用。作者解释并使用Matlab代码,通过举例说明对三种类型的DDE进行的一些测试,所有这些测试都具有单一延迟。值得一提的是,计算的Lyapunov指数收敛于精确指数的严格证明尚不可用,并且超出了当前工作的范围。事实上,ODE近似系统收敛于原始DDE仍然是当前研究的主题。

MSC公司:

37米25 遍历理论的计算方法(不变测度的近似、Lyapunov指数的计算、熵等)
65升03 泛函微分方程的数值方法
65升07 常微分方程解稳定性的数值研究
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] L.Y.Adrianova,《线性微分方程组导论》,第146期,Transl。数学。专著,AMS,普罗维登斯,1995年·Zbl 0844.34001号
[2] H.T.银行;F.Kappel,泛函微分方程的样条逼近,J.Diff.equations,34496-522(1979)·Zbl 0422.34074号 ·doi:10.1016/0022-0396(79)90033-0
[3] A.贝伦;S.Maset,作为抽象柯西问题的常系数线性时滞微分方程的数值解,Numer。数学。,84, 351-374 (2000) ·Zbl 0949.65072号 ·doi:10.1007/s002110050001
[4] A.Bellen和M.Zennaro,《时滞微分方程的数值方法》,《数值数学和科学计算系列》,牛津大学出版社,2003年·Zbl 1038.65058号
[5] G.贝内廷;L.Galgani;A.Giorgilli;J.M.Strelcyn,光滑动力系统和哈密顿系统的Lyapunov指数;一种计算所有这些参数的方法。第一部分:理论,Meccanica,1529-20(1980)·Zbl 0488.70015号 ·doi:10.1007/BF02128236
[6] G.贝内廷;L.Galgani;A.Giorgilli;J.M.Strelcyn,光滑动力系统和哈密顿系统的Lyapunov指数;一种计算所有这些参数的方法。第2部分:数值应用,麦加尼卡,15,21-30(1980)·doi:10.1007/BF02128237
[7] D.Breda,Hilbert空间中的非自治时滞微分方程和Lyapunov指数,Diff.Int.equations,23935-956(2010)·Zbl 1240.34379号
[8] D.布雷达;O.Diekmann;M.Gyllenberg;F.斯卡拉贝尔;R.Vermiglio,非线性时滞方程的伪谱离散化:数值分岔分析的新前景,SIAM J.Appl。动态。系统。,15, 1-23 (2016) ·Zbl 1352.34101号 ·doi:10.1137/15M1040931
[9] D.布雷达;O.迪克曼;D.列西;F.Scarabel,一类非线性更新方程的数值分岔分析,电子。J.资格。理论不同。Equ.、。,65, 1-24 (2016) ·Zbl 1413.65298号
[10] D.布雷达;S.Maset;R.Vermiglio,时滞微分方程特征根的伪谱差分方法,SIAM J.Sci。计算。,27, 482-495 (2005) ·Zbl 1092.65054号 ·doi:10.1137/030601600
[11] D.布雷达;S.Maset;R.Vermiglio,带非局部边界条件导数算子特征值的伪谱逼近,应用。数字。数学。,56, 318-331 (2006) ·Zbl 1099.65064号 ·doi:10.1016/j.apnum.2005.04.011
[12] D.布雷达;S.Maset;R.Vermiglio,线性滞后泛函微分方程演化算子特征值的逼近,SIAM J.Numer。分析。,501456-1483(2012年)·Zbl 1252.34089号 ·doi:10.1137/100815505
[13] D.Breda、S.Maset和R.Vermiglio,《线性时滞微分方程的稳定性——利用MATLAB的数值方法》,《Springer控制、自动化和机器人简讯》,Springer,纽约,2015年·Zbl 1315.65059号
[14] D.布雷达;E.S.Van Vleck,延迟泛函微分方程的近似Lyapunov指数和Sacker-Shell谱,数值。数学。,126, 225-257 (2014) ·Zbl 1284.65082号 ·doi:10.1007/s00211-013-0565-1
[15] M.D.Chekroun;M.Ghil;H.Liu;S.Wang,非线性时滞微分方程的低维Galerkin逼近,离散Contin。动态。美国,36,4133-4177(2016)·Zbl 1336.34088号 ·doi:10.3934/dcds.2016.36.4133
[16] F.Curtain和H.Zwart,《无限维线性系统理论导论》,《应用数学教材》第21卷,Springer-Verlag,纽约,1995年·Zbl 0839.93001号
[17] L.Dieci、M.S.Jolly和E.S.Van Vleck,用于近似Lyapunov指数的数值技术及其实现,J.Compute。非线性动力学6(2010),011003,第7页。
[18] 迭次;罗素(R.D.Russell);E.S.Van Vleck,关于连续动力系统Lyapunov指数的计算,SIAM J.Numer。分析。,34, 402-423 (1997) ·兹伯利0891.65090 ·doi:10.1137/S0036142993247311
[19] 迭次;E.S.Van Vleck,连续和离散动力系统的几个Lyapunov指数的计算,应用。数字。数学。,17, 275-291 (1995) ·Zbl 0834.65069号 ·doi:10.1016/0168-9274(95)00033-Q
[20] 迭次;E.S.Van Vleck,Lyapunov谱区间:理论与计算,SIAM J.Numer。分析。,40, 516-542 (2002) ·Zbl 1021.65067号 ·doi:10.1137/S0036142901392304
[21] 迭次;E.S.Van Vleck,基于Householder和Givens变换的正交积分器,Future Gener。公司。Sy.,19,363-373(2003)·doi:10.1016/S0167-739X(02)00163-2
[22] L.Dieci和E.S.Van Vleck,LESLIS和LESLIL:线性系统Lyapunov指数近似代码,2004,http://www.math.gatech.edu/dieci/software-les.html。
[23] O.Diekmann、S.A.van Gils、S.M.Verduyn Lunel和H.O.Walther,《延迟方程-泛函、复杂和非线性分析》,《应用数学科学》第110期,Springer Verlag,纽约,1995年·Zbl 0826.34002号
[24] J.R.Dormand;P.J.Prince,嵌入式Runge-Kutta公式家族,J.Comput。申请。数学。,6, 19-26 (1980) ·Zbl 0448.65045号 ·doi:10.1016/0771-050X(80)90013-3
[25] K.Engel和R.Nagel,线性发展方程的单参数半群,数学研究生教材,194。Springer-Verlag,纽约,2000年·Zbl 0952.47036号
[26] D.Farmer,无限维动力系统的混沌吸引子,Physica D,4,366-393(1981/82)·Zbl 1194.37052号 ·doi:10.1016/0167-2789(82)90042-2
[27] D.Gottlieb,伪谱切比雪夫方法的稳定性,数学。公司。,36, 107-118 (1981) ·Zbl 0469.65076号 ·doi:10.1090/S0025-5718-1981-0595045-1
[28] D.Gottlieb,M.Y.Hussaini和S.A.Orszag,谱方法的理论和应用,偏微分方程的谱方法,SIAM,费城,弗吉尼亚州汉普顿,1984,1-54·Zbl 0599.65079号
[29] J.K.Hale和S.M.Verduyn Lunel,《泛函微分方程导论》,第2版,《应用数学科学》第99期,Springer Verlag,纽约,1993年·Zbl 0787.34002号
[30] 伊藤(K.Ito);F.Kappel,使用样条的延迟系统的一致可微近似方案,Appl。数学。选择。,23, 217-262 (1991) ·Zbl 0738.65055号 ·doi:10.1007/BF01442400
[31] F.Kappel,《半群和延迟方程》,第152期(Trieste,1984年),《Pitman Res.Notes Math》。序列号。,朗曼科学。哈洛理工学院,1986年·Zbl 0622.47041号
[32] T.H.Koornvinder,带权函数的正交多项式\(1-x)^α(1+x)^-β+Mδ(x+1)+Nδ(x-1)\),Canad。数学。公牛。,27, 205-214 (1984) ·Zbl 0507.33005号 ·doi:10.4153/CBM-1984-030-7
[33] A.M.Lyapunov,运动稳定性的一般问题,国际。《控制杂志》,55,521-790(1992)·Zbl 0786.70001号
[34] M.C.麦基;L.Glass,生理控制系统中的振荡和混沌,《科学》,197287-289(1977)·Zbl 1383.92036号 ·doi:10.1126/science.267326
[35] S.Maset,作为抽象Cauchy问题的延迟泛函微分方程的数值解,J.Compute。申请。数学。,161, 259-282 (2003) ·Zbl 1042.65052号 ·doi:10.1016/j.cam.2003.03.001
[36] A.Pazy,《线性算子半群及其在偏微分方程中的应用》,《应用数学科学》第44期,Springer Verlag出版社,纽约,1983年·Zbl 0516.47023号
[37] A.Prasad,耦合混沌振荡器中的振幅死亡,物理学。版本E,72,056204-10页(2005)
[38] 沙姆平;S.Thompson,在MATLAB中求解DDE,应用。数字。数学。,37, 441-458 (2001) ·Zbl 0983.65079号 ·doi:10.1016/S0168-9274(00)00055-6
[39] D.E.Sigeti,高频功率谱的指数衰减和正Lyapunov指数,Physica D,52,136-153(1995)·Zbl 0888.58024号 ·doi:10.1016/0167-2789(94)00225-F
[40] J.C.Sprott,一个简单的混沌延迟微分方程,Phys。莱特。A、 366397-402(2007)·Zbl 1203.37064号 ·doi:10.1016/j.physleta.2007.01.083
[41] A.斯特凡斯基;A.达布罗夫斯基;T.Kapitaniak,具有时滞的动力系统中最大Lyapunov指数的评估,混沌,孤子和分形。,23, 1651-1659 (2005) ·Zbl 1071.37014号 ·doi:10.1016/S0960-0779(04)00428-X
[42] L.N.Trefethen,MATLAB中的谱方法,软件-环境-工具系列,SIAM,费城,2000年·Zbl 0953.68643号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。