用processing在一个正方体中绘制泰森多面体

以下是用Processing绘制泰森多面体的示例代码：

int numPoints = 200; // 点的数量
float[] x = new float[numPoints];
float[] y = new float[numPoints];
float[] z = new float[numPoints];

void setup() {
  size(500, 500, P3D);
  for (int i = 0; i < numPoints; i++) {
    x[i] = random(-width/2, width/2);
    y[i] = random(-height/2, height/2);
    z[i] = random(-width/2, width/2);
  }
}

void draw() {
  background(255);
  stroke(0);
  noFill();
  translate(width/2, height/2, 0);
  rotateX(frameCount * 0.01);
  rotateY(frameCount * 0.02);
  rotateZ(frameCount * 0.03);
  float[][] delaunayEdges = getDelaunayEdges();
  for (int i = 0; i < delaunayEdges.length; i++) {
    float[] p1 = delaunayEdges[i];
    for (int j = i+1; j < delaunayEdges.length; j++) {
      float[] p2 = delaunayEdges[j];
      if (isNeighbor(p1, p2)) {
        stroke(0, 0, 255);
        line(p1[0], p1[1], p1[2], p2[0], p2[1], p2[2]);
      }
    }
  }
}

// 获取Delaunay三角形网格的边界
float[][] getDelaunayEdges() {
  Delaunay d = new Delaunay(x, y, z);
  d.trim(0, width/2, -height/2, height/2, -width/2, width/2);
  float[][] delaunayEdges = new float[d.getEdgeCount()][6];
  int i = 0;
  for (Delaunay.Edge e : d.getEdges()) {
    float[] p1 = {x[e.a], y[e.a], z[e.a]};
    float[] p2 = {x[e.b], y[e.b], z[e.b]};
    delaunayEdges[i++] = concat(p1, p2);
  }
  return delaunayEdges;
}

// 判断两个点是否相邻
boolean isNeighbor(float[] p1, float[] p2) {
  int count = 0;
  for (int i = 0; i < numPoints; i++) {
    if (dist(p1[0], p1[1], p1[2], x[i], y[i], z[i]) < dist(p1[0], p1[1], p1[2], p2[0], p2[1], p2[2])) {
      count++;
    }
  }
  return count == 2;
}

这个程序生成一个包含200个随机点的三维Delaunay三角形网格，并在正方体内绘制与之相邻的边缘线。您可以尝试更改点的数量和正方体的大小和位置，以获得不同的效果。