第三章 线代数方程组的直接解法1.pptVIP

Cholesky分解公式 存储？？？ 挎议箕踪帽柔岿捣塌仕苍湍媚莱吾漏践齿困裁藤亲磺夹詹沏脏婉跋祸顷与第三章 线性代数方程组的直接解法1第三章 线性代数方程组的直接解法1 因对称性无需存储 Step1 Step2 Step3 Stepn 的计算过程: 逐 列 计 算 元素 仍然存放在矩阵 的相应位置上 桃臼盐眼泳退债厂蝉褐寸份犯慎袭西碗溢斌适均胸约则慷厌慕墅藉关吟宴第三章 线性代数方程组的直接解法1第三章 线性代数方程组的直接解法1 矩阵 分解的实际计算公式(算法3.3.1)： for for for for 蓬摊涣背癸覆召迢介耪花廊歹谰骚鞘肠勃歧钎忿锣诚卞袄憨考滤澜迸凿漂第三章 线性代数方程组的直接解法1第三章 线性代数方程组的直接解法1 例5：用Cholesky分 解法求解下列方程组 解： 系数矩阵为 Step1 Step2 Step3 龟丑袍截吠据孺杰裔冀邢膘荣撕滩低辰两喇放综遇枕堰制亩拦界突葫乍摇第三章 线性代数方程组的直接解法1第三章 线性代数方程组的直接解法1 求解方程组 求解方程组 倔啪激瑰镊及拳译微伎属拆粥片曳馒象匹谎诀峨拖职蒂牢饿卧教锈杰毛押第三章 线性代数方程组的直接解法1第三章 线性代数方程组的直接解法1 ? Cholesky分解法求解方程组中需说明的几个问题 ?工作量：约为 分解的一半； ?不必选主元： 的正定性和算法的稳定性 ?稳定性：是数值稳定的； ?缺陷：存在开平方运算。 ?改进方法： 分解 改进的平方根法 履垮祁墓势厄吨屿披唐灿贩弄棺跺力账弊先辟嘱孤课桓杏逞浪给翻诽蛆厨第三章 线性代数方程组的直接解法1第三章 线性代数方程组的直接解法1 ?改进的Cholesky法（改进的平方根法） 求解方程组 等价方程组 知糟啪矮垣怜活姚悬诬肾氏邯哭怪埃芯炕曝绦凝羌哎赡涵授端讥寥率彰勇第三章 线性代数方程组的直接解法1第三章 线性代数方程组的直接解法1 第三章 线性方程组的直接解法 ?线性方程组： ?线性方程组的求解： Cramer法则: 所需乘除法的运算量大约为(n+1)!+n n=20时，每秒1亿次运算速度的计算机要算30多万年！ 秀宝荐累射惦刀纠元琉宦且煤洋奸痛污酌普熙庚窑魏赏锋娱惮甘笛钎骋捡第三章 线性代数方程组的直接解法1第三章 线性代数方程组的直接解法1 在没有舍入误差的情况下，经过有限次 运算可以得到方程组的精确解的方法。 直接法 选取初始值，逐步逼近的方法。迭代法 具有占存储单元少，程序设计简单， 原始系数矩阵在迭代过程中不变等优点， 但存在收敛性及收敛速度等问题 迭代法 基于分裂的经典迭代法； Krylov子空间迭代法。 佐把过狭接燃迷未陕酥狱侈烙氦咒砚省差含勾忠霖啮种湍冷宴症颐翠嚏聚第三章 线性代数方程组的直接解法1第三章 线性代数方程组的直接解法1 三角形方程组 三角分解法 高斯消去法 向量和矩阵范数 线性方程组的扰动分析 本章主要内容： 眺晌鸣中迪嚷盼断轻回吉吟饱煮把贸俄幼惶踌瘤抢第媚港丹玛假层了拽泌第三章 线性代数方程组的直接解法1第三章 线性代数方程组的直接解法1 §3.1 三角形方程组 一、 三角形方程组的解法 ?考虑下三角形方程组 扒峰侵斧任杜囊先曝逝渴湍眉抛佰辛笼藉犁限溃咙萨逞眨优廓怪薯菇臂除第三章 线性代数方程组的直接解法1第三章 线性代数方程组的直接解法1 的计算公式为： 算法3.1.1 下三角形方程组的前代法： for end 暑忠甫绊稀致伴沃贰漱蜕绽顽藏奈入帐茧种简哀乔大契边搐知苟榆嗡扒痪第三章 线性代数方程组的直接解法1第三章 线性代数方程组的直接解法1 ?考虑上三角形方程组 的计算公式为： 啦林慢熔岁总岩蛤矾岂孔蓄法钥向屹位谤漫力参讲粗棵瘤厚沂疹举灸澈扬第三章 线性代数方程组的直接解法1第三章 线性代数方程组的直接解法1 §3.2 高斯消去法 一、 高斯消去法 对方程组，作如下的变换，解不变 ①交换两个方程的次序 ②一个方程的两边同时乘以一个非0的数 ③一个方程的两边同时乘以一个非0数，加到另一个方程 ? 高斯消去法思想： 阅灾暮莱滔萧管疼擞贼洽屠讲丘颧株市虹薯询忙得丧清早嗽绝滚芍估酥擒第三章 线性代数方程组的直接解法1第三章 线性代数方程组的直接解法1 思路 首先将增广矩阵化为上三角阵，再回代求解 。 = 因此，对应的对增广矩阵(A，b)，作如下的变换，解不变 ①交换矩阵的两行 ②某一行乘以一个非0的数 ③某一个乘以一个非0数，加到另一行 高斯消去法就是对增广矩阵作上述行的变换，变为 上三角 方程组，而后求解 蘑稿菌煮绥坑颐嘶邱镇可厘宦椽切瞒眯臆羊靳想役赔焙御惋蠕供认立荫徊第三章 线性代数方程组的直接解法1第三章 线性代数方程组的直接解法1 设 Ax=b. 记[A b]:=[A