python样条插值的实现代码


Posted in Python onDecember 17, 2018

本文实例为大家分享了python样条插值的具体实现代码,供大家参考,具体内容如下

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
 
plt.rcParams['font.sans-serif']=['SimHei'] #用来正常显示中文标签
plt.rcParams['axes.unicode_minus']=False #用来正常显示负号
 
#导入数据
data1=pd.read_csv('data1.csv',encoding='gbk')
 
#数据准备
X=data1.index #定义数据点
Y=data1['沪深300'].values #定义数据点
x=np.arange(0,len(data1),0.15) #定义观测点
 
#进行样条差值
import scipy.interpolate as spi
 
#进行一阶样条差值
ipo1=spi.splrep(X,Y,k=1) #源数据点导入,生成参数
iy1=spi.splev(x,ipo1) #根据观测点和样条参数,生成插值
 
#进行三次样条拟合
ipo3=spi.splrep(X,Y,k=3) #源数据点导入,生成参数
iy3=spi.splev(x,ipo3) #根据观测点和样条参数,生成插值
 
 
##作图
fig,(ax1,ax2)=plt.subplots(2,1,figsize=(10,12))
ax1.plot(X,Y,label='沪深300')
ax1.plot(x,iy1,'r.',label='插值点')
ax1.set_ylim(Y.min()-10,Y.max()+10)
ax1.set_ylabel('指数')
ax1.set_title('线性插值')
ax1.legend()
ax2.plot(X,Y,label='沪深300')
ax2.plot(x,iy3,'b.',label='插值点')
ax2.set_ylim(Y.min()-10,Y.max()+10)
ax2.set_ylabel('指数')
ax2.set_title('三次样条插值')
ax2.legend()

python样条插值的实现代码

基于《Python 与量化投资 从基础到实战》的内容练习

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

Python 相关文章推荐
Linux中Python 环境软件包安装步骤
Mar 31 Python
Python利用multiprocessing实现最简单的分布式作业调度系统实例
Nov 14 Python
利用Python找出序列中出现最多的元素示例代码
Dec 08 Python
wxPython的安装与使用教程
Aug 31 Python
python+openCV利用摄像头实现人员活动检测
Jun 22 Python
Django 自动生成api接口文档教程
Nov 19 Python
tensorflow通过模型文件,使用tensorboard查看其模型图Graph方式
Jan 23 Python
对pytorch的函数中的group参数的作用介绍
Feb 18 Python
python字典的值可以修改吗
Jun 29 Python
通过代码实例了解Python sys模块
Sep 14 Python
详解python百行有效代码实现汉诺塔小游戏(简约版)
Oct 30 Python
python状态机transitions库详解
Jun 02 Python
解决python3 HTMLTestRunner测试报告中文乱码的问题
Dec 17 #Python
python3爬虫获取html内容及各属性值的方法
Dec 17 #Python
解决Python中list里的中文输出到html模板里的问题
Dec 17 #Python
python实现三次样条插值
Dec 17 #Python
Python命名空间的本质和加载顺序
Dec 17 #Python
对python的unittest架构公共参数token提取方法详解
Dec 17 #Python
Python单元测试unittest的具体使用示例
Dec 17 #Python
You might like
php db类库进行数据库操作
2009/03/19 PHP
PHP实现中文圆形印章特效
2015/06/19 PHP
XHProf报告字段含义的解析
2016/05/17 PHP
php 自定义错误日志实例详解
2016/11/12 PHP
php 读取文件夹下所有图片、文件的实例
2018/10/17 PHP
thinkphp5+layui实现的分页样式示例
2019/10/08 PHP
为Yahoo! UI Extensions Grid增加内置的可编辑器
2007/03/10 Javascript
jQuery 学习入门篇附实例代码
2010/03/16 Javascript
jQuery 动画弹出窗体支持多种展现方式
2010/04/29 Javascript
Jquery公告滚动+AJAX后台得到数据
2011/04/14 Javascript
正负小数点后两位浮点数实现原理及代码
2013/09/06 Javascript
js函数与php函数的区别实例浅析
2015/01/12 Javascript
JavaScript错误处理
2015/02/03 Javascript
JavaScript保存并运算页面中数字类型变量的写法
2015/07/06 Javascript
JS基于FileSystemObject创建一个指定路径的TXT文本文件
2015/08/05 Javascript
JS去除空格和换行的正则表达式(推荐)
2016/06/14 Javascript
js原生map实现的方法总结
2020/01/19 Javascript
JavaScript实现省份城市的三级联动
2020/02/11 Javascript
jquery检测上传文件大小示例
2020/04/26 jQuery
[03:12]TI9战队档案 - Virtus Pro
2019/08/20 DOTA
[02:38]DOTA2 夜魇暗潮2020活动介绍官方视频
2020/11/04 DOTA
python抓取网页中的图片示例
2014/02/28 Python
推荐10款最受Python开发者欢迎的Python IDE
2018/09/16 Python
keras多显卡训练方式
2020/06/10 Python
深入了解Python装饰器的高级用法
2020/08/13 Python
CSS3 实现穿梭星空动画
2020/11/13 HTML / CSS
英语系本科生个人求职信
2013/09/21 职场文书
大学生预备党员自我评价分享
2013/11/16 职场文书
计算机科学系职业生涯规划书
2014/03/08 职场文书
生日宴会主持词
2014/03/20 职场文书
英文推荐信格式范文
2014/05/09 职场文书
大学生迟到检讨书500字
2014/10/17 职场文书
2015年元旦晚会活动总结(学生会)
2014/11/28 职场文书
运动会5000米加油稿
2015/07/21 职场文书
导游词之峨眉山
2019/12/16 职场文书
Python制作春联的示例代码
2022/01/22 Python