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样条插值的实现代码-创新互联
本文地址:http://myzitong.com/article/dhodos.html