Numpy 多项式函数、求导

多项式函数是变量的整数次冥与系数的乘积之和,形如
f ( x ) = a n x n + a n 1 x n 1 + . . . . . . + a 2 x 2 + a 1 x 1 f(x)=a_{n}x^{n}+a_{n-1}x^{n-1}+......+a_{2}x^{2}+a_{1}x^{1}

例: f ( x ) = 3 2 x 2 + 7 1 x 1 + 2 f(x)=3_{2}x^{2}+7_{1}x^{1}+2

# python
a = np.array([3,7,2]) 
p = np.poly1d(a)
print(p)

numpy.poly1d(x)函数指只有一个变量的的多项式,x的维数即为多项式的最高次数

①加减

p = p + np.poly1d([1,2])
# p + [1,2]  亦可

输出: 3 2 x 2 + 7 1 x 1 + 2 1 x 1 2 3_{2}x^{2}+7_{1}x^{1}+2-{1}x^{1}-2


②乘除

a = np.array([3,7,2]) 
p = np.poly1d(a)
p/[1,2]
>>>(poly1d([ 3.,  4.]), poly1d([-2.])) #返回2个多项式除法的结果,分别为商式和余式  

输出:
3 x + 4
2


③微分与积分
#p指对象多项式
p.deriv()     #微分
f ( x ) = 3 2 x 2 + 7 1 x 1 + 2 f(x)=3_{2}x^{2}+7_{1}x^{1}+2

a = np.array([3,7,2]) 
p = np.poly1d(a)
print(p.deriv())

输出:6x + 7


p.integ()     #积分        integ(m=1,k= 0)   m是积几次分,k是常系数是多少


④开根

np.roots§

更多详见:
nupy官方文档

猜你喜欢

转载自blog.csdn.net/HaruStone/article/details/88741583
今日推荐