用迭代法和超松弛法求解线性方程组.docVIP

实验名称 高斯-赛德尔、超松弛迭代法求解线性方程组 实验时间 2012-5-30 姓名 班级 学号 成绩 一 实验目的 [1] 掌握Gauss顺序消去和选列主元消去解线性方程组基本原理和方法 [2] 用MATLAB编写程序实现Gauss顺序消去和选列主元消去解线性方程组 二 实验内容 1.编写MATLAB程序实现雅可比迭代法、高斯-赛德尔、超松弛迭代法求解线性方程组 2.编写MATLAB程序实现超松弛迭代法求解线性方程组 ， 二 程序代码 1.（1）雅可比迭代法： function [x,n]=jacobi(A,b,x0,eps) if nargin==3 eps= 1.0e-6; M = 1000; elseif nargin3 error return elseif nargin ==5 M = 1000; end D=diag(diag(A)); %求A的对角矩阵 L=-tril(A,-1); %求A的下三角阵 U=-triu(A,1); %求A的上三角阵 B=D\(L+U); f=D\b; x=B*x0+f; n=0; %迭代次数 while norm(x-x0)=eps x0=x; x=B*x0+f; n=n+1; if(n=M) disp(Warning: 迭代次数太多，可能不收敛！); return; end end （2）高斯赛德尔法： 1.function x=gaussMethod(A,b) %高斯列主元消去法，要求系数矩阵非奇异的， n = size(A,1); if abs(det(A))= 1e-8 error(系数矩阵是奇异的); return; end for k=1:n ak = max(abs(A(k:n,k))); index = find(A(:,k)==ak); if length(index) == 0 index = find(A(:,k)==-ak); end %交换列主元 temp = A(index,:); A(index,:) = A(k,:); A(k,:) = temp; temp = b(index);b(index) = b(k); b(k) = temp; %消元过程 for i=k+1:n m=A(i,k)/A(k,k); %消除列元素 A(i,k+1:n)=A(i,k+1:n)-m*A(k,k+1:n); b(i)=b(i)-m*b(k); end end %回代过程 x(n)=b(n)/A(n,n); for k=n-1:-1:1; x(k)=(b(k)-A(k,k+1:n)*x(k+1:n))/A(k,k); end x=x; end 2. 超松弛法： function [x,n]=SOR(A,b,x0,w,eps,M) if nargin==4 eps= 1.0e-6; M = 200; elseif nargin4 error return elseif nargin ==5 M = 200; end if(w=0 || w=2) error; return; end D=diag(diag(A)); %求A的对角矩阵 L=-tril(A,-1); %求A的下三角阵 U=-triu(A,1); %求A的上三角阵 B=inv(D-L*w)*((1-w)*D+w*U); f=w*inv((D-L*w))*b; x=B*x0+f; n=1; %迭代次数 while norm(x-x0)=eps x0=x; x =B*x0+f; n=n+1; if(n=M) disp(Warning: 迭代次数太多，可能不收敛！); return; end end 三 数据结果 A=[-0.98 -0.05 -0.02;-0.04 -0.9 0.07;-0.02 0.09