求矩阵特征多项式的Leverrier方法

引用论文：求特征多项式的Leverrier方法

摘要：本文向读者提供一种既简单又常用的求方阵的特征多项式的方法——leverrier方法。此法比起展开|λI-A|来要简单得多，而且运用此法求的特征多项式的同时，还能求出A的行列式|A|及A的逆矩阵A^(-1)（如果的逆阵存在的话）。

 

matlab程序

以4阶方阵矩阵为例

%求矩阵特征多项式的Leverrier方法
clc,clear
A=[-2 0 1 1;1 -1 1 2;1 2 -1 1;1 1 1 2];
B=[1 2;0 1;2 1;1 -1];
C=[1 1 0 1;0 2 1 2];
n=4;
I4=eye(n);
B3=I4;
a3=-trace(A*B3)
B2=a3*I4+A*B3;
a2=-trace(A*B2)/2
B1=a2*I4+A*B2;
a1=-trace(A*B1)/3
B0=a1*I4+A*B1;
a0=-trace(A*B0)/4