时间:2021-05-22
本文实例为大家分享了python样条插值的具体实现代码,供大家参考,具体内容如下
import numpy as npimport pandas as pdimport 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 与量化投资 从基础到实战》的内容练习
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持。
声明:本页内容来源网络,仅供用户参考;我单位不保证亦不表示资料全面及准确无误,也不保证亦不表示这些资料为最新信息,如因任何原因,本网内容或者用户因倚赖本网内容造成任何损失或损害,我单位将不会负任何法律责任。如涉及版权问题,请提交至online#300.cn邮箱联系删除。
本文实例为大家分享了python实现三次样条插值的具体代码,供大家参考,具体内容如下函数:算法分析三次样条插值。就是在分段插值的一种情况。要求:在每个分段区间上
本文实例讲述了Python使用scipy模块实现一维卷积运算。分享给大家供大家参考,具体如下:一介绍signal模块包含大量滤波函数、B样条插值算法等等。下面的
一维插值插值不同于拟合。插值函数经过样本点,拟合函数一般基于最小二乘法尽量靠近所有样本点穿过。常见插值方法有拉格朗日插值法、分段插值法、样条插值法。拉格朗日插值
一维插值插值不同于拟合。插值函数经过样本点,拟合函数一般基于最小二乘法尽量靠近所有样本点穿过。常见插值方法有拉格朗日插值法、分段插值法、样条插值法。拉格朗日插值
本文实例为大家分享了python使用插值法画出平滑曲线的具体代码,供大家参考,具体内容如下实现所需的库numpy、scipy、matplotlib实现所需的方法