ea604.m

来自「一本关于工程线性代数及其matlab编程方面的好书。」· M 代码 · 共 9 行

% 《工程线性代数（MATLAB版）》第六章例题6.4程序ea604
%
A=[3,2,4;2,0,2;4,2,3];      % 输入系数矩阵
f=poly(A),r=roots(f),       % 求特征多项式并求其根（特征值）
for k=1:3
B(k)=real([A-r(k)*eye(3)]);   % 把三个特征值代入方程，去掉产生的虚数误差 
% B=rref(B,1e-12),    % 若因计算误差造成B非奇异,无基础解，加这一句
p=null(B(k),'r')               % 求其基础解系
end