数值分析第章 解线性方程组的直接方法.pptVIP

§1 引言与预备知识 §2 高斯消去法 §3 高斯主元素消去法 §4 矩阵的三角分解法 克劳特分解方法 设A为n×n阶非奇异矩阵, 且各阶主子矩阵为非奇异,则矩 阵A的克劳特(Crout)分解为 A=LU 其中 这样,L、U中的元素都已求出。计算L的各列与U的各行的次序如图所示 。 对方程组Ax=b的系数矩阵A作出LU分解后,方程组便化为 LUx=b 则求解上列方程组就化为依次解方程组 Ly=b Ux=y 由于L为下三角矩阵,U为单位上三角矩阵,故上述方程组的求解极为方便。他们的计算公式分别为 用克劳特分解求解线性方程组Ax=b的计算过程为: ① LU分解过程:对于k=1,2,…,n依次计算 例4 用克劳特分解方法求解下列方程组 利用矩阵乘法可得到 这样原方程组就化为依次求下列两个三角形方程组 §5 向量和矩阵的范数 §6 误差分析 图