python多项式拟合：np.polyfit和np.polyld
python数据拟合主要可采⽤numpy库，库的安装可直接⽤pip install numpy等。
1. 原始数据：假如要拟合的数据yyy来⾃sin函数，np.sin
import numpy as np
import matplotlib.pyplot as plt
xxx = np.arange(0, 1000)  # x值，此时表⽰弧度
yyy = np.sin(xxx*np.pi/180)  #函数值，转化成度
2. 测试不同阶的多项式，例如7阶多项式拟合，使⽤np.polyfit拟合，np.polyld得到多项式系数
z1 = np.polyfit(xxx, yyy, 7) # ⽤7次多项式拟合，可改变多项式阶数；
p1 = np.poly1d(z1) #得到多项式系数，按照阶数从⾼到低排列
print(p1)  #显⽰多项式
3. 求对应xxx的各项拟合函数值
yvals=p1(xxx) # 可直接使⽤yvals=np.polyval(z1,xxx)
4. 绘图如下
plt.plot(xxx, yyy, '*',label='original values')
plt.plot(xxx, yvals, 'r',label='polyfit values')
numpy库需要安装吗plt.xlabel('x axis')
plt.ylabel('y axis')
plt.legend(loc=4) # 指定legend在图中的位置，类似象限的位置
plt.title('polyfitting')
plt.show()
5. np.polyfit函数：采⽤的是最⼩⼆次拟合，numpy.polyfit(x, y, deg, rcond=None, full=False, w=None, cov=False)，前三个参数是必须的
6. np.polyld函数：得到多项式系数，主要有三个参数
A one-dimensional polynomial class.
A convenience class, used to encapsulate "natural" operations on
polynomials so that said operations may take on their customary
form in code (see Examples).
Parameters
----------
c_or_r : array_like
The polynomial's coefficients, in decreasing powers, or if
the value of the second parameter is True, the polynomial's
roots (values where the polynomial evaluates to 0).  For example,
``poly1d([1, 2, 3])`` returns an object that represents
:math:`x^2 + 2x + 3`, whereas ``poly1d([1, 2, 3], True)`` returns
one that represents :math:`(x-1)(x-2)(x-3) = x^3 - 6x^2 + 11x -6`.
r : bool, optional
If True, `c_or_r` specifies the polynomial's roots; the default
is False.
variable : str, optional
Changes the variable used when printing `p` from `x` to `variable`
(see Examples).
参数1表⽰：在没有参数2（也就是参数2默认False时），参数1是⼀个数组形式，且表⽰从⾼到低的多项式系数项，例如参数1为[4,5,6]表⽰：
参数2表⽰：为True时，表⽰将参数1中的参数作为根来形成多项式，即参数1为[4,5,6]时表⽰：(x-4)(x-5)(x-6)=0，也就是：
参数3表⽰：换参数标识，⽤惯了x，可以⽤ t，s之类的
⽤法：
1. 直接进⾏运算，例如多项式的平⽅，分别得到
xx=np.poly1d([1,2,3])
print(xx)
yy=xx**2  #求平⽅，或者⽤ xx * xx
print(yy)
2. 求值：
yy(1) = 36
3. 求根：即等式为0时的未知数值yy.r
4. 得到系数形成数组：
yy.c 为：array([ 1,  4, 10, 12,  9]) 5. 返回最⾼次幂数：
6. 返回系数：
yy[0] —— 表⽰幂为0的系数
yy[1] —— 表⽰幂为1的系数
参考：

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。