数值分析matlab版第八章.docVIP

实验报告八 题目：解线性方程的迭代法 摘要： 1.描述Jacobi迭代法、Gauss-Seidel迭代法和SOR迭代法。 2编写Jacobi迭代法、Gauss-Seidel迭代法和SOR迭代法的M文件。 选取不同的初始向量)0(x及右端面项向量b给定迭代误差要求分别用编写Jacobi迭代法和Gauss-Seidel迭代法程序求解观察得到的序列是否收敛若收敛通过迭代次数分析计算结果并得出你的结论。 (2)用编写的SOR迭代法程序对于所选取的初始向量及右端面项向量b进行求解，松驰系数ω取1ω2的不同值 目的 ： 1．通过本实验加深对Jacobi迭代法、Gauss-Seidel迭代法的构造过程的理解； 2．能对上述三种迭代法提出正确的算法描述编程实现，进一步理解迭代法的改进过程；设计思想 自选线性方程组，编制一个程序，分别用Jacobi迭代法、Gauss-Seidel迭代法求解，比较两种迭代法收敛速度的快慢。 第八章 习题1 Jacobi迭代法 A=[5,2,1,-12;-1,4,2,20;2,-3,10,3] MAXTIME=50; eps=1E-5; [n,m]=size(A); x=zeros(n,1); y=zeros(n,1); k=0; disp(迭代过程X的值情况如下:) disp(X=); while 1 disp(x); for i=1:1:n s=0.0; for j=1:1:n if j~=i s=s+A(i,j)*x(j); end y(i)=(A(i,n+1)-s)/A(i,i); end end for i=1:1:n maxeps=max(0,abs(x(i)-y(i))); end if maxeps=eps for i=1:1:n x(i)=y(i); end return; end for i=1:1:n x(i)=y(i); y(i)=0.0; end k=k+1; if kMAXTIME error(超过最大迭代次数,退出); return; end end Untitled4 A = 5 2 1 -12 -1 4 2 20 2 -3 10 3 迭代过程X的值情况如下: X= 0 0 0 -2.4000 5.0000 0.3000 -4.4600 4.2500 2.2800 -4.5560 2.7450 2.4670 -3.9914 2.6275 2.0347 -3.8579 2.9848 1.8865 -3.9712 3.0922 1.9670 -4.0303 3.0237 2.0219 -4.0139 2.9815 2.0132 -3.9952 2.9900 1.9972 -3.9954 3.0026 1.9960 -4.0002 3.0031 1.9999 -4.0012 3.0000 2.0010 -4.0002 2.9992 2.0002 -3.9997 2.9998 1.9998 -3.9999 3.0002 1.9999 -4.0000 3.0001 2.0000 Gauss-Seidel迭代法 format long; A=[5,2,1,-12;-1,4,2,20;2,-3,10,3] [n,m]=size(A); Maxtime=50; Eps=10E-5; x=zeros(1,n); disp(x=); for k=1:Maxtime disp(x); for i=1:n s=0.0; for j=1:n if i~=j s=s+A(i,j)*x(j); end end x(i)=(A(i,n+1)-s)/A(i,i);