共轭梯度法--MATLAB程序4页.doc

%conjugate gradient methods%method:FR,PRP,HS,DY,CD,WYL,LS%精确线搜索,梯度终止准则function m,k,d,a,X,g1,fv = conjgradme( G,b,c,X,e,method)if nargin6 error(输入参数必须为6);endn=length(G);if n=2format long e %ratsyms x1 x2f=1/2*x1,x2*G*x1;x2+b*x1;x2+c;g=diff(f,x1);diff(f,x2);g1=subs(subs(g,x1,X(1,1),x2,X(2,1);d=-g1;a=-(d*g1)/(d*G*d);% a=-(X(:,1)*G*d+b*d)/(d*G*d); a=g1(:,1)*g1(:,1)/(d(:,1)*G*d(:,1);X(:,2)=X(:,1)+a*d;g1=g1 subs(subs(g