数值分析与算法作业4.docVIP

第四章上机题作业 徐清瑶仪器科学与技术学号：20130308040002 1 (P145)解：(1)用各方法求解方程组A*x=b的程序如下： 生成Hilbert矩阵的程序如下： %hilbert.m function A=hilbert (n) %生成hi 1 bert矩阵八 A二zeros (n); for i=l:n for j=l:n A(i, j)=l/(i+j-l); end end Jacobi迭代法求解方程组A*x=b的程序如下： %myjacobi.m function [x, k]=myjacobi(八，b，xO, N) y(^jaCobi迭代法解线性方程组A*x=b %说明：A为系数矩阵，b为右端向量，xO为初始向量， % N为最大迭代次数，x返回近似解 n二length(b); x=zeros(n, 1);k=0; %判断雅对比迭代法的收敛情况 D=diag(diag(A)); B=inv ⑼ *(D-A); d=eig(B); j=abs(d); m=max (j) if ml dispC对任意的初始值雅可比迭代法收敛’) else dispC对任意的初始值雅可比迭代法不收敛’) end %求解线性方程组的解和迭代次数 while k=N for i二l:n x(i) = (b(i)-A(i, [1:i-1, i+1:n])*xO([l: i-1, i+1 :n]))/A(i, i end if norm(x-xO, inf)0. 0001, break;end %0. 0001 为误差精度 xO=x;k=k+l; end if k〉N dispC雅可比迭代次数已经超过最大迭代次数N’) end 输入： A=hilbcrt(10); b=[l;1/2;1/3;1/4;1/5;1/6;1/7;1/8;1/9;1/10]; x0=zeros(10, 1); N=300; [x, k]=myjacobi (A, b, xO, N) 结果为： m = 7. 7798 对任意的初始值雅可比迭代法不收敛 雅可比迭代次数已经超过最大迭代次数N x = 1.0e+267 * 0. 7445 1.7419 4384 TOC \o 1-5 \h \z 9670 3864 7293 4.0157 2590 4686 4. 6512 k = 301 分析：由8=11^(0)*(04)矩阵的谱半径m =7. 7798〉1可知对任意的初始值雅 可比迭代法不收敛，且雅可比迭代次数己经超过最大迭代次数N=300。 Gauss-Seidel迭代法求解方程组A*x=b的程序如下： %mygauss.m function [x，k]=mygauss (A, b, xO, N) 9(^Gauss-Seidel迭代法求解线性方程组Ax=b %说明：A为系数矩阵，b为右端向量，xO为初始向量， % N为最大迭代次数，x返回近似解 n二length(b); x=zcros(n, 1);k=0; %通过判断A是否对称正定来判断Gauss-Seidel迭代法的收敛情况 %由题意知A为对称矩阵，故只需判断是否正定 dl=A(l, 1); for i=2:n d=det (A(l: i, l:i)); end if (dl0)(max(d)0) dispC对任意的初始值Gauss-Seidel迭代法收敛’) else dispC对任意的初始值Gauss-Seidel迪代法不收敛’) end %求解线性方程组的解和迭代次数 while k〈=N for i=l:n if i==l x(l) = (b(l)-A(l，2:n)*xO(2:n))/A(l, 1); elseif i==n x (n) = (b (n) -A (n, 1 :n-l)*x(l :n-l))/A(n, n); else x(i) = (b(i)-A(i, i+l:n)*xO(i+l:n))/A(i, i); end end if norm(x-xO, inf)〈0. 0001, break;end %0. 0001 为误差精度 xO=x;k=k+l; end if k〉N dispC Gauss-Seidel迭代次数己经超过最大迭代次数N’) end 输入： A=hilbert (10); b=[1;1/2;1/3;1/4;1/5;1/6;1/7;1/8;1/9;1/10]; x0=zeros(10, 1); N二300; [x, k] =mygauss (A, b, xO, N) 结果为： 对任意的初始值Gauss-Seidel逃代法收敛 x 二 1 0 0 0 0 0 0 0 0 0 k = 1 分析：由A为对称正定矩阵可知对任意的初始值Gauss-Seidel迭代法收敛, 且 G