![](https://img-blog.csdnimg.cn/20201014180756754.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
python实现插值法
文章平均质量分 76
小陈皓
自动化小白一枚
展开
-
python样条插值的实现(四)
上一次我们完成了对第二类边界条件的样条插值的函数编写,现在的问题是:如何求函数上的最大曲率以我非常有限的python知识,想到了三种办法:1.先将曲率的表达式给出,再给足x的采样点,得到x对应的曲率的值列表,利用列表的寻找最大值函数来找曲率最大值当然这个方法弊端比较明显,首先是精度问题,精度肯定不是很高,并且计算量也较大,这里就不给出代码了(因为被我删了)当然也可以利用取样比较再迭代的方法求得高精度,不过这是后话。不过max函数有个小技巧,例如一个嵌套列表:[[1,50],[2,40],原创 2021-05-14 10:57:09 · 1321 阅读 · 0 评论 -
2021-05-13python样条插值实现(三)
要实现自己写样条插值的python代码,首先得了解样条插值的计算公式。这里采取的是第二类边界条件的样条函数:第二类边界条件的计算公式如下:上图为《数值分析(第五版)》清华大学出版社的关于第二类边界情况下样条插值的计算公式,要将其转化为python代码,需解决两个问题:1.利用python解线性方程组在这里,我用的是numpy.linalg.solve()函数,依旧先给出官方定义:官方文档:numpy.linalg.solve(a, b)举个栗子:imp..原创 2021-05-14 00:20:47 · 961 阅读 · 3 评论 -
python样条插值(二)
python实现样条插值在MATLB中有自带的库,在python中也有,位于scipy库中,具体定义如下:scipy官方文档:class scipy.interpolate.interp1d(x, y, kind='linear', axis=-1, copy=True, bounds_error=None, fill_value=nan, assume_sorted=False)上为定义,下为参数解释:interp1d ,scipy.interpolate包里有很多的模块可以.原创 2021-05-13 23:39:24 · 6534 阅读 · 1 评论 -
2021-05-13 python样条插值(一)
python在数值分析中的运用最近在做数分大作业,因为对matlab不太熟悉,尝试用python来解大作业。首先,上题:捕捞机器人轨迹规划如下图3所示,水面捕捞机器人以船体的形式在水面上自动行驶,通过摄像机拍摄照片获得水面上漂浮物的坐标,驱动船体前进至所需位置捕捞。摄像机获取的信息为水面上每个区域的中心点(图(b)),由于摄像处理信息有延迟,且机器人船体行进过某区域时不便再次转向获取某水面信息,需要对行驶轨迹提前进行规划,以成功捕捞漂浮物。图(c)表示采用单目标识别的方法获得图片上物体的位坐原创 2021-05-13 20:24:46 · 1034 阅读 · 1 评论