MATLAB 绘制二维平面内的voronoi图

绘制二维平面内的voronoi图，并给出多边形顶点坐标及面积

clear;clc;
m=100;
x=rand（1m）;   % 生成随机点x坐标
y=rand（1m）;   % 生成随机点y坐标
figure（2）;
voronoi（xy）;   % 绘制voronoi多边形
axis（[0101]）;grid;
[vc]=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）];
    [rs]=size（t）;
    a=0;
    for q=2:r
        a1=1/2*（t（q-11）*t（q2）-t（q1）*t（q-12））;   % 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（r1）*t（12）-t（11）*t（r2））;
    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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        398  2014-06-09 15:16  voronoi local.rar

     文件        934  2014-06-05 18:05  voronoi多边形面积.m

     文件        390  2014-06-05 16:57  voronoi多边形顶点坐标.m

----------- ---------  ---------- -----  ----

                 1722                    3