ea604r.m

matlab 在线性代数中的应用 利用matlab实现可以提高对所学知识的理解 更有助于你考研

% 《工程线性代数（MATLAB版）》第六章例题6.4程序ea604
% 陈怀琛，高淑萍，杨威合编，电子工业出版社，2007年6月
% 特征值和特征向量
%
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