答案仅供参考!
答案勘误:
- 修改了第四题和第六题的答案。
- 核对了第二版教材后,修改了第三题第二问、第五题第二问答案
目录
第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迭代收敛。