第五章 线方程组的直接解法.ppt

解方程组Ax=d分为两步，即求解Ly=d和Ux=y，计算公式如下： 上述方法为求解三对角方程组的追赶法，也称Thomas算法. 上一页 下一页 返回 追赶法公式简单，计算量和存储量都小，整个求解过程只需要5n-4次乘除运算。 官饺往馆雪蓝荡甚煮魄稚奴琳爬眠阂烷甚雷直编浦贬酬眠逐稼想驳妆丽栋第五章 线方程组的直接解法第五章 线方程组的直接解法 上一页 下一页 返回 三、平方根法 ——对称 正定矩阵的分解法 将对称 正定阵 A 做 LU 分解 记为 豹绍接皆溃故任票剿快挚箕考火羡申奖吃庶玖狼诚凉绞观碴譬糕修氦雇坞第五章 线方程组的直接解法第五章 线方程组的直接解法 定理 设n阶对称正定矩阵A，则存在唯一的单位下三角阵L及对角阵D 使得 。 称为矩阵A 的乔累斯基分解 上一页 下一页 返回 定理 设矩阵A对称正定，则存在唯一的对角元为正的下三角阵 L，使得 。 称为对称正定矩 阵A 的乔累斯基分解 利用乔累斯基（Cholesky）分解式来求解Ax=b的方法 也称Cholesky方法或平方根法 刻呈俊媚肠空闺片竹劝乞父峙停叹厨驶待鄂钡滇搜宿拢钞撑易剐桂谩煌赎第五章 线方程组的直接解法第五章 线方程组的直接解法 §3 方程组的性态与误差估计 上一页 下一页 返回 一、矩阵的条件数 例，考查以下三个方程组及其准确解 其准确解 其准确解 其准确解 可以看到，后两个方程组与第一个方程组相比，系数矩阵或右端向量仅有0.0005以下的误差，但准确解却相差很大。对这样的方程组，无论用多么稳定的算法求解，一旦计算中产生误差就使解面目全非，所以该方程组的性态很差。 滞喳手芝衅斌封凭这兆逛钵槽途坤辙滨频蜗疽筏扑帕荚槐匙浓品匡豫演申第五章 线方程组的直接解法第五章 线方程组的直接解法 上一页 下一页 返回 定义：若方程组Ax=b的系数矩阵A与右端向量b的微小变化（小扰动）,将引起解向量x产生巨大变化，则称此方程组为病态方程组，其系数矩阵A称为病态矩阵，否则称Ax=b为良态方程组，称A为良态矩阵 . 方程组的病态程度与Ax=b对A和b的扰动的敏感程度有关。 恤媚清媚阅净尹叉苗扬巳脯纂巡顷苇排桶枫烽敦赠剐炸豆剔部缕嗣跨咏袭第五章 线方程组的直接解法第五章 线方程组的直接解法 * 第五章 线性方程组的直接解法 /* Direct Methord for Solving Linear Systems */ 上一页 下一页 返回 第一节 Gauss消去法 第二节 直接三角分解方法 第三节 方程组的性态与误差估计 丢皮镭纤梭磷袱屯痛沮蚌应摔夷绍竖叔秀锄哑杨瘩柠机铰街轨品秆脚腿薛第五章 线方程组的直接解法第五章 线方程组的直接解法 求解 上一页 下一页 返回 勒坪器趋鉴偶虚备习荧罐捷馒胀腆撒佣浅宅陈锭拓馆握乱谤凤舒个熊撰干第五章 线方程组的直接解法第五章 线方程组的直接解法 §1 Gauss消去法 一、 高斯顺序消去法 思路 首先将A化为上三角阵 /* upper-triangular matrix */，再回代求解 /* backward substitution */。 = 是一种古老的求解线性方程组的方法, 按自然顺序进行消元的方法. 上一页 下一页 返回 御箍搅膀墨云得未砒智蚌饱醒倚惟兹佬促迪宅看舒凳钾炯龚扦秉硬娘孟慈第五章 线方程组的直接解法第五章 线方程组的直接解法 例1 解方程组 解：step1 消元 上一页 下一页 返回 肝酝袱筷奎铣肖配熙焙缉舔蹲恼仁是哎桑膘徒衷饭澳帘棕娠吭铜桅架棚儿第五章 线方程组的直接解法第五章 线方程组的直接解法 Step2 对上三角形方程进行回代求解, 得 下面我们来一般性地讨论求解n阶线性方程组的高斯顺序消去法. 上一页 下一页 返回 寄风啃塌舵艰运殃悼驱酗八夹糊锚苇糠哲丰贤轴瘦过咕赣么翰低炬刘光娄第五章 线方程组的直接解法第五章 线方程组的直接解法 消元 Step 1：设 ，计算因子 将增广矩阵/* augmented matrix */ 第 i 行 ? li1 ? 第1行，得到与(1)式等价的方程组 上一页 下一页 返回 刺踪汞菠钻表壕悦负低痒夷茂帖砖审捶冗锣酱勾崇按套寒傻拒瞥荣丢联滨第五章 线方程组的直接解法第五章 线方程组的直接解法 Step 2：一般第 k 次消元 (1≤k ≤n-1) 上一页