使用matlab写一个刀具路径规划算法，实现以下要求：首先，画一个含有多边形岛屿以及一个圆形岛屿的型腔，然后，再通过一系列的平行线与型腔以及岛屿的轮廓边界进行求交，并且，按照单调链法将交点存储起来，（即同一条扫描线上的交点保存在一个数组和链表结构中，同时对链表结构中的各个交点进行编号，编号按照沿水平方向递增的原则增加），然后，用往复的线将这些交点连接起来，并且遍历完除岛屿以外的所有型腔空间，所连接的往复的线一定要避开岛屿，请给我写一个可供参考的代码。

这是一个比较复杂的问题，涉及到几个不同的步骤和算法。以下是一个参考代码，仅供参考，具体实现可能需要根据实际情况进行调整和修改。

首先，我们需要生成一个包含多边形岛屿和圆形岛屿的型腔。这可以通过matlab中的多边形和圆形函数进行实现，例如：

% 多边形岛屿
xpoly = [0.2 0.4 0.6 0.8 0.6 0.4];
ypoly = [0.2 0.4 0.2 0.4 0.6 0.4];
poly = polyshape(xpoly, ypoly);

% 圆形岛屿
[xcirc, ycirc] = circle(0.5, 0.5, 0.1);
circ = polyshape(xcirc, ycirc);

接下来，我们需要通过一系列的平行线与型腔以及岛屿的轮廓边界进行求交，并且按照单调链法将交点存储起来。这可以通过matlab中的intersect函数和monotone_chain函数进行实现，例如：

% 平行线
xline = [0 1];
yline = linspace(0, 1, 11);

% 求交点
[xint, yint] = intersect(poly, xline, yline);
[xint2, yint2] = intersect(circ, xline, yline);
xint = [xint xint2];
yint = [yint yint2];

% 单调链法
[p, q] = monotone_chain(xint, yint);

最后，我们需要用往复的线将这些交点连接起来，并且遍历完除岛屿以外的所有型腔空间。这可以通过matlab中的plot函数进行实现，例如：

% 连接交点
figure;
hold on;
plot(poly);
plot(circ);
for i = 1:length(p)-1
    plot([xint(p(i)) xint(p(i+1))], [yint(p(i)) yint(p(i+1))], 'k');
end

% 遍历型腔空间
[h, v] = size(xline);
for i = 1:h-1
    if mod(i, 2) == 1
        yscan = yline(i+1:-1:i);
    else
        yscan = yline(i:i+1);
    end
    [xint, yint] = intersect(poly, xline(i:i+1), yscan);
    [xint2, yint2] = intersect(circ, xline(i:i+1), yscan);
    xint = [xint xint2];
    yint = [yint yint2];
    [p, q] = monotone_chain(xint, yint);
    for j = 1:length(p)-1
        plot([xint(p(j)) xint(p(j+1))], [yint(p(j)) yint(p(j+1))], 'k');
    end
end

完整代码如下所示：

% 多边形岛屿
xpoly = [0.2 0.4 0.6 0.8 0.6 0.4];
ypoly = [0.2 0.4 0.2 0.4 0.6 0.4];
poly = polyshape(xpoly, ypoly);

% 圆形岛屿
[xcirc, ycirc] = circle(0.5, 0.5, 0.1);
circ = polyshape(xcirc, ycirc);

% 平行线
xline = [0 1];
yline = linspace(0, 1, 11);

% 求交点
[xint, yint] = intersect(poly, xline, yline);
[xint2, yint2] = intersect(circ, xline, yline);
xint = [xint xint2];
yint = [yint yint2];

% 单调链法
[p, q] = monotone_chain(xint, yint);

% 连接交点
figure;
hold on;
plot(poly);
plot(circ);
for i = 1:length(p)-1
    plot([xint(p(i)) xint(p(i+1))], [yint(p(i)) yint(p(i+1))], 'k');
end

% 遍历型腔空间
[h, v] = size(xline);
for i = 1:h-1
    if mod(i, 2) == 1
        yscan = yline(i+1:-1:i);
    else
        yscan = yline(i:i+1);
    end
    [xint, yint] = intersect(poly, xline(i:i+1), yscan);
    [xint2, yint2] = intersect(circ, xline(i:i+1), yscan);
    xint = [xint xint2];
    yint = [yint yint2];
    [p, q] = monotone_chain(xint, yint);
    for j = 1:length(p)-1
        plot([xint(p(j)) xint(p(j+1))], [yint(p(j)) yint(p(j+1))], 'k');
    end
end

需要注意的是，这个代码只是一个简单的示例，实际应用中可能需要进行更多的处理和优化。例如，可能需要考虑交点的顺序问题，以确保连接的线段不会交叉或重叠。同时，可能需要对型腔空间进行更复杂的遍历，以处理更加复杂的岛屿和轮廓边界。