前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
社区首页 >专栏 >matlab interp1 c,SciPy interp1d结果与MatLab interp1不同[通俗易懂]

matlab interp1 c,SciPy interp1d结果与MatLab interp1不同[通俗易懂]

作者头像
全栈程序员站长
发布2022-07-01 13:04:33
发布2022-07-01 13:04:33
2760
举报

大家好,又见面了,我是你们的朋友全栈君。

我正在将一个MatLab程序转换为Python,我很难理解为什么scipy.interpolate.interp1d给出的结果与MatlabInterp1不同。

在MatLab中,用法略有不同:yi = interp1(x,Y,xi,’cubic’)

剪刀:f = interp1d(x,Y,kind=’cubic’)

yi = f(xi)

对于一个简单的例子,结果是相同的:

MatLab软件:interp1([0 1 2 3 4], [0 1 2 3 4],[1.5 2.5 3.5],’cubic’)

1.5000 2.5000 3.5000

Python:interp1d([1,2,3,4],[1,2,3,4],kind=’cubic’)([1.5,2.5,3.5])

array([ 1.5, 2.5, 3.5])

但对于现实世界的例子来说,它们并不相同:x = 0.0000e+000 2.1333e+001 3.2000e+001 1.6000e+004 2.1333e+004 2.3994e+004

Y = -6 -6 20 20 -6 -6

xi = 0.00000 11.72161 23.44322 35.16484… (2048 data points)

Matlab软件:-6.0000e+000

-1.2330e+001

-3.7384e+000

7.0235e+000

7.0028e+000

6.9821e+000

剪刀:array([[ -6.00000000e+00],

[ -1.56304101e+01],

[ -2.04908267e+00],

…,

[ 1.64475576e+05],

[ 8.28360759e+04],

[ -5.99999999e+00]])

关于如何得到与MatLab一致的结果有什么想法吗?

编辑:我知道三次插值算法的实现有一定的自由度,这可能解释了我看到的差异。似乎我正在转换的原始MatLab程序应该使用线性插值,所以这个问题可能是没有意义的。

发布者:全栈程序员栈长,转载请注明出处:https://javaforall.cn/131343.html原文链接:https://javaforall.cn

本文参与 腾讯云自媒体同步曝光计划,分享自作者个人站点/博客。
如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档