CCG 05-02:锥奇点与离散 Gauss-Bonnet

本页只解决一个问题

分片平坦曲面的曲率为何集中在顶点,角亏怎样满足拓扑恒等式?

本页是原长篇讲义的单概念版本。第一次阅读只追踪本页定义的对象;看见未掌握的前置符号时,先沿链接退回,不要跳过定义硬背结论。

开始前

3. 分片平坦度量与锥奇点

每个三角形内部是欧氏平面,沿边按等长方式粘合。曲率只集中在顶点,顶点邻域像一个欧氏圆锥。

内部顶点ii 的总角为

Θi=∑t∋iθit.\Theta_i=\sum_{t\ni i}\theta_i^t.

离散高斯曲率(角度亏损)为

Ki=2π−Θi.\boxed{K_i=2\pi-\Theta_i.}

边界顶点的边界曲率为

ki=π−Θi.\boxed{k_i=\pi-\Theta_i.}

这里KiK_i 是积分曲率,不是除以面积后的曲率密度。若需近似点值曲率,可除以 barycentric 或 Voronoi dual area。

4. 离散 Gauss–Bonnet 的推导

设网格有边界。将所有顶点曲率求和:

∑i∈V∘(2π−Θi)+∑i∈∂V(π−Θi).\sum_{i\in V^\circ}(2\pi-\Theta_i) +\sum_{i\in\partial V}(\pi-\Theta_i).

每个三角形内角和为π\pi,所以

∑iΘi=πF.\sum_i\Theta_i=\pi F.

记内部顶点数ViV_i、边界顶点数VbV_b:

∑Ki=2πVi+πVb−πF.\sum K_i=2\pi V_i+\pi V_b-\pi F.

对三角网格,数面-边关联:

3F=2Ei+Eb,3F=2E_i+E_b,

每条边界分支是闭合多边形,所以Eb=VbE_b=V_b。代数化简得到

∑i∈V∘Ki+∑i∈∂Vki=2π(V−E+F)=2πχ(M).\boxed{ \sum_{i\in V^\circ}K_i+ \sum_{i\in\partial V}k_i =2\pi(V-E+F)=2\pi\chi(M). }

离散 Gauss–Bonnet 对任意三角剖分严格成立,不是收敛意义下的近似。

读完必须能回答

  1. 不看公式,用自己的话回答开头的核心问题。
  2. 指出本页最重要公式中每个符号的类型:标量、向量、矩阵、映射还是集合。
  3. 找出定理成立所需的假设,并举一个删除某条假设后可能失败的例子。
  4. 对最小例子手算一遍;若本页含代码,再比较手算与程序输出。

前后课程

上一课:三角网格与边长度量 | 下一课:有限元基函数与余切权重