凸优化中:单纯形是一种多面体的证明
思路:先构造一个单纯形,然后将这个单纯形用多面体的定义形式表达出来即可。
证明:设v0,v1,⋯ ,vk∈Rn{v_0},{v_1},\cdots,{v_k} \in R^nv0,v1,⋯,vk∈Rn,v1−v0,v2−v0,⋯ ,vk−v0,{v_1}-{v_0},{v_2}-{v_0},\cdots,{v_k}-{v_0},v1−v0,v2−v0,⋯,vk−v0,线性独立,则由这k+1个点构成的单纯形为:
C=conv{v0,v1,⋯ ,vk}={θ0v0+θ1v1,⋯ ,+θkvk∣θ≽0,1Tθ=1}C=conv\left\{ {{v_0},{v_1},\cdots,{v_k}} \right\}=\left\{ {{\theta _0}{v_0}+{\theta _1}{v_1},\cdots,+{\theta _k}{v_k}\left| {\theta \succcurlyeq 0,{1^T}\theta = 1} \right.}\right\}C=conv{v0,v1,⋯,vk}={θ0v0+θ1v1,⋯,+θkvk∣∣θ≽0,1Tθ=1}
令y=[θ1,⋯ ,θk]Ty=\left[ {\theta _1}, \cdots ,{\theta _k} \right]^Ty=[θ1,⋯,θk]T,则y⩾0y \geqslant 0y⩾0且1Ty⩽1{1^T}y \leqslant 11Ty⩽1,令[v1−v0,⋯ ,vk−v0]=B∈Rn×k\left[ {v_1} - {v_0}, \cdots ,{v_k} - {v_0} \right] = B \in R^{n \times k}[v1−v0,⋯,vk−v0]=B∈Rn×k,则单纯形C中的任意一点XXX可以表示为:X=θ0v0+⋯+θkvk=v0+θ1(v1−v0)+⋯+θk(vk−v0)=v0+ByX = {\theta _0}{v_0} + \cdots + {\theta _k}{v_k} = {v_0} + {\theta _1}\left( {{v_1} - {v_0}} \right) + \cdots + {\theta _k}\left( {{v_k} - {v_0}} \right) = {v_0} + ByX=θ0v0+⋯+θkvk=v0+θ1(v1−v0)+⋯+θk(vk−v0)=v0+By
由单纯形的定义我们可知BBB中的kkk个向量线性无关,所以存在rank(B)=k⩽1rank\left(B\right)=k\leqslant1rank(B)=k⩽1,即BBB为列满秩矩阵。
接下来我们将说明任何一个列满秩矩阵BBB都可以通过线性变换变成[Ik0]\left[ \begin{matrix} {I_k} \\0 \end{matrix} \right][Ik0]这种形式。其中,Ik{I_k}Ik表示k×kk \times kk×k的单位矩阵,000表示(n−k)×k\left(n-k\right)\times k(n−k)×k的零矩阵。具体变换方法是我们可以找到一个非奇异矩阵A=[A1A1]∈Rn×nA=\left[ \begin{matrix} {A_1} \\{A_1} \end{matrix} \right]\in R^{n\times n}A=[A1A1]∈Rn×n使得AB=[A1A1]B=[Ik0]AB=\left[ \begin{matrix} {A_1} \\{A_1} \end{matrix} \right]B=\left[ \begin{matrix} {I_k} \\0 \end{matrix} \right]AB=[A1A1]B=[Ik0]。
根据上面的变换,我们可以得到:
X=v0+By⇔AX=Av0+ABy⇔[A1A1]X=[A1A1]v0+[A1BA1B]y⇔{A1X=A1v0+yA2X=A2v0X = {v_0} + By \Leftrightarrow AX=A{v_0} + ABy\Leftrightarrow\left[ \begin{matrix} {A_1} \\{A_1} \end{matrix} \right]X=\left[ \begin{matrix} {A_1} \\{A_1} \end{matrix} \right]{v_0}+\left[ \begin{matrix} {A_1}B \\{A_1}B \end{matrix} \right]y\Leftrightarrow\left\{\begin{matrix} {A_1}X = {A_1}{v_0} + y\\ {A_2}X = {A_2}{v_0}\end{matrix}\right.X=v0+By⇔AX=Av0+ABy⇔[A1A1]X=[A1A1]v0+[A1BA1B]y⇔{A1X=A1v0+yA2X=A2v0,因为y⩾0y \geqslant 0y⩾0且1Ty⩽1{1^T}y \leqslant 11Ty⩽1,所以上述两个等式可以转换成如下的不等式:
{A1X=A1v0+yA2X=A2v0⇔{A1X⩾A1v01TA1X⩽1TA1v0A2X=A2v0\left\{ \begin{matrix} {A_1}X = {A_1}{v_0} + y\\ {A_2}X = {A_2}{v_0} \end{matrix} \right.\Leftrightarrow \left\{ \begin{matrix} {A_1}X \geqslant {A_1}{v_0} \\ {1^T}{A_1}X \leqslant {1^T}{A_1}{v_0}\\{A_2}X={A_2}{v_0} \end{matrix} \right. {A1X=A1v0+yA2X=A2v0⇔⎩⎨⎧A1X⩾A1v01TA1X⩽1TA1v0A2X=A2v0
注意到XXX为单纯形CCC中的任意一点,故上述两个不等式和一个等式构成了一个多面体,即命题得证。
凸优化中:单纯形是一种多面体的证明相关推荐
- 浅谈凸优化中的共轭函数
浅谈凸优化中的共轭函数 函数ff的共轭定义: f∗(y)=sup(yTx−f(x))f^*(y) = \sup (y^Tx - f(x)), x∈domf{x\in {\bf dom} f} 可见,共 ...
- 凸优化中如何改进GD方法以防止陷入局部最优解
在对函数进行凸优化时,如果使用导数的方法(如:梯度下降法/GD,牛顿法等)来寻找最优解,有可能陷入到局部最优解而非全局最优解. 为了防止得到局部最优,可以对梯度下降法进行一些改进,防止陷入局部最优. ...
- 凸优化中凸函数定义、直线与线段、凸集、仿射集合、仿射函数
凸优化(Stephen Boyd)中自学部分 凸函数定义 函数f:Rn→Rf:\mathbf{R}^{n} \rightarrow \mathbf{R}f:Rn→R是凸的,如果fff在定义域(domd ...
- 如何理解凸优化中的共轭函数的定义?
共轭函数的意义主要就是:一个函数即便不是凸函数,但通过共轭法获得一个凸函数,很方便求解全局最优解的问题. 另外,共轭函数亦称对偶函数.极化函数,函数的某种对偶变换.过多的东西我就不再赘述了. 此处我是 ...
- 机器学习中的凸优化理论
凸优化课程 优化问题 定义:从一个可行解中找到一个最好的元素. 通常来说优化问题都可以写成如下的形式: 最小化目标函数: m个约束函数: 凸规划与非凸规划 凸规划满足: 凸规划都是相对容易解决的,非凸 ...
- 凸优化问题中的Bregman迭代算法
对于搞图像处理的人而言,不懂变分法,基本上,就没法读懂图像处理的一些经典文献.当然,这已经是10年之前的事情了. 现在,如果不懂得Bregman迭代算法,也就没法读懂最近几年以来发表的图像处理的前沿论 ...
- 文献翻译__人工智能时代医学图像重建中的凸优化算法(第4、5、6章)
文章下载–我的Gitee Convex optimization algorithms in medical image reconstruction-in the age of AI 人工智能时代医 ...
- 凸优化学习笔记 10:凸优化问题
前面讲了那么多关于凸集.凸函数的知识,然而都是铺垫,现在我们才来到了这门课的重头戏部分--凸优化问题! 文章目录 1. 一般优化问题 2. 凸优化问题 2.1 凸优化问题定义 2.2 凸优化问题的最优 ...
- 凸优化一(基本概念介绍)
背景知识 优化问题 凸优化是一类数学优化问题,介绍凸优化前先简单介绍数学优化问题. 优化问题定义(形式): minimizef0minimizef0 minimize \quad f_{0} subj ...
- 数学之美:凸优化问题
导言 凸优化(convex optimization)是最优化问题中非常重要的一类,也是被研究的很透彻的一类.对于机器学习来说,如果要优化的问题被证明是凸优化问题,则说明此问题可以被比较好的解决.在本 ...
最新文章
- 锋利的jQuery-4--动画方法总结简表
- 企业运维监控平台架构设计与实现
- JVM 的 工作原理,层次结构 以及 GC工作原理
- Algorithm:C++语言实现之分治法相关问题(给定实数x和整数n,分治法求xn)
- 四十八、微信小程序开发系统组件
- html霓虹灯效果图,HTML5 SVG+CSS3霓虹灯文字边框动画特效
- 【读书笔记】--- 《码出高效:java开发手册》
- windows 7/2008系统字体变成黑体解决
- 如何使用JMeter 对Dubbo接口进行测试
- 使用多线程与不使用多线程访问同一个网址实例对比
- POJ 1472 Instant Complexity
- matlab图像处理教学视频,MATLAB图像处理实例详解视频教程
- 如何才能快速拥有一张高含金量的IE证书?
- 2021年12月电子学会青少年软件编程(图形化)等级考试试卷(四级)答案解析
- (六)CRAFT----2019CVPR论文解读
- PowerDesigner一键导出数据库设计表结构
- PDJCAD皮带机设计软件
- mysql保存微信昵称特殊字符
- layui 数字步进器_光音移动设计规范 — 表单类
- Java-JDK下载过慢的问题解决方案
热门文章
- 计算机能否代替老师英语作文,雅思大作文范文:电脑不可取代老师
- 云计算计算机二级,2021计算机二级office用的是哪个版本
- 2022 基于SpringBoot的高校食堂点评系统 大众点评系统
- Spring的bean管理(注解注入属性)
- Scikit-Learn (浅谈Kmeans聚类算法)
- 计算机音乐tfboys手机,TFBoys王源的iPhone挂了 别怕 手贱有得治
- 下两个网段转发的路由设置_一台路由器下如何实现多个网络互通
- python的http请求和应答_python通过get,post方式发送http请求和接收http响应的方法,pythonget...
- NYOJ 2括号配对 栈的运用
- python做数据可视化的优势_用Python进行数据可视化的10种方法