×

高效的基于正交矩阵多项式的矩阵指数计算方法。 (英文) Zbl 1211.65052号

摘要:矩阵指数在不同科学领域出现的微分系统的求解中起着基础性的作用。本文提出了一种基于Hermite矩阵多项式展开式的计算矩阵指数的有效方法。与公认的最佳计算方法相比,Hermite级数截断、缩放和平方以及对中间结果应用浮点算术界提供了极好的精度结果。给出了精确算法中近似值的反向误差分析。该分析用于为矩阵的最佳缩放提供理论估计。
基于此方法的两种算法已实现为MATLAB软件功能。他们已经与MATLAB软件功能有趣矩阵指数函数在大多数测试中获得更高的准确性。仔细的成本比较分析矩阵指数函数结果表明,对于某些矩阵范数区间,所提出的算法具有较低的最大代价。数值试验表明,对中间结果应用浮点算术界可以大大降低计算成本,在数值试验中达到相对较高的平均成本,而对于最终的Hermite选定顺序,平均成本仅为4.43%,并且在77.36%的测试矩阵中获得了更好的精度结果。这个MATLAB软件基于最佳埃尔米特矩阵多项式的算法的实现已经在线提供。

MSC公司:

65楼30 其他矩阵算法(MSC2010)
15甲16 矩阵的指数函数和相似函数
34A30型 线性常微分方程组
65升05 常微分方程初值问题的数值方法
PDF格式BibTeX公司 XML格式引用
全文: 内政部 链接

参考文献:

[1] 弗雷泽,R.A。;Duncan,W.J。;Collar,A.R.,《初等矩阵及其在动力学和微分方程中的应用》(1946),麦克米伦出版社:纽约麦克米伦·Zbl 0061.01403号
[2] 赫希,M.W。;Smale,S.,微分方程。微分方程、动力系统和线性代数(1974),学术出版社:纽约学术出版社·Zbl 0309.34001号
[3] Kailath,T.,《线性系统》(1980年),Prentice Hall出版社:Prentice Hall Englewood Cliffs,NJ·Zbl 0458.93025号
[4] 莫勒,C.B。;Loan,C.V.,计算矩阵指数的十九种可疑方法,SIAM Rev.,20,40801-836(1978)·Zbl 0395.65012号
[5] Smith,G.D.,《偏微分方程的数值解:有限差分方法》(1985),牛津大学出版社·兹比尔0576.65089
[6] 莫勒,C.B。;Loan,C.V.,《二十五年后计算矩阵指数的十九种可疑方法》,SIAM Rev.,45,3-49(2003)·Zbl 1030.65029号
[7] Higham,N.J.,《矩阵指数的缩放和平方法重访》,SIAM J.矩阵分析。申请。,26, 4, 1179-1193 (2005) ·Zbl 1081.65037号
[8] Al-Mohy,A.H。;Higham,N.J.,矩阵指数的一种新的缩放和平方算法,SIAM J.矩阵分析。申请。,31, 3, 970-989 (2009) ·Zbl 1194.15021号
[9] Dunford,N。;Schwartz,J.,线性算子,第一部分(1957),《跨科学:跨科学出版社》,纽约
[10] 萨克斯,S。;Zygmund,A.,《分析函数》(1971),爱思唯尔科学出版社:荷兰阿姆斯特丹
[11] 德菲兹,E。;Jódar,L.,Hermite矩阵多项式级数展开的一些应用,J.Compute。申请。数学。,99, 105-117 (1998) ·Zbl 0929.33006号
[12] Jódar,L。;Company,R.,Hermite矩阵多项式和二阶矩阵微分方程,J.近似理论应用。,12,2,20-30(1996)·Zbl 0858.15014号
[13] 帕特森,M.S。;Stockmeyer,L.J.,《关于计算多项式所需的非标量乘法的数量》,SIAM J.Compute。,2, 1, 60-66 (1973) ·Zbl 0262.65033号
[14] Higham,N.J.,《矩阵的函数:理论与计算》(2008),工业与应用数学学会:美国宾夕法尼亚州费城工业和应用数学学会·Zbl 1167.15001号
[15] Higham,N.J.,《数值算法的准确性和稳定性》(2002),工业和应用数学学会:美国宾夕法尼亚州费城·Zbl 1011.65010号
[16] Golub,G.H。;Loan,C.V.,矩阵计算。矩阵计算,约翰·霍普金斯数学研究。科学。(1996),约翰霍普金斯大学出版社·Zbl 0865.65009号
[17] Ward,R.C.,矩阵指数的数值计算及精度估计,SIAM J.Numer。分析。,14, 4, 600-610 (1977) ·Zbl 0363.65031号
[18] 戴维斯,P.I。;Higham,N.J.,计算矩阵函数的Schur-Parlett算法,SIAM J.矩阵分析。申请。,2544-485(2003年)·Zbl 1052.65031号
[19] S.Blackford,J.Dongarra,《LAPACK安装指南》,《LABACK工作说明411》,田纳西大学计算机科学系,1999年。;S.Blackford,J.Dongarra,《LAPACK安装指南》,LAPACK工作说明411,特内西大学计算机科学系,1999年。
[20] N.J.Higham,《MATLAB测试矩阵工具箱》,第237号数值分析报告,曼彻斯特计算数学中心,英国曼彻斯特,1993年12月。;N.J.Higham,《MATLAB测试矩阵工具箱》,第237号数值分析报告,英国曼彻斯特计算数学中心,1993年12月。
[21] T.G.Wright,Eigtool,版本2.1,2009年<http://web.comlab.ox.ac.uk/pseudospectra/eigtool/>;T.G.Wright,Eigtool,版本2.1,2009年<http://web.comlab.ox.ac.uk/pseudospectra/eigtool/>
[22] 我·纳杰菲尔德。;哈维尔,T.F.,矩阵指数导数及其计算,高级应用。数学。,16, 321-375 (1995) ·Zbl 0839.15004号
[23] 肯尼,C.S。;Laub,A.J.,计算矩阵对数和指数的Schur-Fréchet算法,SIAM J.matrix Ana。申请。,19, 3, 640-663 (1998) ·Zbl 0913.65036号
[24] Westreich,D.,计算矩阵指数及其积分的实用方法,Commun。申请。数字。方法,6375-380(1990)·Zbl 0704.65028号
[25] Lu,Y.Y.,计算指数积分器的矩阵函数,J.Compute。申请。数学。,161, 203-216 (2003) ·兹比尔1033.65028
[26] Dieci,L。;Papini,A.,块三角矩阵指数的Padé逼近,线性代数应用。,308, 183-202 (2000) ·Zbl 0958.65050号
[27] Dolan,E.D。;Moré,J.J.,《性能曲线基准优化软件》,数学。程序。,91, 201-213 (2002) ·邮编:1049.90004
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。