第二章线性方程组的直接解法解析.ppt

下载文档 降价啦

21
0
约5.68千字
约 49页
2016-03-20 发布于湖北
举报
版权申诉
保障服务

第二章线性方程组的直接解法解析.ppt

1、本文档共49页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

* 例4 利用系数矩阵的LU分解, 求解方程组解：LU分解的紧凑格式为上一页下一页返回 * 则求解原方程组等价于求解下面两个方程组Ly=b及Ux=y 上一页下一页返回 * 思考：在科学与工程计算中，经常遇到求解三对角方程组的问题： * 二、三对角方程组追赶法若A满足Gauss消去法可行的条件，则可用LU分解法求解其中：上一页下一页返回 * 解方程组Ax=d分为两步，即求解Ly=d和Ux=y，计算公式如下：上述方法为求解三对角方程组的“追赶法”，也称Thomas算法. 上一页下一页返回追赶法公式简单，计算量和存储量都小，整个求解过程只需要5n-4次乘除运算。 * 上一页下一页返回三、平方根法 ——对称正定矩阵的分解法将对称正定阵 A 做 LU 分解记为 * 定理设n阶对称正定矩阵A，则存在唯一的单位下三角阵L及对角阵D 使得。上一页下一页返回定理设矩阵A对称正定，则存在唯一的对角元为正的下三角阵 L1，使得。利用乔列斯基（Cholesky）分解式来求解Ax=b的方法也称Cholesky方法或平方根法. 称为对称正定矩阵A 的乔列斯基分解 * §3 方程组的性态与误差估计上一页下一页返回一、矩阵的条件数例，考查以下三个方程组及其准确解其准确解其准确解其准确解可以看到，后两个方程组与第一个方程组相比，系数矩阵或右端向量仅有0.0005以下的误差，但准确解却相差很大。对这样的方程组，无论用多么稳定的算法求解，一旦计算中产生误差就使解面目全非，所以该方程组的性态很差。 * 上一页下一页返回定义1 若方程组Ax=b的系数矩阵A与右端向量b的微小变化（小扰动）,将引起解向量x产生巨大变化，则称此方程组为病态方程组，其系数矩阵A称为病态矩阵，否则称Ax=b为良态方程组，称A为良态矩阵 . 方程组的病态程度与Ax=b对A和b的扰动的敏感程度有关。 * 求解时，A 和的误差对解有何影响？ ? 设 A 精确，有误差　，得到的解为，即绝对误差放大因子又相对误差放大因子上一页下一页返回（教材P25-27） * ? 设精确，A有误差　，得到的解为，即 (只要? A充分小，使得是关键的误差放大因子，称为 A的条件数，记为cond (A) , 越则 A 越病态，难得准确解。大上一页下一页返回 * 注: ? cond (A) 的具体大小与 || · || 的取法有关，但相对大小一致。 ? cond (A) 取决于A，与解题方法无关。常用条件数有： cond 1(A) cond ?(A) cond2 (A) 特别地，若 A 对称，则条件数的性质： ? A可逆，则 cond p (A) ? 1； ? A可逆，? ? R 则 cond (? A) = cond (A) ； ? A正交，则 cond 2 (A) =1； ? A可逆，R正交，则 cond 2 (RA) = cond 2 (AR) = cond (A)2 。上一页下一页返回 * 精确解为例5 计算cond 2(A) 。 A?1 = 解：考察 A 的特征根 39206 1 ? 测试病态程度：给　一个小扰动，其相对误差为此时精确解为 2.0102 200% 上一页下一页返回 * 例6 Hilbert 阵 cond ? (H2) = 27 cond ?(H3) ? 748 cond ? (H6) = 2.9 ? 106 cond? (Hn) ? ? as n ? ? 注：一般判断矩阵是否病态，并不计算A?1，而由经验得出。 ? 行列式的值很大或很小（如某些行、列近似相关）； ? 元素间的数量级相差大，且无规则； ? 主元消去过程中出现小主元； ? 特征值相差大数量级。上一页下一页返回 * 二、方程组解的误差估计定理 ? 上一页下一页返回 * 上一页下一页返回解小结高斯消去法，选主元技术；三角分解法：LU分解法、追赶法（三对角方程组）、平方根法（对称正