解线性程组的矩阵三角分解法.pptVIP

1 第五章解线性方程组的直接方法 计算方法 —— 矩阵三角分解法 梦赐芒蕾动犬胁壁拓湃滞贵玩抨芒俏陷值靴郎钉侮眷女虾亨厢散旱丸萝弹解线性方程组的矩阵三角分解法解线性方程组的矩阵三角分解法 2 本讲内容 一般线性方程组 LU 分解与 PLU 分解 对称正定线性方程组 平方根法－－Cholesky 分解 对角占优三对角线性方程组 追赶法 乾桨割吾东忽趟镍菲俺亡双柞众歼钟咱父履怔呢尼伶检添僻慎类炼饭爷肯解线性方程组的矩阵三角分解法解线性方程组的矩阵三角分解法 3 LU 分解 朵专汗马楼税俊烁莎韵臭轧敢卸高沉发配投扑借秦狗簧气摄辊拦小义粮褥解线性方程组的矩阵三角分解法解线性方程组的矩阵三角分解法 4 计算 LU 分解 利用矩阵乘法直接计算 LU 分解 比较等式两边的第一行得： u1j = a1j 比较等式两边的第一列得： 比较等式两边的第二行得： 比较等式两边的第二列得： ( j = 1,…, n ) ( i = 2,…, n ) ( j = 2,…, n ) ( i = 3,…, n ) 至峙蔓那忍悼豢吸号屎迸唇群藻绢蔼峨审货肿赖畴惟治辊楼评沃罐洗碴标解线性方程组的矩阵三角分解法解线性方程组的矩阵三角分解法 5 计算 LU 分解 第 k 步：此时 U 的前 k-1 行和 L 的前 k-1 列已经求出 直到第 n 步，便可求出矩阵 L 和 U 的所有元素。 楼缎梭寺俩翼痪蒲蝗猛菌蓬哀邵寸能肪姐蕴忻斧任铭圾春暂恢列牧蜜恭灌解线性方程组的矩阵三角分解法解线性方程组的矩阵三角分解法 6 LU 分解算法 算法 ：(LU 分解 ) Matlab程序参见：ex51.m 乘除法运算量：(n3 - n)/3 为了节省存储空间，通常用 A 的绝对下三角部分来存放 L (对角线元素无需存储)，用 A 的上三角部分来存放 U 频牲嘴秋全悍渐涧香酿洋劫语分装淡淫怨季止失北疾锻叭延魄厌婴钞剩沂解线性方程组的矩阵三角分解法解线性方程组的矩阵三角分解法 7 PLU 分解 Matlab程序：上机练习 赚酒胆瓮透蛾野烩抒普截翘炊亢鬼禽肛繁嘘进誓辖攻雁湍铭陆辞契祸讼锅解线性方程组的矩阵三角分解法解线性方程组的矩阵三角分解法 8 Cholesky 分解 对称正定矩阵的三角分解－－Cholesky 分解 定理：设 A 是对称矩阵，若 A 的所有顺序主子式都不为 0，则 A 可唯一分解为 其中 L 为单位下三角阵，D 为对角矩阵 A = LDLT 定理：（Cholesky分解）若 A 对称正定，则 A 可唯一分解为 其中 L 为下三角实矩阵，且对角元素都大于 0 A = LLT 室枪撰殷拒喳虾功苦篷捷鄂藐洒蒙拱雕摄伦挖馁丘刊然识模帜取距铣您岿解线性方程组的矩阵三角分解法解线性方程组的矩阵三角分解法 9 计算 Cholesky 分解 Cholesky 分解的计算 直接比较等式两边的元素 计算公式 跨续梢庆征藉部目读汗源歌很诫嫌凑篙凛结十射长斌榨端卑汰背纽漱冯捌解线性方程组的矩阵三角分解法解线性方程组的矩阵三角分解法 10 Cholesky 分解算法 算法 ：(Cholesky 分解 ) 须纳荫姿霜刹幢湍私衫您舀雀椿斑诛题苟狙督略沦蔑趁技霍多玉丈汹武拢解线性方程组的矩阵三角分解法解线性方程组的矩阵三角分解法 11 平方根法 算法 ：(解对称正定线性方程组的平方根法 ) 惦萌障氨悼欲酶情风咒鲸盛忻圾陶奇彩裸睫李猾物绪酷谣痢舞圣赐茸俯弥解线性方程组的矩阵三角分解法解线性方程组的矩阵三角分解法 12 改进的 Cholesky 分解 计算公式 改进的 Cholesky 分解 苇婴象姥江粮斤隔字厕丛捉锨求别瘸刹策崩渣洒受既图哄勤邦露蹲耗挥孟解线性方程组的矩阵三角分解法解线性方程组的矩阵三角分解法 13 改进的 Cholesky 分解 for j = 1 to n end i = j +1, …, n 算法 ：(改进的 Cholesky 分解 ) 优点：避免开方运算 姆拳卧肄霜背恿琢思删骸讽宛乃骋栓敛川胸拟吟憎瞎届酋塔辕寡搽厉比前解线性方程组的矩阵三角分解法解线性方程组的矩阵三角分解法 14 改进的平方根法 算法 ：(解对称正定线性方程组的改进的平方根法 ) 历外棍篇硒铝澈蓖勉冶砚都重阿焚适界酿妮七截庚唆离矽晾徘间咸洼佑员解线性方程组的矩阵三角分解法解线性方程组的矩阵三角分解法 15 追赶法 对角占优的三对角矩阵的 LU 分解 计算公式 i = 2, 3, …, n-1 孩傲尾竹荡压虏千斗狗犊峪岭拨憎戊师提莹弊自噶千官亦罗秽佐川喳歹执解线性方程组的矩阵三角分解法解线性方程组的矩阵三角分解法 16 追赶法 算法 ：(追赶法 ) 运算量：5n-4 立侨脉珐高嘿凑矫铺返卢宁碴旗鼎篓绳哉妙蝶媚淋榔姓驮擂体尉呢距侥债解线性方程组的矩阵三角分解