matlab 矩阵jocobi迭代_雅可比迭代法的MATLAB程序 Myth丶恋晨 2023-01-09 01:06 345阅读 0赞 雅可比迭代法的 MA TLAB 程序: 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 nargin<3 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)) 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) 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 nargin<4 it\_max=100;
还没有评论,来说两句吧...