Python实现线性插值和三次样条插值的示例代码


Posted in Python onNovember 13, 2019

(1)、函数

y = sin(x)

(2)、数据准备

#数据准备
X=np.arange(-np.pi,np.pi,1) #定义样本点X,从-pi到pi每次间隔1
Y= np.sin(X)#定义样本点Y,形成sin函数
new_x=np.arange(-np.pi,np.pi,0.1) #定义差值点

(3)、样条插值

#进行样条差值
import scipy.interpolate as spi
 
#进行一阶样条插值
ipo1=spi.splrep(X,Y,k=1) #样本点导入,生成参数
iy1=spi.splev(new_x,ipo1) #根据观测点和样条参数,生成插值
 
#进行三次样条拟合
ipo3=spi.splrep(X,Y,k=3) #样本点导入,生成参数
iy3=spi.splev(new_x,ipo3) #根据观测点和样条参数,生成插值

(4)、绘图

##作图
fig,(ax1,ax2)=plt.subplots(2,1,figsize=(10,12))

ax1.plot(X,Y,'o',label='样本点')
ax1.plot(new_x,iy1,label='插值点')
ax1.set_ylim(Y.min()-1,Y.max()+1)
ax1.set_ylabel('指数')
ax1.set_title('线性插值')
ax1.legend()

ax2.plot(X,Y,'o',label='样本点')
ax2.plot(new_x,iy3,label='插值点')
ax2.set_ylim(Y.min()-1,Y.max()+1)
ax2.set_ylabel('指数')
ax2.set_title('三次样条插值')
ax2.legend()

(5)、综合代码

import numpy as np
import matplotlib.pyplot as plt
 
plt.rcParams['font.sans-serif']=['SimHei'] #用来正常显示中文标签
plt.rcParams['axes.unicode_minus']=False #用来正常显示负号
 
#数据准备
X=np.arange(-np.pi,np.pi,1) #定义样本点X,从-pi到pi每次间隔1
Y= np.sin(X)#定义样本点Y,形成sin函数
new_x=np.arange(-np.pi,np.pi,0.1) #定义差值点
 
#进行样条差值
import scipy.interpolate as spi
 
#进行一阶样条插值
ipo1=spi.splrep(X,Y,k=1) #样本点导入,生成参数
iy1=spi.splev(new_x,ipo1) #根据观测点和样条参数,生成插值
 
#进行三次样条拟合
ipo3=spi.splrep(X,Y,k=3) #样本点导入,生成参数
iy3=spi.splev(new_x,ipo3) #根据观测点和样条参数,生成插值

 
 
##作图
fig,(ax1,ax2)=plt.subplots(2,1,figsize=(10,12))

ax1.plot(X,Y,'o',label='样本点')
ax1.plot(new_x,iy1,label='插值点')
ax1.set_ylim(Y.min()-1,Y.max()+1)
ax1.set_ylabel('指数')
ax1.set_title('线性插值')
ax1.legend()

ax2.plot(X,Y,'o',label='样本点')
ax2.plot(new_x,iy3,label='插值点')
ax2.set_ylim(Y.min()-1,Y.max()+1)
ax2.set_ylabel('指数')
ax2.set_title('三次样条插值')
ax2.legend()

(6)、结果展示

Python实现线性插值和三次样条插值的示例代码

Python实现线性插值和三次样条插值的示例代码

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

Python 相关文章推荐
Python简单进程锁代码实例
Apr 27 Python
学习python 之编写简单乘法运算题
Feb 27 Python
python生成词云的实现方法(推荐)
Jun 13 Python
Python爬虫DOTA排行榜爬取实例(分享)
Jun 13 Python
Python基于matplotlib绘制栈式直方图的方法示例
Aug 09 Python
简单学习Python多进程Multiprocessing
Aug 29 Python
在python2.7中用numpy.reshape 对图像进行切割的方法
Dec 05 Python
对python:threading.Thread类的使用方法详解
Jan 31 Python
Pandas之MultiIndex对象的示例详解
Jun 25 Python
OpenCV+Python--RGB转HSI的实现
Nov 27 Python
如何理解Python中的变量
Jun 01 Python
基于python SMTP实现自动发送邮件教程解析
Jun 02 Python
Python坐标线性插值应用实现
Nov 13 #Python
python如果快速判断数字奇数偶数
Nov 13 #Python
Python 异步协程函数原理及实例详解
Nov 13 #Python
python文字转语音实现过程解析
Nov 12 #Python
python文字转语音的实例代码分析
Nov 12 #Python
Python上下文管理器全实例详解
Nov 12 #Python
python3-flask-3将信息写入日志的实操方法
Nov 12 #Python
You might like
让回调函数 showResponse 也带上参数的代码
2007/08/13 Javascript
动态表格Table类的实现
2009/08/26 Javascript
JavaScript 定义function的三种方式小结
2009/10/16 Javascript
javascript中类的定义及其方式(《javascript高级程序设计》学习笔记)
2011/07/04 Javascript
jQuery动态添加 input type=file的实现代码
2012/06/14 Javascript
JS阻止冒泡事件以及默认事件发生的简单方法
2014/01/17 Javascript
javascript单引号和双引号的区别和处理
2014/05/14 Javascript
原生javascript实现图片滚动、延时加载功能
2015/01/12 Javascript
深入理解JavaScript系列(31):设计模式之代理模式详解
2015/03/03 Javascript
JavaScript获取表单内所有元素值的方法
2015/04/02 Javascript
jquery马赛克拼接翻转效果代码分享
2015/08/24 Javascript
jQuery随手笔记之常用的jQuery操作DOM事件
2015/11/29 Javascript
JavaScript数组操作函数汇总
2016/08/05 Javascript
node.js利用redis数据库缓存数据的方法
2017/03/01 Javascript
关闭Vue计算属性自带的缓存功能方法
2018/03/02 Javascript
Vue CLI3.0中使用jQuery和Bootstrap的方法
2019/02/28 jQuery
小程序跳转到的H5页面再跳转回跳小程序的方法
2020/03/06 Javascript
详解vue3.0 diff算法的使用(超详细)
2020/07/01 Javascript
从源码角度来回答keep-alive组件的缓存原理
2021/01/18 Javascript
[01:58]DOTA2上海特级锦标赛现场采访:RTZ这个ID到底好不好
2016/03/25 DOTA
[02:04]2018DOTA2亚洲邀请赛Secret赛前采访
2018/04/03 DOTA
Python下的twisted框架入门指引
2015/04/15 Python
python: line=f.readlines()消除line中\n的方法
2018/03/19 Python
对Python中内置异常层次结构详解
2018/10/18 Python
python命令行工具Click快速掌握
2019/07/04 Python
python使用梯度下降和牛顿法寻找Rosenbrock函数最小值实例
2020/04/02 Python
keras分类模型中的输入数据与标签的维度实例
2020/07/03 Python
CSS3截取字符串实例代码【推荐】
2018/06/07 HTML / CSS
CSS3实现伪类hover离开时平滑过渡效果示例
2017/08/10 HTML / CSS
英国的领先快速时尚零售商:In The Style
2019/03/25 全球购物
建龙钢铁面试总结
2014/04/15 面试题
医院后勤自我鉴定
2013/10/13 职场文书
大学生优秀的自我评价分享
2013/10/22 职场文书
售房协议书
2014/08/19 职场文书
党员教师群众路线个人整改措施
2014/10/28 职场文书
2016暑期校本培训心得体会
2016/01/08 职场文书