2.8, 2.13, 2.22, A1.5, A1.9.
2.8
![在这里插入图片描述](https://img-blog.csdnimg.cn/20200427184858825.png)
问以下哪些集合是多面体,如果可能的话,用 S={x∣Ax⪯b,Fx=g} 的形式来表示 S
a
S 是一个多面体,它是由 a1+a2,a1−a2,−a1+a2,−a1−a2 这四个顶点构成的,比如在 R2 中:
![在这里插入图片描述](https://img-blog.csdnimg.cn/20200428222948478.png)
集合 S 可以描述为以下几个点集的交集:
- S1 :a1,a2 所定义的平面上的点
- S2:{z+y1a1+y2a2∣a1Tz=a2Tz=0,−1≤y1≤1},其中 z 是垂直于 a1,a2 平面的向量,对于每一个固定的 y1 ,S2 为一个平行于 a2 且垂直于 a1,a2 所形成的平面的平面,变化 y1,S2 实际上是一组这样的平面。
- S3:{z+y1a1+y2a2∣a1z=a2Tz=0,−1≤y2≤1} ,这个和 S2 同理,S3 平行于 a1 ,垂直于 a1,a2 所形成的平面,是一组平面。
以上三个集合取交集刚好是图中所示的平行四边形(阴影部分)
vkTx=0, k=1,...,n−2
其中 vk 是垂直于 a1,a2 的不相关的向量(矩阵 [a1,a2]T 的零空间(null space)),对于 R2 来说这个面就不存在了,对于 R3 及以上,这个就可以想象成 a1,a2 形成的超平面。
−∣c1Ta1∣≤c1Tx≤∣c1Ta1∣
其中 c1=a1−∥a2∥22a1Ta2a2,这个可以理解为平行四边形的中心到右边的那个边做了一个垂线,中心到垂足形成的向量即为 c1(不理解的话就画一画),x 就是空间内任意一个向量,要满足的约束其实就是该向量和 c1 做内积的绝对值要小于等于 ∥c1∥22 (∣c1Ta1∣=∥c1∥22,额不理解的话就画图),这整个的意思就是限制向量 x 在平行四边形的左右两边之内。
S3 可以描述为:
−∣c2Ta2∣≤c2Tx≤∣c2Ta2∣
其中 c2=a2−∥a1∥22a2Ta1a1,这个和 S2 的道理是一样的,意思是限制 x 在平行四边形的上下两边内。
总和起来就可以得到 2n 个线性不等式的约束:
![在这里插入图片描述](https://img-blog.csdnimg.cn/20200428224703976.png)
b
是一个多面体,x 的每一个分量大于等于零,然后整体由三个等式约束。
c
不是一个多面体,它是单位球({x∣∥x∥2≤1})和非负象限 R+n 的相交部分。根据柯西-施瓦茨不等式:
xTy≤1 for all y with ∥y∥2=1⟺∥x∥2≤1
多面体的定义是多个半空间的相交部分,这里是无穷多个半平面的相交。
d
是一个多面体,S 是集合 {x∣∣xk∣≤1,k=1,...,n} 与非负象限 R+n 的交集,即证明:
xTy≤1 for all y with i=1∑n∣yi∣=1⟺∣xi∣≤1,i=1,...,n
首先假设 ∣xi∣≤1 ,则有:
xTy=i∑xiyi≤i∑∣xi∣∣yi∣≤i∑∣yi∣=1
然后再假设 x 是一个满足 xTy≤1,∑i∣yi∣=1 的一个非零向量,令 ∣xk∣=maxi∣xi∣ ,并令 ykxk=∣xk∣,yi=0(i=k),则有:
xTy=i∑xiyi=ykxk=∣xk∣=imax∣xi∣
即 maxi∣xi∣≤1,因此可以用有限的线性不等式来表示 S,即为非负象限与 {x∣−1⪯x⪯1} 之间的交集,即 2n 个线性不等式的解:
![在这里插入图片描述](https://img-blog.csdnimg.cn/20200429184623337.png)
2.13
![在这里插入图片描述](https://img-blog.csdnimg.cn/20200427184959169.png)
外积的圆锥壳:考虑秩为 k 的外积的集合,定义为 {XXT∣X∈Rn×k,rank X=k}。描述其圆锥壳。
首先可以知道 XXT⪰0,并且 rank(XXT)=k,也就是说题目中这个集合就是秩为 k 的半正定矩阵的集合,先说答案:这个集合的圆锥壳是秩 ≥k 的半正定矩阵集合再并上零矩阵。
一个集合的圆锥壳记为该集合元素的正组合,假设组合的秩 <k,即令 A,B 为秩为 k 的半正定矩阵,rank(A+B)=r<k,令 V∈Rn×(n−r) 为 A+B 零空间中的向量组成的矩阵,即 R(V)=N(A+B),即:
VT(A+B)V=VTAV+VTBV=0
又因为 A,B⪰0,因此 VTAV=VTBV=0,所以 rank A≤r,rank B≤r ,与之前 A,B 的秩为 k 的假设矛盾!
因此该集合的的圆锥壳是秩 ≥k 的半正定矩阵集合再并上零矩阵。
2.22
![在这里插入图片描述](https://img-blog.csdnimg.cn/20200427185015198.png)
分离超平面定理(sepqrating hyperplane theorem)说的是:假设 C,D 是非空的不相交的凸集,即 C∩S=∅,则存在 a=0,b 使得 aTx≤b,x∈C 以及 aTx≥b,x∈D 成立。也就是说仿射函数 aTx−b 在 C 上非正,在 D 上非负。超平面 {x∣aTx=b} 被称作 C,D 的分离超平面,如下图所示:
![在这里插入图片描述](https://img-blog.csdnimg.cn/20200501145503222.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L20wXzM3Nzk1MjQ0,size_16,color_FFFFFF,t_70)
在原始的证明中证明的是假如两个集合中存在某两个点的距离刚好等于两个集合的距离,那么这两个集合之间存在分离超平面,关键在于两个集合间距离的定义为:
dist(C,D)=inf{∥u−v∥2∣u∈C,v∈D}
其中 inf 表示下确界,意思是小于等于集合中所有元素的最大值,对于明确存在最小值的集合来说,这个下确界就等于集合中的最小值,而对于开集来说,比如 {x∣1<x<3} ,该集合的下确界就为 1 ,但是 1 并不在这个集合中。所以书本中的证明只在 C,D 是闭集,并且其中一个集合是有界的情况下才成立,这道题目要求证明的就是更加一般的情况,也就是考虑 C,D 集合为开,在这两个集合中不能找出两个元素之间的距离刚好等于两个集合之间的距离。
要分成两种情况进行讨论:
-
假设 0∈cl S ,即 0 不在集合 S 的闭包上,S 的闭包指的是 S 的内部点连通其边界上的点形成的集合,假如 S 本身就是闭集,那么闭包就等于它自己,假如是开的,那么闭包就等于它自己并上它的边界。此时对集合 {0} 和 cl S 应用已经证明过的部分,即在集合 {0} 和集合 cl S 之间是存在分离超平面的(因为这两个集合都是闭集且其中有一个集合有界!),因此存在 a=0 使得
aT(x−y)>0
对于所有 x−y∈cl S 成立,因此对于所有 x−y∈S,aTx>aTy 对于所有 x∈C,y∈D 成立。
-
假设 0∈cl S ,由于 0∈S (两集合不相交),因此 0 肯定在 S 的边界上,因此这样一来我们就不能用第一种情况里的方法了,因为 cl S 和 {0} 这两个集合是重叠的,无法应用书中已经证明的部分。再分情况进行讨论:假如 S 的内部是空的,也就是说,他在一个超平面 {z∣aTz=b} 中,又因为这个超平面通过原点,因此 b=0,因此对于所有 x∈C,y∈D,我们都有 aTx=aTy,这就是一个很 trivial 的分离超平面(这个例子可以脑补一下,就是两个集合 C,D 是紧紧贴在这个超平面两侧的,他们中的元素无限接近这个超平面,但是永远到不了);然后再假如 S 的内部不空(也就是说这两个集合只有一小部分是相互无限贴近的,其他部分中间可以离得比较远。。脑补脑补),然后考虑集合 S−ϵ={z∣B(z,ϵ)⊆S} ,其中 B(z,ϵ) 是中心为 z ,半径为 ϵ>0 的欧几里得球,直观来讲,这个集合就是把原集合 S 向里面收缩了 ϵ ,这样得到的效果就是 cl S−ϵ 是闭集且凸,并且该集合不包含 0,因此我们就可以再愉快地用课本上已经证明过的部分啦,该集合可以用一个法线向量为 a(ϵ) 的超平面与 {0} 分隔开:a(ϵ)Tz>0 for all z∈S−ϵ ,不失一般性,我们可以假设 ∥a(ϵ)∥2=1。现在令 ϵk,k=1,2,... 为一个正序列,并有 limk→∞ϵk=0,由于 ∥a(ϵk)∥2=1 ,因此序列 a(ϵk) 是收敛的,定义这个序列最终可以收敛到 aˉ ,于是对于所有的 k 有 a(ϵk)Tz>0 for all z∈S−ϵk ,因此对于所有 z∈int S 有 aˉTz>0,也就是对于所有的 z∈S 有 aˉTz≥0 ,即对于所有的 x∈C,y∈D 有 aˉTx≥aˉTy 。
A1.5
![在这里插入图片描述](https://img-blog.csdnimg.cn/20200427185037712.png)
![在这里插入图片描述](https://img-blog.csdnimg.cn/2020042718505546.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L20wXzM3Nzk1MjQ0,size_16,color_FFFFFF,t_70)
本题要证明两个闭凸锥交集的对偶等于各自对偶的和。即令 C,D 为 Rn 上的闭凸锥,证明:
(C∩D)∗=C∗+D∗
其中 + 表示集合加法:C∗+D∗ 为 {u+v∣u∈C∗,v∈D∗} ,证明步骤如下:
- 证明 C∩D 和 C∗+D∗ 为凸锥(实际上,这两个是闭的,但是不要求证明)
- 证明 (C∩D)∗⊇C∗+D∗
- 然后证明 (C∩D)∗⊆C∗+D∗ ,证明思路是:首先证明
(C∩D)∗⊆C∗+D∗⟺C∩D⊇(C∗+D∗)∗
其中可以使用以下结论:如果 K 是一个闭凸锥,那么 K∗∗=K ,之后再证明:C∩D⊇(C∗+D∗)∗ ,最后得到结论 (C∩D)∗=C∗+D∗
- 证明多面体锥 V={x∣Ax⪰0} 得对偶可以表示为:
V∗={ATv∣v⪰0}
解答:
- 证明 C∩D 和 C∗+D∗ 为凸锥:假设 x∈C∩D ,这意味着 x∈C 且 x∈D ,还意味着对于任意 θ≥0 来说, θx∈C 以及 θx∈D ,因此对于任意的 θ≥0,有 θx∈C∩D ,故 C∩D 是一个锥,并且它是凸的,因为两个凸集的交集仍是凸集。为了证明 C∗+D∗ 是一个闭凸锥,注意到 C∗ 和 D∗ 都是凸锥,因此 C∗+D∗ 是 C∗∩D∗ 的锥包,因此是一个凸锥。
- 证明 (C∩D)∗⊇C∗+D∗:假设 x∈C∗+D∗ ,我们可以把 x 写成 x=u+v ,其中 u∈C∗,v∈D∗ ,根据对偶锥的定义,对于所有 y∈C 有 uTy≥0,对于所有 y∈D 有 vTy≥0 ,可以推出对于所有 y∈C∩D 有 xTy=uTy+vTy≥0 ,这意味着 x∈(C∩D)∗ ,因此 (C∩D)∗⊇C∗+D∗
- 证明 (C∩D)∗⊆C∗+D∗:已经证过 C∩D 和 C∗+D∗ 是闭凸锥,这意味着 (C∩D)∗∗=C∩D ,因此有:
(C∩D)∗⊆C∗+D∗⟺C∩D⊇(C∗+D∗)∗
假设 x∈(C∗+D∗)∗,则对于所有 y=u+v 有 xTy≥0,其中 u∈C∗,v∈D∗ ,可以写为对于所有的 u∈C∗,v∈D∗ 有 xTu+xTv≥0 。因为 0∈C∗ 并且 0∈D∗ ,令 v=0,则对于所有的 u∈C∗ 有 xTu≥0,令 u=0,则对于所有的 v∈D∗ 有 xTv≥0。这样可以推出 x∈C∗∗=C 并且 x∈D∗∗=D ,即 x∈C∩D 。因此原式证毕。
- 证明多面体锥 V={x∣Ax⪰0} 的对偶可以表示为:V∗={ATv∣v⪰0}:使用本题证明的式子,可以得到:
V∗={x∣a1Tx≥0}∗+⋯+{x∣amTx≥0}∗
其中 {x∣aiTx≥0} 的对偶是集合 {θai∣θ≥0} (注:其中 θ 为常数,半平面的对偶是法线所在的射线,脑补一下~),因此我们得到:
V∗={θa1∣θ≥0}+⋯+{θam∣θ≥0}={θ1a1+⋯+θmam∣θi≥0,i=1,...,m}={ATθ∣θ⪰0}
即:V∗={ATv∣v⪰0} (妙啊!)
A1.9
![在这里插入图片描述](https://img-blog.csdnimg.cn/20200427185121656.png)
确定以下 Sn 的子集是否为凸集:
(a) 相关矩阵的集合:Cn={C∈S+n∣Cii=1,i=1,...,n}
(b) 非负相关矩阵的集合:{C∈Cn∣Cij≥0,i,j=1,...,n}
(c) 高相关性的相关矩阵的集合: {C∈Cn∣Cij≥0.8,i,j=1,...,n}
(a) 令 X∈Cn,Y∈Cn, 因此有 X∈S+n (也就是说 X 是一个半正定矩阵),并且 Xii=1,i=1,...,n,而且 Y∈S+n,Yii=1,i=1,...,n ,考虑 Z=θX+(1−θ)Y,θ∈[0,1],显然 Z∈S+n (因为 S+n 是凸的),并且也显然有 Zii=1,i=1,...,n ,因此有 Z∈Cn ,故 Cn 是凸集。
(b) 同理,令 X,Y 都在这个集合内,考虑 Z=θX+(1−θ)Y,θ∈[0,1],由于已经证过 Cn 为凸集,因此 Z∈Cn ,只需再验证 Zij≥0,i,j=1,...,n 是否成立,显然成立(因为 Xij≥0,Yij≥0,i=1,...,n) ,因此该集合也是凸集。
(c) 这个集合也是凸集,证明方法类似于 (b),不再赘述(手动狗头)