muchener's blogs

数模笔记:多项式拟合(1)

字数统计: 226阅读时长: 1 min
2017/01/22 Share

多项式拟合

多项式拟合:poly= np.polyfit(x,a,n),拟合点集a得到n级多项式,其中x为横轴长度,返回多项式的系数

多项式求导函数:np.polyder(poly),返回导函数的系数得到多项式的n阶导函数:多项式.deriv(m = n)

多项式求根:np.roots(poly)

多项式在某点上的值:np.polyval(poly,x[n]),返回poly多项式在横轴点上x[n]上的值

两个多项式做差运算: np.polysub(a,b)

多项式方程:np.poly1d(ploy)

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
import numpy as np
import pylab as pl
mass = [50*i for i in range(1,12)]
length = [1.000,1.875,2.750,3.250,4.375,4.875,5.675,6.500,7.250,8.000,8.750]

pl.plot(mass, length,'o')
c = np.polyfit(mass, length,1)#按一次多项式拟合
x_new = np.linspace(-100, 900, 100)##取样
y_new = np.polyval(c,x_new)## 求y

pl.plot(x_new,y_new)
pl.xlim(0.0,600)
pl.ylim(0,9)
pl.show()

结果

CATALOG
  1. 1. 多项式拟合