《矩阵与数值分析》上机大作业matlab.doc

《矩阵与数值分析》上机大作业 1.给定n阶方程组，其中 ， 则方程组有解。对和，分别用Gauss消去法和列主元消去法解方程组，并比较计算结果。 %产生三对角矩阵 n=84; %或n=10; A=zeros(n); b=zeros(1,n); for i=1:n-1 A(i,i)=6;A(i,i+1)=1;A(i+1,i)=8; end A(n,n)=6; for i=2:n-1 b(1)=6; b(i)=15; b(n)=14; end Ab=[A b]; %Gauss消元法 for j=1:n-1 %按列循环 for k=j+1:n %消元 Ab(k,:)=Ab(k,:)-Ab(j,:)*(Ab(k,j)/Ab(j,j)); end end x(n)=Ab(n,n+1)/Ab(n,n); for i=n-1:-1:1 %回代法求x for j=n:-1:i+1 Ab(i,n+1)=Ab(i,n+1)-Ab(i,j)*x(j); end x(i)=Ab(i,n+1)/Ab(i,i); endx=1.000000000000000 x=1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000001 1.000000000000000 0.999999999999998 1.000000000000000 1.000000000000004 1.000000000000000 0.999999999999993 1.000000000000000 1.000000000000012 1.000000000000000 0.999999999999979 1.000000000000000 1.000000000000028 1.000000000000000 当n=84时，Gauss消去法的解是错解 Columns 34 through 39 2147483649.00000 -4294967295.00000 8589934592.99999 9999 34359738368.9998 Gauss列主元消去法x与x=(1,1…1)T偏差不大 Columns 34 through 39 1.000000172108412 0.999999661246936 1.000000655651093 0.999998776117961 1.000002098083496 综上，高斯列主元消去法可以避免小数作除数带来的误差，获得满意的数值解。 2. 设有方程组，其中， 选取不同的初始向量和不同的右端向量，给定迭代误差要求，用Jacobi和Gauss-Seidel迭代法计算，观测得出的迭代向量序列是否收敛。若收敛，记录迭代次数，分析计算结果并得出你的结论。 选取b={1}; x0={0};则Jacobi法 （1）Jacobi迭代法 n=20; A=zeros(n);b=1*ones(n,1); m=3; for i=1:n-2 A(i,i)=m;A(i,i+1)=-0.5;A(i,i+2)=-0.25;A(i+1,i)=-0.5;A(i+2,i)=-0.25; end A(19,19)=m; A(20,20)=m; A(19,20)=-0.5;A(20,19)=-0.5; %e=ones(n,1);b=A*e;%习题2的初始条件 K=1000; U=-triu(A,1); L=-tril(A,-1); D=diag(diag(A).^(-1)); f=D*b; B=D*(L+U); if max(abs(eig(B)))=1 disp(迭代不收敛) return end xp=zeros(n,1); for k=1:K x=B*xp+f; if max(abs(xp-x))0.000001 break; end xp=x; end x= 0.481636071409813, 0.573414961849327, 0.632801573980968, 0.652101940832594 0.660