求解voronoi多边形顶点坐标及面积的matalab程序

关于voronoi多边形顶点坐标及面积求算的matalab仿真程序

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.7517）+（0.5210*0.6966-0.5156*0.7134）+（0.5156*0.4851-0.7327*0.6966））

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

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

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

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

                 1324                    2