如果您觉得以下内容比较生疏,不必过于焦虑,请继续关注本公众号。我们将会详尽解释以下所涉及的概念,定理,算法和应用。在未来岁月中,让我们共同学习,共同成长。

计算共形几何是计算机科学和纯粹数学之间的交叉学科,其目的是将现代几何,经典几何的概念和定理转化为计算机算法,为工程实践服务;同时,用组合离散的方法给出经典几何定理的构造性证明,换言之,用计算机算法语言从头构建经典的共形几何理论,从而绕过了历史上传统的逻辑链条。

用组合离散方法重建光滑流形上的几何结构,许多数学工具无法直接使用,例如经典复变函数论的所有方法,曲面微分几何的方法,因此其理论发展困难重重。但是,绝大多数深刻的几何定理能够被推广到离散情形,其内在的思路和直觉是一脉相承的。许多问题的连续证明非常简单,但是离散证明却异常困难,例如刘维尔定理:曲率非负的空间不存在正的非常数调和函数;许多定理的离散证明直接了当,但其连续证明却曲折晦涩,例如高斯-博纳定理:曲面全曲率被拓扑控制。在连续和离散的范畴之间自由游走,自然会领略到独到的风景。

计算机是人脑的延长物,一切抽象的概念和理论一旦被计算机科学所内化,就会被芸芸众生所感知,所体验,从而真正转化为物质力量,用于改造自然和人类社会。例如共形结构的概念,在历史上只存在于职业数学家的脑海里面,他们坚信这一自然界的结构,用理性的思维去摸索探查,虽然情深意长,心醉神迷,但是无法将其感受直白地传达给其他人。一旦共形结构的算法被发明出来,所有人都可以用计算机来观察体验这一自然结构,其实在性和普适性立刻变得无可辩驳而浅显自明。而这,正是人类思想进步的必由之路。

计算共形几何的核心是共形结构。具有共形结构的曲面是黎曼面。共形结构的研究方法有很多种,复分析,代数几何和微分几何。我们所采用的微分几何的途径,换言之,通过黎曼度量来研究共形结构。可定向曲面上的黎曼度量天然地决定了一个共形结构,其决定方式由等温坐标给出(陈省身证明)。等温坐标的存在性由拟共形映射的存在性所保证(可测黎曼映照定理)。可测黎曼映照定理的证明由辅助度量下的Ricci曲率流给出。由此,所有真实生活中的曲面都具有三维欧式空间所诱导的黎曼度量,因此都具有天然的共形结构。共形几何的普适性由此得以保障。

度量曲面上共形结构的计算方法分为三大类:调和映照,全纯微分和Ricci曲率流。

调和映照和共形映射之间的关系微妙。共形映射必然调和,反之则不尽然。拓扑球面间的调和映射必为共形映射,这是因为球面上的全纯二次微分必为零。给定高亏格曲面,共形映射未必存在,若像曲面曲率处处为负,则拓扑度为一的调和映射必为微分同胚(丘成桐定理)。调和映照的具体算法是非线性热流。

全纯微分的计算方法是基于Hodge理论,本质上是说黎曼流形上椭圆形偏微分方程的解空间的维数被流形的拓扑所控制。一个全纯微分由两个相互共轭的调和微分形式所组成。每一个deRham上同调类中,调和微分形式存在并且唯一。由于所有计算都是在离散流形上展开,光滑流形上的deRham上同调被离散流形上的单纯上同调所逼近。上同调群和下同调群同构,因此我们通过计算下同调群的基底来构造对偶上同调基底。经典代数拓扑中,单纯下同调群的计算方法是代数方法,亦即整数矩阵的标准型(Smith Norm)。但是整数矩阵的阶正比于离散流形的单形个数,动辄数万,现实中无法实用。因此,我们采用组合方法将单纯复形转化为CW复形。组合方法实际上给了我们更多,同伦群对于流形拓扑的刻画更为精密。但是,同伦群的不可交换性使得许多计算拓扑问题成为不可解问题。由同伦群引发的覆盖空间概念,特别是万有覆盖空间可以极大地简化底空间上的拓扑问题。底空间中的圈被提升为万有覆盖空间中的路径,底空间的同伦群被转化为万有覆盖空间的覆盖变换群。全纯微分作用于同伦群基底所得到的周期矩阵是黎曼面的共形不变量。全纯微分之积得到全纯二次微分。全纯二次微分统摄了黎曼面上所有的叶状结构(foliation)。曲面自同构的同伦类构成的曲面映射类群在所有的叶状结构空间上的作用本质地反映了曲面自同胚的特性。

