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

matlab interp1 c,SciPy interp1d结果与MatLab interp1不同[通俗易懂]我正在将一个MatLab程序转换为Python,我很难理解为什么scipy.interpolate.interp1d给出的结果与MatlabInterp1不同。在MatLab中,用法略有不同:yi=interp1(x,Y,xi,’cubic’)剪刀:f=interp1d(x,Y,kind=’cubic’)yi=f(xi)对于一个简单的例子,结果是相同的:MatLab软件:interp1…

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

我正在将一个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.net/131343.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • Windows 定时执行脚本[通俗易懂]

    Windows 定时执行脚本[通俗易懂]Windows定时执行脚本1、参考资料windows设置定时执行脚本2、解决办法由于Windows系统,无法使用Linux下强大的crontab命令,所以该定时任务也是针对Windo系统的具体操作步骤右击【我的电脑】,选择【管理】在【任务计划程序】面板中点击【创建基本任务】输入任务的【名称】和【描述】当然是希望每天自动执行防疫打卡啦~设置每天执行任务的时间,以及每隔几天执行一次该任务选择【启动程序】选择启动程序的startup.bat

    2022年7月17日
    57
  • 图片url地址的生成获取方法

    图片url地址的生成获取方法图片url地址的生成获取方法  在写博客插入图片时,许多时候需要提供图片的url地址。作为菜鸡的我,自然是一脸懵逼。那么什么是所谓的url地址呢?又该如何获取图片的url地址呢?  首先来看一下度娘对url地址的解释:url是统一资源定位符,对可以从互联网上得到的资源的位置和访问方法的一种简洁的表示,是互联网上标准资源的地址。互联网上的每个文件都有一个唯一的url,它包含的信息指出文件的位置以及浏览器应该怎么处理它。  简单来说,url地址是是用来定位、访问网上资源用的。常见的网址也属于url地

    2022年9月22日
    0
  • 遍历Arraylist的三种方法及优缺点简单介绍

    遍历Arraylist的三种方法及优缺点简单介绍集合ArrayList是接口List的一种子类,它的特点是:存储的元素是有序的.底层的数据结构是数组.查询快,增删慢.在众多集合中ArrayList的遍历又是比较特殊的,下面就写一下它的三种遍历方式,代码如下:第一种遍历方式:普通for循环第二种遍历方式:增强for循环第三种遍历方式:迭代器importjava.util.ArrayList;importjava.util.Iterator;/***PACKAGE_NAME*/publicclassDemo.

    2022年7月22日
    6
  • apache 配置跨域

    apache 配置跨域<VirtualHost*:80>DocumentRoot“D:/project/xuanhua_shop/public”ServerName192.168.18.182#servername后的ip为接口所在服务器IP<Directory“D:/project/xuanhua_shop/public”>OptionsIndexesFollowSymLinksMultiViews Options+Indexes+Includes+FollowSymLi

    2022年7月14日
    11
  • 2021最新Python全栈工程师学习路线(初级+高级+大神)「建议收藏」

    2021最新Python全栈工程师学习路线(初级+高级+大神)「建议收藏」2019最新Python全栈工程师学习路线(初级+高级+大神)IT行业,技术要比学历、年龄、从业经验更为重要,技术水平直接决定就业薪资,想要学好python,首先要先了解精通Python语言基础、Pythonweb开发、Python爬虫、Python数据分析这四大方面。零基础学习需要要从如下几个阶段入手:阶段一:熟练掌握Python多线程并发编程技术,可以编写爬虫程序和语音识别软件…

    2022年10月26日
    0
  • NL-Mean和BM3D去噪原理「建议收藏」

    NL-Mean和BM3D去噪原理「建议收藏」本文为转载,原博客地址:http://blog.csdn.net/jacke121/article/details/55215416图像去噪是非常基础也是非常必要的研究,去噪常常在更高级的图像处理之前进行,是图像处理的基础。可惜的是,目前去噪算法并没有很好的解决方案,实际应用中,更多的是在效果和运算复杂度之间求得一个平衡,再一次验证了我老师的一句话:所有的工程问题最后都是最优化问题

    2022年5月16日
    46

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号