离散外微积分(DEC:Discrete Exterior Calculus)基础

原文链接

“若人们不相信数学简单,只因为他们未意识到生命之复杂。”——Johnvon Neumann

DEC主要讨论离散情况下的外积分,它在计算机领域有重要用途。我们知道,使用计算机来处理几何图形的时候是不可能完全光滑的(计算机是只有0和1组成的离散化世界),利用DEC的概念也给我们提供了一种刻画离散几何的更好的工具。比如在几何分析中常用的“有限元分析(Finite Element Method)”中使用基于DEC的方法可以使用未uniform的曲面,更加方便简单。

外代数(Exterior Algebra)

首先谈谈外微积分(Exterior Calculus),它也叫外代数。

外代数是一类关于外积分(Exterior Product)的代数,外积分也叫楔形积分(Wedge Product),因为它的运算符号是楔形的。传统的线性代数涉及代数加法和标量乘法,但是不足以满足实际需求,于是定义了具有特色的外积分和内积分。

外积有独特的运算法则,例如在矩阵乘法中:

A = [a1, a2, a3, a4]^T

B = [b1, b2, b3]

C是A和B的外积,则

C = [b1a1, b2a1,b3a1]

[b1a2, b2a2, b3a2]

[b1a3, b2a3, b3a3]

[b1a4, b2a4, b3a4]

外积分还有其他的特点及性质,可以参阅相关教科书。外积分在不同的情况中有很多不同的用途,所以它就独自成为了一个体系,叫外代数。

虽然有外代数,但是没有内代数,只有内积,它们常常如下定义:

A = [a1, a2, a3,]^T

B = [b1, b2, b3]

C是A和B的内积,则

C = a1 * b1 + a2 * b2 + a3 * b3

DEC的微分形态(Differential Form)

DEC主要讨论离散情况下的外积分,它的微分形态如下:

0 形态:表示该积分在标量场(Scalar Field)中。

1 形态:表示该积分在一维场中,比如向量场(Vector Field)。

2 形态:与1形态相同。

3 形态:与0形态相同。

牛顿-莱布尼兹法则(Newton-LeibnitzFormula)

也叫First Fundamental theorem,是求解积分的一种常用方法。在1形态下有:

Int(a,b)(dF) = int(a,b)(f(x)dx) = F(b) –F(a)

单纯形(Simplex)和Chain

首先解释凸壳(Convex Hull, 也叫Convex Envelope)。给定许多点,凸壳是一个在向量空间的点集(Set),这个点集包含了所有的点而且是最小的(Minimal)。

Simplex是DEC的基本单元,有以下几种类型:

0-Simplex:点。

1-Simplex:边,或两个不同点的凸壳。

2-Simplex:三角形,或三个不同点的凸壳。

3-Simplex:四面体,或四个不同点的凸壳。

到这里我们有一个有趣的发现,一个n-Simplex必定有n+1个点组成。并且一个n-Simplex里面也包含了更小的Simplex。比如,一个三角形是2-Simplex,共有3个点,如果取其中两个点组成一个新的1-Simplex, 它也叫2-Simplex的一个1-面(1-Face)。很明显,三角形有3个0-Face,3个1-Face和1个2-Face(它自己)。

所以当我们考虑某一个n-Simplex时,要想到它其实也是许许多多更小的Simplex。

一个n-Simplex拥有很多(n-1)-Face。所有的(n-1)-Face可以组成一个(n-1)-Chain. 这所有的(n-1)-Face也可以合起来叫这个n-Simplex的(n-1)-Face。

单纯复合体(Simplicial Complex)

SimplicialComplex是很多Simplex的集合。这个集合中维度最高的Simplex的维度就是Simplicial Complex的维度。比如,一个2维SimplicialComplex中必须存在至少一个2-Simplex,也就是三角形,但是不能存在任何更高维度的Simplex例如四面体。

假设K是一个Simplicial Complex,它必定满足以下条件:

1、这个集合中的所有的Simplex的任何一个面都必须也同时在K中。换句话说,没有破损的Simplex;

2、集合中如果任意两个Simplex存在交集的话,这个交集必须是一个公共面而不是别的什么。也就是说,一个2维Simplicial Complex中如果存在两个三角形,这两个三角形要么不相接,如果相接只能是顶点到顶点的连接,不能其中一个的顶点直接安在另外一个的边上。

