数值分析三次作业.docVIP

习题.pdf 第六章 解线性方程组的迭代法 解：由于A为对称矩阵，且a11,a22,a33均大于0，故： Gauss-Seidel方法收敛的充要条件是A为正定矩阵； Jacobi方法收敛的充要条件是A和2D-A均为正定矩阵。 对矩阵A，顺序主子式?1=30,?2=3002=60,?3=30-2021-212=10，故A为正定矩阵 2D-A=600040004-30-2021-212=30202-12-12 对2D-A矩阵，顺序主子式?1=30,?2=3002=60,?3=30202-12-12=10，故2D-A为正定矩阵 综上，对线性方程组Ax=b，若A=30-2021-212，则Gauss-Seidel方法和Jacobi方法均收敛。 由于Gauss-Seidel方法利用了更新值，其收敛速度更快。 解：取x(0)=0，SOR法迭代公式为 x1(k+1)=x1(k)+ω4(1-4x1(k)+x2(k))x2(k+1)=x2(k)+ω4(4+x1k+1-4x2k+x3k)x3(k+1)=x3(k)-ω4(3-x2(k+1)+4x3(k)) 1） 当ω=1.03,迭代到第6次时， x(6)=(0.500000280744456,1.000000091599815,-0.499999984294099) x*-x(6)∞=3.880501726061247*10-75*10-6，满足要求。 2） 当ω=1, 迭代到第7次时， x(7)=(0.500000476837158,1.000000238418579, -0.499999940395355) x*-x(6)∞=7.748603820800781*10-75*10-6，满足要求。 3） 当ω=1.1, 迭代到第7次时， x(7)=(0.499999235607344,0.999999924905289,-0.500000016711780) x*-x(6)∞=8.561991469702157*10-75*10-6，满足要求。 第七章 方程求根 解：(1)∵ φx=1+1x2,φx=-2x-3,φx0=2*1.5-3≈0.592591 ∴ 该迭代公式收敛 (2)∵φx=31+x2,φx=23x1+x2-23，φx0=23*1.5*1+1.52-23 ≈0.455771 ∴ 该迭代公式收敛 (3)∵ φx=121+x，φx=-12x-1-32，φx0=12*1.5-1-32 ≈1.414211 ∴ 该迭代公式发散 由于(2)式与(1)式相比，L更小，故选择(2)式进行计算。 依题意可知，需求具有4位有效数字的近似根，取m=0,n=4故需满足 xk-x*≤xk-xk-112*100-4+1=12*10-3 计算结果如下表所示, k XK 1 1.481248 2 1.472706 3 1.468817 4 1.467048 5 1.466243 6 1.465877 由于x6-x512*10-3，故取x*=x6≈1.466 解：由题目，要求计算结果精确到4位有效数字，取m=0,n=4,即xk-x*≤12*100-4+1=12*10-3 (1)由牛顿法，有迭代公式xk+1=xk-fxkfxk=xk-xk3-3xk-13xk2-3 列表计算如下： k XK 0 2 1 1.888888889 2 1.879451567 由上表可知，当迭代2次时，xk-x*12*10-3,故取x*=x2≈1.879 (2)由弦截法，有迭代公式xk+1=xk-fxkfxk-fxk-1xk-xk-1 列表计算如下： k XK 0 2 1 1.9 2 1.881093936 3 1.879411060 由上表可知，当迭代2次时，xk-x*12*10-3,故取x*=x3≈1.879 (3)由抛物线法，有迭代公式 xk+1=xk-2fxkω+sgnωω2-4fxkfxk,xk-1,xk-2ω=fxk,xk-1+fxk,xk-1,xk-2(xk-xk-1) 列表计算如下: k XK 0 1 1 3 2 2 3 1.893149824 4 1.879135257 由上表可知，当迭代2次时，xk-x*12*10-3,故取x*=x4≈1.879 计算机习题.pdf 解:（1）Jacobi程序如下: function [x,count]=h_jacobi(n,deviation) %输入矩阵的维数n,误差||*||∞ %输出估计值x,迭代次数count h=zeros(n); %Hn for i=1:n for j=1:n h(i,j)=1/(i+j-1); end end x_acc=ones(n,1); %x