如果不用polyshape函数，怎么实现上述算法呢

如果不使用polyshape函数，可以使用另一种方法来计算多边形的交点。具体来说，可以使用边界框法（bounding box method）来计算两个多边形的交点。该方法的基本思路是，首先计算两个多边形的边界框（bounding box），然后将边界框内的点作为采样点，计算采样点与两个多边形的交点。这样就可以得到两个多边形的交点。

下面是一个可能的实现代码，供您参考：

% 计算多边形的边界框
x_min = min(poly(:, 1));
x_max = max(poly(:, 1));
y_min = min(poly(:, 2));
y_max = max(poly(:, 2));

% 计算采样点，并计算采样点与多边形的交点
dx = 0.1;  % 采样点的间距
x_points = x_min:dx:x_max;
y_points = y_min:dx:y_max;
n_points = length(x_points) * length(y_points);
points = zeros(n_points, 2);
k = 1;
for i = 1:length(x_points)
    for j = 1:length(y_points)
        points(k, :) = [x_points(i), y_points(j)];
        k = k + 1;
    end
end
in_poly = inpolygon(points(:, 1), points(:, 2), poly(:, 1), poly(:, 2));
x_intersect = zeros(sum(in_poly), 1);
y_intersect = zeros(sum(in_poly), 1);
k = 1;
for i = 1:length(points)
    if in_poly(i)
        [x_int, y_int] = polyxpoly(poly(:, 1), poly(:, 2), points(i, 1), points(i, 2));
        if ~isempty(x_int)
            x_intersect(k) = x_int(1);
            y_intersect(k) = y_int(1);
            k = k + 1;
        end
    end
end
x_intersect(k:end) = [];
y_intersect(k:end) = [];

注意，这里的poly是多边形的顶点坐标，需要根据实际情况进行修改。另外，这里的计算方法可能会比polyshape函数慢一些，但是可以适用于更一般的多边形。

希望这些信息能对您有所帮助。如果您有更多问题，可以随时向我提问。