eigenvalue_vector.m

来自「matlab基础程序 便于初学者学习」· M 代码 · 共 17 行

% eigenvalue_vector.m
% 求解矩阵的特征值和特征向量
disp('求解矩阵的特征值和特征向量:');
A=magic(5)
[V,D]=eig(A)
detV=det(V)
err=A*V-V*D
% 当矩阵元素和截断误差相当的情况
disp('矩阵元素和截断误差相当:')
B=[ 3     -2      -.9    2*eps
    -2      4       1    -eps
    -eps/4  eps/2  -1     0
   -.5    -.5      .1    1   ];
[VB,DB]=eig(B)
err1=B*VB - VB*DB
[VN,DN]=eig(B,'nobalance')
err2=B*VN - VN*DN