共轭梯度法的matlab源程序.docxVIP

共轭梯度法的 matlab 源程序 function 共轭梯度法的 matlab 源程序 function g=conjugate_grad_2d(x0,t)%please input g=conjugate_grad_2d(x0,t)%please input thisconjugate_grad_2d([2,2],0.05)x=x0;syms xi yi thisconjugate_grad_2d([2,2],0.05)x=x0;syms xi yi af=xi^2-xiyi+3yi^2;fx=diff(f,xi);fy=diff(f,yi);fx=subs(fx,{xi,yi af=xi^2-xiyi+3yi^2;fx=diff(f,xi);fy=diff(f,yi);fx=subs(fx,{xi,yi },x0);fy=subs(fy,{xi,yi},x0);fi=[fx,fy];count=0;while },x0);fy=subs(fy,{xi,yi},x0);fi=[fx,fy];count=0;while double(sqrt(fx^2+fy^2))ts=-fi;if double(sqrt(fx^2+fy^2))ts=-fi;if count=0s=-fi;elses=s1;endx=x+as;f=subs(f,{xi,yi},x);f1=dif count=0s=-fi;elses=s1;endx=x+as;f=subs(f,{xi,yi},x);f1=dif f(f);f1=solve(f1);if f(f);f1=solve(f1);if f1~=0ai=double(f1);elsebreakx,f=subs(f,{xi,yi},x),countend f1~=0ai=double(f1);elsebreakx,f=subs(f,{xi,yi},x),countend x=subs(x,a,ai);f=xi^2-xiyi+3yi^2;fxi=diff(f,xi);fyi=diff(f,yi);f x=subs(x,a,ai);f=xi^2-xiyi+3yi^2;fxi=diff(f,xi);fyi=diff(f,yi);f xi=subs(fxi,{xi,yi},x);fyi=subs(fyi,{xi,yi},x);fii=[fxi,fyi];d=(fxi^ xi=subs(fxi,{xi,yi},x);fyi=subs(fyi,{xi,yi},x);fii=[fxi,fyi];d=(fxi^ 2+fyi^2)(fx^2+fy^2);s1=-fii+ds;count=count+1;fx=fxi;fy 2+fyi^2)(fx^2+fy^2);s1=-fii+ds;count=count+1;fx=fxi;fy =fyi;endx,f=subs(f,{xi,yi},x),count =fyi;endx,f=subs(f,{xi,yi},x),count 史锋 的 共轭梯度法matlab 程序以函数 f=1-(1/(sqrt(2*pi)))*(exp((-(xi-3)^2+yi^2)/2)+0.6*exp((-(xi+3)^2+yi^2)/2)); 为例 求解最小值 function f=conjugate_gradient(x0,eps) x=x0; syms xi yi a f=1-(1/(sqrt(2*pi)))*(exp((-(xi-3)^2+yi^2)/2)+0.6*exp((-(xi+3)^2+yi^2)/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(fx^2+fy^2))eps 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 f1~=0 ai=double(f1); else break x,n end x=subs(x,a,ai); f=1-(1/(sqrt(2*pi)))*(exp((-(xi-3)^2+yi^2)/2)+0.6*exp((-(xi+3)^2+yi^2)/2)); fxi=diff(f,xi); fyi=diff(f,yi); fxi=subs(fxi,{xi,yi},x); fyi=subs(fyi,{xi,yi},x); fii=[fxi,fyi]; d=(fxi^2+fyi^2)/(fx^2+fy^2); s1=-fii+d*s; %搜