逼近矩阵的指数经验(一)exp⁡(A)使用泰勒级数

计算科学 矩阵
2021-12-02 18:18:04

我试图近似矩阵的指数。我想使用一个容差,但我对如何计算误差感到困惑。有什么想法或提示吗?

2个回答

第一:关于这个主题的必读

莫勒、克利夫和查尔斯范贷款。二十五年后计算矩阵指数的十九种可疑方法。” 暹罗评论 45.1 (2003): 3-49。

(如果您想知道,原始论文是

莫勒、克利夫和查尔斯范贷款。“计算矩阵指数的 19 种可疑方法。” 暹罗评论 20.4 (1978): 801-836.)

更多问题:此幻灯片中绑定的以下错误将有所帮助:

对于一个复杂的 n×n矩阵A

Tr,s=(i=0r1i!(As)i)s
然后
eATr,sAr+1sr(r+1)!eA.

我还没有完全考虑清楚,但是您是否考虑过使用最大的特征值?如果是您的矩阵,并且是其特征分解,那么您可以对泰勒级数中的每个项使用以下恒等式: 然后,取最大的特征值,并从中计算拉格朗日余数的上限。显然,这对于 Q 是正交的对称矩阵最有效。但我仍然认为这可能是一个很好的起点。AQΛQ1

An=QΛnQ1.