python樣條插值的實(shí)現(xiàn)代碼
更新時(shí)間：2018年12月17日 14:49:28 作者：lucindawuyi
這篇文章主要為大家詳細(xì)介紹了python樣條插值的實(shí)現(xiàn)代碼，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下
本文實(shí)例為大家分享了python樣條插值的具體實(shí)現(xiàn)代碼，供大家參考，具體內(nèi)容如下
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
plt.rcParams['font.sans-serif']=['SimHei'] #用來正常顯示中文標(biāo)簽
plt.rcParams['axes.unicode_minus']=False #用來正常顯示負(fù)號(hào)
#導(dǎo)入數(shù)據(jù)
data1=pd.read_csv('data1.csv',encoding='gbk')
#數(shù)據(jù)準(zhǔn)備
X=data1.index #定義數(shù)據(jù)點(diǎn)
Y=data1['滬深300'].values #定義數(shù)據(jù)點(diǎn)
x=np.arange(0,len(data1),0.15) #定義觀測(cè)點(diǎn)
#進(jìn)行樣條差值
import scipy.interpolate as spi
#進(jìn)行一階樣條差值
ipo1=spi.splrep(X,Y,k=1) #源數(shù)據(jù)點(diǎn)導(dǎo)入，生成參數(shù)
iy1=spi.splev(x,ipo1) #根據(jù)觀測(cè)點(diǎn)和樣條參數(shù)，生成插值
#進(jìn)行三次樣條擬合
ipo3=spi.splrep(X,Y,k=3) #源數(shù)據(jù)點(diǎn)導(dǎo)入，生成參數(shù)
iy3=spi.splev(x,ipo3) #根據(jù)觀測(cè)點(diǎn)和樣條參數(shù)，生成插值
##作圖
fig,(ax1,ax2)=plt.subplots(2,1,figsize=(10,12))
ax1.plot(X,Y,label='滬深300')
ax1.plot(x,iy1,'r.',label='插值點(diǎn)')
ax1.set_ylim(Y.min()-10,Y.max()+10)
ax1.set_ylabel('指數(shù)')
ax1.set_title('線性插值')
ax1.legend()
ax2.plot(X,Y,label='滬深300')
ax2.plot(x,iy3,'b.',label='插值點(diǎn)')
ax2.set_ylim(Y.min()-10,Y.max()+10)
ax2.set_ylabel('指數(shù)')
ax2.set_title('三次樣條插值')
ax2.legend()
基于《Python 與量化投資 從基礎(chǔ)到實(shí)戰(zhàn)》的內(nèi)容練習(xí)
以上就是本文的全部?jī)?nèi)容，希望對(duì)大家的學(xué)習(xí)有所幫助，也希望大家多多支持腳本之家。
相關(guān)文章
Python使用scrapy采集時(shí)偽裝成HTTP/1.1的方法
這篇文章主要介紹了Python使用scrapy采集時(shí)偽裝成HTTP/1.1的方法,實(shí)例分析了scrapy采集的使用技巧,非常具有實(shí)用價(jià)值,需要的朋友可以參考下2015-04-04
在本篇文章里小編給大家分享了關(guān)于Python中添加自定義模塊的實(shí)例方法，需要的朋友們可以學(xué)習(xí)下。2020-06-06
這篇文章主要介紹了pyinstaller打包找不到文件的問題解決，文中通過示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2020-04-04
這篇文章主要介紹了python多環(huán)境切換及pyenv使用過程詳解,文中通過示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2019-09-09
python 對(duì)給定可迭代集合統(tǒng)計(jì)出現(xiàn)頻率,并排序的方法
今天小編就為大家分享一篇python 對(duì)給定可迭代集合統(tǒng)計(jì)出現(xiàn)頻率,并排序的方法，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過來看看吧2018-10-10
詳解Python如何利用裝飾器實(shí)現(xiàn)重試機(jī)制
重試機(jī)制在編程中是比較常見的場(chǎng)景,主要被用于處理那些可能由于臨時(shí)性故障或網(wǎng)絡(luò)波動(dòng)等原因而失敗的操作,下面我們就來看看如何使用裝飾器來實(shí)現(xiàn)重試機(jī)制吧2024-03-03
基于python代碼實(shí)現(xiàn)簡(jiǎn)易濾除數(shù)字的方法
今天小編就為大家分享一篇基于python代碼實(shí)現(xiàn)簡(jiǎn)易濾除數(shù)字的方法，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過來看看吧2018-07-07
在PDF中添加條碼是一個(gè)常見需求,特別是在需要自動(dòng)化處理、跟蹤或檢索PDF文件時(shí),作為一種機(jī)器可讀的標(biāo)識(shí)符,PDF中的條碼可以包含各種類型的信息,本文介紹了如何使用Python在PDF中添加條形碼或二維碼,需要的朋友可以參考下2024-08-08
Python標(biāo)準(zhǔn)庫之日期、時(shí)間和日歷模塊
這篇文章介紹了Python標(biāo)準(zhǔn)庫之日期、時(shí)間和日歷模塊，文中通過示例代碼介紹的非常詳細(xì)。對(duì)大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值，需要的朋友可以參考下2022-05-05
使用django-guardian實(shí)現(xiàn)django-admin的行級(jí)權(quán)限控制的方法
這篇文章主要介紹了使用django-guardian實(shí)現(xiàn)django-admin的行級(jí)權(quán)限控制的方法，小編覺得挺不錯(cuò)的，現(xiàn)在分享給大家，也給大家做個(gè)參考。一起跟隨小編過來看看吧2018-10-10