用牛顿法matlab程序参数,牛顿法MATLAB程序

______________________________________________________________________________________________________________ 精品资料 问题一

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 是迭代次数.