python样条插值的实现代码.pdf

python样条插值的实现代码

本⽂实例为⼤家分享了python样条插值的具体实现代码，供⼤家参考，具体内容如下

importnumpyasnp

importpandasaspd

importmatplotlib.pyplotasplt

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)#定义观测点

#进⾏样条差值

importscipy.interpolateasspi

#进⾏⼀阶样条差值

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与量化投资从基础到实战》的内容练习

以上就是本⽂的全部内容，希望对⼤家的学习有所帮助，也希望⼤家多多⽀持。