使用matlab的符号计算功能即可。
使用到的函数:
eye 生成单位矩阵
det 求矩阵行列式
simplify 符号量化简
示例代码:
syms x
A=[1 2 0;2 2 -2;0 -2 3]%定义一个矩阵
simplify(det(A-eye(3)*x))%求出并展示其特征多项式
运行结果为:
A =
1 2 0
2 2 -2
0 -2 3
ans =
- x^3 + 6*x^2 - 3*x - 10
举个例子:
>> A=[3 7 -3;-2 -5 2; -4 -10 3];
>> AA=sym(A);
>> poly(AA)
ans =
x^3 - x^2 + x - 1