Skip to content

方法Method

有基链复形的单位主元消去

Elementary chain cancellation · Unit-pivot chain contraction · 代数链消去

在交换环有基自由链复形中删除一个单位边界对,显式交出约化微分、来回链映射及收缩同伦,并核清过滤保持、复合和非单位主元边界。

矩阵里有一个容易消去的主元,是否就能把对应的两条链删掉?单独看一张边界矩阵,这像一次普通消元;放回整个复形,却还要与上下两张边界相容。更重要的是,算出一个更小的同调群并不够:如果原来有一条具体循环,怎样把它送进小复形,又怎样恢复回来?本页把这些映射一起写出。

形式陈述 ​

输入不仅有矩阵,还带次数与系数环 ​

设 R 是交换含幺环,C∙ 是有限有基自由模组成的链复形,边界 d 降低一次数且满足 d2=0。选基向量 a∈Cp、b∈Cp−1,并记

u=[b](da).

这里 [b](v) 表示 v 的 b 坐标。要求 u 是 R 的单位,即确有 u−1∈R。域上非零即可;整数上只有 1,−1 可以这样直接消去。

令 C′ 的自由基由原基删除 a,b 得到。记 j:C′→C 为剩余坐标的包含,π:C→C′ 为删除这两项的坐标投影。它们此时只是分次模映射,还未声称是链映射。定义次数升一的映射

h(b)=u−1a,h(x)=0(x≠b),

以及

(1)P=1−dh−hd,ι=Pj,ρ=πP,d′=πPdj.

这份构造满足

(2)(d′)2=0,dι=ιd′,ρd=d′ρ,ρι=1C′,(3)1C−ιρ=dh+hd,h2=0,hι=0,ρh=0.

所以 ι,ρ 是明确的链同伦逆,而不只是两份同调维数碰巧相等。式(3)额外的三个零式让后续复合和验证更直接。

过滤信息需要另一项检查 ​

若每个基向量 x 有实数级别 f(x),并以级别至多 t 的基张成子复形 FtC,就要求边界不提高级别。本页的收缩保持全部这些子复形,只需再要求

(4)f(a)=f(b).

此时 h、ι、ρ 都不提高级别,式(2)–(3)在每个 t 上同时成立。因此它们还保留整个持续同调对象中的跨时刻映射。只有最终总复形可缩,不能替代这项过滤检查。

直觉

删除的是一对可互相解释的方向 ​

写成

da=ub+w,

其中 w 不含 b。把低一度的基向量 b 换成 u−1da,它与 a 就组成边界为单位的两项可缩复形。其他 p 次链的边界如果也经过 b,要先减掉适当倍数的 a,才能与这对方向分开。

具体地,对剩余 p 次基向量 y,

(5)ι(y)=y−au−1[b](dy).

在其他次数,ι 就是保留原基。反向映射则满足

ρ(a)=0,ρ(b)=−u−1π(da),

并把其他基送到其剩余坐标。因此 b 不是简单扔掉,而是被它与 da 的关系替换。

对剩余 y∈Cp、x∈Cp−1,新的边界系数为

(6)[d′y]x=[dy]x−[da]xu−1[dy]b.

它与Schur消元的“原系数减去经主元的间接通路”同形。不过这里的交付还包含相邻边界与链同伦,不能把线性系统的一次行列删除直接当成全部证明。

为什么这些公式真的组成一个复形 ​

由定义,h2=0;而 hdh=h,因为 d(u−1a) 的 b 系数为一。于是 dh 和 hd 都是幂等映射,二者乘积均为零:一边使用 h2=0,另一边使用 d2=0。因此

(7)P2=P,dP=Pd,Ph=hP=0.

对剩余基,Pj 只可能额外添上 a 的倍数,故 πPj=1。还须核 PjπP=P,不能只从这一条左逆猜出右逆。事实上 Pa=0、Pda=0,且每个 Pv 的 b 坐标为零,因为 hP=0。故 Pv−jπPv 只能是 a 的倍数,再用 Pa=0 与 P2=P,便得所需恒等式。

