4峭弟、非線性共軛梯度法的研究

??上節(jié)我們研究了線性共軛梯度法,線性共軛梯度法的研究對象是二次函數(shù)脱拼,且采取的線搜索為精確線搜索瞒瘸。為此可以產(chǎn)生共軛向量組,具有二次終止性熄浓。所謂的二次終止性情臭,并不是迭代兩次就終止,而是對于二次函數(shù)且采取精確線搜索能夠有限步終止「┰冢基于二次函數(shù)的良好性質(zhì)丁侄,我們將推廣到一般函數(shù),采用一般線搜索朝巫。實(shí)際計(jì)算中,發(fā)現(xiàn)方法是有效的石景。便有了非線性共軛梯度法劈猿,在不引起混淆的情況下,非線性共軛梯度法也被稱為線性共軛梯度法潮孽。
??對共軛梯度法的研究主要集中在參數(shù)~\beta_k~的選擇揪荣,混合共軛梯度法,多項(xiàng)共軛梯度法和譜共軛梯度法等方面往史。

1仗颈、前言

??共軛梯度法是無約束優(yōu)化方法,主要解決如下問題
\min_{x\mathbb{R}^n}~f(x)\tag{1}
解決問題 (1)椎例,我們采用是線搜索的迭代方法挨决,即
x_{k+1}=x_k+\alpha_k d_k\tag{2}
其中~d_k~是搜索方向,~\alpha_k~是搜索步長订歪,無論是混合共軛梯度法脖祈,譜共軛梯度法或者是多項(xiàng)共軛梯度法,只是方向~d_k~不同刷晋。

2盖高、經(jīng)典共軛參數(shù)~\beta_k~的選擇

??一般地,共軛梯度法的搜索方向?yàn)?br> d_k=\begin{cases}-g_k,&k=1,\\ -g_k+\beta_k d_{k-1},&k\ge2. \end{cases}
1952 年眼虱,Hestenes 和 stiefel^{[1]}在線性共軛梯度法中提出
\beta_k^{HS}=\frac{g_k^T(g_k-g_{k-1})}{d_{k-1}^T(g_k-g_{k-1})}
1964 年喻奥,F(xiàn)letcher 和 Reeves^{[2]}首次提出了非線性共軛梯度法
\beta_k^{FR}=\frac{\Vert g_k\Vert^2}{\Vert g_{k-1}\Vert^2}
1969 年,Polak , Ribiere^{[3]} 和 Polyak^{[4]} 提出
\beta_k^{PRP}=\frac{g_k^T(g_k-g_{k-1})}{\Vert g_{k-1}\Vert^2}
1987 年捏悬,F(xiàn)letcher^{[5]} 提出
\beta_k^{CD}=\frac{\Vert g_k\Vert^2}{-g_{k-1}^Td_{k-1}}
1991 年撞蚕,Liu 和 Storey^{[6]} 提出
\beta_k^{LS}=\frac{g_k^T(g_k-g_{k-1})}{-g_{k-1}^Td_{k-1}}
1998 年,戴彧虹 和 袁亞湘^{[7]} 提出
\beta_k^{DY}=\frac{\Vert g_k\Vert^2}{d_{k-1}^T(g_k-g_{k-1})}
2001年邮破,戴彧虹 和 Liao^{[8]} 提出
\beta_k^{DL}=\frac{g_k^T(g_k-g_{k-1}-ts_{k-1})}{d_{k-1}^T(g_k-g_{k-1})}
其中~t> 0~诈豌,s_{k-1}=x_k-x_{k-1}.
2005 年,Hager 和 Zhang^{[9]} 提出
\beta_k^{HZ}=\frac{g_k^T(y_{k-1}-2d_{k-1}\frac{\Vert y_{k-1}\Vert^2}{d_{k-1}^Ty_{k-1}})}{d_{k-1}^Ty_{k-1}}
其中~y_{k-1}=g_k-g_{k-1}~

??以上是八種經(jīng)典的共軛梯度法抒和,其收斂性會在后面詳細(xì)介紹矫渔。

