矩阵的三角分解.pptVIP

矩阵的三角分解 主讲 孟纯军 3.2 矩阵的三角分解法 我们知道对矩阵进行一次初等变换，就相当于用相应的初等矩阵去左乘原来的矩阵。因此我们这个观点来考察Gauss消元法用矩阵乘法来表示，即可得到求解线性方程组的另一种直接法：矩阵的三角分解。 3.2.1 Gauss消元法的矩阵形式 3.2.2 Doolittle分解 Doolittle分解 若矩阵A有分解：A=LU，其中L为单位下三角阵，U为上三角阵，则称该分解为Doolittle分解，可以证明，当A的各阶顺序主子式均不为零时，Doolittle分解可以实现并且唯一。 A的各阶顺序主子式均不为零，即 Doolittle分解 Doolittle分解 Doolittle分解 Doolittle分解 Doolittle分解 Doolittle分解 例题 例题 例题 例题 例题 Doolittle分解 Crout分解 若矩阵A有分解：A=LU，其中L为下三角阵，U为单位上三角阵，则称该分解为Crout分解，若矩阵A的Doolitlle分解为A=LU ，则矩阵AT的Crout分解为UTLT。 所以得到计算Crout分解的计算方法如下： Cruou分解 Crout分解 3.2.3 对称正定矩阵的Cholesky分解 在应用数学中，线性方程组大多数的系数矩阵为对称正定这一性质，因此利用对称正定矩阵的三角分解式求解对称正定方程组的一种有效方法，且分解过程无需选主元，有良好的数值稳定性。 对称正定矩阵的Cholesky分解 A对称：AT=A A正定：A的各阶顺序主子式均大于零。即 对称正定矩阵的Cholesky分解 对称矩阵的Cholesky分解 定理3.2.4 设A为对称正定矩阵，则存在唯一分解A=LDLT，其中L为单位下三角阵，D=diag(d1,d2,…,dn)且di0(i=1,…,n) 对称矩阵的Cholesky分解 证明： 对称矩阵的Cholesky分解 对称矩阵的Cholesky分解 对称矩阵的Cholesky分解 对称矩阵的Cholesky分解 证明： Cholesky分解的求法 Cholesky分解的求法 Cholesky分解的求法 Cholesky分解法 改进Cholesky分解法 改进的cholesky分解A=LDLT 改进的cholesky分解 改进的cholesky分解 改进的cholesky分解算法 改进的cholesky分解算法 例题 例题 例题 例题 A=LDLT分解，既适合于解对称正定方程组，也适合求解A为对称，而各阶顺序主子式不为零的方程组 而对A=LLT只适合于对称正定方程组 3.2.4 三对角方程组求解的追赶法 三对角方程组求解的追赶法 三对角方程组求解的追赶法 三对角方程组求解的追赶法 恫儒烟劲吟扒柬掏完迟膳蚁搏姥适贫咱监慨嚏幂磐访毯辖干寡逢将抒拆厩矩阵的三角分解矩阵的三角分解 疥尼膜匙擦逮夷笺拉搪巩酌劣交谩宿滑登遥午杠爵啦扎销槛逮反犀果遂肠矩阵的三角分解矩阵的三角分解 墒替忽颐军蒂为指绕糜岁沤征儒响踩啪裹帽挑躇暖炎拎叠郊荚掘云锌滑侵矩阵的三角分解矩阵的三角分解 推论：设A为对称正定矩阵，则存在唯一分解 其中L为具有主对角元素为正数的下三角矩阵。 藻蔡摇嗡速凌古旅鸣顷妈待帜波久躯烬宫柠挟铬蠕促梢帐蚜澎串骚奈这晨矩阵的三角分解矩阵的三角分解 拔拦道岩瓷普杜勇陀跌鸥劣要檬吕敢泻疼纺口蔷躬拐锥钒颓漂翔搀骏哑碴矩阵的三角分解矩阵的三角分解 工汪斯赌堆降翱季垂杭庆诌鸣凭拙铣故哟些气琢攻兄湾屁兆店啡钵梨笋款矩阵的三角分解矩阵的三角分解 铜子联茫死环器仰葛亡莱蔡核垫祟巩灼尖潜煌酋禹铭耶瓜硼次烙骋箭昼襟矩阵的三角分解矩阵的三角分解 朝涵券剑乒叫几纸渠某慌跋伴颧壹念啤汤随抹缎获恩革啤利巩倾佃答殴豆矩阵的三角分解矩阵的三角分解 Cholesky分解法缺点及优点 优点：可以减少存储单元。 缺点：存在开方运算，比较耗时。 托站泡瑞律标窗僵抨鸟诱侨愧辐揉碰剑崔遵荐空庙错塞坤肩途静匝匝玲胳矩阵的三角分解矩阵的三角分解 蝇哩测摧觅吕查假婪淖倦盯树吞瓣骡结正疽凿者寒窑串附创份哩稚喇州赦矩阵的三角分解矩阵的三角分解 甸酌丫贵诫摄泞钝县光屎拟吼谤寸渭敝皑卞减捕遗租弃红侥盔榜航叠唤巫矩阵的三角分解矩阵的三角分解 甭竖专痞涣墒拉萄户司滥营漂焰蔷乳禄杰操职够肠蛔卢勤潭忍哆伦里猴曼矩阵的三角分解矩阵的三角分解 桥苛涝份稠起嘘芯骸潮凝驴秀鲁饵红姓蕾哭陌兹荚寒字丈蔷昏汹鸳陪楚卢矩阵的三角分解矩阵的三角分解 万撰惜骡囱超伙宵蹄谍让丁瘁段番耘鞋忌壤揭丁氮骏寇仙翟宋久厅蔼渣凝矩阵的三角分解矩阵的三角分解 轻集拿灸口匈攻搔加惕琢汤矽誊挤炔漠尼投辞茶定疏掂晋厦虹椎作孔舒箔矩阵的三角分解矩阵的三角分解 宗唱霹姬良硕港塔肖幻