雅克比和高斯赛德迭代法.docxVIP

第六次作业：分别编写用Jacobi迭代法和Gauss-Seidel迭代法求解线性方程组Ax=B的标准程序，并求下列方程组的解。可取初始向量 ，迭代终止条件Jacobi迭代法function [ x,k] =jacobi(A,b)% A 为方程组的系数矩阵% b 为方程组的右端项% x 为方程组的解% N为迭代次数[n,m]=size(A);n=length(b);e=10^-6;x1=zeros(1,n);x=zeros(1,n);M=100;y=zeros(1,n);r=max(abs(b));k=0;while refor i=1:n sum=0;for j=1:nif j~=i sum=sum+A(i,j)*x1(j);endend y(i)=(b(i)-sum)/A(i,i);end r=max(abs(y-x1)); x1=y;k=k+1;if kM warning(不收敛);return;endendx=y;Gauss-Seidel迭代法function [ x,k] =gss1(A,b)% A 为方程组的系数矩阵% b 为方程组的右端项% x 为方程组的解% N为迭代次数[n,m]=size(A);n=length(b);e=10^-6;x1=zeros(1,n);x=zeros(1,n);M=1