- 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分解法、 追赶法(三对角方程组)、 平方根法(对称正
文档评论(0)