可以给一个matlab的代码，一个矩形区域内部有障碍物，采用梯形分解法分割成若干个小梯形单元，并展示分割后的图像，输出坐标点。

以下是一个 MATLAB 代码示例，用于将包含障碍物的矩形区域分割成若干个小梯形单元，并展示分割后的图像，输出坐标点。

% 定义矩形区域的四个顶点坐标
x1 = 0; y1 = 0;
x2 = 10; y2 = 0;
x3 = 10; y3 = 10;
x4 = 0; y4 = 10;

% 定义矩形区域中的障碍物坐标
obstacle = [3, 3; 4, 3; 4, 4; 3, 4];

% 将矩形区域分割成若干个小梯形单元
[x, y] = polybool('subtraction', [x1, x2, x3, x4], [y1, y2, y3, y4], obstacle(:,1), obstacle(:,2));
n_rect = length(x) - 1;

% 绘制分割后的图像
figure;
hold on;
for i = 1:n_rect
    patch([x(i), x(i+1), x(i+1), x(i)], [y(i), y(i), y(i+1), y(i+1)], 'r');
end
axis equal;
axis([x1, x3, y1, y3]);
xlabel('x');
ylabel('y');
title('分割后的矩形区域');

% 输出坐标点
disp(['矩形区域中共有 ', num2str(n_rect), ' 个小梯形单元。']);
for i = 1:n_rect
    disp(['梯形单元 ', num2str(i), ' 的坐标点为：']);
    disp(['(', num2str(x(i)), ', ', num2str(y(i)), '),']);
    disp(['(', num2str(x(i+1)), ', ', num2str(y(i+1)), '),']);
    disp(['(', num2str(x(n_rect+i+1)), ', ', num2str(y(n_rect+i+1)), '),']);
    disp(['(', num2str(x(n_rect+i)), ', ', num2str(y(n_rect+i)), ').']);
end

请注意，这只是一个简单的示例代码，可能需要根据实际情况进行修改和调整。