python样条插值的实现代码
更新时间：2018年12月17日 14:49:28 作者：lucindawuyi
这篇文章主要为大家详细介绍了python样条插值的实现代码，具有一定的参考价值，感兴趣的小伙伴们可以参考一下
本文实例为大家分享了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使用scrapy采集时伪装成HTTP/1.1的方法
这篇文章主要介绍了Python使用scrapy采集时伪装成HTTP/1.1的方法,实例分析了scrapy采集的使用技巧,非常具有实用价值,需要的朋友可以参考下2015-04-04
在本篇文章里小编给大家分享了关于Python中添加自定义模块的实例方法，需要的朋友们可以学习下。2020-06-06
这篇文章主要介绍了pyinstaller打包找不到文件的问题解决，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧2020-04-04
这篇文章主要介绍了python多环境切换及pyenv使用过程详解,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下2019-09-09
今天小编就为大家分享一篇python 对给定可迭代集合统计出现频率,并排序的方法，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2018-10-10
重试机制在编程中是比较常见的场景,主要被用于处理那些可能由于临时性故障或网络波动等原因而失败的操作,下面我们就来看看如何使用装饰器来实现重试机制吧2024-03-03
今天小编就为大家分享一篇基于python代码实现简易滤除数字的方法，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2018-07-07
在PDF中添加条码是一个常见需求,特别是在需要自动化处理、跟踪或检索PDF文件时,作为一种机器可读的标识符,PDF中的条码可以包含各种类型的信息,本文介绍了如何使用Python在PDF中添加条形码或二维码,需要的朋友可以参考下2024-08-08
这篇文章介绍了Python标准库之日期、时间和日历模块，文中通过示例代码介绍的非常详细。对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下2022-05-05
使用django-guardian实现django-admin的行级权限控制的方法
这篇文章主要介绍了使用django-guardian实现django-admin的行级权限控制的方法，小编觉得挺不错的，现在分享给大家，也给大家做个参考。一起跟随小编过来看看吧2018-10-10
最新评论