基于matlab的voronoi图绘制补充之输出各个多边形面积

clc;
clear al;
close all;
N=30;%定义随机点的数量
L=20;%定义正方形边框边长
X = L*rand(1,N);
Y = L*rand(1,N); %随机生成N个2D点 
Options.plot=1; %设置1表示画出维诺图 
v=[0 0;0 L;L L;0 L;L 0];
P = polytope(v); %生成边界 
Options.pbound=P;
%axis square; 
Pn=mpt_voronoi([X' Y'],Options);

A=zeros(1,N);
for j=1:N
V= extreme(Pn(j)); %这里的V就是第j个多边形的顶点序列 
[k,A(j)] = convhull(V(:,1),V(:,2));%计算第j个多边形的面积A（j）
end
S=sum(sum(A));%%加和所有多边形的面积来验证计算的正确性

源码如上，A矩阵即为各个多边形面积。