3、混合共軛梯度法

??為了利用各種基本共軛梯度法的不同優(yōu)點(diǎn)摧莽,許多學(xué)者采用了不同共軛梯度法的巧妙組合庙洼。

Gilbert 和 Nocedal^{[10]}為保證算法的收斂性和具有較好的數(shù)值表現(xiàn),取
\beta_k=\begin{cases} &-\beta_k^{FR},&~若~\beta_k^{PRP}<-\beta_k^{FR}\\ &\beta_k^{PRP},&~若~\vert \beta_k^{PRP}\vert\le\beta_k^{FR}\\ &\beta_k^{FR},&~若~\beta_k^{PRP}>\beta_k^{FR} \end{cases}
戴雨虹 和 袁亞湘^{[11]} 提出了混合 DY 和 CD 共軛梯度法
\beta_k=\frac{\Vert g_k\Vert^2}{\max \left\{d_{k-1}^T(g_k-g_{k-1}),-g_{k-1}^Td_{k-1} \right\}}
焦寶聰,陳蘭平 和 潘翠英^{[12]} 提出混合 DY 和 FR 共軛梯度法
\beta_k=\begin{cases} &\beta_k^{DY},~&若~g_k^T d_{k-1}\ge \Vert g_{k-1}\Vert^2\\ &\beta_k^{FR},~&否則 \end{cases}
??以上只是列出幾種混合梯度法而已油够,具體他們有什么性質(zhì)蚁袭,收斂性的證明,后面會有更加全面的介紹石咬。

4揩悄、多項(xiàng)項(xiàng)共軛梯度法

?? 基本的共軛梯度法是負(fù)梯度方向與前一搜索方向的組合,許多學(xué)者在此基礎(chǔ)上鬼悠,研究了負(fù)梯度删性、前一搜索方向或位移、梯度差的各種形式焕窝,得到了多項(xiàng)共軛梯度法蹬挺。多項(xiàng)共軛梯度法中最主要的形式還是三項(xiàng)共軛梯度法。

2006 年它掂,張麗巴帮,周偉軍,李董輝^{[13]}提出了改進(jìn)的 PRP 共軛梯度法虐秋,得到了如下的三項(xiàng)共軛梯度法
d_k=\begin{cases} &-g_k,&~k=1,\\ &-g_k+\beta_k^{PRP}d_{k-1}-\frac{g_k^T d_{k-1}}{\Vert g_{k-1}\Vert^2}y_{k-1},&~k\ge 2. \end{cases}
2011 年榕茧,Narushima,Yabe 和 Ford^{[14]}得到了一般的三項(xiàng)共軛梯度法
d_k=\begin{cases} &-g_k,&k=1,或~g_k^Tp_k=0,\\ &-g_k+\beta_k d_{k-1}-\beta_k\frac{g_k^T d_{k-1}}{g_k^Tp_k}p_k,&其他情形熟妓, \end{cases}
其中~p_k~為任意向量
同年雪猪,Andrei^{[15]} 將 PRP 公式改進(jìn)為
d_k=-\frac{y_{k-1}^Ts_{k-1}}{\Vert g_{k-1}\Vert^2}+\frac{y_{k-1}^Tg_k}{\Vert g_{k-1}\Vert^2}s_{k-1}-\frac{g_k^Ts_{k-1}}{\Vert g_{k-1}\Vert^2}y_{k-1}
其中~s_{k-1}=x_k-x_{k-1},~y_{k-1}=g_k-g_{k-1}~.

5、譜共軛梯度法

