zoukankan      html  css  js  c++  java
  • 离散外微积分(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三种元素。

  • 相关阅读:
    GUI 监听事件 (两个按钮,实现同一个监听)
    GUI 监听事件
    GUI 练习
    GUI 之表格布局
    GUI 之边界布局
    GUI 之流布局
    [转帖]Linux 下解压 rar 文件
    Linux 启动、停止、重启jar包脚本
    关于linux下,ls vi等命令失效的解决方法(配置下环境变量出现问题)
    超好用的UnixLinux 命令技巧 大神为你详细解读
  • 原文地址:https://www.cnblogs.com/liangliangdetianxia/p/4021116.html
Copyright © 2011-2022 走看看