可以给一个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];

% 将矩形区域分割成若干个小三角形
tri = delaunay([x1, x2, x3, x4], [y1, y2, y3, y4]);

% 对每个三角形进行划分，最终得到小矩形单元
n_tri = size(tri, 1);
n_rect = 0;
for i = 1:n_tri
    % 获取当前三角形的三个顶点坐标
    x_tri = [x1, x2, x3, x4];
    y_tri = [y1, y2, y3, y4];
    x_tri(tri(i, :)) = [];
    y_tri(tri(i, :)) = [];
    
    % 判断当前三角形是否与障碍物相交
    if polyxpoly(obstacle(:,1), obstacle(:,2), x_tri, y_tri)
        % 如果相交，则将当前三角形分割成若干个小矩形单元
        [x_rect, y_rect] = polybool('subtraction', x_tri, y_tri, obstacle(:,1), obstacle(:,2));
        n_rect = n_rect + length(x_rect);
    else
        % 如果不相交，则当前三角形为一个矩形单元
        n_rect = n_rect + 1;
    end
end

disp(['矩形区域中共有 ', num2str(n_rect), ' 个小矩形单元。']);

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