Jacobi迭代和Gauss_Seidel迭代在MATLAB编程实现

最新推荐文章于 2022-10-12 18:14:43 发布

努力搬砖挣猫猫奶粉钱

最新推荐文章于 2022-10-12 18:14:43 发布

阅读量1.1k

点赞数 4

文章标签： matlab 线性代数算法

本文链接：https://blog.csdn.net/qq_40265311/article/details/109508411

版权

运用矩阵迭代的实现方法

摘要由CSDN通过智能技术生成

Jacobi迭代

function [ x ] = jacobi( A,b,ep,N,x0 )
         %   JACOBI 用矩阵形式迭代的过程
         %   A是系数矩阵 B是常数项 ep是精度要求 N是迭代次数 
 format long    %计算结果显示小数点后更多位数
 if nargin<3
     disp('输入参数数目不足3个')
     return
    else if nargin>=5
          M=N;  
   end
 end
    
D=diag(diag(A));
L=-tril(A,-1);
U=D-L-A;
B=D\(L+U);
f=D\b;
x=B*x0+f;    %第一次迭代
 disp(['第1次收敛结果为最终收敛结果：'])
 disp(x')
  n=2;
  while n<M
    x0=x;
    x=B*x0+f；

        if  norm(x-x0)<ep
        disp(['第',num2str(n),'次收敛结果为最终收敛结果：'])
        disp(x')
        break
        
    end
    disp(['第',num2str(n),