Python中计算Chebyshev级数的根
要计算多项式的根,请使用Python Numpy中的chebyshev.chebroots()方法。 该方法返回系列根的数组。如果所有根都是实数,则输出也是实数,否则是复数。 参数c是一个1-D系数数组。
根据伴随矩阵的特征值来获得根的估计值,远远离开复平面原点的根可能会产生较大的误差, 因为对于这些值,该级数的数值不稳定。在多重性大于1的根附近,由于该级数在这些点附近的值相对于根的误差比较不敏感, 因此误差会更大。靠近原点的孤立根可以通过几次牛顿迭代法来改善。
步骤
首先,导入所需的库 –
from numpy.polynomial import chebyshev as C
要计算多项式的根,请使用Python Numpy中的chebyshev.chebroots()方法 –
print("Result (roots)...\n",C.chebroots((-1,0,1)))
获取数据类型 –
print("\nType...\n",C.chebroots((-1,0,1)).dtype)
获取形状 –
print("\nShape...\n",C.chebroots((-1,0,1)).shape)
示例
from numpy.polynomial import chebyshev as C
# 要计算多项式的根,请使用Python Numpy中的chebyshev.chebroots()方法。
# 该方法返回系列根的数组。如果所有根都是实数,则输出也是实数,否则是复数。
# 参数c是一个1-D系数数组。
print("Result (roots)...\n",C.chebroots((-1,0,1)))
# 获取数据类型
print("\nType...\n",C.chebroots((-1,0,1)).dtype)
# 获取形状
print("\nShape...\n",C.chebroots((-1,0,1)).shape)
输出
Result (roots)...
[-1. 1.]
Type...
float64
Shape...
(2,)