Ricci曲率流可以从给定的曲率来得到相应的黎曼度量,并且和初始度量共形等价。Ricci曲率流解的存在性,唯一性,以及离散曲率流的解到连续曲率流的解的收敛性是我们研究的重点。

首先,我们将光滑曲面用多面体曲面近似。我们证明存在光滑曲面的三角剖分方法,使得离散曲面的曲率测度收敛到光滑曲面的曲率测度。证明依赖于离散法丛(normal cycle)理论的建立,这一理论的本质是说离散逼近要逼近曲面本身,更要逼近曲面的法丛。多面体曲面的高斯-博纳定理是组合技巧的平庸推论,由我们的逼近理论可以推出光滑曲面的高斯-博纳定理。经典高斯-博纳定理的证明是依赖于曲面光滑外微分和活动标架法。活动标架法导出了联络形式和曲率形式,曲率形式的超渡显示了曲面单位切丛的非平凡性。

经由余弦定理,我们得到微分余弦定理的对称性,由此构造了离散曲面上的熵能量,熵能量的梯度流就是Ricci曲率流。熵能量的凸性蕴含着Ricci曲率流解的唯一性,或曰局部刚性。

离散Ricci曲率流人为被三角剖分所桎捁,动态保持三角剖分的Delaunay特性使得曲率流更为灵活自然。我们将离散曲面上的顶点视为奇异点,考察曲面上所有的以顶点为奇异点的平直度量所构成的Teichmuller空间,和同样拓扑曲面上所有的Penner's 带装饰的双曲度量构成的Teichmuller空间。一个关键的观察是我们可以建立两个Teichmuller空间之间的微分同胚,使得Delaunay三角剖分被保持。这样,离散Ricci曲率流解的存在性得以证明,从而推出离散曲面的单值化定理:所有带度量曲面,其度量可以共形地形变为带有常值高斯曲率的度量。

离散Ricci曲率流解到光滑解的收敛性本质上等价于离散图上的调和函数梯度估计。给定一个无穷图,如果其体积增长率为二次,则其等效电阻为无穷大,其上不存在正的非常数调和函数。另一种解释是如果无穷图和欧式空间粗糙等距,则其上不存在正的非常数调和函数。依随三角剖分的不断加密,其极限图上不存在有界调和函数。离散Ricci曲率流中的离散共形因子函数恰为调和,由此必为常数,因此离散曲面的局部变换为相似变换,必然共形。梯度估计可以推测热核估计,热核决定了黎曼度量。我们证明了离散热核决定离散度量。

模空间和Teichmuller空间固定曲面拓扑,其上所有可能的共形结构形成的空间为模空间,模空间的万有覆盖空间为Teichmuller空间,其覆盖变换群为曲面映射类群。每一个曲面自同胚都诱导Teichmuller空间的一个自映射。但是有时这一映射存在不动点,所以模空间不是流形,而是orbifold。

给定两个共形结构,它们之间的每一个微分同胚都对应着一个Beltrami微分。Beltrami微分的几何解释是微分同胚将无穷小椭圆映成无穷小圆,椭圆的偏心率和长轴方向定义了Beltrami微分。所有微分同胚中,存在唯一的一个最为接近共形映射,被称为Teichmuller映射。这时,在源曲面和目标曲面上,存在全纯二次微分,Teichmuller映射保持水平轨迹的结构,并在全纯二次微分诱导的自然坐标下为仿射变换, 同时其Beltrami微分由全纯二次微分给出。由此,一个共形结构上所有全纯二次微分构成的空间就是此共形结构在Teichmuller空间上的切空间,切空间和整体空间同构,从而Teichmuller空间为单连通的。每条全纯二次微分空间中的射线定义了Teichmuller空间的一个边界点,如此我们得到Teichmuller空间的紧化。

