Question

我为Jacobi方法编写了以下函数，需要对其进行修改才能执行Gauss-Seidel

function [ x,iter] = jacobi( A,b,tol,maxit )
%jacobi iterations
%   

x=zeros(size(b));

[n,m]=size(A);

i=1;

iter=maxit;

for i=1:maxit

    for j=1:n

        y(j)=(b(j)-A(j,1:j-1)*x(1:j-1)-A(j,j+1:n)*x(j+1:n))/A(j,j)

    end

    if max(abs(A*y'-b))<tol
        iter=i;
        break;
    end
    x=y';

end

我知道我需要让x(1:j-1)更新，但不确定如何编写它，谢谢

Answer 1

你只需要摆脱y并用x替换任何出现的y。

for j=1:n

    x(j)=(b(j)-A(j,1:j-1)*x(1:j-1)-A(j,j+1:n)*x(j+1:n))/A(j,j)

end

if max(abs(A*x-b))<tol
    iter=i;
    break;
end

Jacobi从旧计算新的向量，然后立即替换所有变量。

Gauß-Seidel就地计算并始终使用最新值。

Answer 2

在Gauss-Seidel中，将迭代循环更改为：

for j = 1:n
    jj = 1:n; 
    jj(j) = [];  
    Xtemp = x;  
    Xtemp(j) = [];  
    y(j) = (b(j) - sum(A(j,jj) * Xtemp)) / A(j,j);
end

Answer 3

function [x succes iter]=gausssedel(A,b,x0,tol,maxiter)

n=length(A);
succes=0;
iter=maxiter;
x=zeros(n,1);

while maxiter > 0
    maxiter=maxiter-1;
    for i=1:n
       % suma=A(i,1:i-1)*x(i,1-i)+A(i,i+1:n)*x0(i+1:n);
         suma=A(i,1:i-1)*x(1:i-1)+A(i,i+1:n)*x0(i+1:n);
        x(i) = (b(i)-suma)/A(i,i);
    end
    if norm(x-x0) < tol
        succes=1;
        break;
    end
    x0=x;

end
iter=iter-maxiter;
end

Jacobi迭代到Gauss-Seidel

3 个答案: