×

多频振荡二阶微分方程组显式ERKN积分器的误差界。 (英语) 兹比尔1302.65166

总结:数值分析的一个实质性问题是关于误差的调查和估计。本文关注Wu等人(2010)[30]针对多频振荡二阶微分方程组提出的扩展Runge-Kutta-Nyström(ERKN)积分器的误差分析。ERKN积分器是经典Runge-Kutta-Nyström方法的重要推广,在这个意义上,更新和内部阶段都进行了改革,以使ERKN积分器的定量行为适应真实解的振荡特性。通过显式ERKN积分器误差的展开,我们导出了三阶的刚性阶条件,并给出了误差界。{}我们证明了满足刚性阶(p)的显式ERKN积分器与阶收敛,并且对于(M)是对称半正定矩阵的一个重要特殊情况,(||q_n-q(t_n)||\)的误差界与(||M||)无关((||\cdot||)表示欧几里德范数)。误差分析中提供的刚性阶条件允许我们为多频振荡系统设计新的高效显式ERKN积分器。我们提出了一种新的显式三阶多频多维ERKN积分器,具有最小的色散误差和色散误差。数值实验表明,我们的新型显式多频多维ERKN积分器比科学文献中现有的各种有效方法更有效。我们使用第一个问题来证明这些方法对非对称矩阵的性能良好。特别是,对于著名的Fermi-Paca-Ulam问题,我们新的显式ERKN积分器的数值行为支持我们的理论分析。

MSC公司:

