%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第五章例5-5-8程序a
% 對稱實矩陣的特征值和特征向量
% 西安電子科技大學出版社出版  陳懷琛編著 2007年7月
%
clear,clf,format compact
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=p'*A*p