使用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);

% 型腔
xbox = [0 1 1 0];
ybox = [0 0 1 1];
box = polyshape(xbox, ybox);

% 平行线
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];
[xint3, yint3] = intersect(box, xline, yline);
xint = [xint xint3];
yint = [yint yint3];

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

% 连接交点
figure;
hold on;
plot(poly);
plot(circ);
plot(box);
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];
    [xint3, yint3] = intersect(box, xline(i:i+1), yscan);
    xint = [xint xint3];
    yint = [yint yint3];
    [p, q] = monotone_chain(xint, yint);
    for j = 1:length(p)-1
        if ~inShape(poly, xint(p(j)), yint(p(j))) && ~inShape(circ, xint(p(j)), yint(p(j))) && ~inShape(box, xint(p(j)), yint(p(j)))
            plot([xint(p(j)) xint(p(j+1))], [yint(p(j)) yint(p(j+1))], 'k');
        end
    end
end

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