最优化共轭梯度法课件.ppt

1、共轭的非零共轭的非零个个是是阶对称正定矩阵，阶对称正定矩阵，是是设设AkdddnAk,21性性无无关关。向向量量，则则这这个个向向量量组组线线. 1定理定理证明证明，使得，使得设存在实数设存在实数k ,21,01 kiiid ，则有，则有上式两边同时左乘上式两边同时左乘AdTj,01 kiiTjiAdd 可化简为可化简为共轭的向量，所以上式共轭的向量，所以上式个个是是因为因为Akdddk,21.0 jTjjAdd ，是正定矩阵，所以是正定矩阵，所以而而因为因为0,0 jTjjAddAd所以所以。kjj,2,1,0 线性无关。线性无关。因此因此kddd,21几何意义几何意义设有二次函数设有二次函

2、数)()(21)(xxAxxxfT 对称正定矩阵，对称正定矩阵，是是其中其中nnA 是一个定点。是一个定点。x的等值面的等值面则函数则函数)(xfcxxAxxT )()(21为中心的椭球面。为中心的椭球面。是以是以 x由于由于,0)()( xxAxf,0)(2 AxfA所以所以正定，正定，因为因为的极小点。的极小点。是是因此因此)(xfxx,)(2Axf 而而点，点，是在某个等值面上的一是在某个等值面上的一设设)0(x处的法向量为处的法向量为该等值面在点该等值面在点)1(x. )()()1()1(xxAxf o1x2xx)1(d)0(x中的一个方向，中的一个方向，是是nRd)1(。以最优步长搜

3、索得到点以最优步长搜索得到点沿着沿着)1()1()0(xdx所所在在等等值值面面的的切切向向量量。是是点点则则)1()1(xd正交，正交，与与则则)()1()1(xfd , 0)()1()1( xfdT即即,)1()2(xxd 令令)1(x所以所以, 0)2()1( AddT共轭。共轭。小点的向量关于小点的向量关于向量与由这一点指向极向量与由这一点指向极即等值面上一点处的切即等值面上一点处的切A)2(d. 2定理定理，设有函数设有函数cxbAxxxfTT 21)(共轭向量。共轭向量。一组一组是是阶对称正定矩阵。阶对称正定矩阵。是是其中其中AdddnAk)()2()1(,进行搜索，进行搜索，为初

4、始点，依次沿为初始点，依次沿以任意的以任意的)()2()1()1(,kndddRx 上的上的在在是函数是函数则则得到点得到点kkkBxxfxxxx )1()1()1()3()2()(,极小点，其中极小点，其中,|1)(RdxxBikiiik 是是时，时，当当，生成的子空间。特别地生成的子空间。特别地是由是由)1()()2()1(, nkxnkddd上上的的唯唯一一极极小小点点。在在nRxf)(推论推论有有在上述定理条件下，必在上述定理条件下，必。kidxfiTk,2,1,0)()()1( 共轭方向法对于极小化问题对于极小化问题:法法为为共共轭轭方方向向法法是是正正定定矩矩阵阵，称称下下述述算算

5、其其中中 A,21)(mincxbAxxxfTT ；共轭方向共轭方向取定一组取定一组)()2()1(,)1(ndddA,)2()1()()1( kkxxx确定点确定点依次按照下式由依次按照下式由任取初始点任取初始点 )(min)()()()()()()()1(kkkkkkkkkdxfdxfdxx 。满足满足直到某个直到某个0)()()( kkxfx注注至多经过至多经过求解上述极小化问题，求解上述极小化问题，可知，利用共轭方向法可知，利用共轭方向法由定理由定理2。次次迭迭代代必必可可得得到到最最优优解解n2. 共轭梯度法 如何选取一组共轭方向？如何选取一组共轭方向？:共轭梯度法共轭梯度法eeve

6、sRFletcher 代点代点向相结合，利用已知迭向相结合，利用已知迭将共轭性和最速下降方将共轭性和最速下降方基本思想：基本思想：进行搜索，求出进行搜索，求出共轭方向，并沿此方向共轭方向，并沿此方向处的梯度方向构造一组处的梯度方向构造一组函数的极小点。函数的极小点。以下分析算法的具体步骤。以下分析算法的具体步骤。cxbAxxxfTT 21)(min是常数。是常数。，是对称正定矩阵，是对称正定矩阵，其中其中cRbARxnn ,；，第一个搜索方向取为，第一个搜索方向取为任取初始点任取初始点)()1()1()1()1(xfdx ，令令，若若，设已求得点设已求得点)(0)()2()1(1)1()1(

7、kkkkxfgxfx:)1(按如下方式确定按如下方式确定则下一个搜索方向则下一个搜索方向 kd)1()(1)1(kkkkdgd 令令共轭。共轭。关于关于和和要求要求Addkk)()1( ？如何确定如何确定k ，得，得）式两边同时左乘）式两边同时左乘则在（则在（AdTk)(1)()(1)()1()(0kTkkkTkkTkdAdAgdAdd )2()()(1)(kTkkTkkdAdgAd 解得解得:)3(搜搜索索步步长长的的确确定定，步长步长利用一维搜索确定最优利用一维搜索确定最优和搜索方向和搜索方向已知迭代点已知迭代点kkkdx ,)()(。即求解即求解)(min)()(kkdxf , )()(

8、)()(kkdxf 记记，令令0)()()()()( kTkkddxf ，即有即有0)()()()( kTkkdbdxA ，则有，则有令令bAxxfgkkk )()()(，0)()( kTkkdAdg )3()()()(kTkkTkkAdddg 解得解得3定理定理次次算法在算法在，对于正定二次函数对于正定二次函数nmFRcxbAxxxfTT 21)(），下下列列关关系系成成立立（且且对对所所有有的的一一维维搜搜索索后后即即终终止止，并并mii 1;1,2,1,0)1()()( ijAddjTi;1,2,1,0)2( ijggjTi。iTiiTiggdg )()3(注注共共轭轭的的。是是可可知知

9、搜搜索索方方向向）由由定定理理（Adddm)()2()1(,31则则构构造造的的搜搜索索必必须须取取负负梯梯度度方方向向，否否算算法法中中第第一一个个搜搜索索方方向向)2(方方向向不不能能保保证证共共轭轭性性。）可知，）可知，的（的（由定理由定理33)3(,0|2)( iiTiiTigggdg处的下降方向。处的下降方向。是迭代点是迭代点所以所以)()(iixd的的计计算算公公式式可可以以简简化化。算算法法中中，由由定定理理iFR 3)4()()(1)(iTiiTiiAddgAd )()()(1iTiiTiAdddAg / )( / )( )()1()()()1(1iiiTiiiiTixxAdx

10、xAg .)()()(bxAxfgiii )()(1)(11iiTiiiTiiggdggg iTiigdg)(21| )4(|221iigg 算算法法步步骤骤：FR。，令，令精度要求精度要求，任取初始点任取初始点1. 1)1( kx 为所求极小点；为所求极小点；停止，停止，若若令令)1(1)1(1,|, )(. 2xgxfg 。令令，）计算）计算利用公式（利用公式（否则，令否则，令)1(1)1()2(11)1(3,dxxgd 为所求极小点；为所求极小点；停止，停止，若若令令)1(1)1(1,|, )(. 3 kkkkxgxfg ）计算。）计算。用公式（用公式（其中其中否则，令否则，令4,)(1

11、)1(kkkkkdgd 。转转，令，令）计算）计算利用公式（利用公式（3,3. 4)()()1(kkkkkdxx 。令令1: kk算算法法求求解解下下述述问问题题：用用例例FR22212)(minxxxf 。初始点取为初始点取为Tx)2,2()1( 解：解：.)2,4()(21Txxxf 次迭代：次迭代：第第1,)4,8(1)1(Tgd 令令)1()1()1(11AdddgTT ,2004),(21)(2121 xxxxxf.2004 A而而 482004)4,8(48)4,8(185 )1(1)1()2(dxx 所以所以TT)4,8(185)2,2( T)98,92( 次迭代：次迭代：第第

12、2.)916,98(2Tg 21221|gg .81448)916()98(2222 )1(12)2(dgd TT)4,8(814)916,98( T)4,1(8140 )2()2()2(22AdddgTT 412004)4,1()8140(41)916,98(81402209 )2(2)2()3(dxx TT)4,1(8140209)98,92( T)0,0( Tg)0,0(3 即为所求极小点。即为所求极小点。)3(x3. 用于一般函数的共轭梯度法nRxtsxf .)(min共轭梯度法进行修改：共轭梯度法进行修改：对用于正定二次函数的对用于正定二次函数的确确定定。）计计算算，需需由由一一维维

13、搜搜索索不不能能利利用用公公式式（搜搜索索步步长长3)2(i 。速下降方向，即速下降方向，即第一个搜索方向仍取最第一个搜索方向仍取最)()1()1()1(xfd 算算：其其它它搜搜索索方方向向按按下下式式计计,)()()1()1(iiiidxfd 。其中其中2)(2)1(|)(|)(|iiixfxf 此此时时可可采采取取一一定定能能满满足足停停止止条条件件，算算法法在在有有限限步步迭迭代代后后不不)3(如如下下措措施施：没没有有求求成成一一轮轮搜搜索索后后，如如果果还还次次迭迭代代为为一一轮轮，每每次次完完以以n新的初始点，新的初始点，的最后一个迭代点作为的最后一个迭代点作为得极小点，则以上一轮得极小点，则以上一轮一轮搜索。一轮搜索。一个搜索方向，开始下一个搜索方向，开始下取最速下降方向作为第取最速下降方向作为第注注，如如算算采采用用其其它它形形式式的的公公式式计计在在共共轭轭梯梯度度法法中中，也也可可i 。）共轭梯度法共轭梯度法PRPgggggiTiiiTii()(11 1( )|(iii TigDixondg 共轭梯度法） 。