×

矩阵符号函数及其倒数的Padé迭代是最优的。 (英语) Zbl 1233.65037号

计算矩阵符号函数的一种常见方法是通过形式\(z_{k+1}=\phi(z_k)\)的有理迭代,对于在\(1)和\(-1)处具有吸引性不动点的有理函数\(\phi[z)\)。这种有理函数的一类是Padé系列迭代定义如下。
给定表示复数的变量\(z)和\(xi=1-z^2),我们得到:\[\文本{符号}(z)=\frac{z}{(z^2)^{1/2}}=\frac{z}{(1-\xi)^{1/2}}=zf(\xi,\]其中\(f(\xi)=(1-\xi,^{-1/2}\)。让\(frac{P_{m,n}(\xi)}{Q_{m、n}\[z_{k+1}=\压裂{z_kP{(m,n)},\]其中,(P_{(m,n)}(z))和(Q_{。基于倒数的迭代表现类似。
给定一个整数(s>1),并允许非负整数(m,n)在所有值上变化,使得(m+n=2s-1),族(z_{k+1}=\phi_{(m,n)}(z_k)是Padé族及其倒易族的并集。作者证明了函数(phi{(m,n)})定义了在局部收敛到(1)和(-1)且阶数至少为(s)的有理迭代中分子和分母的度数之和最小的唯一有理迭代。虽然不能保证这样的迭代可以以最小的代价进行计算,但这个属性对于计算效率来说仍然很重要,因为在一般情况下,使用霍纳格式,计算(frac{a(z)}{b(z){)}所需的算术运算次数与(a)和(b)的次数之和直接相关。

MSC公司:

65小时04 多项式方程根的数值计算
65楼30 其他矩阵算法(MSC2010)
65埃05 复杂分析中数值方法的一般理论(势理论等)
41A21号机组 帕德近似
30立方厘米 多项式、有理函数和一个复变量的其他分析函数的零点(例如,具有有界Dirichlet积分的函数的零点)

软件:

mf工具箱
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] G.A.Baker Jr.,P.Graves-Morris,PadéApproximants,第二版,《数学及其应用百科全书》,第59卷,剑桥大学出版社,1996年。;G.A.Baker Jr.,P.Graves-Morris,PadéApproximants,第二版,《数学及其应用百科全书》,第59卷,剑桥大学出版社,1996年·Zbl 0923.41001号
[2] 巴菲特,X。;Henriksen,C.,On König’s root-finding algorithms,Nonlinearity,16,3,989-1015(2003)·Zbl 1030.37030号
[3] A.Frommer,V.Simoncini,《矩阵函数:模型降阶:理论、研究方向和应用》,数学。印度,第13卷,施普林格,柏林,2008年,第275-303页。;A.Frommer,V.Simoncini,《矩阵函数:模型降阶:理论、研究方向和应用》,数学。印度,第13卷,施普林格,柏林,2008年,第275-303页·Zbl 1153.65333号
[4] Higham,N.J.,矩阵平方根的稳定迭代,Numer。算法,15,2,227-242(1997)·Zbl 0884.65035号
[5] Higham,N.J.,《矩阵的函数:理论与计算》(2008),美国工业与应用数学学会(SIAM):美国宾夕法尼亚州费城·Zbl 1167.15001号
[6] Howland,J.L.,符号矩阵与矩阵特征值的分离,线性代数应用。,49, 221-232 (1983) ·Zbl 0507.15007号
[7] Iannazzo,B.,有理迭代族及其在矩阵pth根计算中的应用,SIAM J.矩阵分析。应用。,30, 4, 1445-1462 (2008) ·Zbl 1176.65054号
[8] 伊安娜佐,B。;Meini,B.,回文矩阵多项式,矩阵函数和积分表示,线性代数应用。,434, 1, 174-184 (2011) ·Zbl 1211.15027号
[9] 卡兰塔里,B。;卡兰塔里,I。;Zaare-Nahandi,R.,《多项式求根及其特征的基本迭代函数族》,J.Compute。申请。数学。,80, 2, 209-226 (1997) ·Zbl 0874.65037号
[10] 肯尼,C。;Laub,A.J.,矩阵符号函数的有理迭代方法,SIAM J.矩阵分析。应用。,12, 2, 273-291 (1991) ·Zbl 0725.65048号
[11] 肯尼,C。;Laub,A.J.,双曲正切恒等式和帕德符号函数迭代的几何,Numer。算法,7111-128(1994)·Zbl 0812.65035号
[12] 肯尼,C.S。;Laub,A.J.,矩阵符号函数,IEEE Trans。自动化。对照,40,8,1330-1348(1995)·Zbl 0830.93022号
[13] 科奇、切廷·K。;贝坦·巴卡洛卢;Shieh,Leang S.,使用连续分式展开计算矩阵符号函数,IEEE Trans。自动化。控制,39,8,1644-1647(1994)·Zbl 0800.93039号
[14] 拉兹基维茨,B。;Zieötak,K.,矩阵扇区函数和矩阵第(p)个根的Padé族迭代,Numer。线性代数应用。,16, 11-12, 951-970 (2009) ·Zbl 1224.65123号
[15] Laub,A.J.,《Riccati方程数值解的不变子空间方法》,Riccati方程式。Riccati方程数值解的不变子空间方法。(1991),《施普林格:柏林施普林格出版社》,第163-196页
[16] Schröder,E.,Ueber unendlich viele Algorithmen zur Auflösung der Gleichungen,数学。安(1870)
[17] E.Schröder,《关于求解方程的无穷多算法》,技术报告TR-92-121,美国马里兰州帕克学院马里兰大学计算机科学系,1992年11月,57页(G.W.Stewart翻译[16])。;E.Schröder,《关于求解方程的无限多算法》,技术报告TR-92-121,马里兰大学计算机科学系,马里兰州大学帕克分校,美国马里兰州,1992年11月,57页(G.W.Stewart翻译[16])。
[18] J.Stoer,R.Bulirsch,《数值分析导论》,第三版,《应用数学文本》,第12卷,Springer-Verlag,纽约,2002年(由R.Bartels,W.Gautschi和C.Witzgall从德语翻译)。;J.Stoer,R.Bulirsch,《数值分析导论》,第三版,《应用数学文本》,第12卷,Springer-Verlag,纽约,2002年(由R.Bartels,W.Gautschi和C.Witzgall从德语翻译)·Zbl 1004.65001号
[19] van den Eshof,J。;Frommer,A。;里佩尔,Th。;希林,K。;van der Vorst,H.A.,QCD重叠算子的数值方法。符号函数和误差边界,计算。物理学。Comm.,146,2,203-224(2002)·Zbl 1008.81508号
[20] 张,Z。;He,R.,计算矩阵不变子空间的Padé迭代的误差分析,Front。数学。中国,4,2,381-404(2009)·Zbl 1183.65039号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。