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 相关文章推荐
Python中列表、字典、元组数据结构的简单学习笔记
Mar 20 Python
python socket多线程通讯实例分析(聊天室)
Apr 06 Python
利用python爬取散文网的文章实例教程
Jun 18 Python
对python程序内存泄漏调试的记录
Jun 11 Python
利用python脚本如何简化jar操作命令
Feb 24 Python
Python基础之函数的定义与使用示例
Mar 23 Python
PyCharm 创建指定版本的 Django(超详图解教程)
Jun 18 Python
python爬取百度贴吧前1000页内容(requests库面向对象思想实现)
Aug 10 Python
在Mac中PyCharm配置python Anaconda环境过程图解
Mar 11 Python
Python+PyQt5+MySQL实现天气管理系统
Jun 16 Python
python连接手机自动搜集蚂蚁森林能量的实现代码
Feb 24 Python
Jupyter notebook 输出部分显示不全的解决方案
Apr 24 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合并discuz用户脚本的方法
2015/08/04 PHP
php正则表达式基本知识与应用详解【经典教程】
2017/04/17 PHP
php使用array_chunk函数将一个数组分割成多个数组
2018/12/05 PHP
javascript 兼容鼠标滚轮事件
2009/04/07 Javascript
在JS中最常看到切最容易迷惑的语法(转)
2010/10/29 Javascript
JS隐藏参数post传值实例
2013/04/18 Javascript
js中的this关键字详解
2013/09/25 Javascript
jquery实现滑动特效代码
2015/08/10 Javascript
纯HTML5制作围住神经猫游戏-附源码下载
2015/08/23 Javascript
在easyUI开发中,出现jquery.easyui.min.js函数库问题的解决办法
2015/09/11 Javascript
jQuery soColorPacker 网页拾色器
2016/06/22 Javascript
微信小程序 安全包括(框架、功能模块、账户使用)详解
2017/01/16 Javascript
vue登录路由验证的实现
2017/12/13 Javascript
Vue实现一个无限加载列表功能
2018/11/13 Javascript
layui 上传图片 返回图片地址的方法
2019/09/26 Javascript
[02:36]DOTA2英雄基础教程 一击致命幻影刺客
2013/12/06 DOTA
Python实现字符串匹配算法代码示例
2017/12/05 Python
50行Python代码获取高考志愿信息的实现方法
2019/07/23 Python
keras多显卡训练方式
2020/06/10 Python
python re模块常见用法例举
2021/03/01 Python
使用css3背景渐变中的透明度来设置不同颜色的背景渐变
2014/03/31 HTML / CSS
CSS3实现可爱的小黄人动画
2016/07/11 HTML / CSS
html5实现完美兼容各大浏览器的播放器
2014/12/26 HTML / CSS
计算机专业推荐信范文
2013/11/20 职场文书
化工专业个人的求职信范文
2013/11/28 职场文书
自我鉴定写作要点
2014/01/17 职场文书
《猫》教学反思
2014/02/26 职场文书
学徒工职责
2014/03/06 职场文书
高校教师岗位职责
2014/03/18 职场文书
党员领导干部承诺书
2014/05/28 职场文书
大学生求职信
2014/06/17 职场文书
教师党员个人自我剖析材料
2014/09/29 职场文书
2014年物业管理工作总结
2014/11/21 职场文书
2014年财政局工作总结
2014/12/09 职场文书
MySQL Shell的介绍以及安装
2021/04/24 MySQL
Hive常用日期格式转换语法
2022/06/25 数据库