曲面映射类作用在紧化的Teichmuller空间上,由Brower不动点定理,存在不动点。根据不动点的个数可以将映射类进行分类:周期映射,保持某个双曲度量不变;pseudo-anosov映射, 保持两个叶状结构不变;可降解映射,保持一族简单闭曲线不变,将曲面沿着不变闭曲线切开,映射在每一个连通分支上的限制都可归结为以上两类映射之一。

为了将理论从经典光滑流形推广到离散空间,我们应用并发展了如下的理论工具:图上的调和分析,微分余弦理论,离散法丛理论,曲面有限元理论,图论,单纯同调理论,凸几何理论,离散Hodge理论,离散Ricci流理论。

由此,学习了这门课程,我们可以完成如下的计算任务:

给定任意维流形的单纯剖分,我们可以计算其CW复形分解,得到其上下同调群,和一维同伦群。

给定带度量曲面,我们能够计算deRham上同调群,得到调和微分,全纯微分形式和全纯二次微分。

给定带度量曲面,经由全纯微分形式,得到曲面共形结构,进一步得到曲面共形模。

固定共形结构,我们用Ricci曲率流,可以得到同共形等价的所有可能的黎曼度量,特别是曲面单值化度量。

固定曲面拓扑,我们用全纯二次微分,可以得到所有可能的共形结构。

固定两个黎曼面,我们用Beltrami微分和辅助度量方法,可以得到它们之间所有可能的微分同胚。

展望未来,我们希望如下问题在近期内得以彻底解决:

如何计算Teichmuller映射对应的全纯二次微分;

如何精准地表示曲面上所有的叶状结构构成的空间,如何表达曲面映射类群在此空间上的作用;

如何找到曲面映射类的不变度量,不变叶状结构;

粗糙等距下,图的梯度估计如何保持;

度量曲面的共形欧式嵌入浸入。

十五年前这一切概念只存在于抽象的数学之中,现在已经成为计算机科学的中心腹地。计算共形几何已经从一粒种子生长成一棵幼树,虽然矮小单薄,但是生命力强大。她为我们提供了一个崭新的方式去体验自然的真理,用眼(可视化),用手(编程),用脑(推理),用心(审美),更用计算机。。。

【老顾谈几何】邀请国内国际著名纯粹数学家,应用数学家,理论物理学家和计算机科学家,讲授现代拓扑和几何的理论,算法和应用。

转载本文请联系原作者获取授权,同时请注明本文来自顾险峰科学网博客。

链接地址:http://blog.sciencenet.cn/blog-2472277-928927.html

上一篇:道可道,非常道-浅谈代数拓扑

下一篇:高山仰止-从微分几何角度浅谈陈省身示性类

