跳转至

仿射与凸

在本节中,我们简要介绍仿射和凸的概念。

仿射

直线的定义想必初中都已经学过,这里我们简单复习一下:

直线(line)

空间Rn中的两个点x1,x2,x1x2可以唯一确定一条直线。这条直线可以被参数化的描述为

x=θx1+(1θ)x2 s.t. θR

我们很容易发现,x=x1时,θ=1x=x2时,θ=0。由此,我们知道了线段被如何描述:

线段(line segment)

穿过两点x1,x2Rn,x1x2的线段可以被参数化的描述为

x=θx1+(1θ)x2 s.t. 0θ1

通过变化,直线还可以被表达为x=x2+θ(x1x2) s.t. θR。这给了我们另一种解读–x基点 x2方向 x1x2按参数θ缩放之和。 因此,θ给出了从x2x1的分数。当θ从0增大到1时,点xx2移动到x1;当θ在此区间以外时,点x位于x1x2以外的线上。

在复习了直线之后,我们可以给仿射集做一个定义:

仿射集

对于集合CR,若经过该集合上任意两点x1,x2C,x1x2的直线上的点都在该集合中,那么我们称集合C为一个仿射集。

线性方程{xAx=b}的解即是一个仿射集,反过来每一个仿射集也都可以表示为一个线性方程系统的解集。因此,我们有时候将仿射集成为线性(linear)。

有同学可以发现,仿射集与线性子空间非常相似。事实上,我们可以将一个仿射集看作是一个线性子空间根据一个固定向量平移而得到的。对于仿射集C与其上一点x0,集合V=Cx0={xx0xC}是一个子空间。

证明V=Cx0对于向量加法与标量乘法闭合

对于v1,v2V,α,βR,我们有v1+x0,v2+x0C,所以αv1+βv2+x0=α(v1+x0)+β(v2+x0)+(1αβ)x0C,由于C是仿射集、α+β+(1αβ)=1,αv1+βv2+x0C,易得αv1+βv2V

因此,我们可以将仿射集C表述为C=V+x0={v+x0vV}x0C上任意一点。我们将子空间V称作仿射集C相关的子空间,将子空间V的维度称作仿射集C的维度。

简单了解了直线与仿射集之后,我们可以延伸出仿射组合:

仿射组合(affine combination)

对于线性组合i=1nθixi,若其满足

i=1nθi=1

则我们称其为仿射组合。

仿射集C包含了该集合中任意点的仿射组合,也就是说,对于x1,x2,xnC,i=1nθi=1i=1nθixiC

我们可以发现,仿射组合其实可以被看成一个所有系数之和为1的特殊的线性组合。这样的特性使得我们可以抛弃原点,而将任意的点视为原点,我们都可以使用同样的线性组合去描述同一个点。

🌰

对于原点以及空间上任意一点p,我们可以将向量a, b分别表述为a+bp+(ap)+(bp)。由此,我们通过仿射组合去描述另一点时,两个描述分别为θa+(1θ)bp+θ(ap)+(1θ)(bp)=p+θaθp+(1θ)b(1θ)p=θa+(1θ)b;我们很容易发现两个描述与其起始点没有关系。

集合S的仿射包aff(S)指的是包含集合S的最小仿射集,它被定义为:

仿射包(affine hull)

aff(S)={i=1kθixik>0,xiS,θiR,i=1nθi=1}

也就是说,如果仿射集CSC,那么aff(S)C

我们注意到到现在为止我们一直在讨论仿射,而没有进入真正的主题–凸。事实上,我们已经差不多学习了凸的所有概念了–仅仅将仿射的直线换成线段,我们就能得到凸。

凸集

对于集合CR,若经过该集合上任意两点x1,x2C,x1x2线段上的点都在该集合中,那么我们称集合C为一个凸集。

凸组合(convex combination)

对于线性组合i=1nθixi,若其满足

i=1nθi=1并且θi,θi0

则我们称其为凸组合。

凸包(convex hull)

conv(S)={i=1kθixik>0, xi0,θiR,i=1nθi=1}

| 我们很容易想到,常见的三角形、方形、圆形都是凸的;而星形则是非凸的,因为任意两点之间存在有些点不在集合之内。

非凸集五角星的凸包是

凸集五边形。

如何判断一个函数是否是凸的

这个问题看似十分简单。我们只需要根据凸集的定义,选取集合中的每一个点对然后判断这两个点之间的线段上的点是否是凸的即可。对于一个有着n个点的集合,上述算法的时间复杂度仅仅为O(n?n),这甚至是一个多项式时间的算法。但是这明显还有很大的改进空间。

在本节中,我们介绍保留凸性的运算,这些运算有助于判定一个集合是否凸以及构建一个凸集。

保凸运算

交集

对于集合S1,S2,如果他们是凸集,那么S1S2也是凸集。

仿射函数

如果函数f:RnRm是一个放射函数,那么:

  • f下的凸集的像也是凸的:
SRnis convexf(S)={f(x)xS}is convex
  • f下的凸集的逆像也是凸的:
CRmis convexf1(C)={xf(x)C}is convex

透视和线性分数函数

透视函数 P:Rn+1Rn

P(x,t)=x/t,dom(P)={(x,t)t>0}

线性分数函数 f:RnRm

f(x)=Ax+bcT+d,dom(f)={xcTx+d>0}

透视函数与线性分数函数下的凸集的像和逆像也是凸的

这些属性扩展到无限个集合的交集中,也就是说,如果Sα是凸集对于每一个αA都成立,那么αASα也是一个凸集。

affine_convex.pdf