共轭梯度法的matlab代码_精品文档.txt资料文档下载
《共轭梯度法的matlab代码_精品文档.txt资料文档下载》由会员分享,可在线阅读,更多相关《共轭梯度法的matlab代码_精品文档.txt资料文档下载(1页珍藏版)》请在冰豆网上搜索。
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
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);
break
x,f=subs(f,{xi,yi},x),count
x=subs(x,a,ai);
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;