exm2_f.m

来自「matlab7.0从入门到精通(光盘内容)」· M 代码 · 共 22 行

function S = exm2_f(A)
[m,n] = size(A);
if m ~= n
    error('输入矩阵应是方阵');
end
e = eig(A);
%检查输入矩阵的特征值是否各异
same = 0;
for i = 1:m-1
   for j = (i+1):m
       if e(j) == e(i)
          same = 1;
       end
   end
end
%A可以对角化的条件是A具有各异特征值或者A位埃米尔特矩阵。
if any(any((A'-A)))&(same ==1)
     error('矩阵无法对角化');
end
[v,d] = eig(A);
S = v;