消元及矩阵分解(II).pptVIP

*/20 为 待定系数， 三、三对角方程组的追赶法 比较上式两边得 不用再计算，而L和U中的主要元素 和 的计算公式如下 (*) 定理2 如果三对角矩阵A的元素满足条件： 则由公式(*)。可计算出 ,且 求解三对角方程组Ax=f 等价于解两个三角形方程组 具体计算公式如下： （1）求解Ly=f，得向量y； （2）求解Ux=y，得方程组的解x。 在应用数学中，很多矩阵具有对称正定这一性质，因此利用对称正定矩阵的三角分解式求解对称正定方程组的一种有效方法，且分解过程无需选主元，有良好的数值稳定性。 四、对称矩阵的Cholesky分解 A对称：AT=A A正定：A的各阶顺序主子式均大于零。即 对称矩阵的Cholesky分解 定理1 设A为对称正定矩阵，则存在唯一分解A=LDLT，其中L为单位下三角阵，D=diag(d1,d2,…,dn)且di0(i=1,…,n) 证明： 定理2：设A为对称正定矩阵，则存在唯一分解 其中L为具有主对角元素为正数的下三角矩阵。 Cholesky分解法缺点及优点 优点：可以减少存储单元。 缺点：存在开方运算，可能会出现根号下负数。 四、改进Cholesky分解法 A=LDLT分解，既适合于解对称正定方程组