% 指定区域大小
width = 10;
height = 10;
% 生成晶核点
numPoints = 200; % 晶核点数量
x = rand(numPoints, 1) * width;
y = rand(numPoints, 1) * height;
% 绘制晶核点
scatter(x, y, 'filled');
axis([0 width 0 height]);
% 进行Delaunay三角剖分
dt = DelaunayTri(x, y);
% 绘制三角形
triplot(dt);
axis([0 width 0 height]);
% 计算Voronoi镶嵌
[vx, vy] = voronoi(dt);
% 绘制Voronoi镶嵌
hold on;
plot(vx, vy, 'r');
axis([0 width 0 height]);