基于各种迭代算法的无约束优化

[toc]

一、牛顿法

1.1 原理及公式

  • 原理:对原始高维泰勒二阶展开用二阶近似,利用满足一阶必要条件解得一阶的极小点0 = \nabla q(x) = g^{(k)} + F(x^{(k)})(x - x^{(k)}),类似于梯度法可得迭代公式。
  • x^{(k + 1)} = x^{(k)} - F(x^{(k)})^{-1} g^{(k)}

1.2 牛顿法的特点:

1.缺点 需要求Hesse矩阵的逆,也就是求解一个n维非齐次方程组,也可能黑塞矩阵不可逆;
当初始点原理目标函数的极小值时,即使F(x^{(k + 1)})>0也可能收敛的很慢。
如果黑塞矩阵非正定,牛顿法确定的方向也不一定是下降方向。

  1. 优点 牛顿法可以求解多元方程组g(x)=0,这个时候F(x)表示的就是函数在x处的雅可比矩阵;
    牛顿法可以至少以阶数为二的收敛率收敛到x^*
    对于正定二次型问题,牛顿法只需要一步就能走到最小点。

1.3牛顿修正法

  • 加步长
  • Levenberg-Marquardt修正,将Hesse矩阵变成F(x^{(k)} + \mu_k I)^{-1},使其变成正定的。

二、共轭法

Prev

Gram-Schmid向量组共轭方法,从一组已知的线性无关向量组构造共轭向量。
d^{(k + 1)} = p^{(k + 1)} - \Sigma_{i = 0}^{k}\frac{p^{(i + 1)T}Qd^{(i)}}{d^{(i)T}Qd^{(i)}} d^{(i)}

2.1 共轭法的基本概念、原理和公式

2.1.1 定义

共轭定义:设Q是对称矩阵,对于所有方向\{d^{(i)}\}(i=0,1,...,m),如果有i\ne j,都有d^{(i)}Qd^{(j)} = 0,则称他们是关于Q共轭的。
性质:如果Q是n阶对称正定矩阵,如果k个方向非零且关于Q共轭,那么它们是线性无关的。
定理\forall k, 0\leq k \leq n, 0 \leq i \leq k, g^{(k + 1)T}d^{(i)}=0,该性质对最速下降法均成立。
原理:由于n个方向是提前给定的,每一步的公式中只需要求步长即可,那么步长的求法与最速下降相同,只需要求argmin_{\alpha}f(x^{(k)} + \alpha d^{(k)}),解得\alpha = -\frac{g^{(k) T}d^{(k)}}{d^{(k)T}Qd^{(k)}}

2.2 共轭法的特点

  1. 收敛性质:对于n维二次型问题,能够在n步之内得到结果
  2. 不需要计算Hesse矩阵也就不需要求逆了。

2.3 共轭梯度法

解二次型问题:

2.3.1 原理及公式

原理:每一次都利用上一次搜索方向和目标函数在当前迭代点的梯度向量之间的线性组合构造一个新方向,得到第k+1次方向:d^{(k + 1)} = -g^{(k + 1)} + \beta_k d^{(k)},使其与前面已经产生的搜索方向组成Q共轭方向。由定义可以求出:\beta_k = \frac{g^{(k + 1)T}Qd^{(k)}}{ d^{(k)T} Q d^{(k )}}
性质:共轭梯度法的搜索方向是关于Q的共轭方向。

2.3.2解非二次型问题:

利用泰勒展开式的二阶近似,F(x)对应于共轭梯度法的Q矩阵,但是要求函数可微,且每次都需要计算黑塞矩阵,运算量可能很大。

2.3.3 实现注意事项

  • n次迭代不一定到极小点,一般经过n次或n+1次迭代之后将搜索方向重新初始化为梯度方向
  • 一维搜索精度:是影响共轭梯度法的关键因素

三、拟牛顿法

思想:构造一个正定矩H_k来近似牛顿法中的F(x)^{-1},避免了求逆的计算量。
迭代公式
d^{(k)}= - H_k g^{(k)}
\alpha_k = argmin_\alpha f(x^{(k)} + \alpha d^{(k)})
x^{(k + 1)} = x^{(k)} - \alpha H_k g^{(k)}
正确性推导f(x^{(k + 1)}) = f(x^{(k)}) + g^{(k)T}(x^{(k + 1)} - x^{(k)}) + o(||x^{(k + 1)} - x^{(k)}||)
拟牛顿条件:由泰勒公式易得:H_{k + 1}(g^{(k + 1)}-g^{(k)}) = ( x^{(k + 1)} - x^{(k)})
秩一修正公式H_{k + 1} = H_k + \alpha_k z^{(k)}z^{(k)T}

H_k正定的条件是构造出来的右边分母大于零。
性质:H_{k + 1} \Delta g^{(i)} = \Delta x^{(i)}

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

推荐阅读更多精彩内容