数值分析 jacobi迭代法求解线性方程组 MATLAB程序实现

jacobi迭代法 参考数值分析第四版 颜庆津著 P36 

数值分析 jacobi迭代法MATLAB程序实现   运行输入为：

运行结果为：

以下是函数内容（保存为jacobi.m文件，在MATLAB中运行）：

%function [G,d,x,N]=jacobi(A,b)

%Jacobi迭代法求解线性方程组。迭代公式x(k+1)=Gx(k)+d,(k=0,1,…).
%返回迭代矩阵G，d和方程组的解x，以及迭代次数N


%2015.11.6  密密编写  (*^__^*) 嘻嘻……
function [G,d,x2,N]=jacobi(A,b)
n=size(A);
if n(1)~=n(2)
    error('矩阵A不是方阵');
end

%初始化
N=1;%迭代次数
%X=zeros(N,n)%存储每次迭代解出的x，迭代次数N
B=zeros(n);%用B表示L+U
D=zeros(n);
G=zeros(n);
x=ones(n,1);%生成全1的n维列向量
d=ones(n,1);
B=A;
for i=1:n
    B(i,i)=0;
    D(i,i)=A(i,i);
end
%迭代开始
G=-inv(D)*B;
if vrho(G)>1
    error('jacobi迭代法不收敛！');
end
d=inv(D)*b;
x1=x;
x2=G*x+d;
while norm(x2-x1,inf)>10^(-4)
    x1=x2;
    x2=G*x2+d;
    N=N+1;
end
x=x2;
end