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常见文件操作的函数示例代码
Nov 15 Python
用Python进行一些简单的自然语言处理的教程
Mar 31 Python
Python中shutil模块的常用文件操作函数用法示例
Jul 05 Python
Python实现的堆排序算法原理与用法实例分析
Nov 22 Python
python3+PyQt5实现拖放功能
Apr 24 Python
Python遍历文件夹 处理json文件的方法
Jan 22 Python
python实现LBP方法提取图像纹理特征实现分类的步骤
Jul 11 Python
Python实现的远程文件自动打包并下载功能示例
Jul 12 Python
Django ORM 自定义 char 类型字段解析
Aug 09 Python
python树的同构学习笔记
Sep 14 Python
Python3运算符常见用法分析
Feb 14 Python
python爬虫爬取图片的简单代码
Jan 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 解决substr()截取中文字符乱码问题
2016/07/18 PHP
Laravel的throttle中间件失效问题解决方法
2016/10/09 PHP
php计算多个集合的笛卡尔积实例详解
2017/02/16 PHP
详解php实现页面静态化原理
2017/06/21 PHP
PHP parse_ini_file函数的应用与扩展操作示例
2019/01/07 PHP
Laravel框架源码解析之反射的使用详解
2020/05/14 PHP
用javascript实现无刷新更新数据的详细步骤 asp
2006/12/26 Javascript
二级域名或跨域共享Cookies的实现方法
2008/08/07 Javascript
面向对象的编程思想在javascript中的运用上部
2009/11/20 Javascript
常用javascript表单验证汇总
2020/07/20 Javascript
基于BootStrap的Metronic框架实现页面链接收藏夹功能按钮移动收藏记录(使用Sortable进行拖动排序)
2016/08/29 Javascript
jquery replace方法去空格
2017/05/08 jQuery
详解使用create-react-app添加css modules、sasss和antd
2018/07/31 Javascript
浅谈关于vue中scss公用的解决方案
2019/12/02 Javascript
js实现树形数据转成扁平数据的方法示例
2020/02/27 Javascript
JS禁用右键、禁用Ctrl+u、禁用Ctrl+s、禁用F12的实现代码
2020/12/01 Javascript
jQuery实现动态操作table行
2020/11/23 jQuery
python正则表达式判断字符串是否是全部小写示例
2013/12/25 Python
深入解析Python的Tornado框架中内置的模板引擎
2016/07/11 Python
Python爬虫之正则表达式基本用法实例分析
2018/08/08 Python
Windows 64位下python3安装nltk模块
2018/09/19 Python
Python3.5内置模块之os模块、sys模块、shutil模块用法实例分析
2019/04/27 Python
anaconda如何查看并管理python环境
2019/07/05 Python
python3-flask-3将信息写入日志的实操方法
2019/11/12 Python
使用Pyhton集合set()实现成果查漏的例子
2019/11/24 Python
python实现将中文日期转换为数字日期
2020/07/14 Python
python开发一个解析protobuf文件的简单编译器
2020/11/17 Python
Silk’n激光脱毛器官网:silkn.com
2016/10/06 全球购物
大学同学聚会邀请函
2014/01/19 职场文书
大学生求职工作的自我评价
2014/02/13 职场文书
纪律教育学习月活动总结
2014/08/27 职场文书
违纪开除通知书
2015/04/25 职场文书
2015年秋季运动会广播稿
2015/08/19 职场文书
MySQL数据库如何给表设置约束详解
2022/03/13 MySQL
详解Vue3使用axios的配置教程
2022/04/29 Vue.js
Nginx使用ngx_http_upstream_module实现负载均衡功能示例
2022/08/05 Servers