请用MATLAB写一段代码，要求根据这些点（45，20）（75,60），（105,60）（135,20）（165，-60）（225，-100）（255,20）的坐标，希望通过插值获得近处的轨道曲线规划，要求轨迹拐点曲率不可过小，同时近处轨迹规划的误差精度要求不高

以下是一种可能的解法：

% 给出已知点的坐标
x = [45 75 105 135 165 225 255];
y = [20 60 60 20 -60 -100 20];

% 定义插值点密度，这里取每个像素点一个插值点
density = 1;

% 计算插值点的坐标
xi = x(1):density:x(end);
yi = interp1(x, y, xi);

% 计算插值点的一阶导数，用于求曲率
dyi = gradient(yi, xi);

% 计算曲率，注意要避免除以零
curvature = abs(dyi(2:end) .* diff(dyi)) ./ ((1 + dyi(2:end).^2).^1.5);
curvature(isnan(curvature)) = 0;

% 按照曲率将轨迹分段，从而确保拐点曲率不过小
threshold = 0.1; % 这里取 0.1 作为阈值，可以根据实际情况调整
segment_idx = [1 find(curvature > threshold) length(xi)];

% 绘制轨迹图像
figure;
plot(xi, yi, '-');
hold on;
for i = 1:length(segment_idx)-1
    plot(xi(segment_idx(i):segment_idx(i+1)), yi(segment_idx(i):segment_idx(i+1)), 'r-', 'LineWidth', 2);
end
axis equal;
title('轨迹规划结果');

上述代码中，我们首先利用 `interp1` 函数对已知点进行线性插值，得到一系列密集的插值点。然后计算了每个插值点的一阶导数，用于求曲率。为了确保拐点曲率不过小，我们将曲率大于某个阈值的点作为轨迹的分段点，从而将轨迹划分成多个小段。最后，我们绘制了轨迹图像，其中红色线段表示分段后的小段。