MATLAB 如何求取离散点的曲率最大值

我就废话不多说了，大家还是直接看代码吧~

x0 = linspace(0.1,2,100);%x0,y0验证函数离散点，可以非等间隔
y0 = 1./x0;
h1 = abs(diff(x0)) ;
h = [h1 h1(end)];
ht = h;
yapp1 = gradient(y0)./ht; %matlab数值近似
yapp2 = del2(y0)./ht; %matlab数值近似
k2 = abs(yapp2)./(1+yapp1.^2).^(3/2);
figure
plot(k2)
title('曲率曲线')
[~,maxFlag] = max(k2);%曲率最大位置
x_max = x0(maxFlag);
y_max = y0(maxFlag);
%画出图像 标注曲率最大点
figure
plot(x0,y0,'.-');
hold on;
plot(x_max,y_max,'rp')
title('标注最大曲率点')
xlabel('log10((norm(B*Xk-L)))')
ylabel('log10((norm(Xk)))')

补充：MATLAB 插值+计算离散点曲率

思路：点足够密的话直接用 diff、gradient 求曲率，稀疏的话先插值再算曲率。

公式：

点密的情况 输入曲线坐标（1-2）求一、二阶导数（4-9）通过公式求得曲率（10）

x = 0:0.01:7;
y = cos(x*0.5*pi);
h1 = abs(diff(x));
h = [h1 h1(end)];
ht = h;
y1 = gradient(y)./ht;
y2 = gradient(y1)./ht;
curv = abs(y2)./sqrt((1+y1.^2).^3);
plot(x,y,'-',x,curv,'--r);
legend('Raw Data, 'Curvature','Location',"best");
grid on

图像与下文理论值图像相同

点稀疏的情况

1、输入散点坐标（1-2）

2、用样条曲线（B-Spline）等方法插值得到拟合曲线（3-4）

3、diff、gradient 函数求拟合曲线的一、二阶导数（6-11）

4、通过公式求得曲率（12）

例：余弦函数取 8 个点，用 B-Spline 插值

x = 0:1:7;
y = cos(x*0.5*pi);
xx = 0:0.01:7;
yy = spline(x,y,xx);
h1 = abs(diff(xx));
h = [h1 h1(end)];
ht = h;
yy1 = gradient(yy)./ht;
yy2 = gradient(yy1)./ht;
curv = abs(yy2)./sqrt((1+yy1.^2).^3);
plot(xx,yy,'-',xx,curv,'--r',x,y,'o-');
legend('B-Spline', 'Curvature','Raw Data','Location',"best");
grid on

补充用法

求最大曲率并在图中标出

[max_val,max_ind]=max(curv);
hold on
plot(xx(max_ind),yy(max_ind),'*r');

与理论值（余弦函数曲线）对比

曲线对比

曲率对比

几种插值方法对比

列举四种方法，分别为：分段线性插值、三次样条曲线（B-Spline）插值、三次 Hermite 插值（PCHIP）、修正 Akima 分段三次 Hermite 插值（Akima）

Case 1： 三维螺线

三维螺线散点

插值

俯视

侧视

Case 2：二维梯形波

二维梯形波

Case 3：三维不规则折线

三维不规则折线（不等间距）

对比可得：

Case 1：B-Spline>Akima>PCHIP>Linear

Case 2：Linear>PCHIP>Akima>B-Spline

Case 3：Linear≈PCHIP≈Akima>B-Spline

故在插值的时候需要选择适合的计算方法

以上为个人经验，希望能给大家一个参考，也希望大家多多支持我们。如有错误或未考虑完全的地方，望不吝赐教。

赞 (0)

相关推荐

  • python matlab库简单用法讲解

    大部分的库都是对应不同的接口使用,然而本期给大家带来的库有点与众不同,涉及的方向是图形化的界面,使用上并不困难,但是也存在着需要注意的使用的,比如安装方面,不同于其他的第三方库,安装的方式不一样,且还需要测试安装后的效果,不然很难去掌控它,好了,下面就来详细的介绍. matlab库安装: python setup.py install 验证安装成功: matlabengineforpython R2017a 语法: matlab.plot([],'') 使用代码: import matlab.e

  • 解决python调用matlab时的一些常见问题

    为什么要用python调用matlab? 我自己的有些数据结构涉及到hash查找,在python中key是tuple形式,在matlab中支持hash查找的数据结构只有containers.Map(),并不能支持cell作为key. 尝试过把向量转为string,但是num2str和str2mat的效率不高,containers.Map()的查找耗时也非常的长,所以只好作罢. 后来发现可以用python通过matlab的API直接调用matlab的函数参与运算,朋友的经验说矩阵运算都交给matl

  • 如何基于matlab相机标定导出xml文件

    1 参数选择 径向畸变3个参数还是两个参数 默认两个参数 如果是三个参数 2准备转化生成结果 二参数的转化代码 writeExternalandIntrinsicMatrix(cameraParams62,'cameraParams622.xml'); function writeExternalandIntrinsicMatrix(cameraParams,file) %writeXML(cameraParams,file) docNode = com.mathworks.xml.XMLUti

  • windows下python 3.9 Numpy scipy和matlabplot的安装教程详解

    学习python过程中想使用python的matlabplot绘图功能,遇到了一大批问题,然后一路过关斩将,最终安装成功,实为不易,发帖留念. 1 首先打开cmd win+r 2 pip安装 pip3 install --user numpy scipy matplotlib –user 选项可以设置只安装在当前的用户下,而不是写入到系统目录.默认情况使用国外线路,国外太慢,我们使用清华的镜像就可以: pip3 install numpy scipy matplotlib -i https://

  • Matlab求解数组中的最大值及它所在的具体位置

    Matlab求解数组中的最大值及它所在的具体位置

    在使用Matlab肯定会碰到Matlab求解数组中的最大值以及它所在的位置的问题.博主开始用循环的方法找,既浪费时间又消耗资源,后面查找后才发现有简单快速的方法.下面就简单介绍一下这种方法. 1.电脑环境 电脑环境:Windows 10 教育版 MATLAB:MATLAB R2014a 2.方法 1.1.一维数组 在Matlab随机生成一维数组或者手动输入 a = [1,9,24,8,7,16] [m,p]=max(a) 利用max函数,将一维数据放入max(一维数组)中 最终利用Max函数输出

  • 如何用Matlab和Python读取Netcdf文件

    如何用Matlab和Python读取Netcdf文件

    NetCDF(网络通用数据格式)文件后缀为.nc,文件中的数据结构包含维(dimensions).变量(variables)和属性(attributes)三种描述类型,每种类型都会被分配一个名字和一个ID. 从数学上来说,netcdf存储的数据就是一个多自变量的单值函数.用公式来说就是f(x,y,z,...)=value,函数的自变量x,y,z等在netcdf中叫做维(dimension)或坐标轴(axis),函数值value在netcdf中叫做变量(Variables).而自变量和函数值在物理

  • MATLAB 如何求取离散点的曲率最大值

    MATLAB 如何求取离散点的曲率最大值

    我就废话不多说了,大家还是直接看代码吧~ x0 = linspace(0.1,2,100);%x0,y0验证函数离散点,可以非等间隔 y0 = 1./x0; h1 = abs(diff(x0)) ; h = [h1 h1(end)]; ht = h; yapp1 = gradient(y0)./ht; %matlab数值近似 yapp2 = del2(y0)./ht; %matlab数值近似 k2 = abs(yapp2)./(1+yapp1.^2).^(3/2); figure plot(k2

  • Java实现求子数组和的最大值算法示例

    Java实现求子数组和的最大值算法示例

    本文实例讲述了Java实现求子数组和的最大值算法.分享给大家供大家参考,具体如下: 一般C和C++在算法实现中使用较多,下面我们通过java语言实现算法,更有亲切感. 题目: 输入一个整形数组,数组里有正数也有负数. 数组中连续的一个或多个整数组成一个子数组,每个子数组都有一个和. 求所有子数组的和的最大值. 例如输入的数组为1, -2, 3, 10, -4, 7, 2, -5,和最大的子数组为3, 10, -4, 7, 2, 因此输出为该子数组的和18. 实现代码: package arrDe

  • 利用C#版OpenCV实现圆心求取实例代码

    利用C#版OpenCV实现圆心求取实例代码

    前言 OpenCVSharp是OpenCV的.NET wrapper,是一名日本工程师开发的,项目地址为:https://github.com/shimat/opencvsharp. 该源码是 BSD开放协议,BSD开源协议是一个给于使用者很大自由的协议.基本上使用者可以"为所欲为",可以自由的使用,修改源代码,也可以将修改后的代码作为开源或者专有软件再发布或商业化销售. 1.OpenCVSharp的下载 可以直接从上面的github上下载源码,自行编译引用: 也可用vs中的nuget

  • python多维数组分位数的求取方式

    在python中计算一个多维数组的任意百分比分位数,只需用np.percentile即可,十分方便 import numpy as np a = [154, 400, 1124, 82, 94, 108] print np.percentile(a,95) # gives the 95th percentile 补充拓展:如何解决hive同时计算多个分位数的问题 众所周知,原生hive没有计算中位数的函数(有的平台会有),只有计算分位数的函数percentile 在数据量不大的时候,速度尚可.但

  • Opencv求取连通区域重心实例

    我们有时候需要求取某一个物体重心,这里一般将图像二值化,得出该物体的轮廓,然后根据灰度重心法,计算出每一个物体的中心. 步骤如下: 1)合适的阈值二值化 2)求取轮廓 3)计算重心 otsu算法求取最佳阈值 otsu法(最大类间方差法,有时也称之为大津算法)使用的是聚类的思想,把图像的灰度数按灰度级分成2个部分,使得两个部分之间的灰度值差异最大,每个部分之间的灰度差异最小,通过方差的计算来寻找一个合适的灰度级别来划分,otsu算法被认为是图像分割中阈值选取的最佳算法,计算简单,不受图像亮度和对比

  • Python 马氏距离求取函数详解

    Python 马氏距离求取函数详解

    马氏距离区别于欧式距离,如百度知道中所言: 马氏距离(Mahalanobis distance)是由印度统计学家马哈拉诺比斯(P. C. Mahalanobis)提出的,表示点与一个分布之间的距离.它是一种有效的计算两个未知样本集的相似度的方法.与 欧氏距离不同的是,它考虑到各种特性之间的联系(例如:一条关于身高的信息会带来一条关于体重的信息,因为两者是有关联的),并且是尺度无关的(scale-invariant),即独立于测量尺度.对于一个均值为μ, 协方差矩阵为Σ的多变量向量,其马氏距离为s

  • Python Requests爬虫之求取关键词页面详解

    Python Requests爬虫之求取关键词页面详解

    目录 需求:爬取搜狗首页的页面数据 使用UA伪装 求取关键词页面 总结 需求:爬取搜狗首页的页面数据 import requestsif __name__=='__main__': #step 1:搜索Url url='https://123.sogou.com/' #step 2:发起请求 #get方法会返回一个响应对象 response=requests.get(url=url) #step 3:获取响应数据,text返回的是字符串形式的响应数据 page_text=response.tex

  • Java利用Dijkstra和Floyd分别求取图的最短路径

    Java利用Dijkstra和Floyd分别求取图的最短路径

    目录 1 最短路径的概述 2 杰斯特拉(Dijkstra)算法 2.1 原理 2.2 案例分析 3 弗洛伊德(Floyd)算法 3.1 原理 3.2 案例分析 4 邻接矩阵加权图实现 5 邻接表加权图实现 本文详细介绍了图的最短路径的概念,然后介绍了求最短路径的两种算法:Dijkstra算法和Floyd算法的原理,最后提供了基于邻接矩阵和邻接表的图对两种算法的Java实现. 阅读本文需要一定的图的基础,如果对于图不是太明白的可以看看这篇文章:Java数据结构之图的原理与实现. 1 最短路径的概述

  • C++调用Matlab函数求特征值

    最近需要用到C++和Matlab的混编,记录一下学习过程. 要实现的是调用Matlab函数,求矩阵前k个最小的特征值及其特征向量. //C++ #include "engine.h" //使用Matlab引擎需要包含的头文件 #include <iostream> using namespace std; int main() { Engine *m_engine; //创建Matlab引擎 m_engine = NULL; //初始化引擎 if((!m_engine &a

  • MATLAB Delaunay算法提取离散点边界的方法

    最近在项目进行中遇到要提取离散点边界的问题,像我这样的对于matlab不是特别熟练的朋友一开始肯定摸不着头脑,到底选用哪种算法可以有效地提取到所有已知点的轮廓线呢.本人经过大量的文献搜索及代码实验找到了几个效果比较好的轮廓提取代码,在这里做个总结,并且希望能够对遇到同样问题的朋友有所启发. 关于离散点边界提取的三种方法: 1.Convhull 离散点集获得边界 2.Alpha Shape算法检测边缘点 3.Delaunay 三角剖分算法 前两种方法在之前的博客中已经做了总结这里就不展开了,现在主

随机推荐