www.pudn.com > gauss_Seidel.rar > gauss_Seidel.m

% 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