×

计算波核矩阵函数。 (英语) Zbl 1403.65020号

摘要:我们导出了计算任意方阵(A\)的波核函数(cosh\sqrt{A}\)和(mathrm{sinhc}\sqrt}A})的算法,其中(mathrm{sinhc}z=\sin(z)/z\)。该算法基于Padé近似和使用双角度公式。我们证明了对(cosh\sqrt{A})的任何逼近的后向误差都可以用超几何函数显式表示。为了限制向后错误,我们推导并利用了一个新的限制,这个限制比以前通过A.H.Al-Mohy先生N.J.海姆[SIAM J.Matrix Anal.Appl.31,No.3,970-989(2010;Zbl 1194.15021号)]. 选择缩放量和Padé逼近的程度是为了在精确算术中实现\(\cosh\sqrt{A}\)的向后稳定性的情况下最小化计算成本。数值实验表明,该算法在浮点运算中表现出向前稳定的特性,在这方面优于应用于这些函数的通用Schur-Parlett算法。

MSC公司:

65楼30 其他矩阵算法(MSC2010)
15A60型 矩阵范数,数值范围,泛函分析在矩阵理论中的应用
65层60 矩阵指数和相似矩阵函数的数值计算
PDF格式BibTeX公司 XML格式引用
全文: DOI程序

参考文献:

[1] L.V.Ahlfors,复分析:一个复变量的分析函数理论导论第三版,McGraw-Hill,纽约,1979年·Zbl 0395.30001号
[2] A.H.Al-Mohy,计算三角和双曲矩阵函数作用的截断泰勒级数算法,SIAM J.科学。计算。,40(2018年),第A1696–A1713页·Zbl 1391.15031号
[3] A.H.Al-Mohy和N.J.Higham,矩阵指数的一种新的缩放和平方算法,SIAM J.矩阵分析。申请。,31(2009),第970-989页·Zbl 1194.15021号
[4] A.H.Al-Mohy和N.J.Higham,计算矩阵指数的作用,并应用于指数积分器,SIAM J.科学。计算。,33(2011),第488–511页·Zbl 1234.65028号
[5] A.H.Al-Mohy和N.J.Higham,改进的矩阵对数反缩放和平方算法,SIAM J.科学。计算。,34(2012),第C153–C169页·Zbl 1252.15027号
[6] A.H.Al-Mohy、N.J.Higham和S.D.Relton,分别或同时计算矩阵正弦和余弦的新算法,SIAM J.科学。计算。,37(2015),第A456–A487页·Zbl 1315.65045号
[7] M.Aprahamian和N.J.Higham,矩阵反三角函数和反双曲函数:理论和算法,SIAM J.矩阵分析。申请。,37(2016),第1453-1477页·Zbl 1388.15012号
[8] T.克劳森,Ueber die Faálle,wenn die Reihe von der Form(y=1+frac{α}{1}·frac{β}{γ}x+frac}α.α+1}{1。2} ·frac{β.β+1}{γ.γ+1}x^{2}+\text{etc.})在四象限von der形式中(z=1+frac{α'}{1}·frac}β'}{γ'}·。2} ·压裂{β'.β'+1}{γ'.γ'+1}·压裂{δ'.δ'+1}}{ε'.ε'+1}x^{2}+\text{etc.}帽子。J.Reine Angew著。数学。,3(1828年),第89-91页。
[9] E.B.戴维斯,近似对角化,SIAM J.矩阵分析。申请。,29(2008),第1051-1064页·Zbl 1157.65024号
[10] P.I.Davies和N.J.Higham,计算矩阵函数的Schur–Parlett算法,SIAM J.矩阵分析。申请。,25(2003),第464–485页·兹比尔1052.65031
[11] I.P.Gavrilyuk和V.L.Makarov,希尔伯特空间中二阶发展微分方程的显式解和近似解,数字。方法偏微分方程,15(1999),第111–131页·Zbl 0939.65084号
[12] I.P.Gavrilyuk、V.L.Makarov和V.Vasylyk,二阶方程,摘自《抽象微分方程的指数收敛算法》,Springer,Basel,2011年,第127-165页·Zbl 1225.47001号
[13] J.A.Goldstein,线性算子半群及其应用牛津大学出版社,纽约,1985年·Zbl 0592.47034号
[14] P.C.Greiner、D.Holcman和Y.Kannai,与二阶算子相关的波核,杜克数学。J.,114(2002),第329-386页·Zbl 1072.35130号
[15] V.Grimm和M.Hochbruck,三角算子的有理逼近BIT,48(2008),第215-229页·Zbl 1148.65075号
[16] S.Gu¨ttel,矩阵函数的有理Krylov逼近:数值方法和最优极点选择,GAMM米特。,36(2013),第8-31页·Zbl 1292.65043号
[17] N.J.Higham,矩阵函数工具箱, .
[18] N.J.Higham,矩阵的功能:理论与计算,SIAM,费城,2008年·Zbl 1167.15001号
[19] N.J.Higham和P.Kandolf,计算三角和双曲矩阵函数的作用,SIAM J.科学。计算。,39(2017),第A613–A627页·Zbl 1365.65135号
[20] N.J.Higham和L.Lin,矩阵分数次幂及其Frechet导数的改进Schur–Padeí算法,SIAM J.矩阵分析。申请。,34(2013),第1341–1360页·兹比尔1279.65050
[21] N.J.Higham和M.I.Smith,计算矩阵余弦,数字。《算法》,34(2003),第13-26页·Zbl 1033.65027号
[22] N.J.Higham和F.Tisseur,矩阵1-范数估计的分块算法及其在1-范数伪谱中的应用,SIAM J.矩阵分析。申请。,21(2000),第1185-1201页·Zbl 0959.65061号
[23] D.J.Inman,带控制装置的振动Wiley,2006年。
[24] A.Magnus和J.Wynn,关于\(\cos z\)的Pade表,程序。阿米尔。数学。Soc.,47(1975),第361-367页·Zbl 0297.65016号
[25] P.Nadukandi、B.Servan-Camas、P.A.Becker和J.Garcia-Espinosa,半拉格朗日粒子有限元法的耐波性《计算粒子力学》,4(2017),第321-329页。
[27] M.S.Paterson和L.J.Stockmeyer,关于多项式求值所需的非标量乘法数,SIAM J.计算。,2(1973年),第60-66页·Zbl 0262.65033号
[28] J·L·拉米,丢番图-弗罗贝尼乌斯问题,牛津大学出版社,英国牛津,2005年·Zbl 1134.11012号
[29] A.夏皮罗,用于2D隧道环境中运动的自主蜘蛛机器人的设计与控制2003年,以色列理工学院博士论文。
[30] A.夏皮罗,二阶非对称线性机械系统的稳定性及其在机器人抓取中的应用,事务处理。ASME J.应用。机械。,72(2005),第966–968页·Zbl 1111.74628号
[31] G.Strang和S.MacNamara,差分矩阵的函数是Toeplitz和Hankel,SIAM修订版,56(2014),第525–546页·Zbl 1304.65195号
[32] M.Taylor,非对易谐波分析,数学。调查专题。22,美国数学学会,普罗维登斯,RI,1986年·Zbl 0604.43001号
[33] H.齐格勒,结构稳定性原理,Birkha¨user,巴塞尔,1977年·Zbl 0383.70001号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。