______________________________________________________________________________________________________________ 精品资料 问题一
1、问题描述
本次作业中使用共轭梯度法求解232122212142min x x x x x x x +-++-,初始点取为
T x )2,2,2(0=。
2、求解方法及求解程序
(1)共轭梯度法:
对于二次函数的无约束最小化问题:x b x Q x x f T T -=2
1)(,共轭梯度方法: k k k k d x x α+=+1
其中步长k α由最小线性化准则确定:)(min )(k k k k k d x f d x ααα
+=+。 既将函数)(x f 对α求导等于零,得到的α就是我们当前所求步长。 梯度方向:b Qx x f g k k k -=∇=)(
共轭梯度的方向由下式生成:00g d -=
1,...,1,1-=+-=-n k d g d k k k k β
其中k β由下式给出:
1)1(--=k k k
k k g g g g T T β
该方法在最多n 次迭代后,将终止于某个最优解处。
(2)求解程序
frcg.m
function [x,val,k]=frcg(fun,funs,x0)
% 功能: 用FR 共轭梯度法求解无约束问题: min f(x)
%输入: x0是初始点, fun, gfun 分别是目标函数和梯度
%输出: x, val 分别是近似最优点和最优值, k 是迭代次数.