python使用插值法画出平滑曲线


Posted in Python onDecember 15, 2018

本文实例为大家分享了python使用插值法画出平滑曲线的具体代码,供大家参考,具体内容如下

实现所需的库

numpy、scipy、matplotlib

实现所需的方法

插值

  • nearest:最邻近插值法
  • zero:阶梯插值
  • slinear:线性插值
  • quadratic、cubic:2、3阶B样条曲线插值

拟合和插值的区别

简单来说,插值就是根据原有数据进行填充,最后生成的曲线一定过原有点。

拟合是通过原有数据,调整曲线系数,使得曲线与已知点集的差别(最小二乘)最小,最后生成的曲线不一定经过原有点。

代码实现

# -*- coding: utf-8 -*-

# 调用模块
# 调用数组模块
import numpy as np
# 实现插值的模块
from scipy import interpolate
# 画图的模块
import matplotlib.pyplot as plt
# 生成随机数的模块
import random

# random.randint(0, 10) 生成0-10范围内的一个整型数
# y是一个数组里面有10个随机数,表示y轴的值
y = np.array([random.randint(0, 10) for _ in range(10)])
# x是一个数组,表示x轴的值
x = np.array([num for num in range(10)])

# 插值法之后的x轴值,表示从0到9间距为0.5的18个数
xnew = np.arange(0, 9, 0.5)

"""
kind方法:
nearest、zero、slinear、quadratic、cubic
实现函数func
"""
func = interpolate.interp1d(x, y, kind='cubic')
# 利用xnew和func函数生成ynew,xnew的数量等于ynew数量
ynew = func(xnew)

# 画图部分
# 原图
plt.plot(x, y, 'ro-')
# 拟合之后的平滑曲线图
plt.plot(xnew, ynew)
plt.show()

注意事项/p>

  • x, y为原来的数据(少量)
  • xnew为一个数组,条件:x⊆⊆xnew
  •       如:x的最小值为-2.931,最大值为10.312;则xnew的左边界要小于-2.931,右边界要大于10.312。当然也最好注意一下间距,最好小于x中的精度
  • func为函数,里面的参数x、y、kind,x,y就是原数据的x,y,kind为需要指定的方法
  • ynew需要通过xnew数组和func函数来生成
  • 理论上xnew数组内的值越多,生成的曲线越平滑

以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持三水点靠木。

Python 相关文章推荐
python中的字典详细介绍
Sep 18 Python
详解python时间模块中的datetime模块
Jan 13 Python
python爬取cnvd漏洞库信息的实例
Feb 14 Python
浅谈Python大神都是这样处理XML文件的
May 31 Python
PyCharm中代码字体大小调整方法
Jul 29 Python
python实现桌面气泡提示功能
Jul 29 Python
Python + Flask 实现简单的验证码系统
Oct 01 Python
基于Tensorflow的MNIST手写数字识别分类
Jun 17 Python
关于python的缩进规则的知识点详解
Jun 22 Python
python操作ini类型配置文件的实例教程
Oct 30 Python
pycharm最新激活码有效期至2100年(亲测可用)
Feb 05 Python
python numpy中setdiff1d的用法说明
Apr 22 Python
python用fsolve、leastsq对非线性方程组求解
Dec 15 #Python
python实现一组典型数据格式转换
Dec 15 #Python
python判断计算机是否有网络连接的实例
Dec 15 #Python
Django model反向关联名称的方法
Dec 15 #Python
django orm 通过related_name反向查询的方法
Dec 15 #Python
关于Django ForeignKey 反向查询中filter和_set的效率对比详解
Dec 15 #Python
django 外键model的互相读取方法
Dec 15 #Python
You might like
PHP读取ACCESS数据到MYSQL的代码
2011/05/11 PHP
PHP判断访客是否手机端(移动端浏览器)访问的方法总结【4种方法】
2019/03/27 PHP
鼠标移动到一张图片时变为另一张图片
2006/12/05 Javascript
jQuery实现表单input中提示文字value随鼠标焦点移进移出而显示或隐藏的代码
2010/03/21 Javascript
js优化针对IE6.0起作用(详细整理)
2012/12/25 Javascript
判断滚动条到底部的JS代码
2013/11/04 Javascript
使用jQuery实现验证上传图片的格式与大小
2014/12/03 Javascript
jQuery设置和移除文本框默认值的方法
2015/03/09 Javascript
详解Node.js中的事件机制
2016/09/22 Javascript
详谈jQuery Ajax(load,post,get,ajax)的用法
2017/03/02 Javascript
AngularJs定时器$interval 和 $timeout详解
2017/05/25 Javascript
详解Vue一个案例引发「内容分发slot」的最全总结
2018/12/02 Javascript
javascript面向对象三大特征之继承实例详解
2019/07/24 Javascript
python通过openpyxl生成Excel文件的方法
2015/05/12 Python
pymssql数据库操作MSSQL2005实例分析
2015/05/25 Python
详解Django通用视图中的函数包装
2015/07/21 Python
Python中使用OpenCV库来进行简单的气象学遥感影像计算
2016/02/19 Python
Django xadmin开启搜索功能的实现
2019/11/15 Python
wxPython实现绘图小例子
2019/11/19 Python
python运用pygame库实现双人弹球小游戏
2019/11/25 Python
tensorflow实现打印ckpt模型保存下的变量名称及变量值
2020/01/04 Python
Python字符串格式化常用手段及注意事项
2020/06/17 Python
tensorflow基于CNN实战mnist手写识别(小白必看)
2020/07/20 Python
Python新建项目自动添加介绍和utf-8编码的方法
2020/12/26 Python
Miller Harris官网:英国小众香水品牌
2020/09/24 全球购物
static函数与普通函数有什么区别
2015/12/25 面试题
电子商务专业毕业生工作推荐信
2013/11/17 职场文书
往来会计岗位职责
2013/12/19 职场文书
2014年三八妇女节活动方案
2014/02/28 职场文书
法人委托书范本
2014/04/04 职场文书
小区文明倡议书
2014/05/16 职场文书
预备党员自我批评思想汇报
2014/10/10 职场文书
思想道德自我评价2015
2015/03/09 职场文书
音乐研修感悟
2015/11/18 职场文书
单身狗福利?Python爬取某婚恋网征婚数据
2021/06/03 Python
JavaCV实现照片马赛克效果
2022/01/22 Java/Android