gauss_seidel.m

用matlab实现的Gauss-Seidel迭代法。

% Created by colent 12/7
% Copyright (C) 2007-2008 XMU 数值分析 Experiment. All rights reserved.\nAuthors: 廖绮绮(ID:04120049).
% Thanks 鞠颖老师
% Gauss-Seidel迭代法

function[] = gauss_Seidel()
% data input
%B = [10 -1 -2;-1 10 -2;-1 -1 5]
%b = [72 83 42]'

B = importdata('data_B.txt')
b = importdata('data_bb.txt')

pause;  
dim = max(size(B));   % 矩阵A的行、列数的最大值
x = zeros(dim ,1);     % 生成dim维全0矩阵
maxCount = 6;              % 最大容许迭代次数
for j = 2:maxCount
    x(1) = (b(1)-B(1 , 2:dim)*x(2 : dim))/B(1,1);
    for i = 2:dim-1
        x(i) = (b(i)-B(i,1:i-1)*x(1:i-1) - B(i , i+1:dim)*x(i+1 : dim))/B(i , i);
    end
        x(dim) = (b(dim)-B(dim , 1:dim-1)*x(1 : dim-1))/B(dim,dim);
    count = j-1
    x = x
    pause;
end