离散点的曲率的计算
来源:互联网 发布:南京医科大学数据库 编辑:程序博客网 时间:2024/05/11 02:12
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('模拟曲率');
0 0
- 离散点的曲率计算
- 离散点的曲率的计算
- 求取离散点的曲率
- 图像曲率的计算
- 054 曲率的计算
- 离散曲率计算
- 如何计算图像的曲率
- 曲率的计算和曲率图的采样
- vtk计算Polydata的表面曲率
- 曲面的主曲率与高斯曲率计算
- 求曲线某一点处的曲率圆
- 网格离散曲率算法(利用Normal cycle 理论计算)
- 三维图像的曲率
- 曲面的各种“曲率”
- 摇杆的死区和曲率
- 离散型最值的期望计算
- 使用自适应阈值和支持动态区域的曲率尺度空间角点检测(分析)
- 使用自适应阈值和支持动态区域的曲率尺度空间角点检测(理论)
- Qt Stylesheet
- ROS action
- iOS开发 一句话隐藏导航栏下的黑线
- java集合迭代删除,一个很诡异的现象
- LeetCode 557. Reverse Words in a String III
- 离散点的曲率的计算
- 自己做个MP3_软件
- springmvc 参数绑定
- SSD: Single Shot MultiBox Detector 训练KITTI数据集(2)
- windows服务器如何开启防火墙设置
- Elastic Kibana 5.2 中文文档 | 那伊抹微笑
- C++-类与结构体的区别
- 文章标题
- JDK环境变量配置