interpolate.splev 函数 利用B样条和它的导数进行插值,interpolate.splev课展现二维曲线插值

#p100下半部分
#4.1.2scipy
#例子2 interpolate.splev 函数  利用B样条和它的导数进行插值,interpolate.splev课展现二维曲线插值
import numpy as np
import matplotlib.pyplot as plt
import scipy as sp
from scipy.interpolate import splprep #增加该行


t=np.arange(0,2.5,.10)#0到2.5之间生成一个一维数组,差值0.10
x=np.sin(2*np.pi*t)
y=np.cos(2*np.pi*t)
tcktuples,uarray=sp.interpolate.splprep([x,y],s=0)
unew=np.arange(0,1.01,0.01)
splinevalues=sp.interpolate.splev(unew,tcktuples)
          
plt.figure(figsize=(10,10))
#plt.plot(x,y,'x',splinevalues[0],splinevalues[1],np.sin(2*np.pi*unew),np.cos(2*np.pi*unew),x,y,'b')
plt.plot(x,y,'x')#linear#圆--离散点
#plt.plot(splinevalues[0],splinevalues[1])#圆
#plt.plot(np.sin(2*np.pi*unew),np.cos(2*np.pi*unew))#圆
#plt.plot(np.sin(2*np.pi*unew))#线
#plt.plot(np.cos(2*np.pi*unew))#线
#plt.plot(x,y,'b')#多边形


plt.legend(['Linear','Cubic Spline','True'])
plt.axis([-1.25,1.25,-1.25,1.25])
plt.title('Parametric Spline Interpolation Curve')
plt.show()
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值