参考53直接三角分解方法.pptxVIP

5.3 直接三角分解法;5.3 直接三角分解法;1.不选主元的三角分解法;如果U的第1至k-1列和L的第1至k-1列已经算出，则由;可得L的第k列元素;（5.3.5）; ;解下三角方程组 Ly = b，即 解(单位)上三角方程组 Ux= y，即;二 Crout分解（＊＊）; 实现 A=LU，则 Ax=b 可以化为 LUx=b。 令 Ux=y，则 Ly=b。 由 Ly=b解出 yi: 再由 Ux=y 解出 xi: ;例5.5’用Crout分解求解方程组：;解下三角方程组 Ly = b，即 解(单位)上三角方程组 Ux= y，即;在数值求解常微分方程边值问题、热传导方程和建立三次样条函数时，都会要解三对角方程组：AX=b;(5.3.7);计算次序为 ;称为(5.3.9)、(5.3.10)和(5.3.11)为求解三对角形方程组的追赶法,又称为Thomas算法。追赶法求解 Ax=f 仅需5n-4次乘除运算。工作量较小。 　　追赶法能实现的条件是ui≠0，i=1,2,…,n.下面给出追赶法的一个充分条件。;(2)***如果A满足Gauss消去法的条件,也可用LU (Crout)分解法求解. 并且, L和U有如下形式; ;;例 用追赶法求解三对角线性方程组AX=b，其中;作Crout分解;练习 用追赶