voronoi多边形面积求算

clear;clc; m=100; x=rand(1,m); % 生成随机点x坐标 y=rand(1,m); % 生成随机点y坐标 figure(2); voronoi(x,y); % 绘制voronoi多边形 axis([0,1,0,1]);grid; [v,c]=voronoin([x;y]') % v是所有的voronoi边的顶点坐标，c{i}就是第i个点周围的voronoi边的顶点下标 for n=1:m disp(['点(' num2str(x(n)) ',' num2str(y(n)) ')周围的voronoi边的顶点的坐标是:']) disp([v(c{n},1),v(c{n},2)]) end for p=1:m t=[v(c{p},1),v(c{p},2)]; [r,s]=size(t); a=0; for q=2:r a1=1/2*(t(q-1,1)*t(q,2)-t(q,1)*t(q-1,2)); % S=1/2[(x1*y2-x2*y1)+(x2*y3-x3*y2)+...+(xn-1*yn-xn*yn-1)+(xn*y1-x1*yn)] a=a+a1; end a=a+1/2*(t(r,1)*t(1,2)-t(1,1)*t(r,2)); disp(['点(' num2str(x(p)) ',' num2str(y(p)) ')周围的voronoi多边形的面积是:']) a end % a3=1/2*((0.7327*0.4513-0.9894*0.4851)+(0.9894*0.7517-0.8156*0.4513)+(0.8156*0.7134-0.5210*0.7517)+(0.5210*0.6966-0.5156*0.7134)+(0.5156*0.4851-0.7327*0.6966))