65升06 常微分方程的多步、Runge-Kutta和外推方法
65升70 常微分方程数值方法的误差界
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 科恩,D。;海尔,E。;Lubich,C.,《多频振荡微分方程的数值能量守恒》,BIT Numer。数学。,45, 287-305 (2005) ·Zbl 1083.65117号
[2] 科恩,D。;Jahnke,T。;Lorenz,K。;Lubich,C.,《高振荡哈密顿系统的数值积分器:综述》,(Mielke,A.,多尺度问题的分析、建模和模拟(2006),Springer:Springer-Blin),553-576·Zbl 1367.65191号
[3] Conte博士。;Lubich,C.,量子动力学多组态含时Hartree方法的误差分析,数学。国防部。数字。分析。,44, 759-780 (2010) ·Zbl 1192.81125号
[4] Franco,J.M.,Runge-Kutta-NyströM方法,适用于微扰振子的数值积分,计算。物理学。Comm.,147770-787(2002)·Zbl 1019.6500号
[5] Franco,J.M.,基于ARKN方法的振荡系统新方法,应用。数字。数学。,56, 1040-1053 (2006) ·Zbl 1096.65068号
[6] 加西亚·阿奇拉,B。;桑兹·塞尔纳,J.M。;Skeel,R.D.,振荡微分方程的长时间步长方法,SIAM J.Sci。计算。,20, 930-963 (1998) ·Zbl 0927.65143号
[7] Gautschi,W.,基于三角多项式的常微分方程数值积分,Numer。数学。,3, 381-397 (1961) ·Zbl 0163.39002号
[8] González,A.B。;马汀,P。;Farto,J.M.,扰动振荡器数值积分的一个新的Runge-Kutta型方法家族,Numer。数学。,82, 635-646 (1999) ·Zbl 0935.65075号
[9] Grimm,V.,关于应用于振荡二阶微分方程的Gautschi型指数积分器的误差界,Numer。数学。,100, 71-89 (2005) ·兹比尔1074.65097
[10] Grimm,V.,关于振荡二阶微分方程的Gautschi型方法的注释,Numer。数学。,102, 61-66 (2005) ·Zbl 1082.65080号
[11] 格林·V。;Hochbruck,M.,振荡二阶微分方程指数积分器的误差分析,J.Phys。A: 数学。Gen.,39,5495-5507(2006)·Zbl 1093.65078号
[12] 海尔,E。;Lubich,C.,振荡微分方程数值方法的长期能量守恒,SIAM J.Numer。分析。,38, 414-441 (2000) ·Zbl 0988.65118号
[13] 海尔,E。;卢比奇,C。;Roche,M.,通过微分代数方程研究刚性问题的Runge-Kutta方法的误差,BIT,28,678-700(1988)·Zbl 0657.65093号
[14] 海尔,E。;诺塞特,S.P。;Wanner,G.,《求解常微分方程I:非刚性问题》(1993),Springer-Verlag:Springer-Verlag Berlin·Zbl 0789.65048号
[15] 海尔,E。;卢比奇,C。;Wanner,G.,《几何-数值积分:结构保持算法》(2006),Springer-Verlag:Springer-Verlag Berlin,Heidelberg·Zbl 1094.65125号
[16] Hayes,L.J.,Galerkin非矩形区域交替方向方法(使用面片近似),SIAM J.Numer。分析。,18, 627-643 (1987) ·Zbl 0468.65059号
[17] Hochbruck,M。;Lubich,C.,振荡二阶微分方程的Gautschi型方法,数值。数学。,83, 403-426 (1999) ·Zbl 0937.65077号
[18] Hochbruck先生。;Ostermann,A.,半线性抛物问题的显式指数Runge-Kutta方法,SIAM J.Numer。分析。,43, 1069-1090 (2005) ·Zbl 1093.65052号
[19] Hochbruck,M。;奥斯特曼,A。;Schweitzer,J.,指数Rosenbrock型方法,SIAM J.Numer。分析。,47, 786-803 (2009) ·Zbl 1193.65119号
[20] Iserles,A.,《关于高振荡常微分方程离散化方法的全局误差》,BIT,42,561-599(2002)·Zbl 1027.65107号
[21] 李,J。;王,B。;你,X。;Wu,X.,振荡系统的两步扩展RKN方法,计算。物理学。Comm.,182,2486-2507(2011)·Zbl 1261.65078号
[22] Petzold,L.R。;Yen,L.O。;Yen,J.,高振荡常微分方程的数值解,数值学报。,7, 437-483 (1997) ·Zbl 0887.65072号
[23] 范德胡温,P.J。;Sommeijer,B.P.,计算振荡解的减少相位误差的显式Runge-Kutta(-Nyström)方法,SIAM J.Numer。分析。,24, 595-617 (1987) ·Zbl 0624.65058号
[24] 维戈·阿奎尔,J。;Simos,T.E。;Ferrándiz,J.M.,《控制一个或多个频率扰动振荡长期数值积分中的误差增长》,Proc。罗伊。Soc.伦敦Ser。A、 460561-567(2004)·Zbl 1041.65058号
[25] 王,B。;Wu,X.,振荡二阶微分方程组的一种新型高精度保能积分器,Phys。莱特。A、 3761185-1190(2012)·Zbl 1255.70013号
[26] 王,B。;刘凯。;Wu,X.,解高振荡二阶初值问题的Filon型渐近方法,J.Compute。物理。,243, 210-223 (2013) ·Zbl 1349.65219号
[27] 王,B。;吴,X。;Zhao,H.,Novel改进了多维Strömer-Verlet公式,并将其应用于科学计算的四个方面,数学。计算。型号1。,37, 2327-2336 (2013) ·Zbl 1349.65233号
[28] 吴,X。;Wang,B.,适用于振荡系统的多维Runge-Kutta-Nyström方法,计算机。物理学。Comm.,1811955-1962(2010)·Zbl 1217.65140号
[29] 吴,X。;你,X。;Xia,J.,ARKN方法求解振动系统的阶条件,计算。物理学。Comm.,180,2250-2257(2009年)·兹比尔1197.65084
[30] 吴,X。;你,X。;Shi,W。;Wang,B.,振荡二阶微分方程组的ERKN积分器,计算。物理学。Comm.,1811873-1887(2010)·Zbl 1217.65141号
[31] 吴,X。;王,B。;Xia,J.,显式辛多维指数拟合修正Runge-Kutta-Nyström方法,BIT,52,773-795(2012)·Zbl 1258.65068号
[32] 吴,X。;王,B。;Shi,W.,振荡哈密顿系统的高效保能积分器,J.Compute。物理。,235, 587-605 (2013) ·Zbl 1291.65363号
[33] 吴,X。;你,X。;Wang,B.,振荡微分方程的结构保持算法(2013),Springer-Verlag:Springer-Verlag Berlin,Heidelberg·兹比尔1276.65041
[34] Yang,H。;吴,X。;你,X。;Fang,Y.,扰动振子数值积分的扩展RKN型方法,计算。物理学。Comm.,180,1777-1794(2009)·Zbl 1197.65088号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。