資源描述:
《雅可比迭代法的MATLAB程序.doc》由會員上傳分享,免費在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、雅可比迭代法的MATLAB程序:Function[x,k,index]=Jacobi(A,b,ep,it-max)%求線性方程組的雅可比法;%A為方程組的系數(shù)矩陣;%b為方程組的右端項;%x為方程組的解;%ep為精度要求,缺省值為le-5;%it_max為最大迭代次數(shù),缺省值為100;%k為迭代次數(shù);%index為指標變量,index=0表示計算失敗,index=1表示計算成功;ifnargin<4it_max=100;endifnargin<3ep=le-5;endn=length(A);k=0;x=zeros(n,1);y=zeros(n,1);index=1;whilek<=it
2、_maxfori=1:nifabs(A(i,i))3、代次數(shù);%index為指標變量,index=0表示計算失敗,index=1表示計算成功;ifnargin<4it_max=100;endifnargin<3ep=le-5;endn=length(A);k=0;x=zeros(n,1);y=zeros(n,1);index=1;whilek<=it_maxfori=1:nifabs(A(i,i))4、i)=(b(i)-A(i,1:i-1)*y(1:i-1)-A(i,i+1:n)*x(i+1:n)/A(i,i);endendifnorm(y-x,inf)