凸优化(四)凸函数分析

1. 概述

\quad之前简单介绍了凸函数的定义,相信大家对凸函数有了简单的认识,但是这是远远不够的,这次通过一些详细的函数讲解来介绍一下部分常见凸函数的特点。

2. 凸函数的四个定义:

(1)第一个定义:如果X为在实数向量空间的凸集。并且有映射f:X\rightarrow R,如果f被称为,则有\forall x_1,x_2\in X,\forall t\in[0,1]: f(tx_1+(1-t)x_2)\leq tf(x_1)+(1-t)f(x_2).如果F被称为严格凸,那么有:\forall x_1\ne x_2\in X,\forall t\in(0,1): f(tx_1+(1-t)x_2)\ngeq tf(x_1)+(1-t)f(x_2).

(2)第二个定义:有映射f:R^n\rightarrow R为凸\leftrightarrow dom f为凸$,$\forall x,y\in dom f,\forall v为任意方向,g(t)=f(x+tv)为凸dom g=\{t|x+tv\in dom f\}

(3)第三个定义:若f可微,对\forall x,y\in dom f,f(y)\geq f(x)+\nabla^Tf(y-x)

(4)第四个定义二阶条件,若f:R^n\rightarrow R二阶可微,则f为凸$\leftrightarrow dom f为凸,\nabla^2f(x)\geq0,\forall x\in dom f(这里的大于等于号是表示特征值大于等0,表示矩阵半正定) 。

这四个定义在不同地方均有用处,但在判断函数是否为凸函数时最常用的是第四个。其中\nabla^2f(x)Hessian矩阵,表示函数的二阶偏导矩阵。

3. 常见凸函数:

(1)仿射函数:f(x)=Ax+b,显然,其二阶导函数为\nabla^2f(x)=0,所以仿射函数为凸函数

