水水水水水水

2021/12/22 23:23:43

本文主要是介绍水水水水水水,对大家解决编程问题具有一定的参考价值,需要的程序猿们随着小编来一起学习吧!

function [ x , k , err , time ]= mjacobi ( A,b,x,tol ,max_it )
if nargin <5
   max_it =1000; 
end 
 if nargin <4
   tol=1.e-5;
 end 
 if nargin <3
  x = zeros ( length ( b ));
 end
 tic;
bnrm2= norm ( b );
r= b - A *x;
 err = norm ( r )/bnrm2;
 if ( err < tol )
 return ;
 end
 D = diag(diag ( A )); 
 for k =1:max_it
 x = D \(( D - A )*x+ b );
 r = b -A *x ;
 err = norm ( r )/bnrm2; 
 if ( err <= tol ) 
     break ; 
 end 
 end 
 time = toc ;
end



这篇关于水水水水水水的文章就介绍到这儿,希望我们推荐的文章对大家有所帮助,也希望大家多多支持为之网!


扫一扫关注最新编程教程