《matlab高斯-赛德尔迭代程序》由会员分享,可在线阅读,更多相关《matlab高斯-赛德尔迭代程序(2页珍藏版)》请在金锄头文库上搜索。
matlab中应用的高斯-赛德尔迭代程序主程序如下:function X=gsdddy(A,b,X0,P,wucha,max1)D=diag(diag(A);U=-triu(A,1);L=-tril(A, -1);dD=det(D);if dD=0 disp(请注意:因为对角阵D奇异,所以此方程无解)else disp(请注意:因为对角阵距D非奇异,所以此方程有解) iD=inv(D -L);B2=iD*U;f2=iD*b;jX=Ab; X=X0;n m=size(A); for k=1:max1 X1=B2*X+f2;djwcX=norm(X1 -X,P); xdwcX=djwcX/(norm(X,P)+eps); if(djwcXwucha)|(xdwcXwucha) return else k;X1;k=k+1;X=X1; end end if(djwcXwucha)|(xdwcX A=10 3 1;2 -10 3;1 3 10; b=14;11;20;X0=0 0 0; X=gsdddy(A,b,X0,inf,0.001,100)请注意:因为对角矩阵D非奇异,所以此方程组有解。X = 1.2820 -0.2592 1.9496