求高斯赛德尔迭代矩阵 高斯-塞德尔迭代法的迭代矩阵怎么求?
浏览量:2990
时间:2021-03-11 04:28:45
作者:admin
高斯-塞德尔迭代法的迭代矩阵怎么求?
在Jacobi迭代中,a=d-l-u,ax=B,DX=(Lu)x,BX=d^{-1}(Lu)x,d^{-1}B
所以B=d^{-1}(Lu)
在Gauss-Seidel迭代中,a=d-l-u,但是ax=B(d-l)x=UX,BX=(d-l)^{-1}UX(d-l)^{-1}B
所以B=(d-l)^{-1}u
%Gauss-Seidel函数%a是未知量的系数矩阵;B是未知量的系数矩阵方程的右常量列向量,x0取1(m,1)m为未知数,EPS为精度,如果不输入EPS,则默认为1.0e-6函数[x,n]=Gauss-Seidel(a,B,x0,EPS)如果nargin==3 EPS=1.0e-6elseif nargin<3 error returnendD=diag(diag(A))L=-tril(A,-1)U=-triu(A,1)G=(D-L)Uf=(D-L)bx=G*x0 fn=1而norm(x-x0)>=EPS x0=x x=G*x0 f n=n1
求高斯赛德尔迭代矩阵 gs迭代法例题 gs迭代法迭代矩阵计算例题
版权声明:本文内容由互联网用户自发贡献,本站不承担相关法律责任.如有侵权/违法内容,本站将立刻删除。