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()