ea621b.m

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

% 《工程线性代数（MATLAB版）》第六章例题6.2.1程序ea621b
%
A=[3,2,4;2,0,2;4,2,3];      % 输入系数矩阵
f=poly(A),			       % 求特征多项式
r=roots(f),r=real(r)       % 求特征根（特征值）,并去掉误差造成的虚部
B1=r(1)*eye(3)-A;           % 把第一个个特征值代入方程，
B1o=rref(B1)
B1=rref(B1,1e-12)           % 修正误差，使它确实奇异
p1=null(B1,'r')             % 求其基础解系
B2=r(2)*eye(3)-A; p2=null(B2,'r')   % 把第二个个特征值代入方程，求其基础解系
B3=r(3)*eye(3)-A; p3=null(B3,'r')   % 把第三个个特征值代入方程，求其基础解系