matlab中的矩阵指数

1、从矩阵A开始A = [0 1 2; 0.5 0 1; 2 1 0]Asave = A;

2、按“Enter”键。得图1所示。

5、exmdem2使用了矩阵指数的经典定义。A = Asave;% Taylor series for e旌忭檀挢xp(A)E = zeros(size(A));F = eye(size(A));k = 1;while norm(E+F-E,1) > 0 E = E + F; F = A*F/k; k = k+1;endE2 = E

6、按“Enter”键。得图3所示。

9、比较结果对于这个矩阵，它们都做得同样好。E = expm(Asave);err1 = E - E1err2 = E - E2err3 = E - E3

10、按“Enter”键。得图5所示。

13、这是一个没有完整特征向量集合的矩阵，因此，expmdemo3失败了。A = [-1 1; 0 -1];E1 = expmdemo1(A)E2 = expmdemo2(A)E3 = expmdemo3(A)

14、按“Enter”键。得图7所示。