矩阵里有一个容易消去的主元,是否就能把对应的两条链删掉?单独看一张边界矩阵,这像一次普通消元;放回整个复形,却还要与上下两张边界相容。更重要的是,算出一个更小的同调群并不够:如果原来有一条具体循环,怎样把它送进小复形,又怎样恢复回来?本页把这些映射一起写出。
形式陈述
输入不仅有矩阵,还带次数与系数环
设 是交换含幺环, 是有限有基自由模理路自由模Free module具有基、因而每个元素可唯一写为有限线性组合的模。组成的链复形,边界 降低一次数且满足 。选基向量 、,并记
这里 表示 的 坐标。要求 是 的单位,即确有 。域上非零即可;整数上只有 可以这样直接消去。
令 的自由基由原基删除 得到。记 为剩余坐标的包含, 为删除这两项的坐标投影。它们此时只是分次模映射,还未声称是链映射。定义次数升一的映射
以及
这份构造满足
所以 是明确的链同伦逆理路链同伦Chain homotopy用升高一次数的态射族见证两个链映射之差沿边界方向可消去的关系。,而不只是两份同调维数碰巧相等。式(3)额外的三个零式让后续复合和验证更直接。
过滤信息需要另一项检查
若每个基向量 有实数级别 ,并以级别至多 的基张成子复形 ,就要求边界不提高级别。本页的收缩保持全部这些子复形,只需再要求
此时 、、 都不提高级别,式(2)–(3)在每个 上同时成立。因此它们还保留整个持续同调理路有限过滤的持续同调Persistent homology · 持久同调 · 持续同调从有限子复形过滤构造带包含诱导映射的同调序列,用区间记录类的出生与死亡,并区分逐层维数、系数选择、同级事件及最终常值延拓。对象中的跨时刻映射。只有最终总复形可缩,不能替代这项过滤检查。
直觉
删除的是一对可互相解释的方向
写成
其中 不含 。把低一度的基向量 换成 ,它与 就组成边界为单位的两项可缩复形。其他 次链的边界如果也经过 ,要先减掉适当倍数的 ,才能与这对方向分开。
具体地,对剩余 次基向量 ,
在其他次数, 就是保留原基。反向映射则满足
并把其他基送到其剩余坐标。因此 不是简单扔掉,而是被它与 的关系替换。
对剩余 、,新的边界系数为
它与Schur消元理路Schur 补与静态凝聚Schur complement · Static condensation把内部变量的精确响应折入界面矩阵与右端,推导 Schur 补、恢复公式及最小能量性质,并逐项凝聚五节点链。的“原系数减去经主元的间接通路”同形。不过这里的交付还包含相邻边界与链同伦,不能把线性系统的一次行列删除直接当成全部证明。
为什么这些公式真的组成一个复形
由定义,;而 ,因为 的 系数为一。于是 和 都是幂等映射,二者乘积均为零:一边使用 ,另一边使用 。因此
对剩余基, 只可能额外添上 的倍数,故 。还须核 ,不能只从这一条左逆猜出右逆。事实上 、,且每个 的 坐标为零,因为 。故 只能是 的倍数,再用 与 ,便得所需恒等式。
现在 、。 与 交换说明 是子复形,而 把 的坐标同构到该子复形。把 通过这个同构搬回去,就得到式(1)的 及全部式(2)。最后 给式(3)的旁条件。
同样由 ,可把微分公式写成
式(6)只是其中发生非平凡系数更新的一块;上下相邻边界应一起按完整公式保留,不能各自挑一套互不相容的基。
例子与边界
一条有两个余面的边,仍能代数消去
取有向三角形 、,共用对角边 。约定顶点递增给定方向,则
选 、,主元 。删除这对后,剩余二维基记作 ,它的边界变成
这是一条四边环。映射具体为
例如 ,公共对角边正好抵消。又如
边 在完整复形中有两个三角余面,并不是可作单纯自由面塌缩的自由面。这里合法的是代数单位对分离;剩余链复形也不再只是原来若干单形的子复形。
非零主元不一定能除
整数复形 的零次同调是 。若把主元二当作单位删掉,剩余零复形的同调为零,显然丢失了挠信息。改到有理数域后,二可逆,这对确实可缩;改到特征二的域,边界却变成零,两条基各自留下一个同调方向。
因此必须先声明系数环,再判断主元。反过来,“当前没有单位条目”也不意味着复形已达到某种最小形式;换基或成块处理可能暴露新的可逆部分。本页只认证已给单位对,不声称找到了最优删除顺序。
异级消去可能抹掉一段真实寿命
令 , 在级别零出现、 在级别一出现。最终两项复形可缩;但在 , 是尚未被填掉的同调类,给区间 。如果把二者直接删除,早期同调就消失了。
在构造中,收缩把 送到 ,提高了过滤级别,所以失败原因可直接在 上看见。相等条件(4)保证 不提高级别,再由边界同样不提高级别推出 的过滤相容。这正是额外条件的用途。
推论与应用
多步消去要把证书一起复合
设先从 收缩到 ,得到 ;再从 收缩到 ,得到 。合成证书为
将 拆成
即可得到 。旁条件由两步各自的 、、 逐项相消。每一步都保过滤时,合成也保过滤。
交付一个小复形之后,应同时保留这些映射。原循环 被送到 ;若在小复形找到 ,则
因为 。这给出原复形里真正的填充链,而不是只报告“该类为零”。
作为持续计算的预处理
持续边界约化理路持续边界矩阵的列约化证书Persistence boundary matrix reduction · Persistent homology matrix reduction按过滤顺序约化边界列,维护R=DV并构造真实过滤基W与DW=WB,严格恢复条码、出生循环和死亡填充链,核清同级空条与R不是新微分的边界。需要保留跨级别的映射。可以先删除同级单位对,再在剩余复形上计算条码,最后用 恢复循环与填充链。重新排序剩余基时,应继续按过滤值排序,同级先低次数,保证所有非零边界仍在列之前。
一次单位对更新在相邻边界上可能产生新非零系数,称为填充。更小的基数不保证更稀疏的矩阵,选择顺序还影响中间存储。本页的数学保证是映射与同调保真;没有把它包装成任何输入上都更快的算法。
参考资料