matlab數值分析 線性方程組的迭代解法?高斯-賽德爾迭代法

?

Function [x,iter]=gs(A,b,tol) D=diag(diag(A)); L=D-tril(A); U=D-triu(A); x=zeros(size(b)); %從x=[0;0…]T開始 for iter=1:500x=(D-L)\(b+U*x); %此句換為x=(D)\(b+L*x+U*x);即為Jacobi迭代error=norm(b-A*x)/norm(b);if(error<tol)break;end end

主函數(調用程序)

A=[2,-1,0;-1,3,-1;0,-1,2];b=[1;8;-5];tol=1e-4;[x,iter]=gs(A,b,tol)

?

?

總結

以上是生活随笔為你收集整理的2021-01-07 matlab数值分析 线性方程组的迭代解法 高斯-赛德尔迭代法的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。