Python实现线性插值和三次样条插值的示例代码_python三次样条插值 - CSDN文库

三次样条插值

5星 · 超过95%的资源 37 浏览量更新于2023-03-03 评论 4 收藏 59KB PDF 举报

身份认证购VIP最低享 7 折!

领优惠券(最高得80元）

资源详情

资源评论

资源推荐

Python实现线性插值和三次样条插值的示例代码实现线性插值和三次样条插值的示例代码

主要介绍了Python实现线性插值和三次样条插值的示例代码，文中通过示例代码介绍的非常详细，对大家的学

习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧

(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='样本点')

本内容试读结束，登录后可阅读更多

下载后可阅读完整内容，剩余1页未读，立即下载

评论5

陌陌的日记

2023-07-24

阅读这个文件后，我对线性插值和三次样条插值有了更深入的了解，可以应用到我的实际项目中。

weixin_38723105

粉丝: 4
资源: 969

会员权益专享

图片转文字

全年可省5，000元立即开通

最新资源

资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！点击此处反馈