python--多项式拟合示例代码(numpy.polyfit)

# 定义x、y散点坐标
import numpy as np
x = [1.414, 1.600, 1.700, 1.720, 1.732, 1.740, 1.750, 1.800, 1.875, 2.000]
y = [929.092, 930.713, 929.903, 930.713, 931.524, 927.472, 925.446, 927.066, 924.230, 928.282]
y= [yi-930.0 for yi in y]
# print(y)
z1 = np.polyfit(x, y, 7) # 用7次多项式拟合,可改变多项式阶数;
print(z1) #显示多项式系数
p1 = np.poly1d(z1) #得到多项式系数,按照阶数从高到低排列
# print(p1)  #显示多项式


from functools import reduce

def polynomial(coefficient, x):
    return reduce(lambda c1, c2: c1 * x + c2, coefficient)


res=polynomial(z1,1.720)
# res=polynomial([1,2,1],3)
# res=polynomial([10841.5747056069, -24457.2572571241, 22489.1422245109,-9130.76364970551,1380.67465252012], 1.732)
print("result")
print(res+930.0)
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值