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实现忽略大小写对字符串列表排序的方法
Sep 25 Python
Python中使用items()方法返回字典元素对的教程
May 21 Python
Python获取运行目录与当前脚本目录的方法
Jun 01 Python
Python中的zipfile模块使用详解
Jun 25 Python
python实现给数组按片赋值的方法
Jul 28 Python
python的scikit-learn将特征转成one-hot特征的方法
Jul 10 Python
用Python实现读写锁的示例代码
Nov 05 Python
在Python中实现shuffle给列表洗牌
Nov 08 Python
Django框架模型简单介绍与使用分析
Jul 18 Python
使用Python三角函数公式计算三角形的夹角案例
Apr 15 Python
Pycharm 设置默认解释器路径和编码格式的操作
Feb 05 Python
python实现会员管理系统
Mar 18 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从完整文件路径中分离文件目录和文件名的方法
2015/03/13 PHP
PHP+JS实现大规模数据提交的方法
2015/07/02 PHP
浅谈PHP进程管理
2019/03/08 PHP
Thinkphp5框架简单实现钩子(Hook)行为的方法示例
2019/09/03 PHP
转自Jquery官方 jQuery1.1.3发布,速度提升800%,体积保持20K
2007/08/19 Javascript
js模拟类继承小例子
2010/07/17 Javascript
浅析ajax请求json数据并用js解析(示例分析)
2013/07/13 Javascript
JS验证邮箱格式是否正确的代码
2013/12/05 Javascript
悬浮数字的实现案例
2014/02/19 Javascript
JavaScript字符串对象slice方法入门实例(用于字符串截取)
2014/10/16 Javascript
angularjs封装$http为factory的方法
2017/05/18 Javascript
jquery对table做排序操作的实例演示
2017/08/10 jQuery
jQuery实现点击下拉框中的值累加到文本框中的方法示例
2017/10/28 jQuery
JS改变页面颜色源码分享
2018/02/24 Javascript
JS实现图片旋转动画效果封装与使用示例
2018/07/09 Javascript
从vue源码解析Vue.set()和this.$set()
2018/08/30 Javascript
关于vue v-for 循环问题(一行显示四个,每一行的最右边那个计算属性)
2018/09/04 Javascript
vue服务端渲染缓存应用详解
2018/09/12 Javascript
解决angular2在双向数据绑定时[(ngModel)]无法使用的问题
2018/09/13 Javascript
新手快速入门JavaScript装饰者模式与AOP
2019/06/24 Javascript
Python Queue模块详解
2014/11/30 Python
浅析Python3爬虫登录模拟
2018/02/07 Python
基于循环神经网络(RNN)的古诗生成器
2018/03/26 Python
使用python画社交网络图实例代码
2019/07/10 Python
python爬虫中多线程的使用详解
2019/09/23 Python
python中的RSA加密与解密实例解析
2019/11/18 Python
python3的UnicodeDecodeError解决方法
2019/12/20 Python
Python中三维坐标空间绘制的实现
2020/09/22 Python
canvas学习和滤镜实现代码
2018/08/22 HTML / CSS
荷兰手表网站:Watch2Day
2018/07/02 全球购物
Rossignol金鸡美国官网:始于1907年法国百年雪具品牌
2019/03/06 全球购物
广播电视新闻学专业应届生求职信
2013/10/08 职场文书
设备动力科岗位职责范本
2014/02/23 职场文书
2014最新实习证明模板
2014/10/02 职场文书
反腐倡廉观后感
2015/06/08 职场文书
勤俭节约主题班会
2015/08/13 职场文书