fastdown.m

rs 里查森迭代法求线性方程组 crs 里查森参数迭代法求线性方程组 grs 里查森迭代法求线性方程组 jacobi 雅可比迭代法求线性方程组 gauseidel 高斯-赛德尔迭代法求线性方程

function [x,n]= fastdown(A,b,x0,eps)
if(nargin == 3)
    eps = 1.0e-6;
end

x=x0;
n=0;
tol=1;
while(tol>eps)                       %以下过程可参考算法流程
    r  = b-A*x0;
    d  = dot(r,r)/dot(A*r,r);
    x  = x0+d*r;
    tol = norm(x-x0);
    x0 = x;
    n  = n + 1;
end