还剩4页未读,继续阅读
文本内容:
矩阵三角分解开题报告矩阵三角分解开题报告范文 在近代数学、工程技术、经济理论管理科学中,大量涉及矩阵理论的知识,很多问题都可以归结为矩阵并最终通过矩阵来解决经查阅发现,目前关于矩阵三角分解的应用研究不少,但对三角分解缺乏系统的研究 矩阵三角分解法是指高斯消去法解线性方程组的变形解法其实质就是将系数矩阵A分解为两个三角形矩阵L和U相乘,即A=LU
一、矩阵的直接三角分解 矩阵的直角三角分解即可以不经过消元步骤,直接将矩阵进行分解 定义1设A∈Rn×n,若A能分解为一个下三角矩阵L与一个上三角矩阵U的乘积,即A=LU,则称这种分解为矩阵A的三角分解
(1)如果A可分解为A=LDU,其中L是单位下三角矩阵,D是对角矩阵,U是单位上三角矩阵,则称A可作LDU分解;
(2)如果在A=LU中,L是单位下三角矩阵,U为上三角矩阵,则称此三角分解为杜利特(Doolittle)分解;
(3)如果在A=LU中,L是下三角矩阵,U是单位上三角矩阵,则称此三角矩阵为克劳特(Crout)分解 定理1n阶方阵A非奇异的充要条件为(或A经行、列变换后)存在LDU分解其中L为n阶单位下三角矩阵,D为n阶非奇异对角阵,U为n阶单位上三角矩阵 推论1奇异矩阵不能进行LDU分解 推论2若矩阵A有奇异主子矩阵,则A不能直接进行LDU分解 第2章线性代数方程组数值解法I直接法
1.矩阵 事实上,顺序Gauss消去过程对应一个矩阵的三角分解,即对Axb的顺序Gauss消去过程的结果,把矩阵A分解成两个三角矩阵L与U的乘积ALU下面来证实这一点.依次取第k步消元的乘法 kk likaikik1k2n/akk k1kk则直接验证可知第k步消元aij的结果等价于对Ak左乘Lk:aijlikakj Ak1LkAk 于是经过n1步消元应有 u11u12u13 u22u23Ln1L2L1AUU
2.
3.1u33 这里U为上三角矩阵,另外,又容易直接验证Lk有下列两个基本性质 1Lk的逆阵存在,且有 1 11 1lLk1kk
2.
3.2 1lnk 1 2逆阵Lk的乘积 1 1l21111 L1L2Ln1==L(单位下三角矩阵)(
2.
3.3) 1ln1ln1 111 从而对
2.
3.1式两端依次左乘Ln1L2Lk可得111 U=LUAL1L2Ln1 L就是
2.
3.3式所示的单位下三角矩阵这就是矩阵的三角分解或称LU 分解 ALU称为A的doolittle分解 ALULDU=LU称为A的克劳特分解 ALDU称为A的LDU分解 对于于有选主元和换行步骤的Gauss消去过程,也可证明它对应于“A左乘排列矩阵P的LU分解”,即有PA=LU例
2.
3.1用直接三角分解法解方程组(
2.1节中的实例) 232x10 1x12224317x3 解把解法分为3个步骤
①令A=LU,用Doolittle分解,即令 u11u12u132321 122lluu212223 41u3331l31l32 考虑A的第1行,对比右边两矩阵的乘积,有 21u11u112 31u12u12321uu2 1313 此结果即U的.第1行与A的第1行全同,这对一般情形也是适用的,因此,在分解计算中,此结果也可直接写出接着,再依次考虑A的第1列、第2行、第3列除去已考虑过的元素,作同样比较有 l211/21l21u11 3lul3/2311131 2l21u121u22u221/2 2l21u131u23l233 1l31u12l32u22l327 4l31u13l32u231u33u3328 1232 1/211/23即得A 1283/27
②用前推过程解下三角方程组 1y10y10 1/2y1得y1122 13/27714y3y3
③用回代过程解上三角方程组 x1232x102 x1得x11/2322 28141/2x3x3 下面以不包括选主元和换行的Doolittle分解为例,给出解n阶方程组Axb的一般计算公式及整个求解过程(分3个步骤)
①令ALU,即令 a1n1u1na11a12u11u12 aal1auu2121222n222n ll1aaaun1n1n2nnnnn1 利用矩阵乘法规则并对比等式两边对应元素由A的第1行得 a1j1u1jj1,2,,na1ju1jj1,2,,n
2.
3.5 由A的第1列(除第1行元素外)得 ak1lk1u11k2,3,,nlk1ak1/u11k2,3,,n
2.
3.6 依此类推,由A的第k列(1kn)(除前k1列元素外)得 akjlkrurjukj r1 k1 ukjakjlkrurjjk,,1,,n r1 k1
2.
3.7 由A的第k列(1kn)(除前k列元素外)得 aiklirurklikukk r1k1 likaiklirurk/ukkik1,,n r1 k1
2.
3.8
②求解下三角方程组Lyb得 y1b1 i1 3,,nyibiliryri2, r1
③求解上三角方程组Uxy得 xnyn/unn n ,2,1xiyiuirxr/uiiin1, ri1 这就是用直接三角分解法求解方程组的公式,其中第
①步中的前两个公式也可合并入后两个公式;第
②,
③步中的前一公式也可并入后一公式,这时当公式中出现和 r10 rn1 n 时均不执行计算,作零处理 (在列主元Gauss消去法一节中已提过这个附注) n3 可以推出,Doolittle算法的乘除法次数大致为,与Gauss消3 去法大致相同,故就计算量而言,采用Doolittle算法解方程组并无特别优势(因为我们已拥有相当高效的列主元Gauss消去法)应用中,主要借助直接三角分解法的处理方法来处理具有特殊情况的方程组这就是下一节要介绍的解三角方程组的追赶法和解对称正定方程组的平方根法模板内容仅供参考 。