共轭梯度法的matlab代码_精品文档.txt

共轭梯度法的matlab代码_精品文档.txt

《共轭梯度法的matlab代码_精品文档.txt》由会员分享，可在线阅读，更多相关《共轭梯度法的matlab代码_精品文档.txt（1页珍藏版）》请在冰豆网上搜索。

functionf=conjugate_grad_2d（x0,t）

%pleaseinputthis:

conjugate_grad_2d（[2,2],0.05）

x=x0;

symsxiyia

f=xi^2-xi*yi+3*yi^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;

whiledouble（sqrt（fx^2+fy^2））>t

s=-fi;

ifcount<=0

s=-fi;

else

s=s1;

end

x=x+a*s;

f=subs（f,{xi,yi},x）;

f1=diff（f）;

f1=solve（f1）;

iff1~=0

ai=double（f1）;

else

break

x,f=subs（f,{xi,yi},x）,count

end

x=subs（x,a,ai）;

f=xi^2-xi*yi+3*yi^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;

count=count+1;

fx=fxi;

fy=fyi;

end

x,f=subs（f,{xi,yi},x）,count