用maltab写一段程序实现，目前有4段圆弧的4个散点的矩阵，4个矩阵都为[1：n],要从4个矩阵中，每个矩阵找一个点，组成一个矩形，判断条件为对角线互相平分且相等的四边形是矩形，后续计算该矩形的面积，找出面积最大的矩形并将4个顶点坐标求出来，储存在一个数组中。

在MATLAB中，你可以按照以下步骤编写这个程序：

% 假设你已经有了四个矩阵，分别存储了四个散点的x坐标和y坐标
% A = [x1 y1; x2 y2; ... xn yn];
% B = [x1' y1'; x2' y2'; ... xn' yn'];
% C = ... (类似)
% D = ...

% 确定每个矩阵的第一列和最后一行作为可能的对角线顶点
% 对角线1的顶点1
line1A = A(1,:);
line1B = B(end,:);
% 对角线2的顶点1
line2A = A(end,:);
line2B = B(1,:);

% 计算对角线长度
d1 = sqrt(sum((line1A - line1B).^2));
d2 = sqrt(sum((line2A - d2) < eps && abs(d1 + d2) == 2 * max([d1 d2]) % eps是一个很小的数用于精度处理
    % 计算矩形面积
    area = d1 * d2;
    
    % 找出面积最大的矩形
    if ~exist('max_area', 'var') || area > max_area
        max_area = area;
        rect_points = [line1A; line1B; line2A; line2B]; % 存储顶点坐标
    end
end

% 如果找到符合条件的最大矩形
if isfield(max_area, 'value')
    fprintf('最大矩形的面积为: %.2f\n', max_area);
    fprintf('矩形的顶点坐标为:\n');
    disp(rect_points);
else
    fprintf('未找到满足条件的矩形.\n');
end

这里假设`eps`是一个非常小的数值，用于比较浮点数的精度差异。你需要根据实际情况调整`eps`的值。