概述
from：http://www.ilovematlab.cn/thread-67207-1-1.html
曲率公式：
注意：
①曲率半径为曲率的倒数
②如果是离散点，先用polyfit和polyval拟合出曲线
程序：
clc; clear all; close all;
x0 = linspace(0, 1);
y0 = sin(x0).*cos(x0);
h = abs(diff([x0(2), x0(1)]));
% 模拟一阶导
figure; box on; hold on;
ythe1 = cos(x0).^2 - sin(x0).^2; %理论一阶导
yapp1 = gradient(y0, h); %matlab数值近似
plot(x0, ythe1, '.');
plot(x0, yapp1, 'r');
legend('理论值', '模拟值');
title('模拟一阶导');
% 模拟二阶导
figure; box on; hold on;
ythe2 = (-4)*cos(x0).*sin(x0); %理论二阶导
yapp2 = 2*2*del2(y0, h); %matlab数值近似
plot(x0, ythe2,'.');
plot(x0, yapp2,'r');
legend('理论值', '模拟值');
title('模拟二阶导');
% 模拟曲率
syms x y
y = sin(x)*cos(x);
yd2 = diff(y, 2);
yd1 = diff(y, 1);
k = abs(yd2)/(1+yd1^2)^(3/2);
k1 = subs(k, x, x0);
k2 = abs(yapp2)./(1+yapp1.^2).^(3/2);
figure; box on; hold on;
plot(x0, k1, '.');
plot(x0, k2, 'r');
legend('理论值', '模拟值', 'Location', 'NorthWest');
title('模拟曲率');
最后
以上就是震动店员为你收集整理的离散点的曲率的计算的全部内容，希望文章能够帮你解决离散点的曲率的计算所遇到的程序开发问题。
如果觉得靠谱客网站的内容还不错，欢迎将靠谱客网站推荐给程序员好友。
本图文内容来源于网友提供，作为学习参考使用，或来自网络收集整理，版权属于原作者所有。
发表评论 取消回复