DDR爱好者之家 Design By 杰米
本文实例为大家分享了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 与量化投资 从基础到实战》的内容练习
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持。
DDR爱好者之家 Design By 杰米
广告合作:本站广告合作请联系QQ:858582 申请时备注:广告合作(否则不回)
免责声明:本站资源来自互联网收集,仅供用于学习和交流,请遵循相关法律法规,本站一切资源不代表本站立场,如有侵权、后门、不妥请联系本站删除!
免责声明:本站资源来自互联网收集,仅供用于学习和交流,请遵循相关法律法规,本站一切资源不代表本站立场,如有侵权、后门、不妥请联系本站删除!
DDR爱好者之家 Design By 杰米
暂无评论...
更新日志
2024年11月27日
2024年11月27日
- 凤飞飞《我们的主题曲》飞跃制作[正版原抓WAV+CUE]
- 刘嘉亮《亮情歌2》[WAV+CUE][1G]
- 红馆40·谭咏麟《歌者恋歌浓情30年演唱会》3CD[低速原抓WAV+CUE][1.8G]
- 刘纬武《睡眠宝宝竖琴童谣 吉卜力工作室 白噪音安抚》[320K/MP3][193.25MB]
- 【轻音乐】曼托凡尼乐团《精选辑》2CD.1998[FLAC+CUE整轨]
- 邝美云《心中有爱》1989年香港DMIJP版1MTO东芝首版[WAV+CUE]
- 群星《情叹-发烧女声DSD》天籁女声发烧碟[WAV+CUE]
- 刘纬武《睡眠宝宝竖琴童谣 吉卜力工作室 白噪音安抚》[FLAC/分轨][748.03MB]
- 理想混蛋《Origin Sessions》[320K/MP3][37.47MB]
- 公馆青少年《我其实一点都不酷》[320K/MP3][78.78MB]
- 群星《情叹-发烧男声DSD》最值得珍藏的完美男声[WAV+CUE]
- 群星《国韵飘香·贵妃醉酒HQCD黑胶王》2CD[WAV]
- 卫兰《DAUGHTER》【低速原抓WAV+CUE】
- 公馆青少年《我其实一点都不酷》[FLAC/分轨][398.22MB]
- ZWEI《迟暮的花 (Explicit)》[320K/MP3][57.16MB]