【功能简介】生成友矩阵。
【语法格式】
A=compan(u)
u为多项式向量,A返回多项式u的友矩阵。A的第一行元素为-u(2:n)/u(1),即u的第2~n个元素除以第1个元素再取相反数。矩阵A的特征值是多项式的根。
【实例3.16】生成多项式(x-1)(x+2)(x-3)=x3-2x2-5x+6的友矩阵。
>> u=[1,-2,-5,6];
>> A=compan(u) %多项式x3-2x2-5x+6的友矩阵
A =
2 5 -6
1 0 0
0 1 0
>> eig(A) %求矩阵A的特征值
ans =
-2.0000
3.0000
1.0000
【实例分析】eig函数求出矩阵的特征值,可以看出A的特征值等于多项式的根。