矩阵指数:计算 exp(A)*b 其中 A 是实数且对称的-matlab开发

时间:2024-06-21 11:18:26
【文件属性】:

文件名称:矩阵指数:计算 exp(A)*b 其中 A 是实数且对称的-matlab开发

文件大小:3KB

文件格式:ZIP

更新时间:2024-06-21 11:18:26

matlab

该方法基于此处讨论的 CF 近似: http://www.mathworks.com/matlabcentral/fileexchange/22055 如果 A 稀疏且大,则此方法比 MATLAB 的 expm(A)*b 快得多。 矩阵应该是对称的,并且应该没有正特征值(当然,您可以使用 exp(a+b) = exp(a)*exp(b) 来移动矩阵)。 向量 b 可以有几列,但如果它们不能同时可用,则可以存储某些移位系统的 LU 分解,以加速此类矩阵向量乘积的重复计算。 此方法不会显式计算 EXPM(A)。 如果这是您的目标,请使用 expm(A) 或将 b 设置为单位矩阵。 有几篇论文在讨论这个想法,请使用 特雷费森、魏德曼、施梅尔策Talbot 正交和有理近似北京理工大学,2006 作为起点。


【文件预览】:
Matrixexponential.zip

网友评论