Matlab在对称正定矩阵的改进平方根分解法中的应用.pdfVIP

第 23 卷第 2 期 　　 　 　　　　　　怀化学院学报 　　 　 　　　　　　Vol23No2 2 0 0 4 年 4 月 　 　 　 　 　J OURNAL OF HUAIHUA UNIVER SITY 　 　 　 　 　 　Apr , 2 0 0 4 数 学 研 究 Matlab 在对称正定矩阵的改进平方根分解法中的应用 刘亮元 (怀化学院 物理系 , 湖南 怀化 418008) 摘 　要 : 用改进的平方根法分解大型的正定对称矩阵是一种行之有效的方法 , 用 Matlab 语言编写程序求解 , 以 此对这种方法进行进一步探讨 关键词 : Matlab ; 　对称正定矩阵 ; 　分解 中图分类号 : O15121 　　　文献标识码 : A 　　　文章编号 : 1671 - 9743 (2004) 02 - 0022 - 03 用矩阵理论分析高斯消去法, 建立矩阵的三角分解定理, 实质是将系数矩阵A 分解为一个单位下三角矩阵 L 和一个单位上三角矩阵U 的乘积, 即矩阵的L U 分解 在实际问题中常遇到对称正定矩阵方程组的求解问题, 对这种具有特殊性质的系数矩阵, 采用改进平方根法求解是一种行之有效的方法, 它是在高斯消去法基础上建 立起来的, 算法思想较复杂些, 但计算量大约为高斯消去法的一半, 用Matlab 编写程序求解使这类问题变得更 简单 正定矩阵 A 有以下五个性质 : ) - - 1 矩阵 A 非奇异, 即 A 存在, 且 A 仍是正定矩阵 ) T 2 设 P 是任意的 n 阶非奇异矩阵, 则矩阵 B = P A P 仍是正定矩阵 ) 3 A 的特征值为正数 ) 4 A 的各阶主子矩阵仍为正定矩阵 ) 5 A 是正定矩阵的充要条件为A 的各阶主子式均为正 1 　改进平方根法的基本思想 设 A 为对称正定矩阵, 则 A 有三角分解, 且具有唯一性 : 1 u11 u12 Λ u1 n l21 1 u22 Λ u2 n A = L U = M M O O M ln1 ln2 Λ 1 unn 1 u11 1 l21 Λ ln1 l21 1 u22 1 Λ ln2 = M M O O O M ln1 ln2 Λ 1 unn 1