雅可比迭代法的MATLAB程序.docxVIP

雅可比迭代法的 MATLAB 程序： Function[x，k,index] =Jacobi (A,b，ep,it—max) ％ 求线性方程组的雅可比法； % A 为方程组的系数矩阵; ％ b 为方程组的右端项； ％ x 为方程组的解； ％ ep 为精度要求，缺省值为 le-5； ％ it_max 为最大迭代次数,缺省值为 100； % k 为迭代次数； % index 为指标变量， index=0 表示计算失败,index=1 表示计算成功; if nargin〈4 it_max=100; end if nargin3 ep=le—5； end n=length(A) ;k=0;x=zeros (n ，1);y=zeros (n，1)；index=1; while k=it_max for i=1：n if abs (A(i,i)) le— 10 index=0； return; end y(i)=(b(i) —A (i，1:n)*x(1:n)+A(i,i) *x (i))/A (i,i)； end if norm (y-x，inf) ep break； end k=k+1； x=y； end 高斯—赛德尔迭代的 MATLAB 程序 Function [x,k，index]=Gau-seidel(A，b,ep,it—max) ％ 求线性方程组的高斯-赛德尔迭代法； ％ A 为方程组的系数矩阵； % b 为方程组的右端项; ％ x 为方程组的解； % ep 为精度要求，缺省值为 le—5； ％ it_max 为最大迭代次数,缺省值为 100; % k 为迭代次数; ％ index 为指标变量,index=0 表示计算失败,index=1 表示计算成功； if nargin4 it_max=100； end if nargin3 ep=le-5； end n=length(A);k=0;x=zeros(n,1)；y=zeros (n ，1)；index=1; while k〈=it_max for i=1：n if abs (A (i,i)) le— 10 index=0; return; End if i==1 y (i)= (b(i)-A(i，i+1 ：n)＊x (i+1:n)/A (i，i)； elseif i==n y (i)=(b (i) -A(i，1 ：i- 1) *y(1：i— 1) /A(i，i) ; else y (i)= (b(i)—A(i,1 ：i— 1)＊y(1:i— 1) -A (i，i+1:n)*x (i+1:n) /A(i，i)； end end if norm (y-x,inf)〈ep break； end k=k+1; x=y； end