matlab 矩阵和数值分析,《矩阵与数值分析》上机大作业matlab

Gauss 消去法 Gauss 列主元消去法 x=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

(2) 当n=84时，Gauss 消去法的解是错解

Columns 34 through 39

2147483649.00000 -4294967295.00000 8589934592.99999 -17179869182.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. 设有方程组Ax b =，其中2020A R ?∈，

3

0.50.250.5

30.50.250.50.50.250.530.50.250.53A --??

?-- ? ?-- ?

= ? ?

-- ?

-- ? ?--?

?

(a) 选取不同的初始向量(0)

x

和不同的右端向量b ，给定迭代误差要求，用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;