??譜共軛梯度法是由譜梯度法和共軛梯度法發(fā)展而來起愈。譜梯度法又稱 BB 算法前域,最早是由 Barzilai 和 Borwein^{[16]} 于 1988 年為求解無約束優(yōu)化問題 (1) 提出來的镜粤。BB 方法的主要思想是在最小二乘意義下搀矫,生成能夠逼近目標(biāo)函數(shù) Hesse 矩陣的逆矩陣刑峡,其迭代具有以下形式
x_{k+1}=x_k-\alpha_k d_k
其中
\alpha_k=\frac{(x_k-x_{k-1})^T(x_k-x_{k-1})}{(x_k-x_{k-1})^T(g_k-g_{k-1})}
BB 方法可以看成是最速下降法的改進(jìn),優(yōu)點(diǎn)是它的數(shù)值表現(xiàn)遠(yuǎn)遠(yuǎn)好于最速下降法阐污。
??2001 年休涤,Birgin 和 Martinez^{[17]} 將譜梯度和共軛梯度相結(jié)合,提出了譜共軛梯度法笛辟,其搜索方向如下
d_k=\begin{cases} &-g_k,&~k=1,\\ &-\theta_kg_k+\beta_k d_{k-1},&~k\ge 2,\tag{1} \end{cases}
其中
\theta_k=\frac{s_{k-1}^Ts_{k-1}}{s_{k-1}^Ty_{k-1}}
\beta_k=\frac{g_k^T(\theta_ky_{k-1}-s_{k-1})}{d_{k-1}^Ty_{k-1}}
y_{k-1}=g_k-g_{k-1},~~s_{k-1}=x_k-x_{k-1}
我們把 (1) 的方法稱為譜共軛梯度法功氨,但是上式的譜共軛梯度法不能保證是下降算法。
??張麗, 周偉軍^{[18]}提出一種譜共軛梯度法
d_k=\begin{cases} &-g_k,&k=1,\\ &-(1+\beta_k\frac{g_k^T d_k}{\Vert g_k\Vert^2})g_k+\beta_k d_{k-1},&k\ge 2, \end{cases}

6手幢、結(jié)束語

\color{red}{以上的內(nèi)容只是簡單的介紹共軛梯度法捷凄,后面將會對共軛梯度法作進(jìn)一步學(xué)習(xí)}

7围来、參考文獻(xiàn)

