神剑山庄资源网 Design By www.hcban.com
(1)、函数
y = sin(x)
(2)、数据准备
#数据准备 X=np.arange(-np.pi,np.pi,1) #定义样本点X,从-pi到pi每次间隔1 Y= np.sin(X)#定义样本点Y,形成sin函数 new_x=np.arange(-np.pi,np.pi,0.1) #定义差值点
(3)、样条插值
#进行样条差值 import scipy.interpolate as spi #进行一阶样条插值 ipo1=spi.splrep(X,Y,k=1) #样本点导入,生成参数 iy1=spi.splev(new_x,ipo1) #根据观测点和样条参数,生成插值 #进行三次样条拟合 ipo3=spi.splrep(X,Y,k=3) #样本点导入,生成参数 iy3=spi.splev(new_x,ipo3) #根据观测点和样条参数,生成插值
(4)、绘图
##作图 fig,(ax1,ax2)=plt.subplots(2,1,figsize=(10,12)) ax1.plot(X,Y,'o',label='样本点') ax1.plot(new_x,iy1,label='插值点') ax1.set_ylim(Y.min()-1,Y.max()+1) ax1.set_ylabel('指数') ax1.set_title('线性插值') ax1.legend() ax2.plot(X,Y,'o',label='样本点') ax2.plot(new_x,iy3,label='插值点') ax2.set_ylim(Y.min()-1,Y.max()+1) ax2.set_ylabel('指数') ax2.set_title('三次样条插值') ax2.legend()
(5)、综合代码
import numpy as np import matplotlib.pyplot as plt plt.rcParams['font.sans-serif']=['SimHei'] #用来正常显示中文标签 plt.rcParams['axes.unicode_minus']=False #用来正常显示负号 #数据准备 X=np.arange(-np.pi,np.pi,1) #定义样本点X,从-pi到pi每次间隔1 Y= np.sin(X)#定义样本点Y,形成sin函数 new_x=np.arange(-np.pi,np.pi,0.1) #定义差值点 #进行样条差值 import scipy.interpolate as spi #进行一阶样条插值 ipo1=spi.splrep(X,Y,k=1) #样本点导入,生成参数 iy1=spi.splev(new_x,ipo1) #根据观测点和样条参数,生成插值 #进行三次样条拟合 ipo3=spi.splrep(X,Y,k=3) #样本点导入,生成参数 iy3=spi.splev(new_x,ipo3) #根据观测点和样条参数,生成插值 ##作图 fig,(ax1,ax2)=plt.subplots(2,1,figsize=(10,12)) ax1.plot(X,Y,'o',label='样本点') ax1.plot(new_x,iy1,label='插值点') ax1.set_ylim(Y.min()-1,Y.max()+1) ax1.set_ylabel('指数') ax1.set_title('线性插值') ax1.legend() ax2.plot(X,Y,'o',label='样本点') ax2.plot(new_x,iy3,label='插值点') ax2.set_ylim(Y.min()-1,Y.max()+1) ax2.set_ylabel('指数') ax2.set_title('三次样条插值') ax2.legend()
(6)、结果展示
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持。
神剑山庄资源网 Design By www.hcban.com
神剑山庄资源网
免责声明:本站文章均来自网站采集或用户投稿,网站不提供任何软件下载或自行开发的软件!
如有用户或公司发现本站内容信息存在侵权行为,请邮件告知! 858582#qq.com
神剑山庄资源网 Design By www.hcban.com
暂无Python实现线性插值和三次样条插值的示例代码的评论...
更新日志
2024年09月30日
2024年09月30日
- 孙悦.1996-伙伴【正大国际】【WAV+CUE】
- 纪钧瀚《钢琴阅读时光 雨中书店聆听轻音乐》[FLAC/分轨][399.62MB]
- 证声音乐图书馆《走向自然 疗心爵士乐》[320K/MP3][87.4MB]
- 证声音乐图书馆《走向自然 疗心爵士乐》[FLAC/分轨][184.94MB]
- 陈慧娴.2018-Priscilla-Ism演唱会3CD(2024环球红馆40复刻系列)【环球】【WAV+CUE】
- 郑秀文.1999-我应该得到(国)【华纳】【WAV+CUE】
- 陈家慧.2011-钢琴酒吧2CD【龙吟唱片】【WAV+CUE】
- 证声音乐图书馆《雨季 蓝调吉他 Rainy Blues》[320K/MP3][45.01MB]
- 证声音乐图书馆《雨季 蓝调吉他 Rainy Blues》[FLAC/分轨][109.13MB]
- 赞多《序章》[320K/MP3][45.54MB]
- 许巍.2004-每一刻都是崭新的【步升大风】【WAV+CUE】
- 群星.2024-四方馆影视原声带【韶愔音乐】【FLAC分轨】
- 陈雷.1997-安锁咧【金圆唱片】【WAV+CUE】
- 关淑怡.2013-MY.FAVORITE.SK.3CD【环球】【WAV+CUE】
- Sweety.2006-花言乔语【丰华】【WAV+CUE】