矩阵指数:一种计算对称矩阵指数的有理 Chebyshev 方法。-matlab开发

时间:2024-06-19 15:41:54
【文件属性】:

文件名称:矩阵指数:一种计算对称矩阵指数的有理 Chebyshev 方法。-matlab开发

文件大小:4KB

文件格式:ZIP

更新时间:2024-06-19 15:41:54

matlab

该方法基于将负实轴上的指数函数展开为有理切比雪夫级数。 调用rcexpmv(A,v) 返回expm(A)*v 的近似值,即A 的矩阵指数对向量v(或向量块)的作用。 注意这是线性初值问题u'(t) = A*u(t), u(0) = v的解u(1)。当然,你也可以通过设置v来计算全矩阵指数= 眼睛(大小(A))。 方阵 A 是对称的,没有正的特征值! 如果有正特征值,只需使用 A 向左移动足够远 expm(A) = exp(s) * expm( A - s*speye(size(A)) ), s > 0。 如果 A 很大,这种方法比 Matlab 的 expm 快得多。 它需要针对不同的右手边使用 A 的单个移位线性系统的重复解,并且移位是真实的。 这意味着 1) 对于真实数据 b 和 A,要求解的线性系统是真实的,并且2) 直接求解器只会对矩阵进行一次因式分解并重用因式分解(我们只使


【文件预览】:
rcexpmv.zip

网友评论