使用Python拟合函数曲线需要用到一些第三方库:

  • numpy:科学计算的基础库(例如:矩阵)

  • matplotlib:绘图库

  • scipy:科学计算库

如果没有安装过这些库,需要在命令行中输入下列代码进行安装:

pip install numpy matplotlib scipy

拟合多项式

'''
Author: CloudSir
Date: 2021-08-01 13:40:50
LastEditTime: 2021-08-02 09:41:54
LastEditors: CloudSir
Description: Python拟合多项式
https://github.com/cloudsir
'''
import matplotlib.pyplot as plt
import numpy as npx = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10]
y = [2.83, 9.53, 14.52, 21.57, 38.26, 53.92, 73.15, 101.56, 129.54, 169.75, 207.59]
z1 = np.polyfit(x, y, 3) #用3次多项式拟合,输出系数从高到0
p1 = np.poly1d(z1) #使用次数合成多项式
y_pre = p1(x)plt.plot(x,y,'.')
plt.plot(x,y_pre)
plt.show()

函数说明

np.polyfit(x, y, n)

功能:拟合曲线

参数:

  • x,y:x和y的原始数据
  • n:要拟合的次数

返回值:

  • 一个列表,拟合出的系数,顺序为从高到底

    例: n=3时,会利用 a x 3 + b x 2 + c x + d ax3+bx2+cx+d ax3+bx2+cx+d 拟合函数,并返回拟合出的系数 [a, b, c, d]

np.poly1d(li, r=False)

功能:生成多项式函数

参数:

  • li:

    • 当没有r参数或 r=False 时,传入一个系数列表(次数从高到低),利用该列表生成多项式函数并返回

      import numpy as np

      f = np.poly1d([2, 3, 4])

      “”" f ( x ) = 2 x 2 + 3 x + 4 f(x) = 2x^2+3x+4 f(x)=2x2+3x+4"“”

      print(f(2)) # 18

    • 当参数 r=True 时,传入一个根列表,利用该列表生成多项式函数并返回

      import numpy as np

      f = np.poly1d([2, 3, 4], True)

      “”“ f ( x ) = ( x 2 ) ( x 3 ) ( x 4 ) f(x) = (x-2)*(x-3)*(x-4) f(x)=(x2)(x3)(x4)”“”

      print(f(0)) # -24

返回值:

  • 见上

拟合任意函数

'''
Author: CloudSir
Date: 2021-08-03 15:01:17
LastEditTime: 2021-08-03 15:26:05
LastEditors: CloudSir
Description: Python拟合任意函数
https://github.com/cloudsir
'''
# 引用库函数import numpy as np
import matplotlib.pyplot as plt
from scipy import optimize as opplt.rcParams['font.sans-serif'] = ['SimHei']  # 用来正常显示中文
plt.rcParams['axes.unicode_minus'] = False    # 用来正常显示负号# 需要拟合的函数
def f_1(x, A, B, C):return A * x**2 + B * x + C# 需要拟合的数据组
x_group = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10]
y_group = [2.83, 9.53, 14.52, 21.57, 38.26, 53.92, 73.15, 101.56, 129.54, 169.75, 207.59]# 得到返回的A,B值
A, B, C = op.curve_fit(f_1, x_group, y_group)[0]# 数据点与原先的进行画图比较
plt.scatter(x_group, y_group, marker='o',label='真实值')
x = np.arange(0, 15, 0.01)
y = A * x**2 + B *x + C
plt.plot(x, y,color='red',label='拟合曲线')
plt.legend() # 显示labelplt.show()

函数说明

op.curve_fit(f, x, y)

功能:拟合任意函数

参数:

  • f:要拟合的函数类型

    # 构建一个二次函数
    def f(x, A, B, C):return A * x**2 + B * x + Cop.curve_fit(f, x, y) # 进行拟合
    
  • x, y:x和y的原始数据

返回值:一个元组 (popt,pcov)

  • popt是一个一维数组,表示得到的拟合方程的参数。

  • pcov是一个二维数组,是在popt参数下得到的协方差。

