ea621a.m

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

% 《工程线性代数（MATLAB版）》第六章例题6.2.1程序ea621a
%
A=[3,2,4;2,0,2;4,2,3];      % 输入系数矩阵
f=poly(A),r=roots(f),       % 求特征多项式并求其根（特征值）
B1=real([A-r(1)*eye(3)]);   % 把三个特征值代入方程，去掉产生的虚数误差 
p1=null(B1,'r')               % 求其基础解系
B2=real([A-r(2)*eye(3)]);   % 把三个特征值代入方程，去掉产生的虚数误差 
p2=null(B2,'r')               % 求其基础解系
B3=real([A-r(3)*eye(3)]);   % 把三个特征值代入方程，去掉产生的虚数误差 
p3=null(B3,'r')               % 求其基础解系