05-解线性方程组的直接法.ppt

第五章 §5.1 高斯消去法 §5.2 解三对角方程组的追赶法 §5.3 矩阵的三角分解及其在解方程组中的应用 §5.4 平方根法 §5.5 向量和矩阵的范数 §5.6 线性方程组的性态和解的误差分析 本章小结 本章小结 本章作业 习 题 习 题 习 题 习 题 习 题 习 题 习 题 习 题 习 题 习 题 习 题 习 题 习 题 习 题 习 题 习 题 习 题 习 题 习 题 海南大学 土木建筑工程学院 2008年8月 解线性方程组的直接法 本章目录 1. 消去法 (1) Jordan消去法 (2) Gauss消去法 (3) 选主元的Gauss消去法 2. 矩阵三角分解法 (1) 直接三角分解法 (2) 平方根（Cholesky分解）法 (3) 解三对角方程组的追赶法 3. 向量、矩阵的范数及方程组的性态 (1) 向量的范数 (2) 矩阵的范数 (3) 矩阵的谱半经 (4) 矩阵的条件数 P230 第 8 题 P231 第 9 题、第 10 题 本章结束 1. 用列主消元法解方程组 【解】对（A , b）做选主元及消元过程： 做回代过程： 2. 分别用 Doolittle 和 Crout 分解法求解方程组 【解】方阵系数矩阵为 常数项为 （1）Doolittle 分解： 按 U 的第一行、L 的第一列 , U 的第二行、L 的第二列 , U 的第三行次序 , 对比等号两边矩阵元素 , 得矩阵 L , U元素 求解 LY=b , 得 (y1 , y2 , y3)=(4.0 , 4.0 , 0.6) （2）Crout 分解： 求解 UX=Y , 得 (x1 , x2 , x3)=(1.0 , 1.0 , 1.0) 按 L 的第一列、U 的第一行 , L 的第二列、U 的第二行 , L 的第三列次序 , 对比等号两边矩阵元素 , 得矩阵 L , U元素 求解 LY=b , 得 (y1 , y2 , y3)=(2.0 , 1.6 , 1.0) （2）Crout 分解： 求解 UX=Y , 得 (x1 , x2 , x3)=(1.0 , 1.0 , 1.0) 3. 用改进的平方根法解方程组 【解】对 i=1 , d1=a11=4 对 i=2 , t21=a21=1 , l21=t21/d1=1/4=0.25 , d2=a22-t21l21=2.75 对 i=3 , 对 i=4 , 由于求解 AX=b 等价于求解 （1）LY=b 求 Y （2）DLTX=Y 求 X 求 LY=b 得 求解 LTX=D-1Y 求得 4. 用追赶法解三对角矩阵方程组 【解】因 A 是严格对角占优矩阵 , 因此可直接对 A 进行 Crout 分解 计算公式为： 再次令 n=5 , aii=4 , aii-1=ai-1i=－1 代入得 求解 LY=b , b=( 100 , 200 , 200 , 200 , 100 )T 求解 UX=Y , y=( 25 , 60 , 69.6429 , 72.2488 , 46.1538 )T 求得方程组向量为 x = ( 46.1538 , 84.6154 , 92.3077 , 84.6154 , 46.1538 )T 5. 用 LDLT 分解法解方程组 【解】（1）分解 由矩阵乘法解得 d1=3 , l21=1 , l31=5/3 d2=2 , l32=2 , d3=2/3 （2）解方程组 解得 , y （3）解方程组 6. 计算方阵 的三种常用范数 【解】（1） 由 , 计算 特征多项式 解得 ATA 的特征值是 1 , 8 , 32 所以 7. 用直接分解法计算 的逆矩阵。 【解】 按 U 的第一行、L 的第一列 , U 的第二行、L 的第二列 , U 的第三行次序 , 对比等号两边矩阵元素 , 得矩阵 L , U元素