所以Simplicial Complex是一个几何完整并且合理拼接的图形,它有利于我们直接分析它的点、边、面的关系,还原它的拓扑结构。只要我们找到了一个分析方法,任何能化简为单纯复合体的复杂物体都可以适用。

刚才提到过Chain其实是一种有方向性的Simplicial Complex.

离散流形(Discrete Manifold)

也叫做Simplicial Manifold。一个n维的Simplicial Manifold就是一个n维的Simplicial Complex,另外还满足了以下条件:

对于其中的所有的n维Simplex,如果他们是接壤的,则必须构成一个n维“球”(Ball),或“半球”(Half Ball)如果是边界点。如此的话,每个(n-1)-Simplex就最多与两个n-Simplex接壤。比如,一个三角形是一个2维Simplicial Manifold。如果把三角形中间的面挖去,只剩下一个三角架,就不再是Manifold了,因为每个孤立的点不再能够满足球或半球的定义。

这样的网孔(Mesh), 也叫Manifold Mesh。通过以上的介绍,可以直观地理解为这种Mesh(通常是3维的)没有孤立的点、线、面,每条线也仅与一(边界时)或两个面(非边界时)接壤。同时也具有Simplicial Complex的所有特点(无破损,限制公共面)。非常纯净、整齐。

同源(Homology)

Homology是一个基本拓扑概念,当我们谈到这个概念时,设定它必须是Simplicial Complex. Homology就是用来讨论图形特性(比如相似性)的一种工具。

先看一个例子。给定三个图形:圆形、方形和圆环,我们认为圆形和方形比较相似,而圆环比较特别,因为圆形和方形可以通过变形相互转换,圆环则不可能。

如果两个Simplicial Complex有相似的结构,它们的内部就具有相同的“洞”(Hole)的数量,这两个图形就叫做“同源”(Homology)。Homology提供了研究几何不变性(Invariants)的重要方法。

边界(Boundary)

每个n-Simplex都是有方向的,它的方向定义如下:n-Simplex的n+1个(n-1)-Simplex(也可以叫(n-1)-Faces)可以用两种排列次序,其中的任意一种是它的方向。

比如,三角形是2-Simplex,它有3条边(1-Simplex),这三条边既可以顺时针组成这个三角形,也可以逆时针组成这个三角形。

一般而言,为了方便我们在考察一个n-Simplex的时候会假定它只有一个方向(通常是逆时针方向),本系列以后的所有内容都是按这样的规则。

所谓边界(Boundary),就是指某个n-Simplex对应的(n-1)-Faces中的任何一个Face。如果这条Boundary的方向与规定的方向一致,记为“+”的boundary,用“+1”表示;反之则为“-”的boundary,用“-1”表示。

边界运算符

定义运算符P为边界运算符。当此符作用于一个n-Simplex时,取得的结果是其所有正向的边界。

例如:三角形(A,B,C)有边界关系为:A—〉B—〉C—〉A。那么P(A,B,C) = (A,B) + (B,C) + (C,A)。如果是A—〉B—〉C & A—C,则P(A,B,C) = (A,B) + (B,C) - (C,A)。读者可以自行画图理解。

值得指出的是,如果n-Simplex是空集,即0-Simplex,那么P(0-Simplex) = (-1)-Simplex。

前面说到一组Simplex是一个Chain,边界运算符其实是把一个Simplex转化成一个Chain了。如果在Chain上取边界运算,等于把Chain上的每一个Simplex取边界运算然后加起来。有趣的是如果用“+1”“-1”表示边界,那么P(Chain)一定等于零!

因为边界运算符把n-Simplex映射到(n-1)-Simplex空间中,所以可以把它看成一个(N-1)*N的矩阵。这是一个稀疏矩阵,因为边界只于向接壤的部分有关。如果用+1-1表示法,这个矩阵就只含有0, +1和-1三种元素。

时间: 2024-10-12 21:02:08

离散外微积分(DEC:Discrete Exterior Calculus)基础的相关文章

Discrete.Differential.Geometry-An.Applied.Introduction(sig2008)笔记

-------------------------------------------------------------- Chapter 1: Introduction to Discrete Differential Geometry: The Geometry of Plane Curves . A better approximation than the tangent is the circle of curvature. . If the curve is sufficientl

Discrete.Differential.Geometry-An.Applied.Introduction(sig2013) 笔记

