scipy.linalg.fiedler_companion¶
- scipy.linalg.fiedler_companion(a)[源代码]¶
返回Fiedler伴随矩阵
给定一个多项式系数数组
a
,此函数形成具有特殊结构的五对角矩阵,其特征值与的根重合a
。- 参数
- a(n,)类似数组
前导系数非零的多项式系数降序的一维数组。为
N < 2
则返回空数组。
- 退货
- c(N-1,N-1)ndarray
生成的伴阵
参见
注意事项
类似于
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])