数值计算CH2 解线性方程组及直接法—2.ppt

华长生制作 华长生制作 §4 直接三角分解法 顺序主元 2、Crout 分解 3、求解 4、紧凑格式（分解与求解） 5、部分选主元的Doolilttle分解 §5 平方根法 §6 追赶法 * * 第二章 解线性方程组的直接法 张红梅 自动化学院 2010年3月 2.4 直接三角分解法 Gauss 消去法的一种变形。利用初等变换将方阵化为两个三角阵的乘积，进而将原方程组化为两个三角形方程组，即 方法的产生及思想 基本的三角分解法 Doolittle 分解 Crout 分解 部分选主元的Doolittle分解 平方根法 追赶法 …… 方法的分类 ＝ + = = = 对于确定的 A，如果首先对它进行上述三角 分解，那么对于不同的 ，只需计算两个更 简单的方程组 和 。 高斯消元法的矩阵形式 Step 1: 记 L1 = ，则 Step n ? 1: 其中 Lk = 记为 L 单位下三角阵 A 的 LU 分解 记 U = 上三角 ？ A 能进行这样分解的条件是... 第 k 步消元时，在做行交换之前位于主对角线上的元素 称为顺序主元。显然，不带行交换的消元过程完成的条件是： 这也是 A 可以进行直接 LU 分解的条件。 由矩阵分块乘法可知： Ak = LkUk Dk = det Ak = detUk det Lk = 1 1 1 = 1 顺序主子阵 顺序主子式 由归纳法得： 的充要条件是A 的顺序主子式D1 ,D2,…, Dk 均不为零, 对于k =1,2,…,n-1 成立。 定理3.1 若 n阶方阵A 的顺序主子式 均不为0，则 A 的 LU 分解存在且唯一（L 为单位下三角阵, U为上三角阵）。 证明：存在性由上述分析已得，以下采用证明唯一性。 采用反证法，若不唯一，则可设： A = L1U1 = L2U2 推出 上三角 单位下三角 = I 因此U1 =U2 且 L1 = L2。 消元结束时,U 可能是奇异的，即 。这种情况对分解没有影响，但是对回代能否进行起决定作用。 1、Doolittle 分解 通过比较法直接导出L 和 U 的计算公式。 思路 Step 1a: 比较等式两侧第1行元素，得U 的第1行 Step 1b: 比较等式两侧第1列元素(从第2个开始)，得L 的第1列 Step 2a: 比较等式两侧第2行元素(从第2个开始)，得U 的第2行 Step 2b: 比较等式两侧第2列元素(从第3个开始)，得L 的第2列 ··· ··· 直接比较两侧对应元素 固定 i : 对 j = i, i+1, …, n 有 lii = 1 a 将 i ，j 对换，对 j = i, i+1, …, n 有 b Doolittle 分解 Step 1: u1j = a1j; lj1 = aj1 / u11; ( j = 1, …, n ) Step 2: 计算 和 for i = 2, …, n?1; Step 3: a b 一般采用列主元 法增强稳定性。但注意 也必须做相应的 行交换。 类似于推导Doolittle分解的方法，可得Crout分解公式如下： 可证，当A满足定理3.1的条件时，A可分解为一个一般的下三角阵 L 和一个单位上三角阵U 的乘积，这种分解称为Crout 分解。 1. Crout分解的顺序与Doolittle分解的顺序正好相反：先计算 L 后计算 U； 2. 只要考虑 A 的转置即 AT 的 LU 分解 ，则 即是 A 的 Crout 分解。 对系数矩阵实现Doolittle 分解或Crout分解后，原方程组化为两个三角方程组，它们的求解公式分别为： Doolittle 求解 Crout 求解 例1 p.43 电算时, 采用节省空间的方式存储数据, 将矩阵的分解与方程组的求解同时进行。 思想 数据存储： A(1:n; 1:n+1)←(A,b) 分解计算时，求得的L和U的元素覆盖掉A的相应元素，而y覆盖掉b的相应元素，例如，用Doolittle法： A(i, j)←uij (i=1,2,…,n; j=i,i+1,…,n+1) A(i, j)←lij (j=1,2,…,n-1; i=j+