现在 ιρ=P、ρι=1。P 与 d 交换说明 imP 是子复形,而 ι 把 C′ 的坐标同构到该子复形。把 d 通过这个同构搬回去,就得到式(1)的 d′ 及全部式(2)。最后 Ph=hP=0 给式(3)的旁条件。

同样由 d2=0,可把微分公式写成

d′=π(d−dhd)j.

式(6)只是其中发生非平凡系数更新的一块;上下相邻边界应一起按完整公式保留,不能各自挑一套互不相容的基。

例子与边界

一条有两个余面的边,仍能代数消去 ​

取有向三角形 [012]、[023],共用对角边 [02]。约定顶点递增给定方向,则

d[012]=[12]−[02]+[01],d[023]=[23]−[03]+[02].

选 a=[012]、b=[02],主元 u=−1。删除这对后,剩余二维基记作 [023]′,它的边界变成

(8)d′[023]′=[01]+[12]+[23]−[03].

这是一条四边环。映射具体为

ι([023]′)=[023]+[012],ρ([02])=[01]+[12],h([02])=−[012].

例如 ιd′[023]′=d([023]+[012]),公共对角边正好抵消。又如

(dh+hd)[02]=−[12]+[02]−[01]=[02]−ιρ[02].

边 [02] 在完整复形中有两个三角余面,并不是可作单纯自由面塌缩的自由面。这里合法的是代数单位对分离;剩余链复形也不再只是原来若干单形的子复形。

非零主元不一定能除 ​

整数复形 Za→×2Zb 的零次同调是 Z/2。若把主元二当作单位删掉,剩余零复形的同调为零,显然丢失了挠信息。改到有理数域后,二可逆,这对确实可缩;改到特征二的域,边界却变成零,两条基各自留下一个同调方向。

因此必须先声明系数环,再判断主元。反过来,“当前没有单位条目”也不意味着复形已达到某种最小形式;换基或成块处理可能暴露新的可逆部分。本页只认证已给单位对,不声称找到了最优删除顺序。

异级消去可能抹掉一段真实寿命 ​

令 db=a,a 在级别零出现、b 在级别一出现。最终两项复形可缩;但在 0≤t<1,a 是尚未被填掉的同调类,给区间 [0,1)。如果把二者直接删除,早期同调就消失了。

在构造中,收缩把 a 送到 b,提高了过滤级别,所以失败原因可直接在 h 上看见。相等条件(4)保证 h 不提高级别,再由边界同样不提高级别推出 P,ι,ρ 的过滤相容。这正是额外条件的用途。

推论与应用

多步消去要把证书一起复合 ​

设先从 C 收缩到 C′,得到 (ι1,ρ1,h1);再从 C′ 收缩到 C″,得到 (ι2,ρ2,h2)。合成证书为

(9)ι=ι1ι2,ρ=ρ2ρ1,h=h1+ι1h2ρ1.

将 1−ιρ 拆成

(1−ι1ρ1)+ι1(1−ι2ρ2)ρ1

即可得到 dh+hd。旁条件由两步各自的 hi2=0、hiιi=0、ρihi=0 逐项相消。每一步都保过滤时,合成也保过滤。

交付一个小复形之后,应同时保留这些映射。原循环 z 被送到 ρz;若在小复形找到 ρz=d′w,则

z=d(ιw+hz),

因为 dz=0。这给出原复形里真正的填充链,而不是只报告“该类为零”。

作为持续计算的预处理 ​

持续边界约化需要保留跨级别的映射。可以先删除同级单位对,再在剩余复形上计算条码,最后用 ι 恢复循环与填充链。重新排序剩余基时,应继续按过滤值排序,同级先低次数,保证所有非零边界仍在列之前。

一次单位对更新在相邻边界上可能产生新非零系数,称为填充。更小的基数不保证更稀疏的矩阵,选择顺序还影响中间存储。本页的数学保证是映射与同调保真;没有把它包装成任何输入上都更快的算法。

参考资料
  • Dmitry N. Kozlov,Discrete Morse Theory for free chain complexes,2005,§1的单位权与代数配对,§2 Theorem 2.1的逐对换基分离。本文限定单对显式收缩并自行验证过滤及复合公式,不以几何塌缩代替代数条件。
关系图谱10 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系