如何使用 GUDHI 正确计算并显示一维同调群(H₁)

星辰大大_6251

星辰大大_6251

2026-08-15

472人浏览

原创

如何使用 GUDHI 正确计算并显示一维同调群(H₁)

gudhi 默认仅对实际存在的拓扑特征输出同调信息;原始代码中未构造闭合环路,故无非平凡的一维同调类,导致维度 1 的结果不出现。需显式构建 1-单纯形闭链(如四边形循环)才能观测到 h₁ 的持久区间。

gudhi 默认仅对实际存在的拓扑特征输出同调信息;原始代码中未构造闭合环路,故无非平凡的一维同调类,导致维度 1 的结果不出现。需显式构建 1-单纯形闭链(如四边形循环)才能观测到 h₁ 的持久区间。

在拓扑数据分析(TDA)中,GUDHI 库通过单纯复形计算持久同调(persistent homology),其输出的同调群维度严格对应数据中真实存在的拓扑结构:

  • 维度 0(H₀):反映连通分支数量(如孤立点、连通块);
  • 维度 1(H₁):对应“洞”或一维循环(如环、孔);
  • 更高维度(H₂+):对应空腔、球面等更复杂结构。

你最初的代码构建的是三个二维单纯形(三角形):[1,2,3]、[2,3,4]、[3,4,5]。它们共享边 [2,3]、[3,4],整体形成一条“链状”结构(类似锯齿形),没有闭合环路,因此不存在非平凡的一维同调类 —— GUDHI 正确地不输出任何 Dimension 1 区间,这并非 bug,而是数学上的准确结果。

要观测到 Dimension 1 的输出,必须构造一个至少由四个 1-单纯形(边)组成的闭合循环,例如四边形 [1,2] → [2,3] → [3,4] → [4,1]:

import gudhi

# ✅ 正确构造一个 1D 循环(无填充面,确保 H₁ 非零)
simplex_tree = gudhi.SimplexTree()
simplex_tree.insert([1, 2])
simplex_tree.insert([2, 3])
simplex_tree.insert([3, 4])
simplex_tree.insert([4, 1])

# 计算持久同调(默认使用 Vietoris-Rips 框架下的单纯复形)
homology = simplex_tree.persistence()

print("Homology groups:")
for dim, (birth, death) in homology:
    print(f"Dimension {dim}: birth = {birth}, death = {death}")

输出:

Homology groups:
Dimension 0: birth = 0.0, death = inf
Dimension 1: birth = 0.0, death = inf

✅ 解释:

  • Dimension 0:所有顶点连通(因边已构成单连通图),故一个无限长的连通分支;
  • Dimension 1:四条边形成一个无填充的环(即未插入面 [1,2,3,4] 或任意三角形),因此产生一个持久性为 ∞ 的一维洞。

⚠️ 关键注意事项:

  • 不要插入高维单纯形(如 [1,2,3])——它会“填充”环,使 H₁ 消失;
  • 若需多尺度分析,可配合 gudhi.RipsComplex 自动构建邻域复形,但本例中手动插入边已足够;
  • death = inf 表示该特征在复形构建完成时仍存在(无更高尺度将其填满);
  • 使用 simplex_tree.get_filtration() 可验证复形结构,确保无意外添加的面。

总结:GUDHI 的输出始终忠实反映输入复形的代数拓扑结构。观察不到 H₁ 并非库失效,而是输入数据本身不具备一维洞。构造明确的闭合 1-循环是触发 Dimension 1 输出的必要且充分条件。

PHP速学视频免费教程(入门到精通)
PHP速学视频免费教程(入门到精通)

PHP怎么学习?PHP怎么入门?PHP在哪学?PHP怎么学才快?不用担心,这里为大家提供了PHP速学教程(入门到精通),有需要的小伙伴保存下载就能学习啦!

下载

相关标签:

本站声明:本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系admin@php.cn

相关专题

更多
python打包成可执行文件
python打包成可执行文件

本专题为大家带来python打包成可执行文件相关的文章,大家可以免费的下载体验。

2023.07.20

1611

4

python能做什么
python能做什么

python能做的有:可用于开发基于控制台的应用程序、多媒体部分开发、用于开发基于Web的应用程序、使用python处理数据、系统编程等等。本专题为大家提供python相关的各种文章、以及下载和课程。

2023.07.25

3864

7

format在python中的用法
format在python中的用法

Python中的format是一种字符串格式化方法,用于将变量或值插入到字符串中的占位符位置。通过format方法,我们可以动态地构建字符串,使其包含不同值。php中文网给大家带来了相关的教程以及文章,欢迎大家前来阅读学习。

2023.07.31

1609

3

python教程
python教程

Python已成为一门网红语言,即使是在非编程开发者当中,也掀起了一股学习的热潮。本专题为大家带来python教程的相关文章,大家可以免费体验学习。

2023.08.03

22337

23

python环境变量的配置
python环境变量的配置

Python是一种流行的编程语言,被广泛用于软件开发、数据分析和科学计算等领域。在安装Python之后,我们需要配置环境变量,以便在任何位置都能够访问Python的可执行文件。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

2747

5

python eval
python eval

eval函数是Python中一个非常强大的函数,它可以将字符串作为Python代码进行执行,实现动态编程的效果。然而,由于其潜在的安全风险和性能问题,需要谨慎使用。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

2787

5

scratch和python区别
scratch和python区别

scratch和python的区别:1、scratch是一种专为初学者设计的图形化编程语言,python是一种文本编程语言;2、scratch使用的是基于积木的编程语法,python采用更加传统的文本编程语法等等。本专题为大家提供scratch和python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

1123

5

python合并两个列表
python合并两个列表

Python是一种强大的编程语言,具有许多方便的功能和工具。在Python中,有多种方法可以合并两个列表。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.10

596

4

python是前端还是后端
python是前端还是后端

Python属于前端也属于后端,其灵活性和丰富的生态系统使得开发人员能够在不同的领域中灵活运用。本专题为大家提供python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

2163

5

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
热门推荐
/
最新课程
phpStudy极速入门视频教程
phpStudy极速入门视频教程

共6课时 | 54.6万人学习

独孤九贱(4)_PHP视频教程
独孤九贱(4)_PHP视频教程

共89课时 | 133.3万人学习