第五章 解线性数方程组.docVIP

求解线性方程组的直接解法 5.3??????特殊矩阵的三角分解 ①? 实对称矩阵的LDLT分解 设A是实对称阵,且A的所有顺序主子式均不为零，则LDR分解中R=LT, 故可用以作LDLT分解.这就是说,当A的对角元素非零时，我们可以作LU分解,也就得到LDLT分解,L相同,是单位上三角阵,U的对角元素构成D.不过没有利用对称性,存储量运算量都未能节省—预计是一半。试用n=3的计算表格说明如何实现节省。 d1=u11 =a11 u12=a12 l21=u12/d1 u13=a13 l31=u13/d1 ? d2=u22=a22-l21u12 u23=a23-l21u13 l32=u23/d2 ? ? u33=a33-l31u13-l32u23 这样,可用上半部元素逐列计算D,LT。也可用下半部元素逐行计算L,D。引进輔助量t1, t2代替u1j,u2j,并利用对称性得到： d1=a11 ? ? t1=a21 l21= t1/d1 d2= a22-t1l21 ? t1=a31 l31=t1/d1 t2=a32-t1l21 l32=t2/d2 d3=a33-t1l31-t2l32 ? 据此不难写出LDLT分解A=LDLT的计算公式和程序(逐行计算L,D). d1=a11 for i=2:n for j=1:i-1 tj=aij-lj1t1-lj2t2-…-l