离散度计算公式 python_计算离散点的曲率(附Python, MATLAB代码)

在很多学科中的很多计算任务中都需要用到曲线的曲率(或者曲率半径),numpy库里和matlab build-in里都没有现成的能从离散点来算曲率的方法,网上找到的代码又不敢直接用,毕竟是要高频率用到自己科研上的工具,所以决定结合找到的资料自己推一下,并造出python和matlab的轮子,造福后人

公式很简单:

曲率:

equation?tex=%5Ckappa+%3D+%5Cfrac%7B%5Cddot%7B%5Cvec%7Br%7D%7D%5Ctimes%5Cdot%7B%5Cvec%7Br%7D%7D%7D%7B%5Cleft%7C+%5Cdot%7B%5Cvec%7Br%7D%7D%5Cright%7C%5E3%7D

在二维情况下,其标量形式为:

equation?tex=%5Ckappa+%3D+%5Cfrac%7Bx%27%27y%27-x%27y%27%27%7D%7B%5Cleft%28+%28x%27%29%5E2%2B%28y%27%29%5E2%5Cright%29%5E%7B3%2F2%7D%7D

所以对于解析情况非常简单,可以直接对于曲线表达式进行解析求导,但是对于离散的点,情况反倒会比较复杂,因为这里的x和y的一阶和二阶导数如果直接用差分方法来计算的话会造成比较大的误差。

  • 1
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值