欢迎来到天天文库
浏览记录
ID:53314037
大小:12.94 KB
页数:1页
时间:2020-04-03
《共轭梯度法求线性方程组的解matlab程序.docx》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、一、函数代码:functionx=cg(a,b,c,x0)注释:ax=b;c表示误差精度;x0表示初始向量,c和x0均可省略,但如果c省略了,则x0必须省略。n=length(a);ifnargin<4x0=zeros(n,1);endifnargin<3c=10^(-3);endr0=b-a*x0;r=r0;d=r0;m=zeros(n,1);p=zeros(n,1);x=x0;fork=0:n-1m(k+1)=r'*r/(d'*a*d);x=x+m(k+1)*d;r=b-a*x;ifnorm(r,i
2、nf)<=c
3、
4、k+1==nbreak;endp(k+1)=norm(r)^2/norm(r0)^2;d=r+p(k+1)*d;r0=r;end二、在命令窗体中输入参数a、b、c、x0以及x=cg(a,b,c,x0),即可求出线性方程组的解。
此文档下载收益归作者所有