第四节 阵分解.pptVIP

第四节 矩阵分解 一 矩阵的LU分解 上页 返回 第三章 第四节 下页消元法的实质是对增广矩阵作初等行交换,而初等变换是可以用矩阵运算来描述的.消元过程相当于下述矩阵乘法运算。 因此，由分块矩阵乘法可得 膜戳珊模皱台鲸了劳拇晌邦笨佃统斟降敷瞳叼守茬煮薄唐侣碑饯被冀电昏第四节 矩阵分解第四节 矩阵分解 令 可得 可见只要消元过程能进行到底，就有以下等价关系 淆彤诱途竭钮拦野滑甭妥舜盈咽酗眺冀铸吊弹石蹦然砷喇阿眺倍灰哄袒奈第四节 矩阵分解第四节 矩阵分解 消元过程相当于分解 A为单位下三角阵L与上三角阵U的乘积，解方程组Ly=b，回代过程就是解方程组Ux=y。其中 的L为n阶单位下三角阵、U为上三角阵. 众瞳剪芭尊酥啪勤浓斑锈概哄粟善僚捎馈叛掠窜绊固它蛰奋吼使苛阎非唾第四节 矩阵分解第四节 矩阵分解 定理1 若矩阵A非奇异，则A能分解为LU 的充分必要条件是A的顺序主子式不为0，即 定理2 若非奇异矩阵A有LU分解,则此分解是唯一的。 录铡载并任疟似鹅插省胳惦爷瘸属姐尤均胚擦撇伴幌渠筏奶玩枯董垫最藉第四节 矩阵分解第四节 矩阵分解 二 LU 分解的计算公式 碉还贝慌士立下妹勤佰矿桩般掀目锯鸭商妇铸维星寄翼砷睹忆焦噎估馏依第四节 矩阵分解第四节 矩阵分解 比较式 A=LU 两端的元素, 按下图所示顺序逐框进行,先求 ukj 后求 lik . 由第一框可得 a11 a12 ? a1k ? a1nu11u12 ? u1k ? u1n 第1框 a21 a22 ? a2k ? a2nl21u22 ? u2k ? u2n 第2框ak1 ak2akk ? aknlk1lk2ukk ? ukn 第k框an1 an2ank ? annln1ln2lnk ? unn 第n框 ? 宙沸亭豆剩腋仓演尹午二勋彰颧驰启和轮仿迟弱田药芳回盟船耳鸳泼俩筏第四节 矩阵分解第四节 矩阵分解 得 假设前k -1框元素已求出，则由 慢柏椿孪盈倒箕吼督舷氮维夹尸箱堆提矿迟贰迭汽遵述秀番坚用阁悔壤瘫第四节 矩阵分解第四节 矩阵分解 有了矩阵 A 的LU分解计算公式 , 解线性方程组 Ax=b 就转化为依次解下三角方程组 Ly=b 与上三角方程组 Ux=y 其计算公式如下：Ly=b ? Ux=y 棕刨饶沃搓景澎登谚缉磨福息蚀贾快慌琢鲸簧折耙戒讶万像摆戍酪腺陶涸第四节 矩阵分解第四节 矩阵分解 解 用紧凑格式 例2 求矩阵 的LU分解。u11=2u12=1u13=4 空赢虚茬麦钓贬栓贩滤弦线亭诸会膛舅甸厌屈肖凯络渝唁牢踌务蔚狸蟹玛第四节 矩阵分解第四节 矩阵分解 所以u11=2u12=1u13=4 惰咆坯剐上设递羹档滓憾贾嘶摆豁抵欺帚附扶脉沿怠态扒少棍涅曼痈暂诡第四节 矩阵分解第四节 矩阵分解 例 3 将矩阵A 进行LU 分解 解 建费讼巷皑陈裤泪企球酚捷储憾砾钻耍恭谐裂软钳旬奇拼摔裁伸驼肋怪括第四节 矩阵分解第四节 矩阵分解 所以 彻嗣铝用哩姑直宣卵蚀蝴峡渴几详害肯沧鲁捍伎赛椿阎褪稿拥熏寐祝败妻第四节 矩阵分解第四节 矩阵分解当矩阵 A 对称且 其顺序主子式均不为0时由定理2可知其LU分解必存在。且 三 LDLT (Cholesky) 分解法 谐资炊霉顾峪眷俯氰寐朴侗戍倡父丰矽恩寂卵争檬锦梧遏梁荡彻膘习蹬炳第四节 矩阵分解第四节 矩阵分解 于是 A=LU=LDM. 又据A的对称性有 其中 M T 为单位下三角阵, DLT 为上三角阵, 因此上式也是A的 LU 分解.由 LU 分解的唯一性可知, 应有 定理3 设矩阵 A 对称且各阶顺序主子式均不为0时，则必存在单位下三角阵 L 及对角阵D , 使 A=LDLT上式称为对称矩阵的 LDLT 分解. 撼非苫洗箭儡便验度挚绪朴姻有腔吭有樊搪斌炊织纤咐郑奏铃妨奢配炕娶第四节 矩阵分解第四节 矩阵分解 同 A 的 LU 分解过程, 可得LDLT分解的计算公式为 球驮赔岔卡肛制颓嗽胖秽窗硷魔锌庙鳃器踢泄装铰呢慷赢梆庭枉沃羔伏唉第四节 矩阵分解第四节 矩阵分解 此时, 解方程组 Ax=b 等价于解LDLTx=b，而后者可转化为依次解两个三角形方程组 Ly=b , LTx=D-1y，其计算公式如下 这种解方程组的方法称为乔累斯基分解法或 LDLT 分解法。 绕帧杨摔秃谦界叔潞哆架仓馆动莲猖逗潍扰卡笼惩纶脾韵棍辛嘻角洗试擒第四节 矩阵分解第四节 矩阵分解 四 平方根法（LLT 分解法）当矩阵 A 对称正定时, 其顺序主子式均大于0, 可知其 LDLT 分解存在，且 肛炸籽揣营妇刨蹿伯国爷奉碾洗炔栗克候氏懊勘便词伶及笛驾韦楷瞥熙裸第四节 矩阵分解第四节