最优化—基于matlab的共轭梯度法(共3页).docx

精选优质文档-倾情为你奉上共轭梯度算法程序：function f=conjugate_gradient(x0,t)x=x0;syms xi yi af=(xi-2)2+(yi-4)2;fx=diff(f,xi);fy=diff(f,yi);fx=subs(fx,xi,yi,x0);fy=subs(fy,xi,yi,x0);fi=fx,fy;n=0;while double(sqrt(fx2+fy2)t s=-fi; if n=0 s=-fi; else s=s1; end x=x+a*s; f=subs(f,xi,yi,x); f1=diff(f); f1=solve(f1); if fi=0 ai=double(f1); else break x,f=subs(f,xi,yi,x),n end