poly1.m

有关matlab中simulink的应用

function c=poly1(A)
%POLY1 函数用 Fadeev-Faddeva 算法来取矩阵的特征多项式系数。
%   它完全可以取代 MATLAB 本身提供的 POLY 函数，且运算精度更高。

%Designed by Prof D Xue (c) 2000
[nr,nc]=size(A); 
if nc==nr % 给出若为方阵，则用Fadeev算法求特征多项式
   I=eye(nc); R=I; c=[1 zeros(1,nc)]; 
   for k=1:nc 
      c(k+1)=-1/k*trace(A*R);  
      R=A*R+c(k+1)*I;
   end
elseif (nr==1 | nc==1) % 给出为向量时，构造矩阵
   A=A(isfinite(A)); n=length(A); % 除去非数或无界的特征根
   c=[1 zeros(1,n)];
   for j=1:n
      c(2:(j+1))=c(2:(j+1))-A(j).*c(1:j);
   end
else   % 参数有误则给出错误信息
   error('Argument must be a vector or a square matrix.')
end