数模笔记：多项式拟合（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)

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()

结果

Next Post

数模笔记：差分方程
Previous Post

prompt(1)-to-win

CATALOG

1. 多项式拟合

