SOR算法解线性方程组的matlab程序(共4页).docx

编写用SOR方法求解线性方程组Ax=B的标准程序，并求下列方程组的解，并比较松弛因子取1.0、1.25、1.5时所需迭代的次数。可取初始向量 ，迭代终止条件function x,k =sor(A,b,w)% A 为方程组的系数矩阵% b 为方程组的右端项% x 为方程组的解% N为迭代次数n,m=size(A);n=length(b);e=10-8;x1=ones(1,n);x=zeros(1,n);M=100;s=zeros(1,n)y=zeros(1,n);r=max(abs(b);k=0;while re for i=1:n sum=0; for j=1:n if ji sum=sum+A(i,j)*x1(j); end end s(i)=w*(b(i)-sum)/A(i,i);