数值分析matlab方法.docVIP

塞德尔迭代法（matlab）： A=[5 -1 -1 -1; -1 10 -1 -1; -1 -1 5 -1; -1 -1 -1 10]; b = [-4; 12; 8; 34]; x0=[0;0;0;0]; d0=diag(A); d=diag(d0); l=tril(A,-1); r=triu(A,1); b0=inv(d+l); E=-b0*r; F=b0*b; X1=E*x0+F; K=1; while norm(X1-x0,inf)0x0=X1; X1=E*x0+F K=K+1 End 松弛因子迭代法：(注意w的变化对收敛的影响) A=[5 -1 -1 -1; -1 10 -1 -1; -1 -1 5 -1; -1 -1 -1 10]; b = [-4; 12; 8; 34]; x0=[0;0;0;0]; d0=diag(A); d=diag(d0); l=tril(A,-1); r=triu(A,1);w=1.4; b0=inv(d+w*l); E=b0*[(1-w)*d-w*r]; F=w*b0*b; X1=E*x0+F; K=1; while norm(X1-x0,inf)0.001 x0=X1; X1=E*x0+F K=K+1 end 冥法： A=[2,4,6;3,9,15;4,16,36]