共轭梯度法的matlab代码_精品文档.txt
- 文档编号:1520038
- 上传时间:2022-10-22
- 格式:TXT
- 页数:1
- 大小:698B
共轭梯度法的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
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 共轭 梯度 matlab 代码 精品 文档