自动关联目录:黎曼面、Abel 微分与参数域
本章目标 圆盘型曲面可直接映到平面区域,但闭合环面或高亏格曲面不可能用一个无周期的平面坐标表示。本章回答:
全纯函数为什么在紧致黎曼面上过于刚性? 为什么要研究全纯/亚纯 1-形式? 2 g 2g 2 g 条同调基与g g g 维全纯 1-形式如何组成周期矩阵?Abel–Jacobi 映射为何把曲面送入复环面? Riemann–Roch 如何计数允许给定零极点的函数/微分? 这些理论怎样变成高亏格网格参数化算法? 1. 黎曼面上的全纯对象 黎曼面M M M 带有复坐标图( U α , z α ) (U_\alpha,z_\alpha) ( U α , z α ) ,过渡函数全纯。函数f : M → C f:M\to\mathbb C f : M → C 全纯,是指每个局部表示
f ∘ z α − 1 f\circ z_\alpha^{-1} f ∘ z α − 1
在平面复分析意义下全纯。
1.1 紧致性带来的刚性 若M M M 紧致且连通,任意全纯函数f : M → C f:M\to\mathbb C f : M → C 都是常数。
证明:∣ f ∣ |f| ∣ f ∣ 在紧致空间上取得最大值;最大模原理说明非常数全纯函数不能在内部取最大值,而紧致无边界曲面的每一点都是内部点。
所以要在紧致曲面上得到非平凡对象,必须允许:
函数取值到 Riemann 球面,即允许极点(亚纯函数); 研究随坐标变换带权变化的微分形式。 2. 亚纯函数、零点与极点 在点p p p 的局部坐标z ( p ) = 0 z(p)=0 z ( p ) = 0 中,亚纯函数有 Laurent 展开
f ( z ) = z m ( a 0 + a 1 z + ⋯ ) , a 0 ≠ 0 , f(z)=z^m(a_0+a_1z+\cdots), \qquad a_0\ne0, f ( z ) = z m ( a 0 + a 1 z + ⋯ ) , a 0 = 0 ,
其中:
m > 0 m>0 m > 0 :p p p 是m m m 阶零点;m = 0 m=0 m = 0 :既非零点也非极点;m < 0 m<0 m < 0 :p p p 是− m -m − m 阶极点。记阶数为ord p ( f ) = m \operatorname{ord}_p(f)=m ord p ( f ) = m 。
紧致黎曼面上的亚纯函数满足
∑ p ∈ M ord p ( f ) = 0 , \sum_{p\in M}\operatorname{ord}_p(f)=0, p ∈ M ∑ ord p ( f ) = 0 ,
即零点总阶数等于极点总阶数。
3. 除子 除子是点的有限整数线性组合:
D = ∑ p ∈ M n p p . D=\sum_{p\in M}n_pp. D = p ∈ M ∑ n p p .
次数:
deg D = ∑ p n p . \deg D=\sum_pn_p. deg D = p ∑ n p .
亚纯函数的主除子
( f ) = ∑ p ord p ( f ) p . (f)=\sum_p\operatorname{ord}_p(f)p. ( f ) = p ∑ ord p ( f ) p .
若D = ∑ n p p D=\sum n_pp D = ∑ n p p ,定义
L ( D ) = { f 亚纯 : ( f ) + D ≥ 0 } ∪ { 0 } . L(D)=\{f\text{ 亚纯}: (f)+D\ge0\}\cup\{0\}. L ( D ) = { f 亚纯 : ( f ) + D ≥ 0 } ∪ { 0 } .
它表示极点位置和最高阶数由D D D 限制的函数空间,维数记为
ℓ ( D ) = dim L ( D ) . \ell(D)=\dim L(D). ℓ ( D ) = dim L ( D ) .
4. 全纯 1-形式 在局部复坐标中,全纯 1-形式写成
ω = f α ( z α ) d z α , \omega=f_\alpha(z_\alpha)dz_\alpha, ω = f α ( z α ) d z α ,
其中f α f_\alpha f α 全纯。换坐标z β = z β ( z α ) z_\beta=z_\beta(z_\alpha) z β = z β ( z α ) 时
f α ( z α ) = f β ( z β ) d z β d z α . f_\alpha(z_\alpha) =f_\beta(z_\beta) \frac{dz_\beta}{dz_\alpha}. f α ( z α ) = f β ( z β ) d z α d z β .
所以ω \omega ω 作为整体不依赖坐标。
全纯 1-形式空间记为H 0 ( M , K ) H^0(M,K) H 0 ( M , K ) ,其中K K K 是 canonical bundle。对亏格g g g 的紧致黎曼面:
dim C H 0 ( M , K ) = g . \boxed{\dim_{\mathbb C}H^0(M,K)=g.} dim C H 0 ( M , K ) = g .
4.1 与调和 1-形式的关系 若α \alpha α 是实调和 1-形式,在二维中定义
ω = α + i ⋆ α . \omega=\alpha+i\star\alpha. ω = α + i ⋆ α .
则ω \omega ω 是全纯 1-形式(符号依复结构约定可能为α − i ⋆ α \alpha-i\star\alpha α − i ⋆ α )。反过来,全纯 1-形式的实部和虚部都是调和 1-形式。
实调和空间维数2 g 2g 2 g ,复全纯 1-形式空间维数g g g ,两者完全吻合。
5. 零点与 canonical divisor 非零全纯 1-形式的零点总阶数为
∑ p ord p ( ω ) = 2 g − 2. \sum_p\operatorname{ord}_p(\omega)=2g-2. p ∑ ord p ( ω ) = 2 g − 2.
这个除子称为 canonical divisorK K K ,其次数
deg K = 2 g − 2. \deg K=2g-2. deg K = 2 g − 2.
例子:
球面g = 0 g=0 g = 0 :不存在非零全纯 1-形式; 环面g = 1 g=1 g = 1 :d z dz d z 无零点; g ≥ 2 g\ge2 g ≥ 2 :任意全纯 1-形式必须有总阶数2 g − 2 2g-2 2 g − 2 的零点。这些零点会成为参数化中的锥奇点和 foliation 的临界点。
6. canonical homology basis 选取2 g 2g 2 g 条闭曲线
a 1 , b 1 , … , a g , b g a_1,b_1,\ldots,a_g,b_g a 1 , b 1 , … , a g , b g
构成H 1 ( M ; Z ) H_1(M;\mathbb Z) H 1 ( M ; Z ) 的辛基,使交数满足
a i ⋅ a j = 0 , b i ⋅ b j = 0 , a i ⋅ b j = δ i j . a_i\cdot a_j=0, \qquad b_i\cdot b_j=0, \qquad a_i\cdot b_j=\delta_{ij}. a i ⋅ a j = 0 , b i ⋅ b j = 0 , a i ⋅ b j = δ ij .
交数是带符号的横截相交点数量。这个基不是唯一的,不同基之间由整数辛群S p ( 2 g , Z ) Sp(2g,\mathbb Z) S p ( 2 g , Z ) 作用联系。
沿全部a i , b i a_i,b_i a i , b i 切开,曲面变成一个有4 g 4g 4 g 条边的基本多边形。
7. 周期与归一化 Abel 微分 取全纯 1-形式基ω 1 , … , ω g \omega_1,\ldots,\omega_g ω 1 , … , ω g 。定义a a a 周期和b b b 周期:
A i j = ∮ a j ω i , B i j = ∮ b j ω i . A_{ij}=\oint_{a_j}\omega_i, \qquad B_{ij}=\oint_{b_j}\omega_i. A ij = ∮ a j ω i , B ij = ∮ b j ω i .
可线性组合使a a a 周期归一化为
∮ a j ω i = δ i j . \oint_{a_j}\omega_i=\delta_{ij}. ∮ a j ω i = δ ij .
此时b b b 周期组成 Riemann period matrix
τ i j = ∮ b j ω i . \tau_{ij}=\oint_{b_j}\omega_i. τ ij = ∮ b j ω i .
Riemann 双线性关系推出:
τ T = τ , Im τ > 0. \tau^T=\tau, \qquad \operatorname{Im}\tau>0. τ T = τ , Im τ > 0.
所以τ \tau τ 位于 Siegel 上半空间。
8. Riemann 双线性关系 对闭 1-形式α , β \alpha,\beta α , β :
∫ M α ∧ β = ∑ i = 1 g ( ∮ a i α ∮ b i β − ∮ b i α ∮ a i β ) . \int_M\alpha\wedge\beta =\sum_{i=1}^g \left( \oint_{a_i}\alpha\oint_{b_i}\beta -\oint_{b_i}\alpha\oint_{a_i}\beta \right). ∫ M α ∧ β = i = 1 ∑ g ( ∮ a i α ∮ b i β − ∮ b i α ∮ a i β ) .
证明思想:沿 canonical basis 切开M M M 得基本多边形;在单连通多边形上写α = d f \alpha=df α = df ;应用 Stokes 定理把面积积分化成边界积分,再利用成对边上的函数跳跃等于周期。
这一定理把曲面上的积分内积完全转化为有限个周期。
9. Jacobian variety 与 Abel–Jacobi 映射 定义格
Λ = Z g + τ Z g . \Lambda=\mathbb Z^g+\tau\mathbb Z^g. Λ = Z g + τ Z g .
复环面
Jac ( M ) = C g / Λ \operatorname{Jac}(M)=\mathbb C^g/\Lambda Jac ( M ) = C g /Λ
称为M M M 的 Jacobian variety。
固定基点p 0 p_0 p 0 ,Abel–Jacobi 映射
A : M → Jac ( M ) , \mathcal A:M\to\operatorname{Jac}(M), A : M → Jac ( M ) ,
A ( p ) = ( ∫ p 0 p ω 1 , l d o t s , ∫ p 0 p ω g ) m o d Λ . \mathcal A(p)= \left( \int_{p_0}^{p}\omega_1,ldots, \int_{p_0}^{p}\omega_g \right)\mod\Lambda. A ( p ) = ( ∫ p 0 p ω 1 , l d o t s , ∫ p 0 p ω g ) mod Λ.
路径改变会增加某个周期,因此模掉Λ \Lambda Λ 后定义良好。
9.1 环面情形 g = 1 g=1 g = 1 时,M ≅ C / ( Z + τ Z ) M\cong\mathbb C/(\mathbb Z+\tau\mathbb Z) M ≅ C / ( Z + τ Z ) ,Abel–Jacobi 映射本身给出环面的平坦周期坐标。
g > 1 g>1 g > 1 时,曲面嵌入g g g 维复环面,而不是直接等同于 Jacobian。
10. Abel 定理 次数为 0 的除子
D = ∑ i n i p i , q q u a d ∑ i n i = 0 D=\sum_i n_ip_i,qquad\sum_i n_i=0 D = i ∑ n i p i , q q u a d i ∑ n i = 0
是某个亚纯函数的主除子,当且仅当它的 Abel–Jacobi 像为 0:
∑ i n i A ( p i ) = 0 in Jac ( M ) . \sum_i n_i\mathcal A(p_i)=0 \quad\text{in }\operatorname{Jac}(M). i ∑ n i A ( p i ) = 0 in Jac ( M ) .
这把“是否存在具有指定零极点的亚纯函数”转化为周期积分条件。
11. Riemann–Roch 定理 设K K K 为 canonical divisor。Riemann–Roch 定理:
ℓ ( D ) − ℓ ( K − D ) = deg D − g + 1. \boxed{ \ell(D)-\ell(K-D)=\deg D-g+1. } ℓ ( D ) − ℓ ( K − D ) = deg D − g + 1.
11.1 如何读这个公式 ℓ ( D ) \ell(D) ℓ ( D ) :允许D D D 所规定极点的亚纯函数自由度;ℓ ( K − D ) \ell(K-D) ℓ ( K − D ) :对偶障碍空间;deg D − g + 1 \deg D-g+1 deg D − g + 1 :拓扑给出的预期维数。若deg D > 2 g − 2 \deg D>2g-2 deg D > 2 g − 2 ,则ℓ ( K − D ) = 0 \ell(K-D)=0 ℓ ( K − D ) = 0 ,所以
ℓ ( D ) = deg D − g + 1. \ell(D)=\deg D-g+1. ℓ ( D ) = deg D − g + 1.
11.2 三个检验 D = 0 D=0 D = 0 :ℓ ( 0 ) = 1 \ell(0)=1 ℓ ( 0 ) = 1 (常数函数),于是1 − ℓ ( K ) = 1 − g 1-\ell(K)=1-g 1 − ℓ ( K ) = 1 − g ,得到ℓ ( K ) = g \ell(K)=g ℓ ( K ) = g 。球面g = 0 g=0 g = 0 :次数d ≥ 0 d\ge0 d ≥ 0 的除子给ℓ ( D ) = d + 1 \ell(D)=d+1 ℓ ( D ) = d + 1 ,对应有理函数多项式自由度。 canonical divisor:deg K = 2 g − 2 \deg K=2g-2 deg K = 2 g − 2 。 12. Abel 微分的三类 经典上按奇点分类:
第一类:无极点,即全纯 1-形式; 第二类:有极点但所有留数为 0; 第三类:有简单极点,留数非零但总和为 0。 它们分别适合:
构造周期平坦参数; 构造带高阶奇点的几何场; 构造源汇型参数、圆域与边界映射。 13. 二次微分与 foliation 全纯二次微分局部写成
q = φ ( z ) d z 2 . q=\varphi(z)dz^2. q = φ ( z ) d z 2 .
它定义两族方向场:
水平轨线:q ( γ ˙ 2 ) > 0 q(\dot\gamma^2)>0 q ( γ ˙ 2 ) > 0 ; 垂直轨线:q ( γ ˙ 2 ) < 0 q(\dot\gamma^2)<0 q ( γ ˙ 2 ) < 0 。 局部自然坐标
w = ∫ q w=\int\sqrt{q} w = ∫ q
中,q = d w 2 q=dw^2 q = d w 2 ,水平/垂直轨线变成水平/垂直直线。
若q q q 在p p p 有m m m 阶零点,则有m + 2 m+2 m + 2 条水平临界轨线从p p p 发出,对应锥角
( m + 2 ) π . (m+2)\pi. ( m + 2 ) π .
临界轨线图可把曲面分割成圆柱或多边形区域。这是 quad mesh、各向异性网格生成和 surface foliation 的理论来源。
14. 离散算法:从同调基到全纯微分 典型计算流程:
flowchart TD
A[三角网格] --> B[计算 canonical homology basis]
B --> C[为每个上同调类构造闭 1-形式]
C --> D[能量最小化得到调和代表]
D --> E[Hodge 星构造共轭形式]
E --> F[组合成全纯 1-形式]
F --> G[计算 a/b 周期并归一化]
G --> H[切开网格并积分得到参数域] 14.1 构造闭 1-形式 对一条 cut graph 边或对偶循环,可构造具有指定 period 的初始闭 1-上链η \eta η 。它通常不调和。
14.2 求调和代表 在同一上同调类中找最小能量代表:
min f ∥ η + d f ∥ ⋆ 1 2 . \min_f\|\eta+df\|_{\star_1}^2. f min ∥ η + df ∥ ⋆ 1 2 .
法方程:
d 0 T ⋆ 1 d 0 f = − d 0 T ⋆ 1 η . d_0^T\star_1d_0f =-d_0^T\star_1\eta. d 0 T ⋆ 1 d 0 f = − d 0 T ⋆ 1 η .
得到
α = η + d f , \alpha=\eta+df, α = η + df ,
它闭且余闭,是该上同调类唯一调和代表。
14.3 构造全纯形式 连续理论给ω = α + i ⋆ α \omega=\alpha+i\star\alpha ω = α + i ⋆ α 。离散 Hodge 星未必严格满足⋆ 2 = − 1 \star^2=-1 ⋆ 2 = − 1 ,实际算法常在调和空间内求最接近的复结构算子,或同时求调和共轭并做周期归一化。
14.4 积分 沿切开后的单连通网格,从基点开始累计边 1-形式:
z j − z i = ω i j . z_j-z_i=\omega_{ij}. z j − z i = ω ij .
闭性确保同一单连通域内路径无关;跨越切缝时坐标相差一个周期。
15. 数值验证 一个全局微分算法至少检查:
闭性 :∥ d 1 α ∥ \|d_1\alpha\| ∥ d 1 α ∥ ;余闭性 :∥ d 0 T ⋆ 1 α ∥ \|d_0^T\star_1\alpha\| ∥ d 0 T ⋆ 1 α ∥ ;周期 :沿a i , b i a_i,b_i a i , b i 的离散积分;Riemann 矩阵 :∥ τ − τ T ∥ \|\tau-\tau^T\| ∥ τ − τ T ∥ 与Im τ \operatorname{Im}\tau Im τ 最小特征值;积分一致性 :单连通切域中不同路径积分差;切缝粘合 :成对边坐标差是否为常周期。16. 与网格生成的关系 全纯 1-形式或二次微分提供方向场和自然坐标:
Re ω , Im ω \operatorname{Re}\omega,\operatorname{Im}\omega Re ω , Im ω 给局部正交坐标;零点产生必要奇点; 轨线给出切图和 patch 边界; 在每个 patch 中生成规则三角/四边形网格; 按周期与边界对应粘合回原曲面。 拓扑决定奇点总指数,无法通过优化完全消除。
17. 习题 计算题 证明环面上d z dz d z 没有零点,并计算其两个基本周期。 用 Riemann–Roch 计算球面上最多在无穷远有d d d 阶极点的函数空间维数。 对q = z m d z 2 q=z^mdz^2 q = z m d z 2 ,分析原点附近水平轨线数。 验证g = 2 g=2 g = 2 时全纯 1-形式零点总阶数为 2。 证明题 用最大模原理证明紧致黎曼面上的全纯函数为常数。 解释为何 Abel–Jacobi 积分模周期格后与路径无关。 从 Riemann–Roch 推出dim H 0 ( M , K ) = g \dim H^0(M,K)=g dim H 0 ( M , K ) = g 。 写出 Riemann 双线性关系的切多边形证明框架。 实现题 在离散环面上构造两个上同调类并求调和代表。 计算 period matrix,并检查对称性与虚部正定性。 沿同调基切开环面,积分全纯 1-形式得到平行四边形参数域。 在 genus-2 网格上绘制一个调和 1-形式的方向场和零点。 18. 本章结论 高亏格参数化的关键不是寻找一张“更好的 UV 图”,而是承认坐标天然具有周期:
同调基 → 调和 1-形式 → 全纯微分 → 周期矩阵 → 多值/周期参数域 . \boxed{ \text{同调基} \to \text{调和 1-形式} \to \text{全纯微分} \to \text{周期矩阵} \to \text{多值/周期参数域}. } 同调基 → 调和 1- 形式 → 全纯微分 → 周期矩阵 → 多值 / 周期参数域 .
下一篇将采用另一条路线:直接改变离散度量,使顶点曲率达到指定值,从而得到球面、平面或双曲标准结构:离散 Ricci 流与双曲结构 。