[1] Hestenes M R, Stiefel E. Method of conjugate geadient for linear equations[J]. Research of the National Bureau of Standards, 1952, 49(6): 409-436
[2] Flecher R, Reeves C. Fuction minimization by conjugate gradient gradients[J]. Computer Journal, 1964 7(2): 149-154.
[3] Polak E, Ribiere G. Note surla convergence de directions conjugate[J]. Rev Fr Inform Rech Oper, 1969, 16(3): 35-43.
[4] Polyak B T. The conjugate gradient method in extreme problems. USSR Comput Math Phys, 1969, 9(1): 94-112.
[5] Fletcher R. Practical methods of optimization, vol I: unconstrained optimization[M]. New York: John Wiley and Sons, 1987.
[6] Liu Y, Storey C. Efficient generalized conjugate gradient algorithms, I. Theory[J]. J Optim Theorey Appl, 1991, 69(1): 129-137.
[7] Dai Y H, Yuan Y X. A nonlinear conjugate gradient method with a strong global convergence property. SIAM J Optim, 1999, 10(1): 177-182.
[8] Dai Y H, Liao L Z. New conjugacy conditions and related nonlinear conjugate gradient methods[J]. Applied Mathematics and Optimization, 2001, 43 : 87-101.
[9] Hager W W, Zhang H C. A new conjugate gradient method with guaranteed descent and an efficient line search, SIAM Journal on Optimization, 2005, 1(16) : 170-192.
[10] Gilbert J C, Nocedal J. Global convergence proerties of conjugate gradient methods for optimization[J]. SIAM Journal on Optimization, 1992, 2, 21-42.
[11] Dai Y H, Yuan Y X. Some properties of a new conjugate gradient methods[J]. Advances in Nonlinear Programming. 1998, 12, 251-262.
[12] 焦寶聰跺涤,陳蘭平匈睁,潘翠英. Goldstein 線搜索下混合共軛梯度法的全局收斂性[J]. 計(jì)算數(shù)學(xué), 2007,2(29): 137-146.
[13] Zhang L, Zhou W J, LI D H. A descent modified Polak_Ribiere-Polak gradient method and its global convergence[J]. IMA Journal of Numerical Analysis, 2006, 26: 629-640.
[14] Yasushi N, Hiroshi Y, John A F. A three-term conjugate gradient method with sufficient descent property for unconstrined optimization[J]. 2011,
[15] Andrei N. Amodified Polak-Ribiere-Polyak conjugate gradient algorithm for unconstrained optimization. Optimization, 60(12), 1457-1471.
[16] Barzilai J, Borwein M J. Two-point step size gradient methods[J]. IMA Journal of Numerical Analysis. 1988, 1(8): 141-148.
[17] Birgin E G, Martinez J M. A spectural conjugate gradeint method for unconstrained optimization[J]. 2001, 2(42): 117-128.
[18] Zhang L, Zhou W J. Two descent hybrid conjugate gradient methods for optimization[J]. Journal of Computation and Applied Mathematics, 2008, 251-264.

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市桶错,隨后出現(xiàn)的幾起案子航唆,更是在濱河造成了極大的恐慌,老刑警劉巖院刁,帶你破解...
    沈念sama閱讀 216,470評論 6 501
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件糯钙,死亡現(xiàn)場離奇詭異,居然都是意外死亡退腥,警方通過查閱死者的電腦和手機(jī)超营,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,393評論 3 392
  • 文/潘曉璐 我一進(jìn)店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來阅虫,“玉大人,你說我怎么就攤上這事不跟⊥堑郏” “怎么了?”我有些...
    開封第一講書人閱讀 162,577評論 0 353
  • 文/不壞的土叔 我叫張陵窝革,是天一觀的道長购城。 經(jīng)常有香客問我,道長虐译,這世上最難降的妖魔是什么瘪板? 我笑而不...
    開封第一講書人閱讀 58,176評論 1 292
  • 正文 為了忘掉前任,我火速辦了婚禮漆诽,結(jié)果婚禮上侮攀,老公的妹妹穿的比我還像新娘。我一直安慰自己厢拭,他們只是感情好兰英,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,189評論 6 388
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著供鸠,像睡著了一般畦贸。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上楞捂,一...
    開封第一講書人閱讀 51,155評論 1 299
  • 那天薄坏,我揣著相機(jī)與錄音,去河邊找鬼寨闹。 笑死胶坠,一個胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的鼻忠。 我是一名探鬼主播涵但,決...
    沈念sama閱讀 40,041評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼杈绸,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了矮瘟?” 一聲冷哼從身側(cè)響起瞳脓,我...
    開封第一講書人閱讀 38,903評論 0 274
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎澈侠,沒想到半個月后劫侧,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 45,319評論 1 310
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡哨啃,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,539評論 2 332
  • 正文 我和宋清朗相戀三年烧栋,在試婚紗的時候發(fā)現(xiàn)自己被綠了。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片拳球。...
    茶點(diǎn)故事閱讀 39,703評論 1 348
  • 序言:一個原本活蹦亂跳的男人離奇死亡审姓,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出祝峻,到底是詐尸還是另有隱情魔吐,我是刑警寧澤,帶...
    沈念sama閱讀 35,417評論 5 343
  • 正文 年R本政府宣布莱找,位于F島的核電站酬姆,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏奥溺。R本人自食惡果不足惜辞色,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,013評論 3 325
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望浮定。 院中可真熱鬧相满,春花似錦、人聲如沸桦卒。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,664評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽闸盔。三九已至悯辙,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間迎吵,已是汗流浹背躲撰。 一陣腳步聲響...
    開封第一講書人閱讀 32,818評論 1 269
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留击费,地道東北人拢蛋。 一個月前我還...
    沈念sama閱讀 47,711評論 2 368
  • 正文 我出身青樓,卻偏偏與公主長得像蔫巩,于是被迫代替她去往敵國和親谆棱。 傳聞我的和親對象是個殘疾皇子快压,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,601評論 2 353

推薦閱讀更多精彩內(nèi)容