方程组直接3角分解法.pptVIP

4.2 直接三角分解?法;4.2 直接三角分解?法;可得U的第k行元素;用向后回代的方法即可求得?x。设x=（x1 ，x2， ·?·· xn) T, y=（y1?, y2, ··· yn) T?,b= （b1 ，b2， ··?· bn) T, 则有计算公式? ; ? ;2.列选主元的三角分解法 ;当i=k时， si对应于（4.2.3）中的ukk，它可能不宜在（4.2.4）作除法。当i=k+1,k=2,….n， si对应于（4.2.4）中的分子。记 ;例4.6 用列选主元的三角分解法解;由于方车程组的右端参与了消元计算，所以Ly=Pb的解为y=b（3）=（20，14/3，216/39） T 。解Ux=y得x=（1，2，3） T; ? ;?;.另一方面,有;例4.7 用追赶发求解三对角方程组Ax=d，其中 解 由(4.2.9)得; 对另一类方程组，在周期样条插值等为题遇到的循环三对角方程组Ax=d，其中;　　当?A为对称正定矩阵时，对A可直接?作LU分解。由（4.1.8）式?可得下面的定??。 ? 定理4.7 ?设A∈Rn×n， A = AT?且A的顺序主子式Di≠0（I=?1,2,…,n），则存在唯一的?单位下三角阵和三角阵，使 ? ? ? A?= LD LT ? 定理4.8 ? 设A∈Rn×n，当A 为对称?正定矩阵 ，则存在唯一的对角元?素为正的下三角阵L，使 ? ? ? A=? L LT ? ; 令 　，则有 由分解式 的唯一性可得（4.2.3）分解式的唯一性。定理得证。 称（4.2.13）式为矩阵A的Cholesky分解。利用A的Cholesky分解式来求解方程组Ax=b的方法称为Cholesky方法或平方根法，这是因为计算过程含开方运算。 ;; ? 解 不难验证系数矩阵?是对称正定的，按（4.2.14?）和（4.2.15）依次计算得?;则可避免开方根运算，称为改进的平方根法。 它即适合于求接对称正定方程组，也适合于A求解对称且其顺序主子式全不为零的方程组。分解式的计算公式为(j=1,2,…n);解Ly=b得y=(?6,1,-1)T。解LTx=D?-1y得x=(2,1,-1)T?。