命令友矩阵

　　函数compan

　　格式A=compan(u)%u为多项式系统向量,A为友矩阵,A的第1行元素为-u(2:n)/u(1),其中u(2:n)为u的第2到第n个元素,A为特征值就是多项式的特征根.

　　例1-14求多项式的友矩阵和根

　　>>u=[10-76];

　　>>A=compan(u)%求多项式的友矩阵

　　A=07-6100010

　　>>eig(A)%A的特征值就是多项式的根

　　ans=-3.00002.00001.0000

　　友矩阵的特点是主对角线上/下方的元素均为1；最后一行/第一行的元素可取任意值；而其余元素均为零.