标签：dl end 练习 曲率 matlab dx dy xp

曲线参数方程：

的曲率公式如下：

matlab代码如下：

clear all;close all;clc;

x = 0:0.01:2*pi;
y = sin(x);

dx  = 0.5*(x(3:end)-x(1:end-2));
dy  = 0.5*(y(3:end)-y(1:end-2));
dl  = sqrt(dx.^2 + dy.^2);
xp  = dx./dl;
yp  = dy./dl;
% approximate 2nd derivatives of x & y with discrete differences
xpp = (x(3:end)-2*x(2:end-1)+x(1:end-2))./(dl.^2);
ypp = (y(3:end)-2*y(2:end-1)+y(1:end-2))./(dl.^2);

% Compute the curvature
K = (xp.*ypp - yp.*xpp) ./ ((xp.^2 + yp.^2).^(1.5));

subplot(2,1,1);
plot(x,y);
title('曲线')
subplot(2,1,2);
plot(K);
title('曲率')

结果如下：

之前还有写过Mathematica版的三维曲线曲率挠率计算，可以参考这里。

本文代码参考了stackoverflow上的一篇回答。

标签：dl,end,练习,曲率,matlab,dx,dy,xp
来源： https://www.cnblogs.com/tiandsp/p/16536759.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。