【Python】利用Python拟合函数曲线相关推荐

  1. python使用matplotlib可视化函数曲线、设置y轴为对数坐标(log scale)、默认情况下坐标轴为线性坐标

    python使用matplotlib可视化函数曲线.设置y轴为对数坐标(log scale).默认情况下坐标轴为线性坐标 目录

  2. python 利用python处理excel文件

    python 利用python处理excel文件 1.相关包的安装 2.读取.创建excel文件,遍历内容 3.字体样式调整 4.填充样式调整 5.边框样式调整 6.对齐样式调整 7.excel内添加 ...

  3. python爬虫简单实例-Python 利用Python编写简单网络爬虫实例3

    利用Python编写简单网络爬虫实例3 by:授客 QQ:1033553122 实验环境 python版本:3.3.5(2.7下报错 实验目的 获取目标网站"http://bbs.51tes ...

  4. [Python] 利用 Python 制作“飞花令”小工具

    转载请注明:陈熹 chenx6542@foxmail.com (简书号:半为花间酒) 若公众号内转载请联系公众号:早起Python 介绍:飞花令是古时候人bai们经常玩一种"行酒令" ...

  5. 简单的病毒编程代码Python,利用python编写勒索病毒

    电脑紧急求助 你中毒了,建议你在DOS下杀毒,看看能不能启动,如果还是不能,恐怕你要从装系统了,可以联系我,跟你一起解决这个问题联系qq:201716下面是病毒简介"AdobeR" ...

  6. python中高阶函数改写学生信息管理程序_python利用高阶函数实现剪枝函数

    本文为大家分享了python利用高阶函数实现剪枝函数的具体代码,供大家参考,具体内容如下 案例: 某些时候,我们想要为多个函数,添加某种功能,比如计时统计,记录日志,缓存运算结果等等 需求: 在每个函 ...

  7. 用python进行数据分析举例说明_《利用python进行数据分析》读书笔记 --第一、二章 准备与例子...

    第一章 准备工作 今天开始码这本书--<利用python进行数据分析>.R和python都得会用才行,这是码这本书的原因.首先按照书上说的进行安装,google下载了epd_free-7. ...

  8. 算法动画 - 理解函数曲线

    这篇梳理一些有关算法动画的生成思路. 用算法生成动画,大致可分成两类.一类是基于时间( time-based ),一类是基于帧( frame-based ).其中有何区别,我们先通过两段 Proces ...

  9. python画抛物线_在python中利用最小二乘拟合二次抛物线函数的方法

    1.最小二乘也可以拟合二次函数 我们都知道用最小二乘拟合线性函数没有问题,那么能不能拟合二次函数甚至更高次的函数呢?答案当然是可以的.下面我们就来试试用最小二乘来拟合抛物线形状的的图像. 对于二次函数 ...

  10. python函数拟合不规则曲线_python 对任意数据和曲线进行拟合并求出函数表达式的三种解决方案...

    第一种是进行多项式拟合,数学上可以证明,任意函数都可以表示为多项式形式.具体示例如下. ###拟合年龄 import numpy as np import matplotlib.pyplot as p ...

最新文章

  1. 计算几个变量之间的相关系数,计算协方差矩阵时:TypeError: cannot perform reduce with flexible type
  2. java实现图片对比功能_Java 照片对比功能的实现
  3. 20175204 张湲祯 2018-2019-2《Java程序设计》第九周学习总结
  4. openresty获取nginx 请求方法
  5. Python日志模块的两种用法
  6. linux编译安装Nginx1.6.0+MySQL5.6.19+PHP5.5.14
  7. Git SSH key配置
  8. python面试技巧和注意事项_Python 程序员面试须知须会的5个问题
  9. c#后的完整cookie
  10. 苹果Mac备份应用推荐:Get Backup Pro
  11. tensorflow CNN常用函数汇总
  12. 2022年最新《谷粒学院开发教程》:6 - 整合SpringCloud
  13. Golang语言 零基础入门教程
  14. js神秘的电报密码---哈弗曼编码
  15. 解决h5py.h5 ImportError问题
  16. 前端开发:报错Error in... ”SyntaxError:Unexpected token u in JSON at position 0”…解决方法
  17. Matlab小课堂1
  18. 48万的无人共享车,能让百度破局?还是能“拯救”极狐?
  19. 2023 年度 A 类学科竞赛项目清单
  20. 在Vue中Promise.all的使用

热门文章

  1. android桌面小工具,Android桌面小工具 清新实用的桌面小插件推荐
  2. python关于numpy常用函数思维导图
  3. Python+Tensorflow+CNN实现车牌识别
  4. YOLO+PaddleOCR实现车牌检测识别
  5. VGG16和VGG19网络结构图
  6. 制作种子怎么上传服务器,怎么做种子,又怎么发布呢?
  7. vue php跨域,Vue+php处理跨域问题
  8. PS去除图片和PDF中的水印
  9. UWB定位记录二(DWM1000模组介绍)
  10. python对电影进行预测评分_TMDb电影数据分析 电影评分预测