计算共形几何是计算机科学和,科学网—计算共形几何概览 - 顾险峰的博文相关推荐

  1. 计算共形几何是计算机科学和,科学网—《计算共形几何》教程第一章 - 顾险峰的博文...

    [最近老顾等人合著的汉语教程<计算共形几何>已经完成初稿.这里我们将第一章公布,其他章节会在清华暑期课程中讲授.希望大家批评指正,不吝赐教.有兴趣预定者,请联系gu@cmsa.fas.ha ...

  2. 曲率流的计算机应用,科学网—清华笔记:计算共形几何讲义 (23)离散曲面曲率流 (Discrete Surface Ricci Flow)V - 顾险峰的博文...

    前面我们介绍了离散曲面的曲率流理论,曲面上配备着欧氏度量带有奇异点.这次,我们介绍双曲离散曲面的曲率流理论.对于欧拉示性数为负的曲面,其单值化度量自然是双曲度量.双曲度量具有非常多的优点,因此在工程实 ...

  3. 计算机控制采样周期计算公式,科学网—计算步长和采样周期 - 王永林的博文

    采用计算机进行数值模拟或系统仿真时,经常需要采用数值解法求解模型方程.系统模型一般采用微分方程描述(也有其它形式),采用数值方法求解,就是采用迭代法逐步逼近. 一.计算步长 积分与微分是互逆的,对微分 ...

  4. 计算特征数据matlab代码,科学网—MATLAB特征提取代码 - 蒋样明的博文

    PS=imresize(PS,[300,300],'bilinear');%归一化大小 PS=rgb2gray(PS); [m,n]=size(PS);                       % ...

  5. 悖论对计算机科学影响,科学网—基于对角线引理和维特根斯坦思想对于悖论的分析 - 庄朝晖的博文...

    第六届分析哲学会议发言稿 庄朝晖,厦门大学计算机科学系 悖论的定义 ...悖论的出现:在20世纪初期关于数学基础的讨论中,出现了大量的悖论,比如康托尔悖论.罗素悖论.理查德悖论等等. 对角线引理 .. ...

  6. 群论和计算机科学,科学网—魔方和群论(2) - 顾险峰的博文

    魔方方面的研究主要有两个方向:人类复原魔方算法和计算机复原魔方算法.这两类算法非常不同,但是都是基于群论的基本原理. 人类复原魔方算法是将魔方群分解成子群的乘积,然后构造每个子群的生成元,应用交换子. ...

  7. 计算机组装与维护手写笔记,科学网—计算机技术的简单小结 - 熊伟的博文

    我的博客已经搬家到 xiongbox.com  欢迎访问熊伟博士的网站! 1.计算机技术是指计算机领域中所运用的技术方法和技术手段.计算机技术具有明显的综合特性,它与电子工程.应用物理.机械工程.现代 ...

  8. python 画风场 scipy_科学网-Python: 扩展库SciPy-刘洋洋的博文

    1 什么是SciPy? SciPy是基于Python的软件生态系统,开源,主要为数学.科学和工程服务. SciPy包含6个核心库: NumPy, SciPy library, Matplotlib, ...

  9. linux去重文件第一列,科学网—Linux实用命令 - 刘洪波的博文

    1.查看系统访问日志(实时) tail --lines=100 -f /var/log/secure 2.按列合并文件 paste file1 file2 >combind.txt 3.按行合并 ...

最新文章

  1. python任务队列 http_基于Python开发的分布式任务队列:Celery
  2. 旋转质量轮实验对象上的传感器
  3. python交互式命令_在python中运行交互式命令
  4. python null byte_如何以“正确”的方式处理带有nullbytes的Python unicode字符串?
  5. java在所选路径中找到分号,classpath
  6. 聊几个与赚钱相关的小事情
  7. oracle rds 运维服务_从运维的角度分析使用阿里云数据库RDS的必要性–你不应该在阿里云上使用自建的MySQL/SQL Server/Oracle/PostgreSQL数据库...
  8. c语言实现线性表的算法,数据结构算法代码实现——线性表的定义(一)
  9. yield表达式形式的应用
  10. 【Vue.js 3.0】不使用根标签,Error:[vue/no-multiple-template-root] The template root requires exactly one ...
  11. php Immutable,Immutable.js详解
  12. 【Matlab】滤波函数
  13. 微软说:麻将AI系统终获突破,人工智能商学院来啦!
  14. Clark变换与Park(派克)变换
  15. 使用n切换node版本
  16. 微信记账本小程序(代码+功能介绍)
  17. python实现图像二分类特异度(numpy)
  18. 低功耗设计实现中secondary power pin的连接方法汇总
  19. 二进制(通俗简要版,一看就会)
  20. mysql 1062_mysql #1062 报错

热门文章

  1. The History of Tcl/Tk
  2. 一次失败的面试经历:我只想找个工作,你却用面试题羞辱我!
  3. 思维模型 波特五力模型
  4. 忻州----穿上旗袍为家乡代言!
  5. kindeditor html5 播放,kindeditor编辑器,视频不显示,视频格式embed改成video
  6. serializeArray 与 serialize
  7. 关于duilib的字体font
  8. windows 7 active directory用户和计算机,Windows – Active Directory用户和计算机上的Active Directory管理中心...
  9. 8.18 B组 T1 分火腿
  10. 一次 WIN7X64无法定位序数“3906”于动态链接库LIBEAY32.dll处理记录