scipy.linalg.fiedler_companion

scipy.linalg.fiedler_companion(a)[源代码]

返回Fiedler伴随矩阵

给定一个多项式系数数组 a ,此函数形成具有特殊结构的五对角矩阵,其特征值与的根重合 a

参数
a(n,)类似数组

前导系数非零的多项式系数降序的一维数组。为 N < 2 则返回空数组。

退货
c(N-1,N-1)ndarray

生成的伴阵

参见

companion

注意事项

类似于 companion 前导系数应为非零。在前导系数不是1的情况下,在阵列生成之前对其他系数进行重新缩放。为避免数字问题,最好提供一次多项式。

1.3.0 新版功能.

参考文献

1

M·费德勒,“关于伴随矩阵的注记”,“线性代数及其应用”,2003 DOI:10.1016/S0024-3795(03)00548-2

示例

>>> from scipy.linalg import fiedler_companion, eigvals
>>> p = np.poly(np.arange(1, 9, 2))  # [1., -16., 86., -176., 105.]
>>> fc = fiedler_companion(p)
>>> fc
array([[  16.,  -86.,    1.,    0.],
       [   1.,    0.,    0.,    0.],
       [   0.,  176.,    0., -105.],
       [   0.,    1.,    0.,    0.]])
>>> eigvals(fc)
array([7.+0.j, 5.+0.j, 3.+0.j, 1.+0.j])