cggm.m

matlab编程的三个线性方程组求解子程序：cggm.m,fdm.m,febs.m.其中cggm.m为共轭梯度法求解子程序,fdm.m为最速下降法求解自程序,febs.m为三对角追赶法求解子程序.

function x=cggm(A,b,x0,tol)
%x=cggm(A,b,x0,tol)
%线性方程组共轭梯度法求解程序。
%时间：二零零七年五月二十八日
%地点：南华大学
%作者：王礼广
%A为系数矩阵，b为右向量,x0为迭代初值,tol为误差容许

xk=x0;
rk=b-A*xk;
pk=rk;
d=norm(rk);
while d>tol
    uk=A*pk;
    ak=(rk'*rk)/(pk'*uk);
    xk1=xk+ak*pk;
    rk1=rk-ak*uk;
    bk=(rk1'*rk1)/(rk'*rk);  
    pk1=rk1+bk*rk;  
    d=norm(xk1-xk);
    rk=rk1;
    pk=pk1;
    xk=xk1;
end
x=xk;
return;