于是可以对PAPT进行块对角化，即 秦 -g非形e 其中C-EB-1ET是Schur补，对称非奇异 不断重复以上过程，就可以得到A的块LDLT分解： PAPT=LDLT 其中D是拟对角矩阵，即块对角矩阵且对角块的大小为1×1或2×2. http://math.ecnu.edu.cn/~jypan 15/42 于是可以对 P AP⊺ 进行块对角化, 即 P AP⊺ =   I 0 EB−1 I     B 0 0 C − EB−1E ⊺     I B−1E ⊺ 0 I   , 其中 C − EB−1E ⊺ 是 Schur 补, 对称非奇异. 不断重复以上过程, 就可以得到 A 的块 LDLT 分解: P AP⊺ = LDL˜ ⊺ , 其中 D˜ 是 拟对角矩阵 , 即块对角矩阵且对角块的大小为 1 × 1 或 2 × 2. http://math.ecnu.edu.cn/~jypan 15/42