实验三 迭代法解线性方程组.docVIP

实验三 迭代法解线性方程组 实验目的 学会用Jacobi迭代法、Gauss-Seidel迭代法和超松驰迭代法求线性方程组解。 学会对迭代法做收敛性分析，研究求方程组解的最优迭代方法。 学会用共轭梯度法求线性方程组的解，研究共轭梯度法的计算效率。 实验要求 按照题目要求完成实验内容。 写出相应的Matlab程序。 给出实验结果。 对实验结果进行分析讨论。 写出相应的实验报告。 实验步骤 1、研究Jacobi迭代法求解线性方程组的方法和相应的收敛性。 （1）用Jacobi迭代法（Jacobi.m）求解线性方程组 （4.31） 取初始点，精度要求为。请给出满足精度要求的迭代次数和相应的计算结果。 function [x,k]=jc(a,b,x0,ep,max) n=length(a); k=0; if nargin5 max=500; end if nargin4 ep=1e-5; end if nargin3 x0=zeros(n,1); y=zeros(n,1); end x=x0;x0=x+2*ep; while norm(x0-x,inf)epkmax k=k+1;x0=x; for i=1:n y(i)=b(i); for j=1:n if j~=i y(i)=y(i)-a(i,j)*x0(j); end end if abs(a(i,i))1e-10||k==max warning(a(i,i) ì?D?); return; end y(i)=y(i)/a(i,i); end x=y; end end a=[11 -3 -2;-1 5 -3;-2 -12 19]; b=[3 6 -7]; [x,k]=jc(a,b) x = 0.999985953146656 1.999977859006902 0.999978180649185 k = 33 研究相应迭代矩阵的谱半径和Jacobi迭代的渐近收敛速度。 2、研究Gauss-Seidel迭代法求解线性方程的方法和相应的收敛性 用Gauss-Seidel迭代法求解线性方程组，取初始点，精度要求为，请给出满足精度要求的迭代次数和相应的计算结果。 function [x,k]=gsd(a,b,x0,ep,max) n=length(a); k=0; if nargin5 max=500; end if nargin4 ep=1e-5; end if nargin3 x0=zeros(n,1); end x=x0;x0=x+2*ep; while norm(x0-x,inf)epkmax k=k+1; x0=x; z=x; for i=1:n z(i)=b(i); for j=1:n if j~=i z(i)=z(i)-a(i,j)*x(j); end end if abs(a(i,i))1e-10||k==max warning(a(i,i) 太小); return end z(i)=z(i)/a(i,i); x(i)=z(i); end end if k==max warning(已达上限); end end 研究相应迭代矩阵的谱半径和Gauss-Seidel迭代的渐近收敛速度。 3、研究松弛因子对超松弛迭代法的影响 选取用超松弛（SOR）迭代法（sor.m)求解线性方程组（4.31），取初始点，精度要求为，列出不同松弛因子下的迭代次数。 function [x,k]=sor(a,b,om,ep,Nmax) n=length(b); k=0; if nargin6 Nmax=500; end if nargin5 ep=1e-5; end if nargin4 x0=zeros(n,1); end if nargin3 om=1.5; end x=x0; x0=x+2*ep;k=0; l=tri