凸优化:仿射集、凸集、凸锥
凸集
仿射集合和凸集
主要内容:
- 直线和线段
- 仿射集
- 凸集
- 凸锥
直线和线段
定义1(直线和线段)
对于RnR^nRn空间中不同的两点∀x1≠x2\forall x_1\ne x_2∀x1=x2,令
y=θx1+(1−θ)x2=x2+θ(x1−x2),y=\theta x_1+(1-\theta)x_2=x_2+\theta(x_1-x_2),y=θx1+(1−θ)x2=x2+θ(x1−x2),
当θ∈R\theta\in Rθ∈R时,具有上述形式的点,即构成了穿过x1,x2x_1,x_2x1,x2的直线。
当θ∈[0,1]\theta\in [0,1]θ∈[0,1]时,具有上述形式的点,即构成了以x1,x2x_1,x_2x1,x2为端点的线段。
从第二个等式来看,上述形式的点以x2x_2x2为起点,x1−x2x_1-x_2x1−x2为方向,θ\thetaθ衡量了点在方向x1−x2x_1-x_2x1−x2上距离起点x2x_2x2的距离。
仿射集
定义2(仿射集)
通过集合C⊆RnC\subseteq R^nC⊆Rn中任意两点的直线仍旧在集合CCC中,称集合CCC是仿射的。
等价于,∀x1,x2∈C\forall x_1,x_2\in C∀x1,x2∈C,θx1+(1−θ)x2∈C,θ∈R\theta x_1+(1-\theta)x_2\in C,\theta\in Rθx1+(1−θ)x2∈C,θ∈R
注
仿射的拉丁语意为相关,即∑Xi∈C\sum X_i\in C∑Xi∈C
例1
一条直线也是仿射集,因此直线上任意两点构成的直线仍在原先的直线上。
定义3(仿射组合)
扩展至多点,当θ1+...+θn=1\theta_1+...+\theta_n=1θ1+...+θn=1,称具有形式θ1x1+...+θnxn\theta_1x_1+...+\theta_nx_nθ1x1+...+θnxn的点为x1,...,xnx_1,...,x_nx1,...,xn的仿射组合。
例2
一类特殊的仿射集:

对于线1,当x1,x2∈Cx_1,x_2\in Cx1,x2∈C时,由于仿射集要求组合点的系数之和为1,显然x1+x2=x∗∉Cx_1+x_2=x^*\notin Cx1+x2=x∗∈/C,从图像上即组合后的点x∗x^*x∗不落在线1上;但对于线2,当x1,x2∈Cx_1,x_2\in Cx1,x2∈C时,x1+x2=0∈Cx_1+x_2=0\in Cx1+x2=0∈C,即此时C是仿射集,但此时x1x_1x1和x2x_2x2的系数之和并非1。
定义4(仿射包)
包含C的最小仿射集(集合C中的点的所有仿射组合构成的集合):
affC={θ1x1+...+θnxn∣xi∈C,i=1,...,n,θ1+...+θn=1}\text{aff}C=\{\theta_1x_1+...+\theta_nx_n|x_i\in C,i=1,...,n,\theta_1+...+\theta_n=1\}affC={θ1x1+...+θnxn∣xi∈C,i=1,...,n,θ1+...+θn=1}
例3
线性方程组的解集(超平面上的点):
线性方程组的解集:C={x∣Ax=b}C=\{x|Ax=b\}C={x∣Ax=b},A∈Rm×nA\in R^{m\times n}A∈Rm×n,b∈Rmb\in R^mb∈Rm,直接可从仿射集的定义可判断。
假设x1,x2∈Cx_1,x_2\in Cx1,x2∈C,∀θ\forall \theta∀θ,有:
A(θx1+(1−θ)x2)=θb+(1−θ)b=bA(\theta x_1+(1-\theta)x_2)=\theta b+(1-\theta )b=bA(θx1+(1−θ)x2)=θb+(1−θ)b=b
因此C中任意点的仿射组合仍在C中,C是仿射集。
换个角度,对比例1,一条线是仿射集,很容易得出超平面也是仿射集。
凸集
定义5(凸集)
通过集合C⊆RnC\subseteq R^nC⊆Rn中任意两点的线段仍旧在集合CCC中,称集合CCC是凸集。
等价于,∀x1,x2∈C\forall x_1,x_2\in C∀x1,x2∈C,θx1+(1−θ)x2∈C,θ∈[0,1]\theta x_1+(1-\theta)x_2\in C,\theta\in [0,1]θx1+(1−θ)x2∈C,θ∈[0,1]
注
将仿射集合中的θ∈R\theta\in Rθ∈R限制在θ∈[0,1]\theta\in [0,1]θ∈[0,1]即成了凸集。换个角度该定义意为,从点x2x_2x2看向点x1x_1x1,中间无任何阻碍,都能看见。
例4
一个集合CCC是仿射集,必然也是凸集。因为对仿射集CCC而言,任意两点连线的直线属于CCC,那么任意两点连线的线段必然属于CCC,满足凸集的定义。
例5
下图是R2R^2R2空间中一些常见的凸集和非凸集的例子:

