weixin_44028636的博客

曲率公式： 注意： ①曲率半径为曲率的倒数 ②如果是离散点，先用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)