matlab中求二维离散曲线的曲率，写出matlab代码

以下是求解二维离散曲线曲率的Matlab代码：

% 假设二维离散曲线为(x, y)
dx = diff(x);
dy = diff(y);
d2x = diff(x, 2);
d2y = diff(y, 2);

% 计算曲率
k = (dx(1:end-1).*d2y - dy(1:end-1).*d2x) ./ (dx(1:end-1).^2 + dy(1:end-1).^2).^1.5;

% 将首尾点的曲率设为0
k = [0; k; 0];

其中，`diff()` 函数用于计算差分，`.^` 表示对每个元素进行幂运算。