DISCRETE DIFFERENTIAL GEOMETRY : AN APPLIED INTRODUCTION Last updated: November 19, 2015 ======================================== Ch 1. INTRODUCTION 把exterior calculus作为模型处理的语言.the exterior calculus of differential forms is the modern language of dif

Surface curvature

在介绍曲面曲率之前,我们必须引入微分几何的一些基本概念. 如下图所示,$\mathbb{R}^3$空间中的一个连续曲面可以表示为由一个$\mathbb{R}^2$区域$M$通过映射函数$f:M\to\mathbb{R}^3$到$\mathbb{R}^3$中的一个集合,即$f(M)$.其中,$df(X)$为平面上向量$X$在曲面上的映射,也即为$f(M)$在点$p$处的一个切向量. 如果曲面点$p$处的某一向量$u$对点$p$处所有的切向量都满足 \[ df(X) \cdot u = 0 \] 则

java基础知识总结8(数据库篇1)

一. Oracle的安装(windowXP.win7.Linux)和卸载 1.1 Oracle的安装 1.1.1 在WindowsXP.Win7下安装 第一:解压win32_11gR2_database_1of2.win32_11gR2_database_2of2,生成detabase目录 第二:安装oracle A.点击setup图标即可,注意:安装目录不要含有中文 B.在弹出的第一个界面中取消更新选择项,点击下一步 C.在弹出的警告框中选择是 D.选择创建和配置数据库选项,下一步 E.选择桌

【深度分解】听趣拍云产品经理剖析视频基础知识(2)

"随着技术的不断进步,视频技术的制作加工门槛逐渐降低,信息资源的不断增长,同时由于视频信息内容更加丰富完整的先天优势,在近年来已经逐渐成为主流.在基础知识(1)里面已经讲了模拟时代和数字化时代的视频技术.接下来将对视频编码与压缩.画面压缩.运动压缩.互联网视频应用的到来做一个详细的介绍." 视频编码与压缩 视频编码与压缩,是数字化视频非常重要的技术,以至于它直接影响到视频在各个领域的应用.如果没有视频编码技术的不断提高,我们今天也不可能在方方面面享受到视频的便利性. 首先,视频编码是一

数据库 SQL 外键约束 多表查询

多表设计与多表查询 1.外键约束        表是用来保存现实生活中的数据的,而现实生活中数据和数据之间往往具有一定的关系,我们在使用表来存储数据时,可以明确的声明表和表之前的依赖关系,命令数据库来帮我们维护这种关系,像这种约束就叫做外键约束.          定义外键约束  foreign key(本表的列名 ordersid) references orders(id)(引用的表明及列名)     create table dept(         id int primary key

【SQL】SQL中笛卡尔积、内连接、外连接的数据演示

SQL的查询语句中,常使用到内连接.外连接,以及连接的基础--笛卡尔积运算. 在简单的SQL中,也许我们还分辨清楚数据如何连接,一旦查询复杂了,脑子也犯浆糊了,迷迷糊糊的. 本文,简单以数据形式记录连接的数据结果,在迷糊时可翻阅. 以MySQL运行. > 原始的表 select * from t_user u; select * from t_address a; select * from t_phone p; > 笛卡尔积 select * from t_user, t_address;

Oracle内连接、外连接、右外连接、全外连接小总结

数据库版本:Oracle 9i 表TESTA,TESTB,TESTC,各有A, B两列 A B 001 10A 002 20A A B 001 10B 003 30B A B 001 10C 004 40C 连接分为两种:内连接与外连接. A.内连接 内连接,即最常见的等值连接,例: SELECT * FROM TESTA,TESTBWHERE TESTA.A=TESTB.A 结果 A B A B 001 10A 001 10B B.外连接 外连接分为左外连接,右外连接和全外连接. 1.  左外

Matlab信号处理基础

简介 离散傅立叶.离散余弦和离散小波变换是图像.音频信号常用基础操作,时域信号转换到不同变换域以后,会导致不同程度的能量集中,信息隐藏利用这个原理在变换域选择适当位置系数进行修改,嵌入信息,并确保图像.音频信号经处理后感官质量无明显变化. 一维离散傅立叶变换对定义 一维离散傅里叶变换: 一维离散傅里叶逆变换: 一维离散余弦变换对定义 一维离散余弦正变换: 一维离散余弦反变换: 一维连续小波变换对定义 一维连续小波变换,其总h(t)是小波母函数: 一维连续小波逆变换: 二维离散傅立叶变换对定义 二