(2)指数函数:f(x)=e^{ax},x\in R,显然f^{'}(x)=ae^{ax},f^{''}(x)=a^2e^{ax}\geq0,所以指数函数是凸函数

(3)幂函数:f(x)=x^a,x\in R_{++},接着求导啊求导~,f^{'}(x)=ax^{a-1}f^{''}(x)\begin{cases} \geq 0,a\geq1 或a\leq0 (凸函数)\\ \leq 0,0\leq a\leq1 (凹函数)\end{cases},显然啦,当a=1、0时,幂函数就成为了仿射函数,所以即凸又凹。

(4)负熵函数:f(x)=xln{x},x\in R_{++},还是求导,f^{'}(x)=lnx+1,f^{''}(x)=\frac{1}{x}\nleq0,嗯,还是个严格凸函数。(也是个非常重要的函数!!)

(5)极大值函数(重中之重):现在来一个比较复杂却非常常见的函数:f(x) = max\{x_1,x_2,...,x_n\}这个函数显然是不可导的,那么首先根据定义一来看一下是否为凸函数。取两值x,y\in dom f,构造凸组合的新值\theta x+(1-\theta)y,f(\theta x+(1-\theta)y)=max\{\theta x_i+(1-\theta)y_i\}\leq \theta max\{x_i,i=0,1..n\}+(1-\theta)max\{y_i,i=0,1..n\}=\theta f(x)+(1-\theta)f(y),发现满足凸函数定义,所以极大值函数时凸函数。但是啊,由于其无法求导,后续处理会出现各种问题。所以,这里有一个解析逼近,就是用一个解析函数去逼近极大值函数。这个函数是这样的log-sum-exp:f(x)=log(e^{x_1}+...+e^{x_n}),x\in R^n,且有性质max\{x_i\}\leq f(x)\leq max\{x_i\}+ln(n) \quad 那么来证明一下这个函数也是凸函数吧!这里就要用到凸函数的第四个定义了,轮到Hessian矩阵出场了。对上述函数求二次偏导得到如下关系(公式打得累死):
\frac{\partial^2f}{\partial{x_i}\partial{y_i}}=\frac{-e^{x_i}e^{x_j}}{(e^{x_1}+...e^{x_n})^2},i\neq j, \frac{\partial^2f}{\partial{x_i}\partial{y_i}}=\frac{-e^{x_i}e^{x_i}+e^{x_i}(e^{x_1}+...e^{x_n})}{(e^{x_1}+...e^{x_n})^2},i=j\tag{1}
这个式子看上去也很丑,那么定义列向量z=[e^{x_1},...,e^{x_n}]^T,那么(1)式就变成了\frac{\partial^2f}{\partial{x_i}\partial{y_i}}=\frac{-e^{x_i}e^{x_j}}{(1^Tz)^2},i\neq j, \frac{\partial^2f}{\partial{x_i}\partial{y_i}}=\frac{-e^{x_i}e^{x_i}+e^{x_i}(1^Tz)} {(1^Tz)^2},i=j,函数的Hessian矩阵可以写成H(x)=\frac{1}{(1^Tz)^2}((1^Tz)diag(z)-zz^T),另后项矩阵为K,diag(z),为向量z展开的对角矩阵那么大家还记得半正定矩阵如何证明么?就是\forall x\in R^n,x^TAx\geq0成立,那么A则为半正定矩阵。好,那么开始构造!!V^TKV=(1^Tz)V^Tdiag(z)v-V^Tzz^TV=(\sum_iz_i)\sum_i(V_i^2z_i)-(\sum_iv_iz_i)^2\tag{2}a_i=v_i\sqrt{z_i},b_i=\sqrt{z_i},那么(2)式就变成了:(b^Tb)(a^Ta)-(a^Tb)^2\geq0\tag{3}此式成立,用到的性质为柯西-施瓦茨不等式,所以log-sum-exp函数为凸函数。

(6)行列式的对数:f(X)=\log det(X),dom f=S^n_{++},首先说明一下啊,当矩阵X只有一维时,那么原函数则为\log x,显然是凹函数。所以我们是在已经知道其为凹函数的前提下证明它是凹函数的了~根据凸函数的第二个定义当\forall z\in S^n_{++},\forall t\in R^{n \times n},z+tv\in S^n_{++}=dom f,构造凸组合的函数g(t)=f(z+tv)=\log det(z+tv)=\log det(z^{\frac{1}{2}}(I+tv^{-\frac{1}{2}}vz^{-\frac{1}{2}})z^{\frac{1}{2}})\tag{4}继续化简得到为:\log det\{z\}+\log det\{I+tz^{-\frac{1}{2}}vz^{-\frac{1}{2}}\}=\log det\{z\}+\sum^n_i\log(1+t\lambda_i),其中\lambda_i\geq0\tag{5}接着只要分析这个式子就可以,求导即可,得到:g^{'}(t)=\sum_i\frac{\lambda_i}{1+t\lambda_i},g^{''}(t)=\sum_i\frac{-\lambda_i^2}{(1+t\lambda_i)^2}\leq0\tag{6}到这里证明就结束了,原函数为凹函数得证。

4. 总结:

\quad可见啊,分析函数凸性一般都是通过其Hessian矩阵来分析,但是对于部分凸函数的证明也不是简单的,总体的计算过程也在越来越复杂,后面会逐步讲解凸问题的理论与求解。但是在证明的过程中会发现,其理论也是一步一步建立起来的,弄懂了原理之后看问题就会举一反三了。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 212,816评论 6 492
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 90,729评论 3 385
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 158,300评论 0 348
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 56,780评论 1 285
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 65,890评论 6 385
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,084评论 1 291
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,151评论 3 410
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 37,912评论 0 268
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,355评论 1 303
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,666评论 2 327
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 38,809评论 1 341
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,504评论 4 334
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,150评论 3 317
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,882评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,121评论 1 267
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,628评论 2 362
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,724评论 2 351

推荐阅读更多精彩内容