richason.m

matlab 多个计算方程租的算法汇集 主要讲的是一般方法 还有正定矩阵的解法

function [x,n]=richason(A,b,x0,eps,M)
%Richardson法求解线性方程组 Ax=b
%方程组系数矩阵:A
%方程组之常数向量:b
%迭代初始向量:X0
%e解的精度控制:eps
%迭代步数控制：M
%返回值线性方程组的解：x
%返回值迭代步数：n
if(nargin == 3)
    eps = 1.0e-6;
      M = 200;
elseif(nargin == 4)
      M = 200;
end

I =eye(size(A));
x1=x0;
x=(I-A)*x0+b;
n=1;

while(norm(x-x1)>eps)
    x1=x;
    x=(I-A)*x1+b;
    n = n + 1;
    if(n>=M)
        disp('Warning: 迭代次数太多，现在退出！');
        return;
    end
end