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 相关文章推荐
wxpython中利用线程防止假死的实现方法
Aug 11 Python
Python import自定义模块方法
Feb 12 Python
Python变量类型知识点总结
Feb 18 Python
Python3.5面向对象程序设计之类的继承和多态详解
Apr 24 Python
pytorch torchvision.ImageFolder的用法介绍
Feb 20 Python
python实现opencv+scoket网络实时图传
Mar 20 Python
Python插件机制实现详解
May 04 Python
Python常见反爬虫机制解决方案
Jun 01 Python
Python虚拟环境的创建和使用详解
Sep 07 Python
Python3.9最新版下载与安装图文教程详解(Windows系统为例)
Nov 28 Python
python自动化办公操作PPT的实现
Feb 05 Python
python周期任务调度工具Schedule使用详解
Nov 23 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类Class的概念
2012/06/14 PHP
php写的带缓存数据功能的mysqli类
2012/09/06 PHP
php根据年月获取季度的方法
2014/03/31 PHP
PHP+Mysql实现多关键字与多字段生成SQL语句的函数
2014/11/05 PHP
PHP生成word文档的三种实现方式
2016/11/14 PHP
php pdo连接数据库操作示例
2019/11/18 PHP
通过PHP实现用户注册后邮箱验证激活
2020/11/10 PHP
firefox下jquery ajax返回object XMLDocument处理方法
2014/01/26 Javascript
JS通过分析userAgent属性来判断浏览器的类型及版本
2014/03/28 Javascript
jQuery中scrollLeft()方法用法实例
2015/01/16 Javascript
JS+CSS实现淡入式焦点图片幻灯切换效果的方法
2015/02/26 Javascript
AngularJS定时器的使用与移除操作方法【interval与timeout】
2016/12/14 Javascript
Angularjs 实现移动端在线测评效果(推荐)
2017/04/05 Javascript
浅谈js-FCC算法Friendly Date Ranges(详解)
2017/04/10 Javascript
JS实现问卷星自动填问卷脚本并在两秒自动提交功能
2020/06/17 Javascript
基于 flexible 的 Vue 组件:Toast -- 显示框效果
2017/12/26 Javascript
关于vue单文件中引用路径的处理方法
2018/01/08 Javascript
js实现鼠标单击Tab表单切换效果
2018/05/16 Javascript
在vue中使用公共过滤器filter的方法
2018/06/26 Javascript
webpack 如何解析代码模块路径的实现
2019/09/04 Javascript
Vue多选列表组件深入详解
2021/03/02 Vue.js
python 运算符 供重载参考
2009/06/11 Python
python 爬虫出现403禁止访问错误详解
2017/03/11 Python
Python基于回溯法子集树模板解决找零问题示例
2017/09/11 Python
修改python plot折线图的坐标轴刻度方法
2018/12/13 Python
Python实现计算字符串中出现次数最多的字符示例
2019/01/21 Python
python Tkinter的图片刷新实例
2019/06/14 Python
django 数据库 get_or_create函数返回值是tuple的问题
2020/05/15 Python
python 引用传递和值传递详解(实参,形参)
2020/06/05 Python
keras model.fit 解决validation_spilt=num 的问题
2020/06/19 Python
通过css3的filter滤镜改变png图片的颜色的示例代码
2020/05/06 HTML / CSS
幼儿园端午节活动方案
2014/08/25 职场文书
纪检干部个人对照检查材料
2014/09/23 职场文书
社区义诊通知
2015/04/24 职场文书
2016春节放假通知范文
2015/08/18 职场文书
电脑关机速度很慢怎么办 提升电脑关机速度设置教程
2022/04/08 数码科技