3.2典型例题精解.docVIP

3.2典型例题精解 例3.2.1 设 , 解三角矩阵方程组 和 解 由下三角方程组 向前回代方法有 得，所以 再由上三角方程组 向后回代有 得，所以。 例3.2.2 设 ， 试用Gauss顺序消元法解方程组Ax=b。 解 线性方程组的增广矩阵为 则，再由得 再由，计算得 由回代过程，计算得 所以原方程组的解。 例3.2.3 用Gauss列主元消元法解线性方程组 解 当k=1时，由于，因此，第二行为主行，即为第1步消元的主元，交换原方程组的增广矩阵第一行与第二行，即 计算乘法因子 消元得 当k=2时，由于，因而第三行为主行，交换第3行与第2行得 计算乘法因子 消元得 最后回代解得，故解得原方程组的解。 例3.2.4 用Gauss全主元消元法解方程组 解 设数据组，用于记录列交换。 当k=1时，由于，所以取主元为。交换，即。 交换原方程组的增广矩阵第1行与第4行，在交换第1列与第4列得 计算乘法因子 消元得 当k=2时，选主元，交换，即，交换增广矩阵第2行与第3行，在交换第2列与第3列得 计算乘法因子 ， 消元得 当k=3时，选主元，计算乘法因子 消元得 最后回代得 交换得 在交换得 故解得方程组的解为。 例3.2.5 用Doolittle 分解法解线性方程组 解 首先对系数矩阵A做分解A=LU 于是得到 ， 解方程组的解得 解方程组得。 例3.2.6 求矩阵 的Cholesky分解A=LL。 解 由式（3.1.6）,当k=1时 因此 当k=2时 从而 当k=3时 即 故得到Cholesky分解A=LL,其中 例 3.2.7 用改进的Cholesky分解法解方程组 解 由改进的Cholesky分解（3.1.7）式，对i=1,2,3,4计算得 i=1时， i=2时， i=3时， i=4时， 由Ly=b(b=(3,3,3,3))解得 再由解得 故解得原方程组的解。 例3．2．8 用追赶法解三对角方程组 解 首先，对方程组系数矩阵A作LU分解，由追赶法（3.1.9）,对i=1,2,3,4计算得 当i=1时, 当i=2时, 当i=3时, 当i=4时, 所以 由解得 再由解得，故得原方程组的解为。 例3.2.9 设矩阵 试用Gauss-Jordan原地求逆法求解A。 解 设，用于记录选列主元所在的行。由Gauss-Jordan原地求逆法： 当k=1时，计算得 当k=2时，计算得 当k=3时，计算得 由知，交换第2列与第3列后，再交换第1列与第3列得 例3.2.10 设，求向量的范数。 解 由向量范数定义计算得 例3.2.11 设为对称正定矩阵，，证明是中的一种向量范数。 证明 只要验证它满足向量范数的三个条件。 （1）因为对称正定，因此对任意，，有 当且仅当时； （2）对任意有 （3）因为对称正定，所以总可以分解为，其中为非奇异的下三角阵，于是 从而对任意，，恒有 综上所述，是中的一种范数。 例3.2.12 设矩阵 试求，，。 解 由范数定义计算得 例3.2.13 设矩阵 试求,,,。 解 因为，所以。由 即 解得，所以，故。 对于,由 即，解得，所以，故 。 而由 即，解得，所以 例3.2.14 设矩阵 ， 且矩阵的特征值，求。 解 因为，所以,因此有 例3.2.15 设线性方程组 若系数矩阵有微小扰动 试求解得相对误差及相对误差估计。 解 原方程组的解为。把此方程组的系数矩阵记作A，则，而 因此 这个方程组的条件不太好。 现在系数矩阵A做微小扰动，则方程组 的解为。解的相对误差为 由于 因此，，由定理3.7得到解的相对误差估计为 例3.2.16 Hilbert矩阵 讨论的性质。 解Hilbert矩阵是有名的坏条件矩阵，当时 ， 于是 当时，可计算得。一般，当n越大时，的病态越严重。 例3.2.17 设 讨论顺序消去法求解方程组的数值稳定性。 解 方程组的精解解为 经计算得 因此这个方程组的条件并不坏。但是，若应用基本消去法来求解它，把 元素消为零，得到方程组 这时，若以9位十进制数运算，便得到 由回代得解为。 由此可见，该方程组的性态不是很坏，但用某一算法求解问题，也可能产生较大的误差，因此该方程组用顺序消元法数值不稳定。 例3.2.18 设线性方程组 试用迭代法及迭代法求方程组的解 解 该方程组的精确解为。迭代公式为 取初始向量，迭代结果如表3.2.1所示。 表3.2.1 k 0 1 2 3 4 5 6 x 0 0.9 0.97 0.991 0.9973 0.99919 0.999757 x 0 0.7 0.91 0.973 0.9919 0.99757 0.999271 x 0 0.6 0.88 0.964 0.9892 0.99