函数名 功能
rs 里查森迭代法求线性方程组Ax=b的解
crs 里查森参数迭代法求线性方程组Ax=b的解
grs 里查森迭代法求线性方程组Ax=b的解
jacobi 雅可比迭代法求线性方程组Ax=b的解
gauseidel 高斯-赛德尔迭代法求线性方程组Ax=b的解
SOR 超松弛迭代法求线性方程组Ax=b的解
SSOR 对称逐次超松弛迭代法求线性方程组Ax=b的解
JOR 雅可比超松弛迭代法求线性方程组Ax=b的解
twostep 两步迭代法求线性方程组Ax=b的解
fastdown 最速下降法求线性方程组Ax=b的解
conjgrad 共轭梯度法求线性方程组Ax=b的解
preconjgrad 预处理共轭梯度法求线性方程组Ax=b的解
BJ 块雅克比迭代法求线性方程组Ax=b的解
BGS 块高斯-赛德尔迭代法求线性方程组Ax=b的解
BSOR 块逐次超松弛迭代法求线性方程组Ax=b的解
function [x,N]= BGS(A,b,x0,d,eps,M)
if nargin==4
eps= 1.0e-6;
M = 10000;
elseif nargin<4
error
return
elseif nargin ==5
M = 10000;
end
NS = size(A);
n = NS(1,1);
bnum = length(d);
bs = ones(bnum,1);
for i=1:(bnum-1)
bs(i+1,1)=sum(d(1:i))+1;
%获得对角线上每个分块矩阵元素索引的起始值
end
DB = zeros(n,n);
for i=1:bnum
endb = bs(i,1)+d(i,1)-1;
DB(bs(i,1):endb,bs(i,1):endb)=A(bs(i,1):endb,bs(i,1):endb);
%求A的对角分块矩阵
end
LB = -tril(A-DB); %求A的下三角分块阵
UB = -triu(A-DB); %求A的上三角分块阵
N = 0;
tol = 1;
while tol>=eps
invDL = inv(DB-LB);
x = invDL*UB*x0+invDL*b; %块迭代公式
N = N+1;
tol = norm(x-x0);
x0 = x;
if(N>=M)
disp('Warning: 迭代次数太多,可能不收敛!');