溫馨提示×
您好，登錄后才能下訂單哦！
點擊 登錄注冊 即表示同意《億速云用戶服務條款》
您好，登錄后才能下訂單哦！
本文實例為大家分享了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 與量化投資 從基礎到實戰》的內容練習
以上就是本文的全部內容，希望對大家的學習有所幫助，也希望大家多多支持億速云。
免責聲明：本站發布的內容（圖片、視頻和文字）以原創、轉載和分享為主，文章觀點不代表本網站立場，如果涉及侵權請聯系站長郵箱：is@yisu.com進行舉報，并提供相關證據，一經查實，將立刻刪除涉嫌侵權內容。
