9787.0)(618.02=-+=a b a x ,0199.0)(22-==x f f ; 因为ε=<=-5.0438.0a b ,所以算法停止,得到
927.02
146
.1708.02*=+=+=
b a x 。 例 6 用FR 共轭梯度法求解问题()2212min 2f x x x =+,要求选取初始点05,5x ??= ???
610ε-=。
解 ???
?
??=?2142)(x x x f ,???? ??=20100g ,???? ??--=-=201000g d , 2
200)205(2)105(205105)(ααααα-+-=???
? ??--=+f d x f ,
令05001800)(00=-=+αααd x f d d ,则1850=α,于是??
????
??-=+=959200001d x x α;
则?
????
? ??-=?=920940)(11x f g ,95201=g ,8140
0110==g g g g T T β,??????
??-=+-=81100814000011d g d β, 2
211)9
201(8150)9201(81400)(ααα+-+-=+d x f , 令
0)(11=+d x f d d αα,则209
1=α,于是???
? ??=+=001112d x x α; 则???? ??=?=00)(22x f g ,ε<=02g ,故???
? ??==002*x x 为所求。 例7 用外罚函数法求解:
解 即
于是
令
)()01.1min 222
21≥-+-=x t s x x x f ()()222
122,(1)min 0,1P x M x x M x =-++-????
()()22122222
1
222(1),
10,(1)1,10x x x P x M x x M x x ?-+-≥?=?-++--?112(1)
P x x ?=-?()2222222,1221,1x x P x M x x x ≥??=?
+-?02
1=??=??x P x P