左一内部点全能连通,因此是凸集;中间的图,两点间连线有一部分点不在集合内,不是凸集;右一,边界有一部分缺失,当连接缺失两边点,其构成的线段不在集合内,因此不是凸集。
定义6(凸组合)
当θ1+...+θn=1,θi≥0\theta_1+...+\theta_n=1,\theta_i\ge 0θ1+...+θn=1,θi≥0,称具有形式θ1x1+...+θnxn\theta_1x_1+...+\theta_nx_nθ1x1+...+θnxn的点为x1,...,xnx_1,...,x_nx1,...,xn的凸组合。
定义7(凸包)
包含C的最小凸集:
convC={θ1x1+...+θnxn∣xi∈C,θi≥0,i=1,...,n,θ1+...+θn=1}convC=\{\theta_1x_1+...+\theta_nx_n|x_i\in C,\theta_i\ge0,i=1,...,n,\theta_1+...+\theta_n=1\}convC={θ1x1+...+θnxn∣xi∈C,θi≥0,i=1,...,n,θ1+...+θn=1}
例6
按定义,很容易验证,空集ϕ\phiϕ,单点集{x}\{x\}{x},空间RnR^nRn是仿射集,自然也都是凸集。
例7
半空间{x∣aTx≤b}\{x|a^Tx\le b\}{x∣aTx≤b}是凸集,但非仿射集,如下R2R^2R2上的半空间:

半空间(阴影区域)任意两点线段的连线仍在半空间内。
例8
半空间的另一种表述形式:假设RnR^nRn空间中有两个不同的点a,ba,ba,b,则所有距离aaa比bbb近的点构成的集合(2范数):
{x∣ ∣∣x−a∣∣2≤∣∣x−b∣∣2}={x∣xTx−2aTx+aTa≤xTx−2bTx+bTb}={x∣2(b−a)Tx≤bTb−aTa}\{x|\ ||x-a||_2\le ||x-b||_2\}=\{x|x^Tx-2a^Tx+a^Ta\le x^Tx-2b^Tx+b^Tb\}=\{x|2(b-a)^Tx\le b^Tb-a^Ta\}{x∣ ∣∣x−a∣∣2≤∣∣x−b∣∣2}={x∣xTx−2aTx+aTa≤xTx−2bTx+bTb}={x∣2(b−a)Tx≤bTb−aTa}
例9
多面体,有限个线性不等式和等式的解集:
C={x∣ajTx≤bj,j=1,...,m,cjTx=dj,j=1,...,p}C=\{x|a_j^Tx\le b_j,j=1,...,m,c_j^Tx= d_j,j=1,...,p\}C={x∣ajTx≤bj,j=1,...,m,cjTx=dj,j=1,...,p}
可见多面体是有限个半空间和超平面的交集,因此半空间也是多面体,自然多面体也是凸集。
例10
Euclid球和椭球,RnR^nRn空间中的Euclid球有如下形式:
B(x0,r)={x∣ ∣∣x−x0∣∣2≤r}={x0+ru∣ ∣∣u∣∣2≤1}B(x_0,r)=\{x|\ ||x-x_0||_2\le r\}=\{x_0+ru|\ ||u||_2\le1\}B(x0,r)={x∣ ∣∣x−x0∣∣2≤r}={x0+ru∣ ∣∣u∣∣2≤1}
即以x0x_0x0为圆心,rrr为半径的球。根据定义,假设∣∣x1−x0∣∣2≤r,∣∣x2−x0∣∣2≤r||x_1-x_0||_2\le r,||x_2-x_0||_2\le r∣∣x1−x0∣∣2≤r,∣∣x2−x0∣∣2≤r。
∣∣θx1+(1−θ)x2−x0∣∣2=∣∣θ(x1−x0)−(1−θ)(x2−x0)∣∣2≤∣∣θ(x1−x0)∣∣2+∣∣θ(x2−x0)∣∣2≤r||\theta x_1+(1-\theta)x_2-x_0||_2=||\theta(x_1-x_0)-(1-\theta)(x_2-x_0)||_2\le ||\theta(x_1-x_0)||_2+||\theta(x_2-x_0)||_2\le r∣∣θx1+(1−θ)x2−x0∣∣2=∣∣θ(x1−x0)−(1−θ)(x2−x0)∣∣2≤∣∣θ(x1−x0)∣∣2+∣∣θ(x2−x0)∣∣2≤r
凸组合的点仍旧在B(x0,r)B(x_0,r)B(x0,r)内,因此为凸集。(R2R^2R2中为圆,R3R^3R3中为球,其内点的连线仍在其中,自然满足凸集的要求。)
椭球也是凸集,ε={x∣(x−x0)TP−1(x−x0)≤1}={x0+Au∣ ∣∣u∣∣2≤1}\varepsilon=\{x|(x-x_0)^TP^{-1}(x-x_0)\le 1\}=\{x_0+Au|\ ||u||_2\le1\}ε={x∣(x−x0)TP−1(x−x0)≤1}={x0+Au∣ ∣∣u∣∣2≤1},其中PPP是对称正定矩阵。
锥
定义8(凸锥)
∀x1,x2∈C\forall x_1,x_2\in C∀x1,x2∈C,θ1x1+θ2x2∈C,θ1,θ2≥0\theta_1 x_1+\theta_2x_2\in C,\theta_1,\theta_2\ge0θ1x1+θ2x2∈C,θ1,θ2≥0
注
既然θ1,θ2\theta_1,\theta_2θ1,θ2都可取到0,那么0点必然属于集合CCC
下图即是凸锥,根据向量加法的平行四边形法则,θ1x1+θ2x2\theta_1x_1+\theta_2x_2θ1x1+θ2x2即红箭头交点处,改变θ1,θ2\theta_1,\theta_2θ1,θ2的取值,红箭头交点可取遍锥内所有的点。

