共轭梯度法求线性方程组的解matlab程序(总1页).docx

精选优质文档-倾情为你奉上一、函数代码：function x=cg(a,b,c,x0) 注释: ax=b;c表示误差精度;x0表示初始向量,c和x0均可省略,但如果c省略了,则x0必须省略。 n=length(a);if nargin4 x0=zeros(n,1);endif nargin3 c=10(-3);endr0=b-a*x0;r=r0;d=r0;m=zeros(n,1);p=zeros(n,1);x=x0;for k=0:n-1 m(k+1)=r*r/(d*a*d); x=x+m(k+1)*d; r=b-a*x; if norm(r,inf)=c | k+1=n break; end p(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)，即可求出线性方