ea604.m

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

% 《工程线性代数（MATLAB版）》第六章例题6.4程序ea604
% 陈怀琛，高淑萍，杨威合编，电子工业出版社，2007年6月
% 特征值和特征向量
%
A=[3,2,4;2,0,2;4,2,3];      % 输入系数矩阵
f=poly(A),			    % 求特征多项式(步骤1)
r=roots(f),r=real(r)    % 求特征根（步骤2）,并去掉误差造成的虚部
B1= r(1)*eye(3)-A; 		% 后面可能需要插入一条语句
B1=rref(B1,1e-12),		% 为保证矩阵B1奇异性而设的语句
p1=null(B1,'r')   		% 把第一个个特征值代入方程，求其基础解系
B2=r(2)*eye(3)-A; p2=null(B2,'r')   % 求第二个个特征值的基础解系
B3=r(3)*eye(3)-A; p3=null(B3,'r')   % 求第三个个特征值的基础解系