数值分析上机作业3.docVIP

实验一 实验步骤： 1、 得到 则有： Jacobi迭代法： 令则称为雅克比迭代矩阵 可得雅克比迭代的迭代格式如下： Gauss-Seidel迭代法： 令，则称为Gauss-Seidel迭代矩阵 可得Gauss-Seidel迭代的迭代格式如下： SOR迭代法： 令,则有： 令带入的值可有 称为SOR迭代矩阵 可得SOR迭代的迭代格式如下： 2、编写程序： 一、Jacobi迭代法M文件： function [y,n]=Jacobi(A,b,x0,eps) n=length(A); if nargin3 error(输入错误，最少要输入三个参数); return; end if nargin==3 eps=1e-6; end D=diag(diag(A)); L=-tril(A,-1); U=-triu(A,1); M=D; N=L+U; B=M\N; f=M\b; if max(abs(eig(B)))=1 disp(谱半径大于等于1，迭代不收敛，无法进行); return; end n=1; for i=1:1:n if sum(A(i,i)~=n)~=n error(输入的A矩阵的对角线元素不能为0); return; end end y=B*x0+f; while norm(y-x0)=epsn100 x0=