共轭梯度法求线性方程组的解matlab程序.docx

《共轭梯度法求线性方程组的解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)，即可求出线性方程组的解。