定义9(锥组合)
当θi≥0,i=1,...,n\theta_i\ge0,i=1,...,nθi≥0,i=1,...,n,称具有形式θ1x1+...+θnxn\theta_1x_1+...+\theta_nx_nθ1x1+...+θnxn的点为x1,...,xnx_1,...,x_nx1,...,xn的锥组合。
定义10(锥包)
包含C的最小凸锥:
{θ1x1+...+θnxn∣xi∈C,θi≥0,i=1,...,n}\{\theta_1x_1+...+\theta_nx_n|x_i\in C,\theta_i\ge0,i=1,...,n\}{θ1x1+...+θnxn∣xi∈C,θi≥0,i=1,...,n}
例11
过原点的直线是凸锥。
例12
任意子空间是仿射的、凸锥。
例13
半正定锥,
用SnS^nSn表示n维对称矩阵,Sn={X∈Rn×n∣X=XT}S^n=\{X\in R^{n\times n}|X=X^T\}Sn={X∈Rn×n∣X=XT},
用S+nS^n_+S+n表示n维半正定对称矩阵,S+n={X∈Sn∣X⪰0}S^n_+=\{X\in S^n|X\succeq0\}S+n={X∈Sn∣X⪰0},其是凸锥,首先过原点,再者,若θ1,θ2≥0\theta_1,\theta_2\ge 0θ1,θ2≥0,A,BA,BA,B是半正定对称矩阵,
xT(θ1A+θ2B)x=θ1xTAx+θ2xTBx≥0x^T(\theta_1A+\theta_2B)x=\theta_1x^TAx+\theta_2x^TBx\ge 0xT(θ1A+θ2B)x=θ1xTAx+θ2xTBx≥0
满足凸锥的定义。
用S++nS^n_{++}S++n表示n维正定对称矩阵,S++n={X∈Sn∣X≻0}S^n_{++}=\{X\in S^n|X\succ0\}S++n={X∈Sn∣X≻0},当n=1,S++n=R++S^n_{++}=R_{++}S++n=R++,此时点集不过原点,不是凸锥,但其是实数轴上的半空间,因此是凸集。
更多推荐
所有评论(0)