fastdown.m

来自「matlab科学计算的几乎所有数学算法以及仿真。每一个m文件在我学习的时候都验证」· M 代码 · 共 17 行

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

r  = b-A*x0;
d  = dot(r,r)/dot(A*r,r);
x  = x0+d*r;
n=1;

while(norm(x-x0)>eps)
    x0 = x;
    r  = b-A*x0;
    d  = dot(r,r)/dot(A*r,r);
    x  = x0+d*r;
    n  = n + 1;
end