高等工程数学张韵华版第六章课后题

答案仅供参考！

• 答案勘误：修改了第四题和第六题的答案。

目录

第6章 线性方程组数值解

        6.1 高斯列主元消元

                6.1.1 高斯消元法

                6.1.2 列主元消元法

        6.2 直接分解法

                6.2.1 LU分解

                6.2.2 对称正定矩阵的LDLT分解

        6.3 解线性方程组的迭代法

                6.3.1 雅可比迭代

                6.3.2 Gauss-Seidel迭代

                6.3.3 松弛迭代

附录6 用Mathematica求解方程组和矩阵分解

主要内容：

1、Gauss列主元消元法

2、矩阵的Courant分解和Doolittle分解

3、Jacobi迭代格式及其收敛性

4、Gauss-Seidel的迭代格式及其收敛性

课后题：

1. 给岀解下三角方程组的算法:

在解方程组中用到以下三种初等变换:
(1)对换某两个方程的次序;
(2)对其中某个方程的两边同乘一个不为零的数;

(3)把某一个方程两边同乘一个常数后加到另一个方程的两边。

2. 用 Doolittle 分解求解方程组：

3.用柯朗分解求解下面的方程组:

4. 用雅可比方法解方程组：

(1) 请写岀迭代格式；

(2) 请写出迭代矩阵；

(3) 讨论其迭代收敛性；

(4) 以为初值，计算

5.用 Gauss -Seidel 方法解方程组

(1) 请写岀迭代格式；

(2) 请写出迭代矩阵；

(3) 讨论其迭代收敛性；

(4) 以为初值，计算

6.设有系数矩阵

证明:

⑴ 对系数矩阵 A，雅可比迭代收敛，而 Gauss -Seidel 迭代不收敛;

⑵ 对系数矩阵 B，雅可比迭代不收敛，而 Gauss -Seidel 迭代收敛。

注：

雅可比迭代收敛条件：对于方程组AX=y,构造雅可比迭代式X^(k+1)=BX^(k)+g.当迭代矩阵B的谱半径p<1时,迭代收敛,这是收敛的充要条件,迭代矩阵的某范数||B||<1,是迭代收敛的充分条件。当方程组的系数矩阵A具有某些特殊性质时（比如A为严格行对角优矩阵，或者严格列对角优矩阵）,可直接判定由它生成的雅可比迭代矩阵是收敛的。

Gauss-Seidel迭代收敛条件：判断Gauss-Seidel迭代收敛与判断雅可比迭代收敛的方法类似。一方面从Gauss-Seidel迭代矩阵S获取信息,当p(S)<1或S的某种范数‖S‖<1时,迭代收敛;另一方面,有时可直接根据方程组系数矩阵的特点做出判断。若方程组系数矩阵为列或行严格对角优的,则Gauss-Seidel迭代收敛。若方程组系数矩阵为对称正定矩阵,则Gauss-Seidel迭代收敛。