数值分析第四次作业-国科大_国科大数值分析肖良资源.docx

。Ax=b,

。

Ax=b,其中

事

数值分析习题

第四单元

1.设方程组：

(1)考察用Jacobi迭代法，Gauss-Seidel迭代法解此方程组的收敛性；

(2)取x(0)=[0,0,0,用Jacobi迭代法及Gauss-Seidel迭代法解此

方程组，要求当|~-x*L102时终止迭代。

(3)用SOR方法解上述方程组(取o=0.9,初始解为[0,0,0]),要

求当|F-x102

解：令方程组为

终止迭代。

;

;

(1)因矩阵A为严格对角占优矩阵，故Jacobi迭代法，Gauss-Seidel

迭代法均收敛。

(2)Jacobi迭代法的计算公式为

经过11次Jacobi迭代算法收敛，近似解为[4.000242,3.003129,1.999862].

Gauss-Seidel迭代法的计算公式为

经过6次Gauss-Seide迭代算法收敛近似解为[-3.999314,3.000003,1.999864]。

(3)SOR迭代法的计算公式为

经过5次SOR迭代算法收敛，近似解为[-3.996732,3.000028,1.999362]。

实验Matalb代码及迭代过程详见附录表格和pr4.m文件。

2.考虑线性方程组Ax=b其中

(1)a为何值时，A是正定的?

(2)a为何值时，Jacobi迭代法收敛?

(3)a为何值时，Gauss-Seidel迭代法收敛?

解：(1)若A各阶顺序主子式值都大于0,即可判定为正定矩阵。

A的一阶顺序主子式川0满足条件；

A的二阶顺序主子式满足条件；

A的三阶顺序主子式

得到-1a1。

所以当-1α1,A是正定的。

(2)显然，矩阵A为对称实矩阵，且对角线元素q0,(i=1,2…,n),

要判断线性方程组Ax=b是否能Jacobi迭代法收敛，即只需判定A和

2D-A都正定。其中D为A对角元素的对角矩阵。

由(1)可知，A正定时，-1a1。因此

2D-A的一阶顺序主子式川0满足条件；

2D-A的二阶顺序主子式

2D-A的三阶顺序主子式

得到-1a1。

满足条件；

所以当-1α1,2D-A也是正定的。

综上，当-1α1时，线性方程组Ax=bJacobi迭代法收敛。

(3)Gauss-Seidel迭代法的迭代矩阵为

谱半径p(B)=mx{a2,α}1即-1a1。

所以当-1α1,线性方程组Ax=bGauss-Seidel迭代法收敛。

附录

表1用acobi方法得到的迭代解序列

迭代步

X?

X?

X?

1

-2,400000

5.000000

0.300000

2

-4.460000

4.250000

2.280000

3

-4.556000

2.745000

2.467000

4

-3.991400

2.627500

2.034700

5

-3.857940

2.984800

1.886530

6

-3.971226

3.092250

1.967028

7

-4.030306

3.023680

2.021920

8

-4.013856

2.981464

2.013165

9

-3.995218

2.989954

1.997210

10

-3.995423

3.002590

1.996030

11

-4.000242

3.003129

1.999862

表2用Guass-Seidel方法得到的迭代解序列

迭代步

X?

X?

X?

1

-2,400000

4.400000

2.100000

2

-4,580000

2.805000

2.057500

3

-3.933500

2.987875

1.983063

4

-3.991763

3.010528

2.001511

5

-4.004513

2.998116

2.000338

6

-3.999314

3.000003

1.999864

表3用SOR方法得到的迭代解序列

迭代步

X?

X?

x?

1

-2.160000

4.014000

1.742580

2

-4.134704

3.186931

2.048976

3

-4.089581

2.976498

2.014677

4

-4.003139

2.990339

1.999424

5

-3.996732

3.000028

1.999362

Matlab实验代码：

pr4.m

%%第三题

clearall;cle;

formatshort

A=[5