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??matalab仿真程序_voronoi多邊形頂點坐標及面積求算\voronoi面積.m

?????文件????????390??2014-06-05?16:57??matalab仿真程序_voronoi多邊形頂點坐標及面積求算\voronoi頂點坐標.m

?????目錄??????????0??2019-03-26?15:11??matalab仿真程序_voronoi多邊形頂點坐標及面積求算

-----------?---------??----------?-----??----

?????????????????1324????????????????????3