数值代数第五章上机作业.docxVIP

第五章 共轭梯度法程序： function [ output,k ] = congra( A,b,x0 ) %%%%共轭梯度法eps=1e-32; fprintf(共轭梯度法：epslon=%f \n,eps); x=x0; r=b-A*x; k=0; del=norm(r,2); while deleps k=k+1; fprintf(第%4d 次迭代，del(norm(r,2))为：%f\n,k,del); if k==1 p=r; else end  bel=(r*r)/(rk*rk); p=r+bel*p; end alp=(r*r)/(p*A*p); x=x+alp*p; rk=