%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第五章例5-5-8程序
% 對稱實矩陣的特征值和特征向量
% 西安電子科技大學出版社出版  陳懷琛編著 2007年7月
%
clear,clf,format compact
%（1）分步計算特征值和特征向量
A=[2,4,9;4,2,4;9,4,18],		% 輸入矩陣參數
f=poly(A)						% 求其特征多項式的系數向量f
r=roots(f)					% 求其特征根
for i=1:3              % 將三個特征根分別代入特征方程，得齊次方程
    p(:,i)=null(A-r(i)*eye(3)); % 求齊次方程的基礎解，即特征向量
end,
p,d=diag(r)            % 列出特征向量矩陣和特征根（對角）矩陣	
%（2）用eig函數求特征值和特征向量的MATLAB語句如下：
[p1,d1] = eig(A)				% 一步求出特征值和特征向量
% 求原矩陣及特征根矩陣的行列式和跡
det(A), det(d), trace(A), trace(d),